Download - University of Washington, February 2016
![Page 1: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/1.jpg)
Machine Learning on Images: Combining Passive Microwave and Optical Data to Estimate Snow Water Equivalent in Afghanistan’s Hindu Kush
Jeff DozierUniversity of California, Santa Barbara
University of Washington, 2016-02-25
![Page 2: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/2.jpg)
![Page 3: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/3.jpg)
![Page 4: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/4.jpg)
News item from IRIN (UN)http://www.irinnews.org/Report/93781/Analysis-Afghan-drought-conditions-could-spell-disaster
KABUL, 21 September 2011 (IRIN) – “The current dry spell sweeping across Afghanistan’s northern, northeastern and western provinces could lead to a large-scale food crisis and the humanitarian community should act quickly to ensure this does not degenerate into a disaster, government and aid officials warn.”
![Page 5: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/5.jpg)
![Page 6: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/6.jpg)
AMSR2,1.49 km3
Reconstruction,3.54 km3
SNODAS,5.22 km3
10 100
200
300
400
500
600
700
800
900
103Snow water equivalent, Sierra Nevada, mm, 2014-04-01
![Page 7: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/7.jpg)
Question: can we estimate the real snow water equivalent, in places like the Hindu Kush, from information available during the season?
2011-04-01, passive microwave SWE, 27 km3
SWE, mm snow fraction
2011-04-01, MODIS snow cover, 190,000 km22011-04-01, reconstructed SWE, 68 km3
![Page 8: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/8.jpg)
Fractional snow cover from MODIS2015-04-01
𝑅𝜆=𝜖𝜆+∑𝑘=1
𝑁
𝑓 𝑘𝑅 𝜆 ,𝑘 Data from http://snow.jpl.nasa.gov
![Page 9: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/9.jpg)
Comparison of MODIS (500m) and Landsat (30m) snow fraction, in the Sierra Nevada
32 scenes with coincident MODIS and Landsat imagesAverage RMSE = 7.8%Range from 2% to 12%
![Page 10: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/10.jpg)
Energy balance reconstruction
American River basin Snow Pillows
on day𝑛 ,𝑆𝑊𝐸𝑛=𝑆𝑊𝐸0−∑𝑗=1
𝑛
𝑀 𝑗
when𝑆𝑊𝐸𝑛=0 ,𝑆𝑊𝐸0=∑𝑗=1
𝑛
𝑀 𝑗
𝑀 𝑗=𝑀𝑝 𝑗× 𝑓 𝑆𝐶𝐴 𝑗
From energy balance model, driven by NLDAS or GLDAS
𝑀𝑝 𝑗
![Page 11: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/11.jpg)
11
Downscaling
Sola
r rad
iatio
nLo
ngwa
ve ra
diat
ion
![Page 12: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/12.jpg)
Water balance? Over a year so is ? (Only for Reconstruction)
![Page 13: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/13.jpg)
Comparison of reconstructed SWE with Airborne Snow Observatory
![Page 14: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/14.jpg)
Comparison of reconstructed SWE with Airborne Snow Observatory
![Page 15: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/15.jpg)
Sierra Nevada, AMSR-E & Reconstruction
![Page 16: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/16.jpg)
![Page 17: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/17.jpg)
Persistent sources of uncertainty in passive microwave retrieval of SWE:deep snow, vegetation, sub-grid heterogeneity
[Luojus et al., GlobSnow ATBD, http://www.globsnow.info/docs/GS2_SWE_ATBD.pdf]
![Page 18: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/18.jpg)
Question: can we estimate the real snow water equivalent from information available during the season?
2011-04-01, passive microwave SWE, 27 km3
SWE, mm snow fraction
2011-04-01, MODIS snow cover, 190,000 km22011-04-01, reconstructed SWE, 68 km3
![Page 19: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/19.jpg)
The problem• Discover a pattern that uses real-time data—passive
microwave SWE and snow-covered area from optical sensor—to match reconstructed SWE, which is available only after the snow is gone
Approach• Among the choices for machine learning methods, start
with neural networks• Investigate alternatives later
• Develop a training set with 2003-2004 data, and predict 2005, from March 1 to June 1• Then adapt 2005 results, and predict 2006• And so on
![Page 20: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/20.jpg)
Results
(modest improvement in R2 over time, but unpredictable bias)
![Page 21: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/21.jpg)
2009 integrated total SWE (km3)
![Page 22: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/22.jpg)
Insight from comparing statistical distributions (red predicted, blue reconstructed)
![Page 23: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/23.jpg)
Improvements to the reconstruction
![Page 24: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/24.jpg)
Remotely sensed albedo of fractional snow
![Page 25: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/25.jpg)
Dirty snow albedo has a similar spectral shape to fine-grain clean snow
[Warren, Rev Geophys, 1982]
![Page 26: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/26.jpg)
Errors in assimilated energy inputs
![Page 27: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/27.jpg)
“Clouds” in cloud-free scenes (Landsat OLI, Sierra Nevada)a
2013-05-20 2015-03-07
![Page 28: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/28.jpg)
Steps forward (but already better than what we have now)• Improve reconstruction
• Sensitive to albedo, so fix the albedo estimate of fractionally snow-covered pixels
• Consider alternatives to NLDAS/GLDAS• Solve the snow-cloud discrimination
• Incorporate better microwave snow retrieval methods• Improve machine learning
• Examine other ensemble statistical methods, e.g. genetic programming, regression boosted decision trees, support vector machines
• Find another input variable that helps with the bias, e.g. satellite precipitation products, models like WRF
• But . . . We’re encouraged by the results, as a way to estimate snow water equivalent in mountains without existing sensors
![Page 29: University of Washington, February 2016](https://reader036.vdocument.in/reader036/viewer/2022081521/58745fc21a28abab198b4dbf/html5/thumbnails/29.jpg)
“the author of all books” (James Joyce, Finnegan’s Wake)
Slides available at http://www.slideshare.net/JeffDozier/
Finis
Data available at ftp://ftp.snow.ucsb.edu/pub/org/snow/users/dozier/MachineLearning/