1 / 76

Dr Masami Tokuno

Assimilation of QuikSCAT/SeaWinds Ocean Surface Wind Data into the Global Data Assimilation System at JMA. From NASA Home Page. Dr Masami Tokuno. Mr. Yasuaki Ohhashi. Mr. Takeshi Iriguchi. NWP Division , JMA. CONTENTS. Introduction 2. Assimilation system for scatterometer

ajuarez
Télécharger la présentation

Dr Masami Tokuno

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. Assimilation of QuikSCAT/SeaWinds Ocean Surface Wind Data into the Global Data Assimilation System at JMA From NASA Home Page Dr Masami Tokuno Mr. Yasuaki Ohhashi Mr. Takeshi Iriguchi NWP Division , JMA

  2. CONTENTS • Introduction • 2. Assimilation system for scatterometer • 3. Impact study for QuikSCAT wind data • 3.1 Configuration of the experiment • 3.2 Quality of SeaWinds data • 3.3 An example of OSE (1 Dec. 2001) • 3.4 Impacts to GSM( Dec. 2001, July 2002) • 3.5 Impacts to typhoon truck forecasts in GSM • 4. Summary

  3. Numerical Forecast and Satellite Data 衛星観測データ Observation Initial Analysis Numerical Forecast Calculation Grasp the actual state of wind, temperature, humidity Numerical Prediction Model Translate into Weather Forecast

  4. Plan to use scatterometer data for NWP FY2001 2002 2003 2004 2005 GA3DV MA4DV RA4DV GA4DV ADEOS-II/SeaWinds assimilation QuikSCAT Operational (Global) 4D-Var QuikSCAT Apply Regional Model

  5. Assimilation of observation data for GSM 3D-Var 00UTC 3D-Var 18UTC 3D-Var 06UTC Ob.Ob. Ob. 3D-Var Early Analysis Forecast

  6. Coverage of QuikSCAT Data used in the Cycle Analysis 2002/07/08/00UTC 2002/07/08/06UTC 2002/07/08/12UTC 2002/07/08/18UTC

  7. Coverage of QuikSCAT Data used in the Early Analysis 2003/01/29/00UTC 2003/01/29/12UTC

  8. The Number of QuikSCAT Data used in the Cycle and Early Analysis, respectively Cycle Analysis Early Analysis days 2003年

  9. SCAT WIND EXTRACTION SeaWinds ambiguity winds (NESDIS L2B) [QC] Rain check,etc (Wind data in rain are not used) Ambiguity removal (NWP nudging, Median filter) [QC] Speed Check:First Guess [QC] Direction check:Group QC, Save calm wind A scatterometer wind field Thinning out windon 1x1deg grid 3D-VAR Analysis

  10. Experiment of Assimilation and Forecast (OSE) Case 1 Case 2 Start up : 1 ~ 7 Dec. 20011~7 Jul. 2002 Period: 8 ~22 Dec. 2001 8~22Jul. 2002 9 days forecast (every 12UTC) Global Data Assimilation System • Three-dimensional Assimilation System (3D-VAR) Forecast Model T213L40 Global Spectrum Model    (JMA’s Operational Model)

  11. OSE (Observation System Experiment) Control run [Cntl] (December) Conventional data + AMW + GMS-5 Moisture Bogus + NOAA-14 + NOAA-15 (SATEM/TOVS) Control run [Cntl] (July) Conventional data + AMW + NOAA-14 + NOAA-15 (SATEM/TOVS) Test run : [Only Wind] Control + QuikSCAT/SeaWinds (only surface wind) (no ERS-2/AMI data in Cntl & Test)

  12. SCAT WINDS(Speeds) Compared to Fist Guess(Dec., ’01, July ’02) Speed 8-22 Dec. Speed 8-22 July N=537134 N=528941 SCAT WIND RMSE RMSE 1.91 m/s 1.87m/s BIAS BIAS 0.81 m/s 0.66 m /s First Guess (m/s) First Guess

  13. SCAT WINDS(Speeds) Compared to Fist Guess( 8 – 22 Dec., ’01) Number SCAT Winds First Guess First Guess Histogram SCAT Winds Speed(m/s) Reject Number SCAT Winds All Observation RMSE ー First Guess 1.91 m/s BIAS Histogram 0.81 m/s Reject (m/s)

  14. SCAT WINDS(Speeds) Compared to Fist Guess( 8 – 22 July, ’02) Number SCAT Winds First Guess First Guess Histogram SCAT Winds Speed(m/s) Reject 個数 SCAT Winds All Observation RMSE ー First Guess 1.87 m/s BIAS Histogram 0.66 m/s Reject (m/s)

  15. SCAT WINDS(Dir.) Compared to Fist Guess(Dec., ’01, July ’02) Dir. 8-22 Dec. Dir. 8-22 July N=537134 N=528941 SCAT Wind RMSE RMSE 20.18 deg 19.70 deg. BIAS BIAS 0.81 deg. 0.06 deg. First Guess First Guess

  16. SCAT WINDS(Dir.) Compared to Fist Guess( 8 – 22 Dec., ’01) Number SCAT Winds SCAT Winds First Guess Histogram First Guess Dir. (deg.) Reject Number SCAT Winds RMSE All Observation ー First Guess 20.18(deg.) BIAS Histogram 0.81(deg.) Reject

  17. SCAT WINDS(Dir.) Compared to Fist Guess( 8 – 22 July, ’02) Number SCAT Winds SCAT Winds First Guess Histogram First Guess Reject Dir. (deg.) Number SCAT Winds RMSE All Observation ー First Guess 19.70(deg.) BIAS Histogram 0.06(deg.) Reject

  18. The Results of Comparison of SCAT Winds and First Guess Speed : RMSEaround 1.9 m/s BIAS 0.6 ~ 0.8 m/s Direction: RMSE around 20 deg. BIASalmost zero

  19. A Example of OSE (1 Dec. 2001) ①Effect of SCAT Winds 6 hours Forecast (initial :18UTC Nov. 30) (First Guess) altitude pressure temperature  ・  ・ SCAT Winds+Conventional data Assimilation Conventional data Assimilation CNTL TEST Effect of SCAT Winds(TEST – CNTL)

  20. Effect of altitude at 1000hPa(00UTC 1 Dec. 2001) SCAT Winds Change of altitude and wind at 1000hPa(TEST-CNTL)

  21. Effect of altitude at 500hPa(00UTC 1 Dec. 2001) SCAT Winds Change of altitude and wind at 500hPa(TEST-CNTL)

  22. Effect of temperature at 500hPa(00UTC 1 Dec. 2001) SCAT Winds Change of Temperature at 500hPa(TEST-CNTL)

  23. Result of OSE (CASE 1) Start up : 1 ~ 7 Dec. 2001 Period : 8 ~22 Dec. 2001 9 Days Forecast (every 12UTC) Global Data Assimilation System • Three-dimensional Assimilation System (3D-VAR) Forecast Model T213L40 Global Spectrum Model    (JMA’s Operational Model)

  24. P(surface)(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. RMSE Increase E.Q. S.H.

  25. P(surface)(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. BIAS Decrease BIAS Increase E.Q. S.H.

  26. Wind Speed(surface)(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. E.Q. S.H.

  27. Wind Speed(surface)(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. BIAS Increase E.Q. S.H.

  28. 850hPa Altitude(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. RMSE Increase E.Q. S.H.

  29. 850hPa Altitude(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. BIAS Decrease E.Q. S.H.

  30. 850hPa Temperature(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. E.Q. S.H.

  31. 850hPa Temperature(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. E.Q. S.H.

  32. 500hPa Altitude(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. RMSE Increase E.Q. S.H.

  33. 500hPa Altitude(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. BIAS Decrease E.Q. S.H.

  34. 500hPa Temperature(RMSE) (8 – 22 Dec. 2001) test cntl Global N.H. E.Q. S.H.

  35. 500hPa Temperature(BIAS) (8 – 22 Dec. 2001) test cntl Global N.H. E.Q. S.H.

  36. The Summary of Comparison of TEST and CNTL to Initial Value (Dec. 2001) BIAS N.H. T.R. S.H. Psurf Almost Same Negative Decrease Negative Increase Z850 Almost Same Negative Decrease Almost Same Z500 Almost Same Negative Decrease Almost Same Wsurf Almost Same Negative Increase Almost Same T850 Almost Same Almost Same Almost Same T500 Almost Same Almost Same Almost Same

  37. The Summary of Comparison of TEST and CNTL to Initial Value (Dec. 2001) RMSE N.H. T.R. S.H. Psurf Almost Same Almost Same Slightly Increase Z850 Almost Same Almost Same Slightly Increase Z500 Almost Same Almost Same Slightly Increase Wsurf Almost Same Almost Same Almost Same T850 Almost Same Almost Same Almost Same T500 Almost Same Almost Same Almost Same

  38. Z500hPa Anomaly Correlation (8 – 22 Dec. 2001) N.H. S.H. Accuracy decrease test cntl E.Q.

  39. P(surface) Anomaly Correlation (8 – 22 Dec. 2001) N.H. S.H. Accuracy decrease test cntl Accuracy Increase E.Q.

  40. Comparison of the experimental results with radiosonde observations(Z850, December 2001) N.H. RMSE BIAS TEST CNTL Forecast hours E.Q. S.H. Red: TEST Green: CNTL Scatter diagram

  41. Comparison of the experimental results with radiosonde observations(Z500, December 2001) N.H. RMSE BIAS TEST CNTL Forecast hours E.Q. S.H. Red: TEST Green: CNTL Scatter diagram

  42. Comparison of the experimental results with radiosonde observations(T850, December 2001) N.H. BIAS RMSE TEST CNTL Forecast hours E.Q. S.H. Red: TEST Green: CNTL Scatter diagram

  43. Comparison of the experimental results with radiosonde observations(T500, December 2001) N.H. RMSE BIAS TEST CNTL Forecast hours E.Q. S.H. Red: TEST Green: CNTL Scatter diagram

  44. The Summary of Comparison of TEST and CNTL to Radiosonde (Dec. 2001) RMSE N.H. T.R. S.H. Z850 Slightly Decrease Slightly Decrease Almost Same Z500 Slightly Decrease Almost Same Almost Same T850 Almost Same Almost Same Almost Same T500 Slightly Decrease Slightly Decrease Almost Same

  45. Result of OSE (CASE 2) Start up : 1~7 Jul. 2002 Period : 8~22Jul. 2002 9 Days Forecast (every 12UTC) Global Data Assimilation System • Three-dimensional Assimilation System (3D-VAR) Forecast Model T213L40 Global Spectrum Model    (JMA’s Operational Model)

  46. P(surface)(RMSE) (8 – 22 July 2002) test cntl Global N.H. E.Q. S.H.

  47. P(surface)(BIAS) (8 – 22 July 2002) test cntl Global N.H. BIAS Increase E.Q. S.H.

  48. Wind Speed(surface)(RMSE) (8 – 22 July. 2002) test cntl Global N.H. E.Q. S.H.

  49. Wind Speed(surface)(BIAS) (8 – 22 July. 2002) test cntl Global N.H. BIAS Increase E.Q. S.H.

  50. 850hPa Altitude(RMSE) (8 – 22 July. 2002) test cntl Accuracy Increase Global N.H. E.Q. S.H.

More Related