Effects of Ocean Particles on the Upwelling Radiance and Polarized Radiance in the Atmosphere–Ocean System

SHI Chong,WANG Pucai,Teruyuki NAKAJIMA,Yoshifumi OTA,TAN Saichun,and SHI Guangyu

1Key Laboratory of Middle Atmosphere and Global Environment Observation,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

2Earth Observation Research Center,Japan Aerospace Exploration Agency 2-1-1 Sengen,Tsukuba,Ibaraki 305-8505,Japan

3Japan Agency for Marine-Earth Science and Technology,3173-25 Showa-machi,Kanazawa-ku,Yokohama,Kanagawa 236-0001,Japan

4State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

5University of Chinese Academy of Sciences,Beijing 100049

Effects of Ocean Particles on the Upwelling Radiance and Polarized Radiance in the Atmosphere–Ocean System

SHI Chong1,5,WANG Pucai∗1,Teruyuki NAKAJIMA2,Yoshifumi OTA3,TAN Saichun4,and SHI Guangyu4

1Key Laboratory of Middle Atmosphere and Global Environment Observation,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

2Earth Observation Research Center,Japan Aerospace Exploration Agency 2-1-1 Sengen,Tsukuba,Ibaraki 305-8505,Japan

3Japan Agency for Marine-Earth Science and Technology,3173-25 Showa-machi,Kanazawa-ku,Yokohama,Kanagawa 236-0001,Japan

4State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

5University of Chinese Academy of Sciences,Beijing 100049

Based on a vector radiative transfer model of the atmosphere–ocean system,the influence of oceanic components on radiation processes,including polarization effects,was investigated in the wavelength region ranging from 0.380 to 0.865µm. The components considered were phytoplankton,inorganic suspended material(sediment),and colored,dissolved organic matter.Due to their important roles in oceanic radiation processes,the sensitivity of the bidirectional reflectance to the rough ocean surface,represented by the wind velocity 10 m above the ocean surface,and aerosol,were taken into account.The results demonstrated that both radiance and polarized radiance just below the ocean surface were sensitive to the change of the concentrations of the considered components,while the dependence of polarized radiance on the observation geometry was more sensitive than radiance.Significant differences in the specular plane existed between the impacts of the phytoplankton and sediment on the degree of polarization just above the ocean surface at 670 nm.At the top of the atmosphere(TOA), polarization was relatively insensitive to changing concentrations of ocean particles at longer wavelengths.Furthermore,the radiance at the TOA in the solar plane was more sensitive to the aerosol optical thickness than wind velocity.In contrast, wind velocity strongly influenced the radiance at the TOA in the sun glint region,while the polarization degree showed less dependence in that region.Finally,a nonlinear optimal inversion method was proposed to simultaneously retrieve the aerosol and wind velocity using radiance measurement.

ocean particles,atmosphere–ocean system,radiative transfer,polarization,ocean color

1.Introduction

Radiative transfer is one of the most important processes in the atmosphere–ocean system.It plays a critical role in the Earth’s energy budget and material exchange by driving thermodynamic processes(Shi,2007).Radiation is emitted, absorbed and scattered,including the inelastic scattering effect by components in the atmosphere–ocean system,such as molecular gases,aerosols,the sea surface,and ocean components(IOCCG,2006).Absorption reduces the intensity of radiation,and reflection changes the direction;scattering changes the intensity and directionality,and also introduces or modifies the polarization state(Duan et al.,2010).

Radiative transfer in the atmosphere–ocean system is comprised of non-uniformly refracting layered media due to different refractive index values;therefore,refraction by the sea surface exerts very important roles in radiation processes. Assuming a flat ocean surface,the Fresnel–Snell law(Jin and Stamnes,1994)can be used to calculate the reflectivity and transmissivity;however,in reality the ocean surface is windroughened and randomly oriented.This is relevant to wind velocity determining the variance of the density function of a wave slope obeying an isotropical Gaussian distribution.Cox and Munk(1954)estimated the statistical characteristics of the wave slope distribution through measurement of the sea surface;however,their calculations neglected the influence of surface roughness in the angular distribution of upwelling ra-diance and diffuse sunlight.Based on the results of Cox and Munk(1954),Nakajima and Tanaka(1983)further proposed a rough ocean surface model to calculate the reflectivity and transmissivity function,including the wave slope and shadowing effect.

