1 / 39

Background [1] Sources of Precipitation Information

A Prototype Algorithm for Gauge – Satellite Merged Analysis of Daily Precipitation over Land 2009.06.23. Background [1] Sources of Precipitation Information. Gauge measurements Accurate point measurements Insufficient station networks Radar observations

bela
Download Presentation

Background [1] Sources of Precipitation Information

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A Prototype Algorithm for Gauge – Satellite Merged Analysis of Daily Precipitation over Land2009.06.23.

  2. Background [1]Sources of Precipitation Information Gauge measurements Accurate point measurements Insufficient station networks Radar observations reasonable temporal / spatial coverage around radar sites Relatively accurate estimations Blocked by mountains / inavailability over some remote regions Satellite estimates Infrared (IR) observations Indirect and less accurate estimates Frequent sampling (every 30-min) Fine spatial resolution (4kmx4km) Passive Microwave (PMW) observations Physically-based / more accurate estimates of instantaneous rain rates Poor sampling (up to twice a day from a single satellite) Precipitation generated by numerical models Less-than desirable performance for convective precipitation and precipitation ocean Good performance for large-scale stratiform precipitation over extra-tropics during cold season

  3. Background [2]Strategy and Status Basic Strategy Two-step merging approach to combined gauge & sate First, to merge observations from multiple satellites to form an all-satellite merged precipitation estimates (CMORPH) Over land, to combine the all-satellite estimates with in situ observations (gauge); Over ocean, to adjust the CMORPH against a long-term data with coarser resolution (pentad GPCP) to facilitate the definition of anomaly Status Land : finished development of the prototype algorithm (to introduced here) in process of implementation over the globe Ocean : finished development of a test algorithm

  4. Objective To develop a prototype algorithm to define high-resolution analysis of daily precipitation over land by combining gauge observations and satellite estimates 0.25olat/lon over China Tested for May – September, 2007

  5. Input Data Gauge data Gridded analysis of daily precipitation defined by interpolating gauge observations from over 2,400 stations over China Daily reports are available from only 198 Chinese stations through GTS High-resolution satellite precipitation estimates Generated by CPC Morphing Technique (CMORPH, Joyce et al. 2004) Precip estimstes of 8km/30-min resolution regridded into 0.25olat/lon/daily

  6. Gauge Data • Daily precipitation • 0.25olat/lon over China • May 1 – Sep.30, 2007 • Daily precip analysis defined by interpolating station data over ~2400 gauges and used as the truth in the development

  7. Characteristics of Input Data Sources Gauge data Accurate point measurement of precipitation (unbiased) Gauge network of insufficient density over most of the global land areas Random error in gauge-based analysis representing areal mean due to spatial variations in precip and inappropriate gauge network to catch them Satellite estimates Ability to catch spatial pattern and temporal variations of precipitation with reasonable accuracy Regionally, seasonally, and range-dependent BIAS Random error

  8. Basic Strategyto Merge Gauge and Satellite Data First Step --- Removing Satellite Bias Assuming gauge-analysis is unbiased Removing bias in the satellite estimates through comparison against the gauge analysis Second Step --- Combining Gauge and Bias-corrected Satellite Estimates Combination through OI Bias-corrected satellite estimates as first guess Gauge analysis used to improve the first guess over regions with gauge stations

  9. Part I : Bias Correction

  10. Bias Correction Through PDF Matching Principal Match the PDF of the CMORPH against that of daily gauge to define and remove the bias, assuming PDF of the gauge analysis represents that of the truth Implementation Collect co-located pairs of gauge and CMORPH over grid boxes within a spatial window centering at the target grid box and a time period ending at the target date; Define PDF for the CMORPH and gauge analysis

  11. Sample ResultsPrecipitation for August 2, 2007

  12. Cross-Validation TestsHow are they carried out Each time, gauge analysis at 10% randomly selected grid boxes is withdrawn PDF bias correction is performed using gauge data over the remaining 90% grid boxes The above processes are repeated 10 times so that analysis at each grid box is withdrawn once Compare the original and bias-corrected CMORPH with the gauge analysis at the withdrawn grid boxes to ensure independence Only data over grid boxes with at least one gauge are included in the calculation to ensure reliable quality of the ground truth

  13. Cross-Validation TestsCombined space / time domain • Successful correction of the bias; • Substantial improvements in correlation

  14. Cross-Validation TestsTime Series

  15. Cross-Validation TestsPDF of Daily Precip over 0.25olat/lon Grid Boxes • Adjusted CMORPH presents PDF very close to that of the gauge analysis • Frequencies of no-rain and light rain in the adjusted CMORPH is slightly larger / smaller than those in the gauge analysis

  16. Part II : Combining Gauge with Bias Corrected CMORPH

  17. The OI-based Combining Algorithm Analyzed fields of precipitation are defined by combining bias-corrected satellite estimates and the gauge data (gauge-based analysis, in this test) The satellite estimates are used as the first guess The gauge data are used as the observations to correct the first guess where gauge data are available The combined analysis will be the same as the first guess over regions with no gauges nearby Weighting coefficients used in the combining are computed from error structures for the first guess and observation fields following the OI formula

  18. Error Structures See backup slides for details Both gauge data and bias-corrected satellite estimates are bias free Random error in gauge analysis is a function of gauge network density and intensity of precipitation No spatial correlation in gauge random error CMORPH random error is a function of sampling size (number of original satellite data used) and the intensity of precipitation determined empirically through comparison with gauge data Spatial correlation in CMORPH error declines with distance

  19. OI Combined Analysis [1] • Sample for August 2, 2007 • Both the precipitation analysis and the error estimation seem reasonable

  20. OI Combined Analysis [2] • Precipitation distribution looks better than the gauge-based analysis and the CMORPH estimates

  21. OI Combined Analysis [3] • Comparison between the combined analysis and the CMORPH, for August 2, 2007 • Differences over regions with gauge data • Grid boxes with negative/positive differences distributed randomly, implying the bias correction for CMORPH works well

  22. OI Combined Analysis [4] • Comparison between the combined analysis and the gauge data, for August 2, 2007 • Differences over regions with less gauge data • Grid boxes with negative/positive differences distributed randomly

  23. OI Combined Analysis [5] • Mean precipitation (mm/day) for May 1 – Sep.30, 2007, from gauge analysis (top), bias corrected CMORPH (middle), and combined analysis (bottom) • Combined analysis presents better distribution than gauge and CMORPH

  24. OI Combined Analysis [6] • Comparison between the combined analysis - the gauge data, and combined analysis – CMORPH for May 1 – Sep.30, 2007 • RMS differences are computed for grid boxes with different numbers of gauges • Differences between the combined analysis and the gauge / CMORPH are smaller / larger over grid boxes with more gauges • This confirms that our algorithm works as we designed

  25. Part III : Real-Time Operation System

  26. Objective To construct a LINUX-based real-time operational system to create analyses of daily precipitation over East Asia by combining gauge observations and satellite estimates 60oE – 160oE; 0o – 60oN 0.25olat/lon grid Daily precipitation ending at 00Z

  27. Gauge Data CMORPH 8km/30min Interpolation Regridding Gauge Analysis (0.25olat/lon) CMORPH (0.25olat/lon) Merging Binary Data GIF Pictures Structure of the Analysis System

  28. Gauge Analysis (0.25olat/lon) CMORPH (0.25olat/lon) Merging Procedures PDF Matching Bias Correction Bias-Corrected CMORPH OI Combining Merged Analysis

  29. Sample Merged Analysis Over China, June 4, 2008

  30. Sample Merged Analysis Over East Asia, June 4, 2008

  31. Sample Merged Analysis Over China, Animation

  32. Sample Merged Analysis Over East Asia, Animation

  33. A prototype algorithm is developed to define daily precipitation analysis on a 0.25olat/lon over land by merging gauge observations and CMORPH satellite estimates • More work needs to be done to refine the system and to implement the system for global applications • Fine tuning of data pair collecting scheme for the PDF matching procedures • Accurate definition of error structure for the gauge analysis and satellite estimates • Comprehensive examinations of system performance • Modifications of the system structure to better fit the input data availability, quality, and stability Summary

  34. CMORPH Error [1] • Assuming CMORPH error variance is proportional to the estimated precipitation amount • Proportional constant is determined through comparison with gauge analysis over grid boxes with at least one gauge using data over China for summer 2007 • CMORPH error is defined for each grid box of 0.25olat/lon using the curve in the bottom panel of the figure • Further work is required to set the error also as a function of region, season and sampling

  35. CMORPH Error [2] • Sample CMORPH error estimation for August 2, 2007

  36. CMORPH Error [3] • In performing OI, we also need to set the correlation between errors occurred over two separate points • For observations (gauge), we assume that the correlation between errors at two different locations is zero • For CMORPH, we assume the error correlation is a function of the distance between the two points, expressed as an exponential function of the negative distance • In this test, we arbitrarily set the e- folding distance as 30km • Further work is required to accurately define the error function

  37. Gauge Error [1] • Assuming gauge error variance is proportional to the observed precipitation amount and inversely proportional to the number of gauges • Proportional constant is set arbitrarily in this test, to make sure the gauge error is smaller than the CMORPH error when there is at least one gauge • Further work is needed to accurately define the error function, through comparison with gauge data over a region with dense network (e.g. OKALAHOMA MESONET)

  38. Gauge Error [2] • Sample gauge error estimation for August 2, 2007

  39. Cross-Validation TestsSpatial Distribution

More Related