Evidence-based logic chains demonstrate multiple impacts of trace metals on ecosystem services.

Trace metals can have far-reaching ecosystem impacts. In this study, we develop consistent and evidence-based logic chains to demonstrate the wider effects of trace metal contamination on a suite of ecosystem services. They demonstrate knock-on effects from an initial receptor that is sensitive to metal toxicity, along a cascade of impact, to final ecosystem services via alterations to multiple ecosystem processes. We developed logic chains to highlight two aspects of metal toxicity: for impacts of copper pollution in soil ecosystems, and for impacts of mercury in freshwaters. Each link of the chains is supported by published evidence, with an indication of the strength of the supporting science. Copper pollution to soils (134 unique chains) showed a complex network of pathways originating from direct effects on a range of invertebrate and microbial taxa and plants. In contrast, mercury pollution on freshwaters (63 unique chains) shows pathways that broadly follow the food web of this habitat, reflecting the potential for mercury bioaccumulation. Despite different pathways, there is considerable overlap in the final ecosystem services impacted by both of these metals and in both ecosystems. These included reduced human-use impacts (food, fishing), reduced human non-use impacts (amenity value) and positive or negative alterations to climate regulation (impacts on carbon sequestration). Other final ecosystem goods impacted include reduced crop production, animal production, flood regulation, drinking water quality and soil purification. Taking an ecosystem services approach demonstrates that consideration of only the direct effects of metal contamination of soils and water will considerably underestimate the total impacts of these pollutants. Construction of logic chains, evidenced by published literature, allows a robust assessment of potential impacts indicating primary, secondary and tertiary effects.

The effect of advanced treatment of sewage effluents on metal speciation and (bio)availability.

The bioavailability of metals can be strongly influenced by dissolved organic carbon (DOC). Wastewater treatment effluents add considerable quantities of DOC and metals to receiving waters, and as effluent controls become more stringent advanced effluent treatments may be needed. We assessed the effects of two types of advanced treatment processes on metal availability in wastewater effluents. Trace metal availability was assessed using diffuse gradients in thin films and predicted through speciation modelling. The results show little difference in metal availability post-advanced treatment. EDTA-like compounds are important metal complexants in the effluents.

WHAM-FTOXß - An aquatic toxicity model based on intrinsic metal toxic potency and intrinsic species sensitivity.

We developed a model that quantifies aquatic cationic toxicity by a combination of the intrinsic toxicities of metals and protons and the intrinsic sensitivities of the test species. It is based on the WHAM-FTOX model, which combines the calculated binding of cations by the organism with toxicity coefficients (αH, αM) to estimate the variable FTOX, a measure of toxic effect; the key parameter αM,max (applying at infinite time) depends upon both the metal and the test species. In our new model, WHAM-FTOXß, values of αM,max are given by the product αM* × ß, where αM* has a single value for each metal, and ß a single value for each species. To parameterise WHAM-FTOXß, we assembled a set of 2182 estimates of αM,max obtained by applying the basic model to laboratory toxicity data for 76 different test species, covering 15 different metals, and including results for metal mixtures. Then we fitted the log10αM,max values with αM* and ß values (a total of 91 parameters). The resulting model accounted for 72% of the variance in log10αM,max. The values of αM* increased markedly as the chemical character of the metal changed from hard (average αM* = 4.4) to intermediate (average αM* = 25) to soft (average αM* = 560). The values of log10ß were normally distributed, with a 5-95 percentile range of -0.73 to +0.56, corresponding to ß values of 0.18 to 3.62. The WHAM-FTOXß model entails the assumption that test species exhibit common relative sensitivity, i.e. the ratio αM,max / αM* is constant across all metals. This was tested with data from studies in which the toxic responses of a single organism towards two or more metals had been measured (179 examples for the most-tested metals Ni, Cu, Zn, Ag, Cd, Pb), and statistically-significant (p < 0.003) results were obtained.

Soil properties influence the toxicity and availability of Zn from ZnO nanoparticles to earthworms.

