Aggregation kinetics of kaolinite-fulvic acid colloids as affected by the sorption of Cu and Pb

Aggregation kinetics of kaolinite-fulvic acid colloids, as influenced by two strongly sorbing trace metal cations, Cu(II) and Pb(II), was investigated by time-resolved dynamic laser light scattering experiments. The effects of Cu and Pb on the aggregation rate and electrophoretic mobilitywere compared with that of Ca, another major divalent metal cation which is less strongly adsorbed. Kaolinite-fulvic acid suspensions (in 0.01 M NaNO3 at pH 4 and pH 6) were spiked with solutions containing Cu, Pb, or Ca to give total divalent cation concentrations between 10(-5) and 8 x 10(-3) M. The concentration of kaolinite was varied between 25 and 200 mg L(-1), while the concentration of fulvic acid ranged from 0.15 to 1.2 mg L(-1). The mass ratio of kaolinite to fulvic acid was kept constant at 500:3 in all experiments. Relative aggregation rates, expressed as attachment efficiency alpha, were determined from linear increases in average hydrodynamic radius with time during the first 5-8 min of the aggregation experiments, always starting with a well-dispersed suspension at time zero. The corresponding slope for fast aggregation (alpha = 1) was measured for pure kaolinite suspended in 0.01 M NaNO3 at pH 4. Addition of fulvic acid to the suspensions completely inhibited kaolinite aggregation at pH 4 and pH 6. Additions of Cu, Pb, and Ca resulted in strongly increased aggregation rates of the kaolinite-fulvic acid particles. The potential of the three cations to enhance aggregation of the kaolinite-fulvic acid colloids increased in the order Ca < Cu < or = Pb. At pH 4, the relationship between particle electrophoretic mobility and aggregation rate was the same for all three divalent metal cations. In the presence of Ca, an increase in pH from 4 to 6 resulted in decreased aggregation rates. However, in the presence of Cu or Pb, the opposite trend was observed and the relationship between electrophoretic mobility and aggregation rate was different than at pH 4. The effects of Cu, Pb, and Ca on the aggregation rates of kaolinite-fulvic acid colloids are explained by the different sorption behavior of the three divalent metal cations.     

Complexation with dissolved organic matter and solubility control of heavy metals in a sandy soil

The complexation of heavy metals with dissolved organic matter (DOM) in the environment influences the solubility and mobility of these metals. In this paper, we measured the complexation of Cu, Cd, Zn, Ni, and Pb with DOM in the soil solution at pH 3.7-6.1 using a Donnan membrane technique. The results show that the DOM-complexed species is generally more significant for Cu and Pb than for Cd, Zn, and Ni. The ability of two advanced models for ion binding to humic substances, e.g., model VI and NICA-Donnan, in the simulation of metal binding to natural DOM was assessed by comparing the model predictions with the measurements. Using the default parameters of fulvic and humic acid, the predicted concentrations of free metal ions from the solution speciation calculation using the two models are mostly within 1 order of magnitude difference from the measured concentrations, except for Ni and Pb in a few samples. Furthermore, the solid-solution partitioning of the metals was simulated using a multisurface model, in which metal binding to soil organic matter, dissolved organic matter, clay, and iron hydroxides was accounted for using adsorption and cation exchange models (NICA-Donnan, Donnan, DDL, CD-MUSIC). The model estimation of the dissolved concentration of the metals is mostly within 1 order of magnitude difference from those measured except for Ni in some samples and Pb. The solubility of the metals depends mainly on the metal loading over soil sorbents, pH, and the concentration of inorganic ligands and DOM in the soil solution.     

Analysis of the effect of pH on Cu2+-fulvic acid complexation using a simple electostatic model

A simple electrostatic model was used to study the effect of pH on the binding of Cu2+ to fulvic acid in solutions containing similar amounts of dissolved organic carbon (DOC) as natural media, such as aquatic environments and soil solutions. Complexation behavior was affected by increased pH because of changes in the electrostatic interaction resulting from an increase in the negative charge on the fulvic acid molecule. Solutions of soil-extracted fulvic acid (FA), at concentrations of 25 and 35 mg L(-1), ionic strength 0.005 M, and pH 5.0, 5.5, 6.0, 6.5, 7.0, and 7.5, conditions that simulate those of natural freshwaters, were titrated with copper ion using differential pulse anodic stripping voltammetry. Assuming the formation of 1:1 complexes, the conditional binding parameters (stability constant and binding capacity) were calculated for each pH value. Use of a 1:1 electrostatic model allowed estimation of the contribution of the electrostatic effect to the ion binding reaction, at each pH value, as well as calculation of a binding constant that was not dependent on pH and which thus represented the contribution of the chemical heterogeneity. Furthermore, it was found that only a small proportion of the FA acid functional binding sites are involved in the formation of copper complexes.     

