search

CNES / RadiativeTransferCode-OSOAA

Ocean Successive Orders with Atmosphere - Advanced (OSOAA) : Accurate computation of the radiative tranfer into the Ocean-Atmosphere system including polarisation, coupling terms, and the agitated surface.
GNU General Public License v2.0
30 stars 12 forks source link

Comparison of OSOAA and HydroLight #7

Closed ghost closed 3 years ago

ghost commented 4 years ago

I compared the water part of OSOAA with the HydroLight model

and the result is very different from the result of Radiative Transfer Model for the Computation of Radiance and Polarization in an Ocean-Atmosphere System: Polarization Properties of Suspended Matter for Remote Sensing(the relative errors are less than 2%). OSOAA parameter input refers to Influence of the Suspended Particulate Matter on the Satellite Radiance in the Sunglint Observation Geometry in Coastal Waters.

My input file is as follows:

 -OSOAA.ResRoot         :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597
 -OSOAA.Log             :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_osoaa.txt
 -OSOAA.Wa              :      0.66000
 -ANG.Rad.NbGauss       :     -999
 -ANG.Rad.UserAngFile   :     NO_USER_ANGLES
 -ANG.Thetas            :       30.0000
 -ANG.Rad.ResFile       :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/RAD_UsedAngles.txt
 -ANG.Mie.NbGauss       :     -999
 -ANG.Mie.UserAngFile   :     NO_USER_ANGLES
 -ANG.Mie.ResFile       :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/MIE_UsedAngles.txt
 -ANG.Log               :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_ang.txt
 -AER.Waref             :      0.55000
 -AER.AOTref            :       0.10000
 -AER.Tronca            :        1
 -AER.Model             :        2
 -AER.MMD.MRwa          :         -999.000
 -AER.MMD.MIwa          :       -999.00000
 -AER.MMD.MRwaref       :         -999.000
 -AER.MMD.MIwaref       :       -999.00000
 -AER.MMD.SDtype        :     -999
 -AER.MMD.LNDradius     :       -999.00000
 -AER.MMD.LNDvar        :       -999.00000
 -AER.MMD.JD.slope      :       -999.00000
 -AER.MMD.JD.rmin       :       -999.00000
 -AER.MMD.JD.rmax       :         50.00000
 -AER.WMO.Model         :     -999
 -AER.WMO.DL            :       -999.00000
 -AER.WMO.WS            :       -999.00000
 -AER.WMO.OC            :       -999.00000
 -AER.WMO.SO            :       -999.00000
 -AER.SF.Model          :        3
 -AER.SF.RH             :         98.00000
 -AER.BMD.VCdef         :     -999
 -AER.BMD.CoarseVC      :       -999.00000
 -AER.BMD.FineVC        :       -999.00000
 -AER.BMD.RAOT          :       -999.00000
 -AER.BMD.CM.MRwa       :         -999.000
 -AER.BMD.CM.MIwa       :       -999.00000
 -AER.BMD.CM.MRwaref    :         -999.000
 -AER.BMD.CM.MIwaref    :       -999.00000
 -AER.BMD.CM.SDradius   :       -999.00000
 -AER.BMD.CM.SDvar      :       -999.00000
 -AER.BMD.FM.MRwa       :         -999.000
 -AER.BMD.FM.MIwa       :       -999.00000
 -AER.BMD.FM.MRwaref    :         -999.000
 -AER.BMD.FM.MIwaref    :       -999.00000
 -AER.BMD.FM.SDradius   :       -999.00000
 -AER.BMD.FM.SDvar      :       -999.00000
 -AER.ExtData           :     NO_USER_AEROSOLS_PHAZE_FCT
 -AER.MieLog            :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_aermie.txt
 -AER.DirMie            :     /home/qiaofeng/OSOAA/DATABASE/MIE_AER
 -AER.ResFile           :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/PM_AER.txt
 -AER.Log               :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_aer.txt
 -PHYTO.ResFile         :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/PM_PHYTO.txt
 -MLP.ResFile           :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/PM_MLP.txt
 -HYD.Log               :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_hyd.txt
 -HYD.DirMie            :     /home/qiaofeng/OSOAA/DATABASE/MIE_HYD
 -HYD.MieLog            :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_hydmie.txt
 -HYD.Model             :        1
 -PHYTO.JD.slope        :          4.00000
 -PHYTO.JD.rmin         :          0.01000
 -PHYTO.JD.rmax         :        200.00000
 -PHYTO.JD.MRwa         :            1.050
 -PHYTO.JD.MIwa         :         -0.00000
 -PHYTO.JD.rate         :          1.00000
 -PHYTO.LND.SM.SDradius :       -999.00000
 -PHYTO.LND.SM.SDvar    :       -999.00000
 -PHYTO.LND.SM.MRwa     :         -999.000
 -PHYTO.LND.SM.MIwa     :       -999.00000
 -PHYTO.LND.SM.rate     :       -999.00000
 -PHYTO.LND.TM.SDradius :       -999.00000
 -PHYTO.LND.TM.SDvar    :       -999.00000
 -PHYTO.LND.TM.MRwa     :         -999.000
 -PHYTO.LND.TM.MIwa     :       -999.00000
 -PHYTO.LND.TM.rate     :       -999.00000
 -SED.JD.slope          :          4.00000
 -SED.JD.rmin           :          0.01000
 -SED.JD.rmax           :        200.00000
 -SED.JD.MRwa           :            1.150
 -SED.JD.MIwa           :         -0.00100
 -SED.JD.rate           :          1.00000
 -SED.LND.SM.SDradius   :       -999.00000
 -SED.LND.SM.SDvar      :       -999.00000
 -SED.LND.SM.MRwa       :         -999.000
 -SED.LND.SM.MIwa       :       -999.00000
 -SED.LND.SM.rate       :       -999.00000
 -SED.LND.TM.SDradius   :       -999.00000
 -SED.LND.TM.SDvar      :       -999.00000
 -SED.LND.TM.MRwa       :         -999.000
 -SED.LND.TM.MIwa       :       -999.00000
 -SED.LND.TM.rate       :       -999.00000
 -HYD.ExtData           :     NO_USER_HYDROSOLS_PHAZE_FCT
 -AP.MOT                :       -999.00000
 -AP.Pressure           :       1013.00000
 -AP.HR                 :          8.00000
 -AP.HA                 :          2.00000
 -SEA.Depth             :      10000.00000
 -PHYTO.Chl             :          0.00000
 -PHYTO.ProfilType      :        1
 -PHYTO.GP.Chlbg        :       -999.00000
 -PHYTO.GP.Deep         :       -999.00000
 -PHYTO.GP.Width        :       -999.00000
 -PHYTO.Userfile        :     NO_USER_PHYTO_PROFILE                                                                                                                                                                                                                                                                                                                                                                                                                                             
 -SED.Csed              :         15.00000
 -YS.Abs440             :          0.00000
 -YS.Swa                :       -999.00000
 -DET.Abs440            :          0.00000
 -DET.Swa               :       -999.00000
 -PROFILE_ATM.ResFile   :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/PROFILE_ATM.txt
 -PROFILE_SEA.ResFile   :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/PROFILE_SEA.txt
 -PROFILE.Log           :     ,/home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_profile.txt
 -SEA.Ind               :      1.34000
 -SEA.Wind              :      5.00000
 -SEA.Dir               :     /home/qiaofeng/OSOAA/DATABASE/SURF_MATR
 -SEA.Log               :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_sea.txt
 -SEA.SurfAlb           :      0.00000
 -SEA.BotType           :        1
 -SEA.BotAlb            :          0.00000
 -OSOAA.View.Phi        :        180.00000
 -OSOAA.View.Level      :        4
 -OSOAA.View.Z          :       -999.00000
 -OSOAA.ResFile.vsVZA   :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Standard_outputs/LUM_vsVZA.txt
 -OSOAA.View.VZA        :       -999.00000
 -OSOAA.ResFile.vsZ     :     NO_RES_VS_Z_FILE
 -OSOAA.ResFile.Adv.Up  :     NO_ADVANCED_FILE
 -OSOAA.ResFile.Adv.Down:     NO_ADVANCED_FILE
 -SOS.Log               :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/log_sos.txt
 -SOS.IGmax             :      100
 -SOS.ResFile.Bin       :     /home/qiaofeng/OSOAA/results/61b1e956b77a00e3fe3d597340d85597/Advanced_outputs/LUM_SF.bin