To develop models that support site-specific risk assessment for nanoparticles (NPs), a better understanding of how NP transformation processes, bioavailability and toxicity are influenced by soil properties is needed. In this study, the influence of differing soil properties on the bioavailability and toxicity of zinc oxide (ZnO) NPs and ionic Zn to the earthworm Eisenia fetida was investigated. Earthworms were exposed to ZnO_NPs and ionic Zn, between 100 and 4400 mg Zn/kg, in four different natural soils (organic matter content: 1.8-16.7%, soil pH: 5.4-8.3, representing sandy loam to calcareous soils). Survival and reproduction were assessed after 28 and 56 days, respectively. Zn concentrations in soil pore waters were measured while labile concentrations of Zn were measured using an in-situ dynamic speciation technique (diffusive gradient in thin films, DGT). Earthworm Zn tissue concentrations were also measured. Soil properties influenced earthworm reproduction between soil controls, with highest reproductive output in soils with pH values of 6-7. Toxicity was also influenced by soil properties, with EC50s based on total Zn in soil ranging from 694 to >2200 mg Zn/kg for ZnO_NP and 277-734 mg Zn/kg for ionic Zn. Soil pore water and DGT measurements showed good agreement in the relative amount of Zn extracted across the four soils. Earthworms exposed to ZnO_NPs survived higher Zn concentrations in the soils and had higher tissue concentrations compared with ionic Zn exposures, particularly in the high organic content calcareous soil. These higher tissue concentrations in ZnO_NP exposed earthworm could have consequences for the persistence and trophic mobility of Zn in terrestrial systems and need to be further investigated to elucidate if there any longer-term risks associated with sustained input of ZnO_NP to soil.

The use of WHAM-FTOX, parameterized with laboratory data, to simulate zooplankton species richness in acid- and metal- contaminated lakes.

The WHAM-FTOX model quantifies cation toxicity towards freshwater organisms, assuming an additive toxic response to the amounts of protons and metals accumulated by an organism. We combined a parameterization of the model, using data from multi-species laboratory toxicity tests, with a fitted field species sensitivity distribution, to simulate the species richness (nsp) of crustacean zooplankton in acid- and metal-contaminated lakes near Sudbury, Ontario over several decades, and also in reference (uncontaminated) lakes. A good description of variation in toxic response among the zooplankton species was achieved with a log-normal distribution of a new parameter, ß, which characterizes an organism's intrinsic sensitivity towards toxic cations; the greater is ß, the more sensitive is the species. The use of ß assumes that while species vary in their sensitivity, the relative toxicities of different metals are the same for each species (common relative sensitivity). Unbiased agreements between simulated and observed nsp were obtained with a high correlation (r2 = 0.81, p < 0.0001, n = 217). Variations in zooplankton species richness in the Sudbury lakes are calculated to be dominated by toxic responses to H, Al, Cu and Ni, with a small contribution from Zn, and negligible effects of Cd, Hg and Pb. According to the model, some of the Sudbury lakes were affected predominantly by acidification (H and Al), while others were most influenced by toxic heavy metals (Ni, Cu, Zn); for lakes in the latter category, the relative importance of heavy metals, compared to H and Al, has increased over time. The results suggest that, if common relative sensitivity operates, nsp can be modelled on the basis of a single set of parameters characterizing the average toxic effects of different cations, together with a species sensitivity distribution.

Developing understanding of the fate and behaviour of silver in fresh waters and waste waters.

The Windermere Humic Aqueous Model (WHAM) is often used for risk assessment of metals; WHAM can be used to estimate the potential bioavailability of dissolved metals, where metals complexed to dissolved organic matter (DOM) are expected to be less toxic than ionic forms. Silver is a potential metal of concern but WHAM has not been rigorously tested against experimental measurements. This study compares WHAM predictions to measured ionic silver during fixed pH (4, 8 or 10) argentometric titrations of DOM from diverse origins. There were almost two orders of magnitude variation in free silver between sources but, within model uncertainty, WHAM captured this variability. This agreement, between measurements and models, suggests that WHAM is an appropriate tool for silver risk assessment in surface receiving waters when DOM is predominantly in the form of humic/fulvic acids. In sewage samples WHAM dramatically underestimated silver binding by approximately 3 orders of magnitude. Simulations with additional specific strong silver binding sulphide-like binding sites could explain Ag binding at low loadings, but not at higher loadings. This suggests the presence of additional intermediate strength binding sites. These additional ligands would represent components of the raw sewage largely absent in natural waters unimpacted by sewage effluents. A revised empirical model was proposed to account for these sewage-specific binding sites. Further, it is suspected that as sewage organic matter is degraded, either by natural attenuation or by engineered treatment, that sewage organic matter will degrade to a form more readily modelled by WHAM; i.e., humic-like substances. These ageing experiments were performed starting from raw sewage, and the material did in fact become more humic-like, but even after 30 days of aerobic incubation still showed greater Ag+ binding than WHAM predictions. In these incubation experiments it was found that silver (up to 1000 µg/L) had minimal impact on ammonia oxidation kinetics.

