Scattering by pure seawater : Effect of salinity

A theoretical model was developed estimating the scattering by seawater that are due to concentration fluctuation. Combining with the model proposed for density fluctuation (Optics Express, 17, 1671, 2009), we evaluated the overall effect of sea salts on the scattering. The variation of seawater scattering with the salinity is a combination of two factors: decreasing contribution due to density fluctuation and increasing contribution due to concentration fluctuation, with the latter effect dominating. The trend is, however, slightly non-linear and the linear adjustment of scattering with salinity that is frequently used would lead to an underestimate by an average of 2%. The results estimated at S = 38.4‰ agree with the measurements by Morel (Cahiers Oceanographiques, 20, 157, 1968) with an average difference of 1%, well within his experimental error of 2%. The spectral signature also varies with salinity, with the power-law slope increasing from -4.286 to -4.306 for salinity from 0 to 40‰. 2009 Optical Society of America OCIS codes: (010.0010) Atmospheric and oceanic optics; (010.4450) Oceanic optics; (290.5840) Scattering, molecules. References and links 1. A. Einstein, "Theorie der Opaleszenz von homogenen Flüssigkeiten und Flüssigkeitsgemischen in der Nähe des kritischen Zustandes," Annalen der Physik 338, 1275-1298 (1910). 2. X. Zhang and L. Hu, "Estimating scattering of pure water from density fluctuation of the refractive index," Optics Express 17, 1671-1678 (2009). 3. A. Proutiere, E. Megnassan, and H. Hucteau, "Refractive index and density variations in pure liquids: A new theoretical relation," J. Phys. Chem. 96, 3485-3489 (1992). 4. R. S. Farinato, and R. L. Rowell, "New values of the light scattering depolarization and anisotropy of water," J. Chem. Phys. 65, 593-595 (1976). 5. A. Morel, "Etude Experimentale de la diffusion de la lumiere par l'eau, les solutions de chlorure de sodium et l'eau de mer optiquement pures," Journal de Chimie Physique 10, 1359-1366 (1966). 6. A. Morel, "Note au sujet des constantes de diffusion de la lumiere pour l'eau et l'eau de mer optiquement pures," Cahiers Oceanographiques 20, 157-162 (1968). 7. A. Morel and B. Gentili, "Diffuse reflectance of oceanic waters: its dependence on Sun angle as influenced by the molecular scattering contribution," Appl. Opt. 30, 4427-4438 (1991). 8. A. Morel, B. Gentili, H. Claustre, M. Babin, A. Bricaud, J. Ras, and F. Tieche, "Optical properties of the "clearest" natural waters," Limnol. Oceanogr. 52, 217-229 (2007). 9. M. S. Twardowski, H. Claustre, S. A. Freeman, D. Stramski, and Y. Huot, "Optical backscattering properties of the "clearest" natural waters," Biogeosciences 4, 1041-1058 (2007). 10. E. Boss, R. Collier, G. Larson, K. Fennel, and W. S. Pegau, "Measurements of spectral optical properties and their relation to biogeochemical variables and processes in Crater Lake, Crater Lake National Park, OR," Hydrobiologia 574, 149-159 (2007). 11. C. D. Mobley, Light and water: radiative transfer in natural waters (Academic Press, San Diego, 1994). 12. E. Boss and W. S. Pegau, "Relationship of light scattering at an angle in the backward direction to the backscattering coefficient," Appl. Opt. 40, 5503-5507 (2001). 13. A. Morel, "Optical properties of pure water and pure sea water," in Optical Aspects of Oceanography, N. G. Jerlov, and E. S. Nielsen, eds. (Academic Press, New York, 1974), pp. 1-24. 14. P. Debye, "Light Scattering in Solutions," J. Appl. Phys. 15, 338-342 (1944). 15. F. J. Millero, R. Feistel, D. G. Wright, and T. J. McDougall, "The composition of Standard Seawater and the definition of the Reference-Composition Salinity Scale," Deep-Sea Res. I 55, 50-72 (2008). 16. M. Jonasz, and G. R. Fournier, Light Scattering by Particles in Water: Theoretical and Experimental Foundations (Academic Press, New York, 2007). #107650 $15.00 USD Received 17 Feb 2009; revised 16 Mar 2009; accepted 21 Mar 2009; published 25 Mar 2009 (C) 2009 OSA 30 March 2009 / Vol. 17, No. 7 / OPTICS EXPRESS 5698 17. K. S. Shifrin, Physical Optics of Ocean Water (American Institute of Physics, New York, 1988). 18. X. Quan and E. S. Fry, "Empirical equation for the index of refraction of seawater," Appl. Opt. 34, 34773480 (1995). 19. B. A. Pethica and C. Smart, "Light scattering of electrolyte solutions," Trans. Faraday Soc. 62, 1890-1899 (1966). 20. F. J. Millero and W. H. Leung, "The thermodynamics of seawater at one atmosphere," Am J Sci 276, 10351077 (1976). 21. UNESCO, "Background papers and supporting data on the International Equation of State of Seawater 1980," UNESCO Technical Papers in Marine Science No. 38 (1981). 22. P. E. Ciddor, "Refractive index of air: new equations for the visible and near infrared," Appl. Opt. 35, 15661573 (1996). 23. F. J. Millero, C.-T. Chen, A. Bradshaw, and K. Schleicher, "A new high pressure equation of state for seawater," Deep-Sea Res. 27A, 255-264 (1980). 24. H. C. Brinkman and J. J. Hermans, "The Effect of Non-Homogeneity of Molecular Weight on the Scattering of Light by High Polymer Solutions," J. Chem. Phys. 17, 574-576 (1949). 25. J. G. Kirkwood and R. J. Goldberg, "Light Scattering Arising from Composition Fluctuations in MultiComponent Systems," J. Chem. Phys. 18, 54-57 (1950). 26. W. H. Stockmayer, "Light Scattering in Multi-Component Systems," J. Chem. Phys. 18, 58-61 (1950). 27. G. Oster, "The Scattering of Light and its Application to Chemistry," Chem. Rev. 43, 319-365 (1948). 28. F. J. Millero, "The Physical Chemistry of Seawater," Annu. Rev. Earth Planet. Sci. 2, 101-150 (1974). 29. M. R. Wright, An Introduction to Aqueous Electrolyte Solutions (John Wiley & Sons Ltd., West Sussex, 2007). 30. M. Lund, B. Jonsson, and T. Pedersen, "Activity coefficients in sea water using Monte Carlo simulations," Marine Chemistry 80, 95-101 (2003). 31. M. Kerker, The Scattering of Light and Other Electromagnetic Radiation (Academic Press, New York, 1969). 32. Y. Shiokawa and S. Suzuki, "Isotropic Rayleigh Light Scattering in Two-component Systems. The Interdependence of the Concentration and the Density Scattering," Bull. Chem. Soc. Japan 48, 171-174 (1975). 33. D. R. Lide, Handbook of Chemistry and Physics (CRC Press, Boca Raton, 2003). 34. K. S. Pitzer and G. Mayorga, "Thermodynamics of electrolytes. II. Activity and osmotic coefficients for strong electrolytes with one or both ions univalent," J. Phys. Chem. 77, 2300-2308 (1973). 35. P. D. T. Huibers, "Models for the wavelength dependence of the index of refraction of water," Appl. Opt. 36, 3785-3787 (1997). 36. I. L. Fabelinskii, Molecular Scattering of Light (Plenum Press, New York, 1968).


