Competitive Adsorptive Removal of Pb, Cu and Cd from Water Using Water Pennywort Fixed on Alginate

Mahamadi C; Mupa M; Jackson GE

doi:10.4172/2155-6199.1000300

ISSN: 2155-6199

Journal of Bioremediation & Biodegradation

Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.

Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on Medical, Pharma, Engineering, Science, Technology and Business

Competitive Adsorptive Removal of Pb, Cu and Cd from Water Using Water Pennywort Fixed on Alginate

*Mahamadi C^1,, Mupa M² and Jackson GE³**
¹Research and Postgraduate Centre, Bindura University of Science Education, P. Bag 1020, Bindura, Zimbabwe
²Chemistry Department, Bindura University of Science Education, P. Bag 1020, Bindura, Zimbabwe
³Department of Chemistry, University of Cape Town, Rondebosch 7700, South Africa
Corresponding Author :	Mahamadi C Research and Postgraduate Centre Bindura University of Science Education P. Bag 1020, Bindura, Zimbabwe Tel: +26377350801 E-mail: cmahamadi@buse.ac.zw
Received June 07, 2015; Accepted June 24, 2015; Published June 26, 2015
Citation: Mahamadi C, Mupa M, Jackson GE (2015) Competitive Adsorptive Removal of Pb, Cu and Cd from Water Using Water Pennywort Fixed on Alginate. J Bioremed Biodeg 6:300. doi:10.4172/2155-6199.1000300
Copyright: © 2015 Mahamadi C, et al. This is an open-a ccess article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Related article at Pubmed Scholar Google

Visit for more related articles at Journal of Bioremediation & Biodegradation

View PDF Download PDF Tables & Figures

Abstract

The adsorption of Pb (II), Cu (II) and Cd (II) on immobilized water pennywort-Hydrocotyle ranunculoides was investigated in batch and continuous flow systems. Physical characterization of adsorbent was performed using fourier transform infrared (FT-IR) spectroscopy and scanning electron microscopy (SEM). Analysis of batch sorption data using single component and competitive Langmuir models showed that maximum equilibrium sorption capacities for single-metal systems followed the trend (mg/g): Pb(II) (22.25)>Cu(II) (12.55)>Cd(II) (11.45). Mutual interference was probed using normalised equilibrium adsorption capacities, qe'/qe, where the prime indicates the presence of the other metal ion. The results showed that removal of Cu (II) and Cd(II) ions was severely suppressed in the presence of Pb(II) whereas the removal of Pb(II) in the presence of Cu(II) and Cd(II) ions was still effective. Adams-Bohart and Thomas models were used to study the effect of important parameters: bed depth (1.5-2.5 cm), flow rate (2-20 ml/min) and initial metal concentration (10-100 mg/l). It was shown that the Thomas model provided a better fit to the experimental data indicating that mass transfer, i.e. external diffusion was not the rate limiting steps. In conclusion, immobilised H. ranunculoides proved to be a potential biosorbent for the removal of heavy metal ions from aquatic environments.

Keywords

Hydrocotyle ranunculoides; Continuous flow; Binary sorption; Lead; Copper; Cadmium

Introduction

The question of toxicity and health hazards associated with heavy metal accumulation in aquatic systems has generated sustained research over centuries. Biosorption, the process of passive cation binding by dead or living biomass, represents a potentially cost-effective way of eliminating toxic heavy metals from industrial waste waters. The initial incentives of biosorption development for industrial processes include the low cost of biosorbents, their great efficiency for metal removal at low concentration, potential for biosorbent regeneration and metal recovery, high velocity of sorption and desorption, limited generation of secondary residues, and more environmentally friendly life cycle of the material [1].

The complexity of the biochemistry of the biosorbents has led to numerous interpretations on the mechanisms involved in the biosorptive accumulation of metal ions [2]. However, it has generally been agreed that metal binding takes place through electrostatic interaction, surface complexation, ion-exchange, and precipitation which can occur individually or combined [3]. On parameters that affect biosorption, it has generally been observed that temperature does not influence biosorption between 20 and 35°C whereas Ph affect biosorption through its influence on metal chemical speciation, the activity of biomass functional groups, and the ion metallic competition sites [4]. In dilute solutions, biomass concentration has been found to have an influence on biosorption capacity with an increase in biosorption being observed at lower concentrations.