Kinetics of uranium(VI) lability and solubility in aerobic soils.

Uranium may pose a hazard to ecosystems and human health due to its chemotoxic and radiotoxic properties. The long half-life of many U isotopes and their ability to migrate raise concerns over disposal of radioactive wastes. This work examines the long-term U bioavailability in aerobic soils following direct deposition or transport to the surface and addresses two questions: (i) to what extent do soil properties control the kinetics of U speciation changes in soils and (ii) over what experimental timescales must U reaction kinetics be measured to reliably predict long-term of impact in the terrestrial environment? Soil microcosms spiked with soluble uranyl were incubated for 1.7 years. Changes in UVI fractionation were periodically monitored by soil extractions and isotopic dilution techniques, shedding light on the binding strength of uranyl onto the solid phase. Uranyl sorption was rapid and strongly buffered by soil Fe oxides, but UVI remained reversibly held and geochemically reactive. The pool of uranyl species able to replenish the soil solution through several equilibrium reactions is substantially larger than might be anticipated from typical chemical extractions and remarkably similar across different soils despite contrasting soil properties. Modelled kinetic parameters indicate that labile UVI declines very slowly, suggesting that the processes and transformations transferring uranyl to an intractable sink progress at a slow rate regardless of soil characteristics. This is of relevance in the context of radioecological assessments, given that soil solution is the key reservoir for plant uptake.

Systematic analysis of freshwater metal toxicity with WHAM-FTOX.

We used the WHAM chemical speciation model and the WHAM-FTOX toxicity model to analyse the published results of laboratory toxicity experiments covering 52 different freshwater biological test species and 24 different metals, a total of 2037 determinations of EC50 with accompanying data on solution composition. The key extracted parameter was αM, the parameter in WHAM-FTOX that characterises the toxic potency of a metal on the basis of its estimated metabolically active body burden. For 16 data sets applying to metal-test species pairs with appreciable variations in solution composition, values of EC50 back-calculated from averaged values of αM showed significantly (p < 0.001) less deviation from the measured EC50 values than did the simple average EC50, confirming that the modelling calculations could account for some of the dependence of toxicity on chemical speciation. Data for different exposure times permitted a simple parameterisation of temporal effects, enabling values of αM,max (values at infinite exposure time) to be obtained, and the effects of different exposure times to be factored out for further analysis. Comparison of averaged values of αM,max for different metals showed little difference among major taxa (invertebrates, plants, and vertebrates). For Cd, Cu, Ni and Zn (the four metals with most data) there were significant differences among αM,max values for different species, but within-species variabilities were greater. Reasonably similar species sensitivity distributions of standardised αM,max applied to Cd, Cu, Ni and Zn. The average values, over all species, of αM,max increased in the order Al < lanthanides < Zn ∼ UO2 < Ni ∼ Cu < Pb < Cd < Ag. Considering all the αM,max values, there was a strong dependence (r2 = 0.56, p < 0.001) on Pearson's hardness-softness categories, and a slightly stronger relationship (r2 = 0.59) if ionic radius was included in the statistical model, indicating that softer, larger cations are the most effective toxicants.

Simulating the long-term chemistry of an upland UK catchment: heavy metals.

CHUM-AM was used to investigate the behaviours of atmospherically-deposited heavy metals (Ni, Cu, Zn, Cd and Pb) in three moorland sub-catchments in Cumbria UK. The principal processes controlling cationic metals are competitive partitioning to soil organic matter, chemical interactions in solution, and chemical weathering. Metal deposition histories were generated by combining measured data for the last 30 years with local lake sediment records. For Ni, Cu, Zn and Cd, default parameters for the interactions with organic matter provided reasonable agreement between simulated and observed present-day soil metal pools and average streamwater concentrations. However, for Pb, the soil binding affinity in the model had to be increased to match the observations. Simulations suggest that weakly-sorbing metals (Ni, Zn, Cd) will respond on timescales of decades to centuries to changes in metal inputs or acidification status. More strongly-sorbing metals (Cu, Pb) will respond over centuries to millennia.