Radiation processes in the ocean are influenced by four main components:pure seawater,phytoplankton,inorganic suspended material(sediment),and colored,dissolved organic matter(CDOM).In addition,temperature and salinity also have an effect on the absorption of water(R¨ottgers et al.,2010,2014).In general,ocean waters are classified into Case 1 or Case 2 waters.The optical properties of Case 1 waters are dominated by phytoplankton,while Case 2 waters are influenced not just by phytoplankton but also by sediment and/or CDOM(IOCCG,2000).Many studies have reported the impact of ocean components on radiation processes.Tanaka and Nakajima(1977)showed that ocean turbidity and the refractive index of hydrosols both affect the radiation processes in the atmosphere–ocean system.In Case 1 waters,observations and semi-analytic and OGCM model simulations have suggested that the upwelling radiance and variation of irradiance reflectance are influenced by the change of phytoplankton pigment concentration(Gordon et al.,1988;Morel and Maritorena,2001;Lin et al.,2007).In the case of polarized radiation,previous research has shown that the polarization signal is important for satellite retrieval of aerosol properties over the ocean(Mishchenko and Travis, 1997;Goloub et al.,1999);the polarized reflectance in the 490 nm band at the TOA remains insensitive to chlorophyll, even in mesotrophic waters,according to analyses of Polarization and Anisotropy of Reflectance for Atmospheric Sciences coupled with Observation from a Lidar(PARASOL) and Research Scanning Polarimeter(RSP)data(Harmel and Chami,2008;Chowdhary et al.,2012),and the use of parallel polarized radiance effectively reduces sun-glint interference and enhances the ocean color signal at the TOA(He et al.,2014).Moreover,some studies have also indicated that polarization information may be used in remote sensing (Takashima and Masuda,1985;Chowdhary,1999;Chowdhary et al.,2006;Chami,2007;Zhai et al.,2010;Hollstein and Fischer,2012a,2012b).

However,research on the radiation impact of Case 2 waters,aerosols and the rough ocean surface,especially polarization effects,is still rare.In addition,a comprehensive understanding of radiation processes in the atmosphere–ocean system also plays an important role in the remote sensing of ocean color or aerosol monitoring.In this study,the impacts of ocean particles,a rough ocean surface,and aerosols,on upwelling radiance and polarization radiance for Case 2 waters were investigated using a vector radiative transfer model known as Pstar(Ota et al.,2010).The Pstar model uses the discrete ordinate and matrix operator method and has been developed based on the N–T model(Nakajima and Tanaka, 1983,1986,1988;Nakajima et al.,2000).The ocean part of the Pstar model has been recently improved from Case 1 waters to Case 2 waters.

This paper first discusses the radiative transfer model and describes the modified ocean components in detail.Then, the effects of ocean particles on radiation processes,as well as polarization,are analyzed and discussed.Finally,the influence of the rough ocean surface and aerosol on the bidirectional reflectance and polarization degree at the TOA are discussed,and a nonlinear optimal method is also used to investigate the inversion of aerosol and wind velocity simultaneously.

2.Theoretical background

Due to the reflectance of the ocean surface,the source function in the RTE for the diffuse component of the Stokes vector is invoked by the direct solar beam reflected by the sea surface.Thus,the RTE is as follows:

3.Radiative transfer simulations

The atmosphere contains a mixture of molecules and aerosols.In the simulation,the gas absorption and scattering were taken into consideration.The 1976 United States standard atmosphere,with 30 layers,was used for the atmosphere profile,and a more realistic multi-component scattering approach was adopted in the aerosol model.In this approach, instead of using an average refractive index of a single component to represent the mixture of different kinds of aerosols, the refractive index of each type of aerosol particle was calculated by taking hygroscopic growth into account(Shettle and Fenn,1979;Yan et al.,2002).In the simulation,an external mixture aerosol model,comprising sea spray and tropospheric aerosols,was used(Shettle and Fenn,1979).

In the ocean body,we considered four optical propertyaltering components:pure seawater,phytoplankton,sediment,and CDOM,or yellow substance with a homogeneous vertical distribution.The absorption and scattering coefficient of pure seawater were taken from Pope and Fry(1997) and Hale and Querry(1973),respectively.Moreover,the absorption coefficient is influenced by the temperature and salinity according to R¨ottgers et al.[2010,Eq.(4)].The calculation of the phase function of seawater was treated using Rayleigh scattering theory with a depolarization factor of δ=0.039(R¨ottgers et al.,2010),due to their similar volume scattering function(Morel,1974).

Living phytoplankton strongly absorb visible light,and this is the principle component of Case 1 waters.The absorption coefficientofphytoplankton was computed using the empirical formula of Bricaud et al.(1995):

Here,[Chl]is the chlorophyll concentration(mg m−3),and A(λ)and B(λ)are positive,principal wavelength-dependent parameters.The scattering coefficient of phytoplankton bph(λ)was calculated using Huot et al.(2008),based on Morel and Maritorena(2001):