¨Although much research has been carried out on the uptake of heavy metal ions by many researchers, this has been largely limited to single-metal adsorption. However, under most conditions, wastewater streams contain more than one kind of metal ions and examining the interactive effects of divalent cations in various combinations is more representative of the actual environmental problems faced by organisms than the studies of a single metal. As a result, multicomponent or competitive adsorption of metals from water solution has been investigated on various adsorbents in the last decade. The following multicomponent adsorption systems have been studied: Cu(II)-Cd(II) on chitosan cross-linked with epichlorohydrin-triphosphate [5]; Cu(II)-Zn(II) on Myriophyllum spicatum [6] Cu(II)/Pb(II)-Cd(II) on immobilized Lentinus edodes residues [7]; Zn(II)-Cd(II)-Ni(II)-Cu(II)- Cr(III) on clay-bearing mining wastes [8]; Pb(II)-Zn(II) on sulfured orange peel [9]; and Cu(II)-Pb(II) on activated carbon obtained by pyrolysis of cow bone [10]. In a study of the competitive adsorption of Pb²⁺ and Zn²⁺ onto activated carbon derived from van apple pulp, it was found that individual metal ions adsorption decreased whilst more metal ions were adsorbed in multi-metal solutions [11]. Binary and ternary sorption of Pb²⁺, Cd²⁺ and Zn²⁺ ions onto Eicchornia crassipes showed that the combined action of the metals was antagonistic and the metal sorption followed the order Pb²⁺>>Cd²⁺>Zn²⁺ [12]. The same authors also observed that Pb²⁺ ions could still be effectively removed in the presence of Zn²⁺ and Cd²⁺ ions, but the removal of Cd²⁺ and Zn²⁺ ions would be suppressed in the presence of Pb²⁺ ions. Zhu et al. [13], carried out sorption studies of Pb²⁺, Cu²⁺, and Zn²⁺ from ternary solutions onto xanthate-modified magnetic chitosan and reported that the adsorption capacities followed the order: Pb²⁺>Cu²⁺>Zn^²⁺. More recently, Sulaymon et al. [14] observed that the biosorption capacity for single metals onto activated sludge decreased by 11-51%, 53-88% and 79-94% in the binary, ternary and quaternary systems respectively at the optimum agitation speed of 600-800 rpm.

In many cases, biosorption studies have been carried out using native biomasses and products obtained from biomasses, which are generally biopolymers (polysaccharides and glycoproteins). Use of biosorbents in native forms from microbial biomasses (e.g. yeasts, microalgae, bacteria, etc.) has some drawbacks such as difficulties on separation of cells after the biosorption, mass loss during the separation, and the low mechanic resistance of the cells [15]. Several techniques exist which can make biosorbents suitable for process applications. Among these, adsorption on inert supports by preparation of biofilms; encapsulation in polymeric matrices as calcium alginate, polyacrylamide, polysulfone, and polyhydroxyetilmetacrilate; covalent linkage on supports by chemical agents; and crosslinking by chemical agents that form stable cellular aggregated have been found to be practical for biosorption [16].

The present work studies for the first time, the immobilization of H. ranunculoides root and stem biomass in calcium alginate and its potential use as a biosorbent in batch and continuous-flow systems. The perennial plant occurs abundantly at Lake Chivero in Zimbabwe where its extensive mats disrupt ecology and poses a threat to recreational uses of the waterways. The high abundance of the plant makes its possible conversion into chemically and physically stable beads to be an attractive option for removal of toxic metals from aquatic environments. The alginate-H. ranunculoides beads are investigated for their ability to remove Pb(II), Cd(II) and Cu(II) from binary solutions in batch system and from single metal solutions in continuous system.

Materials and Methods

Reagents

Sodium alginate, calcium chloride, copper sulphate, lead nitrate and cadmium chloride were purchased from Sigma Aldrich, Germany. Heavy metal stock solutions (1000 ppm) were separately prepared by dissolving analytical grade CuSO₄.5H₂O, Pb(NO₃)₂, CdCl₂.2.5H₂O in deionised water. Solutions of different concentrations were obtained by diluting the stock solution.

