
==== Front
Environ Sci Pollut Res Int
Environ Sci Pollut Res Int
Environmental Science and Pollution Research International
0944-1344
1614-7499
Springer Berlin Heidelberg Berlin/Heidelberg

39167144
34758
10.1007/s11356-024-34758-w
Research Article
Migration of copper(II) ions in humic systems—effect of incorporated calcium(II), magnesium(II), and iron(III) ions
Klučáková Martina klucakova@fch.vut.cz

Enev Vojtěch
https://ror.org/03613d656 grid.4994.0 0000 0001 0118 0988 Brno University of Technology, Faculty of Chemistry, Purkyňova 118, 612 00 Brno, Czech Republic
Responsible Editor: Kitae Baek

21 8 2024
21 8 2024
2024
31 40 5299653007
26 6 2024
16 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The mobility of heavy metals in natural soil systems can be affected by the properties and compositions of those systems: the content and quality of organic matter as well as the character of inorganic constituents. In this work, the diffusion of copper(II) ions in humic hydrogels with incorporated calcium(II), magnesium(II), and iron(III) ions was investigated. The methods of instantaneous planar source and of constant source were used. Experimental data yielded the time development of the concentration in hydrogels and the values of effective diffusion coefficients. The coefficients include both the influence of the hydrogel structure and the interaction of diffusing particles with the hydrogel. Our results showed that the presence of natural metal ions such as calcium, magnesium, or iron can strongly affect the diffusivity of copper in humic systems. They indicate that the mobility of copper ions depends on their concentration. The mobility can be supported by higher contents of copper in the system. While the incorporation of Ca and Mg resulted in the decrease in the diffusivity of copper ions, the incorporation of Fe(III) into humic hydrogel resulted in an increase in the diffusivity of Cu(II) in the hydrogel in comparison with pure humic hydrogel.

Keywords

Humic acid
Hydrogel
Copper
Calcium
Magnesium
Iron
Diffusion
Immobilization
Brno University of TechnologyOpen access publishing supported by the National Technical Library in Prague.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Humic substances as the major components of soil organic matter play an important role in the bioavailability, toxicity, and transport of different organic and inorganic pollutants in nature (Ding et al. 2019; Hayes et al. 2008, 2012; Sureshkumar et al. 2013). Transport processes are therefore important for understanding the complex behaviour of natural systems and the role of their constituents. The environmental behaviour of metal ions is strongly affected by their binding to humic substances (Klučáková and Pekař 2003b, 2004; Paradelo et al. 2012; Saito et al. 2005).

Humic substances as major constituents of soil organic matter form complexes with metal ions and tend to control their speciation, toxicity, and bioavailability (Gao et al. 2015; Kogut and Voelker 2001; Tipping et al. 2008). Because of the high affinity of metal ions to humic substances, the formed complexes are stable. Martyniuk and Wieckowska (2003) determined the highest affinity to humic acids for Pb, Ag, Hg, Cu, Ba, and Cd. Ding et al. (2019) investigated the following order: Pb > Cu > Ni > Cd > Zn. In contrast, Dinu (2015) demonstrated that Fe and Al ions have higher conditional stability constants than Pb, Cu, and Zn.. van Oort et al. (2020) monitored Zn and Pb contents in 173 soil samples and considered Pb as a less-mobile metal in comparison with Zn. Metal ions may distribute among more types of binding sites. A few of them can be kinetically labile under specific experimental conditions, such as weak bidentate sites for Cu (Shi et al. 2016). The transport of metal–humic complexes can be enhanced when the complexes are mobile in the soil solution or retarded when the organic matter is attached to a solid matrix (Bryan et al. 2005; Paradelo et al. 2012).

The behavior of metal ions and their complexes with humic substances in soils are strongly affected by other different factors such as the pH and ion composition of the solution in pores. The mobility of metal ions as well as their immobilization by complexation with humic substances should thus be influenced by the presence of other ions and substances. Knowledge of the migration behavior of metal ions and their humic complexes under different circumstances is crucial for an understanding of this complex problem.

Copper is an essential nutrient element for plant growth, but its excessive amounts can result in its accumulation in soil and cause risks to the ecosystem as well as to human health (Christl et al. 2005; Wu and Hendershot 2010). Wu and Hendershot (2010) investigated the accumulation of copper in soil, and its toxicity. They showed that increased contents of calcium can suppress the uptake of copper. The effect of calcium and other ions on the complexation of copper in soil and copper toxicity was described in several studies. Kogut and Voelker (2001) studied the copper binding behavior of terrestrial humic substances in seawater. They found that copper binding was weaker than the binding reported for freshwater, possibly indicating the effects of calcium and magnesium, or (in general) ionic strength.

In contrast, Hamilton-Taylor et al. (2002) reported that the presence of seawater concentrations of calcium and magnesium had no detectable effect on the binding of copper, although its binding decreased generally with increasing ionic strength. Their results indicated that the competition for humic binding sites between calcium and magnesium may be important for weakly-binding metals but not for strongly-binding metals such as copper. Kang et al. (2010) studied the adsorption of copper on humic-like substances extracted from activated sludge. Their results confirmed competitive adsorption between free copper and calcium at lower calcium concentrations. Higher calcium concentrations resulted in its precipitation. Gao et al. (2015) observed a decrease in copper binding on dissolved organic matter in the presence of calcium as a result of the decrease in available binding sites.

Saito et al. (2005) analyzed copper binding in a system containing humic acids and goethite. Binding was enhanced in the system containing goethite. Ding et al. (2019) fractionated humic acids by ferrihydrite and studied the binding characteristics of heavy metals. Their results showed that changes in the reactivity of humic acids during their fractionation on ferrihydrite affected the kinetic behavior of heavy metals in solution. Similarly, Lippold et al. (2007) concluded that the competitive influence of iron on the transport of metal ions in water is kinetically controlled. Farina et al. (2018) precipitated ferrihydrite and goethite organo-mineral complexes with humic acids as an analogue of natural organic matter and measured the adsorption of copper. Their results showed that copper adsorption was additive if the mineral fraction predominated in the prepared complexes, and non-additive if the organic fraction predominated. Muller and Cuscov (2017) studied the copper-binding capacity of iron-rich humic substances during transport in marine. The interactions resulted in the decrease in the copper-binding capacity of humic substances. Christl et al. (2001, 2005) demonstrated that adsorption isotherms of copper are independent of concentrations of humic and fulvic acids. In contrast, the binding capacity was affected by pH and ionic strength.

As can be seen, many authors deal with the effects of some ions such as iron, calcium, magnesium, and sodium on the binding of copper or other heavy metals in systems containing humic substances as representatives of natural organic matter. Many authors are interested in adsorption experiments. Studies of interactions of metal ions with humic substances or organic matter directly during their diffusion (or, generally, transport) are relatively scarce. Paradelo et al. (2012) studied the transport of humic substances and copper ions in pore water by means of experiments in water saturated columns filled by quartz sand. They showed that the transport of copper was coupled with humic substances at neutral pH and occurred mostly as free copper ions at pH = 5. Some diffusional characteristics of metal ions in several types of hydrogels were published in connection with the monitoring of metal ions in nature by means of the DGT technique (Scally et al. 2006; Zhang and Davison 1999).

Our approach to the mobility of metal ions is different. We are interested in the diffusion of metal ions through a gelled humic medium, which can be considered as a model of natural soil systems (Klučáková and Pekař 2003b, 2004, 2009). This means that metal ions diffuse through the hydrogel and the front of diffusing ions penetrates into the hydrogel containing non-occupied binding sites. Metal ions can thus interact immediately with humic acids and be (partially) immobilized. These interactions between metal ions and humic binding sites can achieve an equilibrium between immobilized and free mobile ions able to continue in motion.

