Surface Area Variability of a North-Central Tanzanian Crater Lake

: A history of modern (1973–2015) surface area variability for Lake Basotu in north-central Tanzania has been reconstructed using archived Landsat images from the dry season between June and October. This record was compared to local weather data as well as larger scale weather patterns. The lake has been in a state of decline interrupted by major ﬂood events since the beginning of the satellite record. From 1973 to 1997, the lake area was between 0.97 km 2 and 4.28 km 2 . Lake extent abruptly increased to 13.86 km 2 in 1998, when a co-occurrence of El Niño and a positive Indian Ocean Dipole led to extensive ﬂooding. It is hypothesized that local agricultural practices leading to soil erosion and subsequent basin sedimentation have most likely increased the sensitivity of Lake Basotu to climatic ﬂuctuations.


Introduction
With food and water security issues at the forefront of concerns over global climate change, it is imperative to understand the sensitivity of freshwater resources, especially in semi-arid regions.Reconstructing the temporal variability of lakes may be the key to distinguishing which changes are human-made and which are due to natural variability.
The authors present a record of surface area variation for a Tanzanian crater lake created via a simple and effective method using archived Landsat images.Lake Basotu (4 ˝22 1 11 11 S, 35 ˝41 22 11 E) in the Hanang district is an important freshwater resource that has been at the center of social conflict in recent years [1][2][3].Intensive agriculture practices have been cited as responsible for increased siltation, pollution, and reduction in water level.As this area mainly comprises soils with a high clay content (locally known as "mbuga"), it is highly susceptible to extreme weather events and is, therefore, a suitable location for investigating climatic effects in a semi-arid environment.Moreover, logistically, flood prediction is important to infrastructure, the economy, and preventing loss of life.
The Tanzanian government increased budgeted funds for climate change-related activities by 57% between 2009-2010 and 2012-2013, with 70% allocated for mitigation and adaptive strategies [4].In March 2015, dozens of lives were lost and hundreds of people were left homeless due to flooding in Dar es Salaam and the Shinyanga region of Tanzania [5,6].
It is very likely for the frequency of heavy precipitation events to increase over tropical regions in the 21st century [7].During times of drought, soil becomes compacted.This, combined with periodic and extreme precipitation events, increases the chance of flooding and less effective recharge of groundwater, which will negatively impact agricultural viability.Mbuga soils are rich in clay and characterized by developing deep cracks upon drying, which is typical of vertisols.Upon rewetting, the cracks fill and close due to swelling of clays, making the soil nearly impervious to further infiltration [8].This leads to runoff, and vertisols are consequently susceptible to gully formation.Even in low rainfall areas, many Tanzanian storms surpass rainfall intensities required for surficial erosion [9].
Flood events in Tanzania have been linked to the co-occurrence of an abnormally strong El Niño coinciding with a positive phase of the Indian Ocean Dipole (IOD) [10].The effects can be generally positive in Tanzania, with an increase in rain leading to gains in crop yields, although experiences from the 1997-1998 El Niño actually resulted in localized losses and damages to infrastructure.Losses associated with heavy rainfall were greatest in low-lying areas with predominantly heavy clay soils, such as mbuga [11].
An El Niño event is evident when the five-month running means of sea surface temperature (SST) anomalies exceed 0.4 ˝C for six months or more in the region between 5 ˝N-5 ˝S and 120 ˝-170 ˝W [12].One result of this disturbance on global climate is the equator-ward shift of the ITCZ.Due to this, equatorial eastern Africa experiences wetter than normal conditions during its rainy season, with the strongest relationship evident from October to December [13,14].
The IOD signifies the difference in SST between Sumatra and the western Indian Ocean.During a positive phase, warmer waters are present along the eastern coast of Africa, with a reversal of the easterly trade winds, resulting in increased precipitation over tropical East Africa [15].The IOD is commonly represented as the Dipole Mode Index (DMI), which is the difference between these anomalies.Using the Hadley Center Sea Ice and Sea Surface Temperature dataset, with the linear trend of SST removed, monthly DMI was calculated from 1960 to 2014 [16].
A greater understanding of the mechanisms behind flood events is integral to planning efforts for the agricultural sector as well as disaster preparation and mediation.This study aims to provide a record of surface area variation for Lake Basotu in Tanzania.Lake extent data are compared to both indices of large-scale weather patterns and local meteorological data in order to ascertain what is associated with natural variability and what may be an effect of human influence on the environment.

