The Simultaneous Removal of Zinc and Cadmium from Multicomponent Aqueous Solutions by Their Sorption onto Selected Natural and Synthetic Zeolites

: Natural and synthetic aluminosilicate minerals, in particular zeolites, are considered to be very useful in remediation processes, such as purification of waters polluted with heavy metals. That is due to their unique and outstanding physico-chemical properties, rendering them highly efficient, low-cost, and environmentally friendly sorbents of various environmental pollutants. The aim of this study was to examine the sorption capacity of four selected zeolites: A natural zeolite and three synthetic zeolites (3A, 10A, and 13X), towards zinc and cadmium present in multicomponent aqueous solutions, in relation to identified sorption mechanisms. It was stated that synthetic zeolites 3A and 10A were the most efficient in simultaneous removal of zinc and cadmium from aqueous solutions. Additionally, zeolite 10A was demonstrated to be the mineral best coping with prolonged pollution of water with those elements. The mechanism of sorption identified for tested minerals was physisorption.


Introduction
The observed, in the last two centuries, rapid technological development, pronouncing itself in growing industrialization and urbanization, resulted in the rapid increase of pollution of the environment with various chemicals, including trace elements such as heavy metals. Their ongoing accumulation in various environmental compartments, resulting from a systematic deposition, is now considered to be a permanent phenomenon. From that arises a severe danger and an onerous problem to all environmental compartments, including, in the aquatic compartment, that of the permanent pollution, resulting in the inclusion of the toxic elements to the food chains, impacting the proper functioning of the ecosystems and putting into high risk human health [1][2][3][4]. For that reason, the concentrations of many of those pollutants, including heavy metals, are continuously monitored by relevant national and international organizations. In those activities special consideration is given to water, one of the most important renewable resources of the Earth. Water is one of the most precious and most commonly used renewable resources. Abundant in the Earth's crust, it covers two thirds of the planet's surface [5]. However, only 0.3-0.5% of that amount is fresh water accessible to human beings, and that volume diminishes due to pollution from various point and spatial sources. Therefore, so important becomes the issue of water protection by reducing the emissions from identified and potential sources and remediation of polluted surface and ground water using different, efficient techniques.
Heavy metals are listed among the most important environmental pollutants, due to their chemical nature. As chemical elements they are not prone to degradation in the environment, only to change of the form of their occurrence. As a result, they tend to accumulate in various environmental compartments thus becoming permanent pollutants. They enter the environment as a result of various human activities, metallurgy, chemical, and electro-technical and engineering industries, but also agriculture. Among most important of their point sources are listed various untreated industrial effluents as well as the accidental emissions from industrial and municipal wastewater treatment installations and dumping sites. Major spatial sources are indicated as agricultural and municipal effluents as well as dry and wet atmospheric deposition.
Zinc and cadmium are two heavy metals ranked highly among the environmental pollutants, due to their wide abundance and toxicity [6]. Belonging to the same group in the periodic table, IIb, those two elements have similar properties and, hence, display very similar behavior in the environment, for instance, in soil [7][8][9]. Lighter zinc, although an important microelement of well examined roles in living organisms [10], in elevated concentrations is known to adversely affect living organisms. Heavier cadmium has no known biological role, but well demonstrated high toxicity to all living organisms already in low concentrations. Both elements display a tendency to accumulate in living organisms and ecosystems, as well as a long residence time (in the case of cadmium, estimated at 10 to 33 years), which causes severe health problems. For that reason, so important becomes removing those compounds from various environmental matrices, in particular water, the medium in which they are most readily bioavailable.
The possibilities of the efficient removal of heavy metals from polluted waters and of the limiting of the inclusion of heavy metals to the food chain at the level of plants-soils interactions was the subject of numerous studies, also performed by the authors of this work [11][12][13][14][15][16][17][18][19][20][21][22][23][24]. Those findings were at the basis of the experimental hypothesis of this research, which assumed that natural and synthetic aluminosilicate minerals were highly efficient, low-cost, and environmentally friendly materials in removal of heavy metals from water. The possibility of removal of heavy metals from aqueous solution by means of their sorption onto various materials, including mineral sorbents, was extensively examined [22,25,26], but most of those research activities concentrated on single-solute systems. It was also demonstrated that the sorption of metal ions from aqueous multisolute systems onto various sorbents may result in different sorption parameters than those obtained for the same chemicals in single-solute systems [27][28][29].
The main aim of the study was the assessment of the possibility of use of selected porous aluminosilicate minerals having a network of pores and chambers, zeolites, to decontaminate water containing elevated levels of zinc and cadmium, as observed in the case of accidental emissions in highly concentrated effluents bearing various metal ions. This was done with the aim to provide a deeper insight into the problem of simultaneous sorption of metal elements with almost identical chemical behavior from binary system onto sorbents belonging to the same class but varying in sorptive capacities. Also, the results should provide a clear indication for practitioners on which sorbent to use in case of accidental high-concentration releases of zinc and cadmium to the aquatic environment, hitherto seldom done in scientific papers.
To achieve that, four selected zeolites, one natural and three synthetic, were examined for their sorption capacity and selectivity towards zinc and cadmium ions present in multicomponent aqueous solutions. An attempt was also made to identify the possible mechanisms of sorption of Zn 2+ and Cd 2+ by the tested minerals. To obtain a clear and unambiguous answer to the question of the efficiency of each tested mineral in removal of the elements of concern from remediated water, the concentration of the test metals in the study were correlated with the cation exchange capacity (CEC) factor of the minerals undergoing examination.