The scattering phase matrix of phytoplankton cells was calculated using Mie theory.The refractive index relative to seawater was 1.05 and the volume spectrum distribution was assumed to follow a Junge-4 distribution with a radius range of 0.1–50µm(Chami et al.,2001).

The absorption and scattering coefficient of sediment were calculated based on the bio-optical model adopted from Bowers et al.(1988).The scattering phase matrix was also calculated using Mie theory,and determined to have a refractive index of 1.2–0.001i.The size distribution also followed a Junge distribution(He et al.,2014).The yellow substance was treated as a pure absorber:it strongly absorbs blue light, and its absorption coefficient decreases rapidly with increasing wavelength.The wavelength-dependent absorption coefficient of yellow substance ays(λ)can be parameterized by an exponential empirical relation based on the absorption coefficient at 440 nm,ays(440),as proposed by Bricaud et al. (1998):

where Iuand Iupare the upwelling radiance and upwelling polarized radiance,Edis the downwelling irradiance,andθ andφare the viewing zenith and azimuth angles,respectively.This calculation was performed at three“altitudes”: just above the ocean surface,just below the ocean surface, and at the TOA.

4.Results and discussion

4.1.The standard underwater radiative transfer problem

The Pstar model has been shown to simulate radiation processes well,including polarization in the atmosphere (Kokhanovsky et al.,2010).With regard to oceanic radiation processes,the performance of the Pstar model was determined by comparing the radiative processes in the ocean with the standard radiative transfer problems defined by Mobley et al.(1993).Table 1 shows the results of the Pstar model and its comparison with the standard values for problem 5(Mobley et al.,1993),a check of performance of the ocean radiative transfer model considering chlorophyll,sea water and the rough ocean surface.It assumes that the water body is horizontally homogeneous and infinitely deep,with a single scattering albedo of 0.9,and a Petzold phase function(Petzold, 1972),truncated by the delta-M method(Wiscombe,1977). Before the radiative transfer calculation,the real part of the refractive index of water was 1.340.The ocean surface was assumed to be rough,with a wave slope standard deviation of 0.2,corresponding to a wind velocity of 7.23 m s−1,based on the Cox and Munk(1954)capillary-wave spectrum at a solar zenith angle of 80◦.The irradiance perpendicular to the Sun’srays just above the ocean surface was 1 W m−2nm−1,and the influence of the atmosphere was ignored.Good agreement between the Pstar and Mobley-derived average values was demonstrated.The average difference of radiance(1.6%) and irradiance(1%)were very small,indicating that the Pstar model simulated radiative processes in the ocean body well.

4.2.Effects of phytoplankton on radiation processes

A series of sensitivity experiments were performed to study the influence of ocean particles in the upwelling radiance,as well as polarization effects over the ruffled ocean surface.The simulation was performed at ten wavelengths (380,400,412,443,490,510,550,670,750 and 865 nm),i.e. the central wavelengths of the detecting bands of ocean color satellites;a constant salinity of 35 practical salinity units and a temperature of 12◦C was assumed.The solar zenith angle was fixed at 30◦.All results were specified in the principal plane.Atmospheric conditions were taken into account with three aerosol optical thicknesses(AOTs)ofτ=0.01,0.2 and 0.5,and three wind velocities at 10 m height over the ocean surface(2,5 and 8 m s−1).

The numerical results for Case 1 waters,where the nonwater component is dominated by phytoplankton,were investigated first.The concentrations of chlorophyll used in the simulations were 0.03,0.1,1.0,5.0,10.0 and 20.0 mg m−3. Table 2 summarizesthe absorption coefficientsofchlorophyll at different concentrations for two different wavelengths.As expected,aph(λ)increased with concentration and displayed a wavelength dependence:the absorption coefficient at 412 nm was almost twice as high as that at 670 nm.

Research on bidirectionalreflectance justbelow the ocean surface is valuable due to its close relationship with remote sensing.Figures 1a and b show the angular distribution ofρ(λ)for different chlorophyll concentrations at 412 nm and 670 nm,respectively.The results indicate thatρ(412) decreased as the concentration increased due to the strong absorption of blue light by chlorophyll;more chlorophyll added more of an absorption effect.At 490 nm,a pivot point appeared in whichρ(λ)was barely dependent on chlorophyll concentration(not shown).The bidirectional reflectance increased as the chlorophyll concentrations increased at 670 nm,owing to the weak absorption effect and significant scattering(Fig.1b).The distinct relationship between chlorophyll and upwelling radiance on both sides of the pivot point form a theoreticalbasisforthe ratio method ofchlorophyllremote sensing(Gordon and Morel,1983).It is noted thatρ(λ) showed a slight dependence on the geometry of observation just below the ocean surface and its variation with viewing zenith angle curved more uniformly for high chlorophyll concentrations.