Simulating the long-term chemistry of an upland UK catchment: major solutes and acidification.

CHUM-AM was used to investigate changes in soil and water chemical variables in four moorland sub-catchments in Cumbria UK, to which non-marine S deposition has declined by 65% since the 1970s. The principal processes represented in the model comprise N and S uptake and release, water movements, the binding of cations by soil organic matter, chemical interactions in solution, and chemical weathering. CHUM-AM reproduced reasonably well the current soil pH and pools of N and S, and changes in streamwater chemistry over the period 1970-2000, notably decreases in the concentrations of alkaline earth cations and sulphate, and increases in pH. The model also predicts streamwater pH-Al relationships in agreement with observations. Predictive calculations suggest that constant atmospheric deposition of N at present rates will lead to N saturation and re-acidification, whereas a 50% reduction in N would stabilise soil and streamwater pH at about the present levels.

Using isotopic dilution to assess chemical extraction of labile Ni, Cu, Zn, Cd and Pb in soils.

Chemical extractants used to measure labile soil metal must ideally select for and solubilise the labile fraction, with minimal solubilisation of non-labile metal. We assessed four extractants (0.43 M HNO3, 0.43 M CH3COOH, 0.05 M Na2H2EDTA and 1 M CaCl2) against these requirements. For soils contaminated by contrasting sources, we compared isotopically exchangeable Ni, Cu, Zn, Cd and Pb (EValue, mg kg(-1)), with the concentrations of metal solubilised by the chemical extractants (MExt, mg kg(-1)). Crucially, we also determined isotopically exchangeable metal in the soil-extractant systems (EExt, mg kg(-1)). Thus 'EExt - EValue' quantifies the concentration of mobilised non-labile metal, while 'EExt - MExt' represents adsorbed labile metal in the presence of the extractant. Extraction with CaCl2 consistently underestimated EValue for Ni, Cu, Zn and Pb, while providing a reasonable estimate of EValue for Cd. In contrast, extraction with HNO3 both consistently mobilised non-labile metal and overestimated the EValue. Extraction with CH3COOH appeared to provide a good estimate of EValue for Cd; however, this was the net outcome of incomplete solubilisation of labile metal, and concurrent mobilisation of non-labile metal by the extractant (MExtEValue). The Na2H2EDTA extractant mobilised some non-labile metal in three of the four soils, but consistently solubilised the entire labile fraction for all soil-metal combinations (MExt ≈ EExt). Comparison of EValue, MExt and EExt provides a rigorous means of assessing the underlying action of soil chemical extraction methods and could be used to refine long-standing soil extraction methodologies.

Potentially toxic metals in ombrotrophic peat along a 400 km English-Scottish transect.

Four samples of ombrotrophic peat were collected from each of 10 upland locations in a transect from the southern Pennines to the Highland Boundary Fault, a total distance of ca. 400 km. Bulk compositions and other properties were determined. Total contents of Al and heavy metals (Ni, Cu, Zn, Cd, Pb) were determined following digestion with hydrofluoric acid, and concentrations of metals extractable with dilute nitric acid were also measured. Supernatants obtained from aqueous extractions of the peat samples were analysed for pH, major cations and anions, dissolved organic carbon and dissolved metals, and concentrations of free metal ions (Al(3+), Ni(2+), etc.) were estimated by applying a chemical speciation model. Both total and HNO(3)-extractable metal concentrations varied along the transect, the highest values being found at locations close to industrial and former mining areas. The HNO(3)-extractable soil metal contents of Ni, Cu and Cd were appreciably lower than lowest-observed-effect-concentrations (LOEC) for toxicity towards microorganisms in acid, organic rich soils. However, the contents of Zn at two locations, and of Pb at five locations exceeded LOECs, suggesting that they may be exerting toxic effects in the peats. Soil solution concentrations of free heavy metal ions (Cu(2+), Zn(2+), Cd(2+), Pb(2+)) were substantially lower than LOECs for toxicity towards vascular plants, whereas concentrations of Al(3+) were near to toxic levels at two locations.

The solid-solution partitioning of heavy metals (Cu, Zn, Cd, Pb) in upland soils of England and Wales.