Materials and Methods
In this study we used four porous aluminosilicate minerals-zeolites. They were crystalline hydrated minerals displaying high variability of internal structure, containing characteristic empty spaces, channels and chambers, filled with ions and water molecules having a high degree of freedom of movement.
Chemically, zeolites are hydrated aluminosilicates of metals belonging to Ia (Li, Na, K) and IIa (Mg, Ca Sr, Ba) groups of the periodic table, which compensate the negative charge. The simplified chemical formula of any zeolite can be presented as follows: where Me is metal cation, n is valency of the metal cation, and x is molar ratio SiO2/Al2O3.
From the structural point of view, the basic elements of their crystalline network are tetrahedras [SiO4] 2− and [AlO4] − linked by the oxygen bridges to form the polyhedras, e.g., cuboctahedras. The spatial arrangements of cuboctahedras and their binding structures result in different structural types of zeolites, e.g., A, X and Y, all belonging to the faujasite group.
Zeolites are also called molecular sieves, because of their ability of selective sorption of chemical molecules that are smaller than their sorption holes.
In this study we used the natural zeolite, clinoptilolite, obtained from Caucasian deposits, and three synthetic zeolites, 3A, 10A, and 13X.
The natural zeolite, clinoptilolite, has dominant exchangeable cations K + and Ca 2+ and the pores' diameter of 0.44-0.55 nm. To obtain higher homogeneity it was ground to obtain the grain size of <0.2 mm.
Of the three synthetic zeolites, two are the representatives of type A and the third is an X-type zeolite. Their exemplary structural models are presented below in Figure 1.