Introduction
Scattering by pure water or pure seawater is a quantity of fundamental importance in aquatic optics.The contribution by these hypothetically pure liquid, often used as a "blank", must be subtracted from the observed scattering in order to derive the properties of particles that are often of primary interest.According to Einstein-Smoluchowski theory [1], the scattering of light by a particle-free liquid is due to the microscopic fluctuation of its dielectric constant (ε, or equivalently, the refractive index, n, where ε = n 2 ), which in turn is caused by localized fluctuations in density and mixing ratio (i.e., concentration).Representing the scattering with the volume scattering function at 90 degree, β(90), then (90) (90) (90) where β d (90) and β c (90) are the scattering due to density and concentration fluctuation, respectively.
Recently we [2] evaluated theoretically the scattering by pure water, which is due to density fluctuation entirely.where λ is the wavelength, k the Boltzmann constant, and ρ, n, T, β T and f(δ) are the density, the refractive index in vacuum, the absolute temperature, the isothermal compressibility, and the Cabbanes factor of water, respectively.Table 1 lists the details of these parameters including their estimates and the associated errors.Theoretical calculation using the model that we proposed for estimating the density derivative of the refractive index ( ) n ρ ρ ∂ ∂ ) [3] agrees with Morel's measurements within the experiment error (2%) with an average difference ranging from -0.67% to 1.97% for values of δ varying from 0.039 to 0.051 [4].
Presence of sea salts induces additional fluctuations in the concentration leading to additional scattering.Morel [5,6] measured the scattering by natural seawater of salinity S = 38.4‰ at five wavelengths of 366, 405, 436, 546, and 578 nm and found an average increase of 30% over the pure water.To our best knowledge, there have been no other measurements reported evaluating the effect of sea salts on the scattering of seawater.While its contribution to the total scattering is generally small, up to 10% in the visible for very clear ocean waters, seawater can account for a sizable portion of total backscattering, up to 80% in the clearest waters [7][8][9], owing to its nearly isotropic characteristics of angular distribution.In alpine lakes, the contributions to the backscattering by water are also significant (30-60%) [10].
For the past 4 decades, the values observed by Morel at one salinity value have been used throughout global oceans and coastal waters [e.g., 11], though sometimes modified via a linear adjustment as a function of salinity [e.g., 12].Based on the interpretation in Morel [13], Boss and Pegau [12] proposed an empirical model for adjusting the seawater scattering, where β pw (90) is the scattering of pure water and S 0 =37.Given that the reported salinity for the natural seawater that Morel measured was 38.4‰, a value of 38.4 for S 0 may be more exact.
The linearity assumption is only valid when salinity approaching zero, namely, the solution is ideal, under the condition of which Morel [13] used a simple model to estimate the excessive scattering following Debye [14], where N A is the Avogadro number, M ss represents the mean molecular weight of the electrolyte, and S the mass concentration of solute with a unit of g/kg.For seawater, S in Eq. ( 4) differs from the Practical Salinity S (‰) in both definition and values, S(g/kg)=1.0047×S(‰)[15].But we will use them interchangeably, because the results will be barely affected.Equation ( 4) was also used by Jonasz and Fournier [16] and its variation with S having a unit of g/l by Shifrin [17] to estimate scattering due to sea salts.Note, the term 1/ρ in Eq. ( 4) was not present in Morel [13] or Jonasz and Fournier [16]; it is required either from theoretical deduction or from dimensional analysis.However, because ρ would take a unit of kg/l and hence have a value close to unity, their computations are still valid numerically.The term 2 n S ∂ ∂ in Eq. ( 4) can be assumed as a constant (within 0.6%) for a given temperature and wavelength [18], which leads to a prediction of linear variation of β c (90) with S. Apparently, if the effect of ρ is considered, the linearity between β c (90) and S will no longer be warranted.
The development of Eq. ( 4) was based on the assumption that the solution of two components is dilute and ideal.Even though seawater can be deemed as a dilute solution and the composition and proportion of its sea salts remains constant fairly well over global ocean, it is well known that an electrolyte solution such as seawater is not ideal.Morel [5] measured scattering by pure NaCl solution at 366 nm and 546 nm for concentrations from 0.01 to 0.1 g/ml, and his results at 546 nm in terms of percentage increase relative to the scattering by pure water are shown in Fig. 1.Obviously, the excessive scattering due to addition of NaCl varies with the concentration in a non-linear manner.The similar pattern of non-linear variations were also found in other simple electrolyte solutions of KCl, KNO 3 , or NaI [19].
Seawater is a much more complex system than a NaCl solution, even though NaCl is its major sea salt.While we know theoretically that sea salts enhance the scattering, which was also confirmed by Morel's measurement, we still do not know how exactly the scattering varies with salinity, which constitutes the purpose of this study.