Biosorbent preparation

H. ranunculoides samples were collected from Lake Chivero, near Harare, Zimbabwe. Root and stem samples were washed with running tap water to remove sand and grit, rinsed with deionised water and sun dried. The samples were then ground to a fine powder using mortar and pestle and sieved through mesh of size less than or equal to 150 μm. Acid pretreatment of the biomass was then followed by agitating approximately 5 g of the biomass for every 100 ml of 1 M HCl for 12 hr. After HCl treatment the biomass was collected by filtration under vacuum suction. This was followed by washing with de-ionised water to remove excess acid. The process was repeated 5 times before drying the biomass at 40°C in an oven until constant weight.

The dried biomass powder was fixed on alginate using a method described by Vijayaraghavan and Yuan [17]. A 2% (w/v) polymer solution was prepared by dissolving 2 g of Sodium alginate for every 100 g of deionised water and then shaking for 24 hr at 300 rpm to make a homogenous polymer solvent mixture. Two grams of the biomass powder was then blended with 100 ml of the 2% (w/v) sodium alginate solution which was then extruded dropwise through a funnel into 0.1 M calcium chloride solution to make beads with a biomass loading of 2% and dimensions of 3 to 5 mm in diameter. The beads were hardened by placing them in 2% calcium chloride solution for 12 hr. Those that floated on the surface were discarded. The remaining beads were washed by agitation in deionised water at 100 rpm in 250 ml conical flasks for 30 min, discarding the solution, and then repeating the process 5 times until the conductivity the resulting filtrate approached that of the deionised water. After removing residual water with paper towels the beads were dried gradually at room temperature to constant weight and then used for further experiments.

Biosorbent characterization

Scanning electron microscopy (SEM) morphologies were done in an FEI Nova Nano SEM 230 instrument (Holand) and images were acquired at magnification of 500x, 2000x, and 10000x. IR studies were carried out using a Perkin-Elmer Spectrum 100 FT-IR over a range of 4000-4500 cm^-1. Specific surface area was determined according to the Brunaru Emmet-Teller (BET) protocol on a Micromeritics model ASAP 2000. The point of zero charge (pHpzc) of the adsorbent was determined using the method described by Ofomaja and Ho [18].

Adsorption column

The column used was a glass tube with an inner diameter=1.2 cm, outer diameter=1.83 cm, length=30 cm. An adjustable plunger was attached with a 20 μm filter of 1.1 mm thickness at the top of the column, and a stainless sieve with a mesh size of 70 (0.5 mm) was placed at the bottom. The following parameters were kept constant throughout the studies: bulk density (0.25 g/cm³), packing density (1.2 g/cm³), void fraction (0.25), exchangeable calcium (130 mg/g) and biomass loading of 5%.

Batch biosorption experiments

Biosorption experiments onto alginate-H. ranunculoides beads were conducted in 250 ml flasks containing 50 ml of single-component or binary mixtures of heavy metal solutions and 0.1 g of biomass. Initial pH was adjusted to 5.5 using 0.1 M HNO3 or 0.1 M NaOH. Binary metal mixtures [Pb(II)-Cd(II), Pb(II)-Cu(II) and Cd(II)-Cu(II)] were made for combinations of the metals in the ratios; 1:1, 1:2 and 2:1 respectively by mixing appropriate quantities of the metal solutions as reported by Mahamadi and Nharingo [12]. For example, a 1:1 Cu(II)- Pb(II) mixture was prepared by making a composite mixture comprising of 10 mg/l Cu(II) and 10 mg/l Pb(II) respectively and a 1:2 Cu(II)-Pb(II) mixture was prepared by making a composite mixture comprising of 10 mg/l Cu(II) and 20 mg/l Pb(II) respectively. The mixtures were shaken on a rotary shaker (agitation rate, 150 rpm) for a predetermined equilibrium period of 6 h. The biosorbent was then filtered through an acid-cleaned 0.45 μm Millipore filter and the concentrations of heavy metal ions in the filtrate were analysed by flame atomic absorption spectrometry. The amount of heavy metal ions uptake by the beads in each flask was determined using the mass balance equation:

(1)

