EVALUATION OF STATISTICAL MODELS FOR CLIMATIC CHARACTERIZATION OF GKVK STATION BHOOMIKA RAJ, R. PALB 3153

Size: px
Start display at page:

Download "EVALUATION OF STATISTICAL MODELS FOR CLIMATIC CHARACTERIZATION OF GKVK STATION BHOOMIKA RAJ, R. PALB 3153"

Transcription

1 EVALUATION OF STATISTICAL MODELS FOR CLIMATIC CHARACTERIZATION OF GKVK STATION BHOOMIKA RAJ, R. PALB 3153 DEPARTMENT OF AGRICULTURAL STATISTICS, APPLIED MATHEMATICS AND COMPUTER SCIENCE UNIVERSITY OF AGRICULTURAL SCIENCES GKVK, BENGALURU 2015

2 EVALUATION OF STATISTICAL MODELS FOR CLIMATIC CHARACTERIZATION OF GKVK STATION BHOOMIKA RAJ, R. PALB 3153 Thesis submitted to the University of Agricultural Sciences, Bengaluru In partial fulfilment of the requirements For the award of the degree of Master of Science (Agriculture) In Agricultural Statistics BENGALURU JULY, 2015

3 AFFECTIONATELY DEDICATED TO MY BELOVED PARENTS, MY SISTER, MY FRIENDS & MY TEACHERS

4 DEPARTMENT OF AGRICULTURAL STATISTICS, APPLIED MATHEMATICS AND COMPUTER SCIENCES UNIVERSITY OF AGRICULTURAL SCIENCES BENGALURU CERTIFICATE This is to certify that the thesis titled Evaluation of statistical models for climatic characterization of GKVK station. Submitted by Ms. BHOOMIKA RAJ, R., ID NO. PALB 3153 in partial fulfilment of the requirement for the degree of Master of Science (Agriculture) in Agricultural Statistics to the University of Agricultural Sciences, Bengaluru, is a record of bonafide research work done by her during the period of her study in this University under my guidance and supervision and that no part of the thesis has previously formed the basis for the award of any degree, diploma, associate ship, fellowship or other similar titles. Bengaluru, July, 2015 Mr. K.N. Krishnamurthy (MAJOR ADVISOR) Approved by: Chairman: Mr. K.N. Krishnamurthy Members: 1. Dr. M. Gopinath Rao 2. Dr. M.B. Rajegowda 3. Dr. G. B. Mallikarjuna

5 ACKNOWLEDGEMENT Success of any scientific research is not dependent on an individual effort. It is a harvest of seeds sown by group of like-minded people. It is a great juncture to recall all those faces and spirits in the form of teachers, friends and dear ones. I place on record my deep sense of gratitude with at most sincerity and heartfelt respects to the esteemed chairman of my Advisory Committee Mr. K. N. Krishnamurthy, Associate Professor, Dept. of Agricultural Statistics, GKVK, Bengaluru, for his excellent guidance, constant support, close counsel and valuable suggestions throughout the period of my study. I admit that it has been a great fortune for me to be associated from him during my degree programme. I am very much grateful and thankful to my advisory committee Dr. M. Gopinath Rao, Professor and Head, Dept. of Agricultural Statistics,UAS, GKVK, Bengaluru, Dr. M.B. Raje Gowda, The Registrar, UAS, GKVK, Bengaluru and Dr. G.B. Mallikarjuna, Associate Professor, Dept. of Agricultural Statistics, UAS, GKVK, Bengaluru, for their constant supervision, invaluable guidance and all the facilities extended during the course of this investigation. I gratefully record my indebtedness to the staff members Dr. D.M. Gowda, Professor and University Head, Mr. H..S.S Surendra, Associate Professor, Mr. V. Manjunath, Associate Professor, Mr. R Munirajappa,, Associate Professor, Dr. H. Chandrashekhar, Professor, Dr. S..S..S. Surpur, Professor, Mr. K. B. Murthy, Associate Professor, Dr. S. S Patil., Associate Prfessor for their timely help rendered during the course of my research. On my pesronel note, I humbly place before my parents Rajanna,T and Asha,T, my most sincere gratitude. Their blessings have renewed me every day, all the way on

6 the journey through my life. I am also thankful to my sister Yashaswini Raj, R for her constant support. I am fortunate to have a myriad of friends here. I am thankful for the emotional support from my friends Amrutha,T., Bhavana,D., Chethan,A. J., Nishitha,K., Pavithra, H.K, Priyanka Basavaraddi, Seema Raj, Shashi Kumar,N., Sitara Makam, Shivani Dechamma, Supriya,D., S Swapna Shetty, who encouraged me in each and every step of my post graduation. My special thanks to Aditya K. S., Bhavyashree,S., Golive Prashanthi, Navya, H. S., for their valuable suggestions, moral support and cooperation during the degree programme. It gives me great pleasure to thank all my classmates Deepti K.A., Kavitha D.S., Madhushree C.S., Pavithra R., Satish Y., Shilpashree G.S., Sowmya S Desai, and all juniors in the Department of Agricultural Statistics for their affection and cooperation rendered to me in all aspects. Last but not the least, I would like to thank non-teaching faculty Miss. Jayasudha, Mr. Bayanna and Mr. Umesh, Department of Agricultural Statistics for their kind co-operation and encouragement during my study and research. I owe all my success to my beloved father, mother, sister and teachers for their inspiration and support that saw me through all the trials and tribulations. BENGALURU JULY 2015 Bhoomika Raj,R.

7 EVALUATION OF STATISTICAL MODELS FOR CLIMATIC CHARACTERIZATION OF GKVK STATION BHOOMIKA RAJ, R. ABSTRACT Agro-climatic characters play an important role in deciding the cropping pattern of a region. The distribution of rainfall is one such climatic character essential to plan farm activities in a given region. The present study was conducted to know the climatic characterization of GKVK station. The secondary data of weather parameters over a period of 38 years ( ) and Finger millet (GPU 28) crop yield (qtl/ha) data during kharif season for a period of 16 years ( ) was collected from AICRP on Agro Meteorology and AICSMIP, UAS, GKVK, Bengaluru respectively. Among the weather parameters, amount of maximum daily rainfall (mm) was considered to fit appropriate probability distributions. The probability distributions viz., Normal, Log- normal, Gamma (1P, 2P, 3P), Generalized Extreme Value (GEV), Weibull (1P, 2P, 3P), Gumbel and Pareto were used to evaluate the best fit for maximum daily rainfall (mm). Kolmogorov-Smirnov test for the goodness of fit of the probability distributions showed that for majority of the data sets on rainfall at different study periods, Weibull (3P) distribution was found to be the best fit. However, all the data sets were scale dominated which indicated large variation in the distribution of rainfall. Simple and multiple linear regression analysis showed that none of the weather parameters influenced the finger millet crop yield significantly. The path coefficient analysis indicated that among all the weather parameters, direct effect of PET on rainfall was found to be highly negative, while evaporation had the highest positive direct effect on rainfall. Further, vapor pressure and evaporation had positive and negative indirect contributions through the other weather parameters respectively. Bengaluru July 2015 K.N. Krishnamurthy Major Advisor

8 fpé«pé ದ ದ ಹ ನ ರ ಪ ಯ ಅ ಅ ಶಗಳ ದ ಗಳ ಲ ಪನ sàæ«äpá gáeï Dgï ಶ ಕ -ಹ ನ ತ ಗಳ ಒ ದ ಪ ಶದ ಕ ದ ಯನ ಧ ಸ ವ ಮಹತ ದ ತ ವನ ವ ಸ ತ. ಅ ತಹ ಹ ನ ಲ ಣಗಳ, ಮ ತರ ಯ MAzÀÄ ಪ ಶದ ಕ ಚಟ ವ ಗಳ ಜ ಯ ಅತ ಗತ ತ ವನ ªÀ9» ÀÄತ zé. ಪ ಸ ತ ಅಧ ಯನªÀÅ fpé«pé ದ ದ ಹ ನ ತ ಯಲ ನ ಸ ತ. 16 ವಷ ಗಳ ( ) ಲ ಮ ರ ಋತ ನ ( ಯ 28) ಇಳ ವ ( / ) ಮತ 38 ವಷ ಗಳ ( ) ಅವ ಯ ಹ ನ ಯ ಕಗಳ ಯಕ ಯನ JL¹Dgï ಕ ಹ ನ ಮತ JL¹JªÀiïL, AiÀÄÄJJ ï, fpé«pé, ಗಳ j AzÀ ಕ ಮ ಪ ಯ ತ. ಹ ನ ಯ ಕಗಳ ಗ ಷ ನ ನ ಮ ( ) ಪ ಣವನ ಸ ಕ ಸ ಭವ ಯ ತರ ಗಳ ಳ ಲ ಪ ಗ ಸ. ಇ ಪ ಗ ¹gÀĪÀ ಸ ಭವ ಯ ತರ ಗ¼ÉAzÀgÉ, ನ, - ನ, (1P, 2P, 3P), ರಣ ದ ಎ ಲ (fe«), ಬ (1P, 2P, 3P), ಗ ಬ ಮತ. ಸ ಭವ ಯ ಹ ಗಳ, ಧ ಅಧ ಯನ ಲಗಳ ಮ ದ ಶ ಗ ಪ ಗಳ ಬಹ ಕ ಬ (3P) ತರ ಯ ಅತ ತ ಮ JA ÄzÀÄ - ಪ AiÀÄ ªÀÄÆ PÀ ಕ ಡ ಬ ತ. ಆ ಗ, ಎ UÀÄA ÀÄ ಪ ಣದ ಮ ಹ ಯ ವ ಸಗಳ ಸ ಸ ವ ಪ ವ ರ ತ. ವ ಹ ನ ಯ ಕಗ¼ÀÄ ಇಳ ವ ಗಮ ಹ ಪ ವ ÃgÀĪÀÅ ಲ JA ÄzÀÄ ಸರಳ ಮತ ಹಲ ರ ತ ಕ ಜ ತ ಷ ÄAzÀ w½zàä A vàä. ಗ ಗ ಕ ಷ AiÀÄ ಪ ರ, ಎ ಹ ನ ಯ ಕಗಳ, ಮ ಯ ಇ AiÀÄÄ ಚ ಗಮ ಹ ನ ತ ಕ ರ ಪ ಮªÀನ ಗ ಕರಣªÀÅ ಚ ಗಮ ಹ ಸ ತ ಕ ರ ಪ ಮªÀನ ಸ zàªàå. ಇದಲ, ಆ ಯ ಒತ ಡ ಮತ ಕರಣUÀ¼ÀÄ ಕ ಮ ಇತರ ಹ ನದ ಯ ಕಗಳ ಮ ಲಕ ಸ ತ ಕ ಮತ ನ ತ ಕ ಪ ಡ ಯನ ದ ªÀÅ. ಗಳ ರ ಜ 2015 ಕ ಷ ಮ,. ಎ. ಪ ಮ ಖ ಸಲ ರ

9

10 CONTENTS CHAPTER TITLE PAGE No. I INTRODUCTION 1-7 II REVIEW OF LITERATURE 8-20 III MATERIAL AND METHODS IV RESULTS V DISCUSSION VI SUMMARY VII REFERENCES 80-84

11 LIST OF TABLES TABLE No. TITLE PAGE No Description of various probability distribution functions Summary of Statistics for Maximum daily rainfall (mm) Study period wise probability distributions using goodness of fit test Parameters of the best fit probability distributions for Maximum daily rainfall Summary of Statistics for weather parameters Estimates of the simple linear regression model relating the crop yield with weather parameters Multiple linear regression of weather parameters on the yield of finger millet Observed and predicted yield of finger millet crop over a period of 16 years ( ) Correlation matrix between weather parameters and rainfall Path coefficients showing direct and indirect effects of weather parameters on rainfall Direct and indirect partitioning of correlation between rainfall and weather parameters

12 LIST OF FIGURES FIGURE No. TITLE PAGE No. 1 Agro-climatic zones of Karnataka Conceptual regimes described in parameter space with the shape parameter on the x-axis and the scale parameter on the y-axis Cause and effect relationship Variation in Annual maximum daily rainfall (mm) Variation in Seasonal maximum daily rainfall (mm) Plot of best fitted distribution functions for Maximum daily rainfall (mm). Effect of weather parameters on the finger millet crop yield. Plot of PRESS residuals against predicted finger millet yield. Normal Probability Plot for finger millet yield against weather parameters. Path diagram for the direct and indirect effects of weather parameters on rainfall

13 LIST OF ABBREVIATIONS SL. No. ABBREVATIONS 1 SMW - Standard Meteorological Week 2 K-S - Kolmogorov-Smirnov 3 Gamma (1P) - One parameter Gamma distribution 4 Gamma (2P) - Two parameters Gamma distribution 5 Gamma (3P) - Three parameters Gamma distribution 6 Weibull (1P) - One parameter Weibull distribution 7 Weibull (2P) - Two parameters Weibull distribution 8 Weibull (3P) - Three parameters Weibull distribution 9 GEV- Generalized Extreme Value 10 PET- Potential Evapotranspiration

14 I INTRODUCTION Climate is a measure of average pattern of variation in temperature, humidity, atmospheric pressure, wind speed, precipitation, atmospheric particle count and other meteorological variables in a given region over long periods of time. Climate is different from weather, in that weather only describes the short-term conditions of these variables in a given region. The Intergovernmental Panel on Climate Change (IPCC) (Anon., 2007) defines climate change as, changes in the mean and/or the variability of its properties that persists for an extended period, typically decades or longer. A region's climate is generated by the climate system, which has five components: atmosphere, hydrosphere, cryosphere, lithosphere, and biosphere. Climate is usually defined as the "average weather" or as the statistical description in terms of the mean and variability of relevant quantities over a period ranging from months to thousands or millions of years. The classical period is 30 years, as defined by the World Meteorological Organization (WMO). These quantities are most often surface variables such as temperature, precipitation, and wind. Climate in a wider sense is the state, including a statistical description, of the climate system. The difference between climate and weather is usefully summarized by the popular phrase Climate is what we expect, weather is what we get. Climate change brings about changes in the weather conditions. There is a reason to believe that climate change could affect agricultural productivity, and cause increased health hazards and submergence of lands due to rise in the sea level to name a few. Climate change is the net result of many factors caused by continuous evolution of Planet Earth through many geological eras. However, there is growing concern about manmade developments causing, even if partially or insignificantly, the climate change outcomes. The industrialization that started from the late 17 th century is believed to have accelerated the process of climate change by emissions of Greenhouse Gases (GHGs) to the atmosphere. The observed levels of GHGs have perhaps nearly crossed tolerance levels in the atmosphere so that the survival for many animal and human species is at stake, while developmental needs of human race are contributing to factors like deforestation, urbanization etc., that can speed up the process of climate change. Climatic variability is the major factor influencing the agriculture productivity. Global climate change and its impacts on agriculture have become an important issue. Agriculture production is highly dependent on climate and it is also adversely affected by increasing climatic variability. Climate change is one of the most important global environmental challenges facing humanity with implications for food production, natural ecosystems, freshwater supply, health, etc. Climate classification The Koppen Climate Classification System is the most widely used system for classifying the world's climates. Its categories are based on the annual and monthly Evaluation of Statistical Models for Climatic Characterization of GKVK Station 1

15 averages of temperature and precipitation. The Koppen system recognizes five major climatic types, each type being designated by a capital letter. A- Tropical Moist Climates: all months have average temperatures above 18 Celsius. B - Dry Climates: with deficient precipitation during most of the year. C - Moist Mid-latitude Climates with Mild Winters. D - Moist Mid-Latitude Climates with Cold Winters. E - Polar Climates: with extremely cold winters and summers. Climate of India India is home to an extraordinary variety of climatic regions, ranging from tropical in the south to temperate and alpine in the Himalayan north, where elevated regions receive sustained winter snowfall. The Himalayas and the Thar Desert strongly influence the nation s climate. The Himalayas act as a barrier to the frigid katabatic winds flowing down from Central Asia keeping the bulk of the Indian subcontinent warmer than most locations at similar latitudes. Land areas in the north of the country have a continental climate with severe summer conditions that alternates with cold winters when temperatures plunge to freezing point. In contrast are the coastal regions of the country, where the warmth is unvarying and the rains are frequent. The country is influenced by two seasons of rains, accompanied by seasonal reversal of winds from January to July. During the winters, dry and cold air blowing from the northerly latitudes from a north-easterly direction prevails over the Indian region. Consequent to the intense heat of the summer months, the northern Indian landmass becomes hot and draws moist winds over the oceans causing a reversal of the winds over the region which is called the summer or the South-west monsoon. This is the most important feature controlling the Indian climate because about 75 per cent of the annual rainfall is received during a short span of four months (June to September). Variability in the onset, withdrawal and quantum of rainfall during the monsoon season has profound impacts on water resources, power generation, agriculture, economics and ecosystems in the country. The variation in climate is perhaps greater than any other area of similar size in the world. There is a large variation in the amounts of rainfall received at different locations. The rainfall pattern roughly reflects the different climate regimes of the country, which vary from humid in the north east (about 180 days rainfall in a year), to arid in Rajasthan (20 days rainfall in a year). So significant is the monsoon season to the Indian climate, that the remaining season are often referred relative to the monsoon. The rainfall over India has large spatial as well as temporal variability. A homogeneous data series has been constructed for the period based on the uniform network of 1476 stations and analyzed the variability and trends of rainfall. Normal monsoon rainfall more than 150cm is being observed over most parts of northeast India, Konkan and Goa. Normal monsoon rainfall is more than 400cm over major parts of Meghalaya. Annual rainfall is more than 200 cm over these regions. 2 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

16 For the country as a whole, mean monthly rainfall during July (286.5 mm) is highest and contributes about 24.2 per cent of annual rainfall ( mm). The mean rainfall during August is slightly lower and contributes about 21.2 per cent of annual rainfall. June and September rainfall are almost similar and contribute 13.8 per cent and 14.2 per cent of annual rainfall, respectively. The mean South-west monsoon (June, July, August & September) rainfall (877.2 mm) contributes 74.2 per cent of annual rainfall ( mm). Contribution of pre-monsoon (March, April & May) rainfall and postmonsoon (October, November & December) rainfall in annual rainfall is mostly the same (11 per cent). Coefficient of variation is higher during the months of November, December, January and February. India is characterized by strong temperature variations in different seasons ranging from mean temperature of about 10 C in winter to about 32 C in summer season. There are four major types of climate found in India. They are presented as below: 1. The cold weather season (December to February) 2. The hot weather season (March to May) 3. The rainy season (June to September) 4. The Season of Retreating South-west Monsoon (October-November). 1) Cold Weather Season January and February are the coldest months of this season. The temperature is between 10 C to 15 C in Northern India and about 25 C in Southern India. 2) Hot Weather Season The North Indian region experiences a well defined hot weather season during the month of April and May. The temperature starts rising by the middle of March and by mid May, mercury touches 41 to 42 C. Temperature even exceeds 45 C in areas of central and north-west India. 3) The Rainy Season The inflow of South-westerly monsoon in India brings the season of rain in India. The monsoon may burst in the first week of June or even earlier in the coastal areas, while in the interior it may be delayed to the first week of July. With the onset of rains, temperature starts falling. The Indian sub-continent receives the bulk of its rainfall during the south-west monsoon period. 4) The Season of Retreating Monsoon The South-west monsoon begins to retreat from northern India by the second week of September. The pattern of retreat also shows interesting regional variations. The Evaluation of Statistical Models for Climatic Characterization of GKVK Station 3