Structural Model of A-Type Zeolite Structural Model of X-Type Zeolite
A cubeoctahedr A single crystal A cubeoctahedr A single crystal The tested synthetic zeolite 3A is a sodium-and-potassium zeolite with pores' diameter of 0.38 nm, zeolite 10A is a sodium zeolite with pores' diameter of 0.9-1.0 nm, and zeolite 13X is a sodium zeolite with pores' diameter of 0.9-1.0 nm.
The basic properties of all characterized above aluminosilicate minerals are presented below in the Table 1. The examination of the sorption of zinc and cadmium onto four tested minerals was performed using multisolute aqueous solutions containing Zn 2+ and Cd 2+ , in initial concentrations corresponding to 2% CEC, 10% CEC, 20% CEC, 30% CEC, 50% CEC, 75% CEC, and 100% CEC (cation exchange capacity) of the given mineral sorbent. Those solutions were prepared from the stock solutions using the serial dilution method.
Stock solutions were prepared individually for each tested mineral by dissolving the appropriate samples of solid nitrate salts of the test elements (analytical grade) in deionized water to obtain the mixture having the total concentration of test ions equal to 100% CEC of the given mineral sorbent.
The experiment was performed to comply with the methodology outlined in the OECD Guideline 106 [30]. It consisted of two stages, initial, at which were determined the adequate equilibration time and the incubation temperature, and the definitive stage, aimed at the determination of the sorption capacity of each tested mineral towards each of the test metal ions.
The whole experiment was performed using a constant, predefined mineral sorbent with solution ratio of 1:25. For the definitive test, aimed at the determination of the sorption capacity of each mineral towards each cation, seven polymetallic solutions were prepared. The experiment was performed in a batch mode. All samples were equilibrated by shaking at constant speed on a horizontal shaker placed in a water bath set to the constant temperature. After equilibration, the samples were centrifuged, clarified supernatants decanted, and the content of metal ions in them determined by AAS to obtain the equilibrium concentrations in solution-Ce.
The equilibrium concentrations of each metal adsorbed onto sorbent were calculated from a difference between the initial concentration in the solution and that determined experimentally after equilibration.
Each sample was prepared in four replicates. For each mineral and each initial concentration, we determined the percent sorption and the distribution coefficient, Kd.
Those two parameters were calculated using the following equations.
1. For the determination of the % of sorption: where CINI is the initial concentration of the given metal ion in the test solution in [cmol/L] Ceconcentration of the given metal ion in the solution at equilibrium in [cmol/L].
2. For the determination of the distribution constant Kd: where Ce is the concentration of the given metal ion in the solution at equilibrium in [cmol/L] and Csorb is the concentration of the given metal ion adsorbed onto the mineral sorbent at equilibrium in [cmol/kg]. The data obtained for each mineral were also analyzed with the aim to determine the parameters of the sorption isotherms, characterizing sorption capacity, and the mechanisms of the process. Three isotherm models were used: Freundlich isotherm, Langmuir isotherm, and Dubinin-Kaganer-Radushkevich (DKR) isotherm. All three isotherms are commonly used to characterize the process of monolayer sorption from solutions onto heterogeneous surfaces.
In their canonical forms they are represented by the following equations.
Freundlich isotherm: x/m = Kf × Cs 1/n where x/m is the amount of the sorbate sorbed at equilibrium, Cs is the concentration of the sorbate in solution at equilibrium, Kf is the Freundlich isotherm sorption constant, and 1/n is the Freundlich sorption coefficient.
where x/m is the amount of the sorbate sorbed at equilibrium, Cs is the concentration of the sorbate in solution at equilibrium, KL is Langmuir's isotherm sorption constant, and Ns is the number of sorption sites (also named "maximum adsorption capacity").
Although, for convenience, they are usually used in their linearized form, as logarithmically transformed equations, for the purpose of this work the data were fitted to the presented above canonical forms. That is due to several limitations of the linear regression methods in comparison to the nonlinear regression ones [44][45][46][47][48].
The third isotherm used in the study, the DKR (Dubinin-Kaganer-Radushkevich) isotherm, is a model characterizing the adsorption onto micropores [49,50]. It is used to determine two important parameters, the maximum sorption capacity of the tested sorbent towards the given adsorbate (Xm) and the mean adsorption free energy E, an indicator of the mechanism of the examined process [51][52][53][54]. Although generally applicable to single-solute systems, that model was demonstrated to be successful in the characterization of sorption in multicomponent systems [55].
Being in its canonical form the exponential equation [56][57][58], for convenience it is commonly used in linearized form (after [59]): where Cads is the amount of metal ions adsorbed per unit weight of sorbent (here in [cmol/kg]), Xm is the maximum sorption capacity, β is the activity coefficient related to mean sorption energy, and ε is the Polanyi potential. Polanyi potential is expressed using the following equation: where Ce is the amount of metal ions remaining in solution at equilibrium The value of the energy is used as the indicator of the mechanism of the process, physisorption, ion exchange, or chemisorption [37,60,61]. It is calculated from the activity coefficient using the following equation: All three sorption isotherm models used in this study were determined using the tool CurveExpert Professional ver. 2.6.5. In the case of Freundlich's and Langmuir's isotherms, the nonlinear regression method was used, while the linearized DKR isotherm was optimized using the linear regression. For Freundlich's isotherm, the built-in tool, "Power" function, was used, Langmuir's isotherm was optimized using the user's model, and for the DKR isotherm we used the model built-in linear function.\

