1 / 59

Monsoon Rainfall Forecasting

Monsoon Rainfall Forecasting. Pankaj Jain IIT Kanpur. Introduction. Monsoon prediction is clearly of great importance for India One would like to make long term prediction, i.e. predict total monsoon rainfall a few weeks or months in advance

dwightb
Download Presentation

Monsoon Rainfall Forecasting

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. Monsoon Rainfall Forecasting Pankaj Jain IIT Kanpur

  2. Introduction • Monsoon prediction is clearly of great importance for India • One would like to make • long term prediction, i.e. predict total monsoon rainfall a few weeks or months in advance • short term prediction, i.e. predict rainfall over different locations a few days in advance

  3. Predicting total monsoon rainfall (June-September) • predicted by using its correlation with observed parameters • The predictors keep changing with time • Several regression and neural network based models are currently available • Indian Met. Dept (IMD) provides statistical forecast in two stages, • March/April • May/June

  4. Model • IMD uses both linear and non-linear regression for their forecast • use ensemble forecast • large number of models are used for all possible combinations of predictors • only a few models with best skill are selected • The forecast is the weighted average of the outcome of these models • The model error • 5% for April forecast • 4% for June forecast

  5. Short term Forecasting • We have been interested in forecasting daily rainfall over a particular location a few days in advance. • The government agency National Center for Medium Range Weather Forecasting (NCMRWF) provides daily forecasts, mainly to assist farmers.

  6. Numerical Weather Prediction (NWP) Models • Numerical Weather Prediction (NWP) models • Used to make short (1-3 days) and medium (4-10) forecast • Navier Stokes equation is written on a spherical grid covering the entire earth • use spherical polar coordinates • need to account for the earth’s rotation, which makes it a non-inertial frame. This introduces fictitious Centrifugal and Coriolis force

  7. The variables are expanded in spherical harmonics, truncated up to a certain multipole, which determines the resolution of the grid. • For example the current model is T254, which implies a grid size of 0.5ox0.5o • 64 vertical levels • 7.5 min time steps

  8. Atmosphere

  9. Long wavelength radiation Kinetic Energy unequal Potential Energy Solar Radiation heating Frictional Dissipation Earth-Atmosphere System

  10. General Circulation

  11. The inputs to the model are the initial conditions obtained by observations throughout the earth temperature, pressure, wind velocity, humidity etc as a function of position and height

  12. The output of the model is the desired prediction However the model is severely limited: • The outcome, especially rainfall, is strongly dependent on local factors • This is particularly true in tropics where the circulation is primarily driven by convection • It is unfeasible to take all local factors into consideration in a global model • The prediction may change considerably by very small changes in the input parameters

  13. The input data, especially high altitude balloon data, is severely limited • Also in many regions, especially India, the data quality is often not very good

  14. There may also exist some unknown effects. An interesting possibility is effect of galactic cosmic rays This possibility has been studied by Tripathi et al (CE, IITK)

  15. Variation of low-altitude cloud cover, galactic cosmic rays and total solar irradiance (1984-1994). The cosmic ray intensity data is from Huancayo observatory, Hawaii Carslaw et al., 2002, Science The physical links for the correlation is subject of research.

  16. Two Routes to cloud modification: Charged Species is the Key! Condensational growth Thermodynamically Stable clusters Molecules Ion-induced nucleation Neutral or ion clusters Initial growth step Aerosol particles Particle-particle coagulation Tripathi and Harrison, 2001; Tripathi et al, 2006 CONDUCTING WATERDROP Drop charge D c Image charge I Cloud drop s Charged Aerosol Trajectory Further growth CCN CHARGED AEROSOL Modgil, Kumar, Tripathi et al., 2005, JGR Charged Aerosol Collision Ion Mediated Rote

  17. Statistical Interpretation of NWF output • It may be better to statistically correlate the model output with observations • This is the technique used by NCMWRF to predict daily rainfall a few days in advance at a particular station • The rainfall at a particular station is obtained by a rain gauge • We have been trying to determine if neural network based relationship can improve the predictability. Jain et al 1999, Jain and Jain 2002

  18. Scale invariance in daily rainfall • The predictability of the quantity can often be judged by its distribution function. • If the variable shows a normal distribution then large fluctuations from the mean value are improbable • However a power law distribution f(x) = x implies no characteristic scale (scale invariant) indicates an underlying chaotic behaviour

  19. Distribution of daily rainfall: Kanpur Frequency Precipitation (0.1mm)

  20. Distribution of Daily Rainfall: Lucknow Frequency Precipitation (0.1mm)

  21. Hourly rainfall Distribution SSM/I satellite data Frequency 5S-10S Rainfall rate in mm/hour A power law fits the distribution very well at low latitudes

  22. Hourly rainfall Distribution Frequency 15N-20N Rainfall rate in mm/hour

  23. Frequency 0N-5N Rainfall rate in mm/hour

  24. Frequency 50N-55N Rainfall rate in mm/hour

  25. Frequency Rainfall rate in mm/hour

  26. short term rainfall over a localized region shows a scale invariant power law distribution Jain and Jain 2002 Peters et al (2002) show this for individual events at the Baltic coast

  27. Power law exponent as a function of the latitude f(x) = x Exponent  Latitude In tropics =1.130.14 At higher latitudes =1.3-1.6 Jain and Jain, 2002

  28. It is better to define variable which we may have a better chance of predicting. • Rather then using a single rain gauge it may be more appropriate to use the rainfall averaged over many rain gauges. • The NWF has a grid size of order 100x100 Km. Hence its predictions should be interpreted as the average over the grid rather than for a particular location.

  29. Predicting Daily Rainfall • We studied daily rainfall forecast one day in advance • The following stations were considered: Delhi, Pune, Hyderabad, Bangalore,Bhubaneshwar • The output variable (y) is the daily rainfall • 47 input variables (xi), each at 9 grid points surrounding the station. select by quadratic fitting over the 9 grid locations • 6 years data (1994-1999) from June to September for Pune, Hyderabad, Bangalore,Bhubaneshwar For Delhi we consider data from June to August

  30. Neural Networks Neural networks differ from statistical regression techniques since one does not try to fit the output. By fitting we mean minimization of the error Sum runs over the training set yi is the predicted variable and y’i the measured variable. Instead one only tries to learn the behaviour of the predictor.

  31. We terminate the training when the error in the validation becomes minimum • Then the results are checked on an independent set. • While training one has to be careful that the network is not struck in some local minima • genetic algorithms or • simulated anealing

  32. Performance Indices We shall predict (a) the probability of rainfall and (b) the actual rainfall (actually the cube root of rainfall) The skill of the model is tested by suitable performance indices A model is skillful if it performs better than persistence model

  33. Probability of Rainfall Performance Indices: Ratio = number correct / total N(dry) = no. of correctly predicted dry days M(dry) = no. of incorrectly predicted dry days

  34. Amount of RainfallCube root of Precipitation Performance Index: Root Mean Square Error

  35. Results (Pune)

  36. Results (Hyderabad)

  37. Results (Bangalore)

  38. Results (Bhubaneshwar)

  39. Results (Delhi)

  40. Results with Average Rainfall over 3 rain gauges in Delhi

  41. Results with Average Rainfall over 3 rain gauges in Delhi

  42. Results with Average Rainfall over 3 rain gauges in Delhi

More Related