where q is the adsorption amount at equilibrium (mg/g); C0 is the initial metal concentration (mg/L); Ce is the concentration remaining in solution at equilibrium (mg/l); W is biosorbent dose (g/l). Metal ion concentration in the exiting solution was determined using a Varian AA- 1275 Spectrometer.

Column biosorption tests

In all column experiments, packing density was kept at approximately value of 1.0 g/cm3 by measuring the appropriate mass of the beads for a given bed depth. Metal solutions were prepared as described for batch experiments. A flowrate of 5 ml/min was delivered by an EYELA Microtube MP-3 (Japan) peristaltic pump and controlled by a flow-meter (Dwyers) and in-line flow-pulse dampener. Metal ion concentration was determined by flame atomic absorption spectrometry.

Results and Discussion

Characterization of the adsorbent

The FT-IR spectrum of alginate-H. ranunculoides is given in Figure 1. The strong and very broad peak at 3632.47-2963.14 cm^-1 indicates the stretching vibrations of a hydrogenbonded O-H group. The peak at 1600.43 is due to the bending vibrations of an amine N-H group or aromatic C=C. The presence of an aromatic ring is further suggested by the peak occurring at 1420.07. Weak bands at 1319.88, 1263.76, and 1155.55 cm^-1 are due to a C-O stretch and suggest the presence of an acyl C-O, phenol C-O, or alkoxy C-O. The presence of an alkoxy C-O is further confirmed by the strong peak at 1031.3 cm^-1. The band centred at ca. 560 cm^-1 was associated with out-of-plane bending vibrations of C-H or O-H group [19]. The SEM images of alginate-H. ranunculoides are shown in Figure 2 and they clearly reveal the polysaccharide nature of the adsorbent having a network of organized fibers. A pHpzc value of 4.7 was obtained which shows that cations adsorption would be favourable at pH values higher than 4.7, and anions adsorption at pH values lower than pHpzc. Consequently, all adsorption experiments for metal sorption were carried out at an initial pH of 5.5. A BET surface area of 11.2 m²/g was obtained for alginate-H. ranunculoides, which is comparable to that reported for polyaniline grafted chitosan beads [20].

Monocomponent adsorption

The equilibrium of the adsorptive removal of the metals from single-component systems was described by fitting data to Langmuir model expressed in a linearized form as follows:

(2)

where q_e (mg/g) is the equilibrium sorption capacity; q_max(mg/g) is the maximum amount of metal ion per unit weight of the alginate-H. ranunculoides beads required to form a complete monolayer on the surface, and b is the equilibrium adsorption constant which is related to the affinity of the binding sites.

Batch binary sorption

The nature of competitive adsorption in two-metal systems was analysed using the Langmuir competitive model [24]:

(3)

where q_max,i, and K _L,i are physical parameters, and , C _e,i, are equilibrium concentrations in the mixture of the solutes. It can be shown that in a system where the concentrations of the solutes are sufficiently large that surface coverage is substantially complete, the unit term in Equation 1 may be neglected and the equation can be linearized as follows:

(4)

Results of the plots of C_e,1/C_e,2 q_e,1 as a function of C_e,1/C_e,2 for the sorption of Pb(II) in the presence of Cd(II), Cd(II) in the presence of Pb(II), and Cu(II) in the presence of Cd(II) are shown in Figures 5-7 and the model parameters are given in Table 1. The equilibrium sorption capacity of one metal in the presence of another, q_e, was normalized to that in the absence of the other metal, q_e, and q_e/q_e the ratios were used to probe the sorption dynamics.