Study Location
The Basotu plains of the Hanang district contain 13 shallow maar lakes that act as the primary freshwater source for local people, as there are no perennial rivers.The area directly surrounding Lake Basotu (4 ˝22 1 11 11 S, 35 ˝41 22 11 E) is comprised of mainly acacia savanna, with soils derived from volcanic rocks and calcareous tuffs (Figure 1).Precipitation in this area is driven by the movement of the ITCZ, and results in two rainy seasons.The primary rains come from March to June, centered on April-May, with a secondary season from September to December, centered on October-November.Lake Basotu and its surrounding land have been at the center of conflict since the 1960s, when local pastoralists, the Barabaig, were displaced in order to implement mechanized wheat farming [17,18].Since then, the lake has suffered from sediment input, drainage from irrigation, and eutrophication due to increasing nutrient loads [19].As the lake continues to decrease in size, the concentration of pollutants will likely increase, which will also have implications for local fishing and the viability of this lake as a source of drinking water.

Methods
A geologic map compiled by the Survey Division of the Ministry of Lands, Settlement and Water Development of the Government of Tanzania (1963) was obtained from the Geological Survey of Sweden [20].All available georeferenced Landsat images of the Basotu area captured during the long dry season as well as the one-second digital elevation model (DEM) by the Shuttle Radar Topography Mission (SRTM) were obtained from the United States Geological Survey (USGS) EarthExplorer via their bulk download service [21].Landsat products are terrain-corrected using the Level 1 Product Generation System [22].All scenes in this study were processed by Standard Terrain Correction (Level 1T), which is the highest level possible, and no additional co-registration was performed.
From the available Landsat images, 20 years of dry-season captures were suitable for analysis.A gap in coverage from 1974 to 1978 can be explained by a more conservative approach to image acquisition with the launch of Landsat-2, which was restricted to on-board data storage and primarily responsible for the years 1975 to 1977.A second gap in coverage occurred from 1980 to 1983.Landsat-3 (1978Landsat-3 ( -1983) ) was prone to technical problems, and it was not until 1984 that acquisition rates were increased with the launch of Landsat-5 [23].
Due to the strong contrast between land and water reflectance in relatively shallow lakes, use of the Landsat satellite series in tracking extent changes is appropriate [24].Dry season images were chosen to avoid atmospheric interference by cloud cover and to maximize the reflective differences between land and water.Maximum lake extent (5 April 2007) and seasonal in/outflows were determined from cloud-free rainy season images.Lake shorelines were digitized manually using false-color composites alternating between the near-infrared (NIR), shortwave-infrared (SWIR), and blue-light wavelengths of available Landsat images (Figure 2).The NIR bands were chosen to emphasize shorelines, while SWIR offered the potential to discriminate the moisture content of soil.In Landsat 1-3, where the SWIR and blue bands are not available, composites were alternated between the two NIR (land-water boundaries) and green (shallow water) bands.
The surface area of Lake Basotu was calculated by manually digitizing the lake extent from 20 years using ArcGIS 10.3.1 ® (Environmental Systems Research Institute, Inc., Redlands, CA, USA) (Table 1) [25].Due to the unavailability of bathymetric information, calculations of lake volume were not possible.Manual digitization of satellite images is both subjective and time-consuming, and so applied to small study areas.The error from this method was found to be between ±1.5 m for distinct shorelines, and ±3.0 m for diffuse shorelines, when no binary encoding was used to facilitate landwater boundaries [26].Values for uncertainty in manual delineation of shorelines were determined by an experiment of multiple operators with different levels of expertise.Levels of expertise were based on field and remote sensing experience related to the study area, and users were divided into

