Quantifying Land Use Land Cover Changes in the Lake Victoria Basin Using Satellite Remote Sensing: The Trends and Drivers between 1985 and 2014

The Lake Victoria Basin (LVB) is a significant resource for five states within East Africa, which faces major land use land cover changes that threaten ecosystem integrity and ecosystem services derived from the basin’s resources. To assess land use land cover changes between 1985 and 2014, and subsequently determine the trends and drivers of these changes, we used a series of Landsat images and field data obtained from the LVB. Landsat image pre-processing and band combinations were done in ENVI 5.1. A supervised classification was applied on 118 Landsat scenes using the maximum likelihood classifier in ENVI 5.1. The overall accuracy of classified images was computed for the 2014 images using 124 reference data points collected through stratified random sampling. Computations of area under various land cover classes were calculated between the 1985 and 2014 images. We also correlated the area from natural vegetation classes to farmlands and settlements (urban areas) to explore relationships between land use land cover conversions among these classes. Based on our land cover classifications, we obtained overall accuracy of 71% and a moderate Kappa statistic of 0.56. Our results indicate that the LVB has undergone drastic changes in land use land cover, mainly driven by human activities that led to the conversion of forests, woodlands, grasslands, and wetlands to either farmlands or settlements. We conclude that information from this work is useful not only for basin-scale assessments and monitoring of land cover changes but also for targeting, prioritizing, and monitoring of small scale, community led efforts to restore degraded and fragmented areas in the basin. Such efforts could mitigate the loss of ecosystem services previously derived from large contiguous land covers which are no longer tenable to restore. We recommend adoption of a basin scale, operational, Earth observation-based, land use change monitoring framework. Such a framework can facilitate rapid and frequent assessments of gains and losses in specific land cover classes and thus focus strategic interventions in areas experiencing major losses, through mitigation and compensatory approaches.


Introduction
The Lake Victoria Basin (LVB) land mass covers approximately 184,200 square kilometers, traversing five East African countries (Kenya, Uganda, Tanzania, Rwanda, and Burundi) [1]. It is home to an estimated 35 million people, with the registered population growth within the 100-km buffer zone around Lake Victoria reported to be significantly higher than that of the rest of Africa. This is due to the wealth of natural resources, economic benefits, and ecosystem services that the basin offers [1]. The LVB contains the world's largest tropical lake [2] and supports diverse socioeconomic activities, among them agriculture, livestock keeping, and fishing activities, and thus is a significant resource for food security for the five countries. It is also important for energy production, domestic and irrigation water, shelter and transport, and recreation [3].
The demand for land and its associated resources and ecosystem services continues to grow all over in East Africa, which drives human related land use changes and subsequent land cover change over short and long-term scales. Land use land cover (LULC) change in East Africa is a widely studied subject [4][5][6][7], and a number of factors have been identified as key drivers of land use land cover change [1]. There are recent continental scale studies [5,8], national, basin, or sub-basin scale studies [9,10], as well as regional efforts [11]. Some work has focused on relatively small riverine catchments, with a focus on the local drivers of land use change [10,12,13]. Attention has also been leveled at the linkages between climate change and variability in land use change [4,11] while others have investigated the relationships between population growth, economic development, and environmental sustainability [6,14]. In almost all the papers, the authors agree on the key drivers of land cover change, mainly land use conversions to agriculture and urbanization (settlements) [1,13]. These researchers have provided invaluable information on the drivers of land use change over time and in some instances advocated for specific mitigation measures for reducing massive land use changes which impact negatively on ecosystem services and the climatic regime of the region. While the drivers are well recognized, correlations among land cover types at a basin scale are yet to be shown. In addition, for the LVB, gaps still abound in land use land cover changes in recent decades and the functional relationships between the area under natural vegetation and that under croplands and urbanized zones at a basin scale. In a basin shared by five countries, land use policies and business regimes influencing land cover changes vary. Furthermore, interpretations of land cover classes across the countries often differ, thus affecting the computation of harmonized estimates of land cover change. This can also adversely affect the quantification of ecosystem services offered by natural land cover, thus leading to underestimation of the value of these services. In the absence of regionally acceptable land cover classifications, the assessment of gains and losses of important ecosystems would not be feasible, and this hampers efforts in joint land use management and conservation at regional level. To address such concerns and create consistency, it is necessary to create a basin wide, standardized land cover map using an internationally recognized classification system over a period of approximately 30 years. Towards this end, our specific objectives were (i) to quantify land use land cover change from 1985 to 2014 using an internationally recognized classification system in the LVB and (ii) to investigate the spatial patterns and land cover changes per class and assess the relationships between areas under different land cover classes.