Results and Their Discussion
The preliminary experiment showed that the equilibrium state was reached after 1 h. That time was, therefore, selected as equilibration time in further experiments.
The examination of the influence of the temperature on the sorption showed that the optimum sorption was attained at T = 20 °C. Therefore, for the definitive experiment, that temperature was selected, as it was considered to be representative for the average environmental conditions.
For each tested combination Me 2+ -mineral sorbent and each tested concentration, we calculated the percentage of sorption and the values of the distribution coefficient Kd. The results are presented numerically in the Table 2. The comparative analysis of the above results showed that for both sorbed elements there was no one, clear behavior pattern of all tested minerals, which is related to their structural properties. In general, it can be stated that while the sorption capacity of natural zeolite and synthetic zeolite 3A was greater towards Zn 2+ than towards Cd 2+ , that tendency was reversed in the case of the synthetic Zeolites 10A and 13X.
Additionally, it was stated that in the case of the zeolite 13X for both tested elements, the level of sorption steadily decreased with increase of the initial concentration. The same was observed for the sorption of Cd 2+ onto natural zeolite. The sorptive behavior of Zn 2+ onto that zeolite was different: The sorption at low concentrations was at the level of 50%, rising to 60%, and then gradually decreasing with increasing concentration.
In the case of the zeolite 3A, the sorption pattern of both elements was similar: It was at the relatively constant levels in the range of ~30%-~60%, without any clear increase-decrease trend.
A similar picture was obtained for the zeolite 10A, but in this case the sorption of both elements was low at low concentrations tested to sharply increase and then remain relatively stable at the level ~30%-~50%.
The observed similarities and differences in sorption can be explained in terms of the mechanisms involved in the phenomena. It is generally acknowledged that the adsorption preferences of sorbent to sorbate may be due to the following factors [28]: (1) Chemistry of the solution, such as pH and ionic strength; (2) Properties of sorbate concentration, ionic size, ionic charge, ionic weight, and standard redox potential; (3) Character of binding sites, functional groups, surface properties, etc.
At the same time, it remains not fully clear how those factors may influence the preferences of the given sorbent towards the given sorbate, in particular in the multicomponent system.
The explanation of the observed phenomena was related, firstly, to the properties of sorbates. As already indicated, Cd 2+ and Zn 2+ are the elements belonging to the group IIb of the periodic table. Therefore, having some key, from the point of view of their sorption onto zeolites, properties that are the same or similar, they have the same coordination number (6) and charge (+2) and similar diamagnetism and electronegativity, 1.7 for Cd 2+ and 1.6 for Zn 2+ . As a result, it may be expected to have competition for the sorption sites. However, those two elements differ in other properties considered crucial for sorption onto zeolites, namely [37,62,63]: (1) Atomic weight, for Zn 2+ M = 65.38 [g/mol] and for Cd 2+ M = 112.41 [g/mol]; (2) Hydrolysis constant, for Zn 2+ Ka = 9.0 and for Cd 2+ pKa = 10.1; (3) Ionic radius, for Zn 2+ it is 0.074 nm and for Cd 2+ 0.097 nm; (4) Hydrated ion radius, for Zn 2+ it is 0.430 nm and for Cd 2+ 0.426 nm; (5) Hydration energy, also called hydration enthalpy, being −1955 kJ/mol for Zn 2+ and −1807 kJ/mol for Cd 2+ . The generally low levels of sorption of both Cd 2+ and Zn 2+ onto tested zeolites may be explained by the radii of their hydrated ions and hydration energies. The necessary condition for the given cation to be bound to the sorption sites of the given sorbent is the removal of water molecules from the hydration sphere, which is related to the hydration enthalpy in the following manner: The higher the hydration energy, the lower the possibility of the removal of water molecules from hydration sphere and, hence, binding of the given cation to the sorption sites. Zn 2+ and Cd 2+ have high, in comparison to other bivalent metal cations hydration, energies, which explains the observed relatively low sorption levels and sharp decrease of sorption with increasing concentrations of sorbates [64].
The differences in hydration energies, hydrated ion radii, and atomic weights may also explain some preference in sorption of Cd 2+ ions to Zn 2+ ions observed for synthetic zeolites 10A and 13X, in line with findings reported in the literature on this subject [20,65].
The low sorption capacity of natural zeolite, in comparison to the synthetic zeolites, is in line with the other literature-reported results [66,67].
Because the results of the determination of the percentage of sorption and Kd values did not provide a clear answer to the main problem, the identification of the most efficient of the four tested zeolites in simultaneous removal of Zn 2+ and Cd 2+ ions from the multisolute aqueous solutions, the obtained data were analyzed using three isotherm-based models. Freundlich, Langmuir, and DKR, demonstrated to be suitable for examination of sorption in multisolute systems [38].
Freundlich's isotherm provided the information on the strength and extent of sorption, as well as on the nature of the process. From the Langmuir isotherm, we derived the information on the maximum sorption capacity. Finally, the results of the DKR isotherm were used to identify the possible mechanism of sorption as well as sorption capacity.
Below are provided the results of the determination of the Freundlich and Langmuir isotherms. The numerical results are provided in Table 3 for Freundlich's isotherms and Table 4 for Langmuir's isotherms.  In graphical form, the isotherms are presented in Figures 2 and 3. Each figure presents a set of the individual isotherms, generated by the tool CurveExpert Pro 2.5.6. The colored bands represent the confidence bands for a 95% level of confidence, darker pink, and 90% level of confidence, lighter pink.