Our previous work has shown that the binding capacities of humic substances are affected by their used form (Klučáková 2014; Klučáková and Pekař 2003a) and content of active sites (Klučáková 2014; Klučáková and Kalina 2015). Humic substances can be prepared in hydrogel form in several ways. One possibility is the precipitation of humic substances dissolved in alkaline medium by an acid solution, as used in our previous research (Klučáková 2014; Klučáková et al. 2014; Klučáková and Kalina 2015; Klučáková and Pekař 2003b, 2004, 2009). The cross-linking of hydrogels can also be achieved by means of multivalent ions (Klučáková and Pekař 2009). This approach can be used for the preparation of hydrogels containing (besides humic substances) ions characteristic of natural systems (e.g., Fe, Mg, Ca) and better simulate real circumstances in soils.

The main goal of this study was to investigate the diffusion of copper(II) ions in humic hydrogels enriched by calcium(II), magnesium(II), and iron(II) ions; determine the diffusion coefficients; and consider the effect of added ions in comparison with “pure” humic hydrogel. Another aspect investigated in this study was the dissolving medium.

Traditionally, and for many years, sodium or potassium hydroxide have been used as alkaline solutions in the isolation of humic substances from solid-phase source materials, including isolation procedures used by International Humic Substances Society (De Nobili et al. 1990; Thurman and Malcolm 1981; Watanabe and Kuwatsuka 1991). Some researchers (Conte and Piccolo 1999; Gerzabek et al. 1997; Maia et al. 2008; Novotny et al. 1999; Nuzzo et al. 2013; Tatzber et al. 2007) used Na4P2O7 solution as an alternative extraction agent or in combination with the traditional NaOH solution. Tatzber et al. (2007) characterized in detail soil humic acids isolated using three different extractants including sodium hydroxide and tetrasodium pyrophosphate. Extraction using NaOH exhibited lower yields and probably extracted a different fraction of soil organic matter in comparison with the other extraction agents. They noted that sodium hydroxide can induce some reactions in organic matter like the Kolbe reaction, autooxidation, the breakdown of humic (macro)molecules, and the condensation of amino-carbonyls, as well as reactions like the hydrolysis of ester groups or the base-catalysed aldol condensation reaction (Head and Zhou 2000; Rosa et al. 2005; Tatzber et al. 2007).

In general, tetrasodium pyrophosphate is considered a milder agent for the extraction of humic acids. However, it might also lead to changes in the humic structure due to the incorporation of orthophosphate in humic fractions of lower molecular weights and pyrophosphate in humic fractions of higher molecular weights (Francioso et al. 1998; Tatzber et al. 2007). It is unknown whether differences in the analyses of humic acids extracted by different agents arose because of the extractions of different humic fractions or because of alterations to the extracted humic samples caused by the extractants themselves (Francioso et al. 1998; Head and Zhou 2000; Rosa et al. 2005; Tatzber et al. 2007). Although, Tatzber et al. (2007) concluded that NaOH extraction did not influence the analysis, we decided to use both NaOH and Na4P2O7, the solutions most frequently used for the extraction of humic substances, in our study. The dissolution of humic acids in these two different solutions and their precipitation/gelation using different agents such as acidic solutions and multivalent ions can provide several types of humic hydrogels. Incorporated ions (Ca, Mg, and Fe) as well as potential changes in the humic structure caused by different solvents should affect the mobility of copper(II) in hydrogels. We believe that investigation of the interactions of copper(II) ions and their diffusion simultaneously and directly in these hydrogels can provide knowledge useful for predicting the toxicity and mobility of trace metals in nature.

Materials and methods

Preparation of humic acids

Humic acids were isolated from South Moravian lignite (Czech Republic). They were prepared by means of an alkali extraction procedure using a mixture of 0.5 M NaOH and 0.1 M Na4P2O7 (1:1). This initial step was followed by centrifugation, the precipitation of humic acids from the alkali-soluble portion by a solution of HCl (pH < 2), washing, purification, and drying. The procedure is described elsewhere (Klučáková 2014; Klučáková and Kalina 2015; Peuravuori et al. 2006). More details on the chemical structure of the initial lignite matrix, as well as that of the isolated humic sample, can be found in previous works (Peuravuori et al. 2006, 2007). The basic characteristics of the humic sample are listed in Table 1.Table 1 Basic characteristics of humic acids

Ca (at. %)	Ha (at. %)	Oa (at. %)	Na (at. %)	Sa (at. %)	Total acidity (mmol g−1)	COOH (mmol g−1)	
32.3	39.0	27.0	1.4	0.3	4.9	3.7	

Preparation of humic hydrogels

The powdered humic acids were dissolved in 0.5 M NaOH solution (hydrogels H) or 0.1 M Na4P2O7 solution (hydrogels P). The ratio between the solid humic sample and the solution was 8 g per 1 dm3. The solutions were used for the preparations of all hydrogels as described below.

Basic hydrogels (H and P) were prepared by the precipitation method—the hydrogels were created after collapse of the humic sols, when their ionic groups were neutralized, and electrostatic repulsion was suppressed. The solutions of humic acids in hydroxide and pyrophosphate were acidified by concentrated HCl solution up to a pH value close to 1. The precipitated humic hydrogels were washed repeatedly with deionized water and centrifuged until Cl− ions were removed.

Hydrogels with incorporated metal ions were prepared by precipitation of the hydroxide and pyrophosphate solutions by 1 M CaCl2, 1 M MgCl2, and 1 M FeCl3. The volume ratio between humic solution and precipitation solution was 1:1. It is known that cations of weak hydroxides can hydrolyze in water, which causes acidic behavior in solutions of their salts. In the case of prepared hydrogels, we assume that the equilibrium of this reaction is shifted to the left as a result of the liberation of H+ ions from humic acids (Klučáková and Pekař 2009). The precipitations of hydroxides were not observed in the preparation of hydrogels (in contrast to blank experiments, in which NaOH solutions were added to salt solutions). The precipitated hydrogels were washed repeatedly by deionized water and centrifuged until Cl− ions were removed. They were designated as H-Ca, H-Mg, and H-Fe (hydrogels based on NaOH solution); and as P-Ca, P-Mg, and P-Fe (hydrogels based on Na4P2O7 solution).

Hydrogels were characterized by means of FT-IR spectrometry (Nicolet iS5, Thermo Scientific). They were dried and the spectra of xerogel (solid matter) were collected using KBr pellets as described previously (Haberhauer and Gerzabek 1999; Tatzber et al. 2007).

Diffusion from an instantaneous planar source

The hydrogels were pressed gently into glass tubes (length 3 cm and diameter 1 cm). A circular piece of filter paper (diameter 1 cm) was sunk into the 1 M CuCl2 solution for 1 min and then added to one side of the tube filled by hydrogel. The tube was wrapped with parafilm and aluminum foil to prevent the hydrogel drying. The durations of the diffusion experiments were 5, 24, 48, and 72 h, respectively. Then, the glass tube was removed, the hydrogel was sliced, and both the paper and hydrogel slices were extracted separately in 1 M HCl solution (Klučáková et al. 2014; Klučáková and Kalina 2015). The concentration of Cu2+ ions in leachate was determined by means of UV/VIS spectrometry (Hitachi U3900H). The obtained data were used to compute the concentration profiles of metal ions in the tubes and their diffusion fluxes (the total amounts diffused into the hydrogel). All experiments were performed at laboratory temperature (25 ± 1 °C) and in triplicate. Data are presented as average values with standard deviation bars.

Diffusion from a constant source

The hydrogels were packed gently into glass tubes (length 5 cm and diameter 1 cm). The tubes were immersed in a saturated solution of CuCl2 (150 cm3). Vessel lids were sealed with parafilm to prevent evaporation. The durations of the diffusion experiments were 5, 24, 48, and 72 h. Then, each tube was taken out from the solution, the glass tube was removed, the hydrogel was sliced, and both the paper and hydrogel slices were extracted separately in 1 M HCl solution (Klučáková and Pekař 2003b, 2004). The concentration of metal ions in leaches was determined by means of UV/VIS spectroscopy (Hitachi U3900H). The obtained data were used to compute the concentration profiles of metal ions in the tubes and their diffusion fluxes (the total amounts diffused into the hydrogel). All experiments were performed at laboratory temperature (25 ± 1 °C) and in triplicate. Data are presented as average values with standard deviation bars.