Use of diffusive gradients in thin films to measure cadmium speciation in solutions with synthetic and natural ligands: comparison with model predictions

The performance of the technique of diffusive gradients in thin films (DGT) was characterized in well-defined systems containing cadmium with chloride and nitrate ions, simple organic ligands (nitrilotriacetic acid and diglycolic acid), and Suwannee river fulvic acid for the pH range 5-8. Cd was fully labile in all Cd, Cl-, and NO3- solutions tested (I= 0.1 and 0.01 M), even atvery low Cd concentrations (10 nM), consistent with there being no binding of Cd to the diffusive gel. Diffusion coefficients of Cd-nitritotriacetic acid (NTA) and Cd-diglycolic acid (DGA) species were measured and found to be ca. 25-30% lower than the equivalent coefficient for free metal ions. These values were used to calculate concentrations of labile Cd from DGT measurements in solutions of Cd with NTA or DGA. Cd-NTA and Cd-DGA species were found to be fully DGT-labile. DGT devices that used a diffusive gel with a reduced pore size, which retarded the passage of fulvic acid species through the gel, were used to estimate the proportion of Cd complexed by fulvic acid. These results were compared with predictions of the solution speciation from models with default parameter values. ECOSAT, incorporating the NICA-Donnan model, correctly predicted the magnitude of the binding and its pH dependence, while predictions from WHAM V (with humic ion binding model V) and WHAM 6 (with humic ion binding model VI) were less satisfactory at predicting the pH dependence. Reasonable fits to the data could be obtained from WHAM 6 when the effective binding constant log K(MA) was changed from 1.6 to 1.5, the value of deltaLK1 from 2.8 to 1.0 to minimize the dependence on pH, and the value of deltaLK2 from 1.48 to 1.0 to decrease the strength of the strong bidentate and tridentate binding sites.     

Interactions of trace metals with hydrogels and filter membranes used in DET and DGT techniques

Equilibrium partitioning of trace metals between bulk solution and hydrogels/filter was studied. Under some conditions, trace metal concentrations were higher in the hydrogels or filter membranes compared to bulk solution (enrichment). In synthetic soft water, enrichment of cationic trace metals in polyacrylamide hydrogels decreased with increasing trace metal concentration. Enrichment was little affected by Ca and Mg in the concentration range typically encountered in natural freshwaters, indicating high affinity but low capacity binding of trace metals to solid structure in polyacrylamide gels. The apparent binding strength decreased in the sequence: Cu > Pb > Ni approximately to Cd approximately to Co and a low concentration of cationic Cu eliminated enrichment of weakly binding trace metal cations. The polyacrylamide gels also had an affinity for fulvic acid and/or its trace metal complexes. Enrichment of cationic Cd in agarose gel and hydrophilic polyethersulfone filter was independent of concentration (10 nM to 5 microM) but decreased with increasing Ca/ Mg concentration and ionic strength, suggesting that it is mainly due to electrostatic interactions. However, Cu and Pb were enriched even after equilibration in seawater, indicating that these metals additionally bind to sites within the agarose gel and filter. Compared to the polyacrylamide gels, agarose gel had a lower affinity for metal-fulvic complexes. Potential biases in measurements made with the diffusive equilibration in thin-films (DET) technique, identified by this work, are discussed.     

Analysis of trace metal humic acid interactions using counterion condensation theory

Recent extensions of counterion condensation theory, originally developed for well-defined linear polyelectrolytes, enable us to analyze the interaction of trace metals with humic acid. In the present model, the heterogeneity of the macromolecule is taken into account as well as the chemical binding of the considered metal ions to the humic material. Experimentally, potentiometric titrations have been performed for humic acid in solution in the presence of different environmentally important (heavy) metals (Ca, Cd, Cu, Ni, and Pb) at various metal concentrations by titrating with potassium hydroxide without additional salt. From proton release data obtained for the initial point in the titration, it was estimated that the interaction of the different metals with the humic acid in terms of binding strength increased in the order Ca < Cd approximately = Ni < Pb approximately = Cu. These results were confirmed by model analysis. Experimentally obtained apparent dissociation constants were in good agreement for the humic acid systems containing Ca, Cd, and Ni at concentrations ranging from 0 up to 0.75 x 10(-3) mol L(-1) and polymer dissociation degree from about 0.1 up to approximately 0.8. Also for the Cu/humic acid and Pb/humic acid systems, the agreement between experimental data and calculated data was satisfactory atthe lowest metal concentrations over the complete titration curve. For elevated levels of Cu and Pb, the agreement between experimental data and theoretical calculations becomes less satisfactory at low degrees of dissociation of the humic acid. This distortion of the potentiometric curves is probably due to changes in the intrinsic pK of the functional groups due to metal binding. This complex process is not included in present polyelectrolytic models.     