The influence of chlorophyll on the upwelling polarization bidirectional reflectanceρp(λ)just below the ocean surface was also investigated(Fig.2).The overall variation pattern ofρp(λ)was similar to that ofρ(λ),except thatρp(λ) showed obvious dependence on the viewing zenith angle,especially in the large viewing zenith angles of the specular plane.The polarization bidirectional reflectance decreased with an increase in chlorophyll concentration at 412 nm due to the strong absorption of chlorophyll and depolarization of multiple scattering(Fig.2a).At 670 nm,ρp(λ)increased with an increase in oceanic turbidity due to the weaker absorption effect of chlorophyll and the stronger backward scattering of multiple scattering after the competition with its depolarization effects(Fig.2b).

Next,the dependence ofρ(λ)andρp(λ)just above the ocean surface on the chlorophyll concentration were investigated(Fig.3).The bidirectional reflectance was clearly sensitive to chlorophyll concentration,with similar rangeability among viewing zenith angles at 412 nm.It decreased with[Chl]increase,but a significant difference was that the change ofρp(λ)for different chlorophyll concentrations was only obvious in the highest viewing angles of the specular plane,with a maximum relative difference over 44%.Meanwhile,the polarization reflectance was fairly insensitive to [Chl]in the backscattering region,owing to less reflected skylight and a low degree of polarization.In general,the varying amplitude ofρ(λ)andρp(λ)just above the ocean surface,dependent on the chlorophyll concentrations,was smaller than that just below the ocean surface.

Figure 4 shows the angular variation of bidirectional reflectance and the degree of polarization at the TOA as a function of chlorophyll concentration.The radiance was slightlysensitive to the chlorophyll concentration,and decreased as the concentrations increased at 412 nm.There was a relative difference of 11.8%when[Chl]changed from 0.03 to 20 mg m−3at nadir zenith angle(Fig.4a),while the degree of polarization barely varied in the principal plane,except at the complementary angle relative to solar zenith(Fig.4b).Note that greater insensitivity was apparent at longer wavelengths due to more significant absorption of sea water(not shown). This demonstrated that total radiance and polarized radiance at the TOA were predominantly from the atmospheric layer; in particular,the Rayleigh scattering strongly polarized the diffuse radiance.

Table 1.Comparison of Pstar with standard values derived by Mobley et al.(1993)for problem 5(Edis downwelling plane irradiance;E0uis upwelling scalar irradiance;I0uis upward nadir radiance).

Table 2.Absorption coefficients(m−1)of chlorophyll at 412 and 670 nm.

4.3.Effects of sediment and yellow substance on radiation processes

In the coastal ocean or estuarine region,ocean particles are dominated by inorganic suspended material with higher refractive index than phytoplankton.In the sensitivity simulations,the concentrations of sediment were defined as 1.0, 5.0,10.0 and 20.0 mg l−1.The concentrations of yellow substance were represented by the ays(440)in Eq.(7),and was set to 0.01,0.1 and 1.0 m−1.

Figure 5 shows the angular variation ofρ(λ)and the degree of polarization just below the ocean surface for different sediment concentrations at 550 nm.ρ(λ)increased with increasing concentration,and the increasing amplification decreased from 1 mg l−1to 20 mg l−1.The degree of polarization,which is typically characterized by a“bell shape”with a neutral point in the total backscattering region,decreased with increasing sediment concentration owing to the depolarization of multiple scattering.It was also noted that the amplitude of variation was more obvious in the longer wavelengths(figure not shown).

The dependence of the radiance reflectance and polarization degree just above the ocean surface on the concentrations of sediment was similar to that just below the ocean surface.However,the response of the degree of polarization just above the ocean surface at 670 nm to the change of chlorophyll and sediment was obviously different(Fig.6). Similar observations were also made in the specular plane.A significant amount of sediment can induce strong depolarization,while the degree of polarization for chlorophyll changes slightly with the increase of concentration.Such a situation is also similar at the TOA,it has been demonstrated that the degree of polarization in the specular plane at 670 nm can potentially be used to retrieve the sediment due to its more significant sensitivity to sediment than chlorophyll (Chami and McKee,2007;Ibrahim et al.,2012).

The yellow substance only absorbed;therefore,increasing its concentration increased the absorption of radiation in the ocean body.The radiance and polarization reflectance just below the ocean surface decreased by a similar ratio(Fig.7), and the dependence ofρ(λ)andρp(λ)on concentration was more significant in the shorter wavelength region due to the stronger absorption coefficient.