17 weather during the season is characterized by high day temperatures, but nights are pleasant with the mean minimum temperature going down to 20 C or even lower. Climate of Karnataka Climate of Karnataka presents an exceptional diversity. While the hilly and plateau regions demonstrate a different climatic behavior, the plain presents comparatively a warmer atmosphere. The varying geographic and physiographic conditions of the State is responsible for the climatic variation in the State from arid to semi-arid in the plateau region, sub humid to humid tropical in the Ghats and humid tropical monsoon type in the west coast plains. For meteorological purposes, the State has been divided into three subdivisions: Coastal Karnataka, which includes, Dakshina Kannada and Uttara Kannada districts. North Interior Karnataka, which includes, Belagavi, Bidar, Vijapura, Dharwad, Kalaburgi and Raichur districts. South Interior Karnataka, which includes the remaining districts of Bengaluru Rural, Begaluru urban, Ballari, Chikkamamagaluru, Chitradurga, Kodagu, Hassan, Kolar, Mysuru, Mandya, Shivamogga and Tumakuru districts. As per Koppen's classification, the State witnesses three climatic types. The tropical monsoon covers the entire coastal belt and the adjoining areas. The southern half of the State, outside the coastal belt experiences hot, seasonally dry tropical savanna climate. The remaining regions of the southern half of the State experiences hot, semiarid, tropical steppe type of climate. According to the Thronthwaite's classification, the coastal and Malnad regions are pre-humid i.e. those having moisture index of 100 per cent and above. The interior regions are semi-arid (moisture index of minus 66.7 per cent to minus 33.3 per cent). Moist sub-humid and dry-humid zones (moisture index of minus 33.3 to plus 20 per cent) are the transition zones covering the region between the Malnad and the coast. The arid zone in the State is confined to the east of Bellary district, most of Raichur district, east of Chitradurga district and the adjoining Pavagada taluk of Tumakuru district with small areas in Vijapura and the adjoining north-eastern Belagavi district. Very dry areas with moisture indices less than minus 60 per cent occur in Chitradurga, Bellary, Raichur and Vijapura districts, west and south of Kalaburgi district and north Tumakuru district. Semi- arid regions with moisture indices of less than 50 per cent occur in Bidar district in the north, Bengaluru district and adjoining areas of Tumakuru and Mysuru districts in the south. The sub- humid zone in the State exists as a narrow belt east of Western Ghats from Belagavi in the north to the west of Heggadadevanakote taluk of Mysuru district in the south. Adjoining to this in the west is a narrow strip of humid and a wider strip of pre-humid zones. About 77 per cent of the total geographical area of the State, covering interior Karnataka is arid or semi-arid with the State contributing 15 per cent of the total semi-arid or 3 per cent of the total arid areas of the country. 4 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

18 Annual Rainfall Coastal Karnataka: 3456mm North Interior Karnataka: 731mm South Interior Karnataka: 1126mm Summer The average weather of Karnataka is dry and warm over different regions and summers starting from the month of April which last till the month of May. These months are the hottest months in Karnataka; somewhere the humidity percentage is comparatively low, but as the month of June starts to pervade humidity in the air could make you uncomfortable as the monsoon is reaching the state soon. The average temperature remains around 34 C with 75 per cent humidity. Monsoon Monsoon season starts from June and lasts till September, as prominent downfalls in temperature are noted, but at this time the percentage of humidity gets a little higher in the atmosphere. The Tropical Monsoon climate covers the entire coastal belt and adjoining areas. This area experiences heavy rainfall 3456 mm annually, while the North interior Karnataka and its adjoining areas; Vijapura, Bagalkote, Belagavi, Haveri, Gadag, Dharwad, Kalaburgi, Ballari, Koppal and Raichur districts experience mediocre rainfall of 731mm per annum. On the other side, the South interior Karnataka receives a blissful shower of monsoon annually. Winter Winters start from January and last till February. During winters, the climate of Karnataka remains pleasant. Every region of Karnataka experiences a cold winter - wave and average temperature falls between 32 C to 20 C or even below that. Karnataka receives a delightful shower during months of October and November. Agro-Climatic Zones 1. North Eastern Transition Zone 2. North Eastern Dry Zone 3. Northern Dry Zone 4. Central Dry Zone 5. Eastern Dry Zone 6. Southern Dry Zone 7. Southern Transition Zone 8. Northern Transition Zone Evaluation of Statistical Models for Climatic Characterization of GKVK Station 5

19 9. Hill Zone 10. Coastal Zone Need for the study Agro climatic characters play an important role in deciding the ways of life, crop of the region and other developmental plans. The study of climate at different scales and its variations with the topography and the nature of soils can play a fundamental role in energy and/or environmental conscious land use planning. Rainfall is one of the most important natural input resources to crop production. Its occurrence and distribution is erratic, and shows temporal and spatial variations in nature. Hence, there is a need to study the distribution both on long term (annually or seasonally) as well as short term (monthly or weekly). Most of the hydrological events occurring as natural phenomena are observed only once. One of the important problems in hydrology deals with the interpreting past records of hydrological event in terms of future probabilities of occurrence. The average rainfall of the region is generally considered as the basis for deciding irrigation management and cropping pattern. Long term analysis is vital for firm planning and execution of crop cultivation. Rainfall is considered as principle source of water. The success or failure of crops particularly under rainfed conditions is closely linked with the rainfall patterns. The important characteristics of rainfall influencing production from rainfed farming are the duration of wet spells and durations of intervening dry spells. Rainfall during the monsoon is not uniform. The annual and seasonal analysis of rainfall will give general idea about the rainfall pattern of the region, whereas the monthly analysis of rainfall will be of much use as far as agricultural planning is concerned. The spatial and temporal variability of rainfall and its uneven and inadequate distribution is based on other weather parameters such as, maximum and minimum temperatures, relative humidity, evapo-transpiration, sunshine hours, etc. This determines the success or failure of crops especially in drought prone areas. Knowledge of the distribution of rainfall is essential for successful farming. It is also important to know the chances of occurrence during the critical stages of the crops for deciding the sowing date, cropping pattern and planning for protective irrigation and intercultural operations. With this view in mind, the study was conducted for GKVK station with the following specific objectives: 1. To fit appropriate Statistical distributions for rainfall. 2. To study the impact of weather parameters on the crop yield. 3. To study direct and indirect effects of weather parameters on rainfall. 6 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

20 Fig.1: Agro-climatic zones of Karnataka Evaluation of Statistical Models for Climatic Characterization of GKVK Station 7

21 II REVIEW OF LITERATURE This chapter deals with the reviews of related past studies. Keeping in view of the objectives of the study, the reviews have been presented under the following headings: 2.1 To fit appropriate Statistical distributions for rainfall. 2.2 To study the impact of weather parameters on the crop yield. 2.3.To study direct and indirect effects of weather parameters on rainfall. 2.1 To fit appropriate Statistical distributions for rainfall. Jinfan et al. (1998) compared the general two-parameter Weibull distribution with other widely used distributions for describing the distribution of daily precipitation event sizes at 99 sites from the U.S. Pacific Northwest. They found that a one-parameter Weibull model of a wet day precipitation amount based on the Weibull distribution provided a better fit to historical daily precipitation data for eastern U.S. sites than other one-parameter models. The fit of the single-parameter Weibull to the 99 stations included in this study was significantly better than other single-parameter models tested, and performed as well as the widely endorsed, more cumbersome, two-parameter gamma model. Both the one and two parameter Weibull distributions were shown to have L- moments that were consistent with historical precipitation data, while the ratio of L-skew and L-variance in the gamma model was inconsistent with the historical record by this measure. In addition, it was found that the two-parameter gamma distribution was better fit using the method of moment estimators than maximum likelihood estimates. Gregory et al. (2007) demonstrated the feasibility of fitting cell-by-cell probability distributions to grids of monthly interpolated, Continent-wise data. They applied these grids to improve satellite-remote sensing of drought and interpretation of probabilistic climate outlook for forecasts. They showed that the gamma distribution was well suited for these applications because it was capable of representing a variety of distribution shapes. This study tested the goodness of fit using the Kolmogorov Smirnov (KS) test, and compared these results against another distribution commonly used in rainfall events, the Weibull distribution. The gamma distribution was suitable for roughly 98 per cent of the locations over all months. The techniques and results presented in this study provided a foundation for the use of the gamma distribution to generate drivers for various rain related models. These models were used as decision support tools for the management of water and agricultural resources as well as food reserves by providing decision makers with ways to evaluate the likelihood of various rainfall accumulations and assess different scenarios in Africa. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 7

22 Suhaila and Jemain (2007) have tested Gamma, Weibull, Kappa and Mixed exponential distributions to fit the daily amount of rainfall in Peninsular Malaysia. They estimated the parameters for each distribution using maximum likelihood method. The selected model was chosen based on the minimum error produced by seven goodness of fit (GOF) tests, namely, the Median of Absolute Difference (MAD) between the empirical and hypothesized distributions, the traditional Empirical Distribution Function (EDF) statistics which included Kolmogorov- Smirnov, Anderson Darling statistic and Cramer-Von-Mises statistic and the new method of EDF statistics based on likelihood statistics. Based on this goodness of fit tests, the Mixed Exponential was found to be the most appropriate distribution for describing the daily amount of rainfall in Peninsular Malaysia. Lars and Vogel (2008) studied the distribution of wet-day daily rainfall and identified that the 2-parameter Gamma (2P) distribution as the most likely candidate distribution based on traditional goodness of fit tests. They used probability plot correlation coefficient test statistics and L-moment diagrams to examine the complete series and wet-day series of daily precipitation records at 237 U.S. stations. The analysis indicated that the Pearson Type-III (P3) distribution fits the full record of daily precipitation data remarkably well, while the Kappa (KAP) distribution best described the observed distribution of wet-day daily rainfall. It was also shown that the Gamma (2P) distribution performed poorly in comparison to either the P3 or KAP distributions. Deka et al. (2009) determined the best fitting distribution to describe the annual series of maximum daily rainfall data for the period 1966 to 2007 of nine distantly located stations in North East India. Five extreme value distributions viz., generalized extreme value distribution; generalized logistic distribution, generalized Pareto distribution, lognormal distribution and Pearson distribution were fitted for this purpose using the method of L-moment and LQ-moment. The performance of the distributions was evaluated by using three goodness of fit tests namely, relative root mean square error, relative mean absolute error and probability plot correlation coefficient. Further, L- moment ratio diagram was also used to confirm the goodness of fit for the above five distributions. Finally, goodness of fit test results was compared and generalized logistic distribution was empirically proved to be the most appropriate distribution for describing the annual maximum rainfall series for the majority of the stations in North East India. Mohita and Singh (2010) studied the annual daily maximum rainfall of Pantnagar and showed that it received rainfall ranging from 49.32mm (minimum) to mm (maximum) indicating a very large range of fluctuation during the period of study. It was observed that the best probability distributions obtained for the maximum daily rainfall for different data sets are different. The lognormal and gamma distribution were found as the best fit probability distribution for the annual and monsoon season period of study, respectively. Generalized extreme value distribution was observed in most of the weekly period as best fit probability distribution. The best fit probability distribution of monthly 8 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

23 data was found to be different for each month. The scientific results clearly established that the analytical procedure devised and tested in this study may be suitably applied for the identification of the best fit probability distribution of weather parameters. Manikandan et al. (2011) analyzed daily rainfall data of 37 years and annual one day maximum rainfall was sorted to estimate the probable one day maximum rainfall for different return periods by using the probability distribution function at Tamil Nadu Agricultural University (TNAU) Campus. The mean value of annual one day maximum rainfall was found to be mm with standard deviation, coefficient of variation and skewness of 28.56, 0.37 and 0.71 respectively. October month received the highest amount of one day maximum rainfall and the highest number of one day maximum rainy days. Five probability distributions such as Normal, Log Normal, Gumbel, Pearson Type- III and Log Pearson Type-III distributions were tested to determine the best fit probability distribution that described the annual one day maximum rainfall by comparing with the Chi-square value. The results revealed that the log-normal distribution was the best fit probability distribution for annual one day maximum rainfall. Based on the best fit probability distribution, the minimum rainfall of mm in a day can be expected to occur with 99 per cent probability and one year return period and maximum of mm rainfall can be received with one per cent probability and 100 year return period. The results from the study could be used as a rough guide by design engineers and hydrologists for appropriate planning and design of small soil and water conservation structures, irrigation and drainage systems. Bhim Singh et al. (2012) analyzed the daily rainfall data of 39 years ( ) to determine the annual one day maximum rainfall of the Jhalarapatan area of Rajasthan, India. The observed values were estimated by Weibull's plotting position and expected values were estimated by four well known probability distribution functions viz., normal, log-normal, log-pearson type-iii and Gumbel. The expected values were compared with the observed values and goodness of fit was determined by Chi-square (χ 2 ) test. The results showed that the log-pearson type-iii distribution was the best fit probability distribution to forecast annual one day maximum rainfall for different return periods. Based on the best fit probability distribution, the minimum rainfall of mm in a day can be expected to occur with 99 per cent probability and one year return period and maximum of mm rainfall can be received with one per cent probability and 100 year return period. And it was observed that July month received the highest amount of one day maximum rainfall (46 per cent) followed by August (38 per cent) and September (13 per cent). Oseni and Femi (2012) studied several types of statistical distributions to describe the rainfall distribution in the Ibadan metropolis over a period of 30 years. The exponential, gamma, normal and Poisson distributions were compared to identify the optimal model for a daily rainfall amount based on data recorded at rain gauge station at the Forestry Research Institute of Nigeria, Jericho, and Ibadan (FRIN). The models were Evaluation of Statistical Models for Climatic Characterization of GKVK Station 9

24 evaluated based on chi- square and Kolmogorov-Smirnov tests. Overall, this study showed that the exponential distribution was the best model followed by normal and Poisson model that had the same estimated rainfall amount for describing the daily rainfall in Ibadan metropolis. Lala et al. (2013) analyzed daily rainfall data for 28 years ( ) of Central Meghalaya, Nongstoin station for estimating maximum daily rainfall. The annual maximum daily rainfall data was fitted to five different probability distribution functions i.e. Normal, Log-normal, Pearson Type-III, Log Pearson Type-III and Gumbel Type-I extreme. The probable rainfall value for different return periods was estimated. These estimated values were compared with the values obtained by Weibull s Method. The analysis indicated that, the Gumbel distribution gave the closest fit to the observed data. Hence, Gumbel distribution may be used to predict maximum rainfall, which will be a great importance for economic planning and design of small and medium hydraulic structures. Majid et al.(2013) determined the best probability distribution of the rainfall data in the Schoeckelbach basin, which is situated at the Northern Graz in Austria using L- moments method (parameter estimation technique). They used Chi- square, Kolmogorov- Smirnov and Root mean square error (RMSE) tests to test the goodness of fit. They compared between four commonly used rainfall frequency distribution, such as Generalized Extreme Value (GEV), Gumbel, Log- Pearson type II I(LP III) and 3 parameter Log normal (LN III) and showed that Gumbel was the best probability distribution for Schoeckelbach basin. The best probability distribution was used to determine the Intensity- Duration- Frequency (IDF) relation for the Schoeckelbach basin. Mohsen et al. (2013) investigated the spatial and temporal characteristics of rainfall in the United Arab Emirates (UAE). The region is divided into four climate zones (East Coast, Mountains, Gravel Plains and Desert Foreland) of distinguished rainfall distribution. The rainfall patterns, rainfall probability of occurrences, rainfall intensityduration-frequency (IDF) relationship, probable maximum precipitation (PMP) and drought scenarios were investigated. Daily rainfall data from a network of stations across the UAE were used. Standard statistical techniques were applied for data analyses. The Gumbel, log Pearson, generalized extreme value, log normal, Wake by and Weibull probability distributions were tested to fit extreme rainfalls. Both Gumbel and Weibull distributions were found adequate. Measures of dispersion and symmetry of rainfall patterns were found relatively high. The estimated PMP values were found highest in the East Coast region and lowest in the Gravel Plains region. The estimated drought severity index showed that the regions have similar trends of drought patterns over the years. The study is useful for sustainable water resources planning and management in the region. Taofik et al. (2013) carried out a study to evaluate the utility of Weibull probability distribution in the estimation of the distribution of rainfall data of South- 10 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

25 Western Nigeria. The study was executed using the rainfall data of 4 stations (Ibadan, Ikenne, Ilora and Orin Ekiti) of the Institute of Agricultural Research and Training, Ibadan. The data obtained from the database were subjected to Weibull distribution analysis using Peak Fit program and the parameters were estimated using maximum likelihood estimation. The results showed that the parameter estimates ranged between for the Ibadan rainfall and for the Ilora rainfall while ranged from for rainfall at Ibadan to for rainfall at Ikenne. The statistical properties of this distribution were that the estimated mode of the rainfall was zero because α and < 1 and the mean ranged between (Ibadan) and (Ikenne) while the estimated mean for Ilora was and Orin Ekiti was Both the estimated and actual means followed the same trend, indicating the plausibility of the distribution function. Deka (2014) demonstrated the application of the method of Trimmed L-moment (TL-moment) to determine the best fitting distribution for annual series of maximum daily rainfall data of nine distantly located stations of North- East India. For this purpose, three extreme value distributions viz. Generalized Extreme Value distribution (GEV), Generalized Logistic distribution (GLD), Generalized Pareto distribution (GPD) were considered and the parameters for each of the distributions were estimated using the method of TL-moment. The performance of the distributions was evaluated by using three goodness of fit tests, namely, relative root mean square error, relative mean absolute error and probability plot correlation coefficient. Goodness of fit test results were compared and Generalized Extreme Value distribution was empirically proved to be the most appropriate distribution for describing the annual maximum rainfall data for all the stations considered when the parameters were estimated by using TL-moment method. 2.2 To study the impact of weather parameters on the crop yield. Jain et al. (1980) have made an attempt to study the forecasting of rice yield on climatic variables at different stages of crop growth in Raipur District using the yield data of 25 years and weekly weather variables viz., maximum temperature, relative humidity, total rainfall and number of rainy days. Based on the principal component analysis, weather variables were obtained and used as independent variables in the regression equation. Among the different models, two models were found suitable. In the first, weighted averages of weekly weather variables and their interactions using the power of the week as weights were used. The respective correlation coefficients with yield in place of week number were taken in the second model. The study revealed that forecasting of rice yield is possible by weekly climatic variables, 2 ½ months after sowing for crop of 5 months duration. Using the second model, the study showed that the weighted weather index/ principal component could be used in place of weather variables. Kulkarni et al. (2004) studied the crop yield weather relations using a nonparametric approach based on the concept of two-way contingency tables that had been proposed to account for the associations of the combined effect of rainfall pattern on crop Evaluation of Statistical Models for Climatic Characterization of GKVK Station 11