Scattering due to sea salts
Seawater is a multi-component solution with dissolved components of sea salts either in their original form or as dissociated ions.Complex models developed for a multi-components system by Brinkman and Hermans [24], Kirkwood and Goldberg [25], Stockmayer [26] involved the coupling terms among any two solutes, which are difficult, if not impossible, to evaluate or measure.To explain the experimental measurements of light scattering by three simple electrolyte solutions of KCl, KNO 3 , and NaI, Pethica and Smart [19] used an equation by Oster [27], where c is salt concentration in kg/m 3 , 0 V and a 0 are the partial molar volume and the activity of the solvent, respectively.While still based on a two-component solution, Eq. ( 5) applies to the non-ideal solution with the correction term 0 ln a c ∂ ∂ .Strictly speaking, an electrolyte solution is a system of at least three components with dissociated ions and solvent.Therefore, the overall close agreement within 6% between the prediction of Eq. ( 5) and the measurements [19] over a wide range of concentrations suggests that Eq. ( 5) might be applicable to a multi-component solution like seawater as well.
In the appendix we showed that Eq. ( 5) would apply to seawater as long as the molar ratios among dissolved components of sea salts remain constant during the fluctuation of concentration or equivalently sea salts act like one hypothetical compound thermodynamically.This assumption was based on the following considerations.The Principle of Constant Proportion suggests that sea salts through global oceans can be regarded as one "sea salt" with major difference in their total content.In seawater, ions are paired through covalent bond, electrostatic bond, or solvent-shared electrostatic bond (separated by one or more water molecules) [28]; electrostatic interaction among ions impose some degree of order over random thermal motions [29].Actually, it is precisely the electrostatic interactions among ions that make seawater non-ideal [29].Also, Millero and Leung [20] have shown that many thermodynamic properties of seawater can be estimated theoretically by treating the solution as being a mixture of one hypothetical sea salt and water, and their results were in excellent agreement with the measurements.Moreover, the terms of 0 V and 0 ln a in Eq. ( 5) are all properties of the solvent; therefore at least part of the errors or uncertainties  4) Calc.w/ Eq. ( 5) Fig. 1.Scattering by NaCl solution measured by Morel [5] at 546 nm was scaled, following the power-law of a slope of -4.30, to the sodium D line wavelength, 589 nm, at which the thermodynamic parameters for NaCl solution are available.The measured values, shown as % increases relative to the pure water, are compared with theoretical calculations using Eq. ( 4) and Eq. ( 5).The theoretical values for pure water were calculated using Eq. ( 2).associated with this assumption regarding the thermodynamic behavior of sea salts may have already been accounted for in the estimates of the thermodynamic properties of their counterpart -water.
By definitions, the partial molar volume coefficient for water, 0 V and salinity S 0 0 Here we use the same definitions of variables as in Appendix.Deriving ( )  7) and taking ratio of them, we obtain, where M 0 is the molecular weight of water.The prediction of Eq. ( 8) agrees with the measurements of 0 V [20] within 0.05%.Replacing the partial derivative with respect to the sea salts concentration, c, with the partial derivative with respect to the salinity, S, using the relationship, c S ρ = , and inserting Eq. ( 8) to Eq. ( 5) lead to Equation ( 9) summarizes the effect of sea salts on the scattering that are due to concentration fluctuation.Equations ( 5) and ( 9) are the same; the advantage of Eq. ( 9) is that the refractive index, the density, and the activity for seawater have all been measured with high precision as a function of salinity.
The activity for seawater was recently estimated using a Monte Carlo simulation [30], the results agree with the experimental values [20] within 0.01%.Millero and Leung [20] fitted the measured values to a general Debye-Hückel form as a function of temperature and salinity.Combining Eqs. ( 14), ( 22), ( 23), (88), and (107) in Millero and Leung [20], we derived a formula for 0 ln a S ∂ ∂ , and the result is shown in Table 1, which also listed the formulae for the other parameters in Eqs ( 2) and ( 9).
If we assume that seawater is ideal, then the activity of seawater a 0 = (1-X ss ), where X ss is the molar fraction of sea salts and 0 ( ) . Combining this and Eq. ( 7) leads to Inserting Eq. ( 10) into Eq.( 9), we will have Eq.(4).By definition, M ss is the mean molecular weight of ions (dissociated sea salts) and non-dissociated sea salts dissolved in the solution.
For seawater, we estimated M ss = 31.33using the latest molar ratio data [15].For aqueous NaCl solution, Morel [13] used water.However, the theoretical mean ion weight for NaCl-type electrolyte should be where M 1 and M 2 represent the molecular weight of cation and anion, respectively.And mathematically, , with equality occurs for M 1 =M 2 , which seldom holds for natural salt.Similarly, Johnaz and Fournier [16] used to estimate the mean molecular weight of ions, where f i is the molar fraction of sea salts and i denotes different sea salt.They derived a value of 32.41, which is bigger than the value of 31.33 that we derived based on the molar fraction of ions directly.Mathematically, these two numbers will never be equal unless all the sea salts were of the same type, e.g., 1-1 like NaCl.Given these uncertainties associated with estimating the mean molecular weight for ions as well as Eq. ( 4) being only an approximation, we recommend the use of Eq. ( 9) for estimating seawater scattering that are due to concentration fluctuation.
The differences in estimating β c (90) for seawater with Eq. ( 4) and Eq. ( 9) are shown in Fig. 2. Overall, Eq. ( 4),which is only valid for an ideal solution, underestimates the scattering due to concentration fluctuation for salinity from 0 to 31‰ with the largest difference of -8.94% occurring at S=6.1‰ and overestimates at higher salinity values.The applicable range of Eq. ( 9), which was derived theoretically, is determined by the validity ranges of its parameters, most of which (shown in Table 1) were determined empirically as a function of temperature, salinity, and/or wavelength.The validity range for salinity is 0~40‰.The trend shown in Fig. 2 suggests that Eq. ( 4) would increasingly overestimate Eq. ( 9) for waters with extremely high salt content.This can also be expected from a theoretical ground because the higher the concentration of electrolytes, the larger the deviation a solution would be from ideal.
The presence of sea salts also affect the scattering that is due to density fluctuation, β d (90), because all the variables of n, ρ, and β T in Eq. ( 2) are functions of salinity (see Table 1).Since the fluctuations due to density and concentration are independent to each other [19,32], the total scattering of seawater is simply the sum of the two terms, β d (90) and β c (90), as in Eq. ( 1).We should stress here that the concentration fluctuation used throughout this study refers to the change in the mixing ratio [1], because otherwise, the variations in density and concentration will not be independent.

