J magi) WIN WHWHHTIW N 312910087 9140 1 " W manner Michigan 3mm 1 University \ f This is to certify that the dissertation entitled PHOSPHATE EXCHANGE BETWEEN LITTORAL SEDIMENT AND LAKE WATER presented by Robert Paul Glandon has been accepted towards fulfillment of the requirements for _Bh_._D_.__degree in Eisherieis. Wildlife 0%722/170010 *ajor pk ‘ofessor Date October 20, 1982 MS U is an Affirmative Action/Equal Opportunity Institution 042771 MSU RETURNING MATERIALS: Place in book drop to ALW~ LJBRARJES remove this checkout from “ your record. FINES Win be charged if book is returned after the date stamped below. HS 2 6 199? PHOSPHATE EXCHANGE BETWEEN LITTORAL SEDIMENT AND LAKE WATER BY ROBERT PAUL GLANDON A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Fisheries and Wildlife 1982 ABSTRACT PHOSPHATE EXCHANGE BETWEEN LITTORAL SEDIMENT AND LAKE WATER BY Robert Paul Glandon Soluble phosphate is an important form of phosphorus in the environment. Current models, which consider a largely hydrostatic sediment-water interface, hold that aerobic sedi- ments act as a barrier to sediment phosphate flux to over- lying water. However, differences between littoral sediment interstitial and lake water phosphate levels indicate poten- tial for diffusive flux. Further, water movement that re- sults in perfusion of epilimnetic water or resuspension of littoral sediment particulates, can modify the particle- solution phosphate system. Information concerning mechanisms of particle response to changes in ambient phosphate or vol- ume of ambient water can enhance models of phosphate exchange under hydrodynamic conditions. Interstitial phosphate of sediments underlying 0.75- 1.00 m water in a shallow lake was variable in space. Lev- els were generally greater than those of overlying water and gradients could be used to estimate static-water phos- phate eff: solution c water, mes proached I Particulai rected by of associz Sedir sampling 5 of particx ponse of 1 levels imI phosphate ambient p1 phate sor} condition; described sorption 1 concepts « ulate res, depends 0 water, an ting from Robert Paul Glandon 1x through surficial sediment layers. Assuming ilcium activity and pH similar to that of overlying 1 soluble phosphate of surficial sediments ap- cedicted levels in equilibrium with hydroxyapatite. a response to ambient perturbation could be di- phosphate mineral equilibria and limited by kinetics ted phosphate exchange. ent material was taken from the site of interstitial ad a portion was treated to modify characteristics late surfaces. Experimental observations of res- articulates to alterations in lake water phosphate licated controlling mechanisms other than calcium solubility. Langmuir-type analysis suggested that osphate level is governed by the reservoir of phos- ed on sediment particles. Relative to experimental , particle-solution phosphate relationships were using quantified parameters related to phosphate aximum and binding energy. A model, incorporating rawn from experimental data, indicates that partic- onse to perturbation of a sediment-phosphate system phosphate activity of interstitial and infusing particle weight to solution volume ratio resul- iisturbance. DEDICATION To my son, Kevin (5 T3 I ACKNOWLEDGEMENTS I would like to express my sincere appreciation to Dr. Clarence D. McNabb whose open-minded approach has been inspiring. Thanks to Drs. Boyd G. Ellis and Bernard D. Knezek for their direct and indirect contributions to this effort. Many of my fellow students deserve thanks, particularly Ted R. Batterson, Frederick C. Payne, and Mehdi Siami. Each has contributed to my graduate experience. My wife, Nancy, deserves thanks for everything. My parents deserve special recognition for their con- tinous encouragement. This study was supported in part by funds from the U.S. Environmental Protection Agency, and the Michigan Agricultural Experimentation Station at Michigan State University. (-3 V5. «4. TABLE OF CONTENTS Page LIST or TABLES ... ....... . ..... ........................ v LIST OF FIGURES ....................................... vii INTRODUCTION .......................................... 1 MATERIALS AND METHODS ................................. 7 RESULTS ............................................... 16' DISCUSSION ............ ..... ........................... 43 LITERATURE CITED ...................................... 69 APPENDIX 1 ............................................ 74 APPENDIX 2 ............................................ 79 APPENDIX 3 ............................ ..... ........... 81 APPENDIX 4 . ..... ............... ....................... 84 iv Table 1. LIST OF TABLES Page Calculated and measured net phosphate exchange (uM P04 g'l) between whole sediments and in- fusing lake water of different initial phosphate activity (uM P04 1'1). Values assume equilibrium conditions with ionic strength 3 x 10’ M and pH 7.60 . Negative Sign indicates net sediment sorption of phosphate ........................... 63 Calculated and measured net phosphate exchange (UM P04 1’1) between fine sediments and infusing lake water of different initial phosphate act- ivity (uM P04 171). Calculated exchange is based on experimentally estimated ambient cross- over activity of 18.63 UM P04 1‘1 (A), and on a hypothetical crossover activity of 22 EM P04 1' (B). Values assume equilibrium conditions with ionic strength 3 x 10'3 M and pH 7.75 ...... 64 Effect of measured additions of arsenate (added as Na2HASO4-7H20) on determinations of phosphate in standard (KH2P04) preparations ............... 75 Effects of addition of arsenate reductant on determinations of phosphate in distilled water solutions containing measured additions of phos- phate and arsenate .............................. 78 Concentration of soluble reactive and total phos- phate, iron, and calcium in shallow sediment collections and lake water used in processing fine sediment. Corresponding values represent estimates of the amount of each component de- posited during drying (at 800C) of 96 g of sedi- ment material ................................... 80 Filterable iron (uM l—l) in lake water of differ- ent initial phosphate levels following 0.5, 2.0, and 48 hours exposure to whole sediments. Lake water, prior to sediment addition, contained 0 734 uM Fe 1'1 ............................ . ..... 82 Filterable iron (UM l-l) in lake water of differ— ent initial phosphate levels following 0.5, 2.0, Table Page and 48 hours exposure to fine sediments. Lake water prior to sediment addition contained 0.84211MFe1-1 O.......OOOOOOOOOOOOOOO0.0.000... 83 Solubility products of calcium phosphate com- pounds and dissociation constants of phosphoric acid species and water at 25°C used in the development of the solubility diagram of Figure 10 ......OOOOOOOOOIOOOO......OOOOOOOOOOOOOO....O. 85 vi Figure 1. LIST OF FIGURES Page Content of total interstitial phosphorus (mean i 1 standard error) with depth in shallow sedi- ments of the south basin of Lake Lansing in the interval from July - October, 1978 ................ 18 Net reactive phosphate exchange between sedi- ment and lake water following 2 hours (A-A) and 48 hours (x-X) exposure to different initial reactive phosphate levels. Initial ambient levels (uM P0 1‘1) are indicated in parentheses. Horizontal aXIs shows final ambient reactive phosphate level. Panel A results from.whole sediments. Panel B results from fine sediments .. 22 Net Reactive phosphate sorbed by whole sediments (solid curve) and fine sediments (dashed curve) from lake water to different initial reactive phosphate levels, following 2 hours (A-A) and 48 hours (X-X) exposure. Horizontal axis indicates resultant ambient reactive phosphate level. Associated initial levels of ambient reactive phosphate (uM P04 1’1) are indicated in paren- theses ........................................... 27 Relationship between release of soluble nonre- active phosphate and release of reactive phosphate by sediments following 0.5 hours (0-0), 2 hours (A-A), and 48 hours (x-x) exposure to lake water with low initial reactive phosphate levels. Solid line, results from whole sediments; dashed line from fine sediments .............................. 30 Relationship between release of soluble nonreact- ive phosphate and sorption of reactive phosphate by sediments following 0.5 hours (0-0), 2 hours (A-A), 48 hours (X-X) exposure to lake water with high initial reactive phosphate levels. Solid line, result from whole sediments; dash line from fine sediments ........... ..... .............. 32 Ambient calcium (mean and 1 standard deviation) and pH (mean and range) with duration of whole (A-A) and fine (0-0) sediment exposure to lake water phosphate solution. The dashed lines vii Figure Page 10. 11. indicate hypothesized course of change in ambient calcium within initial 30 minutes sediment-lake water exposure. Time 0 indicates point of sediment introduction to lake water......... ...... ........... 35 Net reactive phosphate released from fine sedi- ments under conditions of ambient phosphate accumulation.(A—A), and cumulative phosphate release from fine sediments to ambient lake water replaced at 15 minute intervals (0-0)._l Lake water initially contained 0.10 uM P04 1 . Values are expressed as a mean of six replicates t 1 standard error ................................ 40 Relation between pH and relative proportions of soluble reactive phosphate species. Calcula- tions based on dissociation constants of phos- phoric acid, mono-, and dihydrogen phosphate ions in dilute solution ........................... 46 Influence of ionic strength (u) on activity coefficient (7) of mono— and dihydrogen phos- phate for solutions of ionic strength less than 0.01 molar. V was calculated from: log y = -0.5085 22 /fi_ ................................. 49 Solubility diagram for certain calcium phosphate compounds in sediments at 25°C. Levels of solu- ble reactive phosphate were calculated by consid- ering appropriate solubility equilibria and _3 assuming calcium levels maintained at l x 10 molar (solid lines) and 3 x 10‘3 molar (dashed lines. (-) refer to ambient reactive phosphate levels following 48 hours exposure of whole sediments to high initial levels of ambient phosphate. (X) same, but for fine sediment ex- posure. Phosphate activities were converted to concentrations for comparison with observed levels of total interstitial phosphate of south basin sediments ..... ......... ........................... 52 Calculated relationship between phosphate sorbed by whole and fine sediments and ambient phos- phate activity. Each isotherm was developed from estimations of Langmuir constants B and K. Ionic strength 3 x 10‘3 M and equilibrium with respect to particulate-solution phosphate ex- change were assumed. The pH was taken as 7.60 and 7.75 for whole and fine sediments, respec- tively .................................. . ......... 60 Figure 12. Net exchange of particle associated phosphate as a function of phosphate activity of infusing medium and of interstitial water in an ideal- ized sorption-desorption exchange system. Letters indicate phosphate activity level (As ,. e I -00 0» l° 0" 0.. «IN: 3 103 a) E «'0 SE. .2 D is; L. .... g .3 o O J ng 0X 0 ‘- l_ l l 1 L '0. ‘2 '°. N -: O O O O O (mbgem Mp |.bl'od w ‘0‘) pesoeiea etoquoqd Panacea-mu 29 Figure 4. Relationship between release of soluble nonreactive phosphate and release of reactive phosphate by sediments following 0.5 hours (O-O), 2 hours (A-A), and 48 hours (X-X) exposure to lake water with low initial reactive phosphate levels. Solid line, results from whole sediments; dashed line from fine sediments. 30 4 q 'x .. Ix 022 I. .9 3 4 lo0 '0’ ,i r e I ~00 E» |° °" 0. .1»: 3. .49“: U) .3 .... 0'2. .2 0 “0‘5. 8 a at” g a: o 8 #610: odx 0 ‘- L l J l 1 If! ‘1’ '0. 0! -. o o’ o o o (NOIOM Mp .-b’od w or) posoelea eioquoqd mucosa-non 31 Figure 5. Relationship between release of soluble nonreactive phosphate and sorption of reactive phosphate by sediments following 0.5 hours (O-O), 2 hours (A-A), 48 hours (X-X) exposure to lake water with high initial reactive phosphate levels. Solid line, result from whole sediments; dash line from fine sediments. 32 .222. to .-o 40.. 23 933 2289... 388m E m. N. : O. o m b m n v m N . u u n J d '1 4 d u d 10 u d e. d N x .a m. o a a D mw x LNG M o d m d .3 W x .0... ill 0 a 41x1 4 o \«A m. \ .8 w 0 \\ m. \ .4... \ .3 w \ ..d \ . m x x \ 5.... llll'llll'ik 10.0 p ..A. q 0 m .3 M 33 that mechanisms of sediment uptake of reactive phosphate may involve displacement of phosphate-containing organic molecules. The hypothesis that reactive phosphate added to soil systems displaced sorbed organic phosphate was not supported by the data of Wier and Black (1968). According to Latteral, et a1. (1971), however, reactive phosphate added to and sorbed by sediments resulted in a substantial increase in the proportion of organic to reactive phosphate in equilibrium solution. This was interpreted as being caused by displacement of organic phosphorus from sediment sorption sites by the more strongly sorbed reactive phosphate. Interpretation of changes in solution phosphate levels following sediment introduction partially rests on observations of other system variables that may influence results. Assess- ing potential impact of pH-dependent dissolution/precipitation reactions involving calcium phosphate compounds, or changes in oxidation-reduction potential leading to solubilization of phosphate binding sediment compounds is critical to a discussion of mechanisms of sediment-solution phosphate exchange. Toward this end, solution calcium, pH, and iron were monitored. The general directions of change in ambient pH and calcium were similar in both whole and fine sediment-lake water systems (Figure (H. Further, for each sediment type, trends in ambient calcium and pH were similar at all phosphate levels examined, ie. there were no correlations between the extent of ambient calcium increase or pH decrease and initial ambient phosphate, or extent of inorganic phosphate exchange. Measurements of Figure 6. 34 Ambient calcium (mean and 1 standard deviation) and pH (mean and range) with duration of whole (A-A) and fine (o-o) sediment exposure to lake water phosphate solutions. The dashed lines indicate hy- pothesized course of change in ambient calcium within initial 30 minutes sediment-lake water exposure. Time 0 indicates point of sediment introduction to lake water. 35 (pm 1‘ ..3’ “low ) “WINDS W‘NWV o o o o 9' 0' as as r I I j ,r——--1‘r 1"""’—' onion 5555 L Hd "1°!qu .‘q p. ¢ Duration 01' Sediment-Lake Water Exposure (Hours) 36 soluble calcium and ambient pH were grouped to obtain a mean for each exposure interval. Lake water pH increased from 7.4 to 7.95 and 8.10 upon introduction of whole and fine sediment (time 0). Solution calcium increased over lake water levels by about 20 uM 1-1 within 30 minutes sediment-lake water exposure. The dashed lines indicate hypothesized course of solution calcium change; suggesting that increases in soluble calcium corresponded to observed pH increases (of. Turner and Clark, 1956). These trends may reflect response of calcium carbonate equilibria and as such were largely independent of observed phosphate reactions. Calcium carbonate solubility is a function of pH and carbon dioxide content of water, the relationship can be de- scribed by: 8 pCa = pH + a log PCO2 - 4.93 where pCa is the negative log of the molar concentration of calcium, and FCC is the partial pressure of carbon dioxide. 2 At equilibrium, partial pressure of aqueous CO2 is equal to the partial pressure of atmospheric CO2 (Turner and Clark, 1956). Assuming lake water used for experimentation was in equi- librium with atmospheric CO2 prior to sediment addition, PCO2 = 0.0003 atm, and 5 log PCO = —l.761. At pH 7.4, predicted 2 aqueous calcium level would be 3.8 x 10"2 molar. Measured levels were near 1 x 10-3 M. This indicates a tendency for calcium carbonate, in contact with lake water, to dissolve. Dissolution of CaCO3, and subsequent increase in CO32-, could 37 incur a pH increase. A resultant disequilibrium between aqueous and atmospheric CO2 would promote reinvasion of CO2 from the atmosphere causing a gradual pH decrease in solution. The lab— oratory results suggest that the rate of potential calcium car— bonate dissolution early in the exposure period was greater than the rate of invasion of atmospheric CO The observation 2. that pH-Ca equilibrium was not attained in the duration of the experiments (ie. final calcium levels were lower than predict- ed) may reflect insufficient calcium reserves or reduction in the rate of dissolution of CaCO3 coated with adsorbed organic compounds (Wetzel, 1971, 1972). These considerations suggest that observed calcium increases may not involve dissolution of phosphate containing minerals. Filterable iron content did not change measureably in the sediment-lake water system (Appendix 3), suggesting that redox potentials did not fall below a critial level for chemical re- duction of phosphate-binding iron oxides. While observed levels of soluble iron were low relative to analytical precision, they were high relative to predicted concentrations based on solu- bilities of iron oxides and hydroxides. Assuming particulate iron passing the 0.5 um filter was minimal, the observed iron levels may reflect presence of naturally occurring metal che- lating compounds in solution. The reciprocal nature of phosphate exchange implies, for example, that the extent of particulate release of phosphate, in response to a reduction in ambient levels, is modified by the correspounding increase in solution phosphate. This suggests 38 the importance of water movement characteristics on the extent of phosphate release from sediment particulates. Turbulent displacement of sediment particles to the relatively phosphorus- poor water column could maximize phosphate release from particle surfaces.v Figure 7 shows that phosphate release from fine sediments was greater when ambient phosphate was maintained at relative; ly low levels by lake water replacement. The dashed line of the figure indicates that following 2 hours exposure to 100 ml lake water, 0.5 g sediment released approximately 1.28 uM PO4 g-l; corresponding to an ambient phosphate level of 6.5 uM PO4 1-1. While sediment impact on ambient phosphate level was reduced from that observed with l 9 material per 100 ml (Figure 2) weight-normalized phosphate release was higher. This could indicate inhibiting effects of increased ambient phosphate on rate of net particulate release. Fine sediments released a total of approximately 1.57 uM P04 9.1 in response; to ambient medium renewal. Contrast in sediment response shown in Figure 7 was dom- inated by differences in phosphate release during the 15 - 30 minute exposure interval. Within the initial 15 minute exposure interval, sediments released 1.0 uM P04 9.1 resulting in ambi- 1 ent level near 5.1 uM P04 1. . Materials exposed to renewed medium, containing 0.1 uM P0 1.1, released an additional 0.38 4 0M PO4 g-l, while continued exposure to 5.1 uM P04 1"1 restrict— ed release to 0.13 MM PO4 g_l. The ratio of differences in phosphate released to differences in ambient phosphate within Figure 7. 39 Net reactive phosphate released from fine sediments under conditions of ambient phosphate accumulation (A-A), and cumulative phosphate released from fine sediments to ambient lake water replaced at 15 minute intervals (O-O). Lake water initially contained 0.10 uM P04 1‘1. Values are expressed as a mean of six replicates :1 standard error. ON. 0: 8. 305558386 .883 8.0.. (ESE—com to 522.5 cm on Oh cm on O? on d d I d u I l L U 08: m 0 u. A 0 3... W U. D N 08.. m m s O 0. con. I w w 32% A m . K o8. m 9 m. N 41 the 15 - 30 minute exposure period was: 0.38 - 0.13 5.10 - 0.10 "" 0°05 ' This ratio declined during subsequent exposure intervals, in— dicating effects of diminishing quantities of desorbable phos- phate associated with particles exposed to renewed medium. The cumulative total of phosphate released to renewed medium, 1.57 uM P04 9.1, was considered an estimate of de- sorbable particulate phosphate at time 0; and measured release during exposure intervals as a reduction in that amount. The data conformed to first order kinetics. The general form of phosphate exchange can be expressed as: particulate - PO solution - P04 4 Wtivfi‘ 2 where kl and k are exchange rate constants. The rate of part- 2 iculate -lake water phosphate exchange can be expressed by: dP _ _ jg — klP kZS where P is particulate P04 and S is solution P04. If solution phosphate is reduced to levels below those supporting steady state and maintained low through ambient replacement, then the rate of change of P in approaching equilibrium with phosphate- poor lake water may be approximated from: dP _ -a—t- — klP Applied here, the equation indicates that the rate of change of particulate desorbable phosphate, g; , is equal to a constant times the amount on the particle, and 42 —5 = -kldt . Integrating both sides the equation becomes a linear ex- pression, 1n P = -klt + P . The intercept P is the estimated desorbable particulate phos- phate at t = 0. Specifically, for the data presented here, 1) = -o.059 + 0.324 . ln(uM P04 9' The linearity of the data obtained (r2 = 0.995) implies that observed release rates were a function of particle content. The difference between the calculated intercept of 1.383 uM P04 g-l .324 (e = 1.383) and estimated 1.57 uM may reflect error involved in the implicit assumption that phosphate release within the initial 15 minutes interval was negligibly affected by ambient phosphate levels. The slope, k1 was an estimate of the rate of change of particulate phosphate and calculated as 0.059 min-1. While the release rate was specific to particle and environmental characteristics, the suggestion that it was a function of the amount of desorbable phosphate associated with particles exposed to a phosphate-poor environment may be k generally applicable when Ea is small. 1 DISCUSSION In contrast to overlying water, inorganic phosphorus frequently constitutes a major portion of total phosphorus in sediments (Sommers, et a1. 1970). Phosphate originating from mineralization of sedimented organic matter or flux from deep— er anaerobic sediment layers can result in interstitial phos- phate levels greater than those of overlying water. Relatively high interstitial phosphate can result in gradients, such as the one measured here. Because phosphate gradients represent potential for release of sediment phosphate to overlying water, the nature of mechanisms involved in developing and maintain- ing gradients deserves attention. Estimating the potential of phosphorus liberation from littoral sediments to overlying water can be approached by examining the effects of certain phosphate reactions between solid and solution sediment components. In addition to the impact of phosphate sorption/desorption reactions, equilibria characterized by solubility of solid Ca-phosphate precipitates can influence levels of interstitial phosphate, particularly in sediments of lakes that experience epilimnetic calcium pre- cipitation (Stumm and Leckie, 1971). Evidence of a midsummer decline in Lake Lansing littoral and epilimnetic calcium, sug— gests presence of calcium in Shallow sediments and warrants examination of the potential impact of Ca-phosphate solubility 43 44 equilibria on observed phosphorus levels in sediment solution. Since calcium phosphate equilibria expressions involve interactions between activities of specific ionic phosphate species and other soluble components of the aqueous medium, it is useful to consider factors influencing these variables with respect to the data presented. Soluble reactive phos- phorus is considered to exist as one or more of the forms of phosphoric acid, ie. H3PO4, H2P04—, HPO42-, and PO43-. The forms dominating soluble phosphate levels depends on pH. Under pH conditions of the work discussed here, concern is with H2PO4- and HPO42- (Figure 8). If soluble reactive phosphorus measure- ments used in this study are taken as estimates of ionic phos- phate levels, the amount of each form present can be determined from measurements of ambient pH. An additional level of con- sideration is drawn from observations of physical chemistry that the reactivity of a given ion may not be entirely determined by its molar concentration. The concept is contained in the Law of Mass Action which states that 'at a given temperature, the rate of a chemical reaction is proportional to the active masses of the reactants', indicating that activity of soluble phosphate species is a relevant parameter when considering chemical reactions. Concept of activity accounts for differences between observed and ideal behavior of solutes. When greater than 10’4 molar, the ionic strength of a given solution has been empirically determined to be an important factor modifying the effects of the concentration of the ion of interest. 45 Figure 8. Relation between pH and relative prOportions of soluble reactive phosphate species. Calculations based on dissociation constants of phosphoric acid, mono-, and dihydrogen phosphate ions in dilute solution. Ill ' P0 Hpo: H 2 Po;- HsPO4 46 I.0 0.9 t 0.7 ' I l l w. «o. a O O O ewquoud emooeu 0.3 ' ‘t o e .3 L N ". O O 0. 0 names IO uouoou wow I4 I3 l2 l0 pH 47 Quantification of this modifying factor, expressed in the activity coefficient, takes into account the concentration and valence of all ions in solution. For solutions of ionic strength less than 0.01 molar, the activity coefficient, y, can be calculated from: log Y = -A 22/5— where A is the temperature-dependent Debye-Huckel constant, equal to 0.5085 at 25°C, 2 is the valence of the ion of inter- est, and u is the ionic strength of the solution, u== %(.§ molar concentration of ion 1 times the valence of ion i l-l squared). Solute activity is equal to the concentration times the activity coefficient (Tinoco, et a1. 1978). Activity of soluble reactive phosphate is a function of pH dependent ionic distribution (Figure 8) and solution ionic strength (Figure 9). Determination of ionic strength is most accurately made with knowledge of all ions of solution. However, it can be approximated with measurements of dominant forms (Stumm and Morgan, 1981). Ionic strength was approximated here as 3 x molar calcium levels (after Lindsay and Moreno, 1960). Based on mea- surements of soluble calcium of Lake Lansing epilimnion, ionic strength was taken as 3 x 10-3molar. Griffin and Jurinak (1973), working with 124 river water samples and 27 soil extracts, found an excellent correlation (r2 = 0.992) between ionic strength and electrical conductance, u = 0.013 EC, where EC is expressed as millimhos cm-l. Electrical conductance of Lake Lansing epi- . . . . . -l . . . limnion varied around a mean of 0.26 mllllthS cm , Indicating 48 Figure 9. Influence of ionic strength (u) on activity coefficient (7) of mono- and dihydrogen phosphate for solutions of ionic strength less than 0.01 molar. V was calculated from: log y = -0.5085 22 /E— .49 I.000 r _ 3300' .700 " .600 l l 500 4.00 3.00 2.00 Ionic Strength 0f Solution (- I00 males .9." i 50 an ionic strength estimate near that made with calcium. Incorporating information presented above and appro- priate solubility expressions, the expected levels of total soluble reactive phosphorus in equilibrium with dicalcium phos- phate, octacalcium phosphate, and hydroxyapatite can be cal- culated (Figure 10; also of. Appendix 4). Calculations yield activities of soluble phosphate which have been converted to 4 concentrations (using assumptions presented above) for compari- son with data collected in this study. Although calcium phos- phates other than those presented in Figure 10 may form in sediments, lack of information concerning their solubilities precludes inclusion in the Figure. Fluoroapatite which gener- ally has a lower solubility than hydroxyapatite in the pH range of interest, has been excluded for lack of information concern- ing solution levels of fluoride ion. Significance of the solubility diagram of Figure 10 stems from the potential of sediment interstitial water phosphorus levels to approach equilibrium with the stable calcium phos- phate minerals considered. Measured soluble phosphate values falling above a given line indicate supersaturation with respect to that compound and suggest that precipitation is possible. Soluble phosphate levels below the line indicate undersaturation and that the corresponding compound would dissolve. Figure 10 suggests that dicalcium phosphate and octacalcium phosphate form or remain stable only at high levels of soluble phosphate. If soluble phosphate levels decreased, as a result of diffusion, biogenic uptake, or association with hydroxyapatite, these Figure 10. 51 Solubility diagram for certain calcium phosphate compounds in sediments at 250C. Levels of soluble reactive phosphate were calculated by considering appropriate solubility equilibria an assuming calcium levels maintained at l x 10‘ molar (solid lines) and 3 x 10'3 molar (dashed line). (°) refer to ambient reactive phosphate levels following 48 hours exposure to whole sediments to high initial levels of ambient phosphate. (X) same, but for fine sediment exposure. Phosphate activities were converted to concentrations for comparison with observed levels of total interstitial phosphate of south basin sediments. Interstitial Reactive Phosphate Concentration(-Ioo moles P041") 5 6 7 8 52 dicalcium phosphate dihydrate 00 H P04- 2 H20 octacalcium phosphate Ca4HI (20413-3 H20 hydroxyapatite CaDIPO‘Ié 0H I2 P b l b ‘01. Ambient pH 53 minerals would tend to dissolve. Using measured levels of to- tal interstitial phosphorus of sediments as an estimate of soluble reactive phosphate, surficial sediment solution levels ranged from about 0.581 to 3.871 uM 1-1, corresponding to 6.24 and 5.41, respectively, on the vertical scale. Based on the assumptions used to calculate the isotherms, sediments would ex- perience dissolution of any dicalcium or octacalcium phosphate present and precipitation of hydroxyapatite at pH levels above 7. Calcium increases above 1 x 10.3 molar would have the effect of shifting the lines parallel (to the ones indicated) and downward (cf. predicted hydroxyapatite isotherm under conditions of 3 x 10"3 molar calcium, dashed line of Figure 10). In general, solution pH, calcium, and ionic strength can greatly influence equilibrium phosphate levels. Observed variability in interstitial phosphorus could result from differences in the magnitude of these parameters over short distances within the sediment. Considerations of interstitial water buffering of solu- ble phosphate imposed by calcium phosphate equilibria suggest that equilibrium levels can be highly influenced by physico- chemical parameters. In the dynamic environment of Shallow water sediments, an important variable modifying sediment im- pact on soluble phosphate is the rate of response to ambient changes. In general, only the most soluble phosphate compounds react fast enough to determine solution phosphate activity un- der environmental conditions that may change within short time intervals. Kinetic information concerning the approach of 54 ambient phosphate to equilibria with calcium phosphate minerals in sediments is, at best, qualitative. Dicalcium phosphate is con— sidered the most reactive of the compounds discussed here, with reaction kinetics suggesting that equilibrium could be obtained within hours of ambient phosphate perturbation (Moreno, et a1. 1960). Hydroxyapatite is the most insoluble of the forms presented. The pH and temperature dependent rate of formation involves slow transformation of amorphous calcium phosphate into crystalline form. Reported rates of apatite formation vary from two weeks (Stumm and Leckie, 1971) to over three months (Moreno, et a1. 1960). It appears that at pH levels greater than 7, apatite phosphate may influence the direction and endpoint of ambient phosphate changes. However, attainment of equilibrium levels will likely be restricted to largely undisturbed systems where slow dissolution and precip- itation reactions have adequate time to dominate ambient phos- phate levels. Laboratory work reported here shows that under aerobic conditions a major portion of sediment-induced change from initial phosphate levels occurred within 30 minutes of expo- sure. The rate of change following 2 hours exposure to all initial phosphate levels examined was greatly reduced. Quasi- equilibrium phosphate levels observed in the laboratory follow- ing 2 and 48 hours exposure to high initial levels of ambient phosphate are plotted at the corresponding mean final pH (Figure 10). Solution ionic strength, calcium, and temperature in experiments were similar to those assumed when calculating 55 the isotherms. While a few points fall near predicted phos- phate levels in equilibrium with octacalcium phosphate, viewed as a set, the data appear to be independent of the isotherm. Moreno, et a1. (1960) observed that approximately one month was required for formation of the mineral under a similar pH and temperature regime. In the laboratory work presented here, calcium and phosphate activities did not exceed the solubility product of the more responsive dicalcium phosphate mineral. These observations indicate that 2 and 48 hour ambient phosphate concentrations stabilized at levels other than those predicted by solubility equilibria of calcium phosphate minerals consid- ered and implicate other controlling mechanisms. Sediment induced changes in ambient reactive phosphate have been attributed to phosphate association with pre-existing solids rather than through formation of insoluble precipitates. The Langmuir adsorption isotherm has been commonly used to quantitatively describe particulate induced changes in ambient reactive phosphate for soil (cf. Ellis and Knezek, 1972) and sediment (cf. Syers, et a1 1973) systems. In contrast to sol- ubility equilibria, Langmuir analysis describes the net effect of participating particulate-solution phosphate associations. The Langmuir equation was originally based on the ki- netic theory of gases. The same equation has been applied to adsorption of ions from solution by solids (Olsen and Watanabe, 1957). The analysis stems from the general consideration that observed solution levels result from the net effect of desorp- tion and sorption reactions acting Simultaneously (Shapiro and 56 Fried, 1959). Individually, the rate of release of phosphate from sediments is a function of sorbed particulate phosphate, X m ’1) rate of release = kl(m where k1 is a constant. Rate of phosphate sorption on parti- culate surfaces is a function of solution phosphate activity, (P), and the difference between phosphate sorbed and a sediment- specific sorption maximum, B: rate of sorption = k2(P)(B-%) . At equilibrium, the rate of release equals the rate of sorption: kl(%) = 1020:) (Ia-g) . A linear transformation of the equilibrium expressions yields: (P) = 1 + (P) x KB B m k2 where K is a sediment-specific constant, equal to E— , and is 1 related to the binding of particulate surfaces for sorbed phos- phate. The Langmuir constants, B and K, can be used to quanti- tatively describe the relationship between sorbed and equilib- rium ambient phosphate for a given sediment-solution system. Values for the constants(;?n be obtained from experimental ob- servations if a plot of x versus (P) is linear. The Slope m 1 and intercept are equal to l and respectively. i. B KB ' Since E-of the Langmuir relationship refers to amount of phosphate sorbed, measurements of solution phosphate losses Should be adjusted for the amount sorbed prior to sediment-solu- tion exposure. In laboratory work with whole sediments, this 57 adjustment was taken as the amount released following long-term exposure to phosphorus-poor lake water (0.305 uM P04 9.1). This was added to measured solution phosphate losses to estimate %-. Corresponding ambient phosphate activity was estimated by as- suming ionic strength to be 3 x 10-3 molar (ie. 3 x observed soluble calcium) and pH of 7.6 . Under these conditions, mea- sured concentrations of soluble reactive phosphate can be con- verted to activities by multiplying by 0.82 . A plot of values obtained following long-term whole sediment exposure to the three highest levels of initial ambient phosphate resulted in a straight line with a slope of 0.058, B = 17.179 uM P04 9.1; and an intercept of 2.125, K = 0.027 liters per uM P04. Estimation of the constants K and B for fine sediment material was made by adding 1.38 uM P04 9.1, the quanitity of desorbable phosphate estimated from analysis of P04 released to renewed ambient medium, to measured solution phosphate losses. At an ambient pH of 7.75 and ionic strength of 3 x 10"3 molar, the conversion factor used to obtain ambient phosphate activity was 0.81 . The constant B, indicating sorption maximum, was estimated as 47.091 uM PO4 g-1. The constant K, related to phosphate binding energy, was 0.002 liters per uM P04. The sorption maximum of fine material was over twice that calculated for whole sediments, and may reflect the greater proportion of small particles ( < 227 pm in diameter). The larger surface area per unit weight of material implies a greater number of sites for phosphate sorption. An interesting feature of the analysis is the apparent ten fold reduction in the value 58 of the constant related to the binding energy of processed fine sediments. This may indicate that expected carbonate precipi- tation on particulate surfaces during preparation of fine sedi- ments altered phosphate sorption characteristics. Green, et a1. (1978) found a negative correlation between value of the sorption energy parameter and calcite content of Maumee River sediments. Figure 11 illustrates the relationship between sorbed phosphate and ambient activity, based on experimental estimates of Langmuir constants B and K. The sediment-specific relation- ships suggest that particulate-solution phosphate exchange would result from perturbation of equilibrium ambient phosphate act- ivity (resulting from infusion of epilimnetic water through surficial sediments, for example). The direction and extent of particulate response would depend on phosphate activity of infusing medium, and sediment weight/solution volume ratio. An equation that quantitatively describes sediment re- sponse to perturbation from equilibrium ambient phosphate can be obtained from the Langmuir expression: flitL-I-fl x_ KB B m _1 1 - B ( R + P) 1 — 1 l E‘s—p‘rfl’) m x = B(P) m I E + (P) (Equation 1). If y is taken as the net quantity of phosphate exchanged with sediment, then g-- y is the level of sorbed sediment Figure 11. 59 Calculated relationship between phosphate sorbed by whole and fine sediments and ambient phosphate activity. Each isotherm was developed from esti- mations of Langmuir constants B and K. Ionic strength 3 x 10'3 M and equilibrium with respect to particulate—solution phosphate exchange were assumed. The pH was taken as 7.60 and 7.75 for whole and fine sediments, respectively. 60 ON ...... 40.. .23 3.264 22.32... 2.88.... .835 8.8. 8.2. 8. o... 9.8. 8. o: 8. 8 8 2 8 8 0.. on 8 o. q q 1 1 q u u q 1 d d 1 222.68 2.... 382.com 22.3 O) O h 10 D V '0 N (I-D’Od WV) (39qu aloudsoud mucosa enmzo ..“2 61 phosphate following exchange. Ambient phosphate in equilibrium with g'- y will depend on impact of phosphate mass exchanged on solution activity: B(P1 + AP) E. = m Y 1 ... (P1 + AP) K (Equation 2) where P1 is the phosphate activity of infusing medium, and AP is a function of the solid : solution following perturbation and a conversion factor relating phosphate exchanged to result- ant ambient phosphate activity. Laboratory results obtained with whole sediments can be used to illustrate an application of these considerations. Estimated phosphate activity at crossover, 2.624 uM P04 1-1 (ie. 3.20 x 0.82 = 2.624) reflected an estimated sorbed phos- phate level, %" of 1.137 uM P04 9-1 (from Equation 1). Per- turbation of the sediment system by introduction to phosphate- poor lake water (P1 = 0.13 x 0.82 = 0.107 uM PO 1-1) resulted 4 in net phosphate release from sediments. At the solid : solu- tion of 0.8 g dry weight per 100 m1 lake water, each uM P04 released per gram would increase ambient phosphate activity by 6.56 uM PO 1'1 (ie. 8 um PO 1‘1 x 0.82 = 6.56). For this ex- 4 4 ample, AP of Equation 2 is equal to 6.56 g 2-1 times y; where y is expressed in uM P04 9.1. Expected weight-normalized phos- phate release, y, can be calculated from: = 17.179 (0.107 + 6.56y) 37.037 + (0.107 + 6.56y) ' 1 1.137 - Expected phosphate release is 0.28 uM P04 9- In general, calculated net phosphate exchange resulting 62 from whole sediment exposure to initial phosphate activities examined agree reasonably well with observations of net exchange (Table 1). Employing analogous considerations to results obtained with fine sediment, calculated net release following exposure to initial phosphate activities used is somewhat less than that observed (Table 2). This may reflect error in the estimated crossover value since points used for its estimation were widely separated (cf. Figure 2). Influence of solid : solution on net phosphate exchange is indicated from measurements of phosphate desorption from 0.5 9 fine sediment under conditions of ambient medium renewal. In this experiment, 100 m1 lake water was renewed six times and considered similar to a solid : solution of 0.83 g 1-1. Assum- ing pre-disturbed ambient activity of 18.63 uM PO4 1-1, the calculated net release of 1.58 uM P04 9.1 is comparable to measured cumulative release of 1.57 uM P04 9.1. Similarly, calculated desorption from 0.5 g to 100 m1 lake water is 1.22 pM P04 9.1 and comparable to observed release of 1.26 uM. A few limitations of the preceded analysis should be noted. Care should be exercised when applying constants. Any experimental error experienced when obtaining data could lead to serious deviations between calculated and operational values of B and K. Further, an implicit assumption of Equation 2 is that phosphate sorbed is entirely releaseable. This may not be true. If sorbed phosphate forms metastable associations with the particulate sorbent, hydrated ferric oxides for example, 63 Table 1. Calculated and measured net phosphate exchange (uM P0 9.1) between whole sediments and infusing lake water of different initial phosphate activity (uM P04 1'1). Values assume equilibrium conditions with ionic strength 3 x 10'3M and pH 7.60 . Negative sign indicates net sediment sorption of phosphate. Initial Phosphate Calculated Measured Activity Exchange Exchange 0.107 0.28 0.30 1.246 0.15 0.16 -0.27 -0.31 5.108 64 Table 2. Calculated and measured net phosphate exchange (uM P04 g'l) between fine sediments and infusing lake water of different initial phosphate activity (uM P04 1'1). Calculated exchange is based on experimentally estimated ambient crossover activity of 18.63 uM P04 1'1 (A), and on a hypothetical crossover activity of 22 uM P04 1' (B). Values uilibrium conditions with ionic strength assume gq 3 x 10' M and pH 7.75. Initial Calculated Exchange - Phosphate Measured Activity A B Exchange 0.057 0.96 1.14 1.18 1.284 0.90 1.07 1.08 4.965 0.71 0.87 0.82 65 changes in the sorbent over time may affect its ability to re- lease phosphate following perturbation. While these consider- ations suggest that modifications of Equation 2 may be required for quantitative predictions, the expression may be used to qualitatively envision response of a sediment-solution system to ambient perturbation. A generalized representation of sorption/desorption con- cepts discussed is presented in Figure 12. Based on observa- tions of this study, the figure illustrates hypothetical interr actions between sediment particulates, characterized by phosphate binding energy (K) and sorption capacity (B), and infusing media. The direction and extent of phosphate exchange depends on equilibrium interstitial phosphate activity (Afluo¢ .H om.vH mH.H mo.m no.m oo.~m oo.H H~.H Ho.a oa.~m v>.o om.o Hm.o I wH.o I oo.o mm.a oo.v mm.o mm.v~ ho.m oo.¢a mN.o vm.v Hm.H oH.v~ mH.o m¢.N Hm.o t no.0 i oo.o mm.o «om say a 403 2: son 22v 402 21c mmmouocH m pmuommoz ta vom 23v coflumnucmocoo :oflumuucwocou w CH wmmmHOCH mfium m m< H.mc0wumummmum Avommmxv pumocmum ca mumnmmozm mo mcoHumcfiE unmump so Ao~m5.v0m¢mmmz mm cowboy mumcmmum mo mcoflufioom pmuomwme mo Doommm .Hi< magma 76 Figure A-l. Percent over-estimation of phosphate concentration as a function of the ratio of phosphate to arsenate in solution. 77 ¢N mm mm .N ON m. h. w. .....02 237.16.. 5:. n. V. n. N. 0. mo» #IO ‘1' no d J [ml 1 a u 0. ON on o¢ on ow ON 00 om uo-ginnueauog owquoqd JO uogwwusadano % 78 Table A-2. Effects of addition of arsenate reductant on determinations of phosphate in distilled water solutions containing measured additions of phosphate and arsenate. Measured P04 Concentration Concentration (M 1-1) (M 1'1) w/o with 1 P04 AsO4 Reductant Reductant 0.00 0.00 0.00 0.00 0.33 0.03 0.00 1.33 0.16 0.00 0.81 0.00 0.84 0.74 0.33 1.00 0.77 1.33 1.55 0.77 1.61 0.00 1.61 1.58 0.33 1.84 1.58 1.33 2.68 1.71 8.07 0.00 8.07 8.07 0.33 8.39 8.07 1.33 9.29 8.03 l. 5 ml of arsenate reductant containing Na28205 and Na252O3-5H20 were added to samples and standards (cf. Johnson, 1971). 79 Appendix 2. Estimates of phosphorus, iron, and calcium added to fine sediment particulates during the evaporation and drying treatment. 80 .Hzon H How Do omm um cowumoneoo mcfl30a IHom umoH mm3 unmfim3 who unmfiflpmm moan mo wmm mamumE«xoumm£ .usmflm3 wup unmafipmw maons mo wow mmz AsmumEmao a“ E: bmmvv cofluomum unmaflomm mcflm mo unmflmz who .H mmomfl made some whose was ssfloamo Hmuoe nm.m~ «m.m mo.m mw.mH am.o cosH Hmuoe vm.aa mm.oa ¢H.m m~.> me.o macromona sauce oa.m mm.m am.a ~¢.~ GH.o mumnmmosm o>fluommm 121V 12:3 1 na 21c Axuc land 21V ucmnoosoo mansflom Hmuoa ucmEfipmm no omHHou Hmumz mxmq .HHMHHODME newswpmm mo 6 mm mo AUoom umv mcflhuo mCMMSU wmuflmommw ucmcomeoo comm mo unsoEm may mo mmumeflumm uoommummu mosam> mcflccommwunoo .ucmefloom mcflm mcfimmmooum :a com: Hmum3 mme can mcofluomaaoo ucmfifipmm 3oHHmnm cfl Edfloamo com .couw .mumnmmonm Hmuou pom m>wuommn OHQSHOm mo cowumnuomocou .m|¢ canoe 81 Appendix 3. Total iron content of filtrate following exposure of sediment particulates to different initial levels of ambient phosphate. 82 Table A-4. Filterable iron (uM 1'1) in lake water of different initial phosphate levels following 0.5, 2.0, and 48 hours exposure to whole sedi- mentsl. Lake water, prin to sediment addition, contained 0.734 uM Fe 1' . Initial Lake Water Duration of Exposure (hours) ) Phosphate (UM 1’ 0.5 2.0 48 0.13 0.895 0.985 0.985 1.52 0.698 0.985 0.842 6.23 0.734 0.842 0.842 31.29 0.734 0.842 0.806 122.30 0.734 0.842 0.734 302.90 0.859 0.895 0.698 1. Estimated contribution of filterable iron from aqueous of whole sediment is 0.1 uM Fe 1'1. 83 Table A-5. Filterable iron (uM 1'1) in lake water of different initial phosphate levels following 0.5, 2.0, and 48 hours exposure to fine sedi- ments. Lake water prior to sediment addition contained 0.342 uM Fe 1‘1. Initial Lake Water Duration of Exposure (hours) ) Phosphate (UM 1' 0.5 2.0 48 0.07 0.985 0.734 0.842 1.58 0.913 0.985 0.734 6.13 0.842 0.842 0.985 32.26 0.842 0.734 0.842 122.30 0.985 0.842 0.734 302.60 0.967 0.842 0.842 84 Appendix 4. Solubility and dissociation expressions, and calculations used in the development of the calcium phosphate solubility diagram of Figure 10. 85 .omma .ocmuoz pom >mm©cflq .mo .H oo.¢H mom + mm mm Omm momma mm.NH vommm I ¢Omm + an mmm Imvomm mvmcmmocm :mmoup>£ocoz om.h wommmm I wanna + mm mmm Iwommm mumsmmonm :omoupxcflo ~H.~ eommmm I wommmm + no Hem eommm neon oflnonmnoee oe.nHH zoom + «anew + nosed mm mimoceieoneoano whennensxonesm Hm.m¢ wommm + mm + come Mm Ommm.mavomvmvmo momnmmozm Esfioamomuoo om.e vomen + mom mm o~m~.¢ommno whenesefle mumnmmonm EDHOHMOHQ axm :oflmmmsmxm auflawnoaom manfiuom HMOflEmsu pcsomEou .Hoa wusmfim mo Emummflo muHHHQSHOm on» no unmEQon>mp may :a com: 00mm um nouns cam mmwoomm owom canonmmosm mo mucmumcoo coaumfi00mmwp cam mpcoomaoo mumnmmonm Esfloamo mo muosoonm wuflaanoaom .mIm manna 86 CALCIUM PHOSPHATE SOLUBILITY CALCULATIONS The following calculations of soluble reactive phosphate concentration in equilibrium with calcium phosphate compounds at pH 7 and 25°C illustrate the development of Figure 10. Here, calcium concentration was taken as 1 x 10"3 molar, and ionic strength as 3 x 10"3 M. Activity coefficients, y , were calculated from: log y = -A zz/u where A is a Debye-Huckel constant equal to 0.5085 at 250C, and z is the charge of the ion of interest. Accordingly, yCa2+ = 0.774, ynzpo4’ = 0.938, yHPo42' = 0.774, and ypo43‘ = 0.561. 'p' indicates negative log of the activity. For example, if Ca2+ concentration = l x 10.3 M, Ca2+ activity = 1 x 10"3 (0.774) = 7.74 x 10 '4; and pCa = 3.11. A. Dicalcium phosphate dihydrate - CaHPO4-2HZO 6.56 = pCa + pHPO4 6.56 = 3.11 + pHPO4 3.45 = pHPO4 . (of. Table A-6) Equilibrium levels of H2PO4- in solution at pH 7 can be obtained from: 7.20 = pH + pHPO4 - pH2PO4 7.20 = 7.00 + 3.45 - pH2P04 3.25 = pH2PO4 2- _ _ -4 2- -l - _ pHPO - 3.45 - 3.55 x 10 M HPO 1 , and pH PO - 4 4 2 4 3.25 = 5.62 x 10"4 M H p0 ‘ 1'1. 2 4 87 Molar concentration x y = activity; activity/y = molar concen- tration. inn activity (M 1-1) .1. concentration (M 1-1) HPO42’ 3.55 x 10-4 0.774 4.59 x 10-4 H2P04- 5.62 x 10‘4 0.938 5.99 x 10"4 p(reactive phosphate concentration) = 3(3é99 x 10.4 + 4.59 x 10-4) B. Octacalcium phosphate - Ca4H(PO4)3-3H20 46.91 + 4pCa + pH + 3pPO4 46.91 = 4(3.11) + 7.00 + 3ppo4 9.16 = pPO4 . Equlibrium levels of HPO42- can be obtained from: 12.32 = pH + pPO4 - pHPO4 12.32 = 7.00 + 9.16 - pHPO4 3.84 = pHPO4 . (HZPO4-) can be obtained from: 7.20 = 7.00 + 3.84 — szPO4 3.64 = szPO4 . ion activity (M 1-1) Y concentration (M 1-1) PO43‘ 6.92 x 10-10 0.561 1.23 x 10'9 HPO42' 1.45 x 10-4 0.774 1.87 x 10-4 H2P04- 2.29 x 10-4 0.938. 2.44 x 10"4 p(reactive phosphate concentration) = 3.26 . C. Hydroxyapatite - CalO(PO4)6(OH)2 ion activity (M 1-1) Y concentration (M 1-1) HPO42' 3.55 x 10"4 0.774 4.59 x 10'4 H2P04- 5.62 x 10‘4 0.938 5.99 x 10‘4 p( reactive phosphate concentration) p(5.99 x 10"4 + 4.59 x 10-4) 2.98 "1111111111"1111111111“