1 / 51

Overview

Object-based precipitation analysis: application to tropical cyclones and the Slovenian radar data Mini Workshop on NWP Modelling Research in Slovenia 15.December 2011. Gregor Skok Julio Bacmeister , Joe Tribbia , Benedikt Strajnar , Jože Rakovec , Anton Zgonc , Mark Žagar. Overview.

calida
Download Presentation

Overview

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. Object-based precipitation analysis: application to tropical cyclones and the Slovenian radar dataMini Workshop on NWP Modelling Research in Slovenia15.December 2011 GregorSkokJulio Bacmeister, Joe Tribbia, BenediktStrajnar, JožeRakovec, Anton Zgonc, Mark Žagar

  2. Overview • Object based analysis • Analysis of tropical cyclone precipitation using satellite data • Hail-area tracking algorithm using radar data

  3. Object based analysis • Doctoral thesis:“Object-Based Analysis And Verification Of Precipitation Over Low- And Mid-Latitudes”

  4. Motivation • TRMM 3B42 • 3 hourly precipitation accumulations, 0.25 deg

  5. Object identification method • Based on MODE - Method for Object-based Diagnostic Evaluation developed by Devis et al (2006a,b) • Part of Model Evaluation Tools (MET) verification package developed by NCAR • 3 steps: • Smoothing • Thresholding • Identification of self-enclosed areas as objects • The method tries to simulate what a human forecaster or analyst might infer by a more subjective visual evaluation of a field => (Objective) simulation of a subjective evaluation

  6. Thresholding only Smoothing Thresholding after smoothing Original MODE method

  7. Methodology - Extended method • Time evolution of objects • “tri-dimensional” objects • Enables study of properties: location, size, shape movement, lifespan, total object precipiation, ….

  8. Doctoral thesis:Pacific, 6-years of 3-hourly TRMM 3B42 precipitation data

  9. highest density of objects with a longer lifespan (red) is in the ICTZ and in the low-latitudes in the west • eastern tip of the ITCZ – mainly objects with short lifespan • Central America – mostly shortlived objects • ……. Trajectories for 2001 BLUE – short lifespan, RED – long lifespan

  10. Movement in the northern and southern parts of domain is predominantly eastward • In the ITCZ region, movement in both directions is present although westward movement (green) is more frequent • In the eastern and western part of the ITCZ the westward movement is clearly dominant. Trajectories for 2001 ORANGE – eastward, GREEN – westward

  11. Number of objects vs. lifespan Straight in a Log-Log graph = Power law

  12. Analysis of tropical cyclone precipitation using satellite dataGregor Skok, Julio Bacmeister, Joe Tribbia • TRMM 3B42 precipitation data • The IBTrACS tropical cyclones track database • 11 years - 1998-2008

  13. FiT -“Forward in Time” object identification

  14. FiT -“Forward in Time” object identification

  15. FiT -“Forward in Time” object identification

  16. FiT -“Forward in Time” object identification • Merging! • Forced to check into the past and also perform merger there

  17. FiT -“Forward in Time” object identification

  18. FiT -“Forward in Time” object identification

  19. FiT -“Forward in Time” object identification • Don’t allow merging! • Larger “wins” • No need to check into the past -> only forward in time • Side benefit: faster and less memory consuming

  20. FiT -“Forward in Time” object identification

  21. The problem of “missed” precipitation • Inside objects (threshold 7mm/3h) there is only 50 % of all precipitation. • The other 50 % is located in a dislocated self-enclosed areas of low-intensity precipitation or just outside the borders of objects. • We want to include nearby low-intensity precipitation for TC analysis

  22. Estimation of object precipitation by “grown” objects Precipitation threshold

  23. Estimation of object precipitation by “grown” objects Precipitation threshold Secondary threshold

  24. Estimation of object precipitation by “grown” objects Sequentially grow objects: 1 iteration

  25. Estimation of object precipitation by “grown” objects Sequentially grow objects: 4 iterations

  26. Estimation of object precipitation by “grown” objects Sequentially grow objects: 9 iterations -> end

  27. Estimation of object precipitation by “grown” objects Might be more unattributed low intensity precipitation below secondary threshold Unattributed precipitation In GROWN objects (to 1 mm/3h) now 75 % of all precipitation

  28. IBTrACS database

  29. Identification of TC objects Object MATCH?YES IBTrACS TC center Distance smaller than 2.5 deg MATCH?YES Distance larger than 2.5 deg MATCH?NO

  30. Animation

  31. TC precipitation [mm/day]

  32. Contribution of TC precipitation to all precipitation [%]

  33. Zonal means of TC precipitation GLOBAL SEA LAND

  34. Regions

  35. Regions • TCs contribute about 4 % (on average 40 km3/day) • This percentage is on average higher for oceans than for land (4.8 % vs. 1.4 %). • NH the TCs contribute around 5.1 % and in SH about 2.8 % precipitation • Compared to the oceans, the land sub-regions have much smaller TC precipitation volumes. • some land regions get over 3 %:Australia, Maritime continent with E Asian islands and E Asia • some seasons TCs contribute more precipitation; i.e. N America (6 % in SON), Australia (4 and 5.5 % in DJF and MAM), Maritime continent with E Asian islands (5,5 % in JJA and SON), E Asia (3 and 6 % in JJA and SON) and S Asia (4 % in SON)

  36. Yearly global TC precipitation

  37. Hail-area tracking algorithm using radar dataGregor Skok, Benedikt Strajnar, Jože Rakovec, Anton Zgonc, Mark Žagar • Using volumetric radar data from Lisca – 8 years 2002-2010 • Areas with hail precipitation identified using a combination of two methods: Waldwogel et al. (1979) and Gmoser et al.(2006). • This produces a 2D binary field – hail yes/no. • Radar scan is performed every 10 minutes. A sequence of 2D binary “hail” fields is fed into the object identification algorithm • The movement of objects represent the movement of areas with hail precipitation

  38. Animation

  39. Animation

  40. Hail area tracking • No smoothing/thresholding possible since the field is binary • The hail areas are relatively small and move fast – they often do not overlap in 10 minute intervals • To overcome this problem the objects are artificially grown in all directions • This improves the overlap but can merge nearby objects • The value of parameter describing the “extent” of growth has to be selected carefully – sensitivity analysis

  41. not grown

  42. grown by 1 km

  43. grown by 2 km

  44. grown by 3 km

  45. Number of objects Not an exponential distribution Hail events not a random shortlived process

  46. Results Direction of movement by azimuth (regardless of lifespan) Trajectories longer than 150 min. Red = eastward, blue = westward

  47. Thank you

  48. ORIGINAL

  49. SMOOTHING

  50. OBJECTS AFTER THRESHOLDING

More Related