Results and discussion
For all the calculations, we use Eq. ( 2) to estimate the scattering due to density fluctuation with the term 2 ( ) n ρ ρ ∂ ∂ calculated using PMH model [2] listed in Table 1.The scattering due to concentration fluctuation will be calculated using either Eq. ( 4), (5), or (9).The total molecular scattering for a solution is calculated using Eq.(1).Apparently, for S = 0, the estimates will be the scattering by pure water.Also we use a value of 0.039 for the depolarization ratio parameter, δ [4].
Since NaCl is the dominant sea salt among the constituents of seawater, it is of interest to evaluate how the model performs as compared to the measurements by Morel [5].The values for the isothermal compressibility, the refractive index and the density of NaCl solution are from Handbook of Physics and Chemistry [33] and the values for the activity are from Pitzer and Mayorga [34].The prediction based on the assumption that an electrolyte solution is ideal (Eq.( 4), and green curve in Fig. 1) would overestimate significantly at higher concentrations.On the other hand, the estimate by Eq. ( 5) (red curve) can largely account for the observed variation over the entire range of NaCl concentration that had been measured, with an average agreement within 1.7%.The results shown in Fig. 1 provide us with confidence to proceed with seawater.
For seawater, we used Eq. ( 9) to estimate the scattering due to concentration fluctuation.Fig. 3 shows the spectral scattering estimated for S = 38.4‰,and T c = 20 °C and their comparison with the measurements by Morel [5,6].Also compared in Fig. 3 are Morel's theoretical estimates of scattering for seawater [Table 4 in 13], which has been widely used [e.g., 11].The theoretical prediction of our model agrees with Morel's measurements within 1.0%, well within the experimental error of 2%.It was recently suggested [8,9] that Morel's theoretical estimates (the green curve) might be high; and our theoretical calculation (red curve) confirmed the argument.
The scattering of seawater as a function of salinity calculated at 546 nm and 20°C are shown in Fig. 4. Clearly, the variation of scattering with the salinity is a combination of two factors: decreasing contribution due to density fluctuation (Eq.( 2)) and increasing contribution due to concentration fluctuation (Eq.( 9)).The effect of sea salts on the density term of scattering was briefly discussed by Morel [13], who suggested that the effect might be small due to contradictory influence of sea salts on β T (decreasing with S) and on n and ρ (both increasing with S).Our estimate (blue line in Fig. 4) confirms these early analyses and shows that the density term decreases by ~2.6% as S increases from 0 to 40‰.Given the high precisions associated with the estimates of the parameters used in calculating the density contribution (Eq.( 2) and Table 1), the magnitude of variation of 2.6%, though small, is statistically significant.
As can be expected, the scattering due to concentration fluctuation (dashed blue line in Fig. 4) increases with the salinity.Because of the dominant role of the term S in Eq. ( 9), the variation is primarily linear with small adjustments resulting from variations of ρ and 0 ln a S ∂ ∂ with S. The estimates based on the empirical model of Eq. ( 3) with S 0 =38.4 are also shown in Fig. 4 (red line).The comparison between the theoretical and empirical estimates suggests that the linear empirical model would underestimate the seawater scattering by an average of 2% for S up to 40‰ assuming the values for pure water are estimated correctly.
The minimum scattering occurs at 26°C, which is the same as that for pure water.And the range of variation with temperature increases slightly for seawater, 5.5% vs. 3.7% for pure water for T c from 0 to 26°C.Also, sea salts have a weak effect on the spectral dependence of scattering.With anchor wavelength λ 0 at 450 nm, the slope s as in 0 ( ) s λ λ determined from non-linear fitting varies from -4.286 to -4.306 for S from 0 to 40‰.These values are very close to Morel's estimate of -4.32.We used s = -4.30,roughly the mean, in Fig. 1 to scale the scattering of NaCl solution from 546 nm to 589 nm.Equation ( 9) was derived theoretically from thermodynamic principles, therefore its applicability is only limited by the ranges of validity of its parameters that are determined empirically (Table 1).The empirical model for the refractive index, though originally developed for the visible region [18], has been shown to fit the available data well over an extended range of 200 -1100 nm [35].The valid range for the temperature is 0 -30°C and for the salinity 0 -40‰.While these ranges cover most of oceanic and coastal waters, the formulae for density, refractive index, activity, and isothermal compressibility may not apply for inland seas, whose salt components might have different molar ratios among themselves from those for seawater.Also extrapolation beyond these applicable ranges should be verified with independent measurements.
As shown in Table 1, the relative errors in modeling the parameters used to estimate the seawater scattering are very small.We estimated that the relative error in the model itself is only 0.04%, which is primarily due to the term, 0 ln a S ∂ ∂ in Eq. ( 9).Therefore the accuracy of the model for the concentration fluctuation part is primarily determined by the assumption that we have made regarding the thermodynamic behavior of the dissolved components of sea salts -they maintain the relative molar ratio with each other during the fluctuation.The close agreements in scattering between the theoretical predictions by this study and the measurements by Morel [5,6] for NaCl solution and seawater as shown in Figs. 1 and 3, respectively, do provide justification.Recently, Lund et al. [30] estimated the activity coefficients in seawater using Monte Carlo simulations based on electric neutrality, which is similar to our assumption.They showed that their results are in excellent agreement with the experimental determinations and the estimates based on the Pitzer model.This also partly supports the validity of our assumption.However, if the assumption does not stand, we think it would lead to an overestimate of the contribution by sea salts, because if only part of the ions fluctuate while the rest remains still, the magnitude of concentration change would be smaller.
The other major source of uncertainty is the value adopted for the depolarization ratio.Sea salts are expected to modify the value of depolarization ratio of water with two contradictory effects: reducing δ because ions are isotropic and enhancing δ because their electrostatic field causes increased anisotropy [13,17].Both of the effects have been observed; for example, Pethica and Smart [19] measured the variations of δ with concentration that was increasing in a KNO 3 solution but decreasing in a KCl solution.We do not know yet how δ would vary with salinity and in this study we have assumed that δ remains unchanged.
In our previous study [2], we analyzed the effect of δ on the scattering by pure water and the results suggested that within the range of values of δ from 0.039 to 0.051 that have been measured [4], the model we proposed converges with the measurements by Morel [5,6] from -0.7% to 2%.Results in Fig. 3 were calculated with δ = 0.039 and the average difference is 1%; if δ = 0.051, the difference would be 3.7%.Given the experimental error of 2% reported for the measurements by Morel and assuming the same values of δ apply to seawater, the model developed seems to overestimate a little.This is consistent with our estimate of possible errors introduced by the assumption discussed above.However, in order to further constrain the uncertainty in seawater scattering, we need to improve our understanding of the effect of sea salts on the depolarization ratio.