26 yield. The approach identified both the different rainfall patterns and the levels of crop yields that were likely to occur during the crop season and then measured their associations. The approach was applied for studying the weather (rainfall) effect on the three crops viz., Jowar, Groundnut and Cotton by using 40 years crop yield-rainfall data ( ) of Anantapur district of Andhra Pradesh. The analysis provided useful information about the yield response range corresponding to the different rainfall patterns that were likely to occur during the crop growth season. Das et al. (2007) carried out a preliminary study to test a crop simulator model, ORYZA 2000 in the baro seasons during and in West Bengal. The model predicts the drastic changes in the yield by -10 per cent, per cent and 72.1 per cent with change in temperature, +1 0 C, +2 0 C and +3 0 C respectively. They also revealed that even with warmer climate up to 1 0 C, production may increase by 10 per cent in an atmosphere with doubled CO 2 concentration. Chunqiang et al. (2008) reported that long term changes of reference evapotranspiration and crop water requirement can have great effect on agricultural production, hydrological cycle as well as water resources management and analyzed the effect of climate change and variability on reference evapo-transpiration and crop water requirements. The results showed that reference evapo-transpiration and water requirement had decreased with time. Geethalakshmi et al. (2008) used Info Crop model to evaluate the impact of climate change on crop yield for baseline data and daily weather data created for 2020, 2050 and The CO 2 concentration used in the model for the baseline, 2020, 2050 and 2080 runs were 376, and 682 ppm, respectively as projected in the IPCC (2001) report. Simulated rice production indicated that, generally rice production showed a declining phase in Tamil Nadu by 8.7, 23.6 and 42.2 per cent from the production level of the year 2000 in 2020, 2050 and 2080 respectively. Guruswamy et al. (2008) have assessed the mean annual and seasonal rainfall behaviour and also have identified the influence of weather parameters on the coconut productivity. South West Monsoon and winter rainfall seasons had negative correlation with coconut productivity, whereas summer and North East Monsoons had positive correlation with coconut productivity. The percentage of barren nut production in coconut had positive correlation with summer rainfall and negative correlation with winter, south west and north east monsoons. The weather variables namely maximum temperature, minimum temperature, relative humidity, evaporation and rainfall had positive correlation with coconut productivity. But there was no significant influence on the productivity of the palms. Jadhav et al. (2008) conducted an experiment for 5 years on sorghum by using four different sowing windows to study the relationship between weather parameters and yield 12 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

27 in Rabi season. Rabi sorghum sown at MW 39 (24-30 Sept) produced a maximum grain yield (850 kg/ha), fodder yield (29.06q/ha) and total monetary returns (Rs.13781/ha). However, the crop sown at MW 40 was also produced the on par yields with crop sown at MW 39 and of high degree of sustainability. Weather parameter showed a significant positive association with grain yield. The stepwise multiple regressions were obtained for different finishes wise weather parameters with a yield of Rabi sorghum, which explained 78 per cent variation in the yield. Kersebaum et al. (2008) studied the impact of projected climate change for winter wheat production which was simulated for 9 sites across Germany using the dynamic agro ecosystem model HERMES and downscaled climate change scenarios of the GCM ECHAM5 output for SRESM emission scenario A1B until Yield reductions between 2 and 11 per cent were estimated for 8 sites for the period Decreasing summer precipitation led to a reduced efficiency of nitrogen fertilizers. Khadtare et al. (2008) conducted a study for 5 years on sunflower by using 4 different sowing windows to know the relationship between weather parameters and yield in Rabi season. The minimum temperature, relative humidity and bright sunshine had a highly significant positive association with grain yield at all the phonological phases. Significant negative association with grain yield by pan evaporation indicated that, at early growth stages Rabi sunflower did not favour the moisture stress condition. A significant positive association with grain yield at all stages of growth by minimizing temperature indicates Rabi sunflower responds well to the low temperature condition throughout the growth period. The weather parameter influences their contribution and performance of Rabi sunflower crop sown at different dates were assessed and model on combined effect was developed using stepwise multiple regression for predicting grain yield. Chaudhari et al. (2009) developed the regional models to know the yield response to temperature (minimum, maximum and its diurnal range) and precipitation for meteorological (met) sub-divisions of India and to study the impact of future climate change on major food crops viz. wheat, rice, potato and rapeseed-mustard. The area weighted averages of district-wise crop yield data were computed at met sub-division level for for 9 major wheat (Triticum aestivum L.) producing met subdivisions, 16 major rice (Oryza sativa L.) producing met subdivisions, 6 major potato (Solanum tuberosum L.) growing sub-divisions of and 8 major rapeseed-mustard (Brassica spp.) growing subdivisions. Fortnightly correlation weighted weather parameters like minimum and maximum temperature and precipitation for the respective met sub-division and periods of the crop season were used to develop the empirical relationships. The study showed that, there was a clear negative response of yields to increase minimum temperatures for all three Rabi crops such as wheat, potato and rapeseed-mustard while the mixed response was observed for kharif rice. The reduction impact was high for Rabi crop as compared to kharif rice. The crop yield also showed Evaluation of Statistical Models for Climatic Characterization of GKVK Station 13

28 negative response to increase in maximum temperature. The crops like potato and rapeseed-mustard showed positive response to increased maximum temperature. This might be due to its strong positive correlations with diurnal temperature range (DTR). The mixed impact (increase and decrease both) was observed on rice yield for the increased precipitation. The estimated impacts of diurnal temperature range (DTR) changes on yields were generally less for wheat and rice crops while more for potato and rapeseed-mustard crops with the unit increase in DTR. Rojalin et al. (2009) assessed the impact of temperature and CO 2 on the productivity of the four major cereal food crops (wheat, rice, maize and pearl millet) taking New Delhi as study area. This study also aimed at identifying the most sensitive growth phase of these crops to the rising temperature. A soil-plant growth simulator CropSyst model was used for this purpose. Their study concluded that Rabi (winter) crops are (wheat) more sensitive to high temperature stress than the kharif crops. The positive effect of CO 2 was more in C4 crops (maize) as compared to C3 crops. The present study also indicated that adverse effect of rising temperature could be mitigated through the doubled CO 2 up to 10 0 C in pearl millet, up to 20 0 C in rice and wheat, up to 30 0 C in maize. However, the rise in temperature beyond that temperature could cause significant reduction in yield of all cereal crops under the present management environment. Hence, further research is required to find out the appropriate management practices to sustain the present level yield in a changed climate condition. Ratna and Amrender (2011) studied the application of Neural Networks (NNs) for crop yields (rice, wheat and sugar cane), by forecasting using Multilayer Perceptron (MLP) architecture with different learning algorithm has been attempted. For development of neural network based forecast models, yields of the crop at district level (Uttar Pradesh state, India) was considered as an output variable and indices of weather variables viz. maximum and minimum temperatures, rainfall and morning relative humidity were considered as input variables. Forecasts based on MLP architecture using a Conjugate gradient descent algorithm for learning have been found to be close to the observed ones in most of the cases. The findings of the study substantiate that neural networks possess potential for prediction purposes. Parekh and Suryanarayana (2012) carried out a study to determine the predominance of various meteorological data on the yield of wheat. The meteorological data for the Rabi season are collected and correlated with yield of wheat in Vallabh Vidyanagar for the period using neural network fitting. Then, the model was re-trained until the best coefficient of correlation was obtained and this corresponding model was considered as the best model and this was further used to validate the dataset from 2000 to This whole procedure was repeated for three different Alternatives. In Alternative 1, only maximum and minimum temperatures were correlated with yield data. In Alternative 2, maximum and minimum temperature and relative humidity were correlated with yield data. In Alternative 3, Maximum and minimum temperatures, 14 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

29 sunshine hours and relative humidity were correlated with yield. The correlation between maximum and minimum temperatures & relative humidity and yield, the co-efficient of correlation for training, i.e validations, i.e and testing, i.e were best for 70 per cent, whereas in 30 per cent dataset, R came out to be 0.62, which is good. They evidently concluded from the study that considering three variables, the correlation was achieved as the best. This revealed that the yield of a crop was very much depended on maximum and minimum temperatures and relative humidity. Govinda (2013) examined the effect of precipitation and temperature variation on the yield of major cereals (e.g., rice, wheat, maize, millet and barley) in Dadeldhura district of Nepal. He used the time series information of yield and seasonal meteorological data (e.g., precipitation and temperature) to assess the effect using the exponential growth regression methods. The findings of this study confirmed that maximum temperature and minimum precipitation have had adverse effects on the yield of rice and maize. Temperature is statistically significant for the yield of all cereals whereas precipitation is statistically significant for the yield of maize and barley. Finally, precipitation and temperature have a statistically significant effect on rice and maize yield. The low rainfall and high temperature during the crops growing seasons had severely affected the yield of major cereals in Dadeldhura district. Elbariki et al. (2014) carried out a study to examine the trends of important weather parameters and their effect on the production of cereal crops in Solan district of Himachal Pradesh. The descriptive statistics and regression analysis revealed no significant trend in all the selected weather parameters in annual, seasonal and monthly basis for the period ( ). In regression analysis cubic and quadratic functions (non significant) were found to be the best fit. Sen s estimation analysis indicated increasing trend of maximum temperature by 2.95 o C, minimum temperature by 0.50 o C and decreasing trend of rainfall, relative humidity and sunshine hours. Analysis of relationship of crop yields to time in regression analysis did not show any significant trend. However, Sen s estimate showed increasing trends for selected crops. There was no significant correlation of individual weather parameters with crop yields both in case of multiple regression analysis some effects were observed, indicating thereby that crop yields are influenced by combinations of weather parameters. 2.3 To study direct and indirect effects of weather parameters on rainfall. Kaushik et al.(2003) evaluated thirty-seven genotypes of American cotton at Agricultural Research Station, Sriganganagar during kharif Randomized block design was followed with two replications. A significant positive association was observed for a number of bolls per plant, sympodial branches per plant, plant height and monopodium branches per plant with seed cotton yield per plant. Besides, these components had positive inter-relationships among themselves. Path-coefficient analysis at genotypic level revealed that sympodia per plant, monopodia per plant and boll weight Evaluation of Statistical Models for Climatic Characterization of GKVK Station 15

30 (g) had positive direct effect on seed cotton yield. Selection based on these characters may contribute considerably to the improvement in seed yield. Sharma et al. (2004) studied the effect of meteorological factors on the population build up of green leaf hopper Nephotettix virescens Dist (Cicadellidae, Hemiptera), plant hoppers Cofana spectra Dist (Delphacidae, Hemiptera) and C. Yasumatsui Young (Kolla mimica, Hemiptera) and rice gundhi bug Leptocoriza acuta Thunberg (Alydidae, Hemiptera) in rice growing season (July to November) through light trap collection during ten years ( ). Maximum populations of Nephotettix virescens Dist (Cicadellidae, Hemiptera) and C. Yasumatsui Young (Kolla mimica, Hemiptera) were recorded in the third week of October during all the years. Cofana spectra Dist had maintained peak activity in respect of population in the last week of September and third week of October. Leptocoriza acuta Thunberg (Alydidae, Hemiptera) had a maximum population in second and third weeks of October during the aforesaid period. No meteorological factors had a significant effect on the population build up of Nephotettix virescens Dist, Cofana spectra Dist and C. yasumatsui Young in the month of October. In the case of Leptocoriza acuta Thunberg, no other factor, but rainfall had a positive correlation of order with population build up in the fourth week of September. Anwar et al. (2009) conducted a study at Wheat Research Institute, Faisalabad during where significant genotypic differences existed among the genotypes. Correlation coefficients were computed for grain yield plant-1, tillers plant-1, spikelet s spike-1, 1000 grain weight, spike length, days to heading, days to maturity and plant height from the F1 crosses developed from four lines and three testers including their parents. The results showed that grain yield plant-1 was positively and significantly correlated with number of tillers plant-1 and days to maturity at genotypic level but non significantly correlated at phenotypic level. Days to maturity had positive genotypic correlation with grain yield plant-1, number of tillers plant-1 and 1000-grain weight. Days to maturity and tillers plant-1 had a positive direct effect on grain yield plant-1 also. Therefore, more days to maturity and more tillers plant-1 would be important selection criteria for improved grain yield plant-1 in the breeding material studied. Ramazan (2009) conducted an in and growing seasons with 7 different durum wheat genotypes under West Anatolia conditions. The correlations among plant height, grain number per spike, grain weight per spike, 1000 grain weight, test weight and grain yield as well as direct and indirect effects of those traits on the grain yield were investigated using Path analysis. Grain number per spike, 1000-grain weight, plant height and test weight had significant direct effects on grain yield. It was concluded that these characteristics could be important selection criteria in durum wheat breeding studies. Er et al. (2009) conducted a study to determine the relationships between climatic parameters and sugar beet produced in Konya (Cumra) plain. There were statistically 16 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

31 significant relationships between the monthly mean sugar content of sugar beet and the climatic factors; rainfall, temperature, cloudiness mean and wind velocity (r = 0.630*, r = **, r = **, r = * in 1999; r = 0.528, r =-0.980**, r = 0.673, r = in 2000, respectively). According to path analysis results, temperature, cloudiness mean and rainfall had more direct and indirect effects than other climatic factors on sugar content of sugar beet. Javare Gowda et al. (2009) installed a light trap in Teak plantation to study the effect of various weather parameters. The path analysis of five associated weather parameters indicated that, on the outbreak of moth rainfall had a positive direct effect on the emergence of male and female moth during and ( r = ; r = and r = ; r = , respectively). Indirectly relative humidity in the evening through rainfall influenced the outbreak of male moth ( r= ) and through maximum temperature it influenced the emergence of the female moth ( r = ) during Similarly, relative humidity in the morning through the relative humidity in the evening influenced the outbreak of male and female moth ( r = ; r= , respectively) during Diogo et al.(2013) conducted a study to evaluate the relationships between morphological characteristics and fresh matter yield of cactus pear and identify their direct and indirect effects. Nineteen accessions were evaluated for the following traits: number (NCl), thickness (ThCl), length (LCl) and width (WCl) of cladodes, plant height (PH) and plant width (PW), and green matter yield (GMY), dry matter yield (DMY) and dry matter percentage (DM). The correlations were estimated, and path analysis was performed by the method proposed by Wright. GMY was strongly correlated with DMY, allowing indirect selection for this trait. NCl and ThCl had a direct effect on GMY and can be used for indirect selection or as secondary traits in the selection process. Given the lack of significant correlations between MS and DMY, it is possible to select a palm variety with high DMY and DM. Nevzat et al. (2010) conducted a research to determine characters effecting grain yield in fifty bread wheat (Triticum aestivum L.) cultivars and advanced lines by using simple correlation coefficient and path analysis under 2 locations (high rainfall and low rainfall; 745 and 506 mm, respectively). A total of 50 genotypes, 25 for each independent experiment, were tested for grain yield, test weight, 1000-kernel weight, Zeleny sedimentation, protein content and plant height. Grain yield was significantly correlated with plant height in high rainfall condition. It was significantly correlated with all components except Zeleny sedimentation and protein content in low rainfall condition. Results suggest that plant height and test weight are the primary selection criteria for improving grain yield in bread wheat in high and low rainfall conditions. Hassan (2011) carried out an investigation using 13 alfalfa accessions during 2009 Evaluation of Statistical Models for Climatic Characterization of GKVK Station 17

32 to 2010 and located on the experimental field of East Azarbaijan Agriculture and Natural Resources Research Center (AZARAN), Iran. The main objective of his research was to evaluate positive effect and reliability of yield and quality traits as selection criteria in alfalfa breeding. Significant differences were observed for most of the yield and quality components. Variability coefficients were high for yield components, while quality traits showed relatively low variation. Plant height (PH), number of stems (NS), number of nodes (NN) and leaf size (LS) were positively correlated with plant yield. Crude protein (CP) content was correlated directly with acid detergent fiber (ADF) and natural detergent fiber (NDF) while the correlation with crude fiber (CF) was inverse. The direct effect of the number of stems on yield had the highest value (0.698, P < 0.01). Direct effect of number of nodes on yield was positive (0.508, P < 0.01). Only the plant height had a lower direct effect than the indirect effect on yield. The direct effects of independent traits on CP were significant (P < 0.05), except for trait leaf dry weight (LDW). The direct effect of ADF had the highest value (2.440, P < 0.01), which was positive and significant. LDW trait had a negative direct effect on the CP ( , P < 0.05), while the indirect effect on the NDF trait on CP was neither high nor justifiable. ADF, CF, NDF and LDW traits had the highest indirect effects on CP trait via correlation with in vitro dry matter digestibility (IVDMD). XueMei et al. (2011) collected the data of daily mean air temperature, precipitation and runoff during the period of in the Kaidu River watershed, this paper analyzed the changes in air temperature, precipitation and runoff and revealed the direct and indirect impacts of daily air temperature and precipitation on daily runoff by path analysis. The results showed that mean temperature, time series of the annual, summer and autumn had a significant fluctuant increase during the last 50 years (P < 0.05). Only winter precipitation increased significantly (P < 0.05) with a rate of mm/10a. The annual and winter runoff depths in the last 50 years significantly increased with the rates of 7.11 mm/10a and 1.85 mm/10a, respectively. The driving function of both daily temperature and precipitation on daily runoff in annual and seasonal levels is significant in the Kaidu River watershed by correlation analysis. The result of path analysis showed that the positive effect of daily air temperature on daily runoff depth is much higher than that of daily precipitation in annual, spring, autumn and winter, however, the trend is opposite in the summer. Shweta et al. (2012) conducted correlation and path coefficient analyses in order to identify potential traits determining yield under terminal heat stress during grain filling stages and to study character associations among 25 morphological and physiological traits in a set of genetically divergent 36 bread wheat genotypes under normal and late sown conditions. Grain yield depicted positive association with grain weight (GW), grains/spike (G/S), biological yield (BY), harvest index (HI) under normal (NS) and late sown (LS) conditions, GGR-3 (grain growth rate at 28 days after anthesis) in LS and negative association with DM (days to maturity), HU (heat units) and PTU (photo 18 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

33 thermal unit) in NS and with ChA-1 (CHB a at anthesis), CHB a -1 (total CHB at anthesis) and CHB a -2 (total CHB at 28 days after anthesis) in LS. Path coefficient analysis revealed that, out of eight characters significantly related to yield under normal sown conditions, PTU, harvest index, biological yield and grains per spike directly affected the grain yield in a positive direction. On the basis of strong association with yield and marked direct influence on yield, the numbers of grains were considered to be first order yield components and ought to have top priority in selection under normal sown conditions. Under late sown conditions BY, harvest index (HI), CHB a -1 and grain weight had a direct positive effect, but the grains per spike exhibited negative direct effect despite of the fact that it possessed a positive significant correlation with grain yield. Based on the results, it suggested that higher numbers of grains and high grain weight should be given priority for selection of high yielding genotypes in NS and LS, respectively. Khan and Naqvi (2012) conducted a research with the objectives to identify the utmost traits that may be beneficial for the higher productivity of the grains under stress environment. They used thirteen genotypes which were obtained from different sources as research material. An experiment was carried out in randomized complete block in a split plot design. Water regimes (irrigated and non-irrigated) were allocated to the main plots and genotypes to the subplots. Path coefficient could be used as an important tool to bring about appropriate cause and effect relationship between yield and yield components. According to obtained results the selection on the basis of number of spikes, number of spikelet s and number of grains in this material would likely to be most useful for increasing grain yield because of their direct positive contribution to grain yield under irrigated condition. However a number of spikes, spikelet numbers, spike length and grains number may be used as an effective selection criterion for increasing grain yield of wheat under different irrigation levels. Therefore, they concluded that these traits could be selected for the different stress environments and it would be beneficial for the yield. Sileshi and Temesgen (2013) determined the factors that influence stem radial growth of juvenile Eucalyptus hybrids grown in the east coast of South Africa. Measurement of stem radius was conducted using dendrometers on sampled trees of two Eucalyptus hybrid clones (Eucalyptus grandis Eucalyptus urophylla, GU and E. grandis Eucalyptus camaldulensis, GC). Daily averages of climatic data (temperature, solar radiation, relative humidity and wind speed) were simultaneously collected with total rainfall from the site. They employed path analysis for this study. The joint effect of the climatic variables as well as the direct effect of each climatic variable was studied. Bootstrap estimation procedures, which relax the distributional assumption of the maximum likelihood estimation method, were used. It is found that all variables had a positive effect on stem radial growth. The study showed that tree age is the most important determinant of radial measure. Saied et al. (2013) carried out a study on the relationship among some agronomic Evaluation of Statistical Models for Climatic Characterization of GKVK Station 19