Ninety-eight surface soils were sampled from the uplands of England and Wales, and analysed for loss-on-ignition (LOI), and total and dissolved base cations, Al, Fe, and trace heavy metals (Cu, Zn, Cd, Pb). The samples covered wide ranges of pH (3.4-8.3) and LOI (9-98%). Soil metal contents measured by extraction with 0.43 mol l-1 HNO3 and 0.1 mol l-1 EDTA were very similar, and generally lower than values obtained by extraction with a mixture of concentrated nitric and perchloric acids. Total heavy metal concentrations in soil solution depend positively upon soil metal content and [DOC], and negatively upon pH and LOI, values of r2 ranging from 0.39 (Cu) to 0.81 (Pb). Stronger correlations (r2=0.76-0.95) were obtained by multiple regression analysis involving free metal ion (Cu2+, Zn2+, Cd2+, Pb2+) concentrations calculated with the equilibrium speciation model WHAM/Model VI. The free metal ion concentrations depend positively upon MHNO3 and negatively upon pH and LOI. The data were also analysed by using WHAM/Model VI to describe solid-solution interactions as well as solution speciation; this involved calibrating each soil sample by adjusting the content of "active" humic matter to match the observed soil pH. The calibrated model provided fair predictions of total heavy metal concentrations in soil solution, and predicted free metal ion concentrations were in reasonable agreement with the values obtained from solution-only speciation calculations.

Reversal of acidification in tributaries of the River Duddon (English Lake District) between 1970 and 1998.

Long-term changes in stream water chemistry in the upper Duddon catchment (southwest Lake District, UK) were investigated. Ten streams were sampled and analysed weekly during 1998, and the results compared with data for the early 1970s and 1986. The waters exhibited a range of pH, average values for 1998 being 5.04-7.04. For all the streams, the average pH in 1998 was greater than that during 1971-73. Statistical analysis was carried out, using the 1970s data to estimate the magnitude of inter-annual variation, and taking discharge into account on the basis of antecedent rainfall. The results showed that for two of the streams the pH increase was significant at the 2.5% level, while for a further three it was significant at the 10% level. Comparison of the 1998 concentrations of nitrate and non-marine sulphate with data obtained for five streams in 1973-74 showed that average nitrate concentration had increased from 11 to 20 microeq dm(-3) while that of non-marine sulphate had decreased from 94 to 50 microeq dm(-3). For four of the streams, comparisons were also made between the 1998 data and those for 1986. In three cases, pH in 1998 was generally higher, and Al generally lower, than the values for 1986, but in the fourth case little difference was evident. The present results support observations for five nearby standing waters, strengthening the evidence for a general reversal of acidification in the southwest part of Lake District, due to a decline in the deposition of pollutant sulphur.

Metal mixture toxicity to aquatic biota in laboratory experiments: application of the WHAM-FTOX model.

The WHAM-FTOX model describes the combined toxic effects of protons and metal cations towards aquatic organisms through the toxicity function (FTOX), a linear combination of the products of organism-bound cation and a toxic potency coefficient (αi) for each cation. Organism-bound, metabolically-active, cation is quantified by the proxy variable, amount bound by humic acid (HA), as predicted by the WHAM chemical speciation model. We compared published measured accumulations of metals by living organisms (bacteria, algae, invertebrates) in different solutions, with WHAM predictions of metal binding to humic acid in the same solutions. After adjustment for differences in binding site density, the predictions were in reasonable line with observations (for logarithmic variables, r(2)=0.89, root mean squared deviation=0.44), supporting the use of HA binding as a proxy. Calculated loadings of H(+), Al, Cu, Zn, Cd, Pb and UO2 were used to fit observed toxic effects in 11 published mixture toxicity experiments involving bacteria, macrophytes, invertebrates and fish. Overall, WHAM-FTOX gave slightly better fits than a conventional additive model based on solution concentrations. From the derived values of αi, the toxicity of bound cations can tentatively be ranked in the order: H

Critical Limits for Hg(II) in soils, derived from chronic toxicity data.

