Hydro-meteorological related disasters, particularly from climate extremes are responsible for about 90% of all global natural disasters  . Changes in many extreme weather and climate events have been observed since 1950    . The frequency and magnitude of climate extremes including droughts, floods and heat waves have been found to increase on global scale, and have significantly contributed to the loss of life and properties and severe infrastructure damages in various parts of the world    . Recent IPCC reports    confirm that it is very likely that there has been an overall decrease in the number of cold days and cold nights, and an overall increase in the number of warm days and warm nights. The frequency of heavy precipitation is likely to increase in the 21st century over many areas of the globe, while there is medium confidence that droughts will intensify in some areas due to reduced precipitation and/or increased evapotranspiration  .
Over the last few decades most parts of Tanzania have experienced an increased incidence of extreme climatic events, particularly floods and droughts that have been associated with severe socio-economic and ecological implications. The most recent severe floods include those of 2006, 2009, 2010, 2011, 2012, 2014, 2016 and 2017 that affected many parts of the country associated with many fatalities including loss of life and severe infrastructural destructions   . The country has also experienced severe and recurrent droughts, which have caused devastating impacts, particularly in agriculture, water, energy and livestock sectors, and thus paralyzing socio-economic activities and development. The most recent devastating droughts include those of 2003, 2005, 2011, 2014 and 2016  .
Climate extremes therefore poses a serious threat to the achievement of various socio-economic development strategies and plans including the achievement of national development vision (Vision 2025), National Strategy for Growth and Reduction of Poverty and the Sustainable Development Goals (SDG). Goal number one of the SDG aims at ending poverty in all its forms by 2030, while goal number 2 aims at ending hunger, achieving food security and improving nutrition and promoting sustainable development. The Tanzania development vision (Vision 2025) also set a target of graduating from a least developed country to a middle income country and transforming the economy from a low productivity agricultural economy to a semi-industrialized led by modernized and highly productive agricultural activities. Among the goals of the Vision 2025 is attaining high quality livelihoods through ensuring food self-sufficiency and food security, and ending abject poverty by 2025. Successful achievement of these goals requires effective plans and strategies to enhance and sustain production and productivity, particularly in agricultural and livestock sectors in the context of increasing frequency and magnitude of climate extremes.
For most of the developing countries including Tanzania of which socio-economic development and livelihoods are strongly linked to their dependencies on natural resources and rain-fed agriculture, Initiatives and efforts towards poverty eradication and reaching zero hunger are significantly affected by the trends and patterns of climate extremes. Unfortunately, extreme events are projected to increase and will continue to have greater impacts, particularly on climate sensitive sectors such as water, health, forestry, agriculture and food security.
Understanding the spatial and temporal patterns of rainfall and temperature related extremes is therefore pre-requisite and a critical inputs for effective designing and implementing coping and adaptation options to the impacts of climate change. This study provides a current understanding of the countrywide trends of rainfall and temperature related extreme indices. The trends in the percentage change in rainfall, anomaly, percentage change in number of warm days and warm nights, percentage change in number of cold days and cold nights, and trends in Standardized Precipitation Index (SPI) are analyzed and presented.
2. Description of the Study Area
The United Republic of Tanzania is among the East African countries lying in the south of the equator between 1 - 12˚S and 29 - 41˚E. The country has a total area of 945,087 km2 of which 883,749 km2 is covered by land area, and 59,050 km2 is covered by inland water bodies and part of the Indian Ocean. The country is characterized by diverse and complex topographical features extending from a narrow coastal belt of the western Indian Ocean with sandy beaches to an extensive plateau with altitude ranging from 1000 to 2000 m above mean sea level. The plateau is fringed by narrow belts of highlands, including Mount Meru with an altitude of (4566 m), and the highest mountain in Africa and the second highest mountain in the world: Mount Kilimanjaro with an altitude of (5895 m). Tanzania has several fresh water bodies, including the largest Lake in Africa; Lake Victoria, the longest and deepest Lake in Africa; Lake Tanganyika and Lake Nyasa     . The country also has many large rivers, draining into nine drainage basins. The major rivers include Rufiji, Kagera, Mara, Ruaha, Pangani, Ruvuma and Malagarasi. The country experience Unimodal and Bimodal rainfall regimes. The northern parts of the country including areas around Lake Victoria basin, Northern Coast and areas around Kilimanjaro Mountain experience bimodal rainfall regime, the first occurs in the period of March, April and May (MAM) and the second in the period of October, November and December (OND). Central, South and Western areas are characterized by a Unimodal rainfall regime starting from November and continue to the end of April. The Bimodal and Unimodal rainfall patterns are mainly associated with the southward and northward movement of the ITCZ. Most of the rainfall, especially during October to February is convective in nature and distinctly organized.
3. Data and Methodology
Daily rainfall, maximum and minimum temperatures from 18 meteorological stations (Figure 1) for the period 1961-2015 were used in this study. The geographical coordinate and elevation of the stations used in the study are provided in Table 1. These data was collected from the Tanzania Meteorological Agency.
3.2. Data Quality Control
Assessment of the data quality is pre-requisite step in the analysis of climate indices   . It provides opportunity to identify suspicious data and any erroneous outlier that might influence the analysis including daily values outside the user-predefined threshold. It helps to verify whether a reported data value is representative of what was intended to be measured and has not been contaminated by unrelated factors. Software such as RClimdex and recent developed National Climate Monitoring Products (NCMP) software are often used to examine the data and identify and list those that fail pre-defined tests  . This is then followed by expert examinations and judgments of the possible causes of error and determines appropriate corrections. In a case where maximum temperature is indicated to be less than minimum temperature, that data will be flagged out as a suspicious. In this study, rainfall and temperature data were subjected to
Figure 1. Location of the meteorological stations used in the study.
Table 1. Geographical coordinates and station numbers of the stations used in the study.
rigorous quality control check, using the described approach including user- predefined threshold set as the mean of the values for the day plus or minus five times the standard deviation (std) of the value for the day, that is . Daily temperature values outside of these thresholds are manually checked and corrected accordingly. An assessment of unreasonably long consecutive occurrence of the same value, and times when the daily maximum temperature was less than the minimum temperature was done to identify incorrect data entry during the digitization process. The observed errors were corrected using original source of information contained in monthly returns of daily climatological observations report (form no. 626), Synoptic reporting station (form No. 444) and three hourly observation report (form No. 10) stored in Tanzania Meteorological Agency data archive.
3.3. Methodology for Computing the National Climate Monitoring Products for Tanzania (Climate Indices)
Considerable work has been done to date by the World Meteorological Organization (WMO)      in designing methodology and software for computation and analysis of climate indices to facilitate climate monitoring and detection of trends and change in the climate patterns. The Expert Team on Climate Change Detection and Indices (ETCCDI) defined and recommended a sort of 27 core climate extreme indices to enhance study on climate extremes     . Computation of the climate Indices was done using RClimdex and NMCP software, which are R-Based package developed at Climate Research Branch of Meteorological Service of Canada  . Consistent computation and reporting of the indices by countries provide a comprehensive understanding of the regional and global dynamics of the climate extremes and is an important signal for climate change detection. These indices are also used in preparation of various climate-monitoring reports including the reports on status of the global climate, annually produced by the World Meteorological Organization (WMO), and the Bulletin of American Meteorological Society (BAMS). Recently the WMO Commission for Climatology Expert Team on Climate Monitoring Products recommended six climate indices termed as National Climate Monitoring Products (NCMPs) to be consistently computed and reported by all countries to enhance and harmonize climate monitoring  . The computation of the six NCMP involve calculation of the monthly indices for each station or grid point, followed by the interpolation of the station or grid values for each month using ordinary kriging to obtain a spatially-complete analysis on a regular grid. The spatially complete analysis is then averaged across the area covered by the country to calculate the NCMP for that month. The basic steps in computation of the indices involve quality control of the daily data, computation of the indices, computation of the variogram for each Index, interpolation of the data for each Index, and averaging each index across the country. Detailed description of the methods and approach used in the computation of the climate indices using RClimdex and NCMP software is provided in   and can also be accessed from    . The linear trends presented in this study (Table 3 and Figures 2-10) are computed from the indices series using a Kendall’s tau based slope estimator  . This estimator is robust to the effect of outliers in the series. It has been widely used to compute trends in hydro-meteorological series, see for example      . The significance of the trend was determined using the Mann-Kendall test and was estimated at the 5% level  . The trend was considered to be statistically significant if it was less than or equal to a level of 5%     . This test does not assume an underlying probability distribution of the data series. The description of Climate Indices also classified as the NCMP, computed in this study for Tanzania is provided in Table 2.
4. Results and Discussion
The results of the spatial and temporal patterns of the trends in climate extremes indices averaged over the country are presented in Table 3 and in Figures 2-10. The results describe the sign, magnitudes and significance of the trends in mean temperature anomaly, percentage of warm days, percentage of warm nights, percentage of cold days, percentage of cold nights, rainfall anomaly and Standardized Precipitation Index.
Figure 2. Trend in mean annual temperature anomaly for the period 1961-2015.
Figure 3. Trend in the mean percentage of warm days for the period 1961-2015.
4.1. Trends in Temperature Extremes
Trends in mean monthly temperature anomaly averaged over the country for the period 1961-2015 are presented in Table 3 and Figure 2. The bold values represent statistical significant trends at a level of 5% (P < 0.05), while values in italics represent non-statistical significance trends. Mean monthly temperature anomaly averaged over the entire country exhibited a statistical significant increasing temperature trends from January to December ranging from 0.8˚C (April) to 1.14˚C (January).
Percentage of warm days and warm nights for the period 1961-2015 also exhibits a statistical significant increasing trends ranging from 6.21% in March to 15.60% in July for the warm days and from 10.49% (June) to 16.52% (January) for the cold nights (Table 3). These results indicate that nights are warming much faster than days, which literary implies that minimum temperature is becoming warmer than maximum temperature. The percentage of cold days and cold nights are observed to decrease with the decrease being more pronounced in cold nights. The decrease in percentage of cold days ranges from 5.22%
Figure 4. Trend in mean percentage of warm nights for the period 1961-2015.
Table 2. Description of the climate indices, the National Climate Monitoring Products (NCMP).
Table 3. Countrywide trends in temperature and rainfall related extreme indices for the period (1961-2015). In bold are statistical significant trends at 95% confidence interval, while in Italics are non-significant trends.
Figure 5. Trends in spatial distribution of percentage of warm days for the period 1961-2015. Statistical significant trends at 95% level (P < 0.05) are represented by solid triangles. The red symbol for positive trends and blue for negative trends.
Figure 6. Trend in spatial distribution of the percentage of Warm Nights for the period 1961-2015. Statistical significant trends at 95% level (P < 0.05) are represented by solid triangles. The red symbol for positive trends and blue for negative trends.
Figure 7. Trend in Spatial distribution of the percentage of Cold Nights for the period 1961-2015. Statistical significant trends at 95% level (P < 0.05) are represented by solid triangles. The red symbol for positive trends and blue for negative trends.
(March) to −15.21% (June), while for the cold nights the decrease in percentage ranges from −8.14% (April) to −15.36% (January). In general the trend analysis
Figure 8. Spatial distribution of the trends in percentage change of rainfall for the period 1961-2015. Statistical significant trends at 95% confidence level/5% significant level (P < 0.05) are presented by solid triangles and non-statistical significant trend in open triangles. The green colour symbol are for positivetrends and brown colour symbol for negative trends.
Figure 9. Trend in mean percentage rainfall anomaly for the period 1961-2015.
results from temperature extreme indices indicates a coherent and widespread warming in both days and nights.
Figures 2-4 depicts well-pronounced statistical significant increasing trend in mean annual temperature anomaly, mean annual percentage of warm days, and mean annual percentage of warm nights respectively. These results are consistent with the observed temperature trends in different parts of the world as do
Figure 10. Trend in Standardized Precipitation Index for the period 1961-2015.
cumented in different IPCC reports    , and in a series of WMO publication    ). The fifth Assessment Report (AR5) of the IPCC indicates a global average warming of 0.85˚C over the period 1880 to 2012, and that the number of cold days and nights has decreased, while the number of warm days and nights has increased on global scale  .
Spatial distribution of trends in mean percentage of warm days and warm nights are presented in Figure 5 and Figure 6 respectively. A statistical significant increasing trend in percentage of warm days and warm nights is depicted over the eastern part of the country including areas around Kilimanjaro regions, northern coast, southern coast and south-western highlands. Western and eastern parts of the Lake Victoria basin are also characterized by increasing trends of warm days. However, non-statistical decreasing trends in the percentage of warm days and warm nights are depicted in the western parts of the country and in some parts of central areas and the southern part of Lake Victoria basin. Moreover, the western and southernparts of the lake victoria basin showed non- statistical significant decreasing trends of warm nights. These results indicate a clear dipole pattern in temperature dynamics between the eastern side of the country mainly characterized and influenced by the Indian Ocean and the Western side of the country largely influenced by the Congo air mass and westerlies. Figure 7 represents the spatial distribution of the percentage of cold nights for the same period, whereby a decreasing trends is depicted in the entire country, but being more pronounced and statistically significant over the eastern part of the country.
4.2. Trends in Extreme Rainfall Indices
Trends in extreme rainfall indices are presented in Table 3 and Figures 8-10. Non-statistical significant decreasing trends in rainfall patterns are depicted in both percentage changes in rainfall anomaly and in Standardized Precipitation Index (SPI). The spatial distribution of percentage change in rainfall anomaly (Figure 8) indicates a non-statistical significant decreasing rainfall trends over the large parts of the country, except in Songea and Pemba regions, where statistical significant decreasing rainfall trend is depicted. These results underscore the complexity, mixed pattern of change and un-certainties associated with dynamics in rainfall patterns in the region and are consistent with observed global trends.
This study provides an assessment of the observed trends in key climate extreme indices characterized by maximum and minimum temperature, and rainfall. The analyses were done using daily rainfall and daily temperature datasets from 18 meteorological stations for the period 1961-2015. Five temperature extreme indices including temperature anomaly, mean percentages of warm days, warm nights, cold days and cold nights, and two rainfall extreme indices including mean rainfall anomaly and the Standardized Precipitation Index (SPI) were computed and then subjected to trend analysis. It has been observed that both maximum and minimum temperature exhibit statistical significant increasing trends at monthly and annual time scales. Mean annual temperature anomaly has increased by 0.69˚C, while the percentage of warm days has increased by 9.37%. Comparatively mean percentage of warm night has increased by 12.05%, which is relatively faster trend than the trend in warm days. The percentage number of cold days has decreased by 7.64%, while percentage number of cold nights has decreased by 10.0%. Spatial analysis of the observed temperature trend has indicated a coherent and widespread warming over the large part of the country, particularly over the eastern part, while a non-statistical significant decreasing temperature trends is depicted in some few areas over western part of the country. A non-statistical significant decreasing trend is depicted in percentage change in mean rainfall anomaly and in Standardized Precipitation Index. These results are consistent with observed global warming trends and underscore theunderlying complexity and un-certainties in the assessment of precipitation patterns as extensively depicted many literaturesin most parts of the world.
The authors wishes to thank The Tanzania Meteorological Agency for providing data used in this study and WMO for providing guidance in analysis of climate extreme in climate time series.
 Field, C.B., Barros, V., Stocker, T.F., Qin, D., Dokken, D.J., Ebi, K.L., Mastrandrea, M.D., Mach, K.J., Plattner, G.K., Allen, S.K., Tignor, M. and Midgley, P.M. (2012) Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation—Summary for Policymakers, a special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change. World Meteorological Organization, Geneva, Switzerland, 24 p.
 Field, C.B., Barros, V., Stocker, T.F., Qin, D., Dokken, D.J., Ebi, K.L., Mastrandrea, M.D., Mach, K.J., Plattner, G.-K., Allen, S.K., Tignor, M. and Midgley, P.M. (2012) Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, UK, and New York, NY, USA, 582 p.
 Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V. and Midgley, P.M. (2013) Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change—Summary for Policymakers. Cambridge University Press, Cambridge, United Kingdom; New York, New York, USA. 24 p.
 Changa, L.B., Yanda, P.Z. and Ngana, J. (2010) Spatial and Temporal analysis of recent climatological data in Tanzania. Journal of Geography and Regional Planning, 3, 44-65. ISSN 2070-1845©2010 Academic Journals. http://www.academicjournals.org/JGRP
 Luhunga, P., Botai, J. and Kahimba, F. (2016) Evaluation of the Performance of CORDEX Regional Climate Models in Simulating Present Climate Conditions of Tanzania. Journal of Southern Hemisphere Earth Systems Science, 66, 33-55. https://doi.org/10.22499/3.6601.005
 Folland, C.K., Miller, C., Bader, D., Crowe, M., Jones, P., Plumner, N., Richman, M., Parker, D.E., Rogers, J. and Scholefield, P. (1999) Workshop on Indices and indicators for Climate Extremes, Asheville, NC, USA, 3-6 June 1997, Breakout Group C: Temperature Indices for Climate Extremes. Climatic Change, 42, 31-43. https://doi.org/10.1023/A:1005447712757
 Jones, P.D., Horton, E.B., Folland, C.K., Hulme, M., Parker, D.E. and Basnett, T.A. (1999) The Use of Indices to Identify Changes in Climatic Extremes. Climatic Change, 42, 131-149.
 Peterson, T.C., Folland, C., Gruza, G., Hogg, W., Mokssit, A. and Plummer, N. (2001) Report on the Activities of the Working Group on Climate Change Detection and Related Rapporteurs 1998-2001. World Meteorological Organization, Geneva, Rep. WCDMP-47, WMO-TD 1071, 143 p.
 Frich, P., Alexander, L.V., Della-Marta, P., Gleason, B., Haylock, M., Klein Tank, A.M.G. and Peterson, T. (2002) Observed Coherent Changes in Climatic Extremes during the 2nd Half of the 20th Century. Climate Research, 19, 193-212. https://doi.org/10.3354/cr019193
 New, M., Hewitson, B., Stephenson, D.B., Tsiga, A., Kruger, A., Manhique, A., Gomez, B., Coelho, C.A.S., Masisi, D.N., Kululanga, E., Mbambalala, E., Adesina, F., Saleh, H., Kanyanga, J., Adosi, J., Bulane, L., Fortunata, L., Mdoka, M.L. and Lajoie, R. (2006) Evidence of Trends in Daily Climate Extremes over Southern and West Africa. Journal of Geophysical Research, 111, Article ID: D14102. https://doi.org/10.1029/2005JD006289
 Zhang, X., Alexander, L.V., Hegerl, G.C., Klein-Tank, A., Peterson, T.C., Trewin, B. and Zwiers, F.W. (2011) Indices for Monitoring Changes in Extremes Based on Daily Temperature and Precipitation Data. Wiley Interdisciplinary Reviews. Climate Change, 2, 851-870.
 Zhang, X., Aguilar, E., Sensoy, S., Melkonyan, H., Tagiyeva, U., Ahmed, N., Kutaladze, N., Rahimzadeh, F., Taghipour, A., Hantosh, T.H., Albert, P. and Semawi, M. (2005) Trends in Middle East Climate Extreme Indices from 1950 to 2003. Journal of Geophysical Research, 110, D22104. https://doi.org/10.1029/2005JD006181
 Zhang, X., Vincent, L.A., Hogg, W.D. and Niitsoo, A. (2000) Temperature and Precipitation Trends in Canada during the 20th Century. Atmosphere-Ocean, 38, 395-429.
 Luhunga, P., Mutayoba, E. and Ng’ongolo, H. (2014) Homogeneity of Monthly Mean Air Temperature of the United Republic of Tanzania with HOMER. Atmospheric and Climate Sciences, 4, 70-77. https://doi.org/10.4236/acs.2014.41010