230 likes | 384 Views
Geostatistical Mapping of Mountain Precipitation Incorporating Auto-searched Effects of Terrain and Climatic Characteristics. Huade Guan, John L. Wilson , Oleg Makhnin New Mexico Institute of Mining and Technology American Meteorological Society 85 th Annual Meeting
E N D
Geostatistical Mapping of Mountain Precipitation Incorporating Auto-searched Effects of Terrain and Climatic Characteristics Huade Guan, John L. Wilson, Oleg Makhnin New Mexico Institute of Mining and Technology American Meteorological Society 85th Annual Meeting San Diego, Jan.11, 2005
Why use gauge data for precipitation mapping in mountains? • Problems with NEXRAD • Beam blockage • Snow estimation • 4km pixel size NEXRAD rainfallNew Mexico, July 1999 From Hongjie Xie, 2004
today Four types of mapping approaches (examples) Cokriging (P-Z)& De-trended residual kriging
Orographic lifting, & hindrance Reduction in virga effect wind direction:ω terrain aspect: α T↑ T↓ Elevation (Z) wind P (windward)>P( leeward) P (low Z)<P( high Z) P (low Z)<P( high Z) terrain aspect Physical process (1)Orographic effects on precip. We use cos (α-ω) to approximate terrain aspect effects
Study area Physical process (2)Atmospheric effects on precipitation How does this heterogeneous atmospheric moisture distribution (or gradient in atmospheric moisture) influence precipitation? We use geographic coordinates (Longitude or X, and Latitude or Y) to capture the effect of gradient in atmospheric moisture on precipitation GOES East 4-km, infrared imagery 2001.05.04
moist. flux dir. aspect Auto-search orographic and atmospheric effects Regression: gradient in moist., elevation, aspect & moist. flux direct. Data: Gauge precip: X, Y, P; Elev. DEM: X, Y, Z, a ; But what about moisture flux direction, w ?
Auto-search orographic and atmospheric effects Regression turns to: where b5=b4cosω, and b6=b4 sinω, implicitly contain the moisture flux direction. And b1 and b2include the information of gradient in atmospheric moisture. • For example, if b5 >0 and b6 >0, ω= atan (b6/b5) • Similarly, if b1 >0 and b2 >0, gradient in atmospheric moisture, or the wetter direction = atan (b1/b2)
ASOADeKAuto-Searched Orographic and Atmospheric effects De-trended Kriging • Auto-determine moisture gradient, elevation, & moisture direction effects via regressions • Included in b0, b1, b2, b3, b5, and b6. • Construct regression map from DEM • Find residual at each gauge • Generate residual (or de-trended) map by kriging • Construct the final precipitation map • Regression map + residual map
ASOADeK regression improves estimates aspect + moisture flux direction moisture gradient
ASOADeK inferred moisture flux directions November January April July Winter: Southwesterly Summer: Southeasterly
Two weather • patterns • in Summer • Southwesterly moisture flux related North American Monsoon • (picture to the right) • Easterly moisture flux • ASOADeK: Southeasterly From NOAA Mixture of the two may give apparent southeasterly moisture flux as inferred from ASOADeK
From Sellers and Hill, 1974 Weather pattern related to heaviest winter precipitation Southwesterly moisture flux at the study area, ASOADeK: Southwesterly
ASOADeK regression vs. PRISM • Model estimates from both models v. measured values • Scatter plots and fits (R2) • For three months: Feb, May & Aug. • ASOADek regression only! • No residual kriging
ASOADeK vs. PRISM For ASOADeK let’s now include the residual map • Precipitation maps for both models, and • QQ plots, • for same three months: Feb, May & Aug. • ASOADek regression plus residual map.
Cross validation results:ASOADeK gives better estimatesthankriging &cokriging
Conclusions • ASOADeK detects regional climate settings using only precip. gauge data in mountainous terrain. • ASOADeK vs. PRISM • Precipitation maps: ASOADeK ≈ PRISM • ASOADeK product has higher spatial resolution • ASOADeK vs. other geostat. approaches • Precipitation estimates improved in comparison with krigng and co-kriging.
Future work • Further testing ASOADeK auto-searching capacity • Event cases • Other geographic regions • Applications & Extensions of ASOADeK • Mapping Precipitation in mountainous regions • Studying ENSO/PDO effects on precipitation distribution • Recovering NEXRAD beam-blockage shadow • Downscaling precip. products, e.g., NEXRAD
ain Thank you!