250 likes | 260 Views
Exploring refinements for global oceanic IOPs products regarding phase function variations and angular effects for improved accuracy and reliability.
E N D
Some refinements for global IOPs products ZhongPing Lee IOPs Workshop, Anchorage, AK, Oct 25, 2010.
Outline: • Two not-so-minor issues: • a) Phase function of molecule vs particle scattering • b) Angular variation 2. QAA vs Optimization
1a). Phase function of molecule vs particle scattering Why a 0.0012 m-1 bias for bbp? (Behrenfeld et al 2005)
Rrs model commonly used: g or {g0 & g1} depends on phase function shape G(0+) (sr-1) bb/(a+bb)
G or g is NOT a monotonic function of bb/(a+bb)! Facts: • G or g depends on phase function shape. • Molecular & particle scattering have significantly different phase function shapes. • We can have the same bb/(a+bb) for both molecular and particle scattering phase functions, but different G or g values. • For oceanic waters, blue band is dominated by molecular scattering; green/red bands by particle scattering. Conclusions: Need to account for the phase function shape (molecular vs particle) effects, especially for oceanic waters.
Explicitly separate the effects of molecule and particle phase function effects Lee et al (2004) G from formula (sr-1) G from Hydrolight (sr-1) Caveats: Complex function; cannot invert a&bb algebraically.
A practical model for algebraic/analytical inversion: gw&gp model g model (conventional, widely used)
Effects of separating molecular/particle phase function: MODIS Aqua
Comparison of results bbp555 a443 bbp is reduced by 40%; even more for adg443 Consistent in general with the results (Werdell & Franz) that compare f/Q LUT vs Gordon 88 formula. aph443 adg443
Conclusion: Significantly different IOPs (in clear oceans) retrieved when the phase function shape effects are considered.
1b). Angular variation ψ θv θS Ω(10, 20, 30) measured photons going further away from Sun (~forward scatter) Ω(10, 20, 150) measured photons going closer to Sun (~backscatter)
400 nm 640 nm Water-leaving radiance in the Sun plane, zenith dependence (arrow length indicates radiance value) Bottom line: Water-leaving radiance, or reflectance, is a function of angles.
Current, “standard”, approach to deal with angular variation: empirical Rrs(Ω) Chla Rrs[0] IOPs f/Q based on Case-1 The two semi-empirical steps could be omitted with an Rrs model accounts for the angular effects.
Candidate models for angular Rrs: Albert and Mobley (2003) Lee et al (2004) Park and Ruddick (2005) Van Der Woerd and Pasterkamp (2008) Caveats: 1. Some are not visible/transparent about the physics 2. Some are not easily invertible algebraically
A practical choice for algebraic/analytical inversion: Global distribution of Rrs(443) 1:1 G from model [sr-1] (Ω: 60,40,90) Rrs443 [sr-1] G ~ 0.07 G from HL simulation [sr-1]
Angular-dependent model coefficients for Rrs(Ω): Table ((7x13+1)x4x6) array, 2208 elements) of {G(Ω)} (if based on Chl, it is 6x13x7 = 546 elements per band per Chl)
IOP retrieval from angular Rrs: G[Ω] Rrs(Ω) IOPs empirical QAA, optimization, linear matrix, etc. Rrs(Ω) Chla Rrs[0] IOPs Now we get IOPs straightforwardly (one step) from Rrs(Ω)! f/Q based on Case-1
2. QAA vs spectral optimization All start with: “essential” difference: 1) “Philosophic” difference 2) Difference in measuring “signal” vs “noise” 3) Impact of data and/or model 4) Processing efficiency 5) Stability
1) “Philosophic” difference QAA: Total first, then individuals Optimization: Individuals first, then total or simultaneously Under QAA, the error propagation is visual and easy to quantify Δbbp, Δa will be 0 if ∆a(λ0) and ∆η are 0.
In addition: (Lee and Carder 2004) Spectral shape of aph(λ) is a property we want to obtain from each measured Rrs(λ). Spectral optimization assumes a spectral shape before its derivation. aph(λ) [m-1] Wavelength [nm]
2) Difference in measuring “signal” vs “noise” Optimization algorithm (e.g., GSM01, HOPE) Algebraic algorithm (QAA, LMI) (Lee et al. 2002, Hoge and Lyon 1996) (Roesler and Perry 1996, Lee et al. 1996, Maritorena et al. 2001, Doerffer 1999) Every measurement is perceived as signal. Mis-match is perceived as measurement noise.
3) Impact of data and/or model The increase trend of bbp from Aqua, by GSM01, is considered due to “error” of Aqua Rrs412. QAA will have different patterns, at least for bbp, as it does not use Rrs412 for bbp derivation. (Maritorena et al 2010)
4) Processing efficiency (Lee et al 2002) ‘Resolved’ for multi-band data; hyperspectral data? 5) Stability Optimization: software impact; optimized or closed to be optimized?
The pathway for global IOPs products by GIOP: As QAA and optimization schemes (and other semi-analytical algorithms) initiate from the same physics, then IOPs from both retrievals can serve as a consistence/reliability check: Agreeable results highly reliable! Different results need further diagnose …