The Potential Use of Zeolite, Montmorillonite, and Biochar for the Removal of Radium-226 from Aqueous Solutions and Contaminated Groundwater

: The present work investigated the potential of using zeolite (clinoptilolite), montmorillonite (Swy2), and Conocarpus biochar as adsorbents to remove 226 Ra from aqueous solution. The e ﬀ ect of the initial 226 Ra concentrations on sorbents’ equilibrium activity concentrations and sorbents’ radium removal e ﬃ ciency were investigated. The results showed that zeolite has a higher removal e ﬃ ciency for 226 Ra in comparison with the e ﬃ ciencies of montmorillonite and biochar. In addition to the linear isotherm model, the Freundlich model, followed by Temkin’s model, provided a better description of the adsorption process than the Langmuir model. Kinetic studies indicated that a pseudo-second-order kinetic model could be the best ﬁt for the adsorption of 226 Ra onto the three investigated sorbents, which suggests that the mechanism of adsorption of 226 Ra by sorbents was chemisorption. The intraparticle di ﬀ usion model indicated that adsorption of 226 Ra onto the sorbents involves a multistep process: (i) boundary layer di ﬀ usion and (ii) intraparticle di ﬀ usion. Moreover, the remediation of groundwater samples polluted with 226 Ra was assessed using the investigated sorbents; the results showed that zeolite also has the highest removal e ﬃ ciency among other sorbents. Thus, the low cost, availability, and the high adsorption e ﬃ ciency of zeolite can be a promising sorbent on 226 Ra removal from aqueous solutions and groundwater remediation. Contributions: Conceptualization, M.I.A.-W.; analysis, A.R.A.U., F.I.A., A.S.A.-F.; funding acquisition, A.R.A.U.; investigation, M.I.A.-W., F.I.A., Y.J.A., Z.Q.A. A.R.A.U.; methodology, F.I.A., Y.J.A., and A.S.A.-F.; project administration, A.R.A.U.; writing—original F.I.A., Z.Q.A. A.R.A.U.; writing—review F.I.A., Z.Q.A., A.S.A.-F., M.I.A.-W., and A.R.A.U.