Isotherms for Zn 2+
Isotherms for Cd 2+ Figure 2. The Freundlich isotherms obtained in the study for sorption of Zn 2+ and Cd 2+ onto tested natural and synthetic zeolites.

Isotherms for Zn 2+
Isotherms for Cd 2+ Figure 3. The Langmuir isotherms obtained in the study for sorption of Zn 2+ and Cd 2+ onto tested natural and synthetic zeolites.
On the basis of the analysis of the goodness of fit (visual inspection of plotted isotherms and the analysis of the r and R 2 values) for each isotherm, it can be stated that both models well complied with the experimental data. The comparative analysis of the determined values of the r and R 2 values showed that in general Freundlich's model better characterized the adsorption of Zn 2+ and Cd 2+ onto the tested minerals. For that reason, its parameters were selected to characterize the strength and extent of sorption.
The analysis of the strength of sorption, based on the determined Kf values, showed that it cannot be stated for all four tested zeolites whose element, Zn 2+ or Cd 2+ , was sorbed more strongly, hence, which displayed higher affinity towards the tested aluminosilicates. The stronger sorption of Zn 2+ in comparison to that of Cd 2+ was observed in case of natural zeolite and zeolite 3A, while Cd 2+ was sorbed more strongly than Zn 2+ on zeolite 10A and zeolite 13A. It has to be indicated that in the case of zeolite 13X the difference in sorption strength of both elements was not that big, as observed for the other tested mineral sorbents. The observed differences may be attributed to some specific properties of the sorbed ions, indicated in the available literature, atomic weight, ion radius, hydratation energy, and specific properties of the aqua complexes of each ion.
It was also stated that the strongest sorption of both Zn 2+ and Cd 2+ occurred on synthetic zeolites 3A and 10A, while the weakest was on the natural zeolite.
When arranged in order from the strongest to the weakest the strength and extent of sorption followed the order: The values of 1/n, the parameter informing about the nature of the process, was always below 0.8, indicating that the sorption of both Zn 2+ and Cd 2+ on all tested minerals was generally preferential and either favorable or pseudolinear [68]. When arranged from the highest to the lowest, the values of 1/n, and, hence, the decrease of the linearity of sorption followed the order: This may indicate, in the case of the combined pollution of the aquatic environment with those two elements, synthetic zeolite 10A would remain most efficient either over a much broader range of concentrations or in the case of prolonged exposure to a polluted medium, while the least efficient is expected to be the synthetic zeolite 13X. The sorption strength, characterized by the Freundlich adsorption constant, was well correlated with the maximum sorption capacity N, the parameter determined using the Langmuir sorption isotherm. It was stated that for both Zn 2+ and Cd 2+ this parameter was higher for synthetic zeolites than for the natural zeolite and that the A-type zeolites had much greater sorption capacity than the X-type zeolite.
It may be also stated that using the values of Kf from the Freundlich isotherm and N from Langmuir's isotherm, indicators of the selectivity of the tested zeolites, it may be postulated that while natural zeolite and zeolite 3A displayed higher selectivity towards Zn 2+ ions, the synthetic zeolites 10A and 13X showed preference towards Cd 2+ ions.
When arranged in order from the highest to the lowest, the Langmuir's maximum sorption capacity N followed the order: For both Zn 2+ and Cd 2+ , Zeolite 10A > Zeolite 3A > Zeolite 13X > Natural Zeolite It was also indicated that, while for the natural zeolite the difference in N values was significant, about 2.5 times in favor of Zn 2+ , in the case of the tested synthetic zeolites the differences in N values were not that big, and followed the order: Zeolite 10A (51% difference in favor of Cd 2+ ) > zeolite 3A (11% difference in favor of Zn 2+ ) > zeolite 13X (5% difference in favor of Cd 2+ ).
The results of the determination of the parameters of the DKR isotherm are presented below, in numerical form in Table 5. The analysis of the goodness of fit of each isotherm, based on the analysis of R 2 and r coefficients, showed that at least acceptable fits were obtained for each combination mineral, sorbed element, which indicated, as it was in the case of Freundlich and Langmuir sorption models, that the DKR was correctly selected as an adequate model characterizing sorption in the examined test systems.
The analysis of the constant Xm, characterizing the maximum sorption capacity, showed that that parameter was higher for sorption of Zn 2+ for the natural zeolite and synthetic zeolite 3A, while the opposite situation was observed for the synthetic zeolite 13X. In the case of the synthetic zeolite 10A, the values of Xm for both elements were almost identical.
When arranged in order from the highest to the lowest, the maximum sorption capacity Xm followed the order: For Zn 2+ , that order was very similar, and for Cd 2+ , it was identical to that determined for respective N values in the Langmuir isotherms, but the values were different, which can be explained by the fact that while Langmuir's isotherm examines the process of sorption in a more general way, the DKR isotherm is more focused on the examination of sorption in micropores. When the ratio N:Xm was calculated, N was always higher than Xm. The numerical results of that analysis, the calculated ratios N:Xm, are provided below in the Table 6.
The values of the energy of sorption of Zn 2+ on the tested sorbents were in the range of 2.19-4.52 [kJ/mol] and those for Cd 2+ were in the range 2.29-4.94 [kJ/mol].
That parameter provides the information about the possible mechanism of sorption [37,60,61]. That was done using the following classification of the postulated mechanisms of sorption in relation to the determined energy: For E < 8 kJ/mol, the possible mechanism of adsorption was physical sorption; for E being in the range of 8 and 16 kJ/mol, it can be explained by ion exchange; for E > 16 kJ/mol, the possible mechanism was chemical sorption stronger than the ion exchange.
On its basis, it can be stated for both elements the sorption mechanism for all tested minerals was physisorption. Although, in principle, ion exchange was indicated as a primary mechanism of sorption of metal ions onto aluminosilicate minerals, including zeolites, in the case of sorption from multicomponent solutions the mechanism could switch to the physisorption, which in turn may be attributed to the higher competition for the low-energy sorption sites [27].
When arranged in order from the highest to the lowest the sorption energy E followed the order: Additionally, the values of the maximum adsorption capacity towards Zn 2+ and Cd 2+ ions, determined for each tested mineral using Langmuir's and DKR isotherms, were collected in a single table, Table 6 below. In order to render them comparable with the values reported in the literature, those values were converted to [mmol/g] and to [mg/g], units most commonly used in the scientific literature of the subject. The second conversion was performed using the adequate values of molar weights, 65.38 g/mol for Zn and 112.41 g/mol for Cd. When compared to the values provided in other scientific papers on the same subject [19,[22][23][24]26,[69][70][71], the tested minerals were demonstrated to display similar or greater sorption capacity comparing to the similar sorbents and other novel materials tested to remove heavy metals from wastewater effluents. That makes them a good, costefficient, and environmentally friendly alternative in the processes of purification of water polluted with heavy metal ions.

