Mobility of Zn and Cu in Bentonites: Implications for Environmental Remediation

The aim of this study was to evaluate the mobility of copper (Cu) and zinc (Zn) and their impact on the properties of bentonites and unfrozen water content. Limited research in this area necessitates further analysis to prevent the negative effects of metal interactions on bentonite effectiveness. Tests involved American (SWy-3, Stx-1b) and Slovak (BSvk) bentonite samples with Zn or Cu ion exchange. Sequential extraction was performed using the Community Bureau of Reference (BCR) method. Elemental content was analyzed via inductively coupled plasma optical emission spectrometry (ICP-OES). Unfrozen water content was measured using nuclear magnetic resonance (1H-NMR) and differential scanning calorimetry (DSC). Results showed a significant influence of the main cation (Zn or Cu) on ion mobility, with toxic metal concentrations increasing mobility and decreasing residual fractions. Mobile Zn fractions increased with larger particle diameters, lower clay content, and shorter interplanar spacing, while the opposite was observed for Cu. Zn likely accumulated in larger clay pores, while Cu was immobilized in the bentonite complex. The stability of Zn or Cu ions increased with higher clay content or specific surface area. Residual Zn or Cu fractions were highest in uncontaminated bentonites with higher unfrozen water content, suggesting the potential formation of concentrated solutions in sub-zero temperatures, posing a threat to the clay–water environment, especially in cold regions.