Radiance just above the ocean surface was also more sensitive to the change of the concentrations of yellow substance when compared with the polarized radiance below 670 nm (figure not shown).For example,the variation ofρ(λ)was 54.8%,while it was about 22.6%forρp(λ)at 412 nm when the absorption coefficient of CDOM changed from 0.01 m−1to 0.1 m−1.

At the TOA,the radiance and polarization reflectance were insensitive to the variation in yellow substance concentration,since the main radiance and polarization effects came from the contribution of the atmospheric layer and ocean surface.

4.4.Effects of AOT and wind velocity on upwelling radiation

At longer wavelengths,radiation sensitivity to varying concentrations of oceanic components was smaller,due to the significant absorption effects of seawater,which allowed the derivation of aerosol concentration and wind velocity.A primary objective of atmospheric correction of satellite images is to extract the effects of aerosol on the upwelling radiance at the TOA.

In the following sensitivity experiments,the solar zenith angle was fixed at 30◦and the concentrations of ocean particles were set to 0.03 mg m−3,0.1 mg l−1and 0.01 m−1for phytoplankton,sediment and CDOM,respectively.The variations ofρ(λ)and degree of polarization at the TOA are shown in Fig.8.

The results of the specular plane direction simulation demonstrate that the bidirectional reflectanceρ(λ)increased with AOT when the viewing zenith angle was large.The opposite was true when the viewing zenith angle was small,especially in the sun-glint region(Fig.8a).This may be because,at larger viewing zenith angles,diffused sky radiation is reflected more by the atmospheric layer than direct solar radiation in the turbid atmosphere with high AOT,which increases as AOT increases.At smaller viewing zenith angles, the upwelling radiance mainly comes from the reflection of direct solar radiation,so the direct transmittance decreases when the atmosphere is more turbid with heavy aerosol.As for the degree of polarization,it mainly decreases with increasing AOT,due to the depolarization effect by increased scattering order,and this phenomenon is more significant in the larger viewing zenith angles.

In view of the wind velocity dependence of the radiance reflectance and degree of polarization on the viewing zenith angles(Fig.8),it is worth noting thatρ(λ)was the most sensitive to the wind velocity in the sun-glint region,suggesting that a robust signal might be used to retrieve the wind velocity.ρ(λ)decreased with increasing wind velocity in the sun-glint region,while the degree of polarization showed different variation patterns in that region.Relatively less dependence of the degree of polarization on the sun-glint region was found,indicating that the degree of polarization could effectively reduce sun-glint contamination(Fougnie et al.,1999;Zhou et al.,2013;He et al.,2014).Moreover,the pattern of“W-8-T-0.2”(i.e.a wind velocity of 8 m s−1and an AOT of 0.2),shown in Fig.8a,indicates that the curve of the upwelling radiance was less sharp and its variation with viewing zenith angle curves more uniform for stronger wind velocity.Consequently,the ocean surface tended to be a Lambert reflector as wind velocity increased.

Nevertheless,in the solar plane,the radiance and degree of polarization were almost unchanged with wind velocity,and their dependence on wind velocity was significantly smaller than AOT(Fig.8).Taking the viewing zenith angle of−30◦as an example,the results show that theρ(λ)for AOT=0.5(0.0876 sr−1)increased by～194.0%relative to that for AOT=0.05(0.0298 sr−1),while theρ(λ)only increased by～0.8%(from 0.0503 sr−1to 0.0507 sr−1)when the wind velocity increased from 2 to 8 m s−1.Therefore, the radiance in the back-scattered direction can be useful to retrieve aerosol information.

4.5.Simultaneous retrieval of AOT and wind velocity for CASE 1 Water

Figure 9 shows the results of the numerical simulation carried out to investigate the performance of the optimal method for retrieving aerosol and wind velocity.For all the retrievals,synthetic measurements(input values)for a given atmospheric state were generated using the radiative transfer model.Simulated retrievals(output values)were then performed based on optimal estimation theory.The statistics of the total degrees of freedom for signals(the trace of the averaging kernel matrix)were over 4.2,which was comparable to the number of state parameters and reveals that the retrieval was stable and not seriously affected by measurement noise. It was demonstrated that the relative errors in simultaneously determining the optical thickness of fine and sea salt particles were less than approximately 10%.The relative error for the inversion of wind velocity was less than 1.4%,due to an observation in the sun-glint direction.Therefore,a multi-angle satellite observation covering in-and out-of-sun-glint conditions is probably useful for the simultaneous retrieval of aerosol and wind velocity,which reflects the results reported by Harmel and Chami(2012)using polarized information.

5.Conclusion

