Effects of Ionic Strength on Arsenate Adsorption at Aluminum Hydroxide-Water Interfaces

Adsorption processes at mineral-water interfaces control the fate and transport of arsenic in soils and aquatic systems. Mechanistic and thermodynamic models to describe this phenomenon only consider inner-sphere complexes but recent observation of the simultaneous adsorption of innerand outer-sphere arsenate on single crystal surfaces complicates this picture. In this study, we investigate the ionic strength-dependence of the macroscopic adsorption behavior and molecular-scale surface speciation of arsenate bound to gibbsite and bayerite. Arsenate adsorption decreases with increasing ionic strength on both minerals, with a larger effect at pH 4 than pH 7. The observed pH-dependence corresponds with a substantial decrease in surface charge at pH 7, as -potential measurements. EXAFS spectroscopy finds that the number of second shell Al neighbors around arsenate is lower than that required for arsenate to occur solely as an innersphere surface complex. Together, these observations demonstrate that arsenate displays macroscopic and molecular-scale behavior consistent with the co-occurrence of innerand outersphere surface complexes. This demonstrated that outer-sphere species can be responsible for strong adsorption of ions and suggests that environments experiencing an increase in salt content may induce arsenic release to water, especially under weakly acidic conditions.


Introduction
Adsorption is a critical control on the occurrence and fate of arsenic in soils and aquatic systems [1][2][3][4][5] and is utilized to remove this toxic element during water treatment [6,7].Arsenate (H x AsO 4 x−3 ), the common form of arsenic under oxidizing conditions, strongly adsorbs to many mineral surfaces, especially to aluminum and iron (oxyhydr)oxides [1,4,8,9].The macroscopic uptake of this oxoanion onto oxide minerals has not been observed to display a dependence on ionic strength in prior studies [8][9][10], indicating a lack of outer-sphere complexes [11].Previous spectroscopic studies have widely observed that arsenate forms inner-sphere surface complexes on mineral surfaces [4,8,9,[12][13][14][15][16][17].Surface complexation models for arsenate adsorption predict the partitioning of this oxoanion between water and mineral surfaces using only such complexes [9,[18][19][20].A description of arsenate adsorption behavior based solely on an inner-sphere (ligand exchange) binding mechanism was thus established over the last few decades that is widely accepted today.This picture became more complicated when the simultaneous occurrence of both inner-and outer-sphere arsenate surface complexes was observed on aluminum and iron oxide single crystal surfaces [21].While the transferability of processes on such highly-idealized surfaces to particles representative of natural and engineered sorbents is uncertain, a re-evaluation of past observations is warranted.Detection of inner-sphere complexes using extended X-ray absorption fine structure Soil.Syst.2018, 2, 1 2 of 13 (EXAFS) spectroscopy does not preclude co-existing outer-sphere species, as the latter lack diagnostic features in EXAFS spectra [21].Hydrogen bonding of outer-sphere molecules to surface functional groups may yield infrared (IR) spectra of arsenate similar to what would be expected for inner-sphere species [16].Furthermore, macroscopic uptake of arsenate will not display a dependence on ionic strength where outer-sphere complexes are neutrally-charged or when the surface potential is near zero.Prior work with mineral powders thus often cannot unambiguously rule out the co-occurrence of outer-and inner-sphere arsenate surface complexes.
Our recent work investigating binding site effects on arsenate adsorption to Al(OH) 3 polymorphs, gibbsite and bayerite, observed a lower coordination number of second-shell Al neighbors for gibbsite, indicating an additional surface species was present [22].Other inner-sphere binding configurations were not present, suggesting co-occurring outer-sphere arsenate species.This molecular-scale picture conflicts with previous macroscopic studies that have not observed ionic strength-dependent arsenate adsorption behavior on aluminum and iron oxides [8][9][10].Antelo et al. [10] did not detect observable changes in arsenate adsorption isotherms for goethite between 0.01 and 0.1 M KNO 3 fluids.Goldberg and Johnston [9] observed a possible increase in adsorption with increasing ionic strength in adsorption edge measurements on amorphous iron oxide but saw no ionic strength dependence on amorphous aluminum oxide.Arai et al. [8] observed no ionic strength effect on arsenate adsorption to γ-Al 2 O 3 , but their data contains substantial scatter and high percentages of arsenate adsorption at low pH, making detection of ionic strength effects difficult.In all cases, a substantial fraction clearly occurs as inner-sphere species based on many past spectroscopic studies [4,8,9,[12][13][14][15][16][17] and the ionic strength-dependence of macroscopic adsorption would be smaller in magnitude than for an ion that bind solely as an outer-sphere species.Exploration of the ionic strength dependence of arsenate adsorption to aluminum hydroxides is thus warranted given the contrasting past results.
In this study, we seek to determine whether macroscopic arsenate adsorption onto aluminum hydroxides displays a dependence on ionic strength that is consistent with our recent molecular-scale observations suggesting that inner-and outer-sphere species co-occur.Both gibbsite and bayerite are studied because they differ in morphology and thus have distinct distributions of surface functional group types [23] that may lead to a distinct ionic strength-dependence of arsenate adsorption.Adsorption isotherms were measured on both minerals at pH 4 and 7 with ionic strength of 0.001, 0.01, 0.1, and 0.5 M NaNO 3 .ζ-potential measurements assessed the sign and relative magnitude of gibbsite and bayerite surface charge in the presence and absence of adsorbed arsenate at the different ionic strength value used for the isotherm measurements.EXAFS spectroscopy characterized arsenate adsorption mechanisms on each mineral surface at low and high arsenate surface coverages at low and high ionic strength at both pH 4 and 7.These macroscopic and molecular-scale studies were integrated to assess the influences of ionic strength on arsenate adsorption on aluminum hydroxide surfaces.

Mineral and Reagent Preparation
Gibbsite and bayerite were prepared following previously published syntheses [22], modified from the procedures described in Shen et al. [24] and Lefevre et al. [25].After synthesis, gibbsite and bayerite were washed repeatedly in deionized water (>18.2MΩ•cm) to remove excess electrolytes.Both minerals were resuspended in deionized water and stored in polypropylene bottles.The concentrations of the gibbsite and bayerite mineral suspensions were determined through gravimetric analysis.Both gibbsite and bayerite particles have been previously characterized by Brunauer-Emmett-Teller (BET) analysis, powder X-ray diffraction (XRD), and scanning electron microscope (SEM) [22].

Macroscopic Arsenate Adsorption Isotherms
Arsenate adsorption onto gibbsite and bayerite at pH 4 and 7 at solution ionic strength of 0.001, 0.01, 0.1, and 0.5 M NaNO 3 were investigated by batch adsorption isotherm experiments.Samples contained 4 g/L of gibbsite or bayerite reacted in a solution with an initial arsenate concentration from 10 µM to 500 µM, the desired concentration of NaNO 3 (0.001, 0.01, 0.1, and 0.5 M), and, at pH 7, a pH buffer of 0.001 M MOPS, added as the hemisodium salt (no buffer was used at pH 4).It should be noted that the addition of MOPS hemisodium salt increases the ionic strength values to 0.0015, 0.0105, 0.1005, and 0.5005 M for the pH 7 system.For each series of experiment, four mineral-free blanks with 40, 100, 200, and 400 µM arsenate were prepared to assess systematic errors in the initial arsenate concentrations.The samples were rotated on end-over-end rotators for 24 h at room temperature.During the reaction time, sample pH was checked and adjusted to the target values.At the end of the experiment, the samples were filtered (0.22 µm, MCE; Santa Cruz Biotechnology, Dallas, TX, USA), discarding the first 1 mL of the filtrate, and acidified to 2% HNO 3 (trace metal grade; EMD Millipore, Billerica, MA, USA).The acidified filtrates were diluted appropriately with 2% HNO 3 for inductively coupled-plasma mass spectrometry (ICP-MS, Thermo Scientific iCAP Q, Bremen, Germany) analysis for dissolved arsenate concentrations.

ζ-Potential Analysis
ζ-potential analysis was performed to study the electrokinetic response of gibbsite and bayerite particles to different ionic strength with and without arsenate over pH range of 3-8.Suspensions were prepared in the same way as the batch adsorption samples with 40 and 400 µM arsenate.Samples with no arsenate were prepared as the blank references.Suspensions (1 mL) were withdrawn with syringes and injected into a ζ-Potential analyzer (Zetasizer Nano ZS, Malvern, Southborough, MA, USA).Each data point is the average of three replicate measurements.

EXAFS Spectroscopic Measurements
The coordination environment of arsenate adsorbed to gibbsite and bayerite at different ionic strength was investigated with EXAFS spectroscopy.Separate batch adsorption experiments were conducted to prepare fresh samples before beamline measurements.Scaled-up reactors (50 mL total volume) were used to prepare samples with sufficient solid mass for EXAFS measurements.Two adsorbed arsenate samples for each mineral at the lowest and highest background electrolyte concentrations (I = 0.001 and 0.5 M NaNO 3 ) at each pH were prepared with initial low (40 µM) and high (400 µM) loading of arsenate (Table 1).All samples were centrifuged to separate the supernatant and solid (wet paste).The solids were packed as wet paste into polycarbonate sample holders, sealed with 25 µm Kapton tape, and heat sealed in a polyethylene bag with a damp Kimwipe to prevent drying.
Samples were prepared within 3 days of transport to the Stanford Synchrotron Radiation Lightsource (SSRL) at the SLAC National Accelerator Laboratory for measurements of As K-edge EXAFS spectra at SSRL beamline 4-1.Beamline 4-1 uses a Si (220) double crystal monochromator, which was detuned by 35% to attenuate beam harmonics.Fluorescence yield spectra were measured using a 30-element energy-dispersive Ge detector.Energy was calibrated using a Au metal foil, with the first inflection point in the Au L III -edge set to 11919 eV.
The Athena interface [26] to IFEFFIT [27], together with the SamView component of SixPack [28], were used to average, deadtime correct, and normalize the spectra.The mansfieldite (AlAsO 4 •2H 2 O) structure [29] was used to calculate phase-shift and backscattering amplitude functions for structural analysis of the EXAFS spectra using FEFF 7.02 [30].The k 3 -weighted EXAFS spectra were fit over a k range of 3.5-12.5Å and an R range of 0.8-4 Å, with the amplitude reduction factor (S 0 2 ) fixed to 1.0.All spectra for each sample were fitted using a single starting structural model consisting of As-O The parameters for these paths were determined from geometric considerations [31].In the model fits, As-O coordination numbers were fixed at 4 and the coordination numbers of As-Al shells were allowed to vary.For the As-Al single scattering path, the Debye-Waller factor, σ 2 , was fixed to 0.006 to reduce covariation during fitting.

Properties of Synthesized Materials
The chemical and physical properties of the gibbsite and bayerite used in this research have been previously described [22].Briefly, the gibbsite and bayerite used in this study are aluminum hydroxide polymorphs with different morphology.Both mineral particles are enclosed by basal and edge surfaces with a substantially different ratio of these two types of surfaces.Synthetic gibbsite with a platelet morphology and bayerite with a micro-rod morphology were previously estimated to have ~91% of the basal (001) planes, terminated by doubly coordinated (>Al 2 OH) functional groups, and ~94% of the edge planes, terminated by singly coordinated (>AlOH 2 ) functional groups, respectively [22].The gibbsite and bayerite particles have similar surface areas (21.8 m 2 /g and 23.4 m 2 /g, respectively) determined by BET analysis [22].Both minerals used in this study were previously characterized by XRD and SEM, confirming they are consistent with the expected products of the published mineral synthesis methods, i.e., single phase, relatively monodispersed particles of the target mineral.

Macroscopic Arsenate Adsorption
The macroscopic arsenate adsorption isotherms on both minerals at each ionic strength show similar trends, with arsenate displaying greater adsorption at low pH and on bayerite versus gibbsite (Figure 1), consistent with our prior study [22].The macroscopic arsenate adsorption data are then fitted using Langmuir isotherm with co-optimized maximum adsorption capacity parameters (Γ max,i ) and affinity constants (K i ): where Γ ads is the amount of adsorbate adsorbed, C (ads) is the final adsorbate concentration in solution, n = 1 for single Langmuir fits involving only one adsorption site, and n = 2 for dual Langmuir fits involving two adsorption sites.The adsorption data are better fitted using the dual Langmuir isotherm model for both gibbsite and bayerite (Figures 1 and S1).Comparing the fitting parameters in the dual Langmuir model, the total Γ max values (Γ max,1 + Γ max,2 ) of bayerite are larger than those of gibbsite at same pH and same ionic strength (Table 2), which could be attributed to the abundance of edge surface sites on bayerite particles displaying stronger adsorption capacity.However, under the same experimental condition, the two components Γ max,1 and Γ max,2 do not vary in proportion to the relative areas of basal and edge surface sites on both minerals, suggesting that dual Langmuir behavior is not related to particle morphology but the presence of two different surface complexes (e.g., inner-and outer-sphere species) on both minerals.These macroscopic adsorption behaviors in the present study have also been observed in our previous study [22].It should be noted that the very large K 1 value in the dual Langmuir model fit for gibbsite at pH 7 in 0.5 M NaNO 3 is likely not determined uniquely because of the lack of adsorption data at very low arsenate concentration.This is reflected in the small value of Γ max,1 and the large uncertainty in K 1 .
Soils 2018, 2, 1 5 of 13 where Γads is the amount of adsorbate adsorbed, C(ads) is the final adsorbate concentration in solution, n = 1 for single Langmuir fits involving only one adsorption site, and n = 2 for dual Langmuir fits involving two adsorption sites.The adsorption data are better fitted using the dual Langmuir isotherm model for both gibbsite and bayerite (Figure 1 and Figure S1).Comparing the fitting parameters in the dual Langmuir model, the total Γmax values (Γmax,1 + Γmax,2) of bayerite are larger than those of gibbsite at same pH and same ionic strength (Table 2), which could be attributed to the abundance of edge surface sites on bayerite particles displaying stronger adsorption capacity.However, under the same experimental condition, the two components Γmax,1 and Γmax,2 do not vary in proportion to the relative areas of basal and edge surface sites on both minerals, suggesting that dual Langmuir behavior is not related to particle morphology but the presence of two different surface complexes (e.g., inner-and outer-sphere species) on both minerals.These macroscopic adsorption behaviors in the present study have also been observed in our previous study [22].It should be noted that the very large K1 value in the dual Langmuir model fit for gibbsite at pH 7 in 0.5 M NaNO3 is likely not determined uniquely because of the lack of adsorption data at very low arsenate concentration.This is reflected in the small value of Γmax,1 and the large uncertainty in K1.For each mineral, a change in ionic strength has substantial effects on arsenate uptake, with adsorption decreasing as ionic strength increases (Figure 1).These ionic strength-dependence effects on arsenate adsorption are larger in magnitude at pH 4 and on bayerite (Figure 1).At pH 4 the effect of ionic strength varies with final arsenate concentration (CAs): low CAs (<0.3 μmol/L) regions display little variation in uptake with ionic strength, intermediate conditions (0.3 μmol/L< CAs <60 μmol/L) showed the greatest ionic strength effect, and high CAs (>60 μmol/L) regions showed a weaker ionic strength effect.In contrast, at pH 7 the effect of ionic strength on adsorption appears uniform across all arsenate concentrations.These results display clear macroscopic signatures associated with the For each mineral, a change in ionic strength has substantial effects on arsenate uptake, with adsorption decreasing as ionic strength increases (Figure 1).These ionic strength-dependence effects on arsenate adsorption are larger in magnitude at pH 4 and on bayerite (Figure 1).At pH 4 the effect of ionic strength varies with final arsenate concentration (C As ): low C As (<0.3 µmol/L) regions display little variation in uptake with ionic strength, intermediate conditions (0.3 µmol/L< C As < 60 µmol/L) showed the greatest ionic strength effect, and high C As (>60 µmol/L) regions showed a weaker ionic strength effect.In contrast, at pH 7 the effect of ionic strength on adsorption appears uniform across all arsenate concentrations.These results display clear macroscopic signatures associated with the presence of outer-sphere species, i.e., a decrease in adsorption with increasing ionic strength.The smaller effect seen at pH 7 may imply such species are of lower abundance under neutral conditions, but the magnitude of ionic strength effects on outer-sphere adsorption varies with the Soil.Syst.2018, 2, 1 6 of 13 surface potential.If the surface charge is largely neutralized at pH 7 by arsenate adsorption, then little to no effect of ionic strength, is expected even if outer-sphere species are present.a Goodness-of-fit parameter.The R-factor is the sum of the differences between the data and the fit at each data point, divided by the sum of the data at each corresponding point.Smaller R-factor values reflect better fits.

Ionic Strength Effect on Surface Charge Properties
ζ-potential measurements were made to evaluate how different background electrolyte concentrations affect surface charging, and thus trends in the surface potential of gibbsite and bayerite particles in the absence and presence of arsenate.Over the pH range studied, gibbsite and bayerite suspensions have positive ζ-potential values in the absence of arsenate (Figure 2), consistent with previous studies of aluminum hydroxides with a typical isoelectric point (IEP) around pH 9-10 [23,32-35].At individual pH conditions, ζ-potential for both minerals decreases with increasing ionic strength (Figure 2), consistent with the known effect of counterion screening of surface charge, reducing the potential at the shear plane.After arsenate loading, the ζ-potential values of both minerals shift to more negative values, as expected for anion adsorption.For each mineral system, the ζ-potential value is lower near pH 7 than near pH 4, indicating that the surface potential is also lower at pH 7, and ionic strength-dependence effects on outer-sphere species are thus expected to be substantially greater at pH 4 (Figure 2).The ζ-potential appears to always be larger for bayerite than gibbsite in the presence of arsenate below pH 6 (Figure 2), indicating that the bayerite systems should display stronger ionic strength effects at acidic conditions.These two features are consistent with the observations in adsorption isotherm experiments in this study.In addition, the IEP differs between the two minerals, occurring near pH 7 for gibbsite at both arsenate coverages, but remaining above pH 8 for bayerite except at the highest arsenate coverage (Figure 2).

Arsenate Adsorption Mechanisms
Arsenic surface speciation for a series of adsorption samples (Table 1) was characterized using EXAFS spectroscopy to determine the effect of ionic strength on arsenate binding mechanisms.Conditions are chosen to explore low and high arsenate surface loading at the highest and lowest electrolyte concentrations, 0.5 and 0.001 M NaNO3.Structural model fits of the EXAFS spectra yield interatomic distances of ~3.2 Å for As-Al shell for all samples regardless of pH, ionic strength, arsenate coverage, and mineral sorbents (Figure 3, Table 3).This configuration corresponds to a bridging bidentate inner-sphere surface complex, as has been identified in most previous studies [1,8,17,36,37].All spectra for each sample were fitted using a single starting structural model consisting of only one Al neighbor corresponding to this bidentate binding configuration.As noted in our prior study [22], we also evaluated whether monodentate species are present by adding a second As-Al path at R of 3.5 to 3.6 Å in our initial model.However, this did not improve the quality of fit and resulted in a CN for this longer Al shell within an error of zero for all sample studied.This indicates that monodentate inner-sphere species are negligible on these minerals for all conditions examined.
The EXAFS-derived coordination number of Al neighbors associated with the inner-sphere complex on both minerals is less than 2 (Table 3), the expected value if all the adsorbed arsenate was in a bridging bidentate geometry.This could indicate the presence of outer-sphere arsenate complexes as other inner-sphere binding configurations were not observed under our experimental conditions [22].Also, this coordination number of Al neighbors is lower on average on gibbsite than

Arsenate Adsorption Mechanisms
Arsenic surface speciation for a series of adsorption samples (Table 1) was characterized using EXAFS spectroscopy to determine the effect of ionic strength on arsenate binding mechanisms.Conditions are chosen to explore low and high arsenate surface loading at the highest and lowest electrolyte concentrations, 0.5 and 0.001 M NaNO 3 .Structural model fits of the EXAFS spectra yield interatomic distances of ~3.2 Å for As-Al shell for all samples regardless of pH, ionic strength, arsenate coverage, and mineral sorbents (Figure 3, Table 3).This configuration corresponds to a bridging bidentate inner-sphere surface complex, as has been identified in most previous studies [1,8,17,36,37].All spectra for each sample were fitted using a single starting structural model consisting of only one Al neighbor corresponding to this bidentate binding configuration.As noted in our prior study [22], we also evaluated whether monodentate species are present by adding a second As-Al path at R of 3.5 to 3.6 Å in our initial model.However, this did not improve the quality of fit and resulted in a CN for this longer Al shell within an error of zero for all sample studied.This indicates that monodentate inner-sphere species are negligible on these minerals for all conditions examined.
The EXAFS-derived coordination number of Al neighbors associated with the inner-sphere complex on both minerals is less than 2 (Table 3), the expected value if all the adsorbed arsenate was in a bridging bidentate geometry.This could indicate the presence of outer-sphere arsenate complexes as Soil.Syst.2018, 2, 1 8 of 13 other inner-sphere binding configurations were not observed under our experimental conditions [22].Also, this coordination number of Al neighbors is lower on average on gibbsite than that on bayerite.It should be noted that Al coordination numbers obtained from individual samples are not statistically distinguishable because of the large fitting uncertainties.However, statistical tests (t-tests), as applied in the prior study [22], showed that the difference in the average Al coordination numbers for two mineral sets is significant.The lower average Al coordination number on gibbsite indicates a greater proportion of outer-sphere complexes.For individual mineral system, similar statistical tests were performed to evaluate if the average Al coordination number of high ionic strength set differs from that of low ionic strength set.t-tests show that average Al coordination numbers are statistically invariant with different ionic strength for both minerals, suggesting that ionic strength has little effect on arsenate surface complexation geometry as well as the distribution of different types of arsenate surface species on gibbsite or bayerite.
Soils 2018, 2, 1 8 of 13 that on bayerite.It should be noted that Al coordination numbers obtained from individual samples are not statistically distinguishable because of the large fitting uncertainties.However, statistical tests (t-tests), as applied in the prior study [22], showed that the difference in the average Al coordination numbers for two mineral sets is significant.The lower average Al coordination number on gibbsite indicates a greater proportion of outer-sphere complexes.For individual mineral system, similar statistical tests were performed to evaluate if the average Al coordination number of high ionic strength set differs from that of low ionic strength set.t-tests show that average Al coordination numbers are statistically invariant with different ionic strength for both minerals, suggesting that ionic strength has little effect on arsenate surface complexation geometry as well as the distribution of different types of arsenate surface species on gibbsite or bayerite.1.  1. between the data and theory.e Goodness-of-fit parameters [38].f The estimated standard deviations are listed in parentheses, representing the errors in the last digit.

Implications to Arsenate Adsorption Mechanisms
Our macroscopic and spectroscopic measurements on aluminum hydroxide powders, together with the previous surface scattering study of arsenate adsorbed on aluminum and iron oxide single crystal surfaces [21], now provide a generally self-consistent description of arsenate adsorption, indicating that inner-and outer-sphere arsenate species co-adsorb over a wide range of surface coverages and pH.Arsenate adsorption does show the classical dependence on ionic strength, with a substantial effect when the surface potential is large and a negligible effect when the surface potential is near zero.All of these observations are the expected macroscopic behaviors when outer-sphere arsenate surface complexes are present.
However, these results also reveal unexpected observations.The decrease in arsenate adsorption with increasing ionic strength in the isotherms is expected to reduce the fraction of adsorbed arsenate that occurs as an outer-sphere complex, but no clear trend in the coordination number of second shell Al neighbors with increasing ionic strength was observed in EXAFS spectroscopy.This may be because the uncertainty in EXAFS-derived coordination numbers (typically ± 0.5) is too large to observe the expected trends with the number of analyses available.Alternatively, other aspects of the adsorption process may be affected by ionic strength in ways that do not substantially change the balance of inner-and outer-sphere species.For example, a prior study using the extended triple-layer model suggested a possible reaction stoichiometry involving both inner-and outer-sphere arsenate species in a ratio independent of pH and surface coverage and only sensitive to the activity of water [39].Alternatively, if applying a charge distribution model then inner-sphere complexes have two O ligands extending to the β-plane, producing an ionic strength-dependence on macroscopic uptake [19].While this would be weaker in magnitude than for outer-sphere species, the reduction in surface coverage combined with the EXAFS fitting uncertainties may be adequate to explain the apparent lack of change in the distribution of arsenate surface species.More broadly, the ratio of inner-to outer-sphere arsenate species appears to be stable across a wide range of surface coverages [21,22], pH values [22], and salt contents.Reaction stoichiometries that reproduce this phenomenon need to be investigated in future refinements of surface complexation models.
A separated feature of note in the results is an apparent correlation of the Dual Langmuir model fitting parameters (i.e., K and Γ max values) and ionic strength (Table 2).Γ max,1 values decrease with increasing ionic strength for both minerals at both pH conditions, with a weaker variation at pH 7 (Figure S2A).K 1 values show distinct behaviors at the two pH values studied, increasing with increasing ionic strength at pH 4 but decrease at pH 7 (Figure S2C).In contrast, Γ max,2 and K 2 values display little variation as ionic strength changes for all conditions (Figure S2B,D).If the distinct strong and weak correlations of K 1 and K 2 with ionic strength originated from K 1 attributed to inner-sphere adsorption and K 2 attributed to outer-sphere adsorption, a change in the population of inner-and outer-sphere species with increasing ionic strength and arsenate surface coverage should occur.However, as noted above this is not observed in our EXAFS results.The origin of the effects of ionic strength on the dual Langmuir isotherm components is thus uncertain.Many factors can induce variations in K values, including changes in activity coefficients for the dissolved species, surface site competition with electrolyte ions, Na-AsO 4 complexation in solution, and changes in surface potential as ionic strength and arsenate surface coverage increases.Thus, the observed trends in K with ionic strength are most likely attributed to multiple origins instead of a single origin.The smaller difference in Γ max,1 seen at pH 7 compared to pH 4 may be dominantly controlled by the difference in surface potential (large surface potential at pH 4 and small to zero one at pH 7) but also may be affected by changes in solution speciation and Na + and NO 3 -interactions with mineral surfaces.
However, the spectroscopic results clearly rule out the strongly (Langmuir isotherm 1) and weakly (Langmuir isotherm 2) bound arsenate species corresponding to inner-and outer-sphere complexes.
The coexistence of such species is surprisingly unrelated to the apparent adsorption isotherm behavior.Surface complexation models are needed to identify the aspects of the observed variations in uptake behaviors with increasing ionic strength that are induced by activity corrections, solution complexation, surface site competition, and surface potential changes.

Implications to Arsenate in the Environment
The strong adsorption of arsenate to many aluminum and iron (hydro)oxide minerals widely observed in past studies [8,9,12,13,17] has traditionally been attributed solely inner-sphere surface complexes.However, the co-adsorption of inner-and outer-sphere arsenate species indicates that the latter species are partially responsible for this strong adsorption behavior.While such species appear to have high affinity for the mineral surface, they are expected to be more kinetically labile because of their adsorption-desorption does not involve ligand exchange.This suggests that they may desorb more rapidly from mineral surfaces than inner-sphere species, although experimental confirmation is needed as desorption of outer-sphere species was recently shown to be the rate-limiting step in overall Rb + desorption on muscovite [40].Past macroscopic kinetic studies for arsenate have often observed two-stage (i.e., fast and slow) desorption from powder substrates [41][42][43][44] and the fast component may be attributable to the outer-sphere species suggested in the current work.Arsenate adsorption behavior is thus more complex than previously recognized, with inner-and outer-sphere surface complexes displaying similar adsorption affinities and co-occurring over a wide range of conditions.

Soil. Syst. 2018, 2 , 1 4 of 13 and
As-Al single scattering paths and three multiple scattering paths: a triangular As-O-O path, a collinear As-O-As-O path, and a non-collinear As-O-As-O path.

Figure 1 .
Figure 1.Arsenate adsorption isotherms as a function of ionic strength for bayerite at (A) pH 4 and (B) pH 7, on gibbsite at (C) pH 4 and (D) pH 7. Solids lines represent dual Langmuir isotherm fits to the data, with detailed fitting information provided in Table2.

Figure 1 .
Figure 1.Arsenate adsorption isotherms as a function of ionic strength for bayerite at (A) pH 4 and (B) pH 7, on gibbsite at (C) pH 4 and (D) pH 7. Solids lines represent dual Langmuir isotherm fits to the data, with detailed fitting information provided in Table2.

Figure 2 .
Figure 2. ζ-potential as a function of pH for systems with no arsenate on bayerite (A) and gibbsite (B), 40 μM arsenate on bayerite (C) and gibbsite (D), 400 μM arsenate on bayerite (E) and gibbsite (F).The error bars represent the standard deviation of three replicate measurements.

Figure 2 .
Figure 2. ζ-potential as a function of pH for systems with no arsenate on bayerite (A) and gibbsite (B), 40 µM arsenate on bayerite (C) and gibbsite (D), 400 µM arsenate on bayerite (E) and gibbsite (F).The error bars represent the standard deviation of three replicate measurements.

Figure 3 .
Figure 3. Data (dotted) and structural fits (solids) to the As K-edge spectra (A,D), Fourier transform magnitudes (B,E), and real components of the Fourier transforms (C,F) of the series of adsorption samples.Detailed sample information is provided in Table1.

Figure 3 .
Figure 3. Data (dotted) and structural fits (solids) to the As K-edge spectra (A,D), Fourier transform magnitudes (B,E), and real components of the Fourier transforms (C,F) of the series of adsorption samples.Detailed sample information is provided in Table1. 4 •7H 2 O; Sigma Aldrich, ACS Reagent ≥ 98%) stock solutions were made with deionized water.Two 20 g/L suspensions of gibbsite and bayerite were prepared in deionized water.Experiments

Table 1 .
Conditions for the adsorption samples analyzed by Extended X-ray absorption fine structure (EXAFS) spectroscopy.

Table 2 .
Fitting parameters for dual Langmuir isotherms.

Table 3 .
As K-edge EXAFS fitting parameters.
a Coordination number.b Interatomic distance.c Debye-Waller factor.d Difference in the threshold Fermi level