Strong copper-binding behavior of terrestrial humic substances in seawater

In coastal areas, strong complexation of copper generally reduces its toxicity; our ability to monitor and regulate copper as a toxin therefore depends on our understanding of the sources and sinks of the copper-binding ligands. Terrestrial humic substances (HS) are well-recognized contributors to weak ligand concentrations in aquatic systems. In this work, we show that HS are likely contributors to both stronger and weaker ligand classes controlling copper speciation in coastal areas receiving typical inputs of terrestrial organic matter. We used competitive ligand exchange adsorptive cathodic stripping voltammetry (CLE-ACSV), with the added ligands benzoylacetone and salicylaldoxime, to examine copper binding by terrestrial HS in a seawater matrix, at HS and copper concentrations typical of coastal waters. Copper titration data of 1 mg/L Suwannee River humic acid (SRHA) in seawater could be modeled using conditional stability constants of 10(12.0) and 10(10.0) and total ligand concentrations of 10.4 and 199 nM for a stronger and weaker ligand, respectively. Similar results were obtained for Suwannee River fulvic acid (SRFA). Strong copper binding by SRFA in seawater was weaker than previously reported for a freshwater at similar pH, possibly indicating effects of Ca and Mg competition or ionic strength. Nevertheless,the concentrations and binding strengths of copper ligands we observed are comparable to the range reported in previous coastal speciation studies. In addition, we show that the weaker copper ligands cause internal calibration techniques to significantly underestimate the sensitivity of ACSV in the presence of HS concentrations typical of coastal waters. To address this issue, we demonstrate the use of "overload titrations", using a high enough concentration of added ligand to outcompete all natural ligands as an alternative calibration technique for analysis of coastal samples.     

Modeling salt-dependent proton binding by organic soils with the NICA-Donnan and Stockholm humic models

Models are available for simulations of proton dissociation and cation binding by natural organic matter; two examples are the NICA-Donnan and Stockholm Humic (SHM) models. To model proton and metal binding, it is necessary to properly account for the ionic strength dependence of proton dissociation. In previous applications of the models for soils itwas assumed that the electrostatic interactions for solid-phase humic substances were the same as in solution; this assumption was recently challenged. Therefore, we reanalyzed previously published acid-base titrations of acid-washed Sphagnum peat, and we produced additional data sets for two Sphagnum peats and two Spodosol Oe horizons. For the soil suspensions, the original NICA-Donnan and SHM models, which were developed for dissolved humic substances, underestimated the observed salt dependence considerably. When a fixed Donnan volume of 1 L kg(-1) for humic substances in the solid phase was used, the NICA-Donnan model fits were much improved. Also for SHM, slight changes produced improved model fits. The models also produced acceptable simulations of the dissolved Ca, Mg, and Cd concentrations, provided that cation selectivitywas introduced. In conclusion, the proposed extensions to the NICA-Donnan and SHM models were shown to predict the salt dependence of solid-phase humic substances more satisfactorily than earlier model versions.     

Adsorption of nonionic aromatic compounds to single-walled carbon nanotubes: effects of aqueous solution chemistry

We systematically studied effects of pH, ionic strength, and presence of Cu2+ (50 mg/L) or a dissolved humic acid (HA, Fluka) (50 mg/L) on adsorption of three nonionic aromatic compounds, naphthalene, 1,3-dinitrobenzene, and 1,3,5-trinitrobenzene to single-walled carbon nanotubes. Presence of Cu2+ or variance in the ionic strength between 0.02 and 0.1 M (NaNO3) only slightly affected adsorption affinities. Presence of HA reduced adsorption of the three compounds by 29-57% for CNTs, as measured by change in distribution coefficient(Kd), and by 80-95% for graphite. In contrast to nonporous graphite, whose surface area was completely accessible in adsorption, CNTs formed aggregates with microporous interstices in aqueous solution, which blocked large HA molecules from competing with the surface area. Changing the pH from 2 to 11 did not affect adsorption of naphthalene, while it increased adsorption of 1,3-dinitrobenzene and 1,3,5-trinitrobenzene by 2-3 times. Increasing pH apparently facilitated deprotonation of the acidic functional groups (-COOH, -OH) of CNTs, which promoted the pi-electron-donor ability of the graphene surface, therefore enhancing pi-pi electron-donor-acceptor (EDA) interactions of the two nitroaromatics (pi-electron acceptors).     