Introduction
In a broad sense, water resources have been subject to pollution from the increases in concentrations of both heavy metals and natural radionuclides. Radium is considered one of the most naturally occurring radionuclides that exist in appreciable amounts in water and groundwater in particular. The amount of 226 Ra in groundwater depends on the type of bedrock structure and its interaction with the groundwater and the prolonged contact between these rocks and groundwater. Thus, water hosted within igneous rocks has a higher radium concentration than those hosted within sedimentary rocks [1,2]. Human activities can also result in the elevation of radium in the environment, which ultimately reaches surface water and groundwater. Such activities include the production of the oil and gas industry [3,4], uranium mining [5], and phosphate fertilizers [6].
The presence of radium at a higher level in the water and its associated health impact on humans has received considerable attention from researchers in the last few decades. The importance of radium comes from the fact that it is one of the most natural radiotoxic radionuclides that emit alpha particles with relatively high linear energy transfer; thus, ingestion of drinking water that contains a high level of radium increases genetic damage, which may lead to cancer, particularly bone cancer, in humans [7,8]. Therefore, several national and international authorities set limits of radionuclides in drinking water. For example, the Environmental Protection Agency (EPA) regulates the combined radium ( 228 Ra + 226 Ra) in drinking water so that it is less than 0.185 Bq/L (5 pCi/L) [9].
Several studies have been conducted worldwide, including in Saudi Arabia, to investigate the naturally occurring radioactive materials (NORMs) in groundwater [10][11][12][13][14][15][16]. To a large extent, most of these studies showed that groundwater contains a significant level of natural radionuclides, especially radium. In Saudi Arabia, in addition to desalination water, groundwater is considered one of the primary sources of drinking water. It contributes almost 40% of the drinking water in the country [17]. Thus, the elevated level of natural radionuclides in drinking water can be considered as the primary source of contamination. Moreover, Saudi Arabia is characterized by a semiarid to an arid climate, which results in water scarcity. Therefore, the agricultural crops of the farms are mainly irrigated from groundwater, which contains an appreciable number of NORMs. Consequently, this will increase the accumulation of NORM in the soil, which ultimately reaches the human body via soil to plant transfer [18]. Thus, it has become important to carry out radionuclide remediation for the contaminated soil by applying a fixation process for these radionuclides in the soil by using clay minerals.
Hence, water pollution due to radionuclides is a global issue. Therefore, several methods have been employed for removing radionuclides from water, including radium. Such methods include chemical precipitation, reverse osmosis, ion exchange, and lime softening. In these methods, the radium removal efficiency from water can reach up to 95% [19].
In recent years, other attractive methods have been used to remove radionuclides from aqueous solutions based on adsorption of the radionuclides from water by applying different natural adsorbents, such as clays, bentonite, and zeolites [20][21][22][23]. The low cost of these natural adsorbents, their availability, high removal efficiency, which could reach almost 100%, and high sorption properties make them attractive to use in water purification from radioactive and other heavy elements.
Among the natural materials, zeolite has been widely used as a natural sorbent to absorb radionuclides and heavy elements from aqueous solutions, mine water, and soil [20,22,[24][25][26]. This is due to the unique porous structure of zeolites, which are classified as aluminosilicate minerals with high ion exchange and sorption properties.
Most radioactivity studies in Kingdom of Saudi Arabia (KSA) showed a high concentration of 226 Ra in groundwater, and the crops irrigated from this groundwater. Thus, it is imperative to carry out radionuclide remediation for both contaminated soil and water. Therefore, several available natural materials such as zeolite, montmorillonite, and biochar have been proposed for this issue.
This study aimed to investigate the potential use of zeolite, montmorillonite, and biochar to remove 226 Ra from aqueous solutions and groundwater in terms of adsorption isotherms kinetics. In addition, 226 Ra adsorption by sandy soil amended with zeolite was investigated.

Sample Preparation
In this work, the natural zeolite used was clinoptilolite (Morgan Co., Morgan, Mountain Green, UT, USA). The montmorillonite (SWy2) was purchased from the Clay Mineral Society (West Lafayette, IN, USA), while the biochar was produced after pyrolyzing the chopped Conocarpus wood wastes at 400 • C using a closed system, which is made of stainless steel (height of 22 cm and diameter of 7 cm), for 4 h. The results showed that the Conocarpus biochar has an alkaline pH of 9.67. The detailed characteristics of the produced biochar were presented by Al-Wabel and Al-Omran [27]. The surface area and pore size characteristics of each investigated adsorbent were determined by surface area and a pore size analyzer (NOVA touch LX4, Quantachrome Instruments, Graz, Austria). The surface area, pore volume, and pore radius of each adsorbent are listed in Table 1. A bottle of one liter of distilled water was spiked with a high-purity 226 Ra standard from the National Institute for Standard Technology (NIST), USA, with the known activity concentration to get four standards with activity concentrations of 1, 3, 5, and 7 Bq L −1 .
All the standard sources were homogenized. A total of 25 mL of the standard source was taken, with a mass of 0.5 g, from all sorbent materials used for the removal of 226 Ra from the aqueous solutions and were then shaken for 24 h. Subsequently, the solutions were filtered by Whatman 45 mm.

Samples Measurement
An 8 mL of the filtered solution was transferred in a high-density polyethylene vial (size 20 mL) with a Teflon cap to prevent radon from escaping. Each of the standard solutions were then added with 12 mL of a high-efficiency mineral oil scintillator cocktail (PerkinElmer, Boston, MA, USA). Afterward, all samples were stored in a refrigerator for at least 30 days to establish a secular equilibrium between 222 Rn and 226 Ra and their respective decay products before the measurement of the radioactivity using a Liquid Scintillation Counter, Quantulus 1220 Ultra-low-level background (LSC) (PerkinElmer, Turku, Finland).
The LSC equipped with a pulse-shaped analyzer that has the proper setting of the α/β discrimination parameter. The stability study of the detector was carried out by using two standards sources, supplied by PerkinElmer (Turku, Finland) covering high energy beta 14 C and low energy beta 3 H. The calibration was conducted with an in-house mixed standard source containing pure alpha ( 241 Am) and pure beta 90 Sr emitters, which were obtained from the National Institute for Standard Technology (NIST), USA. The window setting in the LSC of 226 Ra was conducted by a standard source of 226Ra obtained from the NIST. The lower limit detection (LLD) was calculated by Curie equation, and the activity concentration was found to be 0.07, 0.32, and 0.08 Bq L −1 for alpha, beta, and 226 Ra, respectively. The LLD was achieved under the following conditions: 8 mL sample size, 500 min of counting time, and efficiency of 99% [16,28-32].

Adsorption Isotherm Analysis
The 226 Ra equilibrium adsorption was described by using the isotherm models of linear form and the nonlinear form of Langmuir, Freundlich, and Temkin. The equations of isotherms are shown below.

The linear isotherm equation is given below
where the q e is the quantity of 226 Ra adsorbed on different adsorbent at equilibrium (Bq g −1 ), C e is the 226 Ra concentration at equilibrium in solution phase (Bq L −1 ), and k d is the linear adsorption constant.
The Langmuir isotherm equation in its linear form is given below where q m is the maximum quantity of 226 Ra adsorbed on different adsorbent at saturation (Bq g −1 ), and k l (L Bq −1 ) is the Langmuir constant, which is related to enthalpy. The linear Freundlich isotherm equation is as follows: where k f is the Freundlich constant, which represents the adsorption capacity, and 1/n is the adsorption intensity related to surface heterogeneity of the adsorbent. The linear form of Temkin isotherm is expressed as follows: ; B T is constant related to the heat of sorption (J/mol); b T is Temkin isotherm constant; R is the universal gas constant, and T is the temperature at 298 K.
To evaluate the goodness of fit for the applied isotherm models, the sum of the squares of the errors (SSE), residual root mean square error (RSME), and Chi-square (χ 2 ) test were computed as follows: (q e − q e,m ) 2 (6) where q e,m is the equilibrium capacity obtained by calculating from the model (mg g −1 ), and q e is experimental data of the equilibrium capacity (mg g −1 ).

Adsorption Kinetic Analysis
In order to investigate the mechanisms of adsorption of 226 Ra onto the investigated sorbents, several kinetics models were used, including first-order, second-order, pseudo-second-order, Elovich, power function, and intraparticle diffusion kinetic models.
First-order equation: Second-order equation: Pseudo-second-order equation: Power function equation: Elovich equation: Intraparticle diffusion equation: where q t and q e are the amount of 226 Ra adsorbed by different sorbent (Bq g −1 ) at any time t (min), and at equilibrium, respectively. k 1 is the adsorption rate constant in the first-order model; K 2 is the adsorption rate constant in the second-order model (g Bq −1 min −1 ); k 2 is the rate constant of pseudo-second-order absorption (g Bq −1 min −1 ); k f is the rate constant of power function adsorption (Bq g −1 min −1 ); α is the initial adsorption rate (Bq g −1 min −1 ); β is Elovich constant (g Bq −1 ); a and b are constants; k id is the intraparticle diffusion rate constant (Bq g −1 min −0.5 ).

Adsorption Isotherms of 226 Ra onto Sorbents
Studies on the adsorption isotherm of 226 Ra on the investigated adsorbents were carried out at an initial pH value of 5, sorbent dose 20 g L −1 , and for initial concentration of 226 Ra (1,3,5,7) Bq L −1 . The 226 Ra equilibrium adsorption isotherms were shown in Figure 1. The results show that adsorption of 226 Ra by zeolite, montmorillonite, and biochar tended to increase by increasing the initially added activity concentrations of 226 Ra. It was observed that the adsorption of 226 Ra onto zeolite shows a different behavior than its adsorption onto montmorillonite and biochar. In general, adsorption of 226 Ra to zeolite was characterized by the Giles adsorption isotherm classification H-type, indicating high affinity overall added initial concentrations. Meanwhile, adsorption of 226 Ra to montmorillonite and biochar was characterized by the Giles adsorption isotherm classification C-type, indicating constant partition. The isotherm shapes indicate that the zeolite has a higher potential for 226 Ra adsorption in comparison with those of montmorillonite and biochar [33]. It has been previously reported that natural zeolite with a porous structure and molecular sieves is a promising tool to remediate radium isotopes from waste mine water [20].
The experimental equilibrium data obtained were analyzed using regression models of isotherms Linear, Langmuir, Freundlich, and Temkin, and their various parameters are shown in Table 2. Linear regression analysis for the Langmuir and Freundlich isotherms is commonly used for the best isotherm fitting. In the current study, among all models, Langmuir's model did not describe the adsorption of 226 Ra on three sorbents with a lower r 2 and higher values of SSE or RSME or X 2 compared to other three applied models (Linear, Freundlich, and Temkin). However, Linear's model well explains the adsorption of 226 Ra on three sorbents with a higher r 2 of 0.9875, 0.9998, and 0.9898, respectively, for zeolite, montmorillonite, and biochar. The zeolite showed a much higher Linear distribution coefficient (k d ) value (4.44) for 226 Ra adsorption than the other two sorbents (0.086 for montmorillonite; 0.022 for biochar). Additionally, Freundlich's model showed a good explanation for the adsorption of 226 Ra on three sorbents with r 2 of 0.9289, 0.9998, and 0.9955, respectively, for zeolite, montmorillonite, and biochar. The explained data by Freundlich's model suggest that the adsorption was not uniform site rather than multilayer adsorption, and it was rather heterogeneous than homogeneous. Overall, the Freundlich coefficient values (1/n) were below 1 for montmorillonite and biochar (but not for zeolite), indicating favorable adsorption of 226 Ra to the sorbents, especially for montmorillonite and biochar. In addition, the zeolite showed a much higher k f value for 226 Ra adsorption than others. Temkin models with r 2 of 0.9130-0.9970 can also describe the adsorption of 226 Ra onto the investigated sorbents. The fit of 226 Ra adsorption onto the sorbents to Temkin models suggests a chemical adsorption mechanism with physical forces. Among three sorbents, zeolite showed much higher Temkin equilibrium binding (A T ) and adsorption constants (b T ) as well as a constant relative to the heat adsorption (B T ).  The removal efficiency of sorbents for 226 Ra in relation to the initial added concentrations is shown in Table 3. The percentage removal of 226 Ra was affected by sorbent type and initial activity. The highest removal efficiency values were pronounced for 226 Ra adsorption onto zeolite. The percentage removal of 226 Ra onto zeolite increased as the initial 226 Ra activity concentration was increased, but the percentage removal of 226 Ra onto montmorillonite and biochar tended to decrease with increasing initial activity concentration. A similar finding of high removal efficiency for 226 Ra in aqueous media (98.73%) was reported by Erenturk and Kaygun [23] when using adsorbent the zeolite clinoptilolite modified by polyacrylonitrile. It has been reported previously that partition coefficient (PC) can be applied for determining the performance of different sorbents, mainly to minimize the bias of the adsorption capacity concepts [34]. Therefore, to compare the performance of the synthesized sorbents, the PCs were calculated. Table 3 shows the partition coefficient PC values in relation with the initial added activity concentration of 226 Ra. The highest PC values were pronounced for 226 Ra adsorption onto zeolite (PC of 1061-3132 L kg −1 ), indicating a high reduction in 226 Ra activity by the solid phase of zeolite. Overall, the lowest PC values were more pronounced for the biochar (PC = 21.3-27.5 L kg −1 ) and followed by montmorillonite (PC = 90.6-125 L kg −1 ). It was observed that PC values for adsorption of 226 Ra onto zeolite tended to increase with the increasing initial activity concentration. The higher PC values that were observed for higher initial activity concentrations are mainly explained by the available sites of high selectivity and strong bonding energies for zeolite. However, the lower PC values with decreasing initial activity concentrations indicate that energetically less favourable sites with decreasing activity of 226 Ra. Though previous studies demonstrated that biochar has the potential to remove heavy metals from aqueous solutions [35,36], the obtained results from the current study indicate that biochar has less ability to remove radioactive 226 Ra.

Adsorption Kinetics of 226 Ra onto Sorbents
In the current study, adsorption kinetics of 226 Ra onto zeolite, montmorillonite, and biochar were investigated by using the model of the first-order, second-order, pseudo-second-order, Elovich, power function, and intraparticle diffusion kinetic models. Kinetics of 226 Ra removal onto each sorbent were investigated for time interval range from 60 to 1440 min using a Mechanical shaker at 150 rounds per minute (rpm), initial activity concentration 5 Bq L −1 , sorbent dose 20 g L −1 at an initial pH value of 5. The results are shown in Figure 2. Additionally, the determined parameters of the best one models studied (pseudo-second-order) for 226 Ra are presented in Table 4. 226 Ra adsorption by zeolite was faster than montmorillonite and followed by biochar in the kinetic adsorption experiment. For example, in the first 60 min, most adsorbed quantity (98.3% uptake) was observed, and then, the adsorbed quantity showed a little increase (99.0-99.8% uptake) with increasing equilibrium kinetic time of 120-1140 min. However, for montmorillonite, in 60, 120, and 240 min, respectively, the adsorbed quantity amounted to 65.3, 75.6, and 83.4% uptake), and then, the adsorbed quantity showed an increase of 90.2-94.80% uptake with increasing equilibrium kinetic time of 480-1140 min. Among all sorbents, biochar showed the slower 226 Ra adsorption, which accounted for 29.9 (at 60 min), 48.0 (at 120 min), 66.7 (at 240 min), 75.8 (at 480 min), 83.8 (at 960 min), 86.2% uptake (1440 min). Among all applied adsorption kinetic models, the pseudo-second-order kinetic's model well explains the adsorption data of 226 Ra. It was found that the pseudo-second-order kinetic model provided higher correlation coefficients with experimental results (r 2 of 1.000 for zeolite, 0.9999 for montmorillonite, and 0.9996 for biochar; Table 4), speculating chemisorption of the 226 Ra onto the sorbents. The q e values of 226 Ra derived from pseudo-second-order kinetic were 0.250, 0.242 and 0.233 Bq g −1 for zeolite, montmorillonite, and biochar, respectively. Meanwhile, the k 2 values of the pseudo-second-order kinetic amounted to 3.665, 0.121, and 0.039 for zeolite, montmorillonite, and biochar, respectively. Among all sorbents, the highest q e and k 2 values were pronounced for zeolite. The value of q e obtained from the pseudo-second-order kinetic model for zeolite is identical to that of the experimental value. As indicated by Houhoune Khemaissia [37], the pseudo-second-order model has successfully described radioactive metals onto sorbents. Our findings generally suggest that the rate-controlling mechanism might be chemisorption, indicating the sharing or exchange of electrons between 226 Ra and zeolite. However, the pseudo-second-order kinetics might provide many sorption steps, including external film diffusion, sorption, and internal particle diffusion. In addition, this model cannot explain the specific sorption mechanism. Therefore, the intraparticle diffusion kinetic model was applied to the data of this study. Several studies showed that the intraparticle diffusion plot might represent multilinearity, indicating that two or more steps could occur. The results showed that the data points mostly have two straight lines, and the plots did not pass through the origin (Figure 3). Table 4 shows that the ki d2 values were smaller than the ki d1 values, suggesting that the intraparticle diffusion was predicted to be the rate-limiting step for the 226 Ra adsorption onto the three investigated sorbents. In this study, the deflection of the straight lines from the origin indicated that the intraparticle diffusion is not only the rate-limiting step in the process of 226 Ra onto the investigated sorbents, and other kinetic models, such as surface adsorption, may control the rate of adsorption. Generally, the kinetic data indicated that the mechanism of 226 Ra adsorption by zeolite, montmorillonite, or biochar is complex and can probably occur through a combination process including external mass transfer, intraparticle diffusion through the macropores and micropores as well as sorption processes. It should also be noted that the kinetic process fitted by the pseudo-second-order model suggests chemisorption. In the case of the rate-limited chemisorption, the inner-sphere complexation involved the radioactive metal ions' sorption.