This study used a radiative transfer model of the atmosphere–ocean system to investigate the influence of ocean components on upwelling radiance and polarized radiance.Simulations were performed at three distinct“altitudes”:just above and below the ocean surface,and at the TOA.The impact of the ruffled ocean surface,represented by the wind velocity 10 m above the ocean surface,was considered,as was that of aerosol.The polarization radiance showed an obvious dependence on the observation geometry for different chlorophyll concentrations,while radiance curved more uniformly in the higher chlorophyll concentration region,just below the ocean surface.The polarized radiance at the TOA was less sensitive to the influence of ocean particles than just below the ocean surface,due to depolarization by multiple scattering in the atmosphere.The upwelling radiance changed slightly for high sediment concentrations just below the ocean surface.At 670 nm,sediment showed a more significant variation pattern than chlorophyll for the degree of polarization.At longer wavelengths,the dependence of radiance or polarized radiance on the changing concentrations of ocean particles was smaller,due to the highly significant absorption effects of seawater.In the backscattered direction,the change of radiation with wind velocity was hardly variable,which is useful for the retrieval of aerosol.Meanwhile,in the sun-glint direction,radiance generated robust information by wind velocity,which could be used for the inversion of wind velocity,while the degree of polarization seemed to reduce the sun-glint contamination. Based on the sensitivity experiments,an inversion algorithm was constructed to simultaneously retrieve the aerosol and wind velocity from multi-wavelength,multi-directional radiance.Since we used a relatively simple aerosol model in the inversion,it would be important to use polarized information for the retrieval of aerosol,wind velocity and water-leaving radiance in future work.

Acknowledgements.This work was supported by the“Strategic Priority Research Program”ofthe Chinese Academy ofSciences (Grant No.XDA05100300),the National Basic Research Program of China(Grant No.2013CB955801),the National Natural Science Foundation of China(Grant Nos.41175030 and 41475136) and the National Basic Research Program of China(Grant No. 2014CB953703).One of coauthors was supported by funds from MOEJ/GOSAT&GOSAT2,JST/CREST/EMS/TEEDDA,JAXA/ EarthCARE&GCOM-C,MEXT/RECCA/SALSA,MEXT/Kakenhi/ Innovative Areas 2409,and MOEJ/ERTDF/S-12.The authors also wish to thank the anonymous reviewers for their relevant comments and suggestions.

REFERENCES

Bowers,D.G.,S.Boudjelas,and G.E.L.Marker,1998:The distribution of fine suspended sediments in the surface waters of the Irish Sea and its relation to tidal stirring.Int.J.Remote Sens.,19,2789–2805.

Bricaud,A.,M.Babin,A.Morel,and H.Claustre,1995:Variability in the chlorophyll-specific absorption coefficients of natural phytoplankton:Analysis and parameterization.J.Geophys.Res.,100,13 321–13 332.

Bricaud,A.,A.Morel,M.Babin,K.Allali,and H.Claustre, 1998:Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic(case 1)waters: Analysis and implications for bio-optical models.J.Geophys. Res.:Oceans(1978–2012),103,31 033–31 044.

Chami,M.,2007:Importance of the polarization in the retrieval of oceanic constituents from the remote sensing reflectance. J.Geophys.Res.:Oceans,112,doi:10.10292006JC003843.

Chami,M.,and D.McKee,2007:Determination of biogeochemical properties of marine particles using above water measurements of the degree of polarization at the Brewster angle.Optics Express,15,9494–9509.

Chami,M.,R.Santer,and E.Dilligeard,2001:Radiative transfer model for the computation of radiance and polarization in an ocean-atmosphere system:Polarization properties of suspended matter for remote sensing.Appl.Opt.,40,2398–2416.

Chowdhary J.,1999:Multiple scattering of polarized light in atmosphere–ocean systems:Application to sensitivity analysesofaerosolpolarimetry.Ph.D.thesis,Columbia University. 276 pp.

Chowdhary,J.,B.Cairns,and L.D.Travis,2006:Contributionof water-leaving radiances to multiangle,multispectral polarimetric observations over the open ocean:Bio-optical model results for case 1 waters.Appl.Opt.,45,5542–5567.

Chowdhary,J.,B.Cairns,F.Waquet,K.Knobelspiesse,M.Ottaviani,J.Redemann,L.Travis,and M.Mishchenko,2012: Sensitivity of multiangle,multispectral polarimetric remote sensing over open oceans to water-leaving radiance:Analyses of RSP data acquired during the MILAGRO campaign. Remote Sen.Environ.,118,284–308.

Cox,C.,and W.Munk,1954:Measurements of the roughness of the sea surface from photographs of the sun’s glitter.J.Opt. Soc.Am.,44,838–850.

Duan,M.Z.,Q.L.Min,and D.R.L¨u,2010:A polarized Radiative Transfer model based on successive order of scattering.Adv. Atmos.Sci.,27,891–900,doi:10.1007/s00376-009-9049-8.