Methods
A geologic map compiled by the Survey Division of the Ministry of Lands, Settlement and Water Development of the Government of Tanzania (1963) was obtained from the Geological Survey of Sweden [20].All available georeferenced Landsat images of the Basotu area captured during the long dry season as well as the one-second digital elevation model (DEM) by the Shuttle Radar Topography Mission (SRTM) were obtained from the United States Geological Survey (USGS) EarthExplorer via their bulk download service [21].Landsat products are terrain-corrected using the Level 1 Product Generation System [22].All scenes in this study were processed by Standard Terrain Correction (Level 1T), which is the highest level possible, and no additional co-registration was performed.
From the available Landsat images, 20 years of dry-season captures were suitable for analysis.A gap in coverage from 1974 to 1978 can be explained by a more conservative approach to image acquisition with the launch of Landsat-2, which was restricted to on-board data storage and primarily responsible for the years 1975 to 1977.A second gap in coverage occurred from 1980 to 1983.Landsat-3 (1978Landsat-3 ( -1983) ) was prone to technical problems, and it was not until 1984 that acquisition rates were increased with the launch of Landsat-5 [23].
Due to the strong contrast between land and water reflectance in relatively shallow lakes, use of the Landsat satellite series in tracking extent changes is appropriate [24].Dry season images were chosen to avoid atmospheric interference by cloud cover and to maximize the reflective differences between land and water.Maximum lake extent (5 April 2007) and seasonal in/outflows were determined from cloud-free rainy season images.Lake shorelines were digitized manually using false-color composites alternating between the near-infrared (NIR), shortwave-infrared (SWIR), and blue-light wavelengths of available Landsat images (Figure 2).The NIR bands were chosen to emphasize shorelines, while SWIR offered the potential to discriminate the moisture content of soil.In Landsat 1-3, where the SWIR and blue bands are not available, composites were alternated between the two NIR (land-water boundaries) and green (shallow water) bands.
The surface area of Lake Basotu was calculated by manually digitizing the lake extent from 20 years using ArcGIS 10.3.1 ® (Environmental Systems Research Institute, Inc., Redlands, CA, USA) (Table 1) [25].Due to the unavailability of bathymetric information, calculations of lake volume were not possible.Manual digitization of satellite images is both subjective and time-consuming, and so applied to small study areas.The error from this method was found to be between ˘1.5 m for distinct shorelines, and ˘3.0 m for diffuse shorelines, when no binary encoding was used to facilitate land-water boundaries [26].Values for uncertainty in manual delineation of shorelines were determined by an experiment of multiple operators with different levels of expertise.Levels of expertise were based on field and remote sensing experience related to the study area, and users were divided into groups as either having one or both.In this experiment, manual delineation was found to be more accurate than semi-automatic classification when data quality was poor or uncertain.groups as either having one or both.In this experiment, manual delineation was found to be more accurate than semi-automatic classification when data quality was poor or uncertain.Another source of uncertainty may arise from intra-annual variability, as this study does not account for seasonality of the lake.In 2007, the year of greatest lake extent, the wet season lake extent was calculated to be 16.11 km 2 , which is roughly 14% more than the dry season value of 14.12 km 2 .The selection of dry season images offers continuity in analysis, and should provide a reliable assessment of inter-annual variability that is adequate for the scope of this investigation.
As separating shallow water from fen-like vegetation is the most challenging aspect of the manual classification, using false color composites is the most appropriate way possible to delineate these borders.Other sources of errors may stem from differences in resolution between Landsat 1-3 (79 m) and 4-8 (30 m) as well as the 2003 failure of the Scanning Line Corrector in Landsat 7, which resulted in coverage gaps [27].
The Lake Basotu watershed was delineated in ArcGIS 10.3.1 ® using the hydrology toolbox and the SRTM 30 m DEM [28].The historical outflow point for the maximum lake extent was determined based on the flow paths derived from the DEM (Figure 3).Another source of uncertainty may arise from intra-annual variability, as this study does not account for seasonality of the lake.In 2007, the year of greatest lake extent, the wet season lake extent was calculated to be 16.11 km 2 , which is roughly 14% more than the dry season value of 14.12 km 2 .The selection of dry season images offers continuity in analysis, and should provide a reliable assessment of inter-annual variability that is adequate for the scope of this investigation.
As separating shallow water from fen-like vegetation is the most challenging aspect of the manual classification, using false color composites is the most appropriate way possible to delineate these borders.Other sources of errors may stem from differences in resolution between Landsat 1-3 (79 m) and 4-8 (30 m) as well as the 2003 failure of the Scanning Line Corrector in Landsat 7, which resulted in coverage gaps [27].
The Lake Basotu watershed was delineated in ArcGIS 10.3.1 ® using the hydrology toolbox and the SRTM 30 m DEM [28].The historical outflow point for the maximum lake extent was determined based on the flow paths derived from the DEM (Figure 3).Discontinuous time series of precipitation, acquired from the Tanzanian Meteorological Agency, were available for two nearby stations, namely Babati and Singida agricultural stations (Figure 4).The two stations had coverage of 81% and 91%, respectively, for the period of 1970-2010.A distanceweighted average of both stations was used to estimate precipitation for this period.An additional station, Singida meteorological station, was also available for the period 2000-2015, with coverage of 93% for the period.Precipitation for the period 2010-2015 was estimated based on a linear regression between the distance-weighted average of Babati and Singida agricultural stations and the Singida meteorological station for 2000-2010 (r 2 = 0.80).Global Precipitation and Climatology Center v6 (GPCC) (The National Oceanic and Atmospheric Administration -Earth System Research Laboratory, Boulder, Colorado, USA) data were used to fill the missing data.GPCC was used since it had previously been shown to provide good precipitation estimates in Tanzania [29,30].The GPCC data were downscaled using quantile mapping considering the scale of the catchment [31].The downscaled data were then further Discontinuous time series of precipitation, acquired from the Tanzanian Meteorological Agency, were available for two nearby stations, namely Babati and Singida agricultural stations (Figure 4).The two stations had coverage of 81% and 91%, respectively, for the period of 1970-2010.A distance-weighted average of both stations was used to estimate precipitation for this period.An additional station, Singida meteorological station, was also available for the period 2000-2015, with coverage of 93% for the period.Precipitation for the period 2010-2015 was estimated based on a linear regression between the distance-weighted average of Babati and Singida agricultural stations and the Singida meteorological station for 2000-2010 (r 2 = 0.80).Discontinuous time series of precipitation, acquired from the Tanzanian Meteorological Agency, were available for two nearby stations, namely Babati and Singida agricultural stations (Figure 4).The two stations had coverage of 81% and 91%, respectively, for the period of 1970-2010.A distanceweighted average of both stations was used to estimate precipitation for this period.An additional station, Singida meteorological station, was also available for the period 2000-2015, with coverage of 93% for the period.Precipitation for the period 2010-2015 was estimated based on a linear regression between the distance-weighted average of Babati and Singida agricultural stations and the Singida meteorological station for 2000-2010 (r 2 = 0.80).Global Precipitation and Climatology Center v6 (GPCC) (The National Oceanic and Atmospheric Administration -Earth System Research Laboratory, Boulder, Colorado, USA) data were used to fill the missing data.GPCC was used since it had previously been shown to provide good precipitation estimates in Tanzania [29,30].The GPCC data were downscaled using quantile mapping considering the scale of the catchment [31].The downscaled data were then further Global Precipitation and Climatology Center v6 (GPCC) (The National Oceanic and Atmospheric Administration -Earth System Research Laboratory, Boulder, Colorado, USA) data were used to fill the missing data.GPCC was used since it had previously been shown to provide good precipitation estimates in Tanzania [29,30].The GPCC data were downscaled using quantile mapping considering the scale of the catchment [31].The downscaled data were then further corrected based on a linear regression between the distance-weighted gauge average and the downscaled GPCC data for the period 1970-2010 (r 2 = 0.75).