Results and discussion

Humic hydrogels

Prepared hydrogels were characterized by means of FT-IR spectrometry. The spectra shown in Figs. 1 and 2 show that the use of different solvents resulted in changes of spectra. In spite of the fact that humic acids were isolated from lignite by a mixture of NaOH and Na4P2O7, their re-dissolving in one of these solvents can influence the character of the gelled sample. As mentioned above, tetrasodium pyrophosphate is considered a milder agent for the extraction of humic acids, while sodium hydroxide is considered a stronger extractant which can induce some reactions in organic matter (Tatzber et al. 2007). According to previous studies (Francioso et al. 1998; Head and Zhou 2000; Rosa et al. 2005; Tatzber et al. 2007), differences in the analyses of humic acids extracted by different agents result from the extraction of different humic fractions or from alterations of the extracted humic samples by the extractants themselves. In our case, changes to the spectra should be caused only by an alteration of the humic sample. Changes were recorded mainly in the area < 1800 cm−1, which can be considered as the region exhibiting the predominant characteristics of humic acids, where almost all vibrational information on humic backbones can be identified (Chen et al. 2015). The bands around 1620 cm−1 (a C = O band with a possible contribution of nitrates and nitrites) and 1450 cm−1 (aliphatic C-H) were stronger in the case of xerogel based on hydroxide, and the spectrum showed a higher resolved fingerprint area. In this area, the OH deformation bands (1410–1260 cm−1), ester bands (1330 and 1050 cm−1), and fingerprint bands of aromatic structures (1225–950 cm−1) can be found. It should be noted that the area around 1000 cm−1 can be influenced by inorganic impurities in the samples. Similar differences were observed in spectra published by Tatzber et al. (2007). The incorporation of metal ions into humic samples influenced the spectra mainly in the same area (< 1800 cm−1). Stronger bands around 3600 cm−1 in the spectra of the H-Ca and H-Mg hydrogels can be assigned to free OH groups (although precipitations of hydroxides were not observed in the preparation of hydrogels). Changes in spectra differed both for individual metal ions and the solvent used for the preparation of hydrogel. In contrast, the common characteristic of xerogels with incorporated metal ions was the practical disappearance of the band of carbonyl vibration in carboxylic groups, aldehydes, ketones, and esters (~ 1700 cm−1) (or it was observed only as a shoulder). Two strong bands (around 1600 and 1400 cm−1), which were observed in the spectra of xerogels with incorporated metal ions, may be attributed to the asymmetric and symmetric stretches stretching of COO− groups, respectively (Paradelo et al. 2012; Sowers et al. 2018). This can indicate that binding through COO− groups may be a mechanism for the incorporation of metals into humic structure. Although some of the bands overlap with uncertain assignments, it can be seen that the incorporation of metal ions into the humic sample resulted in specific changes in the measured spectra and that the characteristic bands changed with the addition of different metals.Fig. 1 FT-IR spectra (transmittance) of xerogel samples based on sodium hydroxide: H (black), H-Ca (blue), H-Mg (green), and H-Fe (red)

Fig. 2 FT-IR spectra (transmittance) of xerogel samples based on tetrasodium pyrophosphate: P (black), P-Ca (blue), P-Mg (green), and P-Fe (red)

Diffusion from an instantaneous planar source

This method is based on the placement of a small, predefined amount of diffusing substance as a diffusion source on the surface of the diffusion medium. In this study, the source was the filter paper saturated by diffusing particles and applied on the frontal circular surface of the hydrogel in the tube. The particles then diffuse into the hydrogel and their surface concentration decreases (because the source is exhausted relatively soon). At equilibrium, the homogeneous distribution of metal ions in the hydrogel is achieved. This means that the concentration of metal ions is then constant and independent of the distance from the interface. The concentration profile in equilibrium is thus a line parallel with the x-axis. The concentration in the hydrogel is much lower in comparison with the method of constant source discussed below. Therefore, in this study, the potential dependence of diffusivity on the content of diffusing particles in the hydrogel could be investigated. The mathematical model applied to the data was based on the presumption that the diffusion was not close to equilibrium, and that the diffusion trajectories were shorter that the length of the hydrogel. Simultaneously, the amount of metal ions applied on the interface was known and no metal ions were present in the hydrogel at the start of the experiment. The final solution of second Fick’s second law is (Klučáková 2022; Klučáková et al. 2014; Klučáková and Kalina 2015):1 c=nSπDeftexp-x24Deft,

where c is the concentration of diffusing particles, n is the total mass of the diffusing substance applied in the form of a narrow pulse, S is the cross-section available for diffusion, Def is the effective diffusion coefficient, t is time, and x is the distance from the interface between the filter paper and the hydrogel. The equation is usually applied in its logarithmic form:2 lnc=lnnSπDeft-x24Deft,

which can be used for the calculation of the effective diffusion coefficient Def. In Fig. 3, the example of experimental data fitted by Eq. 2 is shown. The effective diffusion coefficients were calculated from the slopes of the lines. The obtained values are listed in Tables 2 and 3. They are average values for all diffusion times (durations of experiments) in three replications and include both the influence of the pore hydrogel structure and the chemical interactions of copper with the hydrogel. The pore structure can be characterized by the porosity φ and the tortuosity τ (Klučáková 2022; Klučáková et al. 2014; Klučáková and Kalina 2015). The porosity is the proportion of the pore volume in the total volume of hydrogel. If we assume that particles can diffuse only through pores in the hydrogel and that diffusion into the (solid) hydrogel network is negligible, the area of hydrogel accessible for them is smaller in comparison with a homogeneous medium.Fig. 3 Example of data fitting for the method of instantaneous planar source—Eq. 2: H hydrogel after 24 h (blue) and P hydrogel after 24 h (red)

Table 2 Effective diffusion coefficients determined by Eq. 2 for hydrogels based on hydroxide

Def, H (m2 s−1)	Def, H-Ca (m2 s−1)	Def, H-Mg (m2 s−1)	Def, H-Fe (m2 s−1)	
(3.78 ± 0.18) × 10−10	(1.62 ± 0.12) × 10−10	(1.51 ± 0.13) × 10−10	(2.74 ± 0.15) × 10−10	

Table 3 Effective diffusion coefficients determined by Eq. 2 for hydrogels based on pyrophosphate

Def, P (m2 s−1)	Def, P-Ca (m2 s−1)	Def, P-Mg (m2 s−1)	Def, P-Fe (m2 s−1)	
(3.45 ± 0.13) × 10−10	(1.30 ± 0.10) × 10−10	(1.18 ± 0.09) × 10−10	(2.60 ± 0.16) × 10−10	

A further effect is that the diffusion rate is dependent on the shape and real length of the diffusion pathway. Because the pores are not straight, the diffusion takes place more effectively over a longer distance than it would in a homogenous material. The tortuosity is a value characterizing the longer distance traversed in the pores. The structural parameter µ (= φ / τ) represents the influence of the structure of humic hydrogel and its local geometry in the diffusion (Klučáková 2022; Klučáková et al. 2014; Klučáková and Kalina 2015). If the hydrogel contains an active substance able to interact with diffusing particles, transport through the hydrogel is affected by these interactions and a proportion of the particles can be immobilized by the chemical reaction. The effective diffusion coefficient thus includes the effect of such interactions (besides the influence of the pore structure of the hydrogel). Interactions between humic acids and metal ions are relatively complex. Humic acids can be described as a mixture of many different organic substances containing different types of binding sites. A detailed description of the interactions is not possible. On the other hand, we can assume that the interactions are much faster in comparison with the diffusion. Therefore, a local equilibrium between immobilized and free movable Cu(II) ions can be attained. This means that the front of diffusing ions is in contact with the hydrogel containing non-occupied binding sites and (simultaneously) that this passing front is equilibrated from the point of view of the chemical reaction between humic acids and diffusing Cu(II) ions. If this simple equilibrium between immobilized and free movable Cu(II) ions is assumed, an apparent equilibrium constant K can be defined as the ratio between the concentrations of these two forms:3 K=cimcfree,