Conclusions
We estimated the effect of sea salts on the scattering by seawater by developing a model calculating the scattering due to concentration fluctuation as a function of salinity.The results agreed with Morel's measurements with an average difference of 1%, which is within his experimental error of 2%.While in general, the scattering by seawater increases with salinity, the linear adjustment frequently used will cause an underestimate by 2% on average.The underlying assumption of the model was that relative molar ratios of components of dissolved sea salts remain constant with each other at microscopic scales during concentration fluctuation.As we have discussed and shown by the results, we think the impacts of the assumption are limited.Another factor contributing to the improved modeling is the recent availability of high quality estimates (either through measurement or modeling) of the key thermodynamic parameters, such as the activity and partial molar volume of seawater.To further test its applicability, the Matlab code of the model developed in this study is available to the community (http://www.und.edu/instruct/zhang/programs/betasw_ZHH2009.m).
where m denotes the total number of moles of electrolytes, M the molecular weight, and f the fractional abundance of different electrolytes in mole/mole ( 1 ).For the subscripts in Eq. (A2) and hereafter, ss denotes the quantities associated with sea salts or electrolytes in general, and i for one of the components of sea salts.Because the composition of sea salts and their concentration relative to each others are approximately constant over the global oceans, the parameters f i as well as the summation The assumption in Jonasz and Fournier [16] that all the sea salts are merged into one hypothetical compound can also lead to Eq. (A3).With this, the collective variation of dielectric constant due to sea salts can be written as: The concentration fluctuation, 2 c ∆ , results from changes in the mixing ratio of pure water and sea salts within a small volume of ∆V [1], while the total mass and hence the density remains constant.This ensures that fluctuations due to density and concentration are independent to each other.Thermodynamically, 2 where A is the Helmholtz free energy [36].The change in Helmholtz free energy associated with the concentration change at constant temperature and volume is given by, where µ denotes the chemical potential.Within a small volume, ∆V, the number of moles of pure water and sea salts are related by, where V denotes the partial molar volumes.Since ∆V is held as a constant, we have from Eq. (A7) that 0 0 The partial molar volumes depend on the concentration too, but the fluctuations that we are concerned with here are too small to be important.For example, from Eq. ( 8) we estimated that Differentiating Eq. (A10) with respect to c, we have 2 0 , 2 0 ( ) ( ) The differentials of chemical potentials for the components in a system are related by the Gibbs-Duhem equation, ) Recalling Eq. (A7) for ∆V and Eq.(A2) for c, the second term of Eq. ( A13) is effectively 1/c, which leads to a simplification of Eq. (A13) to µ is the chemical potential of the species in a standard state and R is the gas constant.
For seawater, the state can be the temperature and pressure under which the scattering is being measured or calculated.Differentiating Eq. (A15) with respect to the concentration of sea salts and substituting the result into Eq.(A14), we obtain, 2 0 , , 2 0 ln ( ) ( ) Though not appearing in Eqs.(A14) and (A16) explicitly, the effects of sea salts were carried through the changes in the chemical potential and activity of the solvent as the concentration fluctuates.Combination of Eqs.(A16), (A5), (A4) and (A1) leads to Eq. ( 5), which was originally developed for a two-component solution.We have shown that it can be used for seawater too, under the assumption that sea salts would fluctuate in synchronization with each other such that their molar ratios among each other remain unchanged.
$15.00 USD Received 17 Feb 2009; revised 16 Mar 2009; accepted 21 Mar 2009; published 25 Mar 2009 (C) 2009 OSA b T b T b T b T c c T c T c T c T S c c T c T S c S b T b T b T c c T c T c T S c c T c T

