Land cover, land use, and climate change impacts on endemic cichlid habitats in Northern Tanzania

: Freshwater ecosystems are among the most threatened on Earth, facing environmental and anthropogenic pressures often surpassing their terrestrial counterparts. Land use and land cover change (LUCC) such as degradation and fragmentation of the terrestrial landscape negatively impacts aquatic ecosystems. Satellite imagery allows for an impartial assessment of the past to determine habitat alterations. It can also be used as a forecasting tool in the development of species conservation strategies through models based on ecological factors extracted from imagery. In this study, we analyze Landsat time sequences (1984–2015) to quantify LUCC around three freshwater ecosystems with endemic cichlids in Tanzania. In addition, we examine population growth, agricultural expansion, and climate change as stressors that impact the habitats. We found that the natural vegetation cover surrounding Lake Chala decreased from 15.5% (1984) to 3.5% (2015). At Chemka Springs, we observed a decrease from 7.4% to 3.5% over the same period. While Lake Natron had minimal LUCC, severe climate change impacts have been forecasted for the region. Subsurface water data from the Gravity Recovery and Climate Experiment (GRACE) satellite observations further show a decrease in water resources for the study areas, which could be exacerbated by increased need from a growing population and an increase in agricultural land use.


Introduction
Freshwater aquatic ecosystems are under increasing threat from human activity and climate change.They are among the most heavily altered ecosystems with disproportionately high biodiversity loss worldwide [1,2].Major systematic drivers of aquatic species loss include land cover and land use change (LUCC), overexploitation, invasive species, and climate change [3][4][5][6][7][8].Direct and indirect competition for water resources with humans (e.g., water abstraction for irrigation) impose stress upon aquatic ecosystems and imperil fauna [9].Aquatic ecosystems are particularly vulnerable because they receive the cumulative effects of stressors within the watershed.Reconstruction of a region's LUCC history can be used to determine the potential stressors that have impacted or degraded habitats over time and thereby can enhance current and future management [10].While some sources of aquatic habitat degradation such as the construction of manmade dams and river impoundments are obvious, others are more obscure.For example, land burning as the primary cause of eutrophication in Lake Victoria was only recently determined [11].
Combining standardized geospatial data of population and climate (e.g., [12][13][14]) with remotely sensed data can produce spatially explicit models of a region's environmental geography, distribution of endangered species and biodiversity at regional and global scales [15].Satellite imagery provides critical information for assessing changes in terrestrial ecosystems that can be linked to effects on aquatic ecosystems [16].In Tanzania, an assessment of the Usangu Plains wetlands and Malagarasi River drainage concluded that these large river catchments suffered substantial land cover changes with declines in woodland and wetland habitats co-occurring with increases in settlements, agriculture and grassland [17,18].In the Usangu Plains wetlands in particular, vegetated swamp cover decreased by 67% between 1984 and 2001 [18].
Relatively little is known about the endemic African fishes outside of the Great Lakes (primarily Lakes Victoria, Tanganyika and Malawi).Habitat specializations, phylogeny, reproduction and diet have been sparsely documented for fishes inhabiting the smaller lakes, wetlands and river systems in Tanzania, while LUCC and climate change effects on the species are almost entirely unknown [19].Many endemic fishes are highly specialized to their local environment (e.g., rheophilic cichlids and thermal-alkaline-salt tolerant cichlids).The survival of others, such as annual killifish that rely on ephemeral pools in the savannah, is critically dependent upon the land use and land cover [20].Even recently, new species of fishes have been discovered in areas affected by LUCC, for example, three new species of sucker-mouth catfishes (genus Chiloglanis) were found in the lower Malagrasi River [21], and a previously unknown rheophilic cichlid, Haplochromis vanheusdeni (Schedel, Friel & Schliewen 2014), was discovered in the Great Ruaha River drainage [22].Others such as the characin Petersius conserialis (Hilgendorf 1894), from the Ruvu and Rufiji Rivers are only known from past sightings.LUCC also affects endemic species in lakes, such as the cichlids inhabiting African crater lakes (e.g., Lakes Ejagham, Bermin, Barombi-Mbo, and Bosumtwi) that evolved in small bodies of water and are sensitive to environmental change.Some of these cichlids are also highly specialized in their diets, for example Pungu maclareni (Trewavas 1962), and Coptodon spongotroktis (Stiassny, Schliewen & Dominey 1992), feed predominantly on freshwater sponges [23].The primary threats to the aquatic fauna of these crater lakes are siltation due to LUCC and water extraction for agriculture and other uses [24].
In the context of this study, we adopt the definition of vulnerability provided by the fourth assessment report from the Intergovernmental Panel on Climate Change (IPCC) as "the degree to which a system is susceptible to, and unable to cope with, adverse effects of climate change, including climate variability and extremes" [25].Climate change projected for East Africa, including decreased precipitation and increased temperature affect freshwater fish vulnerability by lessening the capacity of water bodies to regulate against temperature change, altering the thermal suitability of the habitats [26] and, in extreme cases, the spatial extent of the habitats and resulting in a decrease in the natural ecosystem's ability to mitigate against the effects of climate change.Native riparian vegetation is known to lessen the impacts of LUCC and climate change on aquatic ecosystems by improving water quality [27][28][29][30], moderating water temperature [29,31], providing food and resources [32] and improving biodiversity [29,31,33].Therefore, degradation or conversion of the riparian vegetation to other land cover types such as crops unavoidably results in the deterioration of aquatic ecosystems [34].In addition to the integrity of the vegetation adjacent to aquatic habitats, historical LUCC (i.e., land use legacy) at larger landscape scales has been shown to be important for fish and invertebrate diversity [10,35].Aquatic habitat quality is dependent upon the extent of LUCC at catchment or regional scales [10].Furthermore, the rarity of many individual species within diverse freshwater communities often hinders conservation efforts.A lack of thorough natural history data on rare species results in extreme difficulties in forecasting the identity of taxa affected by climate change or LUCC at large scales [24].There is a fundamental gap in knowledge of the vulnerability of endemic species to LUCC, climate change and socio-economic factors that threaten their habitats.
The overall objective of this study was to examine the LUCC from satellite imagery over the 1984-2015 period for three freshwater habitats of endemic cichlids.We also summarize putative climate and landscape stressors including population growth and forecasted climate change.
The family Cichlidae contains 1693 described species (as of May 2017), with many new species yet to be described.Cichlids are distributed exclusively in tropical fresh water, have colonized most major bodies of water, and share these habitats with a great diversity of other aquatic species.Because many have limited ranges and high levels of specialization and/or endemism, they are good indicators for overall aquatic habitat degradation from LUCC and climate change.

Study Areas
The study was carried out in three inland freshwater habitats: (1) the thermal freshwater springs on the southern end of Lake Natron, (2) Chemka Springs, a freshwater spring near the town of Bomang'ombe, and (3) Lake Chala a crater lake in the foothills of Mount Kilimanjaro (Figure 1, Supplementary Video S1).All three sites are arid to semi-arid with a bimodal pattern of precipitation.The long rainy season spans from March to May with a short rainy season from October to December.The family Cichlidae contains 1693 described species (as of May 2017), with many new species yet to be described.Cichlids are distributed exclusively in tropical fresh water, have colonized most major bodies of water, and share these habitats with a great diversity of other aquatic species.Because many have limited ranges and high levels of specialization and/or endemism, they are good indicators for overall aquatic habitat degradation from LUCC and climate change.