The normalized sorption ratios, q_e/q_e for the sorption of one metal in the presence of another metal are shown in Table 1. For binary systems of equal metal concentrations (1:1), the ratios were all <1, indicating that the adsorption of the metals was depressed by the presence of other metals in the binary solution. These results suggest that the effects of the mixtures where antagonistic. Furthermore, the q_e/q_e ratios for the sorption of Pb(II) ions in the presence of Cu(II) and Cd(II) ions were 0.54 and 0.72, respectively suggesting that Cu(II) ions had a higher depressive effect on the sorption of Pb(II) ions by alginate-H. ranunculoides. Similar observations were made for the Pb(II)-M (2:1 and 1:2) compositions (where M=Cd(II) or Cu(II)). It is noteworthy that when the concentration of Pb(II) was doubled relative to Cd(II) or Cu(II), it became increasingly difficult for the adsorbent to sequester Cd(II) or Cu(II) ions respectively. These findings indicate that the observed selectivity was not only related to metal ionic characteristics of the species present but was also related to the concentration ratio of the ions present in solution considering the all other environmental factors were held constant. It was also observed that sorption of Cu(II) ions was more significantly depressed in the presence of Pb(II) ions than in the presence of Cd(II) ions and similar findings were made for the sorption of Cu(II) ions. It may also be interesting to note that the overall adsorption capacities of alginate-H. ranunculoides for Pb(II)-Cd(II) (9.56 + 3.32) and Pb(II)-Cu(II) (7.15 + 4.98) mg/g were higher than that of Cd(II) (6.48 mg/g) and Cu(II) (6.55 mg/g) but lower than that of Pb(II) (13.21 mg/g) in single-metal sorption. A possible explanation is that there could be partial overlap between the adsorption sites of Pb(II) and those of Cd(II) and Cu(II) ions in binary systems.

Single and binary system adsorption parameters obtained using the Langmuir Monocomponent and Langmuir Competitive Models (Table 1).

Continuous flow sorption

Breakthrough curves: The breakthrough curves for the sorption of Pb(II), Cu(II) and Cd(II) ions by alginate-H. ranunculoides are shown in Figure 8. The breakthrough times are estimated to be 20, 22, and 35 h, whilst saturation times were 75, 90 and 120 h for the adsorption of Cd(II), Cu(II) and Pb(II) ions respectively. Furthermore, it is noteworthy that the curves were steeper for the adsorption of Cu(II) and Cd(II) ions than for Pb(II) ion, which was characterized by a broadened mass transfer zone. The closeness in the breakthrough times for Cu(II) and Cd(II) is consistent with the closeness in the batch monolayer maximum adsorption capacities as shown in Table 1 and Figure 8

Dynamic adsorption models: The concentration profiles in the liquid and solid phases change in space and time during column adsorption. Therefore, from perspective of design and optimization of the column processes, it is important to predict the breakthrough curve for the effluent. Numerous fixed bed models have been formulated at different levels of complexity in the modeling of environmental sorption and biosorption research [25]. These include the mass transfer, Adams-Bohart, Thomas, Yoon Nelson, Bed Depth Service Time model, etc. In this study, Adams-Bohart and Thomas models were applied to analyze the experimental data.

Adams-Bohart model: The Adams-Bohart model is used only for the description of the initial part of the breakthrough curve, i.e. up to 10-50% of the saturation point [26]. In this model, it is assumed that the adsorption process is continuous and that equilibrium is not attained instantaneously and the most commonly used form of the model is shown in Eq. 5.

(5)

where C_o and C_t are the inlet and outlet adsorbate concentrations (mg l^-1) respectively, Z is the bed height (c_m), U_O is superficial velocity (cm min^-1), defined as the ratio of the volumetric flowrate Q (cm³ min^-1) to the cross-sectional area of the bed A (cm²), KABis the Adams-Bohart kinetic constant (l mg^-1 min^-1), N0 is the saturation concentration (mg l^-1). From this equation, values describing the characteristic operational parameters of the column can be determined from a plot of against t at a given bed depth and flow rate.

The Adams-Bohart plots obtained after applying Equation 5 to the experimental data are shown in Figure 9 and linear relationships where obtained for the relative concentration region up to 0.5, i.e. up to 50% breakthrough, for all breakthrough curves. The kinetic constant (K_AB) and saturation concentration (N₀) values were calculated from the slope and intercept of the curves, respectively, and are given in Table 2. In all cases, the values of KAB decreased with increase in initial metal concentration as well as flow rate, and increased with increase in bed depth. Since the Adams-Bohart model assumes that surface reaction is the rate-limiting step, the variation of the kinetic parameter, KAB, with operational conditions indicates that the sorption process for the metal ions was significantly controlled by mass transfer rather than surface reaction [27]. When mass transfer limitation is present, the rate constant KAB is no longer constant, but becomes a lumped rate constant containing the effects of both the intrinsic kinetics and mass transport.

Thomas model: The linearized form of the Thomas model can be described by the following expression [28,29]