Infrared spectroscopic evidence supporting heterogeneous site binding models for humic substances

Infrared spectroscopy was used to corroborate predictions made by newly developed heterogeneous site binding models for humic substances. Experimental conditions to acquire the spectra of soil humic substances (humic and fulvic acid and a polysaccharide fraction) in an aqueous state using horizontal attenuated total reflectance Fourier transform infrared spectroscopy (HATR-FTIR) were established. Elimination of the water spectrum from that of the sample was achieved by spectral subtraction of the water peak at 2020 cm(-1). A KSCN internal standard with an absorption band at 2067 cm(-1) was used to verify the efficacy of the subtraction procedure. Spectral artifacts produced by the water spectrum subtraction and from contaminants within the humic materials have been identified. Three fulvic and one humic acid solution were examined in solutions of varying pH. Results show that the observed proportion of ionized carboxylate in relation to pH is consistent with models that assume electrostatic effects and a continuous distribution of proton association constants (log KH). The spectroscopic data were in accordance with calculations made using the generic humic and fulvic acid NICA-Donnan model parameters.     

Voltammetric analysis of heterogeneity in metal ion binding by humics

The complexation of Cd, Pb, and Cu by fulvic acids at a fixed pH and ionic strength is studied by means of different voltammetric techniques at any metal-to-ligand ratio. When using Reverse Pulse Polarography (RPP) the complex species are electrochemically labile and not subject to significant electrodic adsorption. RPP titrations of fulvic acid with metal ions are interpreted on the basis of a recently proposed analytical expression for limiting currents valid for fully labile heterogeneous complexation. The voltammetric data are transformed into the corresponding binding curve, i.e., the fraction of occupied sites vs free metal concentration. Finally, the competition between metal ions and protons in their interaction with the fulvic binding sites as well as the concomitant polyelectrolytic effects are analyzed in terms of the NICCA-Donnan model. The results show that voltammetric techniques can be applied to the studies of heterogeneous complex systems in a broad range of metal-to-ligand ratios.     

Solid-solution partitioning of organic matter in soils as influenced by an increase in pH or Ca concentration

Organic matter is an important component of soil with regard to the binding of contaminants. Hence, the partitioning of organic matter influences the partitioning of soil contaminants. The partitioning of organic matter is, among other factors, influenced by the ionic composition and ionic strength of the soil solution. This study focuses on the behavior of organic matter after a change in the ionic composition of the soil solution, particularly in Ca concentration and pH. Different amounts of Ca(NO3)2 and NaOH were added to soil suspensions. The dissolved organic carbon (DOC) concentration increased with increasing pH (addition of NaOH), whereas an increase in Ca (addition of Ca(NO3)2) had the opposite effect. A stronger increase in DOC was observed if a single dose of NaOH was added, compared to a gradual addition of the same amount of NaOH. Cation binding by organic matter in the supernatant was calculated using the NICA-Donnan model. The log DOC concentration appeared to be correlated to the Donnan potential, calculated under the assumption that all DOC equals humic acid. This correlation was found for all eight neutral to acidic soils used in this study, although the slopes and elevations of the regression lines varied. The slope varied by a factor of 2 and the elevation appeared to be strongly influenced by the DOC concentration in the untreated soils, which is related to the total organic matter in the soil. Finally, we predicted the Donnan potential on the basis of an extraction of untreated soil with 0.03 M NaNO3, and the total additions of Ca(NO3)2 and NaOH. Comparison of these predictions with speciation calculations in solution showed a good correlation, indicating that a combination of one batch experiment and the presented calculation procedure can provide good estimations of DOC concentrations after addition of chemicals.     

Impact of natural organic matter on floc size and structure effects in membrane filtration