Study Area
The LVB (184,200 km 2 ) ranges in altitude from 1100 to over 4000 m (highlands in Rwanda, Elgon and Mau escarpment) above sea level and has a variety of terrains that are characterized by mountains to gentle plains in some areas ( Figure 1). The steep highlands are often affected by erosion and landslides during the rainy seasons [15,16]. The soil types are variable and are largely influenced by the volcanic activity of the Great Rift Valley. Soil type and fertility is a major determinant of cropping types and other land use types such as livestock keeping. The vegetation types across the basin include montane forests, savannahs, grasslands, wetlands, woodlands, and croplands (tea, cereals, coffee etc.) [1]. Vegetation cover is influenced by soil types, climatic patterns (precipitation and temperatures), and human activity over time [3]. Lake Victoria and its network of rivers and streams constitute the major hydrological features in the basin. These are important sources of water for various uses within the basin (domestic, agriculture, transport, etc.). The lake and rivers are also major conduits of pollutants and sediment transport and thus an important repository for human, agricultural, and industrial waste [1,17]. The basin is also significant for biodiversity and thus plays a critical role in conservation and promotion of tourism activities [3]. The LVB experiences a complex climatic system which is influenced by a number of inter-tropical convergences and the Indian Ocean Dipole [18,19], which creates mainly two rainy and dry seasons [18]. Over the years, land use and land cover has evolved within the basin due to a number of factors. For instance, major cities (Jinja, Kisumu, Mwanza) in each of the East African countries have grown over the years, accompanied by upcoming towns in the outskirts. Land use and land cover change in the basin is one of the key drivers of water quality degradation in Lake Victoria and its tributaries [20].
Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 17 major conduits of pollutants and sediment transport and thus an important repository for human, agricultural, and industrial waste [1,17]. The basin is also significant for biodiversity and thus plays a critical role in conservation and promotion of tourism activities [3]. The LVB experiences a complex climatic system which is influenced by a number of inter-tropical convergences and the Indian Ocean Dipole [18,19], which creates mainly two rainy and dry seasons [18]. Over the years, land use and land cover has evolved within the basin due to a number of factors. For instance, major cities (Jinja, Kisumu, Mwanza) in each of the East African countries have grown over the years, accompanied by upcoming towns in the outskirts. Land use and land cover change in the basin is one of the key drivers of water quality degradation in Lake Victoria and its tributaries [20].