34 and morphological attributes with protein yield of bread wheat cultivars in randomized complete block design. Relationships among measured traits were assessed by the phenotypic correlation coefficient, step-wise regression and path analysis. Correlation analysis showed the significant relation of seed and protein yield with all the traits except correlation of grain filling duration and peduncle length with seed yield as well as spike yield and peduncle length with protein yield. Regression analysis by using step-wise method showed that 83.8 percent of total variation existed in protein yield accounted for by traits entered into the regression model, namely; grain filling rate, grain filling duration, peduncle length and spike harvest index. Grain filling rate accounted for the highest amount of variation about 46.5 per cent. Other traits accounted for 20, 11.3 and 6 percent of variation of protein yield, respectively. Because of that, these traits are recommended as the main components of protein yield in bread wheat cultivars. Path analysis showed that the traits grain filling rate and grain filling duration had the highest and positive direct effects on protein yield. On the other hand, peduncle length and spike harvest index showed low positive and negative direct effect on protein yield. Overall, grain filling rate and grain filling duration were recognized as the best indirect selection criteria to improve protein yield in bread wheat cultivars via selection especially from preliminary generations of breeding programs. Yogesh and Preeti (2014) studied the abundance of the ladybird beetles (LBB), Coccinella septumpunctata L. (Coleoptera- Coccinellideae) in relation to climatic factor in two consecutive cropping seasons and The study revealed that, the cotton LBB was first recorded in the 27 th SMW i.e. the first week of July during both the year of study and remained up and remained active till 50th SMW (2 nd week of December). The peak population was observed (9.76 / 5 plant) during 37 th SMW i.e. 3 rd week of September. The correlation studies between Population of LBB and weather factors showed that the LBB population had a significant positive correlation with maximum temperature (0.542) & minimum temperature (0.560). The multiple coefficient value between the LBB population and group of variable clearly indicated that per cent change in LBB population was affected by maximum temperature, minimum temperature, morning relative humidity, evening relative humidity, sunshine hours, wind velocity, rainfall and rainy days respectively. The data also revealed that per cent variation caused by inexplicable reason or due to error beyond the control of the experiment or due to factors not included in the investigation. Path coefficient analysis of LBB population and abiotic factors revealed that the minimum temperature had positive and high direct effect (1.6592) followed by morning relative humidity (0.1972), rainfall (0.1535), and sunshine hours (0.1519) and evening relative humidity (0.016) respectively. 20 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

35 III MATERIAL AND METHODS This chapter deals with a brief description of the materials which provide the necessary database for the study along with the statistical tools employed in the analysis. The methodology is presented under the following headings: 3.1. Description of the study area Data description Analytical tools and techniques applied Descriptive statistics 3.5. Fitting of probability distributions 3.6. Regression analysis 3.7. Path analysis 3.1 DESCRIPTION OF THE STUDY AREA The present study was conducted to know the climatic characterization of Gandhi Krishi Vignana Kendra (GKVK) station. The station GKVK is located at Bengaluru Urban District of Bengaluru North Taluk. It belongs to the Eastern Dry zone (Zone-V). The geographical co-ordinates of this station are 77 o 35 longitude, 12 o 58 latitude and 930 amsl altitude. The Eastern dry zone includes Kolar, Tumakuru, Bengaluru (Urban), Bengaluru (Rural), Chikkaballapur and Ramanagara. The zone has an average annual rainfall of 944 mm ( ). Major crops grown in this zone are ragi, groundnut, rice and maize. 3.2 DATA DESCRIPTION The present study was based on the secondary data on weather parameters over a period of 38 years ( ) which was collected from AICRP on Agro Meteorology, University of Agricultural Sciences, GKVK, Bengaluru. Weather parameters: 1. Temperature ( o C) 2. Relative Humidity (%) 3. Rainfall (mm) 4. Potential Evapotranspiration (mm) 5. Number of rainy days Evaluation of Statistical Models for Climatic Characterization of GKVK Station 21

36 6. Sunshine hours (hrs/day) 7. Vapor Pressure (pa) 8. Wind Speed (km/hr) 9. Evaporation (mm) 10. Cloud Amount (Oktas/day) 11. Soil temperature ( o C) 1.Temperature ( o C) Temperature is one of the important weather parameters which regulate the growth of a crop in any region. The mean atmospheric temperature is derived using maximum and minimum temperatures observed during 24 hours. Daily maximum and minimum temperatures were recorded using thermometer. 2. Relative Humidity (%) Relative humidity is the moisture content in the atmosphere which is indirectly calculated using the dry bulb and wet bulb thermometer readings recorded at 7.00 A.M. IST (RH I) and 2.00 P.M. IST (RH II). Average of these two time periods was considered in the present study. 3. Rainfall (mm) The available daily rainfall data was used to compute weekly, monthly, seasonal and annual rainfall. 4. Potential Evapotranspiration (mm) It is the maximum amount of soil moisture that soil and a crop canopy jointly evaporate and transpire under no limiting availability of soil moisture. 5. Number of rainy days If 2.5 mm or more rainfall is received in an area in 24 hours, then such a day is considered as rainy day. The number of rainy days in each year was recorded. 6.Sunshine hours(hrs/day) It is the duration of sunshine in a given period for a given location which is expressed as an average of several years. 22 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

37 7.Vapor Pressure (pa) It is an indication of a liquid's evaporation rate and relates to the tendency of particles to escape from the liquid (or a solid). 8. Wind Speed (km/hr) It is an average velocity of the atmosphere recorded at 10 feet height above the ground level. 9. Evaporation (mm) Evaporation is a type of vaporization of a liquid that occurs from the surface of a liquid into a gaseous phase that is not saturated with the evaporating substance. 10. Cloud Amount (Oktas/day) It is the extent of cover by clouds in the portion of the sky visible from a particular location. 11. Soil temperature ( o C) The temperature measured at a given soil depth, typically 5cm, 10cm and 15 cm at 7.00 A.M. (Soil temperature I) and 2.00 P.M. (Soil temperature II). It is the measurement of the warmth of the soil. CROP DATA Finger millet (GPU 28) crop yield (qtl/ha) data during kharif season for a period of 16 years from was collected from ZARS, UAS, Bengaluru. This data was used to know the influence of weather parameters on the crop yield. 3.3 ANALYTICAL TOOLS AND TECHNIQUES APPLIED To assess the climatic characterization of the GKVK station, following statistical tools were used Descriptive Statistics The Descriptive Statistics such as minimum, maximum, mean, standard deviation, skewness and coefficient of variation (CV) were calculated for the weather parameters. Standard formulae were used to compare the above measures. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 23

38 3.3.2 Fitting of probability distributions Rainfall is highly variable in a given period. Hence, there is a need to study the distribution both on long term (annually or seasonally) as well as short term (monthly or weekly) basis. Among the weather parameters, amount of daily maximum rainfall (mm) was considered to fit appropriate probability distributions. The probability distributions viz. normal, log normal, Gamma (1P, 2P, 3P), generalized extreme value (GEV), Weibull (1P, 2P, 3P), Gumbel and Pareto were used to evaluate the best fit probability distribution for rainfall. Description of parameters: Shape parameter A shape parameter is any parameter of a probability distribution that is neither a location parameter nor a scale parameter (nor a function of either or both of these only, such as a rate parameter). Shape parameters allow a distribution to take on a variety of shapes, depending on the value of the shape parameter. These distributions are particularly useful in modeling applications since they are flexible enough to model a variety of data sets. Examples of shape parameters are skewness and kurtosis. Scale parameter In probability theory and statistics, a scale parameter is a special kind of numerical parameter of a parametric family of probability distributions. The larger the scale parameter, the more spread out the distribution. The scale parameter of a distribution determines the scale of the distribution function. The scale is either estimated from the data or specified based on historical process knowledge. In general, a scale parameter stretches or squeezes a graph. The examples of scale parameters include variance and standard deviation. Location parameter The location parameter determines the position of central tendency of the distribution along the x-axis. The location is either estimated from the data or specified based on historical process knowledge. A location family is a set of probability distributions where µ is the location parameter. The location parameter defines the shift of the data. A positive location value shifts the distribution to the right, while a negative location value shifts the data distribution to the left. Examples of location parameters include the mean, median, and the mode. 24 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

39 Table 3.3.2: Description of various probability distribution functions. Distribution Probability density function Range Parameters Gamma (1P) 1 k 1 f ( x) = x exp( x) 0 x <+ Γ( k) k >0 Gamma (2P) and β > 0 k 1 ( k µ ) / β f ( x) = ( x) exp k Gamma (3P) GEV f ( x) = ( x µ ) k 1 1 exp β f ( x) = 1 exp β exp β τ ( k) ( k µ ) / β k β Γ( k) 1 ( 1+ kz) k ( 1+ kz) [ z exp( z) ] 1 1/ k k 0 k = k z > 0 for k 0 - < x <+ for k=0 where, ( µ) z = x β k = Shape parameter β = Scale parameter µ = Location parameter Γ = Gamma function k = Shape parameter β = Scale parameter µ = Location parameter Normal Log- normal Gumbel Pareto f ( x ) f ( x) f ( x) Where, f ( x µ ) = exp 2 σ 2π 1 exp xσ 2π 2σ 2 ( ln( x) µ ) = 2 1 = exp β z = x µ β kβ β k ( x) = k + 1 Weibull (1P) β β f x = β x 1 ( ) exp( x ) Weibull (2P) Weibull (3P) f ( x) = f ( x) = k x β β k 1 k x µ β β 2σ z ( z + e ) x exp β k 1 k x µ exp β k - < x <+ - < µ <+ σ > 0 µ = Mean σ = Standard deviation 0 < x < + σ = Scale parameter µ = Location parameter β >0 - < x <+ 1 x + k, β > 0 x 1 x > 0 β > 0 0 x <+ k, β, > 0 β = Scale parameter µ = Location parameter k = Shape parameter β = Scale parameter β = Shape parameter k = Shape parameter β = Scale parameter µ = Location parameter ( µ 0 yields the two parameter Weibull distribution) Evaluation of Statistical Models for Climatic Characterization of GKVK Station 25

40 A conceptual understanding of the parameters The shape parameter describes the form of the curve. Distributions with a low shape parameter are positively skewed, and as the shape value increases the distribution curve becomes more symmetrical. The mapped parameters provide some spatial context to the rainfall values and distributions. The term shape-dominated refer to the distributions with a larger shape parameter, and scale dominated refer to the distributions with a larger scale parameter. Figure 2 shows a conceptual graphic with the shape and scale parameters on the x- axis and y-axis, respectively. The axes on this graphic are not numerated as the concept of small and large values may vary in time and space. This graphic provides an idea of how the parameters describe regimes. Areas in the low rainfall region of the plot describe study periods that are typically arid during the interval of analysis. The empty area in the graph indicates that areas with at least a minimal amount of rainfall are in either the shape-dominated or scale dominated category. A shape-dominated regime describes a pattern where the rainfall tends to be symmetrically distributed, indicating that drier-than-average events are as common as wetter-than-average events. Scale-dominated rainfall describes locations where the variance is quite large in comparison to the mean. Scale - dominated Scale Small Large Variable rain More extreme events Low rainfall Average is unable to sustain agriculture Shape dominated Consistent rain Fewer extreme events Small Shape Large Fig.2: Conceptual regimes described in parameter space with the shape parameter on the x-axis and the scale parameter on the y-axis. 26 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

41 Testing for goodness of fit: The goodness of fit test measures the discrepancy between observed values and the expected values. Kolmogorov- Smirnov test was used to test for the goodness of fit. In the present investigation, the goodness of fit test was conducted at α = 0.05 level of significance. It was applied for testing the following hypothesis: H 0 : The maximum daily rainfall data follows a specified distribution. H 1 : The maximum daily rainfall data does not follow a specified distribution. Kolmogorov- Smirnov test ( K-S test) This test was used to decide whether a sample comes from a hypothesized continuous PDF. The KS test compares the cumulative distribution functions of the theoretical distribution the distribution described by the estimated shape and scale parameters with the observed values and returns the maximum difference between these two cumulative distributions (Wilks, 1995). This maximum difference in cumulative distribution functions is frequently referred to as the KS-statistic. It is based on the empirical distribution function i.e., on the largest vertical difference between the theoretical and empirical cumulative distribution functions which is given as: i 1 i D = max F i, 1 i n n n ( X ) F( X ) Where, X i = Random sample, i= 1, 2,, n. i CDF 1 = Fn ( X ) = [Number of observations x] n Plot of the best fitted distribution functions for maximum daily rainfall (mm) The frequency density at different levels of rainfall was calculated using the following formulae. Relative Frequency = Class frequency Total frequency Frequency Density = Relative frequency Class width The above analysis will be carried out using XLSTAT software. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 27

42 3.3.3 REGRESSION ANALYSIS The influence of weather parameters on the crop yield was assessed separately by relating finger millet crop yield with the weather parameters through linear functional forms. For this, regression models were used. Regression Data Finger millet crop yield was taken as the dependent variable (Y) and weather parameters (X 1, X 2, X 3, X 4, X 5, X 6, X 7, and X 8 ) were considered as independent variables. X 1 : Minimum temperature (Min Temp.) X 2 : Maximum temperature (Max Temp.) X 3 : Relative Humidity (RH) X 4 : Rainfall X 5 : Number of rainy days X 6 : Potential Evapotranspiration (PET) X 7 : Sunshine hours (SSH) X 8 : Soil temperature (Soil Temp.) Simple Linear Regression Analysis The simple linear regression models were fitted using crop yield as dependent variable (Y) and weather parameters as independent variables (X 1, X 2,, X 8 ) one at a time. Y β β X + e where, i = 1,2,,8. i = o + 1 i i Multiple Linear Regression Analysis The multiple linear regression model for the crop yield (Y) and the weather parameters (X i ) is given as, Y = Xβ + ε Where, Y = ( Y Y,..., ) finger millet. is the vector of values of the dependent variable, crop yield of 1, 2 Y n 28 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

43 X = X X... X n X X X 1 k 2 k... nk, a (n x p) matrix = ( β β ) is the vector of the parameters and ε = ( ε ε,..., ) β,..., 1, 2 β k is the error 2 vector. The error vector ε is assumed to be normally distributed with N [ 0, σ ]. The least square estimator of the parameter vector β is 1 ( X X ) X Y ˆ β = b = And its sampling variance is V 2 ( b) = ( X X ) 1 σ The fitted equation is Y ˆ = Xb The variance of Ŷ at the given values of explanatory variables is ( ˆ ) = σ X 0( X X ) X 0 V Y Where, = ( x, x x ) x ,..., 0k 1, 2 is the vector of given values of the predictor variables. An unbiased estimator of ˆ σ 2 σ ( y y ) n ˆ i i = ( Y Y b X Y ) ( n k) = MSres = i= 1 n k Where, k is the number of parameters to be estimated. Testing significance of multiple linear regression model The significance of multiple linear regression models was tested using F- test. Test Hypothesis: H : β i = β i j 0 j ε k Evaluation of Statistical Models for Climatic Characterization of GKVK Station 29

44 H 1 : β β i j Level of significance: α = 5% Test statistic: F = MS MS Reg Res ANOVA Table Source of variation Degrees of freedom Sum of Squares Mean Sum of Squares F Regression k SS Reg MS Reg MS MS Reg Res Error n-k-1 SS Res MS Res Total n-1 SS T SS Re g = ˆ β X Y ( Yi ) n 2 SS s = Y Y βˆ X Y Re SS T = Y Y ( Yi ) n 2 MS MS Re g = Re s SS Re g k SS Re s = n k 1 Critical region: Reject H 0 at α level of significance if F cal F( k, n k 1) Testing significance of regression coefficient The significance of regression coefficients was tested using t- test. 30 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

45 Test Hypothesis: H : β 0 j = 1,2,..., k 0 j = H 1 : β j 0 Level of significance: α = 0.05 Test statistic: ˆ β j t = ~ tα ( n k 2 ) df SE ( ˆ β j ) Where, SE( ˆ β ) = j MS S Re s xx S xx ( X ) 2 = i X Model Adequacy Checking Adequacy of a model depends on the validity of assumptions underlying the model. The assumptions made in a linear regression model are, apart from lineardependence of Y on predictors, independence and identical distribution (normal) of the predictors with zero mean. Gross violations of the assumptions may yield an unstable model in the sense that a different sample could lead to a totally different model with opposite conclusions. We cannot detect departures from the underlying assumption by examination of the summary statistics such as t or F statistic or R 2 (Draper and Smith, 1981). These are Global model properties and as such they do not ensure model adequacy. Hence, diagnostic methods, primarily based on the study of the model residuals are used. The diagnostics checks for the present study used are for randomness and normality of residuals. Coefficient of Determination (R 2 ) The statistic R 2 indicates the percentage of variation in the dependent variable explained by the independent variables included in the regression model. R 2 is computed using the uncorrected sum of squares. R n 2 i= 1 = n ( yˆ y) ( yi y) i= 1 Where, 0 R 2 1 i 2 2 Evaluation of Statistical Models for Climatic Characterization of GKVK Station 31

46 Assumptions of error term An important assumption of multiple regression models is that the residual, ε follows normal distribution. This assumption is required to test the hypothesis about the regression coefficients. This assumption is verified using, i. Randomness of error terms ii. Normal probability plots. Randomness of error terms (Runs test) Randomness of residuals is tested using the non-parametric one sample runs test. A run is defined as the succession of identical symbols which are followed and preceded by different symbols or no symbols at all. The null and alternative hypothesis is given by H 0 : The sequence of residuals is random. H 1 : The sequence of residuals is not random. Let n 1 be the number of elements of one kind and n 2 be the number of elements of the other kind in a sequence of N = n 1 + n 2 binary events. If both n 1 and n 2 are less than or equal to 20, and if the number of runs (r) falls between the critical values provided in the tables, we cannot reject the null hypothesis. Normal probability plot Suppose e 1, e 2,,e n are the residuals ranked in increasing order as e 1 < e 2 < < e n. If we plot e i against the cumulative probabilities P i =(i-1/2)/n, i=1,2,,n, in the normal probability plot. The residuals can be assumed to be normally distributed if the resulting points lie approximately on a straight line PATH ANALYSIS Path analysis quantifies the direct and indirect effect of each variable and helps to visualize correlation in its totality to determine selection indices. It is an extension of the regression model, which researchers use to test the fit of a correlation matrix with a causal model that has been tested. Path analysis is a standardized partial regression and a technique to elucidate the contribution of each variable to the total correlation coefficient. The aim of the path analysis is to provide estimates of the magnitude and significance of the hypothesized causal connections among sets of variables displayed through the use of path diagrams. 32 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