and included into the value of the effective diffusion coefficient (together with structural factor µ):4 Def=DφτK+1=DμK+1,

where D is the diffusion coefficient of Cu(II) ions in water.

The values of the effective diffusion coefficient obtained by means of the method of instantaneous planar source are listed in Tables 2 and 3. If we compare values obtained for hydrogels based on hydroxide (Table 2) and pyrophosphate (Table 3), we can see that the effective diffusion coefficients are lower for pyrophosphate cross-linked by different agents (HCl, CaCl2, MgCl2, FeCl3). However, the differences are in some cases comparable with the measuring error (e.g., H-Fe and P-Fe). The highest values of Def were obtained for pure hydrogels (without incorporated metals), the lowest ones for hydrogels containing incorporated magnesium. The effective diffusion coefficients obtained for hydrogels without incorporated metal ions are much higher than values obtained for hydrogels enriched by metals but lower than the diffusion coefficient of Cu(II) ions in water (Haynes 2012). This means that both the real length of the diffusion pathway in hydrogel pores and the interactions of Cu(II) ions with humic acids resulted in a decrease in diffusion coefficients.

Differences between Def values obtained for hydrogels based on hydroxide and pyrophosphate can be connected with structural changes to the humic acids. Francioso et al. (1998) reported that pyrophosphate could be incorporated in the structure of humic acids in their extraction process. Similarly, Tatzber et al. (2007) observed differences in the structural characteristics of humic acids extracted by pyrophosphate and other reagents. In contrast to the mentioned studies (Francioso et al. 1998; Tatzber et al. 2007), the humic acids used in this work were extracted by a mixture of NaOH and Na4P2O7 solutions. Then, the final powdered humic sample was used for the preparation of hydrogels. Differences in the properties of hydrogels based on hydroxide and pyrophosphate thus resulted from re-dissolving humic acids in one of the reagent solutions during hydrogel preparation.

As described above, differences in the FT-IR spectra of xerogels based on hydroxide and pyrophosphate were observed. The differences were similar to those published for the extraction of humic acids by means of various reagents (Francioso et al. 1998; Tatzber et al. 2007). This means that the structure and properties of humic acids can be changed by their re-dissolving. It is probable that the structural differences are connected with differences in the diffusion characteristics of hydrogels. The potential incorporation of pyrophosphate into the humic structure can result in changes in the permeability of the hydrogel as well as changes in its affinity to copper(II) ions. Both effects can cause a decrease in diffusion coefficients, and it is not easy to determine which effect predominates. If metal ions such as Ca(II), Mg(II), and Fe(II) were incorporated into hydrogels during their preparation, the decrease in the diffusion coefficient was much more significant. As described above, the presence of the mentioned metal ions can affect the binding of copper to organic matter (e.g., Gao et al. 2015; Hamilton-Taylor et al. 2002; Kogut and Voelker 2001; Muller and Cuscov 2017; Saito et al. 2005). Many authors have observed a decrease in the binding of copper in the presence of calcium and magnesium (Gao et al. 2015; Kang et al. 2010; Kogut and Voelker 2001). Wu and Hendershot (2010) confirmed the positive effect of calcium with respect to the suppression of copper toxicity.

In contrast, Hamilton-Taylor et al. (2002) showed that the presence of calcium and magnesium had no measurable effect. Our results show that the presence of calcium and magnesium suppressed the mobility of copper(II) ions in humic hydrogels. This indicates that calcium and magnesium can promote reactivity of humic hydrogels. It is known that the competition between copper and magnesium (and calcium) can decrease copper binding as a result of there being lower amounts of binding sites available for copper(II) ions (Gao et al. 2015; Kang et al. 2010). The effects observed in our experiments were different. The diffusion rate decreased as a result of chemical interactions and changes in the hydrogel structure. The most effective metal in our study was magnesium, closely followed by calcium. We assume that the metals are incorporated in the humic structure and occupy a proportion of humic binding sites. The stability of humic complexes with magnesium and calcium is lower in comparison with copper (Dinu 2015; Evangelou and Marsi 2001; Wieckowska and Martyniuk 2003). Therefore, metal ions can be replaced by copper(II) ions in the diffusion through humic hydrogel. Rey-Castro et al. (2009) derived the effective distribution of affinities of metal ions to humic substances under natural water conditions. They stated that calcium and magnesium display a greater effective affinity for carboxylic functional groups while the distributions of trivalent iron and copper between carboxylic and phenolic groups are overlapped. They computed that the affinity of Ca and Mg to carboxylic binding sites was higher than the affinity of Cu and Fe(III). These findings show that the presence of calcium and magnesium can suppress the interactions of copper with humic substances. In contrast, our results indicate that the diffusivity of Cu(II) ions significantly decreases if calcium and magnesium are incorporated into humic hydrogel. This means that their presence can suppress the permeability of the hydrogel and support its interactions with copper. The influence of the incorporation of Fe(II) ions was weaker but still significant. It can be connected with the different affinity to humic functional groups described above. In order to confirm or disprove the potential support of interactions between humic hydrogel and Cu(II) ions in the presence of other metals, the method of diffusion from a constant source was used for subsequent experiments.

Diffusion from a constant source

This method is based on a constant concentration at the interface between the hydrogel and the donor solution. In contrast to the previous method (diffusion from a coinstantaneous planar source), a higher amount of Cu(II) ions diffused into the hydrogel. Therefore, the potential dependence of the diffusion coefficient on concentration could be investigated. In this study, the source was a super-saturated solution, i.e., it was a saturated solution with a small amount of undissolved salt on the bottom of vessel. This salt was able to be dissolved when copper(II) ions diffused into the hydrogel, thus resulting in the constant concentration of the donor solution. The mathematical model applied to the data was based on the presumption of a semi-infinite medium. This means that the diffusion trajectories were shorter that the hydrogel dimension. The final solution of Fick’s second law is (Klučáková 2022; Klučáková et al. 2014; Klučáková and Kalina 2015):5 c=cserfcx2Deft,

where c is the concentration of diffusing particles, cs is the concentration at the interface between the hydrogel and the donor solution, S is the cross section available for diffusion, Def is the effective diffusion coefficient, t is time, and x is the distance from the interface. The total amount of Cu(II) ions diffused into the hydrogel (mt) can be expressed as (Klučáková 2022; Klučáková et al. 2014; Klučáková and Kalina 2015):6 mt=2csDeftπ

The effective diffusion coefficient Def can be calculated from the slope of the linear dependence of mt against the square root of time. In Fig. 4, the experimental example fitted by Eq. 6 is shown. The effective diffusion coefficients obtained by the method of constant source are listed in Tables 4 and 5. As can be seen, they are higher in comparison with the method of instantaneous planar source (Tables 2 and 3). This means that the diffusion coefficient is dependent on the concentration of diffusing particles. The diffusion rate is affected not only by the diffusion coefficient and also by the gradient of concentration. It means that lower diffusion coefficient and very high concentration gradient can result in higher diffusion rate than system with very low concentration gradient and high diffusion coefficient. Both parameters are higher in the case of the constant diffusion source; therefore, the diffusion rate should also be higher.Fig. 4 Example of data fitting for the method of constant source—Eq. 6: H hydrogel (red) and H-Ca hydrogel (blue)

Table 4 Effective diffusion coefficients determined by Eq. 6 for hydrogels based on hydroxide

Def, H (m2 s−1)	Def, H-Ca (m2 s−1)	Def, H-Mg (m2 s−1)	Def, H-Fe (m2 s−1)	
(9.91 ± 0.37) × 10−10	(6.89 ± 0.36) × 10−10	(5.63 ± 0.29) × 10−10	(1.68 ± 0.04) × 10−9	