My results is:

TSS results:

Figure_1

Chla results:

Figure_2 The x-axis is the wavelength and the y-axis is the $\rho$. $\rho$ Is the output of OSOAA 'I'/cos($\theta$), $theta$ is the sun zenith angle; in HydroLight it is, $\rho$ = pi*Rrs

My question is,

Why is the difference so big, is my input parameter wrong?

What are your input parameters?

ghost commented 4 years ago

The CDOM results seem to have similar trends, but there are still gaps in values

Figure_1

blafranc commented 3 years ago

After checking your OSOAA inputs, it seems that there are some differences in comparison to the parameters used in Chami et al., 2020 (MDPI), which are listed in table 1 of the latter paper.

Here are the differences: -AER.AOTref: Qiaofeng (Q): 0.1, Chami 2020 (C) : 0.2 -AER.SF.Model : Q: 3 ; Ch: 4 -AER.SF.RH : Q: 98 ; Ch: 90 -SED.JD.MIwa: Q: 0.001 ; Ch: 0 -YS.Abs440 : Q: 0 ; Ch: 0.07 -SEA.Depth : Q: 10 000 ; Ch: 50 -SEA.BotType : Q: 1 (seabed albedo=0) ; Ch: 2 (sand)

You could thus see that your OSOAA simulations concern an infinitely deep water (z=10 000m) while Chami et al. (2020)'s simulation concerns a shallow water case ( 50 m) with a sand bottom type. Since we do not know what inputs you used for Hydrolight simulations, it is difficult for us to identify potential differences between OSOAA and Hydrolight ; we suggest you to verify that your Hydrolight inputs parameters are consistent with OSOAA inputs