Results with Discussion
In maar lakes, sedimentation by diatreme and subsequent lacustrine material seal the explosion crater from its feeder vent [32,33].The closed basin state of Lake Basotu was confirmed by testing of geochemical mechanisms and solute evolution pathways.For this reason, groundwater is not believed to be a current source of variability in lake size [32].There have been three distinct periods of lake extent decline observed since 1973 (A-C in Figure 5).Between these are the major flood events of 1997-1998 and 2006-2007, which were observed in various locations across eastern Africa [10,34,35].corrected based on a linear regression between the distance-weighted gauge average and the downscaled GPCC data for the period 1970-2010 (r 2 = 0.75).

Results with Discussion
In maar lakes, sedimentation by diatreme and subsequent lacustrine material seal the explosion crater from its feeder vent [32,33].The closed basin state of Lake Basotu was confirmed by testing of geochemical mechanisms and solute evolution pathways.For this reason, groundwater is not believed to be a current source of variability in lake size [32].There have been three distinct periods of lake extent decline observed since 1973 (A-C in Figure 5).Between these are the major flood events of 1997-1998 and 2006-2007, which were observed in various locations across eastern Africa [10,34,35].Period A, from 1973 to 1997 shows a general decreasing trend of ~0.11 km 2 /year.After the dramatic increase in lake size from 1997 (0.97 km 2 ) to 1998 (13.86 km 2 ), there was a more rapid decline of ~0.92 km 2 /year until the next flood event from 2006 to 2007.Since 2008, the general trend has been a reduction of ~0.32 km 2 /year; however, a slight increase has been noted since 2013.The surface area level towards the end of Part C may have been affected by local conflicts over irrigation leading to a stabilization of lake levels.In 2011, an organized protest of local women demanded government officials to end the large-scale drawing of lake water, which they depend on for their daily water supplies [36].
The DMI and Multivariate El Niño-Southern Oscillation Index (MEI) from October to December (OND) were averaged and standardized for the period 1960-2014 in order to search for similar patterns when flood events would be expected to occur (Figure 6).As both indices for 1997 and 2006 are above z-scores of +1, this was chosen as the base level for flood events.Based on this, flooding would have been expected in the 1972-1973, 1982-1983, and 1994-1995  To improve flood prediction in East Africa, it would be helpful to understand the driving mechanisms behind enhanced precipitation events.The co-occurrence of El Niño and a positive IOD alone is not enough to accurately predict flood-level precipitation evident in the 1982-1983 and 1994-1995 seasons, when enhanced precipitation was not observed in the reanalysis data.Period A, from 1973 to 1997 shows a general decreasing trend of ~0.11 km 2 /year.After the dramatic increase in lake size from 1997 (0.97 km 2 ) to 1998 (13.86 km 2 ), there was a more rapid decline of ~0.92 km 2 /year until the next flood event from 2006 to 2007.Since 2008, the general trend has been a reduction of ~0.32 km 2 /year; however, a slight increase has been noted since 2013.The surface area level towards the end of Part C may have been affected by local conflicts over irrigation leading to a stabilization of lake levels.In 2011, an organized protest of local women demanded government officials to end the large-scale drawing of lake water, which they depend on for their daily water supplies [36].
The DMI and Multivariate El Niño-Southern Oscillation Index (MEI) from October to December (OND) were averaged and standardized for the period 1960-2014 in order to search for similar patterns when flood events would be expected to occur (Figure 6).As both indices for 1997 and 2006 are above z-scores of +1, this was chosen as the base level for flood events.Based on this, flooding would have been expected in the 1972-1973, 1982-1983, and 1994-1995   The mechanized wheat farming program, developed in response to declining food security, was a cooperation between the Tanzanian Government and the Canadian International Development Agency, and was primarily based in the district surrounding Lake Basotu [36].Between 1973 and 1997, there was a negative correlation (r(23) = −0.902,p = 0.000) between the calculated surface area of Lake Basotu and the number of hectares (ha) seeded in the Hanang Wheat Complex (Figure 7).The values for the seeded area were based on statistics reported by Folster [37].The measurements for 1998 were excluded due to the large flood event, which is considered to be an outlier.Missing values were calculated using linear interpolation before a Pearson correlation was performed.The negative correlation observed contradicts the wheat project being reported as exclusively rainfed [38].Between 1973 and 1997, there is no significant relationship between lake surface area and precipitation.Correlations were calculated between surface area and precipitation totals for both the short wet season of the previous year (November-December) (r(25) = 0.129, p = 0.540) and the sum of the months preceding the dry season of the same year (January-May) (r(25) = -0.316,p = 0.124).Additional information is needed to determine if and how this agricultural program affected the decline in lake extent.Increased soil erosion and sediment transport can be assumed to have occurred during the time period of intensive wheat farming in the area around Lake Basotu, causing further infilling of the lake basin.Financial support for the Canadian-Tanzanian wheat program officially ended in 1993, but wheat farming continued until 1998.Flood events since 1996 may be due to filling of vertisol cracks by agricultural runoff, bringing erosional material into the lake basin and preventing infiltration of water into the already dense soil.This process should be visible in Landsat satellite images by the formation of gilgai, which are alternating micro-depressions and micro-mounds common to vertisols [39].With an altitudinal difference of <1 m, these features To improve flood prediction in East Africa, it would be helpful to understand the driving mechanisms behind enhanced precipitation events.The co-occurrence of El Niño and a positive IOD alone is not enough to accurately predict flood-level precipitation evident in the 1982-1983 and 1994-1995 seasons, when enhanced precipitation was not observed in the reanalysis data.
The mechanized wheat farming program, developed in response to declining food security, was a cooperation between the Tanzanian Government and the Canadian International Development Agency, and was primarily based in the district surrounding Lake Basotu [36].Between 1973 and 1997, there was a negative correlation (r(23) = ´0.902,p = 0.000) between the calculated surface area of Lake Basotu and the number of hectares (ha) seeded in the Hanang Wheat Complex (Figure 7).The values for the seeded area were based on statistics reported by Folster [37].The measurements for 1998 were excluded due to the large flood event, which is considered to be an outlier.Missing values were calculated using linear interpolation before a Pearson correlation was performed.The negative correlation observed contradicts the wheat project being reported as exclusively rainfed [38].Between 1973 and 1997, there is no significant relationship between lake surface area and precipitation.Correlations were calculated between surface area and precipitation totals for both the short wet season of the previous year (November-December) (r(25) = 0.129, p = 0.540) and the sum of the months preceding the dry season of the same year (January-May) (r(25) = -0.316,p = 0.124).Additional information is needed to determine if and how this agricultural program affected the decline in lake extent.The mechanized wheat farming program, developed in response to declining food security, was a cooperation between the Tanzanian Government and the Canadian International Development Agency, and was primarily based in the district surrounding Lake Basotu [36].Between 1973 and 1997, there was a negative correlation (r(23) = −0.902,p = 0.000) between the calculated surface area of Lake Basotu and the number of hectares (ha) seeded in the Hanang Wheat Complex (Figure 7).The values for the seeded area were based on statistics reported by Folster [37].The measurements for 1998 were excluded due to the large flood event, which is considered to be an outlier.Missing values were calculated using linear interpolation before a Pearson correlation was performed.The negative correlation observed contradicts the wheat project being reported as exclusively rainfed [38].Between 1973 and 1997, there is no significant relationship between lake surface area and precipitation.Correlations were calculated between surface area and precipitation totals for both the short wet season of the previous year (November-December) (r(25) = 0.129, p = 0.540) and the sum of the months preceding the dry season of the same year (January-May) (r(25) = -0.316,p = 0.124).Additional information is needed to determine if and how this agricultural program affected the decline in lake extent.Increased soil erosion and sediment transport can be assumed to have occurred during the time period of intensive wheat farming in the area around Lake Basotu, causing further infilling of the lake basin.Financial support for the Canadian-Tanzanian wheat program officially ended in 1993, but wheat farming continued until 1998.Flood events since 1996 may be due to filling of vertisol cracks by agricultural runoff, bringing erosional material into the lake basin and preventing infiltration of water into the already dense soil.This process should be visible in Increased soil erosion and sediment transport can be assumed to have occurred during the time period of intensive wheat farming in the area around Lake Basotu, causing further infilling of the lake basin.Financial support for the Canadian-Tanzanian wheat program officially ended in 1993, but wheat farming continued until 1998.Flood events since 1996 may be due to filling of vertisol cracks by agricultural runoff, bringing erosional material into the lake basin and preventing infiltration of water into the already dense soil.This process should be visible in Landsat satellite images by the formation of gilgai, which are alternating micro-depressions and micro-mounds common to vertisols [39].With an altitudinal difference of <1 m, these features would not be captured by the available DEM.If the distance between the gilgai features were at a range that would be visually distinguishable in Landsat images, the presence of swamp conditions in the basin prevents this observation.
Using Landsat images, the area of cultivation within the Basotu watershed was determined and digitized manually.The percentage of this watershed dedicated to agriculture increased from approximately 23% in 1973 to 43% in 1997 (Figure 8).would not be captured by the available DEM.If the distance between the gilgai features were at a range that would be visually distinguishable in Landsat images, the presence of swamp conditions in the basin prevents this observation.Using Landsat images, the area of cultivation within the Basotu watershed was determined and digitized manually.The percentage of this watershed dedicated to agriculture increased from approximately 23% in 1973 to 43% in 1997 (Figure 8).