47 Path analysis was used to examine the influence of weather parameters on the rainfall. The direct and indirect effects of weather parameters on the rainfall were identified considering the following weather variables: X 1 : Minimum temperature (Min Temp.) X 2 : Maximum temperature (Max Temp.) X 3 : Potential Evapotranspiration (PET) X 4 : Vapor Pressure (VP) X 5 : Relative Humidity (RH) X 6 : Wind Speed (WS) X 7 : Sunshine hours (SSH) X 8 : Cloud Amount (CL AMT) X 9 : Evaporation (EVP) PATH DIAGRAM The graphical or diagrammatical representation of path analysis is known as path diagram. Generally, arrows are used to show the relations. A single headed arrow is meant for showing cause to effect. In path diagram, there are two kinds of variables, independent variables known as the causal variables and the dependent variables known as endogenous variables. If Y is a character, which is linearly determined by correlated variables, X 1, X 2 and X 3 because of the closed system, a path diagram can be formulated as in Fig.3. X 1 a r 12 Y b X 2 r 13 c r 23 Residual X 3 Fig.3: Cause and effect relationship Evaluation of Statistical Models for Climatic Characterization of GKVK Station 33

48 Symbols a, b and c represent path coefficients from causes X 1, X 2 and X 3 respectively. r 12, r 23 and r 13 are the correlation coefficients and R is the variation in Y that is due to the unknown agencies (residuals). A character Y is determined by three correlated variables, viz. X 1, X 2 and X 3. The relationship can be expressed in the form of a partial regression equation and is given by, Y µ + β X + β X + X + R = β 3 3 Where, β 1, β2 and β 3 are the partial regression coefficients. β 1 is the partial regression coefficient of Y on X 1, which measures the amount of change that can be brought about in Y due to one unit change in X 1, when X 2 and X 3 are held constant. β 2 and β 3 have similar meanings. R is the residual component. When the variables are expressed as deviation from their respective means and also in terms standardized units, their means will be zero and variance will be equal to unity. Thus, if we denote ( X 1 X 1 )/ σ Y by a lower case letter X 1, mean of X 1 will be zero and variance will be unity. Replacing the expression β 1σ X 1 / σ Y by term a which is a standardized partial regression coefficient and it is called as the path coefficient. It should be pointed out that the path coefficients have a direction similar to the regression coefficients, i.e., they may bear a positive or negative sign. Unlike the correlation coefficients, the path coefficients could be greater than or lesser than unity like regression coefficients. A set of three simultaneous equations is then formulated from the path diagram as: r = a + br + cr YX 1 X 1 X 2 X 1 X 2 r = ar + b + cr YX 2 X 2 X X 1 2 X 3 ryx = ar + br + c 3 X 2 X1 X 3 X 2 r = Estimates of simple correlation coefficient between X 1 and X 2. X 1 X 2 r = Estimates of simple correlation coefficient between X 2 and X 3. X 2 X 3 r = Estimates of simple correlation coefficient between X 1 and X 3. X 1 X 3 a = Direct effect of variable X 1. b = Direct effect of variable X 2. c = Direct effect of variable X 3. Y = Dependent variable. Information provided in the above equation can be arranged in the matrix form as: 34 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

49 1 rx 1 r X1 X X 2 3 r X X r X X 2 r X X r X X 3 a b c = r r r YX YX YX Which is expressed as A B = C Here the matrix A is symmetrical about the diagonal elements. The diagonal elements in this matrix represent the direct effects of X 1, X 2 and X 3, respectively on Y. All the direct effects have been replaced by unity in this matrix. Elements of the column matrix B specify the path coefficients to be estimated. The column matrix C represents the estimated values of correlation coefficients between dependent variable Y, and the component variables, i.e., X 1, X 2 and X 3. Estimates of unknown are calculated as: B = A -1 C Where, A -1 = Inverse of matrix A. C = Estimates of simple correlation coefficient between dependent variables. The indirect effects are calculated from the three simultaneous equations which provide the estimated path coefficients. This method can permit us to identify relatively important components of a dependent variable, on the basis of their direct and indirect influences. Besides the direct effect of the exogenous variables (X i ) on the endogenous (Y), there is an indirect effect of X i on Y via other X s by virtue of its relationship with others. A change in variable, say X 1, will affect its linear correlation r 12, with another variable, X 2, which invariably affects Y. Let P i be the direct effect of X i on Y, and change is only partial and proportional to r 12. That is, r ij P i is an indirect effect of X i via X j. The percentage direct contributions of any variable and its residual effect are worked out as follows: Direct percentage contribution = (P i ) Where, P i is the direct effect of X i on Y. Once the path coefficient values corresponding to different independent variables are worked out, the next problem lies in estimating the residual (R) i.e., the unexplained part of the model. The residual is given by: Where, R = 1 - r ij P i r ij = Correlation between i th independent variables and the response variable. P i = Direct effect. The regression analysis and path analysis will be carried out using MS Excel. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 35

50 IV RESULTS Keeping in view the specific objectives of the present investigation, the data on the weather parameters and the finger millet crop was subjected to statistical analysis using appropriate techniques. The results of the analysis are presented under the following headings: 4.1 Fitting appropriate probability distributions for rainfall. 4.2 Impact of weather parameters on the yield of finger millet crop. 4.3 Direct and indirect effects of weather parameters on rainfall. 4.1 Fitting appropriate probability distributions for rainfall. Rainfall data was classified into 28 sets viz., 1 annual, 1 seasonal, 5 seasonal (June to October) months and 21 seasonal Standard Meteorological Weeks (23 rd SMW- 43 rd SMW) to study the distribution pattern of rainfall at different levels. The best fit probability distributions for maximum daily rainfall on different sets of data were identified based on the criteria of Kolmogorov-Smirnov test Descriptive statistics The summary of statistics such as maximum, minimum, mean, standard deviation, skewness and coefficient of variation values of maximum daily rainfall are presented in Table The results show that both annual and seasonal maximum rainfall was observed to be 200 mm. Monthly maximum rainfall during monsoon season ranged from 98.6 mm to 200 mm while weekly maximum rainfall was between 43.2 mm to 200 mm. It was also observed that the minimum rainfall was 56.8 mm annually, 43.6 mm seasonally and it ranged from 6.4 mm to 16 mm monthly. For all the weeks, the minimum rainfall was found to be 0.0 mm. Mean annual rainfall was found to be 94.6 mm, whereas for overall seasonal months, it was 90.6 mm. During the seasonal months the mean rainfall ranged from 32.1mm to 67.6 mm while for weekly, it varied from 5.2 mm to 33.5 mm. The value of coefficient of variation annually was observed to be 39.2 per cent which was the lowest among all the data sets. Seasonally, the rainfall varied by about 43.6 per cent whereas monthly, the variation ranged from 58.6 to 66.8 per cent. For Standard Meteorological Weeks, the rainfall variation was found to be very high ranging from 86.9 to per cent. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 36

51 Table 4.1.1: Summary of statistics for maximum daily rainfall (mm) Study Period Range Mean Standard deviation Parameters CV (%) Skewness Maximum Minimum Annual 1Jan 31 Dec Seasonal 1 June- 28 Oct June 1 June-30 June July 1 July-31 July August 1 Aug-31 Aug September 1 Sept-30 Sept October 1 Oct- 28 Oct rd SMW 4 June-10 June th SMW 11 June-17 June th SMW 18 June-24 June th SMW 25 June-1 July th SMW 2 July-8 July th SMW 9 July-15 July th SMW 16 July-22 July th SMW 23 July-29 July st SMW 30 July-5 Aug nd SMW 6 Aug-12 Aug rd SMW 13 Aug-19 Aug th SMW 20 Aug-26 Aug th SMW 27 Aug-2 Sept th SMW 3 Sept-9 Sept th SMW 10 Sept-16 Sept th SMW 17 Sept-23 Sept th SMW 24 Sept-30 Sept th SMW 1 Oct-7 Oct st SMW 8 Oct-14 Oct nd SMW 15 Oct-21 Oct rd SMW 22 Oct-28 Oct Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

52 250 Annual Maximum rainfall (mm) Year Fig. 4: Variation in Annual maximum daily rainfall (mm) Seasonal Maximum rainfall (mm) Year Fig. 5: Variation in Seasonal maximum daily rainfall (mm) Evaluation of Statistical Models for Climatic Characterization of GKVK Station 38

53 The coefficient of skewness measures the asymmetry of the frequency distribution of the data. For the present study, the coefficient of skewness for all the data sets ranged from 0.6 to 4.5 indicating positive skewness. From this result, we can conclude that in all the data sets the average rainfall for the period exceeded the modal value. Figure 4 shows the year wise variation in annual maximum daily rainfall from It ranged from 56.8 mm (in 2006) to 200 mm (in 1997). The minimum of the maximum daily rainfall, i.e., 56.8 mm occurred during 10 th SMW (5 th Mar-11 th Mar) and the maximum of the maximum daily rainfall of 200 mm occurred during the 41 st SMW (8 th Oct-14 th Oct). Figure 5 represents the year wise variation in seasonal maximum daily rainfall from and it ranged from 43.6 mm (in 2009) to 200 mm (in 1997). The minimum of the maximum daily rainfall of 43.6 mm occurred during 30 th SMW (23 rd Jul- 29 th July) and the maximum daily rainfall occurred during the 41 st SMW (8 th Oct-14 th Oct). From both the figures it reveals that the variation in rainfall is stationary over the years as well as during the seasons Study period wise fitting of probability distributions The probability distributions used to evaluate the best fit for maximum daily rainfall were Normal, Log normal, Gamma (1P, 2P, 3P), Generalized Extreme Value (GEV), Weibull (1P, 2P, 3P), Gumbel and Pareto. The goodness of fit for different probability distributions was tested using Kolmogorov- Smirnov test (K-S test). The test statistic D along with the p-values for each data set was computed for 11 probability distributions. Table represents the study period wise probability distributions and their goodness of fit statistics. Annual For the data set on annual maximum daily rainfall, the calculated K-S test statistic value was found to be non-significant for the probability distributions such as, Gamma (2P), Log- normal and Weibull (2P, 3P). This means that the null hypothesis which states that the rainfall data follow specified distributions cannot be rejected and hence it fits the above mentioned probability distributions. The test statistic value for Gamma (2P) distribution was and the p-value was found to be while log- normal distribution had a test statistic value of and a p-value of The test statistic value for Weibull (2P) was found to be and the p-value was whereas Weibull (3P) had a test statistic value of and a p-value of Among the four probability distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. Seasonal (June-October) During the seasonal maximum daily rainfall, the distributions which fitted well were Gamma (2P), Log- normal and Weibull (2P, 3P). Gamma (2P) distribution had a test statistic value of and a p-value of while the test statistic value for lognormal distribution was found to be with a p-value of The test statistic 39 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

54 value and p-value for Weibull (2P) distribution was observed to be and respectively, whereas for Weibull (3P) distribution, the test statistic value and p-value were found to be and respectively. Weibull (3P) had the lowest test statistic value and the highest p-value among all the four probability distributions. June For the seasonal month of June, Gamma (2P), Generalized Extreme Value (GEV), log- normal and Weibull (2P, 3P) distributions fitted well. The test statistic value and p- value for Gamma (2P) distribution was observed to be and respectively. GEV distribution had a test statistic value of and p-value of while the test statistic value for log- normal distribution was observed to be with a p-value of The test statistic value and p-value for Weibull (2P) distribution was observed to be and respectively whereas Weibull (3P) distribution had a test statistic value of with a p-value of Among these five probability distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. July It was observed that for the seasonal month of July, the distributions which fitted well were Gamma (2P), log- normal, Pareto and Weibull (2P, 3P). Gamma (2P) distribution had a test statistic value and p-value of and respectively while it was observed to be and respectively for log-normal distribution and for Pareto distribution, the test statistic value and p-value were found to be and respectively. The test statistic value for Weibull (2P) was observed to be with a p-value of whereas Weibull (3P) had a test statistic value of with a p-value of Among these five probability distributions, log-normal had the lowest test statistic value and the highest p-value. August For the seasonal month of August, Gamma (2P), Generalized Extreme Value (GEV), log- normal and Weibull (2P, 3P) were found to be the best fits. The test statistic value and p-value for Gamma (2P) distribution was observed to be and respectively. GEV distribution had a test statistic value of with a p-value of while the test statistic value for log- normal distribution was observed to be with a p-value of The test statistic value and p-value for Weibull (2P) distribution was observed to be and respectively whereas Weibull (3P) distribution had a test statistic value of and p-value of Among these five probability distributions, Weibull (2P) had the lowest test statistic value and the highest p-value. September During the seasonal month of September, the distributions which fitted well were Gamma (2P), GEV, log- normal, Pareto and Weibull (2P, 3P). Gamma (2P) distribution Evaluation of Statistical Models for Climatic Characterization of GKVK Station 40

55 had a test statistic value of and p-value of respectively, while it was observed to be and respectively for GEV distribution. For log-normal distribution, the test statistic value and p-value were found to be and respectively and it was found to be and for Pareto distribution respectively. The test statistic value for Weibull (2P) was observed to be and the p-value was whereas Weibull (3P) had a test statistic value of and a p- value of Among these six probability distributions, Gamma (2P) had the lowest test statistic value and the highest p-value. Maximum number of distributions fitted well for the study period September when compared to all the other study periods. October It was found that for the seasonal month of October, the distributions which fitted well were Gamma (2P), log- normal, Pareto and Weibull (2P, 3P). Gamma (2P) distribution had a test statistic value and p-value of and respectively while it was observed to be and respectively for log-normal and Pareto distribution, the test statistic value and p-value were found to be and respectively. The test statistic value for Weibull (2P) was observed to be with a p- value of whereas Weibull (3P) had a test statistic value of and a p-value of Among these five probability distributions, log-normal had the lowest test statistic value and the highest p-value. 23 rd SMW (4 June-10 June) For this Standard Meteorological Week, Gamma (2P) was the only distribution which fitted well with a K-S test statistic value of and a p-value of th SMW (11 June-17 June) For the 24 th Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was , while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 25 th SMW (18 June-24 June) During the 25 th Standard Meteorological Week, Gamma (2P) and Weibull (3P) distributions were found to be the best fits. Gamma (2P) had a test statistic value of and a p-value of The K-S test statistic value and p-value for Weibull (3P) distribution was observed to be and respectively. Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p- value. 41 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

56 26 th SMW (25 June-1 July) For the 26 th Standard Meteorological Week, the distributions which fitted well were Pareto and Weibull (3P). The test statistic value for Pareto distribution was observed to be and the p-value was , while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 27 th SMW (2 July-8 July) It was observed that for the 27 th Standard Meteorological Week, the distributions such as Gamma (2P) and Weibull (3P) fitted well. Gamma (2P) had a test statistic value of and a p-value of The test statistic value and p-value for Weibull (3P) distribution was observed to be and respectively. Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p- value. 28 th SMW (9 July-15 July) For the 28 th Standard Meteorological Week, Weibull (3P) was the only distribution which fitted well with a K-S test statistic value of and a p-value of th SMW (16 July-22 July) During the 29 th Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was , while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 30 th SMW (23 July-29 th July) It was found that Gamma (2P) and Weibull (3P) distributions fitted well for the 30 th Standard Meteorological Week. Gamma (2P) hada test statistic value of and a p-value of The test statistic value and p-value for Weibull (3P) distribution was observed to be and respectively. Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. 31 st SMW (30 July-5 Aug) For the 31 st Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). Gamma (2P) had a test statistic value of and a p- value of The test statistic value and p-value for Weibull (3P) distribution was Evaluation of Statistical Models for Climatic Characterization of GKVK Station 42

57 observed to be and respectively. Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 32 nd SMW (6 Aug-12 Aug) It was observed that the 32 nd Standard Meteorological Week, the distributions such as Gamma (2P) and Weibull (3P) fitted well. The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Gamma(2P) distribution had the lowest test statistic value and the highest p- value among the two best fitted distributions. 33 rd SMW (13 Aug-19 Aug) During the 33 rd Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was , while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Among the two best fitted distributions, Gamma(2P) distribution had the lowest test statistic value and the highest p-value. 34 th SMW (20 Aug-26 Aug) For the 34 th Standard Meteorological Week, Gamma (2P) and Weibull (3P) distributions were found to the best fits. The K-S test statistic value for Gamma (2P) was found to be and the p-value was while for Weibull (3P) distribution, test statistic value was observed to be and the p-value was found to be Gamma (2P) distribution had the lowest test statistic value and the highest p-value among the two best fitted distributions. 35 th SMW (27 Aug-2 Sept) It was found that for the 35 th Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 36 th SMW (3 Sept-9 Sept) During the 36 th Standard Meteorological Week, Gamma (2P) and Weibull (3P) distributions were found to be the best fits. The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to 43 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

58 be Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. 37 th SMW (10 Sept-16 Sept) For the 37 th Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 38 th SMW (17 Sept-23 Sept) It was observed that for the 38 th Standard Meteorological Week, Gamma (2P) and Weibull (3P) distributions were found to be the best fits. The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. 39 th SMW (24 Sept-30 Sept) For the 39 th Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 40 th SMW (1 Oct-7 Oct) During the 40 th Standard Meteorological Week, Gamma (2P), GEV and Weibull (3P) distributions were found to be the best fits. The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for GEV distribution, test statistic value was found to be and the p-value was observed to be The test statistic value for Weibull (3P) distribution was found to be and the p- value was Among the three best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. 41 st SMW (8 Oct-14 Oct) It was found that for the 41 st Standard Meteorological Week, Gamma (2P) and Weibull (3P) distributions were found to be the best fits. The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed Evaluation of Statistical Models for Climatic Characterization of GKVK Station 44

59 to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. 42 nd SMW (15 Oct-21 Oct) During the 42 nd Standard Meteorological Week, the distributions which fitted well were Gamma (2P) and Weibull (3P). The K-S test statistic value for Gamma (2P) was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Among the two best fitted distributions, Weibull (3P) had the lowest test statistic value and the highest p-value. 43 rd SMW (22 Oct-28 Oct) For the 42 nd Standard Meteorological Week, Pareto and Weibull (3P) distributions were found to be the best fits. The K-S test statistic value for Pareto was observed to be and the p-value was while for Weibull (3P) distribution, test statistic value was found to be and the p-value was observed to be Weibull (3P) had the lowest test statistic value and the highest p-value among the two best fitted distributions. The above results indicated that Weibull (3P) distribution fitted well for majority of the study periods Parameters of the best fit probability distributions The best fit probability distribution for each data set was found out based on the K- S test statistic and the p- value. Among all the fitted probability distributions for maximum daily rainfall, the distribution with lowest test statistic value and highest p- value was considered as the best fit. The best fit probability distribution for each data set along with the parameters is presented in the Table Annual For the annual maximum daily rainfall, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6(A). The distribution followed reverse J- shaped curve on linear axes which approached each of the orthogonal axes asymptotically. Since the distribution was scale dominated, there was large variation in the distribution of rainfall. Seasonal For the seasonal maximum daily rainfall, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6(B) shows the shape of the distribution. The distribution approximated normal with consistency in the distribution of rainfall. 45 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