Hematite (10 mg of Fe/L) floc-humic acid assemblages have been formed at pH 4 either by first aggregating hematite particles with salt (100 mM KCl) and then adding humic acid (salt-particle-organic or SPO assemblages) or by suspending the hematite particles in humic acid solutions and then adding salt to induce aggregation (organic-particle-salt or OPS assemblages). The behavior of these assemblages upon deposition on microfiltration (MF) membranes has then been investigated. In the OPS case, the fractal dimension (dF) of the assemblages formed varied dramatically depending upon the extent of charge neutralization by added fulvic acid with dF values typical of diffusion-limited cluster aggregates at low (0.1-0.2 mg/L) humic acid concentrations and dF values typical of reaction-limited cluster aggregates either in the absence of humic acid or concentrations greater than 0.4-0.6 mg/L. In the SPO case, dF values on the order of 2.1 were initially observed and were found to decrease to around 1.8-1.9 for humic acid concentrations greater than 0.6-0.8 mg/L. OPS assemblages with low fractal dimensions were found to be highly compressible once deposited on MF membranes with significantly higher specific cake resistances than was the case for SPO assemblages at transmembrane pressures of 50 kPa and above. These results highlight the importance of both the choice of coagulant (e.g., preformed vs formed in situ) and the transmembrane pressure to which a membrane filtration process might be allowed to rise prior to removal of the fouling layer.     

Modeling complexometric titrations of natural water samples

Complexometric titrations are the primary source of metal speciation data for aquatic systems, yet their interpretation in waters containing humic and fulvic acids remains problematic. In particular, the accuracy of inferred ambient free metal ion concentrations and parameters quantifying metal complexation by natural ligands has been challenged because of the difficulties inherent in calibrating common analytical methods and in modeling the diverse array of ligands present. This work tests and applies a new method of modeling titration data that combines calibration of analytical sensitivity (S) and estimation of concentrations and stability constants for discrete natural ligand classes ([Li]T and Ki) into a single step using nonlinear regression and a new analytical solution to the one-metal/two-ligand equilibrium problem. When applied to jointly model data from multiple titrations conducted at different analytical windows, it yields accurate estimates of S, [Li]T, Ki, and [Cu2+] plus Monte Carlo-based estimates of the uncertainty in [Cu2+]. Jointly modeling titration data at low-and high-analytical windows leads to an efficient adaptation of the recently proposed "overload" approach to calibrating ACSV/CLE measurements. Application of the method to published data sets yields model results with greater accuracy and precision than originally obtained. The discrete ligand-class model is also re-parametrized, using humic and fulvic acids, L1 class (K1 = 10(13) M(-1)), and strong ligands (L(S)) with K(S) > K1 as "natural components". This approach suggests that Cu complexation in NW Mediterranean Sea water can be well represented as 0.8 +/- 0.3/0.2 mg humic equiv/L, 13 +/- 1 nM L1, and 2.5 +/- 0.1 nM L(S) with [CU]T = 3 nM. In coastal seawater from Narragansett Bay, RI, Cu speciation can be modeled as 0.6 +/- 0.1 mg humic equiv/L and 22 +/- 1 nM L1 or approximately 12 nM L1 and approximately 9 nM L(S), with [CU]T = 13 nM. In both waters, the large excess (approximately 10 nM) of high-affinity, Cu-binding ligands over [CU]T results in low equilibrium [Cu2+] of 10(-14.5 +/- 0.2) M and 10(-13.3 +/- 0.4) M, respectively.     

In-situ speciation of Ni and Zn in freshwaters: comparison between DGT measurements and speciation models

The technique of DGT (diffusive gradients in thin films) was used for the first time to measure in situ the distribution of Zn and Ni between inorganic species and complexes with fulvic and humic acids in natural waters. With DGT, metals are bound to a resin embedded in a layer of hydrogel after diffusive transport through an adjacent layer of hydrogel. The metal concentrations in the waters can be quantified using simple diffusion equations. By using devices with hydrogels of different pore size, large and small complex species were discriminated. Inorganic species diffuse freely through all gels, but larger organic complexes with fulvic and humic acids diffuse less freely through more restricted gels (gels with smaller pore size). Systematic differences between DGT devices containing gels of different pore size were obtained. Their calibration for the diffusion of fulvic and humic complexes allowed calculation of the concentrations of labile inorganic (Zn, 34.6 +/- 2.5 nM; Ni, 23.5 +/- 0.9 nM) and labile organic (Zn, 43.1 +/- 2.9 nM; Ni, 11.2 +/- 0.7 nM) complexes. The concentration of Zn measured by anodic stripping voltammetry in samples returned to the laboratory lay between the DGT-measured inorganic concentration and the total dissolved concentration, consistent with partial measurement of organic complexes of Zn. The speciation model WHAM successfully predicted the species distribution of Ni, Zn, and Cu, provided that competitive binding by Fe(III) was considered. Using the speciation models WHAM and ECOSAT, free ion activities of Zn and Ni were calculated from (1) the total inorganic species measured by DGT and (2) the total dissolved species and dissolved organic carbon. The calculations confirmed the good model predictions of metal-humic binding but highlighted problems with default databases used for the speciation of inorganic components.     