Fougnie,B.,R.Frouin,P.Lecomte,and P.-Y.Deschamps,1999: Reduction of skylight reflection effects in the above-water measurement of diffuse marine reflectance.Appl.Opt.,38, 3844–3856.

Goloub,P.,D.Tanre,J.-L.Deuze,M.Herman,A.Marchand,and F.-M.Breon,1999:Validation of the first algorithm applied for deriving the aerosol properties over the ocean using the POLDER/ADEOS measurements.IEEE Trans.Geosci.Remote Sens.,37,1586–1596.

Gordon,H.R.,O.B.Brown,R.H.Evans,J.W.Brown,R.C. Smith,K.S.Baker,and D.K.Clark,1988:A semianalytic radiance model of ocean color.J.Geophys.Res.:Atmospheres,93,10 909–10 924.

Gordon,H.R.,and A.Y.Morel,1983:Remote assessment of ocean color for interpretation of satellite visible imagery:A review.Springer-Verlag,114 pp.

Hale,G.M.,and M.R.Querry,1973:Optical constants of water in the 200-nm to 200-microm wavelength region.Appl.Opt.,12,555–563.

Harmel,T.,and M.Chami,2008:Invariance of polarized reflectance measured at the top of atmosphere by PARASOL satellite instrument in the visible range with marine constituents in open ocean waters.Optics Express,16,6064–6080.

Harmel,T.,and M.Chami,2012:Determination of sea surface wind speed using the polarimetric and multidirectional properties of satellite measurements in visible bands.Geophys. Res.Lett.,39,L19611,doi:10.1029/2012GL053508.

He,X.Q.,D.L.Pan,Y.Bai,D.F.Wang,and Z.Z.Hao,2014: A new simple concept for ocean colour remote sensing using parallel polarisation radiance.Scientific Reports,4,10 pp.

Hollstein,A.,and J.Fischer,2012a:Effects of salinity,temperature,and polarization on top of atmosphere and water leaving radiances for case 1 waters.Appl.Opt.,51,8022–8033.

Hollstein,A.,and J.Fischer,2012b:Radiative transfer solutions for coupled atmosphere ocean systems using the matrix operator technique.Journal of Quantitative Spectroscopy and Radiative Transfer,113,536–548.

Huot,Y.,A.Morel,M.S.Twardowski,D.Stramski,and R. A.Reynolds,2008:Particle optical backscattering along a chlorophyll gradient in the upper layer of the eastern South Pacific Ocean.Biogeosciences,5,495–507.

Ibrahim,A.,A.Gilerson,T.Harmel,A.Tonizzo,J.Chowdhary, and S.Ahmed,2012:The relationship between upwelling underwater polarization and attenuation/absorption ratio.Optics Express,20,25 662–25 680.

IOCCG,2000:Remote Sensing of Ocean Colour in Coastal, and Other Optically-Complex,Waters.S.Sathyendranath, Eds.,Reports of the International Ocean-Colour Coordinating Group,No.3,IOCCG,Dartmouth,Canada,140 pp.

IOCCG,2006:Remote Sensing of Inherent Optical Properties: Fundamentals,Tests of Algorithms,and Applications.Reports of the International Ocean-Colour Coordinating Group, Z.Lee,Eds.,No.5,IOCCG,Dartmouth,Canada,126 pp.

Jin,Z.,and K.Stamnes,1994:Radiative transfer in nonuniformly refracting layered media:Atmosphere-ocean system.Appl. Opt.,33,431–442.

Kattawar,G.W.,G.N.Plass,and S.J.Hitzfelder,1976:Multiple scattered radiation emerging from Rayleigh and continental haze layers.1:Radiance,polarization,and neutral points. Appl.Opt.,15,632–647.

Kokhanovsky,A.A.,and Coauthors,2010:Benchmark results in vector atmospheric radiative transfer.Journal of Quantitative Spectroscopy and Radiative Transfer,111,1931–1946.

Lin,P.F.,H.L.Liu,and X.H.Zhang,2007:Sensitivity of the upper ocean temperature and circulation in the equatorial Pacific to solar radiation penetration due to phytoplankton.Adv. Atmos.Sci,24,765–780,doi:10.1007/s00376-007-0765-7.

Mishchenko,M.I.,and L.D.Travis,1997:Satellite retrieval of aerosol properties over the ocean using polarization as well as intensity of reflected sunlight.J.Geophys.Res.,102,16 989–17 013.

Mobley,C.D.,and Coauthors,1993:Comparison of numerical models for computing underwater light fields.Appl.Opt.,32, 7484–7504.