Table 5 Effective diffusion coefficients determined by Eq. 6 for hydrogels based on pyrophosphate

Def, P (m2 s−1)	Def, P-Ca (m2 s−1)	Def, P-Mg (m2 s−1)	Def, P-Fe (m2 s−1)	
(9.61 ± 0.05) × 10−10	(6.53 ± 0.22) × 10−10	(6.17 ± 0.28) × 10−10	(1.67 ± 0.26) × 10−9	

As described above, the diffusion is strongly affected by the pore structure of the hydrogel and the interactions of diffusing particles with active substance incorporated in the hydrogel. In our study, the pore structure in the hydrogel was the same for both diffusion methods used and was determined by the type of hydrogel. In contrast, interactions between diffusing particles and active sites in the hydrogel can be affected by the concentration of reactants as well as the ionic strength. It is known that the rate of the chemical reaction is determined by the rate constant and the concentration of reactants (Monk 2004). An increase in concentration results in a higher reaction rate. The simple model of the local equilibrium between immobilized and free movable Cu(II) ions (Eq. 3 and Eq. 4) assumes that interactions are much faster than diffusion and that the system is (at a given place) equilibrated immediately. On the other hand, the equilibrium should be considered to be dynamic, because new diffusing particles are continuously passing through the given place. Therefore, the concentration of Cu(II) ions at a given place increases in time, as can be seen in Fig. 5, and local equilibrium is continually reestablished. It is known that the presence of other ions can affect reaction rates in aqueous medium. This phenomenon, called the primary salt effect, can cause both an increase and decrease in the reaction rate. The final effect depends on the charge of the reacting ions (Monk 2004). If they have the same charge, the reaction should be promoted and its rate increased.Fig. 5 Examples of concentration profiles in H (left) and H-Ca (right) hydrogels for the method of constant source. Data are fitted by Eq. 5: H hydrogel after 5 h (green) and 24 h (black), H-Ca hydrogel after 5 h (blue) and 24 h (red)

The effect of iron on the reaction kinetics in humic systems was also observed by other authors (Ding et al. 2019; Farina et al. 2018; Lippold et al. 2007; Shi et al. 2016). On the other hand, the primary salt effect can affect only the reaction rate. This means that the equilibrium should remain unaltered, but the system should be equilibrated in a shorter time.

However, the salt effect can influence the diffusivity of copper in hydrogel. The mechanisms of interactions between metal ions and humic substances are much more complex than a simple equilibrium between free movable and immobilized metal ions. There are many different binding sites in the structure of humic substances which differ in their chemical character and bond strength. The simplification which is represented by the assumption of the rapid achievement of local equilibrium can be considered as a useful tool for characterizing the transport properties of metal ions in humic systems. It cannot describe the complex process of diffusion connected with interactions in detail but can provide information about the mobility and amount of free movable metal ions. Figure 5 shows another effect of the incorporation of metal ions into hydrogel. We can see that the surface concentration of copper(II) ions differs for individual types of hydrogels. This confirms that the structure of hydrogels changed with the incorporation of metal ions into them. The concentration ratio at the phase interface between hydrogel and donor solution can thus be considered as resulting from the hydrogel structure and the ability of copper to interact with the hydrogel.

Conclusions

The migration behavior of metal ions in natural systems and their interactions with organic matter are crucial for determining their toxicity and bioavailability. We developed methods of mobility and reactivity mapping (Klučáková 2014; Klučáková et al. 2014; Klučáková and Kalina 2015; Klučáková and Pekař 2003a, b, 2004, 2009) which allow investigation of the transport of diffusing particles and their interactions with active substances directly with respect to their motion through the studied system. This study took into account the real situation in soils, which subsumes the existence of many organic and inorganic constituents. While previous studies were focused on the role of organic matter and humic substances in the migration and bioavailability of metal ions, the novelty of this study is the reactivity and transport mapping in systems containing complexed ions naturally occurring in soils (such as Ca, Mg, and Fe). The presence of Ca, Mg, and Fe in organo-mineral complexes can affect the biotoxicity of heavy metals (such as copper), whose accumulation in soil poses risks to human and ecosystem health (Christl et al. 2005; Wu and Hendershot 2010); for example, increased contents of calcium can suppress the uptake of copper by pea roots (Wu and Hendershot 2010). Our results confirm that the incorporation of Ca, Mg, and Fe into humic systems can strongly affect the mobility of copper in humic hydrogels. While the incorporation of Ca and Mg resulted in the decrease in the diffusivity of copper ions, the incorporation of Fe(III) into humic hydrogel resulted in an increase in the diffusivity of Cu(II) in the hydrogel in comparison with pure humic hydrogel. Diffusion coefficients based on the method of diffusion from constant source were higher in comparison with the instantaneous planar source which confirmed their concentration dependence. The use of pyrophosphate in the preparation of humic hydrogels resulted in moderate decrease in diffusivity of copper.

The effective diffusion coefficients determined by our methods include both the influence of the structure of the diffusion medium and the effects of interactions between diffusing particles and active substances in the medium. Our findings suggest that the soil type and its composition can affect the mobility of toxic pollutants and their uptake by plants. The current study has major implications for future research into the bioavailability of toxic metals in soils. It provides basic characteristics of the diffusivity of copper in humic hydrogels containing natural metals such as calcium, magnesium, and iron. Further investigations are required to determine how the mobility of copper (and other toxic metals) can be affected by the different contents of natural ions and organic matter.

Author contribution

Martina Klučáková: conceptualization, methodology, formal analysis, writing. Vojtěch Enev: conceptualization, methodology, data curation, investigation.

Funding

Open access publishing supported by the National Technical Library in Prague.

Data availability

Data will be made available on request.

Declarations

Ethics approval

Not applicable. Research does not involve humans and/or animals, and the authors declare compliance with ethical standards.

Consent to participate