60 June For the data set on the seasonal month of June, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (C). The distribution was observed to be positively skewed where the average rainfall exceeded the modal rainfall. Since the distribution was scale dominated, there was large variation in the distribution of rainfall. July During the seasonal month of July, log-normal was found to be the best fit probability distribution with shape parameter of and scale parameter of Fig.6 (D) shows the shape of the distribution in which positive skewness was observed with the average rainfall exceeding the modal rainfall. Since the distribution was shape dominated, there was not much variation in the distribution of rainfall. August For the seasonal month of August, Weibull (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter of The shape of the distribution is shown in Fig.6 (E). The distribution showed slight positive skewness with the average rainfall exceeding the modal rainfall. There was large variation in the distribution of rainfall since it was found to be scale-dominated. September During the seasonal month of September, Gamma (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter of Fig.6 (E) shows the shape of the distribution. The shape of the distribution approximated normal and there was consistency in the distribution of rainfall since the distribution was shape dominated. October For the data set on the seasonal month of October, log-normal was found to be the best fit probability distribution with shape parameter of and scale parameter of The shape of the distribution is shown in Fig.6 (F) where the distribution was positively skewed with the average rainfall exceeding the modal rainfall. There was consistency in the distribution of rainfall since it was found to be shape-dominated. 23 rd SMW (4 June-10 June) During this Standard Meteorological Week, Gamma (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter of Fig.6 (a) shows the shape of the distribution where the curve got stretched out to Evaluation of Statistical Models for Climatic Characterization of GKVK Station 46

61 the right and its height has decreased. Since the distribution was scale dominated, large variation was observed in the distribution of rainfall. 24 th SMW (11 June-17 June) During 24 th Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of 0.72, scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (b) where a vertical asymptote near the origin was observed. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 25 th SMW (18 June-24 June) For the data set on 25 th Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (c) shows the shape of the distribution. The distribution got stretched out to the right and its height has decreased. Since the distribution was scale dominated, there was large variation in the distribution of rainfall. 26 th SMW (25 June-1 July) During the 26 th Standard Meteorological Week Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (d) in which the distribution got stretched out to the right and its height has decreased. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 27 th SMW (2 July-8 July) For the 27 th Standard Meteorological Week Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (e) shows the shape of the distribution where the distribution got stretched out to the right and its height has decreased. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 28 th SMW (9 July-15 July) For the data set on the 28th Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter , scale parameter and location parameter The shape of the distribution is shown in Fig.6 (f). The curve had a vertical asymptote near the origin. Since the distribution was scale dominated, there was large variation in the distribution of rainfall. 47 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

62 29 th SMW (16 July-22 July) During the 29 th Standard Meteorological Week Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (g). The distribution showed positive skewnes with the shape parameter value more than 1. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 30 th SMW (23 July-29 th July) For the data set on 30 th Standard Meteorological Week, Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (h) in which positive skewness was observed with shape parameter value more than 1. Since the distribution was dominated by the scale parameter, there was large variation in the distribution of rainfall. 31 st SMW (30 July-5 Aug) During this week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of The shape of the distribution is shown in Fig.6 (i). The distribution showed positive skewness with the shape parameter value more than 1. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 32 nd SMW (6 Aug-12 Aug) For the 32 nd Standard Meteorological Week, Gamma (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter The shape of the distribution is shown in Fig.4 (j). It was observed that the distribution stretched out towards right. Since the distribution was dominated by the scale parameter, there was large variation in the distribution of rainfall. 33 rd SMW (13 Aug-19 Aug) During the 33 rd Standard Meteorological Week, Gamma (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter Fig.6 (k) shows the shape of the distribution where a vertical asymptote near the origin was observed. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 48

63 34 th SMW (20 Aug-26 Aug) For the data set on 34 th Standard Meteorological Week, Gamma (2P) was found to be the best fit probability distribution with shape parameter of and scale parameter of Fig.6 (l) shows the shape of the distribution where a vertical asymptote near the origin was observed. Since the distribution was dominated by the scale parameter, there was large variation in the distribution of rainfall. 35 th SMW (27 Aug-2 Sept) During the 35 th Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6(m) shows the shape of the distribution. The distribution got stretched out to the right and its height has decreased as the shape parameter approached a value 1. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 36 th SMW (3 Sept-9 Sept) For the data set on 36 th Standard Meteorological Week, Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (n) shows the shape of the distribution. Since the shape parameter was near to a value 0.5, the distribution showed vertical asymptote near the origin. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 37 th SMW (10 Sept-16 Sept) During this week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (o) shows the shape of the distribution where the curve stretched out towards right. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 38 th SMW (17 Sept-23 Sept) For the 38 th Standard Meteorological Week, Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (p) shows the shape of the distribution where the curve stretched out towards right. Large variation was observed in the distribution of rainfall as it was dominated by the scale parameter. 49 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

64 39 th SMW (24 Sept-30 Sept) During the 39 th Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (q) shows the shape of the distribution. The distribution showed vertical asymptote near the origin. As the distribution was dominated by the scale parameter, there was large variation in the distribution of rainfall. 40 th SMW (1 Oct-7 Oct) For the data set on 40 th Standard Meteorological Week, Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (r) shows the shape of the distribution. The distribution was observed to be positively skewed where the average rainfall exceeded the modal rainfall. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 41 st SMW (8 Oct-14 Oct) During this Standard Meteorological Week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6(s) shows the shape of the distribution. The distribution showed vertical asymptote near the origin. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 42 nd SMW (15 Oct-21 Oct) For the data set on 42 nd Standard Meteorological Week, Weibull (3P) was observed to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (t) shows the shape of the distribution. Since the shape parameter approached a value 1, the distribution got stretched out to the right and its height has decreased. There was large variation in the distribution of rainfall as it was dominated by the scale parameter. 43 rd SMW (22 Oct-28 Oct) During this week, Weibull (3P) was found to be the best fit probability distribution with shape parameter of , scale parameter of and location parameter of Fig.6 (u) shows the shape of the distribution. The distribution showed vertical asymptote near the origin. The distribution was dominated by the scale parameter which showed large variation in the distribution of rainfall. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 50

65 Table 4.1.2: Study period wise probability distributions using goodness of fit test. Study Range Kolmogorov Smirnov Period Distribution Statistic p- value Gamma (2P) Annual 1Jan 31 Dec Log- normal Weibull (2P) Weibull (3P) Gamma (2P) Seasonal 1 June- 28 Oct Log- normal Weibull (2P) Weibull (3P) Gamma (2P) GEV June 1 June-30 June Log- normal Weibull (2P) Weibull (3P) Gamma (2P) Log- normal July 1 July-31 July Pareto Weibull (2P) Weibull (3P) Gamma (2P) GEV August 1 Aug-31 Aug Log- normal Weibull (2P) Weibull (3P) Gamma (2P) GEV September 1 Sept-30 Sept Log- normal Pareto Weibull (2P) Weibull (3P) Gamma (2P) Log- normal October 1 Oct- 28 Oct Pareto Weibull (2P) Weibull (3P) rd SMW 4 June-10 June Gamma (2P) th SMW 11 June-17 June Gamma (2P) Weibull (3P) Contd. 51 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

66 Table 4.1.2: Study period wise probability distributions using goodness of fit test. Study Kolmogorov Smirnov Range Period Distribution Statistic p- value 25 th SMW Gamma (2P) June-24 June Weibull (3P) th SMW 25 June-1 July Pareto Weibull (3P) th SMW Gamma (2P) July-8 July Weibull (3P) th SMW 9 July-15 July Weibull (3P) th SMW Gamma (2P) July-22 July Weibull (3P) th SMW Gamma (2P) July-29 July Weibull (3P) st SMW Gamma (2P) July-5 Aug Weibull (3P) nd SMW Gamma (2P) Aug-12 Aug Weibull (3P) rd SMW Gamma (2P) Aug-19 Aug Weibull (3P) th SMW Gamma (2P) Aug-26 Aug Weibull (3P) th SMW Gamma (2P) Aug-2 Sept Weibull (3P) th SMW Gamma (2P) Sept-9 Sept Weibull (3P) th SMW Gamma (2P) Sept-16 Sept Weibull (3P) th SMW 17 Sept-23 Sept Gamma (2P) Weibull (3P) th SMW 24 Sept-30 Sept Gamma (2P) Weibull (3P) Gamma (2P) th SMW 1 Oct-7 Oct GEV Weibull (3P) st SMW 8 Oct-14 Oct Gamma (2P) Weibull (3P) nd SMW 15 Oct-21 Oct Gamma (2P) Weibull (3P) rd SMW 22 Oct-28 Oct Pareto Weibull (3P) Evaluation of Statistical Models for Climatic Characterization of GKVK Station 52

67 Table 4.1.3: Parameters of the best fit probability distributions for Maximum daily rainfall. Study Period Range Best fit Shape paramete r (k) Parameters Scale paramete r (β) Location paramet er (µ) Annual 1Jan 31 Dec Weibull (3P) Seasonal 1 June- 28 Oct Weibull (3P) June 1 June-30 June Weibull (3P) July 1 July-31 July Log-normal August 1 Aug-31 Aug Weibull (2P) September 1 Sept-30 Sept Gamma (2P) October 1 Oct- 28 Oct Log-normal rd SMW 4 June-10 June Gamma (2P) th SMW 11 June-17 June Weibull (3P) th SMW 18 June-24 June Weibull (3P) th SMW 25 June-1 July Weibull (3P) th SMW 2 July-8 July Weibull (3P) th SMW 9 July-15 July Weibull (3P) th SMW 16 July-22 July Weibull (3P) th SMW 23 July-29 July Weibull (3P) st SMW 30 July-5 Aug Weibull (3P) nd SMW 6 Aug-12 Aug Gamma (2P) rd SMW 13 Aug-19 Aug Gamma (2P) th SMW 20 Aug-26 Aug Gamma (2P) th SMW 27 Aug-2 Sept Weibull (3P) th SMW 3 Sept-9 Sept Weibull (3P) th SMW 10 Sept-16 Sept Weibull (3P) th SMW 17 Sept-23 Sept Weibull (3P) th SMW 24 Sept-30 Sept Weibull (3P) th SMW 1 Oct-7 Oct Weibull (3P) st SMW 8 Oct-14 Oct Weibull (3P) nd SMW 15 Oct-21 Oct Weibull (3P) rd SMW 22 Oct-28 Oct Weibull (3P) Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

68 4.2 Impact of weather parameters on the yield of finger millet crop Descriptive statistics Descriptive statistics of different weather parameters is given in the Table It was observed that minimum temperature had a mean value of C whereas standard deviation was found to be C. The coefficient of variation was found to be 2.08 per cent and negative skewness was observed. Minimum temperature ranged from C to C. The mean of the maximum temperature was found to be C and standard deviation was C. It had a lesser coefficient of variation of 1.79 per cent and the skewness was observed to be positive. The maximum temperature ranged from C to C. Relative humidity had a mean value of per cent and standard deviation was observed to be 2.07 per cent. The coefficient of variation of the relative humidity was observed to be 2.88 per cent and it was negatively skewed. The relative humidity ranged from per cent to per cent. The mean value of daily rainfall was found to be 4.08 mm while the standard deviation was 1.32 mm. It had a higher value of coefficient of variation i.e., per cent and skewness was found to be positive. The daily rainfall ranged from 2.10mm to 6.28 mm. Number of rainy days had a mean value of 30 and standard deviation was observed to be 7 days. The coefficient of variation was found to be per cent and the skewness was negative. The number of rainy days ranged from 19 to 40. The mean value of PET was observed to be mm and standard deviation was 7.63 mm. It was observed that the coefficient of variation was per cent and PET was found to be positively skewed. PET ranged from mm to mm. It was observed that mean value of the sunshine hours was 4.95 hrs/day while it had a standard deviation of 0.72 hrs/day. The coefficient of variation was observed to be per cent and it was negatively skewed. The sunshine hours ranged from 3.56 hrs/day to 6.18 hrs/day. Soil temperature had a mean value of 26 0 C whereas the standard deviation was observed to be C. The coefficient of variation was observed to be 6.23 per cent and soil temperature was negatively skewed. The soil temperature ranged from C to C. Evaluation of Statistical Models for Climatic Characterization of GKVK Station 54

69 Frequency Density Frequency Density Frequency Density Rainfall (mm) Rainfall (mm) Rainfall (mm) Weibull (3)(0.813,37.164,56.350) Weibull (3)(2.150,74.924,12.075) Weibull (3)(1.269,29.740,4.882) Fig.6(A): Distribution of Annual rainfall Fig.6(B): Distribution for Seasonal rainfall Fig.6(C): Distribution of rainfall during June Frequency Density Frequency Density Frequency Density Rainfall (mm) Log-normal(3.514,0.620) Rainfall (mm) Weibull (2)(1.780,48.727) Rainfall (mm) Gamma (2)(2.915,23.178) Fig.6(D): Distribution of rainfall during July Fig.6(E): Distribution of rainfall during August Fig.6(F): Distribution of rainfall during September Fig.6: Plot of the best fitted distribution functions for maximum daily rainfall (mm) Contd.

70 Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Log-normal(4.017,0.618) Gamma (2)(0.947,19.053) Weibull (3)(0.723,8.067,-0.200) Fig.6(G): Distribution of rainfall during October Fig.6(a): Distribution of rainfall during 23 rd SMW Fig.6(b): Distribution of rainfall during 24 th SMW Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Weibull (3)(0.720,7.399,-0.300) Weibull (3)(0.870,5.270,-0.451) Weibull (3)(0.820,14.931,-0.450) Fig.6(c): Distribution of rainfall during 25 th SMW Fig.6(d): Distribution of rainfall during 26 th SMW Fig.6(e): Distribution of rainfall during 27 th SMW Fig.6: Plot of the best fitted distribution functions for maximum daily rainfall (mm) Contd.

71 Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Frequency Density Rainfall (mm) Weibull (3)(0.821,13.060,-0.450) Weibull (3)(1.028,17.958,-0.451) Weibull (3)(1.175,14.167,-1.560) Fig.6(f): Distribution of rainfall during 28 th SMW Fig.6(g): Distribution of rainfall during 29 th SMW Fig.6(h): Distribution of rainfall during 30 th SMW Frequency Density Frequency Density Frequency Density Rainfall (mm) Rainfall (mm) Rainfall (mm) Weibull (3)(1.149,19.743,-0.694) Gamma (2)(0.943,15.018) Gamma (2)(0.850,20.169) Fig.6(i): Distribution of rainfall during 31 st SMW Fig.6(j): Distribution of rainfall during 32 nd SMW Fig.6(k): Distribution of rainfall during 33 rd SMW Fig.6: Plot of the best fitted distribution functions for maximum daily rainfall (mm) Contd.

72 Frequency Density Frequency Density Frequency Density Rainfall (mm) Rainfall (mm) Rainfall (mm) Weibull (3)(0.938,22.466,-0.678) Gamma (2)(0.677,31.916) Weibull (3)(0.668,24.647,-0.450) Fig.6(l): Distribution of rainfall during 34 th SMW Fig.6(m): Distribution of rainfall during 35 th SMW Fig.6(n): Distribution of rainfall during 36 th SMW Frequency Density Frequency Density Frequency Density Rainfall (mm) Rainfall (mm) Rainfall (mm) Weibull (3)(0.890,31.438,-0.450) Weibull (3)(0.856,24.088,-1.024) Weibull (3)(0.644,24.726,-0.200) Fig.6(o): Distribution of rainfall during 37 th SMW Fig.6(p): Distribution of rainfall during 38 th SMW Fig.6(q): Distribution of rainfall during 39 th SMW Fig.6: Plot of the best fitted distribution functions for maximum daily rainfall (mm) Contd.

73 Frequency Density Frequency Density Frequency Density Rainfall (mm) Rainfall (mm) Rainfall (mm) Weibull (3)(1.455,34.983,-5.123) Weibull (3)(0.874,24.822,-1.527) Weibull (3)(0.864,19.599,-1.018) Fig.6(r): Distribution of rainfall during 40 th SMW Fig.6(s): Distribution of rainfall during 41 st SMW Fig.6(t): Distribution of rainfall during 42 nd SMW 14 Frequency Density Rainfall (mm) Weibull (3)(0.518,9.995,-0.200) Fig. 6(u): Distribution of rainfall during 43 rd SMW Fig.6: Plot of the best fitted distribution functions for maximum daily rainfall (mm)

74 Table 4.2.1: Summary of statistics for weather parameters (Period: ) Weather parameters Mean Standard deviation CV (%) Skewness Maximum Minimum Minimum temperature Maximum temperature Relative Humidity Rainfall No. of rainy days PET SSH Soil temperature

75 4.2.2 Simple Linear Regression Analysis Simple linear regression analysis was carried out separately for each climatic factor to know the effect of weather parameters on the finger millet crop yield and the results are presented in Table All the regression coefficients pertaining to the weather parameters were found to be non- significant. The result also revealed that there was a positive influence of minimum temperature, relative humidity, rainfall and number of rainy days on the yield of finger millet whereas maximum temperature, PET, sunshine hours and soil temperature had a negative influence on the yield of finger millet. The fitting of the regression equation of individual weather parameter on the yield are depicted in Fig Multiple Linear Regression Analysis Multiple linear regression analysis was carried out initially involving all the eight weather parameters viz., minimum temperature, maximum temperature, relative humidity, rainfall, number of rainy days, PET, soil temperature and sunshine hours to know the overall effect of weather parameters on the finger millet crop yield. However, there was a problem of multicollinearity in the data set (with a value of VIF exceeding 10) showing high correlation between predictors namely, sunshine hours and maximum temperature. Hence, the predictor sunshine hours was removed from the data set and revised regression analysis was carried out with the remaining seven weather parameters. The results are presented in the Table The results showed that minimum temperature, relative humidity and number of rainy days showed positive influence on finger millet where as maximum temperature, rainfall, PET and soil temperature showed negative influence on finger millet crop yield. None of the weather parameters had significant influence on finger millet crop yield. The model F ratio (0.662) was not significant with R 2 value of The model was free from multicollinearity with all the Variance Inflation Factor (VIF) showing a value below 10.

76 Table 4.2.2: Estimates of the simple linear regression model relating the crop yield with weather parameters Weather parameters Minimum temperature Maximum temperature Relative Humidity Intercept Slope t- value SE R NS NS NS Rainfall NS No. of rainy days NS PET NS SSH NS Soil temperature NS NS Non-significant

77 39 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015 Table 4.2.3: Multiple linear regression of weather parameters on the yield of finger millet Regression Min. Temp Max. Temp RH Rainfall NS Non-significant No. of rainy days PET Soil Temp. Coefficients SE t- value NS NS NS NS NS NS NS VIF F R NS 0.367

78 Yield (qtl/ha) Minimum temperature ( 0 C) Yield(qtl/ha) Maximum temperature ( 0 C) Yield(qtl/ha) Y = X Relative Humidity (%) Yield (qtl/ha) Rainfall (mm) Fig.7: Effect of weather parameters on the finger millet crop yield (qtl/ha) Evaluation of Statistical Models for Climatic Characterization of GKVK Station 40

79 Contd Yield (qtl/ha) Yield (qtl/ha) Number of Rainy days PET (mm) Yield (qtl/ha) Sunshine hours (hrs/day) Y = X Soil Temperature ( 0 C) Fig.7: Effect of weather parameters on the finger millet crop yield (qtl/ha) 41 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015 Yield (qtl/ha)

80 4.2.4 Test for Randomness and Normality of error terms In order to test the adequacy of the model, randomness of errors was verified using runs test. The results are shown in Table From the analysis, it was confirmed that the error terms were randomly distributed (with the total number of runs being 10) for the model. The distribution of predicted residuals against the response variable (yield) indicated a horizontal pattern confirming that the data is free from model defects (as shown in Fig.8). Further, the normality of error terms was checked by using normal probability plot (Fig.9). The plot indicated near normal type with points scattered on either side of the line. Table 4.2.4: Observed and predicted yield of finger millet crop over a period of 16 years ( ). Year Y Y ) = ( Y Y ) e i ) Run test statistic No. of runs (r) =10 NS, Critical values: r 4, r = 13 L = U Evaluation of Statistical Models for Climatic Characterization of GKVK Station 42

81 4.3 Direct and indirect effects of weather parameters on rainfall. Path coefficient analysis was carried out to know the direct and indirect effects of weather parameters on rainfall. The detailed results are presented below: Correlation between weather parameters and rainfall. Correlations between weather parameters and rainfall were useful because they can indicate a predictive relationship that can be exploited in practice. Weather parameters are believed as the main factors that influence the rainfall. Therefore, it becomes important to know the correlation between the weather parameters and rainfall. The weather parameters considered were Maximum temperature (Max temp.), Minimum temperature (Min temp.), Potential Evapotranspiration (PET), Vapor Pressure (VP), Relative Humidity (RH), Wind Speed (WS), Sunshine hours (SSH), Cloud amount (CL AMT), and Evaporation (EVP). Table represents the Correlation matrix between weather parameters and rainfall Maximum temperature Maximum temperature showed a positive significant correlation (0.3575) with evaporation while non-significant positive correlation was observed with minimum temperature (0.1654), PET (0.0856), sunshine hours (0.2325) and rainfall (0.1676). On the other hand, it showed a significant negative correlation with relative humidity ( ) and non-significant negative correlation with vapor pressure ( ), wind speed ( ) and cloud amount ( ) Minimum temperature Minimum temperature had positive non-significant correlation with maximum temperature (0.1654), vapor pressure (0.2306) and rainfall (0.3196) and positive significant correlation was not observed with any of the weather parameters. The negative significant correlation was observed with sunshine hours ( ) while it showed negative non-significant correlation with PET ( ), wind speed ( ), cloud amount ( ) and evaporation ( ) Potential Evapotranspiration (PET) Potential Evapotranspiration showed positive significant correlation with wind speed (0.8040), sunshine hours (0.5551), cloud amount (0.3677) and evaporation (0.6429) while a positive non-significant correlation was observed with maximum temperature (0.0856). PET had negative significant correlation with vapor pressure ( ), relative humidity ( ) and rainfall ( ) whereas a negative nonsignificant correlation was seen with minimum temperature ( ). 43 Bhoomika Raj, R., M.Sc. (Ag. Stat.), 2015

82 Fig. 8:Plot of PRESS residuals against predicted finger millet yield Fig.9: Normal Probability Plot for finger millet yield against weather parameters Evaluation of Statistical Models for Climatic Characterization of GKVK Station 44

INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY

INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY A PATH FOR HORIZING YOUR INNOVATIVE WORK FITTING STATISTICAL DISTRUBTIONS FOR MAXIMUM DAILY RAINFALL AT GKVK STATION K.

More information

Summary and Conclusions

Summary and Conclusions 241 Chapter 10 Summary and Conclusions Kerala is situated in the southern tip of India between 8 15 N and 12 50 N latitude and 74 50 E and 77 30 E longitude. It is popularly known as Gods own country.

More information

Probability models for weekly rainfall at Thrissur

Probability models for weekly rainfall at Thrissur Journal of Tropical Agriculture 53 (1) : 56-6, 015 56 Probability models for weekly rainfall at Thrissur C. Laly John * and B. Ajithkumar *Department of Agricultural Statistics, College of Horticulture,

More information

3. HYDROMETEROLOGY. 3.1 Introduction. 3.2 Hydro-meteorological Aspect. 3.3 Rain Gauge Stations

3. 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 information

Our climate system is based on the location of hot and cold air mass regions and the atmospheric circulation created by trade winds and westerlies.

Our climate system is based on the location of hot and cold air mass regions and the atmospheric circulation created by trade winds and westerlies. CLIMATE REGIONS Have you ever wondered why one area of the world is a desert, another a grassland, and another a rainforest? Or have you wondered why are there different types of forests and deserts with

More information

Lecture Outlines PowerPoint. Chapter 20 Earth Science 11e Tarbuck/Lutgens

Lecture Outlines PowerPoint. Chapter 20 Earth Science 11e Tarbuck/Lutgens Lecture Outlines PowerPoint Chapter 20 Earth Science 11e Tarbuck/Lutgens 2006 Pearson Prentice Hall This work is protected by United States copyright laws and is provided solely for the use of instructors

More information

Climate Chapter 19. Earth Science, 10e. Stan Hatfield and Ken Pinzke Southwestern Illinois College

Climate Chapter 19. Earth Science, 10e. Stan Hatfield and Ken Pinzke Southwestern Illinois College Climate Chapter 19 Earth Science, 10e Stan Hatfield and Ken Pinzke Southwestern Illinois College The climate system A. Climate is an aggregate of weather B. Involves the exchanges of energy and moisture

More information

Analysis of Rainfall and Other Weather Parameters under Climatic Variability of Parbhani ( )

Analysis 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 information

Rainfall is the major source of water for

Rainfall is the major source of water for RESEARCH PAPER: Assessment of occurrence and frequency of drought using rainfall data in Coimbatore, India M. MANIKANDAN AND D.TAMILMANI Asian Journal of Environmental Science December, 2011 Vol. 6 Issue

More information

Energy Systems, Structures and Processes Essential Standard: Analyze patterns of global climate change over time Learning Objective: Differentiate

Energy Systems, Structures and Processes Essential Standard: Analyze patterns of global climate change over time Learning Objective: Differentiate Energy Systems, Structures and Processes Essential Standard: Analyze patterns of global climate change over time Learning Objective: Differentiate between weather and climate Global Climate Focus Question

More information

The Global Scope of Climate. The Global Scope of Climate. Keys to Climate. Chapter 8

The Global Scope of Climate. The Global Scope of Climate. Keys to Climate. Chapter 8 The Global Scope of Climate Chapter 8 The Global Scope of Climate In its most general sense, climate is the average weather of a region, but except where conditions change very little during the course

More information

Key Finding: Long Term Trend During 2014: Rain in Indian Tradition Measuring Rain

Key Finding: Long Term Trend During 2014: Rain in Indian Tradition Measuring Rain Chapter 34 RAINFALL Key Finding: Long Term Trend Despite of theories suggesting increase in rainfall in Asian Region due to global warming, no significant trend has been observed at all India level (confirmed

More information

BANGALORE URBAN SPS / CPS / CSS ANEKAL. Sl. No. Plan Non plan Plan Non plan Plan Non plan Plan Non plan Plan Non plan. Released by

BANGALORE URBAN SPS / CPS / CSS ANEKAL. Sl. No. Plan Non plan Plan Non plan Plan Non plan Plan Non plan Plan Non plan. Released by 1 2202-00-101-0-19 Sarva Shiksha Abhiyana Society ಸವ ಣ ಅ ನ ಸ SPS FD 362.00 0.00 886.54 0.00 856.00 0.00 250.00 0.00 2354.54 0.00 2 2202-00-101-0-63 Elementary Schools GIA ಥ ಕ ಗ SPS FD 25.00 171.00 56.00

More information

Weekly Rainfall Analysis and Markov Chain Model Probability of Dry and Wet Weeks at Varanasi in Uttar Pradesh

Weekly Rainfall Analysis and Markov Chain Model Probability of Dry and Wet Weeks at Varanasi in Uttar Pradesh 885 Environment & Ecology 32 (3) : 885 890, July September 2014 Website: environmentandecology.com ISSN 0970-0420 Weekly Rainfall Analysis and Markov Chain Model Probability of Dry and Wet Weeks at Varanasi

More information

Geographical location and climatic condition of the

Geographical location and climatic condition of the Geographical location and climatic condition of the study sites North eastern region of India is comprised of eight states namely; Arunachal Pradesh, Assam, Manipur, Meghalaya, Mizoram, Nagaland, Sikkim

More information

DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION

DROUGHT 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 information

Advanced Hydrology. (Web course)

Advanced Hydrology. (Web course) Advanced Hydrology (Web course) Subhankar Karmakar Assistant Professor Centre for Environmental Science and Engineering (CESE) Indian Institute of Technology Bombay Powai, Mumbai 400 076 Email: skarmakar@iitb.ac.in

More information

Chapter-1 Introduction

Chapter-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 information

Best Fit Probability Distributions for Monthly Radiosonde Weather Data

Best Fit Probability Distributions for Monthly Radiosonde Weather Data Best Fit Probability Distributions for Monthly Radiosonde Weather Data Athulya P. S 1 and K. C James 2 1 M.Tech III Semester, 2 Professor Department of statistics Cochin University of Science and Technology

More information

Study of Hydrometeorology in a Hard Rock Terrain, Kadirischist Belt Area, Anantapur District, Andhra Pradesh

Study 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 information

Dry spell analysis for effective water management planning

Dry spell analysis for effective water management planning www.ijaser.com 2012 by the authors Licensee IJASER- Under Creative Commons License 3.0 editorial@ijaser.com Research article ISSN 2277 9442 Dry spell analysis for effective water management planning Kandasamy

More information

Study of Changes in Climate Parameters at Regional Level: Indian Scenarios

Study of Changes in Climate Parameters at Regional Level: Indian Scenarios Study of Changes in Climate Parameters at Regional Level: Indian Scenarios S K Dash Centre for Atmospheric Sciences Indian Institute of Technology Delhi Climate Change and Animal Populations - The golden

More information

Frequency analysis of rainfall deviation in Dharmapuri district in Tamil Nadu

Frequency analysis of rainfall deviation in Dharmapuri district in Tamil Nadu Indian Journal of Geo Marine Sciences Vol.46 (8), August 217, pp. 171-1714 Frequency analysis of rainfall deviation in Dharmapuri district in Tamil Nadu V. Rajendran 1* & R. Venkatasubramani 2 1 Department

More information

Effect of rainfall and temperature on rice yield in Puri district of Odisha in India

Effect of rainfall and temperature on rice yield in Puri district of Odisha in India 2018; 7(4): 899-903 ISSN (E): 2277-7695 ISSN (P): 2349-8242 NAAS Rating: 5.03 TPI 2018; 7(4): 899-903 2018 TPI www.thepharmajournal.com Received: 05-02-2018 Accepted: 08-03-2018 A Baliarsingh A Nanda AKB

More information

Rainfall variation and frequency analysis study in Dharmapuri district, India

Rainfall variation and frequency analysis study in Dharmapuri district, India Indian Journal of Geo Marine Sciences Vol. 45 (11), November 216, pp. 156-1565 Rainfall variation and frequency analysis study in Dharmapuri district, India V. Rajendran 1*, R. Venkatasubramani 2 & G.

More information

WHAT CAN MAPS TELL US ABOUT THE GEOGRAPHY OF ANCIENT GREECE? MAP TYPE 1: CLIMATE MAPS

WHAT CAN MAPS TELL US ABOUT THE GEOGRAPHY OF ANCIENT GREECE? MAP TYPE 1: CLIMATE MAPS WHAT CAN MAPS TELL US ABOUT THE GEOGRAPHY OF ANCIENT GREECE? MAP TYPE 1: CLIMATE MAPS MAP TYPE 2: PHYSICAL AND/OR TOPOGRAPHICAL MAPS MAP TYPE 3: POLITICAL MAPS TYPE 4: RESOURCE & TRADE MAPS Descriptions

More information

SEASONAL CLIMATE OUTLOOK VALID FOR JULY-AUGUST- SEPTEMBER 2013 IN WEST AFRICA, CHAD AND CAMEROON

SEASONAL CLIMATE OUTLOOK VALID FOR JULY-AUGUST- SEPTEMBER 2013 IN WEST AFRICA, CHAD AND CAMEROON SEASONAL CLIMATE OUTLOOK VALID FOR JULY-AUGUST- SEPTEMBER 2013 IN WEST AFRICA, CHAD AND CAMEROON May 29, 2013 ABUJA-Federal Republic of Nigeria 1 EXECUTIVE SUMMARY Given the current Sea Surface and sub-surface

More information

Projected Change in Climate Under A2 Scenario in Dal Lake Catchment Area of Srinagar City in Jammu and Kashmir

Projected 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 information

Climate Change Impact on Air Temperature, Daily Temperature Range, Growing Degree Days, and Spring and Fall Frost Dates In Nebraska

Climate Change Impact on Air Temperature, Daily Temperature Range, Growing Degree Days, and Spring and Fall Frost Dates In Nebraska EXTENSION Know how. Know now. Climate Change Impact on Air Temperature, Daily Temperature Range, Growing Degree Days, and Spring and Fall Frost Dates In Nebraska EC715 Kari E. Skaggs, Research Associate

More information

World Geography Chapter 3

World Geography Chapter 3 World Geography Chapter 3 Section 1 A. Introduction a. Weather b. Climate c. Both weather and climate are influenced by i. direct sunlight. ii. iii. iv. the features of the earth s surface. B. The Greenhouse

More information

Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC

Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC This threat overview relies on projections of future climate change in the Mekong Basin for the period 2045-2069 compared to a baseline of 1980-2005.

More information

NAKSHATRA BASED RAINFALL ANALYSIS AND ITS IMPACT ON CROPS DURING MONSOON SEASON AT MANDYA DISTRICT

NAKSHATRA BASED RAINFALL ANALYSIS AND ITS IMPACT ON CROPS DURING MONSOON SEASON AT MANDYA DISTRICT International Journal of Current Agricultural Research, Vol. 4, No. 1, PP. 139-144 October 215 http://www.journalagri.com Full Length Research Article NAKSHATRA BASED RAINFALL ANALYSIS AND ITS IMPACT ON

More information

Mozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1

Mozambique. 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 information

Climate Variability in South Asia

Climate 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 information

By: J Malherbe, R Kuschke

By: J Malherbe, R Kuschke 2015-10-27 By: J Malherbe, R Kuschke Contents Summary...2 Overview of expected conditions over South Africa during the next few days...3 Significant weather events (27 October 2 November)...3 Conditions

More information

Rainfall variation and frequency analysis study of Salem district Tamil Nadu

Rainfall variation and frequency analysis study of Salem district Tamil Nadu Indian Journal of Geo Marine Sciences Vol. 46 (1), January 217, pp. 213-218 Rainfall variation and frequency analysis study of Salem district Tamil Nadu Arulmozhi.S 1* & Dr. Prince Arulraj.G 2 1 Department

More information

RESEARCH NOTE Changing Dew Patterns in Anantapur District, Andhra Pradesh: A Generalistic Observation INTRODUCTION

RESEARCH NOTE Changing Dew Patterns in Anantapur District, Andhra Pradesh: A Generalistic Observation INTRODUCTION RESEARCH NOTE Changing Dew Patterns in Anantapur District, Andhra Pradesh: A Generalistic Observation R.V. Rama Mohan* I INTRODUCTION Ind. Jn. of Agri. Econ. Vol.67, No.4, Oct.-Dec. 212 Anantapur district

More information

Elements of weather and climate Temperature Range of temperature Seasonal temperature pattern Rainfall

Elements of weather and climate Temperature Range of temperature Seasonal temperature pattern Rainfall Climate Earth Science Chapter 20 Pages 560-573 Elements of weather and climate Temperature Range of temperature Seasonal temperature pattern Rainfall Overall rainfall Seasonal distribution of rainfall

More information

The analysis of climatic variability/weather trends (past and future) in eastern U.P.

The analysis of climatic variability/weather trends (past and future) in eastern U.P. 2018; SP1: 1092-1096 E-ISSN: 2278-4136 P-ISSN: 2349-8234 JPP 2018; SP1: 1092-1096 Ajit Singh AK Singh Department of Agriculture Meteorology, Narendra Deva Technology Narendra Nagar, Kumarganj, Faizabad,

More information

GLOBAL CLIMATES FOCUS

GLOBAL CLIMATES FOCUS which you will learn more about in Chapter 6. Refer to the climate map and chart on pages 28-29 as you read the rest of this chapter. FOCUS GLOBAL CLIMATES What are the major influences on climate? Where

More information

KUALA LUMPUR MONSOON ACTIVITY CENT

KUALA LUMPUR MONSOON ACTIVITY CENT T KUALA LUMPUR MONSOON ACTIVITY CENT 2 ALAYSIAN METEOROLOGICAL http://www.met.gov.my DEPARTMENT MINISTRY OF SCIENCE. TECHNOLOGY AND INNOVATIO Introduction Atmospheric and oceanic conditions over the tropical

More information

SEASONAL RAINFALL FORECAST FOR ZIMBABWE. 28 August 2017 THE ZIMBABWE NATIONAL CLIMATE OUTLOOK FORUM

SEASONAL RAINFALL FORECAST FOR ZIMBABWE. 28 August 2017 THE ZIMBABWE NATIONAL CLIMATE OUTLOOK FORUM 2017-18 SEASONAL RAINFALL FORECAST FOR ZIMBABWE METEOROLOGICAL SERVICES DEPARTMENT 28 August 2017 THE ZIMBABWE NATIONAL CLIMATE OUTLOOK FORUM Introduction The Meteorological Services Department of Zimbabwe

More information

Tropical Moist Rainforest

Tropical Moist Rainforest Tropical or Lowlatitude Climates: Controlled by equatorial tropical air masses Tropical Moist Rainforest Rainfall is heavy in all months - more than 250 cm. (100 in.). Common temperatures of 27 C (80 F)

More information

Keys to Climate Climate Classification Low Latitude Climates Midlatitude Climates High Latitude Climates Highland Climates Our Changing Climate

Keys to Climate Climate Classification Low Latitude Climates Midlatitude Climates High Latitude Climates Highland Climates Our Changing Climate Climate Global Climates Keys to Climate Climate Classification Low Latitude Climates Midlatitude Climates High Latitude Climates Highland Climates Our Changing Climate Keys to Climate Climate the average

More information

Development of regression models in ber genotypes under the agroclimatic conditions of south-western region of Punjab, India

Development of regression models in ber genotypes under the agroclimatic conditions of south-western region of Punjab, India Indian J. Agric. Res., 49 (3) 2015: 260-264 Print ISSN:0367-8245 / Online ISSN:0976-058X AGRICULTURAL RESEARCH COMMUNICATION CENTRE www.arccjournals.com/www.ijarjournal.com Development of regression models

More information

Weather and Climate of the Rogue Valley By Gregory V. Jones, Ph.D., Southern Oregon University

Weather and Climate of the Rogue Valley By Gregory V. Jones, Ph.D., Southern Oregon University Weather and Climate of the Rogue Valley By Gregory V. Jones, Ph.D., Southern Oregon University The Rogue Valley region is one of many intermountain valley areas along the west coast of the United States.

More information

1 Ministry of Earth Sciences, Lodi Road, New Delhi India Meteorological Department, Lodi Road, New Delhi

1 Ministry of Earth Sciences, Lodi Road, New Delhi India Meteorological Department, Lodi Road, New Delhi Trends in Extreme Temperature Events over India during 1969-12 A. K. JASWAL, AJIT TYAGI 1 and S. C. BHAN 2 India Meteorological Department, Shivajinagar, Pune - 4105 1 Ministry of Earth Sciences, Lodi

More information

Long Range Forecast Update for 2014 Southwest Monsoon Rainfall

Long Range Forecast Update for 2014 Southwest Monsoon Rainfall Earth System Science Organization (ESSO) Ministry of Earth Sciences (MoES) India Meteorological Department PRESS RELEASE New Delhi, 9 June 2014 Long Update for 2014 Southwest Monsoon Rainfall HIGHLIGHTS

More information

16 Global Climate. Learning Goals. Summary. After studying this chapter, students should be able to:

16 Global Climate. Learning Goals. Summary. After studying this chapter, students should be able to: 16 Global Climate Learning Goals After studying this chapter, students should be able to: 1. associate the world s six major vegetation biomes to climate (pp. 406 408); 2. describe methods for classifying

More information

Seasonal Climate Watch January to May 2016

Seasonal Climate Watch January to May 2016 Seasonal Climate Watch January to May 2016 Date: Dec 17, 2015 1. Advisory Most models are showing the continuation of a strong El-Niño episode towards the latesummer season with the expectation to start

More information

Statistical Analysis of Temperature and Rainfall Trend in Raipur District of Chhattisgarh

Statistical Analysis of Temperature and Rainfall Trend in Raipur District of Chhattisgarh Current World Environment Vol. 10(1), 305-312 (2015) Statistical Analysis of Temperature and Rainfall Trend in Raipur District of Chhattisgarh R. Khavse*, R. Deshmukh, N. Manikandan, J. L Chaudhary and

More information

Rainfall Analysis in Mumbai using Gumbel s Extreme Value Distribution Model

Rainfall Analysis in Mumbai using Gumbel s Extreme Value Distribution Model Rainfall Analysis in Mumbai using Gumbel s Extreme Value Distribution Model S. Sangita Mishra 1, M. Solapure 2, R. Satpathy 3, K. Phadtare 4, S. Pawar 5 1 Assistant Professor, 2,3,4,5 Undergraduate students

More information

Climate and the Atmosphere

Climate and the Atmosphere Climate and Biomes Climate Objectives: Understand how weather is affected by: 1. Variations in the amount of incoming solar radiation 2. The earth s annual path around the sun 3. The earth s daily rotation

More information

Suppressed 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 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 information

Long-Term Trend of Summer Rainfall at Selected Stations in the Republic of Korea

Long-Term Trend of Summer Rainfall at Selected Stations in the Republic of Korea Long-Term Trend of Summer Rainfall at Selected Stations in the Republic of Korea Il-Kon Kim Professor, Department of Region Information Rafique Ahmed Professor, Geography and Earth Science Silla University

More information

Weather and Climate Summary and Forecast Winter

Weather and Climate Summary and Forecast Winter Weather and Climate Summary and Forecast Winter 2016-17 Gregory V. Jones Southern Oregon University February 7, 2017 What a difference from last year at this time. Temperatures in January and February

More information

PART II. Physical Landscape Chapters 2 5 CLIMATE CLIMATE STUDYING CLIMATE R E M I N D E R S. PART II: People and their Physical Environment 10/26/2017

PART II. Physical Landscape Chapters 2 5 CLIMATE CLIMATE STUDYING CLIMATE R E M I N D E R S. PART II: People and their Physical Environment 10/26/2017 R E M I N D E R S Two required essays are due by Nov. 13, 2017. (A third may be used for extra credit in place of a Think Geographically essay.) ESSAY TOPIS (choose any two): ontributions of a noted geographer,

More information

Seasonal Climate Outlook for South Asia (June to September) Issued in May 2014

Seasonal Climate Outlook for South Asia (June to September) Issued in May 2014 Ministry of Earth Sciences Earth System Science Organization India Meteorological Department WMO Regional Climate Centre (Demonstration Phase) Pune, India Seasonal Climate Outlook for South Asia (June

More information

A Study on Storm Water Drainage System of Annanagara and Ashokanagara of Shimoga City Karnataka India

A Study on Storm Water Drainage System of Annanagara and Ashokanagara of Shimoga City Karnataka India A Study on Storm Water Drainage System of Annanagara and Ashokanagara of Shimoga City Karnataka India Veena D Savanth Lecturer, Department of Civil Engineering Jain Institute of Technology, Davanagere,

More information

Prentice Hall EARTH SCIENCE

Prentice Hall EARTH SCIENCE Prentice Hall EARTH SCIENCE Tarbuck Lutgens Chapter 21 Climate 21.1 Factors That Affect Climate Factors That Affect Climate Latitude As latitude increases, the intensity of solar energy decreases. The

More information

Chapter- Three GEOMORPHOLOGY AND RAINFALL PATTERN

Chapter- Three GEOMORPHOLOGY AND RAINFALL PATTERN Chapter- Three GEOMORPHOLOGY AND RAINFALL PATTERN 3.1 INTRODUCTION Geomorphology and its slope as a basic and applied science in general and as tool for searching groundwater resources in efficient geomorphological

More information

1 What Is Climate? TAKE A LOOK 2. Explain Why do areas near the equator tend to have high temperatures?

1 What Is Climate? TAKE A LOOK 2. Explain Why do areas near the equator tend to have high temperatures? CHAPTER 17 1 What Is Climate? SECTION Climate BEFORE YOU READ After you read this section, you should be able to answer these questions: What is climate? What factors affect climate? How do climates differ

More information

Geography 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. 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 information

CLIMATE DEFINED CLIMATE INFLUENCES STUDYING CLIMATE R E M I N D E R S. PART II: People and their 10/19/2018

CLIMATE DEFINED CLIMATE INFLUENCES STUDYING CLIMATE R E M I N D E R S. PART II: People and their 10/19/2018 R E M I N D E R S Two required essays are due by Oct. 30, 2018. (A third may be used for extra credit in place of a Think Geographically essay.) ESSAY TOPIS (choose any two): ontributions of a noted geographer,

More information

UNITED STATES AND SOUTH AMERICA OUTLOOK (FULL REPORT) Wednesday, April 18, 2018

UNITED STATES AND SOUTH AMERICA OUTLOOK (FULL REPORT) Wednesday, April 18, 2018 T-storm Weather Summary Coolness continues over the next week in much of the central U.S., most likely producing the coldest April since 1907 in the Corn Belt, but followed by seasonable to mild weather

More information

Trend and Variability Analysis and Forecasting of Wind-Speed in Bangladesh

Trend and Variability Analysis and Forecasting of Wind-Speed in Bangladesh J. Environ. Sci. & Natural Resources, 5(): 97-07, 0 ISSN 999-736 Trend and Variability Analysis and Forecasting of Wind-Speed in Bangladesh J. A. Syeda Department of Statistics, Hajee Mohammad Danesh Science

More information

Climate. Earth Science Chapter 20 Pages

Climate. Earth Science Chapter 20 Pages Academic Excellence Showcase Thursday, May 27 Two points credit for each talk you attend in Earth science 8:30 to 10:30 AM, in NS 101 General Science video productions: Energy 9:00 to 11:00 AM, in NS 218

More information

1 What Is Climate? TAKE A LOOK 2. Explain Why do areas near the equator tend to have high temperatures?

1 What Is Climate? TAKE A LOOK 2. Explain Why do areas near the equator tend to have high temperatures? CHAPTER 17 1 What Is Climate? SECTION Climate BEFORE YOU READ After you read this section, you should be able to answer these questions: What is climate? What factors affect climate? How do climates differ

More information

DROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE

DROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE DROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE K. Prathumchai, Kiyoshi Honda, Kaew Nualchawee Asian Centre for Research on Remote Sensing STAR Program, Asian Institute

More information

November 2018 Weather Summary West Central Research and Outreach Center Morris, MN

November 2018 Weather Summary West Central Research and Outreach Center Morris, MN November 2018 Weather Summary Lower than normal temperatures occurred for the second month. The mean temperature for November was 22.7 F, which is 7.2 F below the average of 29.9 F (1886-2017). This November

More information

Analysis 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. 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 information

CERTAIN ASPECTS OF THE NUMERICAL SIMULATION OF CIRCULATION IN THE BAY OF BENGAL AND THE ARABIAN SEA

CERTAIN ASPECTS OF THE NUMERICAL SIMULATION OF CIRCULATION IN THE BAY OF BENGAL AND THE ARABIAN SEA CERTAIN ASPECTS OF THE NUMERICAL SIMULATION OF CIRCULATION IN THE BAY OF BENGAL AND THE ARABIAN SEA by INDU JAIN A THESIS SUBMITTED IN FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY

More information

CHAPTER 1: INTRODUCTION

CHAPTER 1: INTRODUCTION CHAPTER 1: INTRODUCTION There is now unequivocal evidence from direct observations of a warming of the climate system (IPCC, 2007). Despite remaining uncertainties, it is now clear that the upward trend

More information

Chapter 2 Variability and Long-Term Changes in Surface Air Temperatures Over the Indian Subcontinent

Chapter 2 Variability and Long-Term Changes in Surface Air Temperatures Over the Indian Subcontinent Chapter 2 Variability and Long-Term Changes in Surface Air Temperatures Over the Indian Subcontinent A.K. Srivastava, D.R. Kothawale and M.N. Rajeevan 1 Introduction Surface air temperature is one of the

More information

Environmental and Earth Sciences Research Journal Vol. 5, No. 3, September, 2018, pp Journal homepage:

Environmental and Earth Sciences Research Journal Vol. 5, No. 3, September, 2018, pp Journal homepage: Environmental and Earth Sciences Research Journal Vol. 5, No. 3, September, 2018, pp. 74-78 Journal homepage: http://iieta.org/journals/eesrj Analysis of rainfall variation over northern parts of Nigeria

More information

variation analysis of Vattamalaikarai sub basin, Tamil Nadu, India.

variation analysis of Vattamalaikarai sub basin, Tamil Nadu, India. A Study on the Rainfall Changes of Orathanadu Taluk Thanjavur District Tamil Nadu Using GIS P. Sujatha 1, Dr. R Baskaran 2 1 Research Scholar, Department of Industries and Earth Sciences, Tamil University,

More information

Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014

Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014 Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014 Short Term Drought Map: Short-term (

More information

Analysis of Meteorological drought condition for Bijapur region in the lower Bhima basin, India

Analysis of Meteorological drought condition for Bijapur region in the lower Bhima basin, India Analysis of Meteorological drought condition for Bijapur region in the lower Bhima basin, India Mamatha.K PG Student Department of WLM branch VTU, Belagavi Dr. Nagaraj Patil Professor and Head of the Department

More information

World geography 3200/3202 Unit 2 review

World geography 3200/3202 Unit 2 review World geography 3200/3202 Unit 2 review 1. Does this statement use the terms revolve & rotate correctly? "Saturn revolves on its axis while several moons rotate around it." 2. Does this statement use the

More information

UPDATE OF REGIONAL WEATHER AND SMOKE HAZE (December 2017)

UPDATE OF REGIONAL WEATHER AND SMOKE HAZE (December 2017) UPDATE OF REGIONAL WEATHER AND SMOKE HAZE (December 2017) 1. Review of Regional Weather Conditions for November 2017 1.1 In November 2017, Southeast Asia experienced inter-monsoon conditions in the first

More information

Weather and Climate Summary and Forecast Summer 2017

Weather and Climate Summary and Forecast Summer 2017 Weather and Climate Summary and Forecast Summer 2017 Gregory V. Jones Southern Oregon University August 4, 2017 July largely held true to forecast, although it ended with the start of one of the most extreme

More information

Water Stress, Droughts under Changing Climate

Water Stress, Droughts under Changing Climate Water Stress, Droughts under Changing Climate Professor A.K.M. Saiful Islam Institute of Water and Flood Management Bangladesh University of Engineering and Technology (BUET) Outline of the presentation

More information

J8.4 TRENDS OF U.S. SNOWFALL AND SNOW COVER IN A WARMING WORLD,

J8.4 TRENDS OF U.S. SNOWFALL AND SNOW COVER IN A WARMING WORLD, J8.4 TRENDS OF U.S. SNOWFALL AND SNOW COVER IN A WARMING WORLD, 1948-2008 Richard R. Heim Jr. * NOAA National Climatic Data Center, Asheville, North Carolina 1. Introduction The Intergovernmental Panel

More information

MARKOV CHAIN MODEL FOR PROBABILITY OF DRY, WET DAYS AND STATISTICAL ANALISIS OF DAILY RAINFALL IN SOME CLIMATIC ZONE OF IRAN

MARKOV CHAIN MODEL FOR PROBABILITY OF DRY, WET DAYS AND STATISTICAL ANALISIS OF DAILY RAINFALL IN SOME CLIMATIC ZONE OF IRAN 30T r30t MARKOV CHAIN MODEL FOR PROBABILITY OF DRY, WET DAYS AND STATISTICAL ANALISIS OF DAILY RAINFALL IN SOME CLIMATIC ZONE OF IRAN N. SHAHRAKI 1, B. BAKHTIARI 2, M. M. AHMADI 2 ABSTRACT- Markov chain

More information

CLIMATE READY BOSTON. Climate Projections Consensus ADAPTED FROM THE BOSTON RESEARCH ADVISORY GROUP REPORT MAY 2016

CLIMATE READY BOSTON. Climate Projections Consensus ADAPTED FROM THE BOSTON RESEARCH ADVISORY GROUP REPORT MAY 2016 CLIMATE READY BOSTON Sasaki Steering Committee Meeting, March 28 nd, 2016 Climate Projections Consensus ADAPTED FROM THE BOSTON RESEARCH ADVISORY GROUP REPORT MAY 2016 WHAT S IN STORE FOR BOSTON S CLIMATE?

More information

National Wildland Significant Fire Potential Outlook

National Wildland Significant Fire Potential Outlook National Wildland Significant Fire Potential Outlook National Interagency Fire Center Predictive Services Issued: September, 2007 Wildland Fire Outlook September through December 2007 Significant fire

More information

Prentice Hall EARTH SCIENCE

Prentice Hall EARTH SCIENCE Prentice Hall EARTH SCIENCE Tarbuck Lutgens Chapter 21 Climate 21.1 Factors That Affect Climate Factors That Affect Climate Latitude As latitude increases, the intensity of solar energy decreases. The

More information

Chapter 3 Section 3 World Climate Regions In-Depth Resources: Unit 1

Chapter 3 Section 3 World Climate Regions In-Depth Resources: Unit 1 Guided Reading A. Determining Cause and Effect Use the organizer below to show the two most important causes of climate. 1. 2. Climate B. Making Comparisons Use the chart below to compare the different

More information

DETECTION OF TREND IN RAINFALL DATA: A CASE STUDY OF SANGLI DISTRICT

DETECTION OF TREND IN RAINFALL DATA: A CASE STUDY OF SANGLI DISTRICT ORIGINAL ARTICLE DETECTION OF TREND IN RAINFALL DATA: A CASE STUDY OF SANGLI DISTRICT M. K. Patil 1 and D. N. Kalange 2 1 Associate Professor, Padmabhushan Vasantraodada Patil Mahavidyalaya, Kavathe- Mahankal,

More information

Climate outlook, longer term assessment and regional implications. What s Ahead for Agriculture: How to Keep One of Our Key Industries Sustainable

Climate outlook, longer term assessment and regional implications. What s Ahead for Agriculture: How to Keep One of Our Key Industries Sustainable Climate outlook, longer term assessment and regional implications What s Ahead for Agriculture: How to Keep One of Our Key Industries Sustainable Bureau of Meteorology presented by Dr Jeff Sabburg Business

More information

Current and future climate of the Cook Islands. Pacific-Australia Climate Change Science and Adaptation Planning Program

Current and future climate of the Cook Islands. Pacific-Australia Climate Change Science and Adaptation Planning Program Pacific-Australia Climate Change Science and Adaptation Planning Program Penrhyn Pukapuka Nassau Suwarrow Rakahanga Manihiki N o r t h e r n C o o k I s l a nds S o u t h e Palmerston r n C o o k I s l

More information

CLIMATE CLIMATE STUDYING CLIMATE. PART II: People and their 3/20/2018. GEOG 101 Part II People and their Physical Environment R E M I N D E R S

CLIMATE CLIMATE STUDYING CLIMATE. PART II: People and their 3/20/2018. GEOG 101 Part II People and their Physical Environment R E M I N D E R S Two required essays are due by March 27, 2018. (A third may be used for extra credit in place of a Think Geographically essay.) ESSAY TOPIS (choose any two): ontributions of a noted geographer, earth scientist

More information

Evidence 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 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 information

Meteorology. Chapter 15 Worksheet 1

Meteorology. Chapter 15 Worksheet 1 Chapter 15 Worksheet 1 Meteorology Name: Circle the letter that corresponds to the correct answer 1) The Tropic of Cancer and the Arctic Circle are examples of locations determined by: a) measuring systems.

More information

1. Evaluation of Flow Regime in the Upper Reaches of Streams Using the Stochastic Flow Duration Curve

1. Evaluation of Flow Regime in the Upper Reaches of Streams Using the Stochastic Flow Duration Curve 1. Evaluation of Flow Regime in the Upper Reaches of Streams Using the Stochastic Flow Duration Curve Hironobu SUGIYAMA 1 ABSTRACT A stochastic estimation of drought evaluation in the upper reaches of

More information

Sudan Seasonal Monitor

Sudan Seasonal Monitor Sudan Seasonal Monitor Sudan Meteorological Authority Federal Ministry of Agriculture and Forestry Issue 5 August 2010 Summary Advanced position of ITCZ during July to most north of Sudan emerged wide

More information

THE CANADIAN CENTRE FOR CLIMATE MODELLING AND ANALYSIS

THE CANADIAN CENTRE FOR CLIMATE MODELLING AND ANALYSIS THE CANADIAN CENTRE FOR CLIMATE MODELLING AND ANALYSIS As Canada s climate changes, and weather patterns shift, Canadian climate models provide guidance in an uncertain future. CANADA S CLIMATE IS CHANGING

More information

Climate Classification Chapter 7

Climate Classification Chapter 7 Climate Classification Chapter 7 Climate Systems Earth is extremely diverse No two places exactly the same Similarities between places allow grouping into regions Climates influence ecosystems Why do we

More information

Meteorology. Circle the letter that corresponds to the correct answer

Meteorology. Circle the letter that corresponds to the correct answer Chapter 3 Worksheet 1 Meteorology Name: Circle the letter that corresponds to the correct answer 1) If the maximum temperature for a particular day is 26 C and the minimum temperature is 14 C, the daily

More information

Climate Classification

Climate Classification Chapter 15: World Climates The Atmosphere: An Introduction to Meteorology, 12 th Lutgens Tarbuck Lectures by: Heather Gallacher, Cleveland State University Climate Classification Köppen classification:

More information