Morel,A.,1974:Optical properties of pure water and pure sea water.Optical Aspects of Oceanography,1,22 pp.

Morel,A.,and S.Maritorena,2001:Bio-optical properties of oceanic waters:A reappraisal.J.Geophys.Res.:Oceans,106, 7163–7180.

Nakajima,T.,and M.Tanaka,1983:Effect of wind–generated waves on the transfer of solar radiation in the atmosphereocean system.Journal of Quantitative Spectroscopy and Radiative Transfer,29,521–537.

Nakajima,T.,and M.Tanaka,1986:Matrix formulations for the transfer of solar radiation in a plane-parallel scattering atmosphere.Journal of Quantitative Spectroscopy and Radiative Transfer,35,13–21.

Nakajima,T.,and M.Tanaka,1988:Algorithms for radiative intensity calculations in moderately thick atmospheres using a truncation approximation.Journal of Quantitative Spectroscopy and Radiative Transfer,40,51–69.

Nakajima,T.,M.Tsukamoto,Y.Tsushima,A.Numaguti,and T. Kimura,2000:Modeling of the radiative process in an atmospheric general circulation model.Appl.Opt.,39,4869–4878.

Ota,Y.,A.Higurashi,T.Nakajima,and T.Yokota,2010:Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere-ocean system.Journal of Quantitative Spectroscopy and Radiative Transfer,111,878–894.

Petzold,T.J.,1972:Volume scattering functionsforselected ocean waters,DTIC Document.77 pp.

Plass,G.N.,G.K.Kattawar,and S.J.Hitzfelder,1976:Multiple scattered radiation emerging from Rayleigh and continental haze layer.2:Ellipti-city and direction of polarization.Appl. Opt.,15,1003–1011.

Pope,R.M.,and E.S.Fry,1997:Absorption spectrum(380–700 nm)of pure water.II.Integrating cavity measurements.Appl. Opt.,36,8710–8723.

Rodgers,C.D.,2000:Inverse methods for atmospheric sounding: Theory and Practice,Series on Atmospheric,Oceanic and Planetary Physics–Vol.2.World Scentic,Singapole,238 pp.

R¨ottgers,R.,R.Doerffer,D.McKee,and W.Sch¨onfeld,2010: Pure water spectral absorption,scattering,and real part of refractive index model.ESA algorithm technical basis document,20 pp.

R¨ottgers,R.,D.McKee,and C.Utschig,2014:Temperature and salinity correction coefficients for light absorption by water in the visible to infrared spectral region.Optics Express,22,25 093–25 108.

Shettle,E.P.,and R.W.Fenn,1979:Models for the aerosols of the lower atmosphere and the effects of humidity variations on their optical properties.DTIC Document,94 pp.

Shi,G.-Y.,2007:Atmospheric Radiation.Science Press,402 pp. (in Chinese)

Takashima,T.,and K.Masuda,1985:Degree of radiance and polarization of the upwelling radiation from an atmosphereocean system.Appl.Opt.,24,2423–2429.

Tanaka,M.,and T.Nakajima,1977:Effects of oceanic turbidity and index of refraction of hydrosols on the flux of solar radiation in the atmosphere-ocean system.Journal of Quantitative Spectroscopy and Radiative Transfer,18,93–111.

Wiscombe,W.J.,1977:The delta-M method:Rapid yet accurate radiative flux calculations for strongly asymmetric phase functions.J.Atmo.Sci.,36,1408–1422.

Yan,B.H.,K.Stamnes,W.Li,B.Q.Chen,J.J.Stamnes,and S.-C.Tsay,2002:Pitfalls in atmospheric correction of ocean colorimagery:How should aerosolopticalproperties be computed?Appl.Opt.,41,412–423.

Zhai,P.-W.,Y.X..Hu,J.Chowdhary,C.R.Trepte,P.L.Lucker, and D.B.Josset,2010:A vector radiative transfer model for coupled atmosphere and ocean systems with a rough interface.Journal of Quantitative Spectroscopy and Radiative Transfer,111,1025–1040.

Zhou,G.,W.Xu,C.Niu,and H.Zhao,2013:The polarization patterns of skylight reflected off wave water surface.Optics Express,21,32 549–32 565.

:Shi,C.,P.C.Wang,T.Nakajima,Y.Ota,S.-C.Tan,and G.-Y.Shi,2015:Effects of ocean particles on the upwelling radiance and polarized radiance in the atmosphere–ocean system.Adv.Atmos.Sci.,32(9),1186–1196,

10.1007/s00376-015-4222-8.

10 October 2014;revised 9 March 2015;accepted 25 March 2015)

∗Corresponding author:WANG Pucai

Email:pcwang@mail.iap.ac.cn