Study Areas
The study was carried out in three inland freshwater habitats: (1) the thermal freshwater springs on the southern end of Lake Natron, (2) Chemka Springs, a freshwater spring near the town of Bomang'ombe, and (3) Lake Chala a crater lake in the foothills of Mount Kilimanjaro (Figure 1, Supplementary Video S1).All three sites are arid to semi-arid with a bimodal pattern of precipitation.The long rainy season spans from March to May with a short rainy season from October to December.Lake Natron, a thermal hypersaline lake, is part of a Ramsar Wetland site but is unprotected except for hunting regulations [36].The main water source for this shallow (<3 m) lake is the Southern Ewaso Ng'iro River originating in Kenya.The lake belongs to the Southern Eastern Rift ecoregion classified as having 'globally outstanding' aquatic biological distinctiveness with a 'critical' conservation status [19].The surrounding anthropogenic biomes are predominantly rangeland, with small areas considered 'semi-natural' and 'cropland' [37].Located on the Southern flank of the Kenyan dome in the Gregory Rift at relatively low elevation (600 m above sea level) Lake Natron is rain shadowed from the Indian Ocean moisture and experiences high air temperatures, resulting in an arid condition with evapotranspiration exceeding precipitation [38].Due to a combination of its aridity and topography, lake levels are not amplified by precipitation events [38], therefore, future water abstraction is a severe threat [36].A planned soda ash mine would severely impact this ecosystem.The notion of soda mining at Lake Natron is not new: in the mid-1820s, 400-500 camel loads amounting to approximately 400 lbs each were reported to have been extracted on an annual basis [39].Saline ecosystems are keystone habitats for specialized endemic fishes such as the Alcolapia species flock (Figure 2A-E, Table 1, Supplementary Video S1), which are confined to the alkaline hydrothermal springs that feed the lake.The lake itself is uninhabitable for fishes.The springs are unique but hostile environments; pH 10.5, salt > 20 ppt, water temperature of 30-42.8• C, and low dissolved oxygen (0.08-6.46 mg/L) [40].The springs provide 25% of the water influx to the lake [41].Lake Natron, a thermal hypersaline lake, is part of a Ramsar Wetland site but is unprotected except for hunting regulations [36].The main water source for this shallow (<3 m) lake is the Southern Ewaso Ng'iro River originating in Kenya.The lake belongs to the Southern Eastern Rift ecoregion classified as having 'globally outstanding' aquatic biological distinctiveness with a 'critical' conservation status [19].The surrounding anthropogenic biomes are predominantly rangeland, with small areas considered 'semi-natural' and 'cropland' [37].Located on the Southern flank of the Kenyan dome in the Gregory Rift at relatively low elevation (600 m above sea level) Lake Natron is rain shadowed from the Indian Ocean moisture and experiences high air temperatures, resulting in an arid condition with evapotranspiration exceeding precipitation [38].Due to a combination of its aridity and topography, lake levels are not amplified by precipitation events [38], therefore, future water abstraction is a severe threat [36].A planned soda ash mine would severely impact this ecosystem.The notion of soda mining at Lake Natron is not new: in the mid-1820s, 400-500 camel loads amounting to approximately 400 lbs each were reported to have been extracted on an annual basis [39].Saline ecosystems are keystone habitats for specialized endemic fishes such as the Alcolapia species flock (Figures 2a-e, Table 1, Supplementary Video S1), which are confined to the alkaline hydrothermal springs that feed the lake.The lake itself is uninhabitable for fishes.The springs are unique but hostile environments; pH 10.5, salt > 20 ppt, water temperature of 30-42.8°C,and low dissolved oxygen (0.08-6.46 mg/L) [40].The springs provide 25% of the water influx to the lake [41].Both Lake Chala and Chemka Springs are part of the Pangani River Basin, which covers approximately 43,650 km 2 [42].Ninety-five percent of the basin is in Tanzania, with five percent in Kenya.Both sites are part of the Pangani ecoregion classified as having 'bioregionally outstanding' aquatic biological distinctiveness with an overall conservation status of 'endangered' [19].At the confluence of the Ruvu and Kikuletwa Rivers (tributaries of the Pangani River) is the Nyumba ya Mungu Dam (NYM) with a 140 km 2 reservoir (Figure 1) [43].Despite one third of the known fishes in the Pangani River Basin being endemic, the ecoregion is considered to be data deficient requiring more comprehensive taxonomic and ecological surveys [44].The Pangani River Basin has been water stressed since the 1940s due to high water losses from irrigation (up to 85%) [42].The ice caps on Mt.Kilimanjaro have decreased by over 75% since 1912; with projected climate change, the Pangani River Basin will no longer receive water input from glacial melt [42].
Chemka Springs is a large volume ground water spring (10 m 3 /s) south of Mount Kilimanjaro near the town of Bomang'ombe (Figure 1, Supplementary Video S1).While the origins of the springs are unknown, the flow remains constant throughout the year.The springs are important water sources for the NYM hydroelectric dam [43].The anthropogenic biomes surrounding the springs consist of rangeland, cropland, villages and dense settlements [37].The high population density in the region has resulted in fires of anthropogenic origins, illegal timber extraction, and overutilization of natural resources [42].The clear, warm, neutral water (28 • C, pH 7.3) of the spring [45] is home to a yet to be described cichlid belonging to the genus Ctenochromis (Figure 2F-J, Table 1, Supplementary Video S1).
Lake Chala is a 98 m deep international volcanic crater lake shared by Tanzania and Kenya on the Eastern side of Mt.Kilimanjaro with a 4.2 km 2 surface area (16.23 km 2 catchment) and an estimated volume of 300-350 Mm 3 [46,47] (Figures 1 and 2K, Supplementary Video S1).The anthropogenic biomes surrounding the lake include residential woodland, remote rangeland, residential rangeland, remote croplands, residential rain fed crops, irrigated villages, rain fed villages, pastoral villages, and mixed settlements [37].Arable land has been increasing in the region since 1973 [48].The lake is fed by groundwater from Mt. Kilimanjaro and receives on average 565 mm of precipitation, exceeded by surface evaporation estimated at 1735 mm per year [46].The water depth is maintained by seepage from precipitation in the forested slopes of Mt.Kilimanjaro and minimal outflow [49]; only 2.5% of the lake volume [50].Permanent anoxic conditions exist below 60 m [49], therefore, only the upper two thirds of the lake are inhabitable by fishes.Due to the high water quality, 7 Mm 3 per year are projected to be abstracted for Kenya's growing population [47,51].Lake Chala is inhabited by Oreochromis hunteri (Günther, 1889), the type species of the genus Oreochromis (Table 1).Fossil evidence from lake sediments suggest that the lake was inhabited by closely related cichlids dating back to ~25 kyr BP [52].Two other cichlids were introduced over the latter half of the 20th century: O. korogwe (Lowe, 1955), by the lake's owner during the colonial era, and Coptodon rendalli (Boulenger, 1896) [52] (Figure 2L).There is uncertainty whether a third cichlid referred to as Haplochromis spec.Chala (Figure 2M,N, Supplementary Video S1) was also introduced [52,53].This is the type species of the genus Oreochromis.It is known from only a few specimens, and little is known about its reproduction and breeding habits.Young O. hunteri feed on algae and organic debris [54].
Tooth morphology suggests aufwuchs (small animals, plants and algae that adhere to open surfaces) form a part of their diet [52].IUCN Red List: Critically endangered.