The Sorbent Potential for Groundwater Remediation
The activity concentrations of 226 Ra in collected groundwater samples ranged from 1.12 to 7.67 Bq L −1 (Table 5), which are higher than their corresponding permissible activity concentrations 0.5 Bq L −1 in water for drinking purposes [38]. After the treatment of groundwater samples with the investigated sorbents, the activity concentrations of 226 Ra in groundwater samples reduced from 1.12-7.68 Bq L −1 to 0.04-0.23 Bq L −1 for zeolite, to 0.32-2.68 Bq L −1 for montmorillonite, and 0.74-5.63 Bq L −1 for biochar ( Figure 4).  Among all sorbents, zeolite showed the highest removal efficiency, and thus reduced the activity concentrations of 226 Ra in all investigated groundwater samples to less than 0.5 Bq L −1 (permissible concentrations for drinking water) ( Figure 4). Previously, it has been reported that natural zeolites have already been identified as a significant tool for environmental remediation [39]. Mostly, it is focused on the ion exchange properties of zeolite. Zeolites with many exchangeable sites and small pores are suitable for processes of sorption [26]. Therefore, it has been reported that zeolites are useful for the removal of radioisotopes from aqueous nuclear wastes [40,41].

226 Ra Adsorption by Sandy Soil Amended with Zeolite
After mixing zeolite in sandy soil for 21 days, the 226 Ra adsorption was carried out with an initial activity concentration of 5.0 Bq L −1 . In Figure 5, considering the control of sandy soil (without the addition of zeolite), the removal efficiency of 226 Ra accounted for 18%. Meanwhile, after mixing zeolite in sandy soil, the removal efficiency for 226 Ra showed significant increases. The results showed that the removal efficiency of 226 Ra increased from 18% in control soil to 73, 88, 94, 96, 98% in soil treated with zeolite at an application rate of 1, 5, 10, 15 and 20% (w/w), respectively.