Conclusions
According to the evidence presented, it is concluded that fluctuations in Lake Basotu appear to be tied to both climatological factors and anthropogenic influences resulting from agricultural practices.In order to use the flooding of Lake Basotu as an indicator of weather patterns with the potential for large flood events, it would be beneficial to understand the reasons behind non-flooding years in the past.Flood events occurring after 1996 may be an effect of increased sedimentation within the lake basin from the mechanized wheat farming program implemented in the Hanang district.This is feasible, as the local news has stated the lake is "overwhelmed with sediments" [40].
There is high confidence that the strongest ocean warming due to climate change will be experienced in the tropics, and it is "very likely" that this will increase the intensity of El Niño events [7] (p.60).As Lake Basotu is an important freshwater resource, further investigation into past lake fluctuations is required to determine if these modern fluctuations are beyond the scope of normal environmental variability.
Lake surface area, Figure S7: Area of wheat seeded (thousands of ha) in the Hanang district versus Lake Basotu surface area calculations (km 2 ), Figure S8: Agricultural area within the Basotu watershed.

Figure 1 .
Figure 1.Map of Lake Basotu surroundings (adapted from a 1963 topographic map published by the Survey Division of the Ministry of Lands, Settlement and Water Development of the Government of Tanzania [20]).

Figure 1 .
Figure 1.Map of Lake Basotu surroundings (adapted from a 1963 topographic map published by the Survey Division of the Ministry of Lands, Settlement and Water Development of the Government of Tanzania [20]).