Comparison of Cd(II), Cu(II), and Pb(II) biouptake by green algae in the presence of humic acid

The present study examines the role of humic acid, as a representative of dissolved organic matter, in Cd(II), Cu(II), and Pb(II) speciation and biouptake by green microalgae. Cellular and intracellular metal fractions were compared in the presence of citric and humic acids. The results demonstrated that Cd and Cu uptake in the presence of 10 mg L(-1) humic acid was consistent with that predicted from measured free metal concentrations, while Pb biouptake was higher. By comparing Cd, Cu, and Pb cellular concentrations in the absence and presence of humic acid, it was found that the influence of the increased negative algal surface charge, resulting from humic acid adsorption, on cellular metal was negligible. Moreover, the experimental results for all three metals were in good agreement with the ternary complex hypothesis. Given that metal has much higher affinity with algal sites than humic acid adsorbed to algae, the contribution of the ternary complex to metal bioavailability was negligible in the case of Cd (II) and Cu (II). In contrast, the ternary complex contributed to over 90% of total cellular metal for Pb(II), due to the comparable affinity of Pb to algal sites in comparison with humic acid adsorbed to algae. Therefore, the extension of the biotic ligand model by including the formation of the ternary complex between the metal, humic acid, and algal surface would help to avoid underestimation of Pb biouptake in the presence of humic substances by green algae Chlorella kesslerii.     

Copper(II) Binding by Dissolved Organic Matter: Importance of the Copper-to-Dissolved Organic Matter Ratio and Implications for the Biotic Ligand Model

The ratio of copper to dissolved organic matter (DOM) is known to affect the strength of copper binding by DOM, but previous methods to determine the Cu(2+)-DOM binding strength have generally not measured binding constants over the same Cu:DOM ratios. In this study, we used a competitive ligand exchange-solid-phase extraction (CLE-SPE) method to determine conditional stability constants for Cu(2+)-DOM binding at pH 6.6 and 0.01 M ionic strength over a range of Cu:DOM ratios that bridge the detection windows of copper-ion-selective electrode and voltammetry measurements. As the Cu:DOM ratio increased from 0.0005 to 0.1 mg of Cu/mg of DOM, the measured conditional binding constant ((c)K(CuDOM)) decreased from 10(11.5) to 10(5.6) M(-1). A comparison of the binding constants measured by CLE-SPE with those measured by copper-ion-selective electrode and voltammetry demonstrates that the Cu:DOM ratio is an important factor controlling Cu(2+)-DOM binding strength even for DOM isolates of different types and different sources and for whole water samples. The results were modeled with Visual MINTEQ and compared to results from the biotic ligand model (BLM). The BLM was found to over-estimate Cu(2+) at low total copper concentrations and under-estimate Cu(2+) at high total copper concentrations.     

Combining spectroscopic and potentiometric approaches to characterize competitive binding to humic substances

In an area that contains high concentrations of natural organic matter, it is expected that it plays an important role on the behavior of rare earth elements (REE), like europium, and of trivalent actinides. Competitive interactions with H+, inorganic species, major cations, e.g. Ca(II) or Mg(II), could influence these metals transport and bioavailability. Competitive experiments between cations, which can bind differently to humic substances and Eu3+, will bring an improved understanding of the competitive mechanisms. The aim of this study is to acquire data for Eu(III)/Cu(II) and Eu(III)/Ca(II) competitive binding to a sedimentary originated humic acid (Gorleben, Germany). The NICA-Donnan parameters for Ca2+, Cu2+, and Eu3+ obtained from competitive binding experiments using Ca2+ or Cu2+ ion selective electrodes were used to model time-resolved laser fluorescence spectroscopy (TRLFS) measurements. Eu3+ and CU2+ are in direct competition for the same type of sites, whereas Ca2+ has an indirect influence through electrostatic binding.