The following sets of images and plots contain 2 sets of comparisons
between Corps of Engineers data and SSM/I-derived SWE. For each date,
the upper left subimage contains Army SWE data
(converted from inches to mm) for the Oahe basin. Each subsequent
subimage is one of the passive microwave algorithm outputs for the
same date, masked for the basin area only. SSM/I algorithms include:
NSIDC1: Chang SMMR algorithm (combination of horizontally
polarized channels) modified for use with SSM/I channels
NSIDC3: Chang's AMSR algorithm, a decision tree built on top of
the basic horizontal channel difference algorithm
GW: Goodison & Walker algorithm (combination of verticlaly polarized channels)
Nagler: Thomas Nagler algorithm, H-pol channels, uses 85 GHz
GB: Grody & Basist, snow extent only
The red outline in the passive microwave images is the outline of the
Army data SWE extent. The plots contain scatter plots of the Army
data vs. each algorithm output. Only SWE values are considered
(i.e. "wet" snow flags and pixels with missing TBs are eliminated).
The comparison data for Survey x SCA for SSM/I data from March 1
was dated March 4_3, indicating that survey data for March 1 were
crossed with SCA data from the closest dates possible, March 4-6
(ref. e-mail from Emily to Mary Jo).
Scatter plots were created by regridding the Army data to a 25-km
EASE-Grid, using drop-in-the-bucket average interpolation. Resulting
Army SWE pixels were compared with non-missing SSM/I-derived SWE for
that pixel.
Nagler algorithm outputs depth. I converted this to SWE (mm)
using a factor of 3, so that images are compared quickly visually.
Statistical correlation values for SWE improved from ~0.25 for
the comparison to the survey data to ~0.75 for the comparison to the
survey data multiplied by the AVHRR-derived snow-covered area
fractions. Note that the prominent vertical artifact at about 50 mm
Army SWE in the Survey comparisons is removed in the Survey X SCA
comparisons.
For the two dates compared (March 1 and 15), the SWE algorithms
performed consistently, relative to one another. In both cases, the
decreasing order of correlation (best to worst) was
GW (Goodison & Walker)
Nagler
NSIDC1
NSIDC3
Due to the differences and theoretical improvements in NSIDC3
versus NSIDC1, Richard and I expected NSIDC3 to improve the
correlation when compared with NSIDC1, but this was not the case. I
don't yet know why. Then again, difference in r values of 0.755
vs. 0.743 probably isn't significant?
GW was better than all of the rest, not particularly surprising,
since this area is very similar to the Canadian prairies, so that
algorithm really ought to do well, here.
Hard to say much about GB (Grody-Basist) since it's only a snow
cover algorithm. It puts down too much snow when compared to Survey
SWExSCA.