(6)

where K_TH is the Thomas rate constant (ml/min mg), q₀ is the adsorption capacity (mg/g) and t stands for total flow time (min). From Eq. (6), the values of KTH and q0 can be determined from the linear plot of

against t.

The Thomas model is suitable for adsorption process where the external and internal diffusion will not be the limiting step. From Figure 10, it was shown that the Thomas model (range from 0.95 to 0.99) provided a better fitting compared to Adams-Bohart model (R2 range from 0.83 to 0.93). An inadequacy of the Adams-Bohart model arises when the isotherm under study exhibits Langmuir equilibrium characteristics which, under certain conditions cannot be approximated by a rectangular isotherm [28,29]. It can be seen from Table 3 that the value of q₀ increased with increasing bed depth and flow rate. However, the value of KTH decreased with bed depth for all the metals studied. Furthermore, the value of KTH decreased while the value of q0 increased as the initial metal concentration increased. These observations were attributed to the driving force for adsorption in the concentration difference. The lower flow rate, higher influent concentration and higher bed depth would increase the adsorption of Pb(II), Cd(II) and Cu(II) onto the H. ranunculoides column. It can also be seen in Table 3 that the trend in q₀ values obtained in column sorption, (Pb(II)>Cu(II)>Cd(II)) is similar to the trend in batch sorption experiments (Table 1). It can further been seen from comparison of batch and continuous flow data that the sorption capacities obtained followed the trend: qmax values- Pb(II) (22.25 mg/g)>Cu(II) (12.55 mg/g)>Cd(II) (11.45 mg/g); N0 values- Pb(II) (6601 mg/l)>>Cu(II) (215 mg/l)>Cd(II) (152 mg/l); and q₀ values- Pb(II) (61321.28∼125602.31 mg/g)>>Cu(II) (1235.45∼6489.71 mg/g)>Cd(II) (989.36∼3756.85 mg/g) for the batch and continuous-flow systems, respectively. These results clearly indicate the higher sorption capacity for Pb(II) compared to the other two metal ions. However, q₀ the values are significantly different from qmax and this observation was further probed by calculating values of q0 at breakthrough (q0.5 ) and saturation points (qS). As an illustration, was calculated as follows:

At breakthrough times of 20, 22 and 35 h for Cd(II), Cu(II) and Pb(II) respectively, flow rate of 10 ml/min, initial metal concentration of 10 mg/L and adsorbent mass of 35.5 g, q_0.5 values of 3.38, 3.72 and 5.92 mg/g were obtained for Cd(II), Cu(II) and Pb(II) respectively. Saturation adsorption capacities (qS) of 17.26, 21.2 and 28.24 mg/g were obtained for Cd(II), Cu(II) and Pb(II) at saturation times of 75, 90 and 120 h, respectively. It is interesting to note that these values of q_S were comparable to q_max obtained in single metal system from Langmuir model.

Conditions: concentration of metal ions 10 mg/L, bed depth 2 cm, flowrate 10 ml/min. C_o and C_t are the inlet and outlet metal ions concentrations respectively.

Conclusions

This paper presented findings from a detailed study of the adsorption of Pb(II), Cu(II) and Cd(II) ions using an abundantly available material: H. ranunculoides. Experiments were carried out in batch and continuous flow systems. Results from single-metal sorption experiments showed that the adsorption trend followed the order: Pb(II)>Cu(II)>Cd(II). In binary systems, the highest inhibition capacity was for Pb(II) ions in the presence of Cd(II) ions, whereas the lowest was for Cd(II) ions in the presence of Pb(II) ions. Furthermore, it was shown that Pb(II) ions could still be effectively removed in the presence of either Cu(II) or Cd(II) ions and the competition in sorption of the heavy metals was affected by the relative concentration of the interfering metal. Continuous flow studies revealed that the Thomas model provided a better fit to the experimental data than the Adams- Bohart model and it was concluded that mass transfer, i.e. external diffusion was not the rate limiting step during the removal of Pb(II), Cu(II) and Cd(II) by immobilised H. ranunculoides.

Acknowledgements

The International Foundation for Science, Stockholm, Sweden supported this research, through Grant Number W/4266-2 to Courtie Mahamadi