a
The numbers in parentheses are the relative errors of estimates in percentage b Tc (0 ~ 30) is the temperature in Degree Celsius, S (0 ~ 40) is the salinity in ‰, λ (220 -1100 nm) is the wavelength in nm cThe index, n, here is defined with respect to the air and therefore it has to be scaled by the refractive index of the air, na, to obtain the index in vacuum.#107650 -$15.00USD Received 17 Feb 2009; revised 16 Mar 2009; accepted 21 Mar 2009; published 25 Mar 2009 (C) 2009 OSA MeasuredCalc.w/ Eq. ( #107650 -$15.00USD Received 17 Feb 2009; revised 16 Mar 2009; accepted 21 Mar 2009; published 25 Mar 2009 (C) 2009 OSA

η
to approximate its value, where M NaCl is the molecular weight of NaCl and η=2, which is number of ions into which NaCl dissociates in #107650 -$15.00USD Received 17 Feb 2009; revised 16 Mar 2009; accepted 21 Mar 2009; published 25 Mar 2009 (C) 2009 OSA

Fig. 4 .
Fig. 4. Scattering of seawater due to density fluctuation (blue line), due to concentration fluctuation (blue dashed line and right Y-axis) and their total are estimated for λ = 546 nm and Tc = 20°C as a function of salinity.The empirical variation estimated using Eq.(3) is shown as red.
(A2) can be assume to be constant.Let M ss denote i i f M ∑.If we further assume that on a microscopic scale, f i also remains constant with each other during fluctuation, then the concentration fluctuation, dc, can safely assume that the partial molar volumes are constant during the concentration fluctuations (∆m ss ).Combining Eqs.(A6) and (A8), we obtain,

Table 1 .
Parameters used in estimating seawater scattering and their formulae for calculation