Conclusions
On the basis of the obtained results it may be stated that: (1) The tested synthetic zeolites, 3A, 10Am and 13X, were more efficient in the simultaneous removal of Cd 2+ and Zn 2+ ions from aqueous solutions than the tested natural zeolite, which may be related to their higher mineralogical homogeneity; (2) The analysis of the performance of the tested zeolites, based on the analysis of the parameters of sorption isotherm models, showed that the most efficient in simultaneous removal of Cd 2+ and Zn 2+ ions from aqueous solutions were zeolite 3A and zeolite 10A. Therefore, those two sorbents should be recommended to be used for rapid reduction of the level of pollution with those two elements and their spreading in the environment; (3) Zeolite 10A displayed a relatively high and constant sorption capacity over a broader range of concentrations, which indicates that it will be efficient in coping with prolonged low-and medium-level pollution of aquatic environment with cadmium and zinc; (4) The identified, on the basis of the adsorption energies, mechanism of sorption was, for all tested zeolites, physisorption, which may be explained by the fact that the sorption was examined in a multisolute system; (5) The synthetic zeolites 10A and 13X displayed higher selectivity to Cd 2+ than to Zn 2+ , while in the case of the natural zeolite and synthetic zeolite 3A that relationship was reversed. Funding: This research received no external funding.