Remotely Sensed Data
Remotely sensed data (i.e., satellite and unmanned aerial vehicle imagery) were used to determine the LUCC and subsurface water storage for the study areas.The primary remotely sensed data consisted of Landsat 5 TM, 7 ETM+, and 8 OLI imagery downloaded from the USGS Earth Explorer archives (https://earthexplorer.usgs.gov/).Images with less than 15% cloud cover were downloaded for the period spanning 1984-2015.In addition to total cloud cover, each scene was also examined for minimal cloud coverage over the specific study areas.For Landsat 7 ETM+, only images collected prior to 2003 (SLC-ON) were considered to avoid missing data from the scan line corrector malfunction.In total, 19 scenes were downloaded for Lake Chala (7 TM5, 8 ETM+, 4 OLI), 19 for Chemka Springs (9 TM5, 3 ETM+, 7 OLI) and 41 for Lake Natron (15 TM5, 11 ETM+, 15 OLI).
To determine the natural land cover classes at Lake Natron and Lake Chala, a DJI Phantom 2 quadcopter (UAV) was used for 4K aerial videography with a GoPro Hero 4™ (stock fisheye lens 2.92 mm, equal to 16.6 mm on a standard sensor, f/2.8 123 • FOV) in March 2015.Every 15 th frame from the video files was extracted and saved as a full resolution TIFF (Figure 3A,C).For Chemka Springs, a Pleides satellite image from March 1, 2015 was acquired consisting of four 2 m multispectral bands and a 50 cm panchromatic band.The multispectral bands were pan-sharpened with the nearest neighbor diffusion algorithm [55] in ENVI 5.3 (Figure 3B).
The spatial distribution of ground water availability was also assessed from satellite data.The main water influx for all three habitats is from groundwater sources.A study of the world's largest aquifers found that the two largest in Eastern Africa (Ogaden-Juba and the Sudd Basin) have a negative mean annual recharge indicating flow out of the system [56].Both of these large aquifers have also been classified as overstressed systems in terms of their renewable groundwater stress ratio (groundwater use to groundwater availability).During 2002-2006 period, the Eastern African Rift Valley region from Lake Victoria to Lake Malawi experienced a substantial decrease in water levels resulting in the Great Lakes becoming shallower by up to 2 m [57].This period of water depletion was followed by large precipitation events caused by the 2006 Indian Ocean Dipole (IOD) returning the lakes to pre-drought levels by 2007 [58].The IOD is an ocean-atmosphere interaction over the Indian Ocean with alternating positive and negative sea surface temperature anomalies; during periods of positive anomalies, flooding is common in East Africa [58].In order to determine the spatial distribution of ground water availability in the study areas, monthly mass grids (Global mascons) from the Gravity Recovery and Climate Experiment (GRACE) [59,60] were retrieved from the Colorado University GRACE-Interactive Visualization Tool and the GRACE Tellus data portal.The twin satellites launched in 2002 measure monthly changes in total terrestrial water storage through the analysis of observed gravity anomalies [61].The mascons are discrete cells of mass concentration (3-degree resolution) over the surface of the Earth, each with its own gravity signal.Small changes in the gravity signal of the cells are due primarily to changes in the amount of water located in the cell [59].The spatial distribution of ground water availability was also assessed from satellite data.The main water influx for all three habitats is from groundwater sources.A study of the world's largest aquifers found that the two largest in Eastern Africa (Ogaden-Juba and the Sudd Basin) have a negative mean annual recharge indicating flow out of the system [56].Both of these large aquifers have also been classified as overstressed systems in terms of their renewable groundwater stress ratio (groundwater use to groundwater availability).During 2002-2006 period, the Eastern African Rift Valley region from Lake Victoria to Lake Malawi experienced a substantial decrease in water levels resulting in the Great Lakes becoming shallower by up to 2 m [57].This period of water depletion was followed by large precipitation events caused by the 2006 Indian Ocean Dipole (IOD) returning the lakes to pre-drought levels by 2007 [58].The IOD is an ocean-atmosphere interaction over the Indian Ocean with alternating positive and negative sea surface temperature anomalies; during periods of positive anomalies, flooding is common in East Africa [58].In order to determine the spatial distribution of ground water availability in the study areas, monthly mass grids (Global mascons) from the Gravity Recovery and Climate Experiment (GRACE) [59,60] were retrieved from the Colorado University GRACE-Interactive Visualization Tool and the GRACE Tellus data portal.The twin satellites launched in 2002 measure monthly changes in total terrestrial water storage through the analysis of observed gravity anomalies [61].The mascons are discrete cells of mass concentration (3-degree resolution) over the surface of the Earth, each with its own gravity signal.Small changes in the gravity signal of the cells are due primarily to changes in the amount of water located in the cell [59].

Population Growth, Agriculture Expansion and Climate Change Data
The Hyde model [12] was used to determine historical cropland, rainfed agriculture and irrigated agriculture extent and population density from 1950 to 2015.While the model reports the variables at a scale of 5 arc minutes, they were summarized at the GRACE mascon scale (i.e., 5 degrees).Cropland extent refers to arable land and permanent crops [12].Population growth rate projections were extracted from the World Bank for the years 1950-2050 [62].Because the spatial distribution of the population in Tanzania is heterogenous (i.e., there are large cities with high population density and remote locations with little to no population), the World Bank "high" population growth rate estimate for 2050 was scaled according to the local growth rates estimated for the 1950-2015 period calculated from the Hyde model.
Current and projected temperature and precipitation changes for the study areas were obtained from WorldClim v1. 4 [13].For the projected temperature and precipitation for the year 2050, downscaled and calibrated NASA Goddard Institute for Space Studies (GISS) global climate model projections (RCP 8.5) using WorldClim v1.4 as the baseline were chosen [63,64].Representative Concentration Pathway 8.5 represents high future greenhouse gas emissions.Characteristics of a future under RCP 8.5 include three times the present-day CO 2 emissions by 2100, rapid methane emissions, an increased use of croplands and grasslands due to an increase in the human population, and heavy reliance on fossil fuels [65].Some key differences between RCP 8.5 and intermediate scenarios are that the intermediate scenarios assume either a decrease in the use of croplands and grasslands due to dietary changes or a decrease in CO 2 emission following 2060 [65].While not a forecast, RCP 8.5 represents a business as usual scenario and currently, emissions have been shown to be consistent with this scenario [66].Contemporary surface temperature anomalies calculated from monthly instrumental temperature records were retrieved from the HadCRU4 global dataset for the 1931-2014 period [14,67].

Image Analysis
Cloud and cloud shadow masks were calculated with Fmask 3.2 [68] implemented in Matlab v. 2015a (Mathworks, Natick, MA, USA).Because the study areas are substantially smaller than an entire Landsat scene, each image and corresponding cloud mask were subset to an area of 1278 km 2 for Lake Natron, 491 km 2 for Lake Chala, and 38 km 2 for Chemka Springs.Total cloud area was recalculated for the subsets and only images with less than 3% cloud cover were retained, resulting in 33 scenes for Lake Natron, 19 for Lake Chala and 15 for Chemka Springs (Figure 4).Individual cloud masks were combined to create one comprehensive mask spanning the entire period  for each study area.From the comprehensive cloud masks, there was 89.5% cloud free area for Lake Natron, 97.3% for Chemka Springs, and 85.3% for Lake Chala on which the remainder of the analyses were based (Figure 5).
The three Landsat datasets (5 TM, 7 ETM+, 8 OLI) have different radiometric and spectral characteristics, therefore, to make the outputs comparable over the multiple image dates and sensors, a tasseled cap transformation was applied to the cloud-free pixels of the images calibrated to top-of-the-atmosphere reflectance in ENVI 5.3 [69].The tasseled cap transformation compresses and decorrelates the multispectral bands into three bands orthogonally distributed along new axes associated with brightness, greenness, and wetness [70,71].The brightness component is associated with the presence of soil and highly reflective natural and manmade features.The greenness component exploits the contrast between the red and near infrared bands commonly seen with vegetation.The third component, wetness, is associated with soil moisture and other 'wet' features in the landscape.
Because several natural and anthropogenic land use classes in the study areas have similar land cover components (e.g., various proportions of exposed soil and vegetation) a time sequence of the tasseled cap transformed image subsets was made for the 1984-2015 period for each study area (Figure 6).The premise of examining the land use trajectory of the pixels over time in the time sequence rather than classifications conducted on single images is that pixels representing natural vegetation will have a spectral phenology that is different from pixels representing land cover modified by human intervention (e.g., irrigated agriculture) during the time sequence.By examining a time sequence of several years the goal was to locate areas of natural vegetation within a complex land cover matrix with a variety of small-scale anthropogenic modifications (Figure 6).
Remote Sens. 2017, 9, 623 8 of The Hyde model [12] was used to determine historical cropland, rainfed agriculture and irrigated agriculture extent and population density from 1950 to 2015.While the model reports the variables at a scale of 5 arc minutes, they were summarized at the GRACE mascon scale (i.e., degrees).Cropland extent refers to arable land and permanent crops [12].Population growth rate projections were extracted from the World Bank for the years 1950-2050 [62].Because the spatial distribution of the population in Tanzania is heterogenous (i.e., there are large cities with high population density and remote locations with little to no population), the World Bank "high" population growth rate estimate for 2050 was scaled according to the local growth rates estimated for the 1950-2015 period calculated from the Hyde model.
Current and projected temperature and precipitation changes for the study areas were obtained from WorldClim v1. 4 [13].For the projected temperature and precipitation for the year 2050, downscaled and calibrated NASA Goddard Institute for Space Studies (GISS) global climate model projections (RCP 8.5) using WorldClim v1.4 as the baseline were chosen [63,64].Representative Concentration Pathway 8.5 represents high future greenhouse gas emissions.Characteristics of a future under RCP 8.5 include three times the present-day CO2 emissions by 2100, rapid methane emissions, an increased use of croplands and grasslands due to an increase in the human population, and heavy reliance on fossil fuels [65].Some key differences between RCP 8.5 and intermediate scenarios are that the intermediate scenarios assume either a decrease in the use of croplands and grasslands due to dietary changes or a decrease in CO2 emission following 2060 [65].While not a forecast, RCP 8.5 represents a business as usual scenario and currently, emissions have been shown to be consistent with this scenario [66].Contemporary surface temperature anomalies calculated from monthly instrumental temperature records were retrieved from the HadCRU4 global dataset for the 1931-2014 period [14,67].

Image Analysis
Cloud and cloud shadow masks were calculated with Fmask 3.2 [68] implemented in Matlab v. 2015a (Mathworks, Natick, MA, USA).Because the study areas are substantially smaller than an entire Landsat scene, each image and corresponding cloud mask were subset to an area of 1278 km for Lake Natron, 491 km 2 for Lake Chala, and 38 km 2 for Chemka Springs.Total cloud area was recalculated for the subsets and only images with less than 3% cloud cover were retained, resulting in 33 scenes for Lake Natron, 19 for Lake Chala and 15 for Chemka Springs (Figure 4).Individual cloud masks were combined to create one comprehensive mask spanning the entire period  for each study area.From the comprehensive cloud masks, there was 89.5% cloud free area for Lake Natron, 97.3% for Chemka Springs, and 85.3% for Lake Chala on which the remainder of the analyses were based (Figure 5).The three Landsat datasets (5 TM, 7 ETM+, 8 OLI) have different radiometric and spectral characteristics, therefore, to make the outputs comparable over the multiple image dates and sensors, a tasseled cap transformation was applied to the cloud-free pixels of the images calibrated to top-ofthe-atmosphere reflectance in ENVI 5.3 [69].The tasseled cap transformation compresses and decorrelates the multispectral bands into three bands orthogonally distributed along new axes associated with brightness, greenness, and wetness [70,71].The brightness component is associated with the presence of soil and highly reflective natural and manmade features.The greenness component exploits the contrast between the red and near infrared bands commonly seen with In order to determine the natural land cover in each study area, the high spatial resolution UAV images and pan sharpened Pleides image were examined in addition to GPS tagged field photographs collected in March 2015.For Lake Natron the natural land cover classes include exposed soil, trona, grassland, aquatic grasses, and bushland (Figure 3A).For Chemka Springs, 'closed canopy tree stand' (Riparian forest) is the natural land cover class (Figure 3B).The primary natural land cover class around Lake Chala is bushland (Figure 3C).
(Figure 6).The premise of examining the land use trajectory of the pixels over time in the time sequence rather than classifications conducted on single images is that pixels representing natural vegetation will have a spectral phenology that is different from pixels representing land cover modified by human intervention (e.g., irrigated agriculture) during the time sequence.By examining a time sequence of several years the goal was to locate areas of natural vegetation within a complex land cover matrix with a variety of small-scale anthropogenic modifications (Figure 6).In order to determine the natural land cover in each study area, the high spatial resolution UAV images and pan sharpened Pleides image were examined in addition to GPS tagged field photographs collected in March 2015.For Lake Natron the natural land cover classes include exposed soil, trona, grassland, aquatic grasses, and bushland (Figure 3a).For Chemka Springs, 'closed canopy tree stand' (Riparian forest) is the natural land cover class (Figure 3b).The primary natural land cover class around Lake Chala is bushland (Figure 3c).
The high spatial resolution data sets (Pleides and UAV), geotagged field photographs and Google Earth were used to select training and validation points for classification.All training and validation points were located at least 90 m apart.Prior to classification the bands in the time sequence were reduced with the BandMax algorithm with the natural class(es) as the target class(es).This algorithm provides a subset of bands that maximize the separability between target(s) and the background.The reduced time sequences were then classified with Spectral Angle Mapper [72].This algorithm determines the similarity between two spectra by calculating the angle between them; smaller angles represent spectra that are a better match to the reference provided.In this case, the reference was the reduced time sequence for the 'natural' land cover class(es).The first time sequence analyzed for all sites spanned 1984-2015.Areas classified as 'natural' in this time sequence are expected not to have changed over the 31 years and represent the 2015 extent of the natural land The high spatial resolution data sets (Pleides and UAV), geotagged field photographs and Google Earth were used to select training and validation points for classification.All training and validation points were located at least 90 m apart.Prior to classification the bands in the time sequence were reduced with the BandMax algorithm with the natural class(es) as the target class(es).This algorithm provides a subset of bands that maximize the separability between target(s) and the background.The reduced time sequences were then classified with Spectral Angle Mapper [72].This algorithm determines the similarity between two spectra by calculating the angle between them; smaller angles represent spectra that are a better match to the reference provided.In this case, the reference was the reduced time sequence for the 'natural' land cover class(es).The first time sequence analyzed for all sites spanned 1984-2015.Areas classified as 'natural' in this time sequence are expected not to have changed over the 31 years and represent the 2015 extent of the natural land cover class(es).Additional time sequences spanning 1984-1987, 1984-1995, and 1984-2000 were analyzed for Lake Chala and 1984-1995 for Chemka Springs.There were not enough images available for Chemka Springs to examine the 1984-1987 and 1984-2000 time sequences.For Lake Natron, no additional time sequences beyond the 1984-2015 period were examined because the entire study area consisted of natural land cover classes with the exception of an agricultural area and Maasai villages (~90 km 2 ) that remained consistent over time (confirmed during the field visit in March 2015).
For an expected accuracy of 85% and allowable error of 5% the minimum number of validation points required to assess the classification accuracy was calculated to be 203 [73].A stratified random sampling scheme was used to assign the points to the two classes of interest, 'natural' and 'modified' for each study area separately.

Fragmentation Analysis
Fragmentation analysis describes the spatial configuration of the landscape classes of interest.Fragstats 4.1 was used to calculate patch and class level statistics [74].Evidence suggests that reporting solely the total area of natural vegetation in a catchment is insufficient and that the spatial configuration and quality of the patches of natural vegetation are also important considerations [30,[75][76][77].At the patch level, we analyzed the distribution of the number of natural vegetation fragments for each period as a proxy for measuring landscape composition and configuration.At the class level the following metrics were calculated: total class area (CA), percentage of the landscape (PLAND), largest patch index (LPI), landscape shape index (LSI), as well as the clumpiness (CLUMPY) index for the 'natural' class.The indices CA and PLAND measure the total amount of the class in hectares and the percentage of that class in the landscape, respectively.LPI is a measure of the percentage of the landscape comprised by the largest patch.LSI quantifies the perimeter-to-area ratio for the landscape as a whole, therefore, it is a proxy for the overall geometric complexity of the 'natural' class; moreover, it can be interpreted as a measure of patch dispersion.The CLUMPY index measures the degree of distribution of patches in the landscape as a measure of landscape configuration.It ranges from −1 for maximally disaggregated patches to 1 for maximally computed patches, while a value of zero indicates a random distribution.Because the fragmentation statistics are based on a single land cover class (i.e., 'natural'), we did not address landscape level metrics.

Land Cover Change
Based on the classification of the satellite imagery for Lake Chala, a decrease in natural vegetation cover (bushland) from 15.5% (1987) to 3.5% in 2015 was observed (Figure 7).Calculation of the bushland extent within the smaller catchment area (16.23 km 2 ) revealed a decline in total extent from 39.2% (1987) to 16% (2000) and a slight increase to 18.6% in 2015.User's accuracy for the bushland class ranged from 83% to 88% for the four time sequence classifications.Producer's accuracy ranged from 88% to 95% for the four classifications.At the patch level, results indicate a large number of small bushland fragments (<10 ha), two medium fragments (10-100 ha) and a very large fragment (6200 ha) in the landscape in 1987 (Table 2, Figures 7 and 8).As land degradation increased in the area, by 1995 we found a reduction in the number of small bushland fragments and an increase in medium sized fragments as a result of the reduction of the largest fragment (Figures 7 and 8).This also resulted in a substantial reduction in the area of the core bushland fragment to 1520 ha.A large increase in the number of bushland fragments from 203 to 526 (Table 2), and a decrease in the mean patch area from 18.2 ha to 4.93 ha (Table 2) was found by the year 2000.For the most recent time period (2015), we found a steep reduction in the number of bushland fragments to 62 indicating that only a few fragments are left in the landscape around the lake.
At the class level, the fragmentation statistics are consistent with patch level statistics, showing a consistent reduction in CA and PLAND.For instance, the total bushland cover decreased from 7633 ha from the initial time period (1987) to 1726 ha by 2015, representing a change in PLAND from 15.5% to 3.5% (Table 2).The LPI, which represents the percentage of the landscape comprised by the largest patch, consistently decreases in the core area around Lake Chala from 12.73 in 1987 to 2.26 in 2000, with a slight increase to 2.88 in 2015.The LSI for Lake Chala shows a peak value in 2000 (28.81) with a decrease to 11.49 by 2015.Finally, the CLUMPY index ranges from 0.83 in 2000 to 0.90 in 1987, indicating a decreased aggregation of the bushland class.For Chemka Springs, the natural vegetation (closed canopy forest) decreased to 3.5% (2015) from an initial extent of 7.4% (1995).The user's accuracy for the natural vegetation ranged from 88% to 94% for the two classifications.Similarly, the producer's accuracy ranged from 80% to 96%.A change in the landscape structure based on patch size distribution can be seen (Figure 8).The results from 1995 illustrate a large number of fragments smaller than 15 ha, however, by 2015 most of the fragments are less than 5 ha, with a lower number of patches overall (194 vs. 104, respectively) and a reduction in the mean patch area (1.45 ha to 0.92 ha).At the class level, there is a decrease from 281.8 ha to 95.7 ha (CA) as well as a decrease in PLAND from 7.4% to 2.5%.This pattern is more dramatic in the LPI, which decreases from 1.19% to 0.25%.Overall, the CLUMPY values are lower than Lake Chala for both periods, with 0.61 and 0.52 for 1995 and 2015, respectively, indicating a more disaggregated landscape than that found around Lake Chala.

GRACE Satellite Subsurface Water Data
The assessment of the GRACE satellite data indicates a sharp decreasing trend of water availability for both the Lake Natron (−4.4 cm/yr) mascon and the Lake Chala and Chemka Springs (−1.69 cm/yr) mascon over the 2002-2006 period followed by a rapid increase in 2007 to a maximum value of 14.8 cm of water equivalent height for Lake Natron and 10.8 cm for Lake Chala and Chemka Springs (Figures 9a,b).The overall trend over the observation period of the GRACE satellites from 2002 to 2015 is 0.82 cm/yr for the Lake Natron mascon and −0.81cm/yr for the Lake Chala and Chemka

GRACE Satellite Subsurface Water Data
The assessment of the GRACE satellite data indicates a sharp decreasing trend of water availability for both the Lake Natron (−4.4 cm/yr) mascon and the Lake Chala and Chemka Springs (−1.69 cm/yr) Remote Sens. 2017, 9, 623 14 of 25 mascon over the 2002-2006 period followed by a rapid increase in 2007 to a maximum value of 14.8 cm of water equivalent height for Lake Natron and 10.8 cm for Lake Chala and Chemka Springs (Figure 9A,B).The overall trend over the observation period of the GRACE satellites from 2002 to 2015 is 0.82 cm/yr for the Lake Natron mascon and −0.81cm/yr for the Lake Chala and Chemka Springs mascon.The mascon for Lake Natron also includes the southeastern edge of Lake Victoria, the largest African lake by surface area (68,800 km 2 ).Therefore, any decrease in subsurface water near the Lake Natron springs specifically may be obscured by the larger gravity signal from Lake Victoria (Figure 9C).Springs mascon.The mascon for Lake Natron also includes the southeastern edge of Lake Victoria, the largest African lake by surface area (68,800 km 2 ).Therefore, any decrease in subsurface water near the Lake Natron springs specifically may be obscured by the larger gravity signal from Lake Victoria (Figure 9c).

Population and Climate Model Results
The changes in cropland extent and rainfed and irrigated agriculture extents summarized at the mascon scale (Figure 10) illustrate how a larger region around the catchments has changed over the 1980-2015 period providing insight into potential future stressors caused by the conversion of natural land cover to agricultural activities.An increase in the area of cropland and rainfed agriculture can be seen for both mascons.Irrigated agriculture also increased for the Chala/Chemka mascon.While covering the smallest total area of the three categories, the irrigated agriculture in the Chala/Chemka mascon for the entire 1980-2015 period is concentrated around Chemka Springs, increasing from 16.7 to 34.0 km 2 immediately surrounding the Ctenchromis habitat.A large increase in cropland and rainfed agriculture can also be seen for the Natron mascon between 2010 and 2015.While there is no increase in agriculture within the Natron basin, encroachment can be seen east of the Gelai volcano

Population and Climate Model Results
The changes in cropland extent and rainfed and irrigated agriculture extents summarized at the mascon scale (Figure 10) illustrate how a larger region around the catchments has changed over the 1980-2015 period providing insight into potential future stressors caused by the conversion of natural land cover to agricultural activities.An increase in the area of cropland and rainfed agriculture can be seen for both mascons.Irrigated agriculture also increased for the Chala/Chemka mascon.While covering the smallest total area of the three categories, the irrigated agriculture in the Chala/Chemka mascon for the entire 1980-2015 period is concentrated around Chemka Springs, increasing from 16.7 to 34.0 km 2 immediately surrounding the Ctenchromis habitat.A large increase in cropland and rainfed agriculture can also be seen for the Natron mascon between 2010 and 2015.While there is no increase in agriculture within the Natron basin, encroachment can be seen east of the Gelai volcano (located SE of Lake Natron) and on the western slopes of the Gregory Rift (western border of the Natron basin).Based on the World Bank high population growth model, Tanzania will experience an increase in population density from 54.9 people/km 2 (1950) to 170.7 people/km 2 by 2050.Regional population density estimates project 102 people/km 2 for Lake Natron, 245 people/km 2 for Chemka Springs, and 249 people/km 2 for Lake Chala.Over the period of this study, the population density tripled for both of the two mascon areas (Figure 10).The total number of inhabitants increased from 2.6 to 7.2 million and 2.3 to 6.2 million for the 1980-2015 period for the Chala/Chemka and Natron mascon regions respectively.
For the three sites the HadCRUT4 data illustrate an increasing trend in the 20-year average surface temperature anomalies since 1968 by 1.6°C.By 2050, the GISS (RCP 8.5) model predicts an increase in temperature of 2.5°C and precipitation decrease of 15 mm/yr for Lake Natron, +1.8°C and −4.6 mm/yr for Chemka Springs and +1.7°C and −3.9 mm/yr for Lake Chala.Based on the World Bank high population growth model, Tanzania will experience an increase in population density from 54.9 people/km 2 (1950) to 170.7 people/km 2 by 2050.Regional population density estimates project 102 people/km 2 for Lake Natron, 245 people/km 2 for Chemka Springs, and 249 people/km 2 for Lake Chala.Over the period of this study, the population density tripled for both of the two mascon areas (Figure 10).The total number of inhabitants increased from 2.6 to 7.2 million and 2.3 to 6.2 million for the 1980-2015 period for the Chala/Chemka and Natron mascon regions respectively.
For the three sites the HadCRUT4 data illustrate an increasing trend in the 20-year average surface temperature anomalies since 1968 by 1.6 • C. By 2050, the GISS (RCP 8.5) model predicts an increase in temperature of 2.5 • C and precipitation decrease of 15 mm/yr for Lake Natron, +1.8 • C and −4.6 mm/yr for Chemka Springs and +1.7 • C and −3.9 mm/yr for Lake Chala.

Discussion
The LUCC history around aquatic habitats has been shown to have an important influence on the present-day fish diversity and habitat conditions [10,78].However, the classification of satellite imagery in the arid ecosystems examined here is inherently challenging.In an assessment of LUCC for the slopes of Mt.Kilimanjaro, it was indicated that from single date imagery, bushland, grassland and cropland cannot be readily separated due to the transitional nature of the classes and their spectral similarity [48].Furthermore, vegetation dynamics in the region are strongly influenced by cultivation and land conversion [16].Small-scale agriculture, with variable irrigation practices can be difficult to separate from the natural vegetation phenological cycles in arid regions [79].However, separating the natural landscape from modified areas is needed in order to understand trends in LUCC dynamics [69].The tasseled cap time sequence of the historical Landsat imagery allowed us to separate out the natural land cover classes from the modified areas and reduce the confusion between spectrally similar classes at single points in time.The results also illustrate an isolation and fragmentation of the natural vegetation around Lake Chala and Chemka Springs.In both cases the total native vegetation area was found to be below the recommended 70% for the maintenance of water quality [30].The width of the Riparian buffer around Chemka springs is below the recommended thresholds for aquatic wildlife protection (30 m), and is at the lower limit of the thresholds for stream temperature stabilization (10 m), nutrient retention (15 m), sediment control (10 m), and pesticide retention (15 m) [80].However, when land use intensity is considered, the minimum buffer recommended to moderate water temperature and quality increases to 30 m [81].In contrast, for Lake Chala the vegetated slopes of the crater and the bushland vegetation on the Tanzanian side afford some protection but this is greatly reduced on the Kenyan side of the lake where, outside of the crater the natural vegetation has been removed.If the rate of natural vegetation loss were to remain constant, the remaining fragments of bushland in the Lake Chala catchment would be converted to other land covers by 2040 and the remaining closed canopy forest around Chemka Springs would be lost by 2033.Mitigation and management efforts could in the long-term to protect the remaining natural vegetation and potentially also to allow it to regrow.
The rapid decline of groundwater in the region followed by the increase in water availability can be seen in the mascon grid time series for our study areas (Figure 9A,B).With the rainfall patterns (and associated vegetation phenology) in East Africa closely associated with IOD and El Nino Southern Oscillation (ENSO) events, further depletion events in terrestrial water storage can be expected [16,82].Depletion events would be important for small aquatic habitats where subsurface sources (e.g., springs) are the primary source of incoming water.Irrigation and other water use such as drinking water in the region almost entirely draws water from sources originating from groundwater [42,47].With precipitation being a significant driver of the water storage behavior of the region encompassing our study areas [83], climate change scenarios for the region suggest an increase in water stress in the future.In addition, interannual high precipitation events associated with large-scale circulations (i.e., IOD, ENSO) will play increasingly important roles in mitigating accumulated water storage deficits [83].A global assessment by Reference [84] found high correlations between threats to human water security and biodiversity, with the strongest correlations in heavily settled regions indicating the negative influence a growing population has on aquatic ecosystems.
The continent-wide approach by Reference [19] assessed ecoregions based on land degradation and threats to aquatic habitats.Our approach focuses on a finer scale by considering current and future stressors to the specific aquatic habitats of cichlid species taking into account the global and regional drivers of ecosystem vulnerability recommended by Reference [85].Of the habitats examined here, we found that the Alcolapia species flock restricted to the Lake Natron Springs was vulnerable to forecasted climate change and potential future landscape alteration.The water in the lake itself is uninhabitable for fish to survive, even a short period of time.Therefore, if the temperature of the shallow springs surpasses their thermal tolerance or the water conditions are altered by future water abstraction [86,87], these fish have no migratory recourse.While the total area of the springs is estimated to be less than 10 km 2 , it is unknown how much of that is inhabited by the Alcolapia.A planned contemporary soda ash mine would result in the construction of dikes and other infrastructure, as well as an increase in the local population.Such mines require large amounts of water, and any impact on water levels could impact the Alcolapia populations that rely on the shallow springs [36].A hydropower dam proposed for the Ewaso Ngiro River in Kenya would further decrease the main source of water input to the lake from the northern side [36], putting additional pressure on the southern freshwater springs.
The highest upper critical temperature (45.6 • C) recorded for a freshwater fish was for Alcolapia grahami (Boulenger 1912), from nearby Lake Magadi in Kenya [87], which is less than 3 • C warmer than the current upper range of the water temperature for the Natron Springs Alcolapia.The only other cichlids known to naturally inhabit similar temperatures are the endemic cichlids from the genus Danakilia, found in the shallow (<10 m) waters of Lake Abaeded in the Danakil Depression of Northern Eritrea and Lake Afrera in Ethiopia [88], and Iranochicla from Southern Iran [89].The water influx to the lakes inhabited by Danakilia is also mainly from springs and ranges in temperature from 29-45 • C [88].Nevertheless, the Danakilia, while abundant in lakes, are generally found in the coolest areas, and rarely in the upper extremes of water temperature, suggesting that while these cichlids may be thermal tolerant, if conditions allow the fish will reside in the cooler temperature range of their habitat.The projected air temperature increase of 2.5 • C (by 2050) from the GISS model for the Natron Springs Alcolapia is consistent with contemporary surface temperature increases observed from ground monitoring stations (HadCRUT4 data).While the high specific heat capacity of water will protect the fish from abrupt temperature changes, we believe that projected climate change will subject the Alcolapia to prolonged exposures to warmer conditions.As ectotherms, fish have poor insulation and their body temperature is predominantly controlled by the water body in which they live [90].Short-term positive impacts on growth of increases in environmental temperature are negated if the temperature surpasses the fish's optimal range, beyond which growth and other necessary biological processes are compromised [91].Tropical fish, in general, experience narrower temperature fluctuations than temperate species and many already live close to their maximal thermal tolerance [92,93].With climate change, fish are expected to be exposed to thermal regimes outside their thermal tolerance windows [94].
Based on the 2006 the IUCN Red List assessment A. alcalicus is considered endangered with a decreasing population trend.One of the primary causes listed for the imperilment of A. alcalicus is climate change.Given that A. ndalalani and A. latilabris are more specialized than A. alcalicus (Table 1), we expect them to be more threatened by changes in their habitat [95,96].Both A. ndalalani and A. latilabris are categorized as vulnerable by the IUCN Red List; the primary threat listed for both is habitat degradation.A projected decrease in precipitation and a potential decrease in the subsurface aquifer could lead to an overall loss of habitat extent or depth.All three species graze on algae and biocover growing on the rocks and stems of marginal plants.The marginal grasses also act as a nursery for juvenile fishes.A decrease in water depth in the springs could have devastating effects on the available habitat including a reduced amount of marginal vegetation (Juncus sp.) that acts both as a vital habitat for the young cichlids, once released from the mouth of the female, cover from cannibalistic larger specimens, and extended surface area for algae and biocover.The projected population growth (average 102 inhabitants/km 2 ) by 2050 is further expected to stress the water resources of the arid region.
Our results also highlight the importance of the management of the local water resources at Chemka Springs.Chemka Springs contribute a major influx of water to the NYM dam reservoir.All natural vegetation except the Riparian zone immediately next to the springs have been removed (Figures 2B and 7F).Palm plantations were observed to be encroaching on the springs.With an increase in population density, cropland, and agriculture (Figure 10), the need for water abstraction also increases.Irrigated agriculture doubled in extent in the immediate vicinity of the springs; this is expected to further increase as rainfall and subsurface water may become inadequate for the growing needs of the population.Pollution from agricultural runoff may further adversely affect the limited habitat of undescribed Ctenochromis sp.Due to potential fluvial habitat fragmentation and breaks in connectivity, it is unclear how much area would be available for the fishes to migrate if their current habitat at the mouth of the springs becomes uninhabitable.The site (springhead) is a popular swimming spot for both locals and tourists, but no formal management of the site exists to protect it from degradation, refuse left behind by visitors and other activities that are potentially detrimental to the habitat (e.g., washing vehicles in the shallow portions of the stream).
Arguably, the habitat of O. hunteri in Lake Chala may undergo the least amount of change unless there is large water abstraction from Lake Chala, increased agricultural runoff, or introduction of additional species.The large volume of Lake Chala provides some protection against climate change.However, increased water abstraction due to climate change, introduced species and water quality degradation from sedimentation and agricultural runoff are important detrimental factors.Based on a 2006 assessment, the IUCN Red List categorizes O. hunteri as critically endangered with a decreasing population.The primary threats listed by the IUCN are climate change and habitat degradation.Another threat to O. hunteri is the introduction of species, especially C. rendalli that compete directly for food and territory.The high reproduction rate of C. rendalli also outnumbers the slower reproduction rate of the mouth brooding Oreochromis.The wide range of tolerated environmental conditions, trophic adaptability and high reproductive rates have been shown to be predominant factors that have allowed tilapias (such as C. rendalli) to become one of the groups of exotics fishes most widely distributed worldwide, exacerbating damage to native fish, and making them a primary cause of species endangerment and extinction [97].The voracity of these tilapias, especially during times of year where food resources are scarce leads to their ability to outcompete native fishes [98].
One of the limitations in this study is the lack of historical or contemporary population data for the species examined.However, as described earlier, such data on African fishes outside the Great Lakes in small aquatic systems are non-existent for most of the endemic species.Therefore, our conclusions on the putative effects of environmental and local human induced stressors are based upon studies of other fishes worldwide undergoing similar stressors [99][100][101][102][103][104][105] among many others.A well-studied example of a fish critically endangered by climate change and water abstraction for human uses is the Devil's Hole pupfish (Cyprinodon diabolis Wales, 1930).It lives in an underground cave with water from a single aquifer-fed thermal pool in the Desert National Wildlife Refuge Complex in Nye County, Nevada [106].By 2013, the wild population was recorded to have decreased to less than 40 individuals [107].It was among the first species to be listed on the U.S. Endangered Species Act in 1967 and was the focus of a legal battle over water rights, the culmination of which was a U.S. Supreme Court decision to cease groundwater abstraction for development.Despite this protection, populations of C. diabolis have declined since 1995; one of the potential causes being increased water temperatures caused by climate change, and therefore, decreased dissolved oxygen.Projected climate change by 2050 was found to compromise the recruitment window for juvenile fish, decrease egg viability and increase larval mortality [107].
While not examined here, increasing atmospheric deposition of sulphur, nitrogen, mercury and other compounds has been shown to have adverse effects on a variety of ecosystems including lakes and streams [108][109][110][111]. Future analyses should incorporate modeled rates of atmospheric pollutant deposition on the habitats.
Another recommendation is to increase the scope of the analyses to the continental scale with a focus on small aquatic habitats outside the Great Lakes.In a systematic review of studies examining predicted and observed effects of climate change on freshwater fishes, Reference [112] found that tropical arid and semi-arid regions have been grossly understudied.Future predictions were heavily biased to the Nearctic (~75%), while the Palaearctic dominated studies examining the effect of recent climate change (~46%).Over half the published studies included at least one salmonid species, indicating a taxonomic bias.The recently recognized Ovalentaria [113,114] includes a number of the fish groups previously classified within the "perciformes" of previous authors.As currently recognized, the Ovalentaria includes an array of marine, and a few freshwater groups, among which the Cichlidae represent a major component of tropical freshwater diversity.Despite this taxonomic importance, Reference [112] highlighted a general lack of studies from tropical freshwaters.Furthermore, less than 10% of the studies focused on fishes categorized as 'threatened with extinction' by the IUCN.Based on the IUCN Red List version 3.1, of the species assessed, 1045 freshwater fish are considered critically endangered or endangered.Of those, climate change is a primary threat for 185 species; six are cichlids, with Cyprinidae accounting for 48.6%.When considering human-induced stressors such as agriculture, energy production, logging/wood harvesting, fishing, human intrusions, natural system modification, and pollution, 838 species have been listed as affected.Of these 143 (17.1%) are cichlids.When invasive species are considered as a major threat, 336 species have been assessed as critically endangered or endangered, of which 76 (22.6%) are cichlids.Native African species account for 33.1% of critically endangered or endangered species for which climate change, human-induced stressors, or invasive species are primary threats.

Conclusions
Humanity's aggregate resource consumption has grown by 78% since 1970, leading to a sustained increase in pressure on global biodiversity [115].For all types of freshwater ecosystems, strong negative relationships have been shown between biodiversity intactness and stressors such as human impact (e.g., LUCC) and climate change [116].Rapid population growth in the regions surrounding the habitats studied here will lead to increased stress on the aquatic ecosystems through water abstraction and pollution of the surface waters.The minimal buffers of natural vegetation surrounding Lake Chala and Chemka Springs are insufficient to moderate water quality degradation from agricultural runoff in the long term.For Chemka Springs with the high intensity agricultural land use surrounding the springs the vegetative buffer is further insufficient to moderate water temperature changes due to projected climate change scenarios.The region encompassed by all three study areas is likely to become increasingly reliant on interannual high precipitation events to mitigate accumulated water storage deficits.
Without protection, the observed trends of decreasing natural vegetation surrounding the aquatic habitats are likely to intensify throughout Tanzania due to large population growth and resulting increases in food demands [117].Deducing the main stressors acting upon the habitats of aquatic species may allow for the development of mitigation measures, the refinement of regional development plans, or the adoption of new paradigms in the management and protection of aquatic habitats [118].With proper management and conservation strategies these habitats could become hydrologic refugia for endemic aquatic species [119].
We reiterate the importance of considering the unique landscape conditions surrounding relatively small ecosystems when conducting remote sensing analyses.Global models of land cover are often unable to take into account the challenges of mapping land cover at such a small scale.

Figure 1 .Figure 1 .
Figure 1.Map of the study areas in Northern Tanzania.Blue circles represent the locations of the study areas at Lake Natron Springs, Chemka Springs, and Lake Chala.Other important points of Figure 1.Map of the study areas in Northern Tanzania.Blue circles represent the locations of the study areas at Lake Natron Springs, Chemka Springs, and Lake Chala.Other important points of references such as Mt.Kilimanjaro, Mt.Meru, Lake Manyara, the town of Bomang'ombe, Nyumba ya Mungu Dam and the city of Arusha are also illustrated.

Figure 4 .
Figure 4. Dates of Landsat TM, ETM+, and OLI used in the analysis.Figure 4. Dates of Landsat TM, ETM+, and OLI used in the analysis.

Figure 4 .
Figure 4. Dates of Landsat TM, ETM+, and OLI used in the analysis.Figure 4. Dates of Landsat TM, ETM+, and OLI used in the analysis.

Figure 6 .
Figure 6.Example of the tasseled cap time sequence for the Lake Chala study area.The differences between the natural land cover class (bushland) and modified classes of agriculture and mixed land use are illustrated.

Figure 6 .
Figure 6.Example of the tasseled cap time sequence for the Lake Chala study area.The differences between the natural land cover class (bushland) and modified classes of agriculture and mixed land use are illustrated.

Figure 8 .
Figure 8. Natural vegetation patch size distributions for Lake Chala and Chemka Springs.

Figure 8 .
Figure 8. Natural vegetation patch size distributions for Lake Chala and Chemka Springs.

Figure 9 .
Figure 9. GRACE satellite derived time series of water equivalent height (cm) for the mascon encompassing (a) Lake Natron and; (b) Chemka Springs and Lake Chala; (c) Mascon grid cells representing water equivalent height from March 2015.

Figure 9 .
Figure 9. GRACE satellite derived time series of water equivalent height (cm) for the mascon encompassing (A) Lake Natron and; (B) Chemka Springs and Lake Chala; (C) Mascon grid cells representing water equivalent height from March 2015.
Remote Sens. 2017, 9, 623 16 of 26(located SE of Lake Natron) and on the western slopes of the Gregory Rift (western border of the Natron basin).

Figure 10 .
Figure 10.Changes in cropland extent and agriculture (rainfed and irrigated) for the (a) Natron and (b) Chemka/Chala mascons.Mean population density for the two mascons is also shown.The area covered by the mascons is illustrated in Figure 9.

Figure 10 .
Figure 10.Changes in cropland extent and agriculture (rainfed and irrigated) for the (A) Natron and (B) Chemka/Chala mascons.Mean population density for the two mascons is also shown.The area covered by the mascons is illustrated in Figure 9.
such as Mt.Kilimanjaro, Mt.Meru, Lake Manyara, the town of Bomang'ombe, Nyumba ya Mungu Dam and the city of Arusha are also illustrated. references

Table 1 .
Descriptions of endemic cichlids from the three study areas.

Table 1 .
Cont.Potential area of occupancy <1 km 2 .Closely related to C. pectoralis (Pfeffer, 1893), described from the Pangani River further downstream.Up to 12 cm TL, mouth brooding, feeds on snails and micro invertebrates.Dominant males excavate large crater nests (up to 60 cm in diameter) under large boulders.IUCN Red List: Not assessed.

Table 2 .
Natural vegetation fragmentation statistics for Lake Chala and Chemka Springs.Values for total class area (CA), percentage of the class in the landscape (PLAND), largest patch index (LPI), landscape shape index (LSI), and clumpiness (CLUMPY) are shown.The 'natural' land cover classes examined are bushland for Lake Chala and closed canopy tree stand for Chemka Springs.