Precipitation forecast over western Himalayas using k-nearest neighbour method
|
|
- Julia Blake
- 5 years ago
- Views:
Transcription
1 INTERNATIONAL JOURNAL OF CLIMATOLOGY Int. J. Climatol. 28: (2008) Published online 25 February 2008 in Wiley InterScience ( Precipitation forecast over western Himalayas using k-nearest neighbour method A. P. Dimri, a * P. Joshi b and A. Ganju b a School of Environmental Sciences, Jawaharlal Nehru University, New Delhi, India b Research and Development Center, Snow and Avalanche Study Establishment, HimParisar, Sector 37A, Chandigarh, India ABSTRACT: Knowledge of precipitation over the western Himalayas is of utmost importance during winter season (December, January, February, and March DJFM) due to extreme weather conditions. During winter, eastward moving synoptic weather systems called western disturbances (WDs) yield enormous amount of precipitation over this region. This amount of precipitation by a number of WDs keeps on accumulating and poses an avalanche threat to the habitat. Not only this, low temperature conditions along with precipitation form extremely hostile winter conditions to live with. In the present study, the nearest neighbour (NN) method is used to forecast probability of precipitation (PoP) occurrence/nonoccurrence and its quantity. The method of NN is introduced to visualize the results. The NN forecast technique attempts to compare similar situations in the past with current data and assumes that similar events are likely to occur under similar conditions. At present, only nine important weather variables are considered for generating a 3-day advance forecast of PoP occurrence/non-occurrence and quantity at eight representative observatories in Jammu and Kashmir (J&K), the northmost state of India. This state receives maximum amount of solid precipitation in the form of snow during winter. Past data of 8 14 years (between 1988 and 2003) is considered, from which the nearest days are looked for and tested with data of 2 4 years (between 2003 and 2005). PoP occurrence/non-occurrence is well predicted by k-nearest Neighbour (k-nn) sampling method. It is evident from the results that the model is able to give the projections 3 days in advance, with an accuracy of 71 88%. Probability forecast of occurrence of precipitation is predicted well by the present model setup, but the quantitative precipitation forecast (QPF) is an important issue that still needs improvement, keeping in view the topographical and land-use heterogeneity of the area under study. Copyright 2008 Royal Meteorological Society KEY WORDS precipitation; k-nearest neighbour; probability Received 3 November 2006; Revised 4 December 2007; Accepted 24 December Introduction Northwest India has complex mountain ranges. This region receives high amount of precipitation, in the form of snow, during winter (December, January, February, and March (DJFM). This winter precipitation is mainly attributed to the passage of weather systems called western disturbances (WDs). These are eastward moving low pressure synoptic weather systems that originate over the Mediterranean Sea or mid Atlantic Ocean and travel eastward over Iran, Afghanistan, Pakistan, and India. These weather systems take their southernmost tracks during winter when they pass over northwest India. Further, the strength and vertical extent of these systems are influenced by the position of the mid and upper air tropospheric troughs in the zonal westerlies. A detailed discussion of the synoptic features of the WDs is given by Rao and Srinivasan (1969); Singh (1979) and others. As pointed out by earlier studies (Pisharoty and Desai, 1956; Kalsi, 1980; Kalsi and Haldar, 1992), WDs * Correspondence to: A. P. Dimri, School of Environmental Sciences, Jawaharlal Nehru University, New Delhi, India apdimri@hotmail.com, apdimri@yahoo.com are manifestations of the interaction between the tropics and mid latitude systems, and are associated with extensive sheets of mid- and high-level clouds. Surface weather elements, like precipitation and temperature, associated with these WDs are highly dependent upon local topography and local atmospheric circulations. Also, northern India is composed, in part, of Himalayan mountain ranges having different altitudes and orientations, causing the prevailing weather conditions to be complex. Topography steers atmospheric circulations by thermodynamical and dynamical forcings. Further, the surface heterogeneity of northern India generates numerous meso/microscale circulations in the narrow valleys and rugged hills, which dominate in determining the precipitation and temperature pattern over the region. Topography exerts a strong dynamical forcing on atmospheric circulations, land surface exchanges of momentum, energy, water, and chemical constituents with the atmosphere. In view of these factors, large-scale general circulation models with coarse horizontal resolution cannot properly simulate the horizontal distribution of the precipitation and temperature (Mohanty and Dimri, 2004). In recent years, with the establishment of the National Centre Copyright 2008 Royal Meteorological Society
2 1922 A. P. DIMRI ET AL. for Medium Range Weather Forecasting (NCMRWF) in New Delhi, India, a direct numerical model output operationally provides precipitation forecasts over the western Himalayas with a horizontal model resolution of 1.5 latitude 1.5 longitude. The performance of numerical weather (NWP) models over this region is severely constrained due to complex topography as well as lack of adequate observations. In view of this; large-scale general circulation models with coarse horizontal resolution cannot properly simulate the quantity of precipitation and its horizontal distribution. Hence it is very difficult to predict/simulate such surface weather elements over complex mountainous region by even sophisticated state-of-the-art NWP models. Forecast of probability of precipitation (PoP) and Quantitative Precipitation Forecast (QPF) at a specific site and time provide important guidance for human activity, preparedness for natural hazards (such as avalanches and floods), and forecast management. Various statistical techniques are available to predict PoP and QPF. Glahn and Lowry (1969, 1972) used outputs from NWP models to develop regression models to forecast PoP over different parts of the USA. Paegle (1974) compared the forecast of PoP over different parts of the USA derived from equations stratified with respect to the synoptic weather patterns and equations that were not stratified. It was found that the stratified methods were more accurate. Kriplani and Singh (1986) developed composite charts of probabilities of 24-h rainfall amounts exceeding 2.5 mm and 65 mm, when there is a monsoon depression over India. Upadhyay et al. (1986) developed a method to forecast precipitation by considering the fact that the precipitation rates are directly proportional to the large-scale vertical velocity. Using this method, precipitation rates were computed for specific monsoon depression situations over central parts of India. Kruzinga (1989) compared the forecasting of PoP over the Netherlands using an analogue technique and logistic regression. For a 1 3 days lead time, the regression method performed better than the analogue technique. Kumar and Ram (1995) developed a technique for QPF over the Rapti catchment region in Uttar Pradesh, India. This is a synoptic-analogue method in which synoptic systems are classified according to the observed rainfall rates in the ranges 11 25, and >50 mm. Mohanty et al. (2001) developed objective methods to forecast PoP and QPF at Delhi using classical multivariate regression and discriminant analysis. Maini et al. (2002) employed the perfect prognostic method for precipitation and temperature forecasts during monsoon season. Mohanty and Dimri (2004) and Dimri and Mohanty (2007) have presented performance of statistical downscaling on numerical model outputs of various models and shown enhanced skills by implementing statistical techniques over the complex Himalayan region. Thus, the use of objective techniques to forecast precipitation for a specific location in the mountainous region of northwest India is limited. It may be noted that in mountainous region mainly synoptic, persistence, climatological methods are used to predict the occurrence of the precipitation. Forecasts based on persistence show very poor results, as precipitation has a strong temporal and spatial variability over the region. Therefore, there is a need for more skilful objective methods to predict precipitation over mountains. It is also desirable to introduce statistical relations between precipitation and other observed meteorological parameters using the Nearest Neighbour (NN) method to produce PoP and QPF guidance over the region. Development of NN models for PoP and QPF is a difficult task over rugged mountainous regions like the western Himalayas. This is mainly due to highly heterogeneous terrain and non-availability of adequate observational data sets. Thus, the of occurrence/non occurrence of precipitation is a challenging task over the western Himalayas. Hence, the objective of the present study is to generate PoP occurrence/nonoccurrence and its quantity (QPF) for day 1, day 2 and day 3. In the present study, the k-nearest Neighbour (k-nn) technique is used, which is a flexible approach, developed here based on past observations. Nonetheless, performance of k-nn with historical data sets will depend upon the quality of the observations. In the present work, siteand time- specific of precipitation is carried out by developing an analogue model over the western Himalayas using the k-nn concept. Section 2 describes the methodology of k-nn technique for formulation of the precipitation forecast model. The data source, quality control, and data preparation for the model development over Jammu and Kashmir (J&K), the northernmost region of India (Figure 1), are presented in Section 3. Results and discussion based on this are given in Section 4. Finally, Section 5 presents the main conclusions of the study. 2. The k-nn sampling methodology A continuous NN query retrieves the NN at every point on a line segment (Tao et al., 2002). The basic assumption of NN technique is as follows: if a day like today can be found in the past, we will have a similar precipitation distribution pattern. Certainly, such a likeness must be based on parameters related to the precipitation occurrence and its quantity. In principle, this method makes forecast based on the assumption that for the same weather conditions precipitation occurrence and its quantity will be similar and therefore the method seeks to identify similar situations in the past. Accadia et al. (2003) have used a simple NN average method, also known as remapping of budge interpolation, which maintains total precipitation to a desired degree of accuracy. Buser (1989) and Martin et al. (2000) have extensively used the NN techniques for formulating site specific avalanche forecasting. An annotated algorithm for sampling precipitation value considered with day-to-day dependence on the (nine) variables is presented below. In the present work, the k-nn model is developed, in which nine daily weather elements are considered
3 PRECIPITATION FORECAST OVER WESTERN HIMALAYAS 1923 Figure 1. Schematic representation of geography and topography (m) of the region under study. This figure is available in colour online at as members of the daily weather vector, which is represented by the vector time series of weather variables x t,t = 1,...,n (n = 9 in the present study), as listed in Table I. This weather vector synoptically defines the future state of the precipitation for that particular day. Being a dependent structure, it explains which and how many lags the future state of precipitation will depend on and the number of NN k touse. A vector x t 1 is then defined. This vector of the (nine) variables of interest of the previous day is considered to assess the next day s precipitation (say ppt t ). Now the current day s vector, x i, is used to Table I. Meteorological parameters considered for the study. SNo Parameter Time of observation 1200 UTC (previous day) 0300 UTC 1 Maximum temperature (T x ) 2 Minimum temperature (T n ) 3 Dry bulb 0300 UTC temperature (T ) 4 Av wind speed (ws avg ) 5 Wind direction (wd) 0300 UTC 6 Wind speed (ws) 0300 UTC 7 Surface Pressure (P) 0300 UTC 8 Cloud amount (cla) 0300 UTC 9 Cloud type (clt) 0300 UTC Average of last 24 h determine/select the k-nns among the historical state vectors (x 1,x 2,...,x m,...,x i,...,x t ) using the weighted Hamming distance. To determine this, the distance between individual parameters is calculated and one threshold value is set for each parameter. If the distance is less than this threshold value, then that particular day would be considered as a nearest day. This method is described in greater detail below. N(x) = (y 1,y 2,...,y 9 )/(x 1 θ 1 y 1 x 1 + θ 1 ), (x 2 θ 2 y 2 x 2 + θ 2 ),...,(x 9 θ 9 y 9 x 9 + θ 9 ) (1) where N(x) is the overall vector distance; y 1,y 2,...,y 9 are the vector distances corresponding to variables(say) x 1,x 2,...,x 9 and θ 1, θ 2,...,θ 9 are the threshold values for each variable. These threshold values are chosen in such a way that variability in each variable is accounted for and is schematically described by Figure 2, in which the dependency is shown with two variables. Threshold value is considered in a way that significant contribution of both the variables is looked into. Therefore, ±θ i is considered within the limits of the shaded portion, where simultaneous contribution of both the variables is seen. Similarly, other threshold values are considered in the study. To elaborate further, say, n is the total number of data points, m is the number of actual days for validation and p is the number of days to be tested and θ 1,θ 2,...,θ 9
4 1924 A. P. DIMRI ET AL. +θ 2 is taken into account for the forecast. The number of NNs selected is based on the criteria that evaluate the mean square forecast error (corrected for degree of freedom of the model) to choose model parameters of time series simulation models. e i = x i x f i (5) k x f 1 i = xi,j k j=1 1m (6) x 2 -θ 2 -θ 1 +θ 1 m=1 x 1 Figure 2. A schematic representation of selection of threshold values for individual parameters. ±θ 1 or ±θ 2 represents the most suitable limits of variables, such as relative humidity and temperature, under which maximum occurrence of precipitation is likely (the shaded portion represents the most suitable region of variables). are the threshold values corresponding to each of the parameters. The differences between parameters of an actual day with the same parameters of the past years, i.e. p days, are calculated. Below is shown the distance, y 1 (k) and y 2 (k), calculated for maximum (T x ) and minimum (T n ) temperature respectively. y 1 (k) = T x (j) T x (k) (2) y 2 (k) = T n (j) T n (k) (3) where j varies from m, m + 1,...,n and k is varying from 1, 2,...,p. Similarly, distances y 3 (k), y 4 (k),...,y 9 (k) are calculated for rest of the parameters considered in this study. Now if the distance y 1 (k) is less than θ 1, threshold value for the first parameter, the value is set to be 0 otherwise 1. Let this value be z 1 (k). Thus if y 1 (k) < θ 1 then z 1 (k) = 0otherwisez 1 (k) = 1 Similarly, If y 2 (k) < θ 2 then z 2 (k) = 0otherwisez 2 (k) = 1, and so on. Then the final Hamming distance is calculated as D = z 1 (k) + z 2 (k) + z 3 (k) z 9 (k) (4) Among these days, the days with zero distance will be nearest to the conditioning day and it will be assumed that the synoptic situations of that particular nearest day will be closest and hence corresponding precipitation state of next 24 h can be estimated as that of the nearest days. Finally, a set of k-nn states is selected in which an element of this set records the time t associated with closest historical state with the current vector (x i ). In this manner, sampling is continued to achieve the similar state and hence corresponding precipitation state where, e i is the error vector of the ith vector x i ; x i is the recorded value that is to be forecast; x f i is the k-nn forecast vector corresponding to x i ; xi,j s is the vector corresponding to the j th nearest neighbor of the vector on which we are conditioning to sample a vector corresponding to x i. 3. Model application and performance measure In general, in most of the statistical methods, a statistical dependence or correlation of daily sequence of weather variables with each other on the same day or on the previous day, persistence, is considered. For example, on rainy days, solar radiation and maximum temperature are likely to be lower, whereas wind speed and minimum temperature may be higher than on dry days. Here, precipitation is chosen as the driving variable to be predicted, as it is one of the requirements over the western Himalayas, which has complex variability within it. A precipitation occurrence model based on NN is developed to generate the sequence of dry and wet days and precipitation amount. The precipitation amount on a rainy day may also depend on the wind, the humidity, and the temperature. Consequently, there is a reason to consider dependence of the daily weather processes on more than just precipitation, as has traditionally being done. Typically, daily precipitation is generated independently, and the other variables are generated by conditioning on precipitation events, i.e. whether a day is wet or dry (Jones et al., 1972). Serial dependence of weather variables is defined by fitting auto regression models of order 1 (AR-1) independently to each variable for each period (Nicks and Harp, 1980). Further, a multivariate autoregressive model of order 1 (MAR-1), adds the consideration of dependence across variables (Richardson, 1981). Mohanty and Dimri (2004) have also shown such dependence using multivariate statistical models. The nine important weather variables considered for the study are maximum temperature (T x ), minimum temperature (T n ), dry bulb temperature (T ), 24 h average wind speed (ws avg ), wind direction (wd), wind speed (ws), surface pressure (P), cloud amount (cla) and cloud type (clt) (Table I). Though these variables are recorded twice a day at 0300 UTC and 1200 UTC, in the present study, the variables observed at the time shown in Table I
5 PRECIPITATION FORECAST OVER WESTERN HIMALAYAS 1925 are considered, whereas for 24 h observations, average of these two records are calculated and used. All these parameters are assessed to determine the state of the precipitation in the next 24 h. Therefore, cumulative precipitation amount recorded in 24 h is considered for study. Precipitation is treated as a binary variable here. If measurable precipitation is observed, the binary value of precipitation is set to 1; if no measurable precipitation is observed, the value is set to 0. The threshold value of precipitation is taken as 0.1 cm, which is the least measurable snow depth. The 24-h accumulated snow depths are found to be highly variable. The precipitation reported on a particular day is the accumulated snow depth in the 24 h ending at the reporting time i.e UTC. A schematic diagram pertaining to observation and forecast time is depicted in Figure 3. Snowfall depths are classified firstly into two classes: non-occurrence and occurrence. Further the occurrence class is classified into four categories: Category I: cm; Category II: cm; Category III: cm and Category IV: 48.1 cm. Snow depth (in cm) is used to verify the results in this method because they are the only observations available. Eight representative observatories (road axes/sectors) considered for study are Kanzalwan (Bandipur Gurez axis/bg), Haddantaj (Nogaon Kaiyan axis/nk), Drass (Srinagar Leh axis/sl), StageII (Chowkibal Tangdhar axis/ct), Banihal (National Highway No. IA/NHIA), Pharkian (Keran sector), Gulmarg (Gulmarg sector), and ZGali (Macchal sector) (Figure 4). These observatories are considered since they represent the climatic and weather conditions of the region in which they fall and are also nodal centres for assessing/disseminating warnings of any natural hazards in time. Further, these sectors are considered because they represent different geographical and climatic conditions of the J&K, the northmost state of India. Data from December 1988 to March 2005 is considered for the present study for predicting probability of occurrence/non-occurrence and category of the precipitation for 3 days in advance. Out of this, December 1988 March 2003 data is used as the historical data set from which nearest days are selected Observation Time N-1 N day N+1 Variable considered Precipitation occurrence/ nonoccurrence and amount Figure 3. Reference time, in UTC, for the variables considered in the study. for the current days which are from December 2003 to March Experimental design The PoP occurrence/non-occurrence model is initiated at 0300 UTC to generate a forecast for the next 3 days. The probability of quantity of precipitation model is initiated at the same time only if the first model indicates the precipitation occurrence as yes. The algorithm defined in Section 2 is applied to data from the eight observatories and selected statistics of the simulated state of precipitation are computed. The statistics of interest, as described in Section 3.2, is computed for each simulation and compared with the actual observation Performance measure The following statistics are considered to be of interest in comparing historical record and the simulated record of precipitation state. None of these statistical measures are explicitly specified in fitting the k-nn model. Consequently, their successful reproduction can be considered a sign of success for the method. For the purpose of verification of categorical forecast, contingency tables are prepared and the verification parameters and skill scores are evaluated as defined in the appendices (Wilks, 1995). 4. Results and discussion PoP occurrence/non-occurrence and probability of occurrence of certain category are assessed. Various statistical scores and skill scores are computed on the basis of simulation. Comprehensive analysis of these scores is carried out for day 1, day 2 and day 3 projections Performance of PoP occurrence/non-occurrence model The model for forecasting performance of probability of occurrence/non-occurrence of precipitation at the eight sites is evaluated using a 2 2 contingency table (Wilks, 1995). For the purpose of verification, measures and skill scores are computed (Appendix A). Table II illustrates that the probability of detection (POD) of occurrence of precipitation is higher in case of day 1 as compared to day 2 and day 3 at all the eight stations. This is attributed to the fact that the relation of precipitation events with synoptic situations decreases as the time lag increases. This also occurs because over the western Himalayas, precipitation has large variability owing to topographic variation. False alarm rate (FAR) for all these stations does not exceed 0.56 for day 1, whereas for day 2 and day 3, its value is less than 0.74 and 0.73, respectively. These extreme values are at Drass, which is situated in Great Himalayan range. Drass, as per the geographic location (Figure 4) receives lot of modified weather as WDs which travel eastward first interact/cross Pir Panjal range and then reach the Drass region. By virtue of this, precipitation patterns over Drass
6 1926 A. P. DIMRI ET AL. Figure 4. Observatory and sector network in Jammu and Kashmir. This figure is available in colour online at region are very different from that at the rest of the stations, which are mostly located in the Pir Panjal range. These values of FAR also illustrate that the model has better POD for non-occurrence of precipitation (C-NON) for all the 3 days (ranging from 0.70 to 0.94; and for day 1, day 2, and day 3, respectively) than it does for occurrence of precipitation (POD ranging from 0.50 to 0.92; and for day 1, day 2, and day 3, respectively). This may be attributed to the fact that the number of no-precipitation days is generally higher compared to the number of days with precipitation (Table III). The model bias shows the stability in the model performance with values ranging from 0.90 to 1.23; , and for day 1, day 2 and day 3, respectively at all the stations except at Banihal, where it shows low values. Also, the model performance is degraded as the lead time increases. Now, this clearly shows the fact that the model was not able to capture the synoptic weather observations satisfactorily for the forecast for day 3. Nonetheless, the model
7 PRECIPITATION FORECAST OVER WESTERN HIMALAYAS 1927 Table II. Verification measures (Appendix A) for the probability of precipitation occurrence/non-occurrence model. Measure Day 1 Day 2 Day 3 Kanzalwan POD FAR MR C-NON CSI TSS BIAS PC 74% 54% 44% Haddantaj POD FAR MR C-NON CSI TSS BIAS PC 88% 68% 58% Drass POD FAR MR C-NON CSI TSS BIAS PC 74% 64% 64% Stage2 POD FAR MR C-NON CSI TSS BIAS PC 76% 65% 60% Banihal POD FAR MR C-NON CSI TSS BIAS PC 82% 73% 69% Pharkian POD FAR MR C-NON CSI TSS BIAS PC 71% 59% 52% Table II. (Continued). Measure Day 1 Day 2 Day 3 Gulmarg POD FAR MR C-NON CSI TSS BIAS PC 76% 62% 57% Zgali POD FAR MR C-NON CSI TSS BIAS PC 81% 49% 49% performance is reasonable for day1 and day2 s. The other skill scores/evaluation indices for the model show reasonably high values towards perfect forecast criteria (Appendix A) particularly for day 1 s. It is interesting to note that the overall performance of the model, which is measured by percent correct (PC), is very high for day 1 at all the eight stations ranging from 71 to 88%. All the forecast verification indices estimated from the contingency table clearly demonstrate that at all the eight stations the model provides very satisfactory performances for day 1, whereas definite improvements for day 2 and day 3 are needed. The general deterioration in day 2 and day 3 values may be because the synoptic scale conditions considered in the study are taken at a 24 h interval. However, sometimes it may happen that complex topography play an important role in modifying the mountain weather within the timescale of 24 h. Owing to non-consideration of this variability, occurrence of precipitation events is overlooked. Also, looking into the philosophy of the NN sampling technique, similar synoptic conditions in the past can yield different precipitation states. So, overall, these aspects and variability in the nature of precipitation occurrence will definitely affect model performance Performance of probability of quantitative precipitation forecast We refer to Section 3 for the definition of snowfall categories. This model, is based on the NN sampling method, is sued to predict probability of a categorical quantity of precipitation. The skill scores and the other verification measures are calculated using 4 4 contingency (Table BI, Wilks, 1995). The skill scores/evaluation indices and other verification measures of this model are illustrated in Table IV for day 1, day 2 and day 3 for the eight stations.
8 1928 A. P. DIMRI ET AL. Table III. Comparison of the number of days with no precipitation with the number of days with precipitation at each station. Station Total no of days considered for study No of days with no precipitation No of days with precipitation Kanzalwan 1715 ( ) Haddantaj 1647 ( ) Drass 1399 ( ) Stage ( ) Banihal 763 ( ) Pharkian 1313 ( ) Gulmarg 2111 ( ) Zgali 1107 ( ) Table IV. Verification measures (Appendix B) for the probability of Quantitative Precipitation Forecast (QPF) model. Measure Day 1 Day 2 Day 3 Kanzalwan CSI I II III IV HSS Percentage Correct 47% 55% 55% % for non-occurrence 72% 62% 62% Haddantaj CSI I II III IV HSS Percentage Correct 79% 49% 49% % for non-occurrence 77% 75% 75% Drass CSI I II III IV HSS Percentage Correct 70% 63% 84% % for non-occurrence 79% 78% 80% Stage2 CSI I II III IV HSS Percentage Correct 45% 35% 55% % for non-occurrence 76% 74% 76% Banihal CSI I II III IV Table IV. (Continued). Measure Day 1 Day 2 Day 3 HSS Percentage Correct 76% 62% 57% % for non-occurrence 94% 91% 85% Pharkian CSI I II III IV HSS Percentage Correct 49% 46% 37% % for non-occurrence 70% 68% 61% Gulmarg CSI I II III IV HSS Percentage Correct 49% 41% 56% % for non-occurrence 83% 76% 78% Zgali CSI I II III IV HSS Percentage Correct 39% 54% 40% % for non-occurrence 77% 51% 52% The Critical Success Index (CSI) of the model in assessing the PoP is higher in category I as compared to other categories for all three days at all the eight stations. This illustrates the fact that the model could capture the precipitation amounts in lower snowfall categories better than the higher categories. In the present model, some of the cases show higher CSI for day 2 and day 3 as compared to day 1. This is attributed to the fact that once
9 PRECIPITATION FORECAST OVER WESTERN HIMALAYAS 1929 the same synoptic situations are found correctly in the past, a category forecast upto 3 days in advance can be made with satisfactory results. Heidke skill score (HSS) ranges from 0.07 to 0.52; and for day 1, day 2 and day 3 respectively at all the eight stations. It is interesting to note that the overall performance in terms of PC ranges from 39 to 79%; 35 79% and 37 84% for day 1, day 2 and day 3, respectively, at all the stations. All the forecasts estimated from NN sampling method for of probability of category of precipitation amount at all the eight stations are able to provide reasonable performance over this data-sparse region. The complexity associated with accurate precipitation occurs mainly for the following reasons. First, the difference may lie in the fact that topography and vegetation cover is very fast changing within a kilometres over the Himalayan region. This variability may be attributed to the fact that although in general topographically induced cold season precipitation maxima are well captured by the model, but the corresponding peak precipitation values are somewhat overestimated. This is evidently a problem in the analogue model relative to numerical models where resolution over the complex topographical Himalayan region does come into play. Secondly, spatial variability of precipitation increases substantially with the effect of land surface heterogeneity due to temperature. Therefore, precipitation can be in the form of snowfall over the higher peaks and rainfall over the valleys. As a result, snow tends to accumulate over the high resolution peaks and melt more effectively over the corresponding valleys over complex Himalayan topography. This can be attributed to the inherent nonlinear nature of snow forming processes. As the temperature threshold for snow formation is reached, say at the high elevation of a peak, snow starts accumulating. Because snow has a higher albedo than a bare soil or vegetation, the overall surface albedo increases and this causes a decrease in absorption of solar radiation at the surface, further bringing down temperatures and allowing for an increase in snow accumulation. 5. Conclusions In this article a k-nn method is used and tested for precipitation forecast over the Himalayan region where both topography and land-use variability is large. In the present multivariate setting, neighbours of the conditioning point correspond to data patterns that are similar to the patterns at the conditioning point. And it was found that for a day with no rain, which is warm, with little wind, and no humidity, neighbours established by calculating the vector distance between the observations will be similar days. The values for the precipitation for the next day are sampled as a vector from a historically similar day. Clearly, this method has shown strong utility in giving a probability forecast to a day that is more similar to the conditioning day than the other neighbours. Using a weighted Hamming distance function, sensitivity in the developed model has increased for selecting the number of NNs for precipitation. The model could reproduce the day 1 forecast with reasonable accuracy as a short range statistical tool. This technique mainly focuses on short range statistical properties and suffers from the drawback of assuming the data to be normally distributed. Therefore, only linear dependence between surface variables and precipitation states from one day to the next is reproduced. Also this model is unlikely to reproduce the variance and related attributes at longer periods (e.g. the inter-annual variance and dependence of seasonal precipitation). In this article, only a winter (DJFM) is considered for study. The results might depend on the specific simulated year, which are not considered in the study. In addition to this, fine-scale heterogeneity of topography and land use strongly affect the large-scale circulations and surface features, and hence suitable interpolations schemes are necessarily to be considered in the study to obtain data over the data-sparse region. Therefore, some of these aspects are proposed for introduction into the k-nn model for future study. Also, considering more variables like wind fields, vorticity, and divergence from the numerical model output will be more helpful to assess larger variability of the method. Apart from this, suitable assumption and selection of threshold values for each parameter will take care of regional and geographical precipitation distribution. Keeping in view the operational use of NWP models over the Himalayan region, the k-nn technique can potentially be used for integration on to the NWP numerical model outputs for statistical downscaling. Since the Himalayan region has a large topographical and land-use variability, simulating NWP models with very fine resolution needs extensive computational power and data assimilation in initial and boundary conditions. This integration will be able to consider effects of local meso/microscale circultation that can be considered as future study. Acknowledgements The authors acknowledge Dr R.N. Sarwade, Director, SASE, for encouraging the study. Thanks are also due to various persons of Avalanche Forecasting Group (AFG), SASE, who have done a tremendous job in the complex Himalayan region during winter. Appendix A The values in the contingency table are defined as follows, 1. When an event is predicted to occur (forecast occurrence) and in reality it does occur (observed occurrence), then it is classified as A, otherwise (observed non-occurrence) it is classified as C.
10 1930 A. P. DIMRI ET AL. Table AI. Verification measures used for forecast evaluation. Observed Yes Forecast No Yes A B No C D 2. When an event is predicted not to occur (forecast non-occurrence) and in reality it does occur (observed occurrence), then it is classified as B, otherwise (observed non-occurrence) it is classified as D. 3. A + B: Total number of cases of occurrence of precipitation as observed. 4. C + D: Total number of cases of non-occurrence of precipitation as observed. 5. A + B + C + D: Total number of forecasts. For a best/perfect forecast series: B = 0andC= 0 and hence POD = 1, FAR = 0, MR = 0, C-NON = 0, Bias = 1, CSI = 1, TSS = 1, HSS = 1, PC = 100% Appendix B Table BI. Categorical verification of forecasts (four category events). Observed Forecast I II III IV Total I a b c d J II e f g h K III i j k l L IV m n o p M Total N O P Q T Probability of Detection (POD) POD = A A + B False Alarm Rate (FAR) FAR = C C + A Miss Rate (MR) MR = B B + A Correct non-occurrence (C-NON) C NON = D D + C Critical Success Index (CSI) A CSI = A + B + C True Skill Score (TSS) TSS = A A + B + D D + C 1 Heidke Skill Score (HSS) 2(AD BC) HSS = B 2 + C 2 + 2AD + (B + C)(A + D) Bias (BIAS) for occurrence BIAS = A + C A + B Percentage Correct (PC) PC = A + D A + B + C + D 100% (A1) (A2) (A3) (A4) (A5) (A6) (A7) (A8) (A9) (Category I: cm; Category II: cm; Category III: cm and Category IV: 48.1 cm) Total number of observed events in category I is J = a + b + c + d Total number of forecast events in category I is N = a + e + i + m (B1) (B2) In the similar wayo,k,p,l,q and M are computed. Then the total numbers of events are T = J + K + L + M = N + O + P + Q Percentage Correct (PC) (B3) PC = a + f + k + p 100% (B4) T Critical Success Index (CSI) CSI = Heidke Skill Score (HSS) a J + N a, f K + O f, k L + P k, p M + Q p (B5) JN + KO + LP + MQ a + f + k + p HSS = T JN + KO + LP + MQ T T (B6) References Accadia C, Mariani S, Casaioli M, Lavagnini A Sensitivity of precipitation forecast skill scores to Bilinear interpolation and
11 PRECIPITATION FORECAST OVER WESTERN HIMALAYAS 1931 a simple Nearest-Neighbor Average Method on High-resolution Verification Grids. Weather and Forecasting 13: Buser O Two years experience of operational avalanche forecasting using nearest neighborhood method. Annals of Glaciology 13: Dimri AP, Mohanty UC Location-specific of maximum and minimum temperature over the Western Himalayas. Meteorological Applications 14: Glahn HR, Lowry DA An operational method for objectively forecasting probability of precipitation. ESSA Technical Memorandum WBTM 27: 24. Glahn HR, Lowry DA The use of model output statistics (MOS) in objective weather forecasting. Journal of Applied Meteorology 11: Jones W, Rex RC, Threadgill DE A simulated environmental model of temperature, evaporation, rainfall and soil moisture. Trans ASAE 15: Kalsi SR On some aspects of interaction between middle latitude westerlies and monsoon circulation. Mausam 38: Kalsi SR, Haldar SR Satellite observations of interaction between tropics and mid latitude. Mausam 43: Kriplani RH, Singh SV Rainfall probabilities and amounts associated with monsoon depressions over India. Mausam 37: Kruzinga S Statistical interpretation of ECMWF products in Dutch weather service. Seminar/Workshop on Interpretation of NWP Products, Reading, European Centre for Medium Range Weather Forecasting, Kumar A, Ram LC Semi-quantitative forecasts for Rapti catchment by synoptic analogue method. Mausam 46: Maini P, Kumar A, Singh SV, Rathore LS Statistical interpretation of NWP products in India. Meteorological Applications 9: Martin G, Hans JE, Karl B, Tom L NXD 2000 an improved avalanche forecasting programme based on nearest neighbor method. ISSW, Mohanty UC, Dimri AP Location specific of probability of occurrence and quantity of precipitation over Western Himalayas. Weather and Forecasting 19(3): Mohanty UC, Ravi N, Madan OP Forecasting precipitation over Delhi during the southwest monsoon season. Meteorological Applications 8: Nicks AD, Harp JF Stochastic generation of temperature and solar radiation data. Journal of Hydrology 48: 1 7. Paegle JN Prediction of precipitation probability based on 500 mb flow types. Journal of Applied Meteorology 13: Pisharoty P, Desai BN Western disturbances and Indian weather. Indian Journal of Meteorology and Geophysics 7: Rao YP, Srinivasan V Forecasting manual, Part II: Discussion of typical synoptic weather situation: winter western disturbances and their associated feature. FMU report no III-1.1. India meteorological Department: India,, issued by. Richardson CW Stochastic simulation of daily precipitation, temperature and solar radiation. Water Resources Research 17(1): Singh MS Westerly upper air troughs and development of western disturbances over India. Mausam 30(4): Tao Y, Papadias D, Shen Q Continuous Nearest neighbor search. In Proceedings of the 28 th VLDB Conference, Hong Kong. Upadhyay DS, Apte NY, Kaur S, Singh SP A dynamical approach to quantitative precipitation forecast. Mausam 37: Wilks DS Statistical Methods in Atmospheric Sciences. Academic Press; San Diego, 466.
Prediction of western disturbances and associated weather over Western Himalayas
Prediction of western disturbances and associated weather over Western Himalayas H. R. Hatwar*, B. P. Yadav and Y. V. Rama Rao India Meteorological Department, Lodi Road, New Delhi 110 003, India Two cases
More informationAnalysis of Heavy Precipitation Events Over Western Himalayan Region
INTERNATIONAL JOURNAL OF EARTH AND ATMOSPHERIC SCIENCE Journal homepage: www.jakraya.com/journal/ijeas Analysis of Heavy Precipitation Events Over Western Himalayan Region B.P. Yadav, Naresh Kumar and
More informationSite-specific analog weather-forecast system for northwest Himalaya, India
224 Annals of Glaciology 49 2008 Site-specific analog weather-forecast system for northwest Himalaya, India Dan SINGH, Amreek SINGH, Ashwagosha GANJU Research and Development Centre, Snow and Avalanche
More informationStatistical interpretation of NWP products in India
Meteorol. Appl. 9, 21 31 (2002) Statistical interpretation of NWP products in India Parvinder Maini, Ashok Kumar, S V Singh and L S Rathore, National Center for Medium Range Weather Forecasting, Department
More informationClimate-change studies in the western Himalaya
Annals of Glaciology 51(54) 2010 105 Climate-change studies in the western Himalaya M.S. SHEKHAR, H. CHAND, S. KUMAR, K. SRINIVASAN, A. GANJU Research and Development Centre, Snow and Avalanche Study Establishment
More informationSuppressed Surface Heating over northwestern parts of the Country and the ensuing Monsoon Rainfall
Suppressed Surface Heating over northwestern parts of the Country and the ensuing Monsoon Rainfall During this year, the temperatures have been below normal over north and northwestern parts of India.
More informationOrography and Monsoons: winter snow-storms over the Himalayas
Orography and Monsoons: winter snow-storms over the Himalayas Massimo Bollasina1 and Laura Bertolani2 1Experimental 2Epson Climate Prediction Center, CA Meteo Center - Milan CEOP/IGWCO Joint Meeting, Tokyo,
More informationARPS MODEL SIMULATION OF SYNOPTIC SCALE WEATHER EVENTS
ARPS MODEL SIMULATION OF SYNOPTIC SCALE WEATHER EVENTS Satish Prakash 1, *R.K. Giri 2 and Satish Chand Agarwal 3 1,3 Department of Physics, Meerut College Meerut-UP-25001-India 2 India Meteorological Department,
More informationUsefulness of mesoscale weather forecast for avalanche forecasting
Usefulness of mesoscale weather forecast for avalanche forecasting K. Srinivasan*, A. Ganju and S. S. Sharma Snow and Avalanche Study Establishment, Research and Development Centre, Chandigarh 160 036,
More informationOutgoing long wave radiation (OLR) a proxy of convection
Outgoing long wave radiation (OLR) a proxy of convection Vijay Garg M.M.College, Modinagar, Ghaziabad Abstract: Outgoing long wave radiation (OLR) is a proxy indicator of convection. It is measured through
More informationNEW SCHEME TO IMPROVE THE DETECTION OF RAINY CLOUDS IN PUERTO RICO
NEW SCHEME TO IMPROVE THE DETECTION OF RAINY CLOUDS IN PUERTO RICO Joan Manuel Castro Sánchez Advisor Dr. Nazario Ramirez UPRM NOAA CREST PRYSIG 2016 October 7, 2016 Introduction A cloud rainfall event
More informationApplication and verification of ECMWF products 2008
Application and verification of ECMWF products 2008 RHMS of Serbia 1. Summary of major highlights ECMWF products are operationally used in Hydrometeorological Service of Serbia from the beginning of 2003.
More informationApplication and verification of the ECMWF products Report 2007
Application and verification of the ECMWF products Report 2007 National Meteorological Administration Romania 1. Summary of major highlights The medium range forecast activity within the National Meteorological
More informationApplication and verification of ECMWF products 2015
Application and verification of ECMWF products 2015 Hungarian Meteorological Service 1. Summary of major highlights The objective verification of ECMWF forecasts have been continued on all the time ranges
More informationNumerical Simulation of a Severe Thunderstorm over Delhi Using WRF Model
International Journal of Scientific and Research Publications, Volume 5, Issue 6, June 2015 1 Numerical Simulation of a Severe Thunderstorm over Delhi Using WRF Model Jaya Singh 1, Ajay Gairola 1, Someshwar
More informationApplication and verification of ECMWF products 2016
Application and verification of ECMWF products 2016 Icelandic Meteorological Office (www.vedur.is) Bolli Pálmason and Guðrún Nína Petersen 1. Summary of major highlights Medium range weather forecasts
More informationCombining Deterministic and Probabilistic Methods to Produce Gridded Climatologies
Combining Deterministic and Probabilistic Methods to Produce Gridded Climatologies Michael Squires Alan McNab National Climatic Data Center (NCDC - NOAA) Asheville, NC Abstract There are nearly 8,000 sites
More informationForecast Verification Analysis of Rainfall for Southern Districts of Tamil Nadu, India
International Journal of Current Microbiology and Applied Sciences ISSN: 2319-7706 Volume 6 Number 5 (2017) pp. 299-306 Journal homepage: http://www.ijcmas.com Original Research Article https://doi.org/10.20546/ijcmas.2017.605.034
More informationApplication and verification of ECMWF products 2010
Application and verification of ECMWF products 2010 Icelandic Meteorological Office (www.vedur.is) Guðrún Nína Petersen 1. Summary of major highlights Medium range weather forecasts issued at IMO are mainly
More informationRomanian Contribution in Quantitative Precipitation Forecasts Project
3 Working Group on Physical Aspects 29 Romanian Contribution in Quantitative Precipitation Forecasts Project Rodica Dumitrache, Victor Pescaru, Liliana Velea, Cosmin Barbu National Meteorological Administration,
More informationClimate Modeling: From the global to the regional scale
Climate Modeling: From the global to the regional scale Filippo Giorgi Abdus Salam ICTP, Trieste, Italy ESA summer school on Earth System Monitoring and Modeling Frascati, Italy, 31 July 11 August 2006
More informationApplication and verification of ECMWF products 2014
Application and verification of ECMWF products 2014 Israel Meteorological Service (IMS), 1. Summary of major highlights ECMWF deterministic runs are used to issue most of the operational forecasts at IMS.
More informationSurface energy balance of seasonal snow cover for snow-melt estimation in N W Himalaya
Surface energy balance of seasonal snow cover for snow-melt estimation in N W Himalaya Prem Datt, P K Srivastava, PSNegiand P K Satyawali Snow and Avalanche Study Establishment (SASE), Research & Development
More information3. HYDROMETEROLOGY. 3.1 Introduction. 3.2 Hydro-meteorological Aspect. 3.3 Rain Gauge Stations
3. HYDROMETEROLOGY 3.1 Introduction Hydrometeorology is a branch of meteorology and hydrology that studies the transfer of water and energy between the land surface and the lower atmosphere. Detailed hydrological
More informationONE-YEAR EXPERIMENT IN NUMERICAL PREDICTION OF MONTHLY MEAN TEMPERATURE IN THE ATMOSPHERE-OCEAN-CONTINENT SYSTEM
71 4 MONTHLY WEATHER REVIEW Vol. 96, No. 10 ONE-YEAR EXPERIMENT IN NUMERICAL PREDICTION OF MONTHLY MEAN TEMPERATURE IN THE ATMOSPHERE-OCEAN-CONTINENT SYSTEM JULIAN ADEM and WARREN J. JACOB Extended Forecast
More informationApplication and verification of ECMWF products: 2010
Application and verification of ECMWF products: 2010 Hellenic National Meteorological Service (HNMS) F. Gofa, D. Tzeferi and T. Charantonis 1. Summary of major highlights In order to determine the quality
More information(Statistical Forecasting: with NWP). Notes from Kalnay (2003), appendix C Postprocessing of Numerical Model Output to Obtain Station Weather Forecasts
35 (Statistical Forecasting: with NWP). Notes from Kalnay (2003), appendix C Postprocessing of Numerical Model Output to Obtain Station Weather Forecasts If the numerical model forecasts are skillful,
More informationEnhancing Weather Information with Probability Forecasts. An Information Statement of the American Meteorological Society
Enhancing Weather Information with Probability Forecasts An Information Statement of the American Meteorological Society (Adopted by AMS Council on 12 May 2008) Bull. Amer. Meteor. Soc., 89 Summary This
More informationVerification of the Seasonal Forecast for the 2005/06 Winter
Verification of the Seasonal Forecast for the 2005/06 Winter Shingo Yamada Tokyo Climate Center Japan Meteorological Agency 2006/11/02 7 th Joint Meeting on EAWM Contents 1. Verification of the Seasonal
More informationRAL Advances in Land Surface Modeling Part I. Andrea Hahmann
RAL Advances in Land Surface Modeling Part I Andrea Hahmann Outline The ATEC real-time high-resolution land data assimilation (HRLDAS) system - Fei Chen, Kevin Manning, and Yubao Liu (RAL) The fine-mesh
More informationMozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1
UNDP Climate Change Country Profiles Mozambique C. McSweeney 1, M. New 1,2 and G. Lizcano 1 1. School of Geography and Environment, University of Oxford. 2.Tyndall Centre for Climate Change Research http://country-profiles.geog.ox.ac.uk
More informationAllison Monarski, University of Maryland Masters Scholarly Paper, December 6, Department of Atmospheric and Oceanic Science
Allison Monarski, University of Maryland Masters Scholarly Paper, December 6, 2011 1 Department of Atmospheric and Oceanic Science Verification of Model Output Statistics forecasts associated with the
More informationAtmospheric processes leading to extreme flood events
Atmospheric processes leading to extreme flood events A. P. Dimri School of Environmental Sciences Jawaharlal Nehru University New Delhi, India apdimri@hotmail.com Observational and modelling limitations
More informationECMWF 10 th workshop on Meteorological Operational Systems
ECMWF 10 th workshop on Meteorological Operational Systems 18th November 2005 Crown copyright 2004 Page 1 Monthly range prediction products: Post-processing methods and verification Bernd Becker, Richard
More informationApplication and verification of ECMWF products in Norway 2008
Application and verification of ECMWF products in Norway 2008 The Norwegian Meteorological Institute 1. Summary of major highlights The ECMWF products are widely used by forecasters to make forecasts for
More informationSNOW COVER MAPPING USING METOP/AVHRR AND MSG/SEVIRI
SNOW COVER MAPPING USING METOP/AVHRR AND MSG/SEVIRI Niilo Siljamo, Markku Suomalainen, Otto Hyvärinen Finnish Meteorological Institute, P.O.Box 503, FI-00101 Helsinki, Finland Abstract Weather and meteorological
More informationApplication and verification of ECMWF products 2009
Application and verification of ECMWF products 2009 Icelandic Meteorological Office (www.vedur.is) Gu rún Nína Petersen 1. Summary of major highlights Medium range weather forecasts issued at IMO are mainly
More informationProjected Change in Climate Under A2 Scenario in Dal Lake Catchment Area of Srinagar City in Jammu and Kashmir
Current World Environment Vol. 11(2), 429-438 (2016) Projected Change in Climate Under A2 Scenario in Dal Lake Catchment Area of Srinagar City in Jammu and Kashmir Saqib Parvaze 1, Sabah Parvaze 2, Sheeza
More informationComparison of Interpolation Methods for Precipitation Data in a mountainous Region (Upper Indus Basin-UIB)
Comparison of Interpolation Methods for Precipitation Data in a mountainous Region (Upper Indus Basin-UIB) AsimJahangir Khan, Doctoral Candidate, Department of Geohydraulicsand Engineering Hydrology, University
More informationFog Hazards in Punjab
Pakistan Journal of Meteorology Vol. 8, Issue 16: Jan, 2012 Fog Hazards in Punjab Aslam, M. 1 Abstract Using 30-year (1976-2005) real meteorological data of 10-weather stations, four dense fog years (1998-1999,
More informationA High Resolution Daily Gridded Rainfall Data Set ( ) for Mesoscale Meteorological Studies
National Climate Centre Research Report No: 9/2008 A High Resolution Daily Gridded Rainfall Data Set (1971-2005) for Mesoscale Meteorological Studies M. Rajeevan and Jyoti Bhate National Climate Centre
More informationVerification of the operational NWP models at DWD - with special focus at COSMO-EU
Verification of the operational NWP models at DWD - with special focus at COSMO-EU Ulrich Damrath Ulrich.Damrath@dwd.de Ein Mensch erkennt (und das ist wichtig): Nichts ist ganz falsch und nichts ganz
More informationEvaluation of the Version 7 TRMM Multi-Satellite Precipitation Analysis (TMPA) 3B42 product over Greece
15 th International Conference on Environmental Science and Technology Rhodes, Greece, 31 August to 2 September 2017 Evaluation of the Version 7 TRMM Multi-Satellite Precipitation Analysis (TMPA) 3B42
More informationUse of the models Safran-Crocus-Mepra in operational avalanche forecasting
Use of the models Safran-Crocus-Mepra in operational avalanche forecasting Coléou C *, Giraud G, Danielou Y, Dumas J-L, Gendre C, Pougatch E CEN, Météo France, Grenoble, France. ABSTRACT: Avalanche forecast
More informationDiagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS)
Diagnosing the Climatology and Interannual Variability of North American Summer Climate with the Regional Atmospheric Modeling System (RAMS) Christopher L. Castro and Roger A. Pielke, Sr. Department of
More informationDownscaling in Time. Andrew W. Robertson, IRI. Advanced Training Institute on Climate Variability and Food Security, 12 July 2002
Downscaling in Time Andrew W. Robertson, IRI Advanced Training Institute on Climate Variability and Food Security, 12 July 2002 Preliminaries Crop yields are driven by daily weather variations! Current
More informationPeter P. Neilley. And. Kurt A. Hanson. Weather Services International, Inc. 400 Minuteman Road Andover, MA 01810
6.4 ARE MODEL OUTPUT STATISTICS STILL NEEDED? Peter P. Neilley And Kurt A. Hanson Weather Services International, Inc. 400 Minuteman Road Andover, MA 01810 1. Introduction. Model Output Statistics (MOS)
More informationPrecipitation verification. Thanks to CMC, CPTEC, DWD, ECMWF, JMA, MF, NCEP, NRL, RHMC, UKMO
Precipitation verification Thanks to CMC, CPTEC, DWD, ECMWF, JMA, MF, NCEP, NRL, RHMC, UKMO Outline 1) Status of WGNE QPF intercomparisons 2) Overview of the use of recommended methods for the verification
More informationSevere Freezing Rain in Slovenia
Severe Freezing Rain in Slovenia Janez Markosek, Environmental Agency, Slovenia Introduction At the end of January and at the beginning of February 2014, severe and long-lasting freezing rain affected
More informationApplication and verification of ECMWF products 2016
Application and verification of ECMWF products 2016 RHMS of Serbia 1 Summary of major highlights ECMWF forecast products became the backbone in operational work during last several years. Starting from
More informationGeography Class XI Fundamentals of Physical Geography Section A Total Periods : 140 Total Marks : 70. Periods Topic Subject Matter Geographical Skills
Geography Class XI Fundamentals of Physical Geography Section A Total Periods : 140 Total Marks : 70 Sr. No. 01 Periods Topic Subject Matter Geographical Skills Nature and Scope Definition, nature, i)
More informationHigh spatial resolution interpolation of monthly temperatures of Sardinia
METEOROLOGICAL APPLICATIONS Meteorol. Appl. 18: 475 482 (2011) Published online 21 March 2011 in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/met.243 High spatial resolution interpolation
More informationVerification and performance measures of Meteorological Services to Air Traffic Management (MSTA)
Verification and performance measures of Meteorological Services to Air Traffic Management (MSTA) Background Information on the accuracy, reliability and relevance of products is provided in terms of verification
More informationMesoscale meteorological models. Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen
Mesoscale meteorological models Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen Outline Mesoscale and synoptic scale meteorology Meteorological models Dynamics Parametrizations and interactions
More informationSynoptic Analysis of Total Rainfall Patterns at Azerbaijan District.
Synoptic Analysis of Total Rainfall Patterns at Azerbaijan District Samad Vahdati 1, Shahrokh Shahrokhi Shirvani 2, Abolfazl Nazari Giglou 3 1,3 Department of Civil Engineering, Parsabad Moghan Branch,
More informationEvidence for Weakening of Indian Summer Monsoon and SA CORDEX Results from RegCM
Evidence for Weakening of Indian Summer Monsoon and SA CORDEX Results from RegCM S K Dash Centre for Atmospheric Sciences Indian Institute of Technology Delhi Based on a paper entitled Projected Seasonal
More informationAnalysis of Rainfall and Other Weather Parameters under Climatic Variability of Parbhani ( )
International Journal of Current Microbiology and Applied Sciences ISSN: 2319-7706 Volume 7 Number 06 (2018) Journal homepage: http://www.ijcmas.com Original Research Article https://doi.org/10.20546/ijcmas.2018.706.295
More information8-km Historical Datasets for FPA
Program for Climate, Ecosystem and Fire Applications 8-km Historical Datasets for FPA Project Report John T. Abatzoglou Timothy J. Brown Division of Atmospheric Sciences. CEFA Report 09-04 June 2009 8-km
More informationThe Australian Operational Daily Rain Gauge Analysis
The Australian Operational Daily Rain Gauge Analysis Beth Ebert and Gary Weymouth Bureau of Meteorology Research Centre, Melbourne, Australia e.ebert@bom.gov.au Daily rainfall data and analysis procedure
More informationPrecipitation in climate modeling for the Mediterranean region
Precipitation in climate modeling for the Mediterranean region Simon Krichak Dept. of Geophysics Atmospheric and Planetary Sciences, Tel Aviv University, Israel Concepts for Convective Parameterizations
More informationConvective scheme and resolution impacts on seasonal precipitation forecasts
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 20, 2078, doi:10.1029/2003gl018297, 2003 Convective scheme and resolution impacts on seasonal precipitation forecasts D. W. Shin, T. E. LaRow, and S. Cocke Center
More informationChapter-1 Introduction
Modeling of rainfall variability and drought assessment in Sabarmati basin, Gujarat, India Chapter-1 Introduction 1.1 General Many researchers had studied variability of rainfall at spatial as well as
More informationTC/PR/RB Lecture 3 - Simulation of Random Model Errors
TC/PR/RB Lecture 3 - Simulation of Random Model Errors Roberto Buizza (buizza@ecmwf.int) European Centre for Medium-Range Weather Forecasts http://www.ecmwf.int Roberto Buizza (buizza@ecmwf.int) 1 ECMWF
More informationDevelopment of Pakistan s New Area Weighted Rainfall Using Thiessen Polygon Method
Pakistan Journal of Meteorology Vol. 9, Issue 17: July 2012 Technical Note Development of Pakistan s New Area Weighted Rainfall Using Thiessen Polygon Method Faisal, N. 1, 2, A. Gaffar 2 ABSTRACT In this
More informationInfluence of 3D Model Grid Resolution on Tropospheric Ozone Levels
Influence of 3D Model Grid Resolution on Tropospheric Ozone Levels Pedro Jiménez nez, Oriol Jorba and José M. Baldasano Laboratory of Environmental Modeling Technical University of Catalonia-UPC (Barcelona,
More informationProceedings, International Snow Science Workshop, Banff, 2014
SIMULATION OF THE ALPINE SNOWPACK USING METEOROLOGICAL FIELDS FROM A NON- HYDROSTATIC WEATHER FORECAST MODEL V. Vionnet 1, I. Etchevers 1, L. Auger 2, A. Colomb 3, L. Pfitzner 3, M. Lafaysse 1 and S. Morin
More informationPrecipitation processes in the Middle East
Precipitation processes in the Middle East J. Evans a, R. Smith a and R.Oglesby b a Dept. Geology & Geophysics, Yale University, Connecticut, USA. b Global Hydrology and Climate Center, NASA, Alabama,
More informationdisturbances in the subtropical jetstream & mechanisms associated with extreme rainfall in South Asia
disturbances in the subtropical jetstream & mechanisms associated with extreme rainfall in South Asia Kieran M R Hunt, Andrew G Turner & Leonard C Shaffrey WMO 6 th International Workshop for Monsoons,
More informationDROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION
DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION Researcher: Saad-ul-Haque Supervisor: Dr. Badar Ghauri Department of RS & GISc Institute of Space Technology
More informationStudy of Hydrometeorology in a Hard Rock Terrain, Kadirischist Belt Area, Anantapur District, Andhra Pradesh
Open Journal of Geology, 2012, 2, 294-300 http://dx.doi.org/10.4236/ojg.2012.24028 Published Online October 2012 (http://www.scirp.org/journal/ojg) Study of Hydrometeorology in a Hard Rock Terrain, Kadirischist
More informationImpact of ATOVS data in a mesoscale assimilationforecast system over the Indian region
Impact of ATOVS data in a mesoscale assimilationforecast system over the Indian region John P. George and Munmun Das Gupta National Centre for Medium Range Weather Forecasting, Department of Science &
More informationGeneral Circulation. Nili Harnik DEES, Lamont-Doherty Earth Observatory
General Circulation Nili Harnik DEES, Lamont-Doherty Earth Observatory nili@ldeo.columbia.edu Latitudinal Radiation Imbalance The annual mean, averaged around latitude circles, of the balance between the
More informationDEPARTMENT OF EARTH & CLIMATE SCIENCES Name SAN FRANCISCO STATE UNIVERSITY Nov 29, ERTH 360 Test #2 200 pts
DEPARTMENT OF EARTH & CLIMATE SCIENCES Name SAN FRANCISCO STATE UNIVERSITY Nov 29, 2018 ERTH 360 Test #2 200 pts Each question is worth 4 points. Indicate your BEST CHOICE for each question on the Scantron
More informationCHAPTER 13 WEATHER ANALYSIS AND FORECASTING MULTIPLE CHOICE QUESTIONS
CHAPTER 13 WEATHER ANALYSIS AND FORECASTING MULTIPLE CHOICE QUESTIONS 1. The atmosphere is a continuous fluid that envelops the globe, so that weather observation, analysis, and forecasting require international
More informationSHORT COMMUNICATION EXPLORING THE RELATIONSHIP BETWEEN THE NORTH ATLANTIC OSCILLATION AND RAINFALL PATTERNS IN BARBADOS
INTERNATIONAL JOURNAL OF CLIMATOLOGY Int. J. Climatol. 6: 89 87 (6) Published online in Wiley InterScience (www.interscience.wiley.com). DOI:./joc. SHORT COMMUNICATION EXPLORING THE RELATIONSHIP BETWEEN
More informationSanjeev Kumar Jha Assistant Professor Earth and Environmental Sciences Indian Institute of Science Education and Research Bhopal
Sanjeev Kumar Jha Assistant Professor Earth and Environmental Sciences Indian Institute of Science Education and Research Bhopal Email: sanjeevj@iiserb.ac.in 1 Outline 1. Motivation FloodNet Project in
More informationAnalysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia.
Analysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia. 1 Hiromitsu Kanno, 2 Hiroyuki Shimono, 3 Takeshi Sakurai, and 4 Taro Yamauchi 1 National Agricultural
More informationMonitoring and evaluation of seasonal snow cover in Kashmir valley using remote sensing, GIS and ancillary data
Monitoring and evaluation of seasonal snow cover in Kashmir valley using remote sensing, GIS and ancillary data HSNegi, N K Thakur, Rajeev Kumar and Manoj Kumar Snow and Avalanche Study Establishment,
More informationA Preliminary Climatology of Extratropical Transitions in the Southwest Indian Ocean
A Preliminary Climatology of Extratropical Transitions in the Southwest Indian Ocean Kyle S. Griffin Department of Atmospheric and Environmental Sciences, University at Albany, State University of New
More informationTHE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS
THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS Brian J. Soden 1 and Christopher S. Velden 2 1) Geophysical Fluid Dynamics Laboratory National Oceanic and Atmospheric Administration
More informationThe Climatology of Clouds using surface observations. S.G. Warren and C.J. Hahn Encyclopedia of Atmospheric Sciences.
The Climatology of Clouds using surface observations S.G. Warren and C.J. Hahn Encyclopedia of Atmospheric Sciences Gill-Ran Jeong Cloud Climatology The time-averaged geographical distribution of cloud
More informationPrediction of Snow Water Equivalent in the Snake River Basin
Hobbs et al. Seasonal Forecasting 1 Jon Hobbs Steve Guimond Nate Snook Meteorology 455 Seasonal Forecasting Prediction of Snow Water Equivalent in the Snake River Basin Abstract Mountainous regions of
More informationApplication and verification of ECMWF products 2011
Application and verification of ECMWF products 2011 National Meteorological Administration 1. Summary of major highlights Medium range weather forecasts are primarily based on the results of ECMWF and
More informationConvective-scale NWP for Singapore
Convective-scale NWP for Singapore Hans Huang and the weather modelling and prediction section MSS, Singapore Dale Barker and the SINGV team Met Office, Exeter, UK ECMWF Symposium on Dynamical Meteorology
More informationClimate Variability in South Asia
Climate Variability in South Asia V. Niranjan, M. Dinesh Kumar, and Nitin Bassi Institute for Resource Analysis and Policy Contents Introduction Rainfall variability in South Asia Temporal variability
More informationINVESTIGATION FOR A POSSIBLE INFLUENCE OF IOANNINA AND METSOVO LAKES (EPIRUS, NW GREECE), ON PRECIPITATION, DURING THE WARM PERIOD OF THE YEAR
Proceedings of the 13 th International Conference of Environmental Science and Technology Athens, Greece, 5-7 September 2013 INVESTIGATION FOR A POSSIBLE INFLUENCE OF IOANNINA AND METSOVO LAKES (EPIRUS,
More informationThe ECMWF Extended range forecasts
The ECMWF Extended range forecasts Laura.Ferranti@ecmwf.int ECMWF, Reading, U.K. Slide 1 TC January 2014 Slide 1 The operational forecasting system l High resolution forecast: twice per day 16 km 91-level,
More informationDeterministic and Probabilistic prediction approaches in Seasonal to Inter-annual climate forecasting
RA 1 EXPERT MEETING ON THE APPLICATION OF CLIMATE FORECASTS FOR AGRICULTURE Banjul, Gambia, 9-13 December 2002 Deterministic and Probabilistic prediction approaches in Seasonal to Inter-annual climate
More informationPresentation Overview. Southwestern Climate: Past, present and future. Global Energy Balance. What is climate?
Southwestern Climate: Past, present and future Mike Crimmins Climate Science Extension Specialist Dept. of Soil, Water, & Env. Science & Arizona Cooperative Extension The University of Arizona Presentation
More informationStability in SeaWinds Quality Control
Ocean and Sea Ice SAF Technical Note Stability in SeaWinds Quality Control Anton Verhoef, Marcos Portabella and Ad Stoffelen Version 1.0 April 2008 DOCUMENTATION CHANGE RECORD Reference: Issue / Revision:
More informationMonitoring of Temporal Variation of Snow Depth Using Remote Sensing in Western Himalaya
Monitoring of Temporal Variation of Snow Depth Using Remote Sensing in Western Himalaya Ruby Nanchahal*, H. S. Gusain**, Darshan Singh Sidhu* and V. D. Mishra** *PTU GZS Campus Bathinda, Punjab, 151 001,
More informationSPI: Standardized Precipitation Index
PRODUCT FACT SHEET: SPI Africa Version 1 (May. 2013) SPI: Standardized Precipitation Index Type Temporal scale Spatial scale Geo. coverage Precipitation Monthly Data dependent Africa (for a range of accumulation
More informationLand Surface: Snow Emanuel Dutra
Land Surface: Snow Emanuel Dutra emanuel.dutra@ecmwf.int Slide 1 Parameterizations training course 2015, Land-surface: Snow ECMWF Outline Snow in the climate system, an overview: Observations; Modeling;
More informationGIS based Estimation of snow depth in mountainous Himalayan region: A case study
GIS based Estimation of snow depth in mountainous Himalayan region: A case study Chander Shekhar *1, S K Dewali *, Snehmani * * Snow and Avalanche Study Establishment, Chandigarh, India 160036, India 1
More informationDenver International Airport MDSS Demonstration Verification Report for the Season
Denver International Airport MDSS Demonstration Verification Report for the 2015-2016 Season Prepared by the University Corporation for Atmospheric Research Research Applications Division (RAL) Seth Linden
More informationModel Output Statistics (MOS)
Model Output Statistics (MOS) Numerical Weather Prediction (NWP) models calculate the future state of the atmosphere at certain points of time (forecasts). The calculation of these forecasts is based on
More informationJoint Meeting of RA II WIGOS Project and RA V TT-SU on 11 October 2018 BMKG Headquarter Jakarta, Indonesia. Mrs. Sinthaly CHANTHANA
Joint Meeting of RA II WIGOS Project and RA V TT-SU on 11 October 2018 BMKG Headquarter Jakarta, Indonesia Mrs. Sinthaly CHANTHANA Lao PDR Background Department of Meteorology and Hydrology ( DMH ) in
More informationApplication and verification of ECMWF products 2013
Application and verification of EMWF products 2013 Hellenic National Meteorological Service (HNMS) Flora Gofa and Theodora Tzeferi 1. Summary of major highlights In order to determine the quality of the
More informationForecasting of thunderstorms in the pre-monsoon season at Delhi
Forecasting of thunderstorms in the pre-monsoon season at Delhi N Ravi 1, U C Mohanty 1, O P Madan 1 and R K Paliwal 2 Meteorol. Appl. 6, 29 38 (1999) 1 Centre for Atmospheric Sciences, Indian Institute
More informationMemorandum. Höfundur: Halldór Björnsson, Nikolai Nawri, Guðrún Elín Jóhannsdóttir and Davíð Egilson.
EBV-007-1 Memorandum Date: 17.12.2015 Title: Estimation of evaporation and precipitation in Iceland Höfundur: Halldór Björnsson, Nikolai Nawri, Guðrún Elín Jóhannsdóttir and Davíð Egilson. Ref: 2015-69
More information