Introduction
In recent years, there has been a growing interest in the topic of toxic metals, their mobility in the soil-water environment, methods of their immobilization, and removal from wastewater [1][2][3][4][5][6].Metals accumulating in soil indirectly affect human and animal health through the food chain, dust, and contamination of surface and groundwater, potentially leading to serious damage [6][7][8].Zinc (Zn) is an essential trace element supporting the immune system, but its excess can lead to immune system disorders and copper (Cu) deficiency; Cu aids enzymatic functions, but in excess, it causes liver and kidney damage [6,9].The use of fertilizers containing high concentrations of Cu and Zn, industrial activities, and flotation tailings storage are the most common sources of these metals [10][11][12].Both Cu and Zn are trace elements that play vital roles as essential nutrients for plants.However, their excessive accumulation in the soil can lead to harmful effects.For instance, elevated levels of these metals can decrease plants' nutrient uptake capacity and reduce the population of soil microorganisms.Consequently, this can disrupt vital soil biological processes, including organic matter decomposition and the nitrogen cycle [12][13][14].The heavy metals nickel (Ni), lead (Pb), and chromium (Cr) pose a significant concern for water and soil environments.Their leaching from municipal or industrial waste can result in soil contamination and toxicity to organisms inhabiting it.Moreover, Cr and Pb can also undergo bioaccumulation in living organisms [15].Liao et al. [16] reported that soil pollution by Ni and Cu is an extremely serious problem.Once in the soil, these metals move within the soil and are not degraded by microorganisms.They are susceptible to transformation into various chemical forms such as dissolved ions, adsorbed forms, exchangeable forms, and complexes with organic matter through processes such as dissolution, adsorption, coagulation, and complexation [17].
Interactions of potentially toxic metals with soils are a complex and dynamic process that is still not fully understood.Feng et al. [18] and Li et al. [2] believe that significant influences on these processes are interactions of metals with microorganisms, organic matter, and clay minerals.The main components of clay minerals are silicate minerals and oxides, which affect the binding and migration of metals in soils [19].Clay minerals, such as montmorillonite, which is the main component of bentonites, are often used to improve soil quality for various applications due to their high specific surface area, ion exchange capacity, and strong stability.One of these applications is their use as a mineral barrier in waste landfills, including hazardous ones [3,[20][21][22].Bentonites in contact with various chemical substances, including ions of potentially toxic metals, are not inert substances and may undergo degradation in the interactions due to clay dissolution, cation exchange between layers, or ion competition [3,23].As a result, the engineering properties might change, which were initially meant to keep the facility and the environment safe [24].In cold regions, the effects of toxic metal interactions in bentonites can be more hazardous.The remaining liquid phase in frozen soils, known as unfrozen water, forms a concentrated solution, whose ions actively interact with the surface ions of minerals, often leading to changes in its volume [25,26].Therefore, assessing the mobility of potentially toxic metals in the clay-water environment, especially in bentonites, is an important engineering and environmental issue.Understanding the interactions of ions of potentially toxic metals in bentonites can help explain their migration mechanisms, which may minimize the risk of their spreading into the soil-water system through appropriate protective measures in the future.
To study the dynamics of changes occurring in the soil-water system contaminated with ions of potentially toxic metals, it is important to determine not only the total concentration of metals present in the soil, but also their chemical form, by which they become mobile.Sequential extraction, which separates metals into pools with different tendencies to migrate in the soil medium, should be the basis of studies in this field [2,27,28].One of the generally accepted and most widely used sequential extraction methods, next to the Tessier method, is the BCR technique [28,29].This method consists of four steps, each using chemical extraction solutions of increasing chemical aggressiveness.This allows the division of metals into fractions (F): F I , exchangeable and associated with carbonates; F II , associated with Fe and Mn oxides; F III , associated with organic matter and sulfides; and F IV , associated with primary and secondary minerals.Metals in fractions I and II are generally considered to pose the greatest threat to the aquatic soil environment, whereas those in fraction IV are considered unavailable [28].A literature review has established that certain factors can influence the mobility of toxic metals in clayey soils, many of which are directly related to changes in the soil's physicochemical properties due to chemical interactions with the solution.However, these mechanisms are highly complex, and the situation is challenging because metal ions interact with substances dissolved in the pore solution during migration, which can affect the efficiency of metal migration [30].Various factors influence the mobility of toxic metals in clayey soils, including the anthropogenic origin of metals, bentonite quality, metal concentration, temperature, presence of organic matter, proportion of clay minerals, and pH [24,[31][32][33].Anthropogenic contamination of soils with cadmium (Cd), Zn, and Pb increases their mobility [24].Dutta et al. [34] and Nartowska et al. [35] suggest that at low concentrations of metal salts in bentonites, bentonite quality affects clay permeability, while at higher concentrations, it becomes the dominant factor.He et al. [36] propose that, in low-permeability soils, solute transport results from the interaction of adsorption and diffusion forces, with solution molarity playing a significant role.For instance, increased Zn retention was observed with increasing ionic strength of the solution.Nartowska et al. [26] observed lower unfrozen water content in bentonites contaminated with Zn or Cu ions compared to uncontaminated bentonites.They also identified statistically significant relationships between unfrozen water content and the clay fraction in bentonites contaminated with Zn or Cu.However, under positive temperature conditions, the behavior of bentonites contaminated with Cu ions is more influenced by changes in their structural parameters, including specific surface area [35,37], while in bentonites contaminated with Zn ions, it is related to changes in their granulometric parameters [24].Huang et al. [38] found that the size and characteristics of soil aggregates significantly affect the mobility and bioavailability of heavy metals in soils.Some authors suggest that the content of Cu, Pb, and Zn increases with decreasing soil particle size, while others attribute this relationship to the type of metal [2].The presence of clay minerals and organic matter generally promotes metal immobilization and reduces their mobility [39].The potential of hydrogen (pH) is cited as one of the most significant parameters affecting the total concentration of metals in soils [33,40] and is crucial for assessing metal mobility [41,42].However, the observed variability in metal mobility cannot be solely explained by changes in pH [31].Variations in metal adsorption capacity can be attributed to acidic modification of clay surface properties [43].In turn, potentially toxic metals can exist as particles adsorbed onto clay surfaces, highlighting their undeniable relationship with the physicochemical properties of the soil-water system [44,45].
The characteristics of the migration and transformation of contaminants in altered soils require further study [30].There is a lack of research in the literature regarding the following issues: (i) The mobility of potentially toxic metals in homogeneous (model) clays devoid of organic matter, or in bentonites, was evaluated.Studies on source clays will facilitate the interpretation of the processes responsible for the mobility of potentially toxic metal ions so that the results can be more easily related to natural conditions.In addition, it may not be possible to identify the mechanisms governing sorption in systems with many metals in the clay [46]; (ii) The relationships between the mobility of potentially toxic metals and the physicochemical properties of soils can be interpreted.The authors emphasized that this is especially important in clay soils, where a significant portion of the metals may be bound to the clay surface.Li et al. [2] emphasized that the interactions between potentially toxic metal ions and clay minerals require further research.Ma et al. [19] emphasized that understanding the mechanisms and suitability of different types of clay minerals in immobilizing toxic metals is crucial to achieving remediation goals; (iii) The mobility of toxic metals, in the context of their physicochemical properties at subzero temperatures, was evaluated.Understanding the interactions of these processes under freezing conditions is important because the concentrations of mineral salts and other substances are highest in very-cold winters [45].Because water can be a carrier for the migration of potentially toxic metals in cold regions, it is important to know the relationship between these factors to take appropriate steps to minimize the risk.It is particularly important to understand the relationships between the chemical fractions of metals that determine their mobility and changes in unfrozen water content.To date, no such studies have been conducted.
Thus, the main objective of this study was to evaluate the mobility of potentially toxic metals, Zn or Cu, in bentonites using the BCR procedure in the temperature range 20 • C to −32 • C. Specific objectives included the following: (i) determination of the total concentration of potentially toxic metals (Zn, Cu, Cr, Ni, and Pb) and their contribution to chemical fractions determined via the BCR procedure in uncontaminated bentonites and in bentonites contaminated with Zn or Cu ions of different origins; (ii) determination of the relationship between the total concentration of Zn or Cu, microstructural properties, physicochemical properties of bentonites, and the proportion of chemical fractions of these metals in bentonites contaminated with Zn or Cu ions; and (iii) determination of the relationship between Zn or Cu chemical fractions, microstructural and physicochemical properties, and the content of unfrozen water in bentonites contaminated with Zn or Cu ions at temperatures below 0 • C. Achieving the objectives allowed for a comprehensive analysis of the impact of contaminating bentonites with Zn or Cu ions on their physicochemical properties, which is applicable in both environmental protection and materials engineering.The urgency of the research hypothesis stems from the need to understand the mobility of these metals in bentonites under various environmental conditions, including low temperatures, which is crucial for assessing environmental hazards in cold climate regions and for effective waste landfill management.The main innovations of the present work are as follows: (i) the use of model bentonites, whose main exchangeable cation is Zn or Cu, to analyze the mobility of potentially toxic metals.Because bentonites are widely used worldwide as adsorbents for toxic metals, knowledge of their mobility in these clays is of particular importance; (ii) exploration of the relationship between the chemical fractions of metals and properties of bentonites.This knowledge is also useful when planning remedial works at waste landfills.An example of an engineering solution is the study by Koda et al. [20][21][22], who designed a cut-off bentonite barrier for groundwater protection at the sanitary landfill in Radiowo, Poland; and (iii) capturing the relationships between the physicochemical parameters of bentonites and the unfrozen water content at subzero temperatures.This knowledge is particularly important in cold regions, and the global literature on this subject is still limited.

Materials
Three model clays were used to assess the mobility of potentially toxic metals in bentonites: SWy-3 (Na-Ca bentonite from Crook County, WY, USA), STx-1b (Ca bentonite from Gonzales County, TX, USA), and BSvk (Ca bentonite from Jelsovy Potok, Slovakia).The American bentonites were obtained from the Source Clays Repository of the Clay Mineral Society, while the Slovak bentonite was obtained from ZGM Zębiec S. A., Starachowice, Poland.Each of the three model bentonites was tested in its natural form and previously prepared homoionic forms (Zn or Cu).Nine homoionic forms of bentonites were obtained, designated as follows: SWy-3 1M Cu, STx-1b 1M Cu, and BSvk 1M Cu, copper forms of bentonites formed via saturation with 1 mol/dm 3 copper (II) chloride solutions; BSvk 0.1M Cu, BSvk 0.25M Cu, and BSvk 0.5M Cu, copper forms of Slovak bentonite formed via saturation with copper (II) chloride solutions (Chempur, Poland, CAS No. 10125-13-0) of 0.1-0.5 mol/dm 3 ; and SWy-3 1M Zn, STx-1b 1M Zn, and BSvk 1M Zn, zinc forms of bentonites formed via saturation with zinc (II) chloride (Chempur, Piekary Śl ąskie, Poland, CAS No. 7646-85-7) solutions of 1 mol/dm 3 .A significant disparity in the content of mobile Cu fractions in bentonites of different origins was observed, leading to a study of the effect of Cu ion concentration on the proportions of copper BCR fractions in selected BSvk bentonite at various concentrations.
In the BCR studies, the following reagents (Chempur, Poland) with their respective CAS numbers were used:  The average perimeter of micropores [μm] 7.17 The pore circularity 0.283 0.315 0.245 Each bentonite contained at least 70% montmorillonite in its mineral composition, in addition to silica minerals and biotite, and in the case of bentonites contaminated with Cu 0.1-0.5 M ions, 3-6% of atacamite (copper hydroxychloride) minerals were observed [47].Source bentonites, according to the diffractograms (Figure 2), contained the following: 70% montmorillonite, 30% quartz (SWy-3); 75% montmorillonite, 20% opal-CT, 5% quartz (STx-1b); and 92% smectite, 5% quartz and 3% biotite (BSvk) [47].The results are nearly consistent with the data presented in the literature [48,49].Each bentonite contained at least 70% montmorillonite in its mineral composition, in addition to silica minerals and biotite, and in the case of bentonites contaminated with Cu 0.1-0.5 M ions, 3-6% of atacamite (copper hydroxychloride) minerals were observed [47].Source bentonites, according to the diffractograms (Figure 2), contained the following: 70% montmorillonite, 30% quartz (SWy-3); 75% montmorillonite, 20% opal-CT, 5% quartz (STx-1b); and 92% smectite, 5% quartz and 3% biotite (BSvk) [47].The results are nearly consistent with the data presented in the literature [48,49].(c) SWy-3 bentonite Preparation of Homoionic Forms Bentonite (50 g) was flooded in buckets with 10 liters of zinc (II) chloride solutions of 1 mol/dm 3 and copper (II) chloride solutions at concentrations of 0.1-0.25-0.5-1mol/dm 3 to obtain homoionic forms of clay.Flooding was performed at intervals of 48 hours, after which the solution was decanted from the sediment and the flooding process was repeated with a new portion of the solution.The operation was repeated three times.Subsequently, the remaining suspension was transferred to chloride ion-permeable membranes, which were immersed in a bucket with a constant circulation of distilled water.The water was changed every 48 hours for approximately 35 days.Excess chlorine was flushed out until the characteristic reaction in the water to Cl − ions (from 0.1 mol/dm 3 AgNO3, Chempur, Poland) disappeared.The monoionic bentonite samples were air-dried and sequentially transferred to sealed jars.The details of the methodology are included in the authors' previous work [26].In addition, the absence of free chloride ions in all 12 clays was confirmed with analytical methods.Solutions consisting of 2 g clay and 20 ml of distilled water were shaken for 24 hours on a rotary shaker at room temperature.The samples were successively filtered through Whatman blotting paper, grade 1.6 μm, 50 mm diameter circle (Maidstone, Kent, UK).The prepared solutions were titrated with silver nitrate in the presence of chromate as an indicator (PN-ISO 9297:1994 [50]).

Determination of the Content of Potentially Toxic Metals
A sequential extraction procedure was used to fractionate the metals according to their mobility, providing insight into the ionic dynamics of potentially toxic metals in bentonites.This is crucial because the hazard of a metal depends on its specific form.Understanding the toxic and chemical forms of metals is often more important than knowing their total content, because their bioavailability is closely determined by their mobility [28,51].Zimmermann and Weindrof [29] confirmed that the total content of elements in soils does not correspond to the processes in which they participate in the Preparation of Homoionic Forms Bentonite (50 g) was flooded in buckets with 10 liters of zinc (II) chloride solutions of 1 mol/dm 3 and copper (II) chloride solutions at concentrations of 0.1-0.25-0.5-1mol/dm 3 to obtain homoionic forms of clay.Flooding was performed at intervals of 48 h, after which the solution was decanted from the sediment and the flooding process was repeated with a new portion of the solution.The operation was repeated three times.Subsequently, the remaining suspension was transferred to chloride ion-permeable membranes, which were immersed in a bucket with a constant circulation of distilled water.The water was changed every 48 h for approximately 35 days.Excess chlorine was flushed out until the characteristic reaction in the water to Cl − ions (from 0.1 mol/dm 3 AgNO 3 , Chempur, Poland) disappeared.The monoionic bentonite samples were air-dried and sequentially transferred to sealed jars.The details of the methodology are included in the authors' previous work [26].In addition, the absence of free chloride ions in all 12 clays was confirmed with analytical methods.Solutions consisting of 2 g clay and 20 ml of distilled water were shaken for 24 h on a rotary shaker at room temperature.The samples were successively filtered through Whatman blotting paper, grade 1.6 µm, 50 mm diameter circle (Maidstone, Kent, UK).The prepared solutions were titrated with silver nitrate in the presence of chromate as an indicator (PN-ISO 9297:1994 [50]).

Determination of the Content of Potentially Toxic Metals
A sequential extraction procedure was used to fractionate the metals according to their mobility, providing insight into the ionic dynamics of potentially toxic metals in bentonites.This is crucial because the hazard of a metal depends on its specific form.Understanding the toxic and chemical forms of metals is often more important than knowing their total content, because their bioavailability is closely determined by their mobility [28,51].Zimmermann and Weindrof [29] confirmed that the total content of elements in soils does not correspond to the processes in which they participate in the environment, and that only studies using sequential extraction can provide a basis for investigating the dynamics of changes in the soil-water environment.
This study used a four-step BCR sequential extraction with a modification introduced by European researchers [52].This method is considered one of the widely accepted methods in soil chemistry, the results of which are easy to reproduce [29].The residual fraction was mineralized using aqua regia.Potentially toxic metals were separated into individual mobility fractions that showed different tendencies to migrate in the clay medium.To determine the subsequent fractions, the test clay was dissolved in a reagent specific to that fraction [53].The mixtures were then centrifuged to separate the solid phase from the solution that contained the extracted metals.The contents of potentially toxic metals (Cu, Zn, Cr, Ni, and Pb) in the extracts were determined using a Perkin Elmer Optima 8000 inductively coupled plasma ICP-OES optical emission spectrometer (PerkinElmer, Waltham, MA, USA).The BCR procedure involves the following steps (Table 1).

F IV -Metals bound to primary and secondary minerals
In all the bentonites tested, the total concentration of metals was additionally determined as an independent test, which is useful to assess the recovery of metals from individual BCR fractions and to establish the relationship between the total concentration of Zn or Cu and the proportion of their mobile and immobile fractions.The procedure was the same as that in step 4. The results of these analyses are summarized in Tables A1 and A2 in the Appendix A.

Determination of the Physicochemical Properties
The granulometric composition was analyzed using a HELOS/BF SUCELL laser diffractometer from Sympatec GmbH (Clausthal-Zellerfeld, Germany), equipped with a 'wet' test attachment, according to the ISO 13320:2020 standard [54] and Sympatec guidelines [55].The clay samples were prepared 48 h prior to testing by processing into clay pastes.Immediately prior to testing, 3 g of clay paste were mixed with 50 ml of distilled water.A sample of a few milliliters was taken for analysis to ensure that the optical concentration of particles in the test chamber did not exceed 25%.A reference measurement was made before each new measurement [32].The diameters d 10 and d 32 were obtained from the generated granulometric composition curves.These values denote grain diameters, the mass proportion of which, together with smaller grains in a given clay, is 10% and 32%, respectively.Effective diameters are used in many soil filtration models, so they may prove useful in analyzing the dynamics of potentially toxic metals in porous clays such as bentonites [24].Additionally, the d 32 diameter, known as the Sauter diameter, is a commonly used measure in fluid dynamics as the mean surface volume diameter of a particle.It is thought to characterize many important processes, including chemical processes [56].
The specific surface area was determined, which is useful in the analysis of clay-water systems, using the water vapor sorption (WST) method [26].The specific surface area was determined according to Equation ( 1): where w 50 is the relative water vapor pressure (p/p 0 = 0.50) in a desiccator over a saturated solution of magnesium nitrate (V) for 10 days.
The mineral composition and the interpacket distances in the bentonites, in the plane parallel to the sample surface (d 001 ), were determined based on X-ray diffractograms.Changes in inter-packet distances serve as a useful tool for assessing the metal dynamics, ion exchange, and adsorption properties of bentonites [57].For this study, a Bruker D8 Advance X-ray diffractometer with a Johansson-type monochromator was employed to determine the CuKα1 radiation (λ = 1.5406Å) with a LynxEye position-sensitive detector (Bruker devices, Berlin, Germany).The measurements were conducted in the 2θ range from 4.51 • to 70 • with 0.02 • steps.The applied voltage was 3.540 kV with a current of 530 mA [57].

Determination of the Changes in Unfrozen Water Content
In soils, there is a fraction of liquid water that does not freeze even at very-low subzero temperatures (e.g., up to 20% of this fraction can be expected in bentonites at T = −30 • C).This is referred to as unfrozen water content.Its amount depends on the temperature, type of soil, physical and chemical parameters, and soil chemistry.
For the analyses, the authors used the results of changes in the unfrozen water content of bentonites published in their previous works [26,47].The mean unfrozen water content and standard deviation are summarized in Table A3.To study the changes in unfrozen water content in the Zn-bentonites, the technique of 1 H NMR (nuclear magnetic resonance) was employed.This study was conducted with a minispec mq20 spectrometer (Bruker, Berlin, Germany).The content of unfrozen water at particular freezing temperatures was determined by changes in the intensity of the free induction decay (FID) signal, which is directly proportional to the number of 1 H hydrogen nuclei and thus to the volume of water and ice in the clay sample.Unfrozen water was calculated according to Equation (2).The highest FID signal value for a given sample above 0 • C was determined according to Nartowska et al. [58]: where w u is the unfrozen water content (% of clay dry weight), SI FID is the signal intensity of the free inductive decay at a given temperature (in this study, the highest signal value for a given sample above 0 • C was determined [58]), and w n is the natural water content of the clay (% clay dry matter).
Based on preliminary studies of changes in unfrozen water content in Cu bentonites using the 1 H NMR technique [26,58], it was suspected that the determination results might be affected by ion hydration, which requires further research [58].Therefore, DSC was used to evaluate the changes in unfrozen water content in the Cu bentonites.The tests were performed using a DSC Q200 instrument (New Castle, DE, USA, TA Instruments) with a liquid nitrogen cooling system (RSC 90).The samples were cooled at a scan rate of 2.5 K/min to −90 • C, followed by heating at a scan rate of 5 K/min to 20 • C. Due to the possible phenomenon of overcooling, the heating thermograms were analyzed.The unfrozen water content at a given temperature was calculated using Equation (3): where w u (T i ) is the unfrozen water content at temperature T i as a percentage of dry mass, w is the water content expressed as a percentage of dry mass, m s is the mass of dry clay in the sample (g), and L(T j ) is the latent heat of fusion of ice at temperature T j , calculated according to the empirical equation L(T) = 7.3T + 334.

Data Quality Control
The methods were based on the ISO/IEC 17025 guidelines [59].Source clays with one main exchangeable cation, which determines the physicochemical properties of bentonites, were used in this study.Bentonites are materials with a grain size of d < 63 µm, which eliminates the possible influence of granulometric composition on the mobility of metals due to high homogeneity [51].The metal speciation analysis was carried out step by step according to the strict BCR procedure developed by the European Commission within the framework of the "Certified Reference Materials for the Quality Control of Chemical Speciation" project, improved by European researchers [29].Reagents free of contaminants that could affect the analytical results (pure for analysis or chemically pure, Sigma Aldrich) were used.The validation of the BCR method involved various steps.To reduce errors, four determinations were made for each bentonite sample, and the reagent blanks were established.Coefficients of variation (CV) of <5% ensured the accuracy and reproducibility of the method.Recovery control was performed by comparing the total metal content (determined in a separate test) with the total metal content in all fractions.The efficiency was 95-105% and met the quality control standards.Metals were determined using inductively coupled plasma optical emission spectrometry (ICP-OES).Before the test, the instrument was calibrated using a multielement standard for ICP-OES (Instrument Calibration Standard 2, No. N9301721, PerkinElmer).At this stage, the linearity of the calibration curve, representing the relationship between the optical signal and the element concentration, was determined.The fit of the calibration curve was assessed using a linear regression statistical model with a significance level (p = 0.01).The background signal level was monitored during the experiments.The results were statistically analyzed to exclude gross errors using the Grubbs' test.

Statistical Treatments
The mobility of potentially toxic metals in bentonites contaminated with Zn or Cu ions was analyzed using Statistica 12 software (StatSoft Inc., Tulsa, OK, USA), and the interpretation of the results was based on the guidelines provided by Rabiej [60].Normality of the distribution was checked using the Kolmogorov-Smirnov test, while homogeneity of variance was assessed using Levene's test.To assess how the independent variables, the origin of the clay and its contamination, affect the treated set of dependent variables, which were the fractions of potentially toxic metals (F I + F II ) Cu, Zn, Cr, Ni, and Pb, MANOVA analysis was used.In addition, interaction effects were tested between the bentonite origin and contamination to test more complex hypotheses.The discriminatory ability of the model was assessed using Wilks' lambda values.Scheffe's conservative post-hoc test was performed to determine the significance of differences between group averages of 'origin of land'.Regression analysis was used to assess whether changes in the BCR Zn or Cu chemical fractions could be predicted from the independent variables (total Zn or Cu concentration and changes in unfrozen water content).The value of the coefficient of determination R 2 describes the percentage of the linear model describing the variation in a specific metal chemical fraction (F I -F IV ).The standard error of the estimate makes it possible to assess the deviation of the actual concentration of a given metal chemical fraction from the predicted value.The unstandardized coefficient (B) represents the values of the model coefficients.The standardized coefficient (ß) compares the strength of the effect of the independent variable on the dependent variable.The higher the absolute value of the beta coefficient, the stronger the effect.A test probability value of p < 0.05 allows the conclusion that the independent variables are statistically significantly related to the dependent variable.

The Properties of Bentonites
The physicochemical parameters of bentonites are summarized in Table 2.All potentially toxic elements are listed in Tables A1 and A2 in the Appendix A.

Effect of Bentonite Origin and Main Exchangeable Cation Cu or Zn on Mobility of Potentially Toxic Metals
The percentages of individual BCR Cu or Zn fractions (F I -Exchangeable and weak acid soluble fractions; F II -Reducible fraction, F III -Oxidizable fraction, and F IV -Residual fraction) and the direction of their changes differed depending on the origin of the clay and the type of the main exchangeable cation.The proportion of the fraction of potentially mobile metals (F I + F II ) in each bentonite after contamination with Cu or Zn ions differed from those in the uncontaminated clays (Tables A1 and A2).MANOVA showed a statistically significant effect of clay type, contamination with Zn or Cu ions, and the interaction of these factors on the entire set of potentially mobile metal fractions (F I + F II ) of Cu, Zn, Cr, Ni, and Pb (Table 3).Thus, the results do not exclude the possibility that the mobility of metals in model clays depends on both the origin of the clay and the type of the main exchangeable cation (Figure 1).Studies by Nartowska [24] and Nartowska et al. [57] confirmed the significant influence of the origin of the clay and the type of main exchangeable cation on shaping the physicochemical and structural properties of bentonites.Changes in these properties can result in different metal mobilities.
Significantly higher mobility of metals was observed in the Slovakian BSvk bentonite than in the American STx-1b and SWy-3 bentonites.It is most likely that the quality of bentonite is important here (Figure 1).For example, Dutta et al. [34] confirmed that bentonite quality affects its hydraulic properties.This, in turn, can be related to the transport of solutes into the clay pores.Scheffe's post hoc test showed that the average content of potentially mobile metals varied depending on the origin of the bentonite, with the most significant differences occurring in BSvk bentonite contaminated with Cu ions (p = 0.001).Therefore, this clay was selected for further analysis.The highest potential mobility of these metals was achieved in bentonites contaminated with Cu or Zn ions (Figure 3).Contamination of the clay with a particular metal increases mobility.Kashem et al. [31] confirmed that the anthropogenic origins of Cd, Zn, and Pb affect the increase in the proportion of their mobile fractions in contaminated soils in Japan.In the present study, a five-fold increase in the mobility rates of Cd and Zn and a two-fold increase in the mobility rate of Pb in clays contaminated with these metals were observed compared to uncontaminated clays.Similarly, Ma et al. [19] observed the highest concentrations of anthropogenic Cr and Cd in the mobile fraction.In this fraction, metals such as Pb, Co, Zn, and Ni were not observed, suggesting that these ions may form stable complexes with Fe and Mn oxides, and Co, Ni, and Cu may additionally be complex with humic substances.This phenomenon requires further investigation.Figure 4 shows that the mobility of metals in bentonites contaminated with Zn and Cu ions differs.The lowest mobilities of Pb and Ni were observed in the clays contaminated with Zn ions.Therefore, further analyses were conducted separately for clays contaminated with Cu and Zn ions.

The Effect of the Total Concentration of Zn or Cu on the Proportion of These Metals in Mobile and Non-Mobile Fractions
The proportion of mobile and stable fractions of Zn and Cu varied depending on the type of bentonite as well as the metal concentration (Figure 4).The Zn mobility in uncontaminated bentonites followed the pattern F IV > F II > F III > F I (BSvk) or F IV > F III > F II > F I (Stx-1b, SWy-3).The stable fraction dominated in all clays, confirming the suitability of bentonites for accumulating potentially toxic metals.These results are consistent with those reported in the literature [61,62].After the bentonite was contaminated with Zn ions, the proportion of individual chemical fractions changed significantly.The dominant directions were F III > F I > F II > F IV (BSvk, Stx-1b) or F I > F III > F II > F IV .In any bentonite contaminated with Zn ions, the smallest share was the immobile fraction, F IV , and the mobile fraction, F I , had a high share.The smallest proportion of non-mobile fractions after Zn ion contamination was observed in Wyoming bentonite, originally sodium bentonite.
Copper mobility in uncontaminated bentonites followed the pattern F IV > F II > F III > F I and for samples contaminated with Cu ions, F I > F II > F III > F IV or F II > F III > F I > F IV .In each case, the proportion of stable fractions decreased significantly after the samples were contaminated with Cu ions.Similarly, Liu et al. [63] observed that the availability of Cu and Cd was negatively correlated with the percentage of the residual fraction of metals.
The results of the regression analysis confirmed a statistically significant relationship between the total Zn ion content in the clay and its content in the individual Zn fractions determined with the BCR method (F I -F IV ) (Table 4).
An increase in Zn concentration had a statistically significant (p < 0.05) effect on the increase in mobile fractions (F I , F II ) and a decrease in stable fractions (F III , F IV ), regardless of the origin of the bentonite.In view of the marked disparity in the content of mobile Cu fractions in bentonites of different origins (Figure 4B) and the significant Scheffe test, a study of the effect of Cu ion concentration on the proportion of copper BCR fractions was carried out on BSvk bentonite.Statistically significant relationships were obtained only for the mobile fractions (F I and F II ).With an increase in the concentration of Cu ions, the proportion of fraction I decreased, indicating a different behavior from that of bentonites contaminated with Zn ions.Most likely, in clays contaminated with Cu ions, the higher the concentration of Cu, the greater the proportion of ions trapped in the structure of bentonites, while simultaneously reducing their mobility under exchangeable fraction conditions.This was confirmed in an earlier study by Nartowska and Kozłowski [47], which showed that, in bentonites contaminated with Cu ions, the presence of a new mineral phase, atacamite, increased in proportion with increasing Cu ion concentration.The possible formation of metal compounds, resulting in a reduction in metal mobility, was also reported by Zhang et al. [64].According to the authors, higher soil pH may promote the adsorption of positively charged heavy metal ions by soil colloids and may increase ion binding with OH − groups.In addition to assessing the mobility of the main exchangeable ions Zn and Cu, certain dependencies in the mobility of co-occurring ions in the bentonite were observed.In both Zn-and Cu-ion-contaminated bentonites, relationships were observed between the concentration of metals co-occurring in these clays, such as Cu, Cr, Ni, and Pb, or Cr, Ni, Pb, and Zn, and their mobility.Moreover, their behavior in both clays was the opposite.Similarly, Nartowska [24] showed different behaviors of bentonites contaminated with Zn and Cu ions.Wang et al. [65] confirmed that the presence of co-occurring ions in the structure of bentonites is one of the factors, in addition to specific surface area, pore size, organic matter content, pH, and temperature, that affect the rate of ion release from the clay structure, in this case, phosphates.

The Effect of the Total Concentration of Zn or Cu on the Proportion of These Metals in Mobile and Non-Mobile Fractions
The proportion of mobile and stable fractions of Zn and Cu varied depending on the type of bentonite as well as the metal concentration (Figure 4).The Zn mobility in uncontaminated bentonites followed the pattern FIV > FII > FIII > FI (BSvk) or FIV > FIII > FII > FI (Stx-1b, SWy-3).The stable fraction dominated in all clays, confirming the suitability of bentonites for accumulating potentially toxic metals.These results are consistent with those reported in the literature [61,62].After the bentonite was contaminated with Zn ions, the proportion of individual chemical fractions changed significantly.The dominant directions were FIII > FI > FII > FIV (BSvk, Stx-1b) or FI > FIII > FII > FIV.In any bentonite contaminated with Zn ions, the smallest share was the immobile fraction, FIV, and the mobile fraction, FI, had a high share.The smallest proportion of non-mobile fractions after In the case of bentonites contaminated with Zn ions, the following statistically significant relationships (p < 0.05) were confirmed: an increase in Cu concentration with an increase in Cu mobility index, an increase in Cr concentration with an increase in Cr mobility index and the content of potentially mobile Cr fractions (F I + F II ), an increase in Ni and Pb concentration affected the decrease in potentially mobile fractions of these metals (F I + F II ), and an increase in Cu, Cr, and Ni concentration affected the increase in the Pb mobility index.In the case of bentonites contaminated with Cu ions, statistically significant relationships (p < 0.05) were confirmed: an increase in Cr content decreased the mobility index, and an increase in Ni content increased it.An increase in Pb content is related to an increase in the potentially mobile fraction (F I + F II ) of this metal.The total concentrations of Cr, Ni, and Zn affected the decrease in the mobility of Pb.The observed correlations, as they involve ions co-occurring in the soil structure with theoretically low concentrations of these metals, they do not pose a threat to the environment; but may provide a basis for further research using bentonites in which Cr, Ni, and Pb ions predominate.

The Relationship between the Physicochemical Properties of Bentonites and the Mobility of Potentially Toxic Metals
Using regression analysis, it was assessed whether changes in physicochemical properties of bentonites, such as clay and silt fraction content, effective diameter (d 10 ), Sauter mean diameter (d 32 ), and interplane distance d 001 , are correlated with the content of individual metals (Zn, Cu, Cr, Ni, and Pb) in potentially mobile (F I + F II ) and mobile (F I ) fractions.
The analysis showed statistically significant relationships (p < 0.05) in bentonites contaminated with Zn ions between the content of potentially mobile Cr, Zn, and Ni fractions and the clay fraction content (R Cr,Zn = 0.99, R Ni = −0.99).The content of potentially mobile Cu fractions increased with an increasing silt fraction.The content of potentially mobile Pb fractions increased with a decrease in the specific surface area of the clay and a decrease in the silt fraction.Therefore, Pb may accumulate in coarser fractions.The higher the clay fraction content, the lower the proportion of mobile fractions (F I ) of Cu, Cr, Pb, and Zn.For Ni, the proportion of mobile fractions increased with increasing Ni concentration.Metal mobility (F I ) is also strongly influenced by the granulometric parameters d 10 and d 32 .The proportion of mobile fractions of Cu, Cr, Pb, and Zn increased with an increase in these parameters and a decrease in Ni.One could assume that the leachability of these metals is related to the ions trapped in the clay pores.Generally, the larger the effective diameters, the larger the pores.The higher mobility of Cu, Cr, Pb, and Zn metals with a decrease in d 001 may suggest that ions located in larger pores are released at the F I stage.This occurs because a decrease in d 001 leads to a decrease in microporosity [57].Conversely, this is the case for Ni.Similarly, Angelaki et al. [66] observed an increase in filtration rate, particle aggregation, and porosity in clay soils contaminated with Zn or Cu ions.
The relationship between the Cu mobile fraction (F I ) and clay fraction, d 32 , and interpacket distance was significant, but in the opposite direction from that of bentonites contaminated with Zn ions.Changes in the behavior of metals that co-occurred in the clay structure (Cr, Ni, and Pb) were also observed.The proportion of the Cr mobile fraction (F I ) increased with decreasing specific surface area and d 001 .The stability index of Ni (Ir) increased with an increase in the specific diameter and a decrease in the d 001 .The Pb mobility index increased with increasing effective diameter and decreasing d 001 , specific surface area, and clay fraction content.
Due to the different behaviors of metals in bentonites contaminated with Zn or Cu ions, it seems reasonable to learn about the influence of individual metals in bentonites where they occur as the main exchangeable cation.The data presented herein can provide a basis for further analysis.

The Impact of the Predominant Exchangeable Cation in Bentonites, Whether Zinc or Copper
Because the physicochemical properties of bentonites are determined by the type of main exchangeable cation [57], the most significant relationships between physicochemical parameters and the mobility of potentially toxic metals can be expected by evaluating the influence of Zn or Cu ions.
It was established that the mobility of Zn or Cu is related to the changes in the physicochemical properties of the bentonites resulting from their introduction into the structure of bentonites due to ion exchange (Figure 5).Statistically significant correlations (p < 0.05) were established in bentonites contaminated with Zn ions (Figure 5A) between F IZn and F IIZn and clay fraction content, d 001 , d 10 , and d 32 .Significant (p < 0.05) correlations in bentonites contaminated with Cu ions (Figure 5B) were obtained for F ICu and F IICu with clay fraction content, d 001 and d 32 , and for F IICu additionally with d 10 .For F IVCu with specific surface area, the correlations were not significant.
Figure 6 shows the relationship between the different chemical fractions and physicochemical properties, depending on the type of metal, Zn or Cu.The results indicate different behavior of these bentonites and clarify the results shown in Individual stages in the BCR method correspond to metal ions present in different chemical forms in the soil-water environment: Stage 1, metals in exchangeable form or bound to carbonates (F I ); Stage 2, metals bound to amorphous iron and manganese oxides (F II ); Stage 3, metals bound to organic matter or sulfides (F III ); and Stage 4, metals bound to primary and secondary minerals (F IV ).Metals released in stages 1 and 2 are considered to pose a threat to the environment due to their high mobility, whereas metals released in stages 3 and 4 are considered to have no significant impact on the environment due to their high stability.
parameters and the mobility of potentially toxic metals can be expected by evaluating the influence of Zn or Cu ions.
It was established that the mobility of Zn or Cu is related to the changes in the physicochemical properties of the bentonites resulting from their introduction into the structure of bentonites due to ion exchange (Figure 5).Statistically significant correlations (p < 0.05) were established in bentonites contaminated with Zn ions (Figure 5A) between FIZn and FIIZn and clay fraction content, d001, d10, and d32.Significant (p < 0.05) correlations in bentonites contaminated with Cu ions (Figure 5B) were obtained for FICu and FIICu with clay fraction content, d001 and d32, and for FIICu additionally with d10.For FIVCu with specific surface area, the correlations were not significant.Based on the correlation coefficients assessment (Figure 5) and the diagram (Figure 6), the likely mechanism of Zn or Cu mobility in bentonites depends on clay properties, and the type of chemical fraction determined via the BCR method was established.
For Zn bentonites, it was observed that in stage 1 of the BCR, Zn ions accumulate in larger pores where free water is present.These observations are confirmed by the decrease in interpacket distances, which is closely related to the decrease in microporosity [57].In turn, a decrease in microporosity indicates an increase in the proportion of pores with larger diameters.At stage 2 of the BCR method, higher clay fraction contents and shorter d 001 distances resulted in reduced mobility of Zn ions.This suggests that microporosity may have further decreased, as indicated by the lower correlation coefficient for d 001 compared to F I .The increase in Zn ions with increasing effective diameters indicates that Zn is still released from free water located in larger pores.At stage 3, it was observed that the stability of Zn ions increases with increasing specific surface area, increasing silt fraction content, and increasing interpacket distances.Most likely, Zn accumulates in the interpacket space and on the mineral surface.The association with the silt fraction may result from particle aggregation after contamination with Zn ions.Zn ions are likely to be present in bound water (interlayer water and diffuse double-layer water).At stage 4, the stability of Zn ions increases with increasing clay fraction content and increasing interpacket distances.Larger interpacket distances probably lead to a higher proportion of micropores [57].Zn ions accumulating in micropores are bound to bentonite particles and may be difficult to access for ion exchange.Copper fractions in bentonites contaminated with copper ions.Bentonite parameters tested: Clay-d < 0.002 mm, Siltd (0.002-0.063 mm), d10-effective diameter is that size below which 10% of the material is contained, SMD-Sauter Mean Diameter is also referred to as d32 (Laser diffraction method); d001interpacket distance (XRD method); S-specific surface (Water sorption test method).
Figure 6 shows the relationship between the different chemical fractions and physicochemical properties, depending on the type of metal, Zn or Cu.The results indicate different behavior of these bentonites and clarify the results shown in

Zn
The proportion of the Zn mobile fractions (F I ) increases with an increase in effective diameter (d 10 , SMD) and a decrease in clay fraction and interpacket distances.
The proportion of the Zn mobile fractions (F II ) increases with an increase in d 10 , SMD, and a decrease in clay fraction and interpacket distances.
The content of the Zn stable fractions (F III ) increases with an increase in specific surface area, the increase in silt fraction content, and the increase in interpacket distances.
The content of the Zn stable fractions (F IV ) increases with an increase in the interpacket distances and an increase in the clay fraction.The role of specific surface area decreases.

Cu
The proportion of the Cu mobile fractions (F I ) increases with an increase in the clay fraction content and an increase in the interpacket distances.
The proportion of the Cu mobile fractions (F II ) increases with an increase in d 10 , SMD, and a decrease in clay fraction and interpacket distances.
The content of the Cu stable fractions (F III ) showed no significant statistical relationship with physicochemical parameters.
The content of the Cu stable fractions (F IV ) increases with an increase in the specific surface area.For Cu bentonites, at stage 1 of the BCR, it was observed that the proportion of the Cu mobile fraction increases with the increase in clay fraction content and interpacket distances.Cu ions are likely released from the interpacket spaces of the bentonite.The larger the space between packets, the more Cu ions can be contained and released more easily.The highest clay fraction content and the largest interpacket distances were observed in the bentonite 'BSvk 1M Cu'.Lower values were found in the bentonites 'BSvk 0.1-0.5 M Cu'.It can be assumed that the high proportion of the mobile copper fraction is associated with the mineralogical composition of the bentonites.Bentonites 'BSvk 0.1-0.5 M Cu' contain a few percent of copper-binding minerals (atacamites) in their structure [47]; hence, in such bentonites, the mobility of Cu may be impended.At stage 2 of the BCR, Cu ions are released from water located in larger pores of the bentonite, as evidenced by the increase in effective diameters with the increase in the amount of released ions.Stage 3 of the BCR did not show statistically significant relationships between the proportion of potentially stable Cu fractions and the physicochemical parameters of the bentonite.The content of the stable Cu fraction at stage 4 of the BCR increased with the increase in the specific surface area of the bentonites.Conversely, the specific surface area increased with the increase in the ionic strength of the copper chloride (II) solution used for ion exchange (0.1M-1M).It is likely that a large portion of Cu ions is located on the surface of bentonite particles and is bound to them by attractive forces, which increase with higher molar concentrations of the solution.Hao et al. [43] confirmed that at low ionic strength of the solution, the removal capacity of Cu from sodium bentonite is the highest.
The results obtained were partially consistent with those reported by Ma et al. [19].The authors studied river sediments from mining areas in China and concluded that the Cu, Zn, Cd, and As residual fractions dominated in sediments with grain sizes larger than 150 µm, while the other three BCR fractions inversely dominated in fractions smaller than 150 µm.Although the research carried out in this article deals with bentonites with grain sizes < 63 µm, it can be observed that the contribution of fractions I-III of Zn or Cu is the largest and the residual fraction is the smallest.However, the results of this study suggest that the relationships between granulometric composition and BCR chemical fractions are more complex, as described above, and depend strongly on the type of Zn or Cu metal and the type of BCR fraction.A different mechanism of Zn or Cu mobility was also observed by Alloway et al. [67], who noted that the mobility of Zn is higher than that of Cu because Cu is strongly retained by the soil solid phase.Similarly, a strong association between Cu ions and the specific surface of the soil is evident in the presented results.

Mobility of Zn or Cu at Freezing Temperatures
Bentonites subjected to Zn or Cu ion contamination, as well as additionally exposed to freezing-thawing cycles, constitute a complex clay-water system.Barker et al. [68] suggest that the mobilization of trace metals is seasonally dependent.At subzero temperatures, there exists a fraction of liquid water acting as a carrier for potentially toxic metals, and its quantity may be associated with the interaction of the ions themselves and changes in the properties of bentonites due to their introduction to the structure of bentonite [25,26].Rui et al. [69], on the other hand, argue that freezing-thawing can lead to the disaggregation of soil particles, resulting in a higher proportion of Cd and Pb metals in the thawing front, although they emphasize the significance of soil particle size and specific surface area.Consequently, the correlational relationships between the unfrozen water content at specific negative temperatures with BCR Zn or Cu chemical fractions and clay properties were assessed (Figure 7).
Additionally, to statistically confirm significant relationships between individual BCR Zn or Cu fractions for bentonites contaminated with Zn or Cu ions and unfrozen water content at specific negative temperatures, regression analysis was conducted (Table 5).
The results of changes in unfrozen water content at selected temperatures are presented in Table A3.
Figure 7 shows the best correlations of unfrozen water content at a given freezing temperature with chemical fractions determined using the BCR method: F IZn , F IIZn , F IVZn , and F IVCu , clay, S, d 001 , d 10 , and d 32 .
The proportion of exchangeable and reducible fractions of Zn (F IZn and F IIZn ) correlated significantly and in the same direction as the d 10 and d 32 in bentonites contaminated with Zn ions (Figure 7A).An increase in unfrozen water content was observed with a decrease in F IZn content, F IIZn content at a given negative temperature, effective diameter, and Sauter diameter.The decrease in Zn concentration in the fractions (F IZn and F IIZn ) of Zn bentonites is probably a consequence of a decrease in the d 10 due to contamination of bentonites with Zn ions.The increase in microporosity causes some Zn ions to be trapped in small pores, making them more difficult to elute.
suggest that the mobilization of trace metals is seasonally dependent.At subzero temperatures, there exists a fraction of liquid water acting as a carrier for potentially toxic metals, and its quantity may be associated with the interaction of the ions themselves and changes in the properties of bentonites due to their introduction to the structure of bentonite [25,26].Rui et al. [69], on the other hand, argue that freezing-thawing can lead to the disaggregation of soil particles, resulting in a higher proportion of Cd and Pb metals in the thawing front, although they emphasize the significance of soil particle size and specific surface area.Consequently, the correlational relationships between the unfrozen water content at specific negative temperatures with BCR Zn or Cu chemical fractions and clay properties were assessed (Figure 7)..002mm, d 10 -effective diameter is that size below which 10% of the material is contained, SMD-Sauter Mean Diameter is also referred to as d 32 (Laser diffraction method); d 001 -interpacket distance (XRD method); and S-specific surface area (Water sorption test method).The unfrozen water content was taken from the authors' published results for the same bentonite samples (Nartowska et al. [26], Figures 2-4; Nartowska and Kozłowski [47], Table 6).The unfrozen water content was determined using the method of 1HNMR (bentonites contaminated with Zn ions) in the following temperatures: The content of Cu stable fractions (F IV ) in the Cu bentonites increased with increasing unfrozen water content and specific surface area at a given negative temperature (Figure 7B).The trend of changes in correlation coefficients with unfrozen water content for F IV and S was similar at temperatures ≤ −12 • C. It can be assumed that the specific surface area of the clay had a significant effect on the changes in the stable fractions in the given temperature range.The effect of specific surface area on the proportion of Cu stable fractions (F IVCu ) seems less significant at temperatures ≥ −6 • C.This may be because, near 0 • C, unfrozen water is mainly water trapped in mesopores whose degree of filling is less related to the specific surface area of the bentonite.An antagonistic trend was observed in the relationship between the correlation coefficients of unfrozen water content with specific surface area and the correlation coefficients of unfrozen water with effective diameter at individual negative temperatures (Figure 7B).& Unfrozen water content in given temperatures determined using DSC method based on Nartowska and Kozłowski [32].
As the relationship between unfrozen water content and Zn content in the F IZn mobile fraction is statistically significant at temperatures ≤ −14 • C (Table 5), it can be assumed that the relationships observed in Figure 4A for F IZn and d 10 are most significant in this temperature range.Most likely, the decrease in total zinc concentration, which is statistically significant in this temperature range, leads to a decrease in the proportion of F IZn and an increase in unfrozen water content.In turn, the proportion of Zn in the F I mobile fraction decreases with a decrease in the effective diameter, which can be explained by the fact that the volume of connected pores filled with water decreases owing to thawing [70].Thus, Zn ions trapped in the smaller pores become more difficult to leach.At the same negative temperatures, the increase in the proportion of Zn stable fractions with increasing unfrozen water content can also be explained by a decrease in pore size and by trapping some Zn ions in the bentonite interpacket space.This was evidenced by the significant correlation between F IVZn and d 001 .The phenomenon where a lower proportion of mobile Zn fractions occurs at higher unfrozen water contents does not imply the absence of risk to the groundwater environment under negative temperature conditions.Higher levels of unfrozen water are observed in bentonites with lower total Zn ion concentrations.This means that the proportion of mobile Zn fractions is lower in bentonites with lower total Zn concentrations.
In bentonites contaminated with Cu ions, there was a statistically significant relationship (p < 0.05) between the increase in unfrozen water content at individual negative temperatures and the increase in the proportion of stable Cu fractions (F IVCu ) (Table 5).At temperatures ≤ −12 • C, these changes were consistent with changes in specific surface area.The content of unfrozen water at a given negative temperature is higher when the total concentration of Cu ions in the clay is low [47].Thus, higher contamination of bentonites with Cu ions leads to a decrease in the proportion of Cu-stable fractions with a decrease in the amount of unfrozen water.These results confirm previous findings [24] that granulometric parameters are more responsible for the behavior of bentonites contaminated with Zn ions, and that the behavior of bentonites contaminated with Cu ions is influenced by specific surface area parameters.It seems that, at subzero temperatures, relationships with the same clay properties are important.The results presented suggest, under freezing temperature conditions, the presence of concentrated Zn or Cu solutions in bentonites.A study by Wu et al. [71] confirmed that cyclic freezing-thawing enhanced the contact reaction between the solid phase and the liquid, leading to an increase in the removal efficiency of Cd and Pb metals from the front of the freeze-thaw and thus the formation of concentrated solutions.On the other hand, experiments were carried out with the addition of an eluent, EDTA with tartaric acid or citric acid, whose contact with heavy metals in low-permeability clay soils can be limited, which can reduce the effectiveness of remediation.
The main limitations of this study concern the BCR procedure, which is labor-intensive and sensitive to slight changes in the procedure.The extracted fractions, defined operationally, roughly characterize the forms of metal occurrence in the clay.In addition, there is a lack of standardized sample preparation procedures for testing.Despite its many drawbacks, the sequential extraction method is an important source of information on the bioavailability and stability of potentially toxic metal ions in the environment [2].

Conclusions
Based on the analysis of the results, the following conclusions were drawn: 1.
The mobility of potentially toxic metals such as Cu, Zn, Cr, Ni, and Pb in bentonites is primarily determined by the type of main exchangeable cation, whether Zn or Cu, its total concentration, and, to a lesser extent, the origin of the bentonite.Further studies on bentonites containing Cr, Ni, and Pb as the main exchange cations are necessary to assess their mobility; 2.
High total concentrations of Zn ions in Zn bentonites led to increased mobility of Zn ions while simultaneously decreasing the proportion of stable fractions.Conversely, in Cu bentonites, increasing Cu concentrations resulted in decreased Cu mobility in stage 1 of the BCR method, likely due to the binding of Cu ions to the mineral phase of atacamite; 3.
There is a statistically significant relationship (p < 0.05) between Zn or Cu mobility and the structural and physicochemical properties of bentonites.The proportion of Zn mobile fractions (F IZn ) increased with increasing effective diameter, decreasing clay fraction content, and interpacket distances, while the proportion of Cu mobile fractions (F ICu ) increased with decreasing effective diameter and increasing clay fraction content and interpacket distances.Zn ions likely accumulate in larger soil pores, while Cu ions are immobilized in the sorption complex of bentonite; 4.
The stability of potentially toxic Zn or Cu metals in bentonites increases with increasing clay fraction content or, respectively, with increasing specific surface area; 5.
It is possible that, at subzero temperatures, unfrozen water forms a concentrated solution of Zn or Cu, which may pose a threat to the groundwater environment, especially in cold regions.

Data Availability Statement:
The original contributions presented in this study are included in the article; further inquiries can be directed to the corresponding author.

Conflicts of Interest:
The authors declare no conflicts of interest.
Appendix A   & Samples were rinsed from free chloride ions, * Metal content with standard deviation calculated for 4 samples using Grubbs' statistical tests; F I -Exchangeable and weak acid soluble fraction, F II -Reducible fraction, F III -Oxidisable fraction, F IV -Residual fraction determined BCR method with using ICP-OES; Potentially mobile fractions-sum of F I and F II ; Non mobile fractions-sum of F III and F IV ; M f -mobility index; I r -stability index acc.to Formula (A1); Clays were obtained from the Source Clays Repository of The Clay Minerals Society and included STx-1b (a natural Ca-bentonite from Texas) and SWy-3 (a natural Na-Ca-bentonite from Wyoming), BSvk (a natural Ca-bentonite from Jelsovy Potok, Slovakia).
Table A3.Average unfrozen water content of Zn and Cu bentonites.

F III -Metals bound to organic matter and sulfides 4 Dry
remaining sediment at 105-110 • C. Place 1 g of clay in a 300 cm 3 conical volumetric flask.Add 7.5 cm 3 concentrated HCl and 2.5 cm 3 concentrated HNO 3 .Heat until the solution evaporates.Add 25 cm 3 of HCl (1 + 5) to the cooled flasks, transfer the dissolved precipitate to a volumetric flask and fill up to 50 cm 3 with distilled water.Filter the solution and determine the metals in fraction IV.

Materials 2024 , 29 Figure 3 .
Figure 3.The average value of potentially mobile fractions (FI + FII) of metals, depending on the type of bentonite and contamination.The results for BSvk ''Cu'' were shown independently of Cu concentration as the average of all samples tested.

Figure 3 .
Figure 3.The average value of potentially mobile fractions (F I + F II ) of metals, depending on the type of bentonite and contamination.The results for BSvk "Cu" were shown independently of Cu concentration as the average of all samples tested.

Figure 3 .
Figure 3.The average value of potentially mobile fractions (FI + FII) of metals, depending on the type of bentonite and contamination.The results for BSvk ''Cu'' were shown independently of Cu concentration as the average of all samples tested.

Figure 4 .
Figure 4. Contents of potentially toxic metals in BCR chemical fractions, depending on the type of bentonite and type of contamination: (A) clays contaminated with Zn ions and (B) clays contaminated with Cu ions.

Figure 4 .
Figure 4. Contents of potentially toxic metals in BCR chemical fractions, depending on the type of bentonite and type of contamination: (A) clays contaminated with Zn ions and (B) clays contaminated with Cu ions.

Figure 5 .
Different colors indicate the four successive stages of the BCR sequential extraction.Dark blue represents stage 1, while red represents stage 4.

Figure 5 .
Figure 5. Correlation coefficients of clay parameters with individual chemical fractions of BCR metals (F I -F IV ).(A) Zinc fractions in bentonites contaminated with zinc ions.(B) Copper fractions in bentonites contaminated with copper ions.Bentonite parameters tested: Clay-d < 0.002 mm, Silt-d (0.002-0.063 mm), d 10 -effective diameter is that size below which 10% of the material is contained, SMD-Sauter Mean Diameter is also referred to as d 32 (Laser diffraction method); d 001 -interpacket distance (XRD method); S-specific surface (Water sorption test method).

Materials 2024 , 29 Figure 5 .
Figure 5. Correlation coefficients of clay parameters with individual chemical fractions of BCR metals (FI-FIV).(A) Zinc fractions in bentonites contaminated with zinc ions.(B)Copper fractions in bentonites contaminated with copper ions.Bentonite parameters tested: Clay-d < 0.002 mm, Siltd (0.002-0.063 mm), d10-effective diameter is that size below which 10% of the material is contained, SMD-Sauter Mean Diameter is also referred to as d32 (Laser diffraction method); d001interpacket distance (XRD method); S-specific surface (Water sorption test method).

Figure 5 .
Figure 6 shows the relationship between the different chemical fractions and physicochemical properties, depending on the type of metal, Zn or Cu.The results indicate different behavior of these bentonites and clarify the results shown in Figure 5. Different colors indicate the four successive stages of the BCR sequential extraction.Dark blue represents stage 1, while red represents stage 4.

Figure 6 .
Figure 6.The relationship between the content of BCR chemical fractions and the physicochemical parameters of bentonites, depending on the type of main exchange cation, Zn or Cu.

Figure 6 .
Figure 6.The relationship between the content of BCR chemical fractions and the physicochemical parameters of bentonites, depending on the type of main exchange cation, Zn or Cu.
F I -Exchangeable and/or carbonate-bound metals 2 Flood clay samples with 40 cm 3 of 0.5 mol/dm 3 NH 2 OH•HCl, pH = 2. Use HNO 3 to correct reaction.Perform procedure as in step one.F II -Metals combined with amorphous Fe and Mn oxides 3 Add 10 cm 3 of 8.8 mol/dm 3 H 2 O 2 to remaining clay samples after step 2. Place in water bath at 85 • C until solution evaporates.Pour cooled samples into 50 cm 3 ammonium acetate solution (1 mol/dm 3 , pH = 2 after HNO 3 correction).Shake for 16 h and centrifuge to separate clay precipitate from extract.
& Total content of major elements in the dry clay matrix determined with ICP-OES.All potentially toxic elements are contained in TablesA1 and A2; * clay (d ≤ 0.002 mm), silt (0.002 mm < d < 0.063 mm), d 10 (arithmetic mean diameter), and d 32 (Sauter mean diameter) determined with Laser diffraction method; # Specific surface area determined with Water sorption test method; + interpacket distance determined with XRD method; ± Standard deviation.

Table 3 .
Multivariate tests of the significance of clay, contamination, and their interactions for the potentially mobile metal fractions (F I + F II ) of Cu, Zn, Cr, Ni, and Pb.

Table 4 .
Regression analysis between the total concentration of Zn 2+ ions and the amount of Zn 2+ ions in each BCR fraction (F I -F IV ).

Table 5 .
Regression analysis between the unfrozen water content in a given temperature and the amount of Zn 2+ and Cu 2+ ions in a given BCR fraction (F I -F II -F IV ).
Funding:Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.

Table A1 .
Metal content of model clay soils: total and their percentage in chemical soil fractions.Metal content with standard deviation calculated for 4 samples using Grubbs' statistical tests; F I -Exchangeable and weak acid soluble fraction, F II -Reducible fraction, F III -Oxidisable fraction, F IV -Residual fraction determined BCR method with using ICP-OES; Potentially mobile fractions-sum of F I and F II ; Non mobile fractions-sum of F III and F IV ; M f -mobility index; I r -stability index acc.to Formula (A1); Clays were obtained from the Source Clays Repository of The Clay Minerals Society and included STx-1b (a natural Ca-bentonite from Texas) and SWy-3 (a natural Na-Ca-bentonite from Wyoming), BSvk (a natural Ca-bentonite from Jelsovy Potok, Slovakia).: i-the next step of sequential extraction; k = 4; F I -percentage of metal in the i-th chemical form. *where

Table A2 .
Metal content of model clay soils after contamination with Zn or Cu ions: total and their percentage in chemical soil fractions.