Earth Observation Data Analysis
We analyzed satellite remote sensing data for several epochs in the LVB. A total of 118 Landsat 2, 3, 4, 5, 7, and 8 (MSS, TM, ETM+, and OLI sensors) scenes were downloaded from the USGS Global Visualization Viewer (https://glovis.usgs.gov/) for the years 1985, 1990, 1995, 2000, 2010, and 2014, with approximately 19-20 scenes covering the area of interest defined as the LVB. The specific scenes covered the following paths and rows on the world reference system (WRS): 169/059-63, 170/059-63, 171/061-63, 172/060-63, and 173/060-63. The Landsat data stream was deemed appropriate for mapping basin scale land use changes at 30 m resolution, with minimal bias from seasonality. Image selection was based on data availability within the dry season of the area under consideration, which presented the best chance of minimizing cloud contaminated images. In cases where images from a target year were not available, the images from the previous or next year were used (no more than 2-3 image substitutions for the target year). We adopted the land use classification approach described by the Food and Agriculture Organization's Land Cover Classification System (FAO-LCCS), which uses the dichotomous and modular phases in succession to define land cover classes [21]. Slight modifications were made to accommodate classes that are applicable to the study area

Earth Observation Data Analysis
We analyzed satellite remote sensing data for several epochs in the LVB. A total of 118 Landsat 2, 3, 4, 5, 7, and 8 (MSS, TM, ETM+, and OLI sensors) scenes were downloaded from the USGS Global Visualization Viewer (https://glovis.usgs.gov/) for the years 1985, 1990, 1995, 2000, 2010, and 2014, with approximately 19-20 scenes covering the area of interest defined as the LVB. The specific scenes covered the following paths and rows on the world reference system (WRS): 169/059-63, 170/059-63, 171/061-63, 172/060-63, and 173/060-63. The Landsat data stream was deemed appropriate for mapping basin scale land use changes at 30 m resolution, with minimal bias from seasonality. Image selection was based on data availability within the dry season of the area under consideration, which presented the best chance of minimizing cloud contaminated images. In cases where images from a target year were not available, the images from the previous or next year were used (no more than 2-3 image substitutions for the target year). We adopted the land use classification approach described by the Food and Agriculture Organization's Land Cover Classification System (FAO-LCCS), which uses the dichotomous and modular phases in succession to define land cover classes [21]. Slight modifications were made to accommodate classes that are applicable to the study area and reflect particular land use changes with respect to ongoing human activities. This approach was found necessary to define regionally consistent land uses in a basin whose geographical coverage traverses 5 countries, where national definitions of various land use categories may vary considerably. The land use land cover classes considered in this work were bare soils, urban areas, wetlands, waterbodies, small scale farmlands, large scale farmlands, large scale farmlands (plantation), open grasslands, closed grasslands, open woodlands, closed woodland, and indigenous forests. We adopted definitions of land use land cover classes described by the Intergovernmental Panel on Climate Change (IPCC) (https://www.ipcc-nggip.iges.or.jp/public/gpglulucf/gpglulucf_files/Chp2/Chp2_Land_Areas.pdf) and expanded the definitions of farmlands, grasslands, and woodlands to make the classifications more responsive to land use in the region and also communicate variations of these land use types in the Lake Victoria Basin. Large scale farmlands are mainly tea, coffee, rice, and sugarcane plantations, while small scale farmlands are characterized by small holder farms practicing subsistence agriculture for food crops. Grasslands and woodlands are an important type of land cover, especially for the areas south east of Lake Victoria, and the cover can vary depending on the level of disturbance. For consistency across the basin, maps were georeferenced to the WGS84 projection.
Image processing and classification was conducted in ENVI version 5.1 (Exelis Visual Information Solutions, Boulder, CO, USA), where geometric and radiometric correction, band stacking, visual interpretation of images, and creation of regions of interest (ROI) useful for classification was done. The analysts conducting image classification had a good understanding of the region, and their knowledge was useful in visual feature recognition, aided by high resolution imagery available in Google Earth. This work applied the maximum likelihood classification algorithm (supervised classification), which is a pixel-based classification that classifies images based on the homogeneity of image pixel spectral information. Classification was done in ENVI, where the ROI polygons picked by the analysts were used as training areas for image classification. The classified images were then taken through post classification processing and assignment of consistency codes (using a decision tree approach) for respective classes before mosaicking of individual scenes. A post classification smoothing was done by applying a 3 × 3 kernel filter to minimize the "salt and pepper" effect associated with pixel-based methods.

Reference Data Collection and Accuracy Assessment
A team of reference data collectors visited parts of the basin that were accessible in Kenya, Uganda, and Tanzania ( Figure 1). Prior to the field visits, the analysts discussed and built a consensus around the identification and documentation of land cover classes in the field, following principles of a response design [22]. The number of reference points collected per class was also determined by the accessibility of sites and cost of field visits [22]. The ArcMap software was used to prepare a shapefile of ROIs from which random points would be selected for validation and accuracy assessment. This shapefile was then loaded in ENVI and overlaid on the classified image, from which selected areas of the classified image were chosen as candidate sites for the creation of a representative sample of points to be used for accuracy assessment. A stratified random sample of selected classes was then processed and exported to a shapefile, from which they were loaded into a GPS unit, ready for use in the field. The data collectors used handheld GPS gadgets to locate and mark distinct land cover classes within the basin. Concurrent with the GPS points, the data collectors also took photographs of land cover, which were further used to help in interpretation of the classified images for the 2014 epoch. This was deemed necessary to minimize geolocation errors and interpreter uncertainty [22]. A total of 124 GPS points covering major land cover types in the basin were collected for accuracy assessment as shown in Figure 1. The overall accuracy and kappa statistics of agreement for the 2014 land cover map were computed. A change detection analysis was conducted between selected epochs using ArcGIS version 10, and land use cover change matrices were extracted to determine land use cover transitions in the basin. Our analysis workflow is illustrated in Figure 2. We derived the hectarage for each land cover class for each epoch (1985, 1990, 1995, 2000, 2010, and 2014) and subsequently computed the percent increase or decrease in area per land cover class for each epoch year, using 1985 as the base year. We also regressed the area under indigenous forests, open woodlands, wetlands, closed woodlands, largescale farmlands (plantation), largescale farmlands, and closed grasslands against either urban areas or small scale farmlands. Regression analysis enabled us to measure the statistical relationships between land cover patterns elucidated in our work. We hypothesized that human activity (e.g., urban growth, expansion of small scale agriculture, deforestation) in the LVB could be correlated with a decline in natural vegetation, which is a phenomenon witnessed globally [23].

Results
Our results show that there have been remarkable land cover land use changes in the Lake Victoria Basin, detectable from satellite remote sensing and in-situ observations over the last 30 years. The land cover maps for 1985,1990,2000,2010, and 2014 are shown in Figure 3. The overall accuracy for the 2014 epoch for which the field data were available was 71% (Table 1)  We derived the hectarage for each land cover class for each epoch (1985, 1990, 1995, 2000, 2010, and 2014) and subsequently computed the percent increase or decrease in area per land cover class for each epoch year, using 1985 as the base year. We also regressed the area under indigenous forests, open woodlands, wetlands, closed woodlands, largescale farmlands (plantation), largescale farmlands, and closed grasslands against either urban areas or small scale farmlands. Regression analysis enabled us to measure the statistical relationships between land cover patterns elucidated in our work. We hypothesized that human activity (e.g., urban growth, expansion of small scale agriculture, deforestation) in the LVB could be correlated with a decline in natural vegetation, which is a phenomenon witnessed globally [23].

Results
Our results show that there have been remarkable land cover land use changes in the Lake Victoria Basin, detectable from satellite remote sensing and in-situ observations over the last 30 years. The land cover maps for 1985,1990,2000,2010, and 2014 are shown in Figure 3. The overall accuracy for the 2014 epoch for which the field data were available was 71% (Table 1), with a kappa score of 0.56. The proportion of each land cover class as percentage of the total area in 1985 and 2014 is shown in Figure 4, indicating increases in the area under small scale farmlands. The general trend over the years indicates that, while the area under small scale farmlands and urban areas continues to grow, that under indigenous forests, open and closed woodlands, wetlands, and open grasslands is declining ( Figure 4). Figure 5 emphasizes snapshots of selected areas (whose locations are shown in Figure 1) across the land cover classification spectrum where change is remarkably visible.     Our area computations were consistent with published area estimates for the Lake Victoria water mass and the LVB land mass, approximately 6 million and 18 million hectares, respectively (Table 2). From 1985, land cover land use change in the LVB is noticeable, with increases in the hectarage under small scale farmlands, urban areas, and bare soil ( Table 2). Between 1985 and 2014, the urban areas grew by over 800%, the bare soil by over 140%, the small scale farmlands by over 90%, and the large scale farmlands by 55%. Furthermore, the open grasslands, open woodlands, closed woodlands, and   Our area computations were consistent with published area estimates for the Lake Victoria water mass and the LVB land mass, approximately 6 million and 18 million hectares, respectively (Table 2). From 1985, land cover land use change in the LVB is noticeable, with increases in the hectarage under small scale farmlands, urban areas, and bare soil ( Table 2). Between 1985 and 2014, the urban areas grew by over 800%, the bare soil by over 140%, the small scale farmlands by over 90%, and the large scale farmlands by 55%. Furthermore, the open grasslands, open woodlands, closed woodlands, and indigenous forests declined by 41%, 43%, 75%, and 44%, respectively ( Table 2). The mean area per class (shown as percentage) and the proportions associated with each year are shown in the Sankey diagram in Figure 6. We also noticed highly significant correlations between the decline in indigenous forests, open and closed woodlands, wetlands, and open grasslands and the growth in small scale farmlands. In addition, we observed significant correlations between indigenous forests, wetlands, and open grasslands and the growth in urban areas (Table 3). However, correlations between areas under open and closed woodlands and urban areas were not significant. Moreover, there is no significant correlation between the areas under large scale farmlands and closed grasslands and small scale farmlands ( Table 3). The land cover maps presented in this work are hosted at the RCMRD Geoportal (http://geoportal.rcmrd.org/layers/?limit=100&offset=0&title__icontains=Lake%20Victoria% 20Basin%20Land%20Cover) and are publicly available, which facilitates data sharing as well as other studies focusing on similar objectives.

Discussion
The goal of this work was to assess basin scale land use land cover changes over the last 30 years (1985-2014) and thus determine trends in land use change over time and subsequent deductions in changes to ecosystem services. We used Earth observation data to look at basin scale changes in land use land cover because remote sensing data provide synoptic coverage of global landscapes and also present systematic historical records of changes for consistent comparisons. Our findings indicate that land use land cover changes in the Lake Victoria Basin are mainly driven by small scale agricultural activities and urban development (settlements). The changes in land cover and associated trends in the LVB are similar to those reported in other studies [12,13], with minor discrepancies in spatial and temporal scales of previous studies, a shortcoming that we counterbalanced in our basin wide, 30-year study. While comparing our maps with the maps produced by the European Space Agency Climate Change Initiative (ESA-CCI) land cover project (http://maps.elie.ucl.ac.be/CCI/viewer/index.php) for the LVB in 2014, we also noted similarities, particularly in the classification of farmlands (croplands). However, the ESA-CCI maps were produced at 300-m resolution and present far more disaggregated land cover classes than our maps. Our results indicated overall accuracy of 71%, which is slightly below the 80-85% reasonable limit recommended for overall accuracy of land cover estimates derived from remotely sensed data [24], but similar to results obtained by previous studies in the region using Landsat and MODIS datasets [4,25]. A number of factors can potentially affect the overall accuracy of land cover classifications. Key among them are classifier selection and class distribution in the training dataset [26], as well as the accuracy of the ground reference dataset [27], sampling issues, errors in classification, and accuracy data sampling. Other factors also include heterogeneity of landscapes with complex land cover mosaics and misrepresentation and misinterpretation of the temporal sequence of remotely sensed imagery [28]. We suppose that the accuracy of our classifications could have been affected by the sample size of reference data and the heterogeneity of the land cover mosaic in the LVB, where fine land cover classes can be challenging to disaggregate. Indeed, other researchers within the basin have argued that classification error rates are lower for datasets with fewer land cover classes [25]. However, we feel that our results represent an important piece of knowledge, given that basin wide land cover classifications for the classes that we considered, over a period of 30 years, are rare. Other regional and/or global land cover maps do not provide the same temporal scale, standardized method across time, or spatial resolution to effectively assess LC dynamics in LVB [4,23].
The drivers of land use change are complex in the LVB, just like elsewhere in developing countries [7]. Land use land cover changes in the LVB have implications for water quality degradation in Lake Victoria, as well as causing a decline in the quality of ecosystem services within the catchment and the lake [17]. Our findings on land cover changes and trends can be explained largely by population growth and expansion of socioeconomic activities in the basin (e.g., agriculture, livestock keeping), demand for natural resources to support economic growth (timber, papyrus, and minerals), climate variability and natural disasters, and infrastructure development [7]. Governance, too, is also a key factor around these issues. The Lake Victoria Basin is one of the regions in East Africa with the highest population density [29], averaging approximately 500 persons per km 2 , and, in some parts of Kenya, densities of up to 1200 persons per km 2 have been reported [1]. The basin is home to over 30 million people [30] whose livelihoods depend mainly on agriculture and livestock production, beekeeping, mining, trade, and fishing [9]. Consequently, agricultural expansion and intensification, land degradation, and other anthropogenic activities are key agents of land use change in the basin. As the demand for food grows, land in the basin continues to be subdivided for small scale agriculture and settlement, leading to an increase in small farms and urban centers [9]. Expansion and intensification of agriculture (mechanization and fertilizer inputs) and increases in livestock numbers also render the land more prone to erosion [9]. With the influx of people in urban areas, there is increased demand for services and construction materials such as timber and building stones and fuel wood [31], which further drives conversion of forests and woodlands to either grasslands or farmlands [32] and quarries to bare ground [33]. Interestingly, we found significant inverse correlations between expansion of small scale farmlands and contraction of indigenous forests, grasslands, woodlands, and wetlands, consistent with previous findings [9] relating to a section of the LVB. This confirms that expansion of small scale agricultural activities is contributing to loss of indigenous vegetation. These dynamics are challenging to land use planners throughout Africa, where the urban population is expected to reach around 1.3 billion people by 2050, corresponding to 21% of the world's projected urban population [34].
The LVB has one of the highest rates of urban expansion in Africa [5]. While the land under urban areas in Africa is projected to increase by nearly 600% between 2000 and 2030 [5], our results show that, for the LVB, the area under urban land has increased [35] by over 800%, which may indicate that the expansion of urban areas within the basin is higher than the continental average. The growth in urban areas from small towns to cities comes with major needs for infrastructure to service not only the urban populace but also the urban-rural interactions. The needs include higher demands for food (irrigation, livestock, fisheries, and aquaculture projects), housing, energy (electricity and fuel wood), transport, healthcare, schools (primary to tertiary), security and governance, communication, flood control structures, and solid and liquid waste treatment plants, etc. All this infrastructure (or part of it) constitutes the modern-day infrastructural layouts in urban centers. The development of such infrastructure which is intended to support current and future populations is competing for the same amount of land resources. Our results indicate highly significant (p < 0.05) and strong inverse correlations between the area under urban centers and area under indigenous forests, wetlands, and open grasslands over the years. This implies that urban growth is contributing to the loss of forests, wetlands, and grasslands, either through increased demand for space, materials, or repositories for waste, which leads to land use change, fragmentation, or degradation. Furthermore, we found a highly significant (p < 0.05) and strong positive correlation between the area under urban centers and small scale farmlands. While this may appear an unusual finding, it can be explained by the need to substitute the land taken up by urban areas with new farmlands, hence conversions from other classes to farmlands. This could imply that both the expansion of urban centers and small scale farmlands is taking place at the expense of other land cover classes, and the scale of such conversions can certainly be guided by data-based government policy [36]. It is also true that the demand for land in these two classes is often driven by very similar or the same dynamics. From these findings, it is important for the five states sharing the LVB to encourage drafting, implementation, and frequent review of spatial land use plans across the basin. Spatial land use plans would help governments and stakeholders to have a consolidated and long-term plan for the smaller administrative units that constitute the basin in each of the countries. Where such plans exist and have been implemented with a reasonable degree of success, they could be used as blueprints for drafting spatial plans for other areas. At present, medium to high resolution Earth observation data would be an important complimentary dataset to other sets of information required to compile a comprehensive spatial land use plan. Ultimately, solving the problem of ad hoc land use changes would go a long way in aligning infrastructure development with development plans for the basin and minimize conversions of land use from one use to another without considerations for long-term sustainability.
The LVB is well endowed with natural capital such as forests, woodlands, grasslands, minerals, large rivers, wetlands, fisheries, and fertile farmlands whose demand continues to rise. Efforts to quantify and document the value of the natural capital within the basin have shown promising results [36]. These resources support livelihoods and thriving international trade [1] within the five states that share the LVB. The implications of land use land cover change on different sectors of the economy, such as agriculture and livestock production and mining [9,17], are now evident and can be discussed at length. The basin's economic growth and development is to a large extent dependent on the availability of natural resources, land being one of the key resources in this region. The fundamental question for the LVB (like other regions in Africa) is this: can the basin manage positive economic growth and development and still preserve the ecological services from its natural landscapes? Besides growth related pressures on natural resources, it is also important to consider external factors which affect the natural resource base and drive land cover change. Natural disasters such as wild fires and landslides can cause significant changes in land cover, and, without sufficient time for the ecosystem to recover, it can lead to permanent changes, which exposes the land to further degradation [37]. Forest and savannah fires are reported in the basin especially during prolonged droughts [11]. Furthermore, instability in governance in any one of the countries could trigger civil conflict and cause an influx of refugees or internally displaced persons who exert pressure on natural resources such as woodlands for fuelwood and construction of temporary shelters, fresh water, and land for cultivation [5]. In addition, foreign direct investments on large scale agriculture could also spur land use land cover change [5].
Climate variability and climate related natural disasters have also contributed to land use land cover change in the basin over the years [38]. Droughts in the basin affect food production, availability of water, and generation of hydroelectric power [39]. Severe droughts also affect vegetation and thus increase the percentage of bare ground exposed to agents of erosion. After extreme drought, some vegetation types take much longer to recover [38] and could even give way for colonization of the landscapes by invasive species. According to one study [32], the LVB is a significant hotspot where people suffer from exceptionally high exposure to climate change impacts, primarily due to the collocation of population increase on lands with decreasing rainfall. Besides droughts, parts of the basin also experience flooding, often with catastrophic consequences. Extreme flooding is known to inundate productive land, thus affecting the production cycles, as well as causing massive erosion, destruction of property, and loss of lives [40]. On the contrary, partial flooding could be beneficial for vegetation recovery in drylands [37]. Parts of the basin in Rwanda and Burundi are also highly susceptible to landslides during periods of extreme rainfall [41]. The land related issues surrounding drought and flooding cycles in the basin could partly explain the demand for land cover conversions from indigenous forests, woodlands, grasslands, and wetlands to small scale farming and settlement.
Mitigating the impacts of land use land cover change on ecosystem services of the Lake Victoria Basin has been ongoing work for many years, particularly under the umbrella of the East Africa Community, the Lake Victoria Basin Commission, and other national level agencies. The details of these efforts are beyond the scope of this paper and are adequately addressed elsewhere [3]. However, it is worth pointing out that our findings indicate that the major drivers of land use change in the LVB have a strong human dimension, particularly in connection to the demand for food production and settlement [42]. Any future efforts to address the impacts of land use change would need to intensify measures to address needs around these two factors.
Our work was subject to certain limitations that are worth highlighting. Firstly, we used 30-m Landsat data, which might have challenges in resolving fine changes in land cover change and can introduce errors. Secondly, optical data are subject to interference by clouds, and even though the percentage of cloud cover relative to the entire basin is small, it could still influence the computation of total area under a land cover class. We made efforts to select imagery with as little cloud cover as practically possible. However, we acknowledge that algorithms to generate cloud free mosaics of Landsat imagery are available (e.g., under Google Earth Engine platform). Using a near cloud free mosaic of Landsat data could change the results slightly, but probably not the trends. Thirdly, our results are at basin scale and the land cover class definitions could vary slightly when interpretations are made at national or sub-national levels. The overall accuracy (71%) was slightly below the recommended threshold of 85% for land cover classifications from remote sensing data [24]. Previous studies have come up with similar results for work within the basin and attributed those results to classification error rates being lower for datasets with fewer land cover classes [25]. In addition, the accuracy could improve if a bigger sample of reference data was obtained. In our work, this was constrained by project timelines, reference data collection logistics, and funding.

Conclusions
Satellite remote sensing data analyzed for the LVB over a period of 30 years show significant impact of human activities on land use land cover for the basin. Key among the drivers of land use land cover for the basin are urban growth and expansion of agriculture at the expense of forests, grasslands, and woodlands. This predisposes the land to severe erosion and pollution of aquatic environments, which has been witnessed in the basin over the years.
The utility and value of information from our work can be two-fold: (i) for monitoring land use land cover changes over time and therefore providing empirical estimates of the changes and drivers of land cover change for decision making; (ii) design of a targeted and prioritized ecosystem restoration system whose investments, achievements, and gains can be measured using Earth observation data and other socioeconomic metrics. In addition to the elaborate regional strategic action plans addressing landscape restoration, conservation, waste management, soil conservation, and improved farming practices, we recommend adoption of a basin scale, operational, Earth observation-based, land cover and land use change monitoring framework. Such a framework can facilitate rapid and frequent assessments of gains and losses in specific land cover classes and thus focus strategic interventions in areas experiencing major losses, through mitigation and compensatory approaches. Cognizant of the relentless fragmentation and modification of natural land cover (forests, grasslands, wetlands) through conversions to other land uses, and the realization that the LVB land use plans fall under various jurisdictions, we are of the opinion that successful basin scale land use management efforts will benefit immensely from coordinated national, sub-national, and community level inputs. Here, we propose that Earth observation (EO) data analysis, outputs, and tools should move beyond quantifying the magnitude of the problem, to facilitating design and refinement of community level targeting, prioritization, monitoring, and restoration efforts. In areas where land cover has changed from one class to another, but the land use has not changed officially, assessing the societal cost of restoration is important. In other areas, it may be argued that some of the land use changes witnessed in the basin are radical (e.g., 100-hectare forest converted to farmlands which are eventually turned into fragmented settlements) and therefore cannot be reversed through restoration. While this is true, through spatial analysis and community level consultations, authorities could explore compensatory approaches to identify and prioritize the restoration and connectivity of fragmented units (forests, grasslands, woodlands, etc.) whose total area could equal that lost from the conversion of large units of homogenous land cover classes. The fragmented sets of such units could be designed to offer complimentary ecosystem services lost from destruction of larger homogenous and contiguous natural land cover, where economic, environmental, and social considerations are optimized. Such an approach can be practically supported by an operational EO-based land use land cover assessment tool which regularly provides estimates of losses and gains in key land cover classes for the basin and thus influences short to long-term decision making. Dissemination of such information on web-based platforms which can be accessed via mobile computing devices can greatly improve accessibility of the data and information and thus foster evidence-based decision making.