Mechanism of 226 Ra Adsorption onto Natural Zeolite (Clinoptilolite)
The kinetic study showed that the obtained data of 226 Ra adsorption onto zeolite (clinoptilolite) were described well by the pseudo-second-order model (Table 4), whereby the correlation coefficient is 1.000. This type of kinetic model was used to describe the process of adsorption as chemical adsorption (chemosorption) occurs, and it has been widely applied to describe the pollutant adsorption from aqueous solutions [42].
It has been reported that the adsorption capacity of the zeolite depends on the porosity as well as the chemical reaction between the pollutants and the active functional groups [43]. In the current study, the infrared spectra (FT-IR) for zeolite were determined in a range of 400 to 4000 cm −1 . The bands appeared at 3634-3780, 1625, 1024, 784, and 709 cm −1 (Figure 6). From these obtained results, it is clear that the bands at 3634-3780 cm −1 and 1625 cm −1 confirm the presence of O-H stretching. The FT-IR spectra of clinoptilolite also show a stretching band at 1024 cm −1 , which belongs to an asymmetry stretch (AlO or SiO stretching in AlO4 or SiO4). Meanwhile, the appeared bands at 784 and 709 cm −1 indicate the presence of symmetry stretch (O-Al-O or O-Si-O stretching). The FT-IR results for clinoptilolite showed a change in the band sites or its intensity after the adsorption of 226 Ra ( Figure 6). Following 226 Ra adsorption onto clinoptilolite, the stretching vibration of the AlO or SiO at 1024 cm −1 shifted to 1054 cm −1 with a lower intensity compared to its corresponding band appeared on clinoptilolite samples before adsorption. Similarly, the stretching vibration of the O-Al-O or O-Si-O stretching at 784 cm −1 showed a lower intensity compared to its corresponding band which appeared on clinoptilolite samples before adsorption. Additionally, the OH stretching vibrations at 1625 cm −1 and 3634-3780 cm −1 were diminished following 226 Ra adsorption onto clinoptilolite. This suggests that many interactions with the functional groups of the zeolite may be responsible for contributing to the adsorption process of 226 Ra onto natural zeolite.
In this context, different mechanisms could be proposed for 226 Ra adsorption on zeolite (Figure 7). It could suggest that the zeolite contains various functional groups of hydroxyl and oxygen, which have the capacity to bind 226 Ra via surface complexation, electrostatic attraction, and/or ion exchange.

Conclusions
In this work, an investigation of the adsorption isotherms of 226 Ra onto zeolite, montmorillonite, and biochar equilibrium sorption was carried out. Four adsorption isotherm models were studied: Linear, Langmuir, Freundlich, and Temkin isotherms. In addition to the linear model, the Freundlich adsorption model was found to be the best fit. Among the investigated sorbents, zeolite showed a better performance for 226 Ra removal from aqueous solutions and groundwater and soil remediation. Kinetic adsorption for the removal of 226 Ra from aqueous solution was obtained for the investigated sorbents and fitted to different kinetic models. The results of the kinetic models reveal that the pseudo-second-order model has the best fit for the adsorption of 226 Ra onto the three investigated sorbents, especially zeolite, indicating the chemisorption of the 226 Ra onto the sorbents. However, the mechanism of 226 Ra adsorption by zeolite, montmorillonite, or biochar is complex and can probably occur by a combination process including sorption processes and/or external mass transfer and intraparticle diffusion through the macropores and micropores. It could be concluded that zeolite is an effective sorbent for the removal of 226 Ra from aqueous solution and groundwater and soil remediation.