Published chronic toxicity data for Hg(II) added to soils were assembled and evaluated to produce a data set comprising 52 chronic end-points, five each for plants and invertebrates and 42 for microbes. With end-points expressed in terms of added soil Hg(II) contents, Critical Limits were derived from the 5th percentiles of species sensitivity distributions, values of 0.13 microg(g soil)(-1) and 3.3 microg(g soil organic matter)(-1) being obtained. The latter value exceeds the currently recommended Critical Limit, used to determine Hg(II) Critical Loads in Europe, of 0.5 microg(g soil organic matter)(-1). We also applied the WHAM/Model VI chemical speciation model to estimate concentrations of Hg(2+) in soil solution, and derived an approximate Critical Limit Function (CLF) that includes pH; log [Hg(2+)](crit)=-2.15 pH -17.10. Because they take soil properties into account, the soil organic matter-based limit and the CLF provide the best assessment of toxic threat for different soils. For differing representative soils, each predicts a range of up to 100-fold in the dry weight-based content of mercury that corresponds to the Critical Limit.

Metal accumulation by stream bryophytes, related to chemical speciation.

Metal accumulation by aquatic bryophytes was investigated using data for headwater streams of differing chemistry. The Windermere Humic Aqueous Model (WHAM) was applied to calculate chemical speciation, including competitive proton and metal interactions with external binding sites on the plants. The speciation modelling approach gives smaller deviations between observed and predicted bryophyte contents of Cu, Zn, Cd and Pb than regressions based on total filtered metal concentrations. If all four metals, and Ni, are considered together, the WHAM predictions are superior at the 1% level. Optimised constants for bryophyte binding by the trace metals are similar to those for humic substances and simple carboxylate ligands. Bryophyte contents of Na, Mg and Ca are approximately explained by binding at external sites, while most of the K is intracellular. Oxide phases account for some of the Al, and most of the Mn, Fe and Co.

Solid-solution metal partitioning in the Humber rivers: application of WHAM and SCAMP

The solid-solution partitioning of five trace metals (Co, Ni, Cu, Zn and Pb) in four LOIS rivers has been modelled using a predictive chemical speciation code (WHAM-SCAMP). Observed (log K(D,obs)) and predicted (log K(D,pred)) log K(D) values were similar for Zn. For Co and Ni, the log K(D,pred) values were typically greater than the log K(D,obs) values, while for Cu and Pb the reverse was seen. Removing modelled competition by Ca and Mg for binding sites on particulate organic matter increased the log K(D,pred) values for all the metals except Pb, and gave better agreement between the observations and the predictions for Co, Ni and Zn. Modelling solution iron as iron oxide particles decreased the log K(D,pred) values for Pb. The Cu predictions were very sensitive to the metal binding strength of dissolved and particulate organic matter. Overall, the model shows promise for the prediction of solid-solution metal partitioning in aquatic systems. The modelling exercise has identified the following uncertainties: the extent of Ca and Mg competition for binding sites, the chemical nature of the measured particulate metal, the efficacy of the solid-solution separation method and the strength of copper binding to organic matter.

Modelling the chemical speciation of trace metals in the surface waters of the Humber system

Calculations have been performed to estimate the chemical speciation at equilibrium of six divalent trace metals (Co, Ni, Cu, Zn, Cd, Pb) in riverine, estuarine and marine surface waters of the Humber system. The Windermere Humic Aqueous Model (WHAM) was used to compute distributions of dissolved metals. In the rivers, the free aquo ion (M2+) is a major part of dissolved Co, Ni, Zn and Cd, but accounts for less than 1% of Cu and Pb. The main complexes are formed with carbonate ligands and dissolved natural organic matter, represented by fulvic acid. In the low-salinity region of the estuary and in seawater, complexation with fulvic acid is less significant, although most of the Cu is still in this form, while the speciation of Cd is dominated by chloride complexes. Adsorption of metals by suspended particulate matter was calculated with a simple model involving the concentrations of the free aquo ions (M2+) and H+, together with a constant for each metal estimated from laboratory adsorption data. Calculated adsorbed concentrations were used to predict the partition coefficient (KD) for each metal under different circumstances. The values can vary by an order of magnitude or more, depending upon solution conditions. Typical values for rivers, low-salinity water and seawater are within one order of magnitude of observations. However, there is a general tendency to underestimate KD, possible reasons being (1) neglect of electrostatic enhancement of adsorption at low ionic strengths; and (2) analytical overestimation of particulate metal in equilibrium with the solution phase. There is a strong case for the development of a more sophisticated adsorption model.