Figure 2 .
Figure 2. False -color composite Landsat image (July 2015) of Lake Basotu with manual shore de line ation shown in white .

Figure 2 .
Figure 2. False-color composite Landsat image (July 2015) of Lake Basotu with manual shore delineation shown in white.

Figure 5 .
Figure 5. Calculated lake surface area record for Lake Basotu.
seasons.Although there are no satellite images from 1980 to 1983, it is unlikely that a large flood event would have occurred due to the similarly low values in 1979 and 1984.Of all monthly precipitation values, the indices for OND, DMI, and MEI were found to be most significantly correlated to December precipitation during the time period of 1970 to 2014 (MEI: r(43) = 0.425, p = 0.004: DMI: r(43) = 0.561, p = 0.000).To facilitate comparison, precipitation values were standardized to the period of 1970-2015.For both 1997 and 2006, above average precipitation is observed in conjunction with the suggested base level for flooding (+1 STD, DMI, and MEI).

Figure 5 .
Figure 5. Calculated lake surface area record for Lake Basotu.
seasons.Although there are no satellite images from 1980 to 1983, it is unlikely that a large flood event would have occurred due to the similarly low values in 1979 and 1984.Of all monthly precipitation values, the indices for OND, DMI, and MEI were found to be most significantly correlated to December precipitation during the time period of 1970 to 2014 (MEI: r(43) = 0.425, p = 0.004: DMI: r(43) = 0.561, p = 0.000).To facilitate comparison, precipitation values were standardized to the period of 1970-2015.For both 1997 and 2006, above average precipitation is observed in conjunction with the suggested base level for flooding (+1 STD, DMI, and MEI).

Figure 6 .
Figure 6.(Lower graph): Standardized indices for El Niño and the Indian Ocean Dipole and regional precipitation records.(Upper graph): Lake surface area.

Figure 7 .
Figure 7. Area of wheat seeded (thousands of ha) in the Hanang district versus Lake Basotu surface area calculations (km 2 ).

6 .
(Lower graph): Standardized indices for El Niño and the Indian Ocean Dipole and regional precipitation records.(Upper graph): Lake surface area.

Figure 6 .
Figure 6.(Lower graph): Standardized indices for El Niño and the Indian Ocean Dipole and regional precipitation records.(Upper graph): Lake surface area.

Figure 7 .
Figure 7. Area of wheat seeded (thousands of ha) in the Hanang district versus Lake Basotu surface area calculations (km 2 ).

Figure 7 .
Figure 7. Area of wheat seeded (thousands of ha) in the Hanang district versus Lake Basotu surface area calculations (km 2 ).