All authors agreed with the content of article.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Bryan ND Barlow J Warwick P Stephens S Higg JJW Griffin D The simultaneous modelling of metal ion and humic substance transport in column experiments J Environ Monit 2005 7 196 202 10.1039/B413241F 15735778
Bryan ND, Barlow J, Warwick P, Stephens S, Higg JJW, Griffin D (2005) The simultaneous modelling of metal ion and humic substance transport in column experiments. J Environ Monit 7:196–202. 10.1039/B413241F15735778 10.1039/B413241F
Chen W Habibul N Liu X-Y Sheng G-P Yu H-Q FTIR and synchronous fluorescence heterospectral two-dimensional correlation analyses on the binding characteristics of copper onto dissolved organic matter Environ Sci Technol 2015 49 2052 2058 10.1021/es5049495 25611159
Chen W, Habibul N, Liu X-Y, Sheng G-P, Yu H-Q (2015) FTIR and synchronous fluorescence heterospectral two-dimensional correlation analyses on the binding characteristics of copper onto dissolved organic matter. Environ Sci Technol 49:2052–2058. 10.1021/es504949525611159 10.1021/es5049495
Christl I Milne CJ Kinniburgh DG Kretzschmar R Relating ion binding by fulvic and humic acids to chemical composition and molecular size. 2 Metal Binding Environ Sci Technol 2001 35 2512 2517 10.1021/es0002520 11432556
Christl I, Milne CJ, Kinniburgh DG, Kretzschmar R (2001) Relating ion binding by fulvic and humic acids to chemical composition and molecular size. 2. Metal Binding Environ Sci Technol 35:2512–2517. 10.1021/es000252011432556 10.1021/es0002520
Christl I Metzger A Heidmann I Kretzschmar R Effect of humic and fulvic acid concentrations and ionic strength on copper and lead binding Environ Sci Technol 2005 39 5319 5326 10.1021/es050018f 16082962
Christl I, Metzger A, Heidmann I, Kretzschmar R (2005) Effect of humic and fulvic acid concentrations and ionic strength on copper and lead binding. Environ Sci Technol 39:5319–5326. 10.1021/es050018f16082962 10.1021/es050018f
Conte P Piccolo A Conformational arrangement of dissolved humic substances. Influence of solution composition on association of humic molecules Environ Sci Technol 1999 33 1682 1690 10.1021/es9808604
Conte P, Piccolo A (1999) Conformational arrangement of dissolved humic substances. Influence of solution composition on association of humic molecules. Environ Sci Technol 33:1682–1690. 10.1021/es980860410.1021/es9808604
De Nobili M Bragato G Alcaniz JM Puigbo A Comellas L Characterization of electrophoretic fractions of humic substances with different electrofocusing behavior Soil Sci 1990 150 763 770 10.1097/00010694-199011000-00002
De Nobili M, Bragato G, Alcaniz JM, Puigbo A, Comellas L (1990) Characterization of electrophoretic fractions of humic substances with different electrofocusing behavior. Soil Sci 150:763–770. 10.1097/00010694-199011000-0000210.1097/00010694-199011000-00002
Ding Y Liu M Peng S Li J Liang Y Shi Z Binding characteristics of heavy metals to humic acid before and after fractionation by ferrihydrite Chemosphere 2019 226 140 148 10.1016/j.chemosphere.2019.03.124 30925406
Ding Y, Liu M, Peng S, Li J, Liang Y, Shi Z (2019) Binding characteristics of heavy metals to humic acid before and after fractionation by ferrihydrite. Chemosphere 226:140–148. 10.1016/j.chemosphere.2019.03.12430925406 10.1016/j.chemosphere.2019.03.124
Dinu MI (2015) Interaction between metal ions in waters with humic acids in gley-podzolic soils. Geochem Int 53”265–276. 10.1134/S0016702915030052
Evangelou VP Marsi M Composition and metal ion complexation behaviour of humic fractions derived from corn tissue Plant Soil 2001 229 13 24 10.1023/A:1004862100925
Evangelou VP, Marsi M (2001) Composition and metal ion complexation behaviour of humic fractions derived from corn tissue. Plant Soil 229:13–24. 10.1023/A:100486210092510.1023/A:1004862100925
Farina AO Peacock CL Fiol S Antelo J Carvin B A universal adsorption behaviour for Cu uptake by iron (hydr)oxide organo-mineral composites Chem Geol 2018 479 23 35 10.1016/j.chemgeo.2017.12.022
Farina AO, Peacock CL, Fiol S, Antelo J, Carvin B (2018) A universal adsorption behaviour for Cu uptake by iron (hydr)oxide organo-mineral composites. Chem Geol 479:23–35. 10.1016/j.chemgeo.2017.12.02210.1016/j.chemgeo.2017.12.022
Francioso O Ciavatta C Tugnoli V Sanchez-Cortes S Gessa C Spectroscopic characterization of pyrophosphate incorporation during extraction of peat humic acids Soil Sci Soc Am J 1998 62 181 187 10.2136/sssaj1998.03615995006200010024x
Francioso O, Ciavatta C, Tugnoli V, Sanchez-Cortes S, Gessa C (1998) Spectroscopic characterization of pyrophosphate incorporation during extraction of peat humic acids. Soil Sci Soc Am J 62:181–187. 10.2136/sssaj1998.03615995006200010024x10.2136/sssaj1998.03615995006200010024x
Gao Y Yan M Korshin G Effects of calcium on the chromophores of dissolved organic matter and their interactions with copper Water Res 2015 81 47 53 10.1016/j.watres.2015.05.038 26043370
Gao Y, Yan M, Korshin G (2015) Effects of calcium on the chromophores of dissolved organic matter and their interactions with copper. Water Res 81:47–53. 10.1016/j.watres.2015.05.03826043370 10.1016/j.watres.2015.05.038
Gerzabek M Pichlmayer F Kirchmann H Haberhauer G The response of soil organic matter to manure amendments in a long-term experiment at Ultuna, Sweden Eur J Soil Sci 1997 48 273 282 10.1111/j.1365-2389.1997.tb00547.x
Gerzabek M, Pichlmayer F, Kirchmann H, Haberhauer G (1997) The response of soil organic matter to manure amendments in a long-term experiment at Ultuna, Sweden. Eur J Soil Sci 48:273–282. 10.1111/j.1365-2389.1997.tb00547.x10.1111/j.1365-2389.1997.tb00547.x
Haberhauer G Gerzabek MH Drift and transmission FT-IR spectroscopy of forest soils: an approach to determine decomposition processes of forest litter Vibr Spectr 1999 19 413 417 10.1016/S0924-2031(98)00046-0
Haberhauer G, Gerzabek MH (1999) Drift and transmission FT-IR spectroscopy of forest soils: an approach to determine decomposition processes of forest litter. Vibr Spectr 19:413–417. 10.1016/S0924-2031(98)00046-010.1016/S0924-2031(98)00046-0
Hamilton-Taylor J Postill AS Tipping E Harper MP Laboratory measurements and modeling of metal-humic interactions under estuarine conditions Geochim Cosmochim Acta 2002 66 403 415 10.1016/S0016-7037(01)00777-3
Hamilton-Taylor J, Postill AS, Tipping E, Harper MP (2002) Laboratory measurements and modeling of metal-humic interactions under estuarine conditions. Geochim Cosmochim Acta 66:403–415. 10.1016/S0016-7037(01)00777-310.1016/S0016-7037(01)00777-3
Hayes TM Hayes MHB Skjemstad JO Swift RS Compositional relationships between organic matter in a grassland soil and its drainage waters Eur J Soil Sci 2008 59 603 616 10.1111/j.1365-2389.2007.01007.x
Hayes TM, Hayes MHB, Skjemstad JO, Swift RS (2008) Compositional relationships between organic matter in a grassland soil and its drainage waters. Eur J Soil Sci 59:603–616. 10.1111/j.1365-2389.2007.01007.x10.1111/j.1365-2389.2007.01007.x
Hayes TM Hayes MH Swift RS Detailed investigation of organic matter components in extracts and drainage waters from a soil under long term cultivation Org Geochem 2012 52 13 22 10.1016/j.orggeochem.2012.07.011
Hayes TM, Hayes MH, Swift RS (2012) Detailed investigation of organic matter components in extracts and drainage waters from a soil under long term cultivation. Org Geochem 52:13–22. 10.1016/j.orggeochem.2012.07.01110.1016/j.orggeochem.2012.07.011
Haynes WM Handbook of chemistry and physics 2012 93 Boca Raton CRC Press
Haynes WM (2012) Handbook of chemistry and physics, 93rd edn. CRC Press, Boca Raton
Head MJ Zhou WJ Evaluation of NaOH leaching techniques to extract humic acids from paleosols Nucl Instr Meth Phys Res B 2000 172 434 439 10.1016/S0168-583X(00)00221-4
Head MJ, Zhou WJ (2000) Evaluation of NaOH leaching techniques to extract humic acids from paleosols. Nucl Instr Meth Phys Res B 172:434–439. 10.1016/S0168-583X(00)00221-410.1016/S0168-583X(00)00221-4
Kang F Hamilton PB Long J Wang Q Influence of calcium precipitation on copper sorption induced by loosely bound extracellular polymeric substance (LB-EPS) from activated sludge Fundam Appl Limnol Arch Hydrobiol 2010 176 173 181 10.1127/1863-9135/2010/0176-0173
Kang F, Hamilton PB, Long J, Wang Q (2010) Influence of calcium precipitation on copper sorption induced by loosely bound extracellular polymeric substance (LB-EPS) from activated sludge. Fundam Appl Limnol Arch Hydrobiol 176:173–181. 10.1127/1863-9135/2010/0176-017310.1127/1863-9135/2010/0176-0173
Klučáková M Complexation of metal ions with solid humic acids, humic colloidal solutions, and humic hydrogel Environ Eng Sci 2014 31 612 620 10.1089/ees.2013.0487
Klučáková M (2014) Complexation of metal ions with solid humic acids, humic colloidal solutions, and humic hydrogel. Environ Eng Sci 31:612–620. 10.1089/ees.2013.048710.1089/ees.2013.0487
Klučáková M The effect of supramolecular humic acids on the diffusivity of metal ions in agarose hydrogel Molecules 2022 27 1019 10.3390/molecules27031019 35164280
Klučáková M (2022) The effect of supramolecular humic acids on the diffusivity of metal ions in agarose hydrogel. Molecules 27:1019. 10.3390/molecules2703101935164280 10.3390/molecules27031019
Klučáková M Kalina M Diffusivity of Cu(II) ions in humic gels - influence of reactive functional groups of humic acids Colloid Surface A 2015 483 162 170 10.1016/j.colsurfa.2015.05.041
Klučáková M, Kalina M (2015) Diffusivity of Cu(II) ions in humic gels - influence of reactive functional groups of humic acids. Colloid Surface A 483:162–170. 10.1016/j.colsurfa.2015.05.04110.1016/j.colsurfa.2015.05.041
Klučáková M Pekař M Study of structure and properties of humic and fulvic acids. III. Study of complexation of Cu2+ ions with humic acid in sols J Polym Mater 2003 20 145 154
Klučáková M, Pekař M (2003a) Study of structure and properties of humic and fulvic acids. III. Study of complexation of Cu2+ ions with humic acid in sols. J Polym Mater 20:145–154
Klučáková M Pekař M Study of structure and properties of humic and fulvic acids. IV. Study of interactions of Cu2+ ions with humic gels and final comparison J Polym Mater 2003 20 155 162
Klučáková M, Pekař M (2003b) Study of structure and properties of humic and fulvic acids. IV. Study of interactions of Cu2+ ions with humic gels and final comparison. J Polym Mater 20:155–162
Klučáková M Pekař M Ghabbour EA Davies G Study of diffusion of metal cations in humic gels Humic substances: Nature`s most versatile materials 2004 New York Taylor & Francis 263 273
Klučáková M, Pekař M (2004) Study of diffusion of metal cations in humic gels. In: Ghabbour EA, Davies G (eds) Humic substances: Nature`s most versatile materials. Taylor & Francis, New York, pp 263–273
Klučáková M Pekař M Transport of copper(II) ions in humic gel – new results from diffusion couple Colloid Surface A 2009 349 96 101 10.1016/j.colsurfa.2009.08.001
Klučáková M, Pekař M (2009) Transport of copper(II) ions in humic gel – new results from diffusion couple. Colloid Surface A 349:96–101. 10.1016/j.colsurfa.2009.08.00110.1016/j.colsurfa.2009.08.001
Klučáková M Kalina M Sedláček P Grasset L Reactivity and transport mapping of Cu(II) ions in humic hydrogels J Soil Sediment 2014 14 368 376 10.1007/s11368-013-0730-2
Klučáková M, Kalina M, Sedláček P, Grasset L (2014) Reactivity and transport mapping of Cu(II) ions in humic hydrogels. J Soil Sediment 14:368–376. 10.1007/s11368-013-0730-210.1007/s11368-013-0730-2
Kogut M Voelker B Strong copper-binding behavior of terrestrial humic substances in seawater Environ Sci Technol 2001 35 1149 1156 10.1021/es0014584 11347927
Kogut M, Voelker B (2001) Strong copper-binding behavior of terrestrial humic substances in seawater. Environ Sci Technol 35:1149–1156. 10.1021/es001458411347927 10.1021/es0014584
Lippold H Evans NDM Warwick P Kupsch H Competitive effect of iron(III) on metal complexation by humic substances: characterisation of ageing processes Chemosphere 2007 67 1050 1056 10.1016/j.chemosphere.2006.10.045 17140629
Lippold H, Evans NDM, Warwick P, Kupsch H (2007) Competitive effect of iron(III) on metal complexation by humic substances: characterisation of ageing processes. Chemosphere 67:1050–1056. 10.1016/j.chemosphere.2006.10.04517140629 10.1016/j.chemosphere.2006.10.045
Maia CMBF Piccolo A Mangrich AS Molecular size distribution of compost-derived humates as a function of concentration and different counterions Chemosphere 2008 73 1162 1166 10.1016/j.chemosphere.2008.07.069 18778847
Maia CMBF, Piccolo A, Mangrich AS (2008) Molecular size distribution of compost-derived humates as a function of concentration and different counterions. Chemosphere 73:1162–1166. 10.1016/j.chemosphere.2008.07.06918778847 10.1016/j.chemosphere.2008.07.069
Martyniuk H Więckowska J Adsorption of metal ions on humic acids extracted from brown coals Fuel Process Technol 2003 84 23 36 10.1016/S0378-3820(02)00246-1
Martyniuk H, Więckowska J (2003) Adsorption of metal ions on humic acids extracted from brown coals. Fuel Process Technol 84:23–36. 10.1016/S0378-3820(02)00246-110.1016/S0378-3820(02)00246-1
Monk P (2004) Physical chemistry. Understanding our chemical world, first ed, John Wiley & Sons Ltd, Chichester
Muller FLL Cuscov M Alteration of the copper-binding capacity of iron-rich humic colloids during transport from peatland to marine waters Environ Sci Technol 2017 51 3214 3222 10.1021/acs.est.6b05303 28218520
Muller FLL, Cuscov M (2017) Alteration of the copper-binding capacity of iron-rich humic colloids during transport from peatland to marine waters. Environ Sci Technol 51:3214–3222. 10.1021/acs.est.6b0530328218520 10.1021/acs.est.6b05303
Novotny EH Blum WEH Gerzabek MH Mangrich AS Soil management system effects on size fractionated humic substances Geoderma 1999 92 87 109 10.1016/S0016-7061(99)00022-1
Novotny EH, Blum WEH, Gerzabek MH, Mangrich AS (1999) Soil management system effects on size fractionated humic substances. Geoderma 92:87–109. 10.1016/S0016-7061(99)00022-110.1016/S0016-7061(99)00022-1
Nuzzo A Sanchez A Fontaine B Piccolo A Conformational changes of dissolved humic and fulvic superstructures with progressive iron complexation J Geochem Explor 2013 129 1 5 10.1016/j.gexplo.2013.01.010
Nuzzo A, Sanchez A, Fontaine B, Piccolo A (2013) Conformational changes of dissolved humic and fulvic superstructures with progressive iron complexation. J Geochem Explor 129:1–5. 10.1016/j.gexplo.2013.01.01010.1016/j.gexplo.2013.01.010
Paradelo M Perez-Rodriguez P Fernandez-Calvino D Arias-Estevez M Lopez-Periago JE Coupled transport of humic acids and copper through saturated porous media Eur J Soil Sci 2012 63 708 716 10.1111/j.1365-2389.2012.01481.x
Paradelo M, Perez-Rodriguez P, Fernandez-Calvino D, Arias-Estevez M, Lopez-Periago JE (2012) Coupled transport of humic acids and copper through saturated porous media. Eur J Soil Sci 63:708–716. 10.1111/j.1365-2389.2012.01481.x10.1111/j.1365-2389.2012.01481.x
Peuravuori J Žbánková P Pihlaja K Aspects of structural features in lignite and lignite humic acids Fuel Process Technol 2006 87 829 839 10.1016/j.fuproc.2006.05.003
Peuravuori J, Žbánková P, Pihlaja K (2006) Aspects of structural features in lignite and lignite humic acids. Fuel Process Technol 87:829–839. 10.1016/j.fuproc.2006.05.00310.1016/j.fuproc.2006.05.003
Peuravuori J Simpson AJ Lam B Žbánková P Pihlaja K Structural features of lignite humic acid in light of NMR and thermal degradation experiments J Mol Struct 2007 826 131 142 10.1016/j.molstruc.2006.04.048
Peuravuori J, Simpson AJ, Lam B, Žbánková P, Pihlaja K (2007) Structural features of lignite humic acid in light of NMR and thermal degradation experiments. J Mol Struct 826:131–142. 10.1016/j.molstruc.2006.04.04810.1016/j.molstruc.2006.04.048
Rey-Castro C Mongin S Huidobro C David C Salvador J Garcés JL Galceran J Mas F Puy J Effective affinity distribution for the binding of metal ions to a generic fulvic acid in natural waters Environ Sci Technol 2009 43 7184 7191 10.1021/es803006p 19848120
Rey-Castro C, Mongin S, Huidobro C, David C, Salvador J, Garcés JL, Galceran J, Mas F, Puy J (2009) Effective affinity distribution for the binding of metal ions to a generic fulvic acid in natural waters. Environ Sci Technol 43:7184–7191. 10.1021/es803006p19848120 10.1021/es803006p
Rosa HA Oliveira LC Bellin IC Rocha JC Romão LPC Filho NLD Influence of alkaline extraction on the characteristics of humic substances in Brasilian soils Thermochim Acta 2005 433 81 86 10.1016/j.tca.2005.02.001
Rosa HA, Oliveira LC, Bellin IC, Rocha JC, Romão LPC, Filho NLD (2005) Influence of alkaline extraction on the characteristics of humic substances in Brasilian soils. Thermochim Acta 433:81–86. 10.1016/j.tca.2005.02.00110.1016/j.tca.2005.02.001
Saito T Koopal LK Nagasaki S Tanaka S Analysis of copper binding in the ternary system Cu2+/humic acid/goethite at neutral to acidic pH Environ Sci Technol 2005 39 4886 4893 10.1021/es0500308 16053088
Saito T, Koopal LK, Nagasaki S, Tanaka S (2005) Analysis of copper binding in the ternary system Cu2+/humic acid/goethite at neutral to acidic pH. Environ Sci Technol 39:4886–4893. 10.1021/es050030816053088 10.1021/es0500308
Scally S Davison W Zhang H Diffusion coefficients of metals and metal complexes in hydrogels used in diffusive gradients in thin films Anal Chim Acta 2006 558 222 229 10.1016/j.aca.2005.11.020
Scally S, Davison W, Zhang H (2006) Diffusion coefficients of metals and metal complexes in hydrogels used in diffusive gradients in thin films. Anal Chim Acta 558:222–229. 10.1016/j.aca.2005.11.02010.1016/j.aca.2005.11.020
Shi Z Wang P Peng L Lin Z Dang Z Kinetics of heavy metal dissociation from natural organic matter: roles of the carboxylic and phenolic sites Environ Sci Technol 2016 50 10476 10484 10.1021/acs.est.6b01809 27580051
Shi Z, Wang P, Peng L, Lin Z, Dang Z (2016) Kinetics of heavy metal dissociation from natural organic matter: roles of the carboxylic and phenolic sites. Environ Sci Technol 50:10476–10484. 10.1021/acs.est.6b0180927580051 10.1021/acs.est.6b01809
Sowers TD Adhikari D Wang J Yang Y Sparks DL Spatial associations and chemical composition of organic carbon sequestered in Fe, Ca, and organic carbon ternary systems Environ Sci Technol 2018 52 6936 6944 10.1021/acs.est.8b01158 29770687
Sowers TD, Adhikari D, Wang J, Yang Y, Sparks DL (2018) Spatial associations and chemical composition of organic carbon sequestered in Fe, Ca, and organic carbon ternary systems. Environ Sci Technol 52:6936–6944. 10.1021/acs.est.8b0115829770687 10.1021/acs.est.8b01158
Sureshkumar MK Das D Mary G Nuwad J Adsorption of Pb(II) ions using humic acid coated chitosan-tripolyphosphate (HA-CTPP) beads Sep Sci Technol 2013 48 1132 1139 10.1080/01496395.2012.724500
Sureshkumar MK, Das D, Mary G, Nuwad J (2013) Adsorption of Pb(II) ions using humic acid coated chitosan-tripolyphosphate (HA-CTPP) beads. Sep Sci Technol 48:1132–1139. 10.1080/01496395.2012.72450010.1080/01496395.2012.724500
Tatzber M Stemmer M Spiegel H Katzlberger C Haberhauer G Mentler A Gerzabek MH FTIR-spectroscopic characterization of humic acids and humin fractions obtained by advanced NaOH, Na4P2O7, and Na2CO3 extraction procedures J Plant Nutr Soil Sci 2007 170 522 529 10.1002/jpln.200622082
Tatzber M, Stemmer M, Spiegel H, Katzlberger C, Haberhauer G, Mentler A, Gerzabek MH (2007) FTIR-spectroscopic characterization of humic acids and humin fractions obtained by advanced NaOH, Na4P2O7, and Na2CO3 extraction procedures. J Plant Nutr Soil Sci 170:522–529. 10.1002/jpln.20062208210.1002/jpln.200622082
Thurman EM Malcolm RL Preparative isolation of aquatic humic substances Environ Sci Technol 1981 15 463 466 10.1021/es00086a012 22248415
Thurman EM, Malcolm RL (1981) Preparative isolation of aquatic humic substances. Environ Sci Technol 15:463–466. 10.1021/es00086a01222248415 10.1021/es00086a012
Tipping E Vincent CD Lawlor AJ Lofts S Metal accumulation by stream bryophytes, related to chemical speciation Environ Pollut 2008 156 936 943 10.1016/j.envpol.2008.05.010 18541353
Tipping E, Vincent CD, Lawlor AJ, Lofts S (2008) Metal accumulation by stream bryophytes, related to chemical speciation. Environ Pollut 156:936–943. 10.1016/j.envpol.2008.05.01018541353 10.1016/j.envpol.2008.05.010
van Oort F Monna F Garnier S Zn/Pb concentration ratios emphasize spatiotemporal airborne metal dynamics in soils under different land use Water Air Soil Pollut 2020 231 109 10.1007/s11270-020-04478-1
van Oort F, Monna F, Garnier S (2020) Zn/Pb concentration ratios emphasize spatiotemporal airborne metal dynamics in soils under different land use. Water Air Soil Pollut 231:109. 10.1007/s11270-020-04478-110.1007/s11270-020-04478-1
Watanabe Kuwatsuka S Fractionation of soil fulvic acids using polyvinyl-pyrrolidone and their ionization difference spectra Soil Sci Plant Nutr 1991 37 611 617 10.1080/00380768.1991.10416929
Watanabe, Kuwatsuka S (1991) Fractionation of soil fulvic acids using polyvinyl-pyrrolidone and their ionization difference spectra. Soil Sci Plant Nutr 37:611–617. 10.1080/00380768.1991.1041692910.1080/00380768.1991.10416929
Wieckowska J Martyniuk H Adsorption of metal ions on humic acids extracted from brown coals Fuel Process Technol 2003 84 23 36 10.1016/S0378-3820(02)00246-1
Wieckowska J, Martyniuk H (2003) Adsorption of metal ions on humic acids extracted from brown coals. Fuel Process Technol 84:23–36. 10.1016/S0378-3820(02)00246-110.1016/S0378-3820(02)00246-1
Wu Y Hendershot WH Effect of calcium and pH on copper binding and rhizotoxicity to pea (Pisum sativum L.) root: Empirical relationships and modeling Arch Environ Contam Toxicol 2010 59 109 119 10.1007/s00244-009-9450-4 20052462
Wu Y, Hendershot WH (2010) Effect of calcium and pH on copper binding and rhizotoxicity to pea (Pisum sativum L.) root: Empirical relationships and modeling. Arch Environ Contam Toxicol 59:109–119. 10.1007/s00244-009-9450-420052462 10.1007/s00244-009-9450-4
Zhang H Davison W Diffusional characteristics of hydrogels used in DGT and DET techniques Anal Chim Acta 1999 398 329 340 10.1016/S0003-2670(99)00458-4
Zhang H, Davison W (1999) Diffusional characteristics of hydrogels used in DGT and DET techniques. Anal Chim Acta 398:329–340. 10.1016/S0003-2670(99)00458-410.1016/S0003-2670(99)00458-4
