How Climate Change and Land Use/Land Cover Change Affect Domestic Water Vulnerability in Yangambi Watersheds (D. R.

: In the tropics, the domestic water supply depends principally on ecosystem services, including the regulation and puriﬁcation of water by humid, dense tropical forests. The Yangambi Biosphere Reserve (YBR) landscape is situated within such forests in the Democratic Republic of Congo (DRC). Surprisingly, given its proximity to the Congo River, the YBR is confronted with water issues. As part of its ecosystem function, the landscape is expected to reduce deterioration of water quality. However, environmental consequences are increasing due to conversion of its dense forest into other types of land use/land cover (LULC) in response to human activities. It is therefore important to check how the physicochemical quality parameters of water resources are inﬂuenced by landscape parameters—and to know if the population can adapt to this water vulnerability. To do this, we analyzed the watershed typology (including morphometric and LULC characteristics) and the physical and chemical parameters of water within the principal watershed’s rivers. We also analyzed data from surveys and the Yangambi meteorological station. We found that some landscape indices related to LULC signiﬁcantly inﬂuence water quality deterioration in Yangambi. On average, each person in the Yangambi landscape uses 29–43 liters of water per day. Unfortunately, this falls short of World Health Organization standards regarding some parameters. The best ﬁtted simple linear regression model explains the variation in pH as a function of edge density of perturbed forest, edge density of crop land and patch density of dense forest up to 94%, 92% and 90%, respectively. While many researchers have identiﬁed the consequences of climate change and human activities on these water resources, the population is not well-equipped to deal with them. These results suggest that water management policies should consider the speciﬁcities of the Yangambi landscape in order to develop better mitigation strategies for a rational management of water resources in the YBR in the context of climate change.


Introduction
Many scientists are increasingly noting the influence of land use patterns on the degradation or improvement of water quality [1,2]. Among other types of land use/land cover (LULC), dense tropical forest has been found to help improve the physicochemical quality of water at the watershed level. In addition to other services, these forests also offer vital ecosystem services related to the regulation and purification of water resources. They also improve the supply and quality of water used for the different needs of forest households [3]. In 2006, for example, one study estimated that at least 77 million people depended directly or indirectly on these ecosystems for water supply for their different domestic demands in the Congo Basin [4].
Given that dense tropical rainforests cover three-quarters of the Democratic Republic of Congo (DRC) [5], the country should not experience water problems [6]. Indeed, DRC accounts for around 11% of total tropical rainforests worldwide and 60% in Africa [7]. This constitutes a major asset in terms of the regulation and renewal of its water resources. DRC alone has the most important hydrological resource in Africa-and in the world after the countries of the Amazon basin [6]. Its surface water represents about 35% and 52% of water reserves of the world and Africa, respectively [8]. Due to DRC's low population density, its per capita water availability potential is very high: 15,000 to 70,000 m 3 per person per year [9].
Unfortunately, as in most sub-Saharan African countries, the rate of access to safe water and sanitation remains paradoxically low [4]. Despite DRC's large hydrological reserves, only 29% of its rural population and 79% of its urban population have access to quality water to meet their various basic needs, particularly those related to domestic use [9].
Mismanagement is often cited to explain water issues in sub-Saharan countries in general-and in DRC in particular [10]. In the past, water management in DRC focused mainly on the demands for urban supplies, irrigation and hydropower. This approach did not consider the potential impacts on the environment. Indeed, a good proportion of the population in rural areas lacked access to drinking water [4].
In Yangambi, an urban-rural area of Tshopo province in DRC, the National Rural Hydraulic Service (NRHS) is struggling to fulfill its mission because of a lack of funds and insufficient capacity to monitor the water resource. The few boreholes for standpipes in certain settlements (Bangala, Lumumba and Ekutsu), drilled in collaboration with partners, are no longer operational. Hence, access to quality water from the water points remains problematic. The population obtains water for its various domestic needs from springs and rivers, which are mostly unprotected.
A survey of 29 water suppliers in the United States by the Trust for Public Lands and the Spring Protection Committee showed that treatment costs are inversely related to the proportion of the watershed protected by forests, wetlands and other open space [11]. In many regions, loss of forest land leads directly to water quality degradation. Therefore, land conservation and pollution prevention have proven to be cost-effective strategies [11]. Unfortunately, research in DRC rarely aims to quantify the site-specific benefits and costs of land conservation. This is the case when it comes to YBR's forest massif. It should function as a buffer zone to limit deterioration of water quality for domestic use, at a lower cost, as is the case in several places around the world [4]. However, it is getting more and more degraded. This massif is under intense pressure as a result of climate change [12] and the various livelihood activities linked to the demographic explosion in the zone [13]. Different land uses, mainly small-scale slash-and-burn agriculture, were reportedly responsible for an estimated 84% of forest disturbance in the Congo Basin from 2000 to 2014, fragmenting it into other new LULC types [14,15]. Most alarmingly, predictive models of forest cover in this area show the situation is likely to worsen [14].
This deterioration of primary forest cover into other LULC types is thus likely a danger to the sustainability of watersheds in terms of water quality regulation [16]. Integrated watershed management practices and their conservation strategies-including The study site is in the watersheds (irrigated by the Lusambila, Isalowe, Bonde, Lubilaie and Loweo rivers) of Yangambi in Tshopo province of DRC. It extends over an area of about 430.2 km 2 . The geographical position is between 0 • 45 58.1 north latitude and 24 • 4 6.13 east longitude, approximately 100 km west of Kisangani. Five main rivers irrigate this study area, of which three (Lusambila, Isalowe and Bonde) are direct tributaries to the Congo River, while the two others (Lubilaie and Loweo) initially flow into the Lobilo stream before flowing into the Congo River ( Figure 1). These watersheds are part of the larger Congo River watershed in the equatorial region, consisting mainly of dense rainforests [22].
watershed management practices and their conservation strategies-including social, economic and ecological components-are important ways to improve water quality production and reduce environmental degradation in the territory [17]. Thus, it is getting more difficult for countries to promote development strategies in protected areas such as the YBR without demonstrating their importance for local communities and society in general, especially the links between protected areas and water resources. For example, in many countries, the linkage between domestic water use and drinking water must be integrated into the management of these areas [18].
Many reports by development agencies have noted the degradation of water quality in the region. In response, a careful assessment of the state of domestic water resources at the YBR(Yangambi Biosphere Reserve) is needed. This was proposed in a report by the Center for International Forestry Research (CIFOR) in 2018 [19]. Such an assessment must consider how different types of LULC at the watershed level can influence the quality of the resource. An accurate assessment of the water resource, which considers its various components, is vital to ensure successful water management [20][21].
Hence, the objectives of this study are (i) to typologize the watersheds in the Yangambi area in connection with the water physicochemical quality of the rivers in Yangambi; (ii) to verify whether the water physicochemical quality of the rivers in Yangambi meets World Health Organization (WHO) standards; and (iii) to analyze and assess the link between water and climate change adaptation (including vulnerability to climate change).

Study Site
The study site is in the watersheds (irrigated by the Lusambila, Isalowe, Bonde, Lubilaie and Loweo rivers) of Yangambi in Tshopo province of DRC. It extends over an area of about 430.2 km 2 . The geographical position is between 0°45'58.1'' north latitude and 24°4'6.13'' east longitude, approximately 100 km west of Kisangani. Five main rivers irrigate this study area, of which three (Lusambila, Isalowe and Bonde) are direct tributaries to the Congo River, while the two others (Lubilaie and Loweo) initially flow into the Lobilo stream before flowing into the Congo River ( Figure 1). These watersheds are part of the larger Congo River watershed in the equatorial region, consisting mainly of dense rainforests [22].   To delineate the five main watersheds in the study area, we used a shuttle radar topography mission (SRTM) file to obtain a digital terrain model. From there, we generated flow directions, flow accumulation raster and talweg segments. With the coordinates of the outlet, these tools helped us produce the raster of the watershed limits. This was then polygonized, to transform it into vectors to bring out the watershed boundaries.
(b) Land use at watershed level We also produced landscape characteristics related to LULC in the watersheds from QGIS version 2.18.18 software. To detect differences in land cover within the watersheds, we downloaded two Sentinel_2B images for the complete study area from http://www. usgs.gov/ (accessed on 7 June and 1 September 2018) (United States Geological Survey). The processing of these two images consisted primarily in their correction. We assigned any clouds perceptible on the images-which could bias the quality of the classificationto the non-data class, i.e., without any appropriate theme. We then created virtual files of these two images. On the latter image, we carried out the random forest supervised classification to generate the different LULC in the Yangambi watersheds. Six training zones were created for each thematic class. The main land uses in Yangambi guided the choice of different training zones. These training zones and the images to be classified were imported into Monteverdi/OTB-6.2.0-Win64 software to generate the classification model and the confusion matrix. Based on this classification model, the Sentinel_2B image was classified into six different land use classes related directly to the different types of LULC in Yangambi. The confusion matrix allowed verification only of the accuracy of the classification performed by the computer. It was used to ensure no overlap between thematic classes. The latter had provided a precision rate of 96%, which was considered acceptable. Before carrying out the analyses, we checked the accuracy of the classification. To that end, we used 50 samples points produced randomly over the entire set of each polygon vector of the different types of LULC in the virtual classified image. We then imported these points into Google Earth to verify whether the thematic classes generated by the supervised random forest classification truly corresponded to the information that characterizes the study environment.
Finally, the LecoS 2.0.7 extension in QGIS version 2.18.18 allowed us to automatically produce landscape composition and configuration indices from the classified image. The following landscape indices were produced: Composition indices: dense forest (DF), perturbed forest (PF), crop land (CL), grass land (GL), and bare soil and residential (BSR). II. Configuration indices: edge density (ED), patch density (PD), great patch area (GPA), mean patch area (MPA) and number of patches for the different composition indices (Table 1).

Indices Descriptions
Edge density Total length of all edge segments per unit of area for the thematic class under consideration. Patch density Number of patches per unit of area, for a given thematic class Great patch area The area of the largest patch for a given thematic class Mean patch area Average area of patches for a given thematic class Number of patches Number of patches for a given theme class (c) Analysis of physicochemical quality of water Analysis focused mainly on the measurement of some physicochemical quality parameters, such as temperature, pH, conductivity, turbidity and dissolved oxygen. We recognized these parameters as basic or synthetic. However, parameter selection was constrained by lack of sufficient funds and the unavailability of materials to carry out more detailed analyses.
We carried out on-site analyses of pH, conductivity, dissolved oxygen and turbidity at the site of sample collection following the standard protocols and methods of the American Public Health Organization [23]. The temperature, conductivity, pH, dissolved oxygen and turbidity were measured using multiparameter HACH ® Test Kits. We calibrated the instruments daily using standard solutions prior to taking samples. Each sample was poured in the sample holder and kept inside for a few minutes. After achieving reading stability, we recorded the value. We collected water samples at the springs, at the cisterns' level 50 cm below the surface and at the rivers' outlets, preferably in the morning (to test the samples under homogeneous conditions), over 17 days, from 20 June 2019 to 6 July 2019.
(d) Study of the water resource's vulnerability in relation to climate change We used two different approaches to analyze the climate change vulnerability of water resources in Yangambi: a qualitative approach (based on survey results) and a biophysical approach (based on the meteorological data analysis from the Yangambi station).
The qualitative approach consisted of a survey on the state of water resources (including access, availability, capacity, use and environment) in Yangambi. This was based on the Water Poverty Index (WPI) as proposed by Sullivan and colleagues in 2002 [20]. A sample size of 200 individuals was used for this purpose, based in the different settlements (Yakako, Yanguma, Isalowe, Lusambila, Likango, Obiloto and Lumumba) in the watersheds. The biophysical approach consisted of analyzing meteorological parameters such as wind speed, mean annual temperature and mean annual rainfall.

Data Processing and Statistical Analysis
The fieldwork data were encoded in an Excel spreadsheet. The processing and statistical analysis of the data were done using the same Excel spreadsheet, R software (version R x 64 3.5.1) and QGIS software (version 2.18.18, extension LecoS 2.0.7).
We generated descriptive statistics, including mean, mode, standard deviation and coefficient of variation from the survey data. We performed the Shapiro and Kolmogorov tests to verify the normality and distribution of the data. We then compared the averages of the physicochemical parameters of the water points in the different watersheds based on an analysis of variance (ANOVA). To understand the influence of LULC parameters on water physicochemical quality, we used the Pearson correlation test and simple linear regression (SLR) models. Data with non-normal distribution were log transformed.

Morphometric Characteristics
The drainage density for all the watersheds varies between 0.5 and 0.6 km/m 2 , with a Gravellius compacity index (KG) value ranging from 1.8 to 2.5 (see Table 2), conferring an elongated shape for most of these watersheds. Only the Isalowe watershed tends toward a Land 2021, 10, 165 6 of 21 circular shape, with a KG value of 1.8. The watershed mean slope varies between 7.3-9.2% (Table 2), with lowest and highest values at Loweo and Isalowe watersheds, respectively. All five watersheds are geomorphologically mature as is shown by the general aspect of the hypsometric curves ( Figure 2)-i.e., with soft plains around their main stream channels, where the variations in altitude are small despite their large surface area. This could also indicate that at basin scale, the erodibility is negligible. However, the Lusambila and Bonde watersheds are steeper than those of Isalowe, Loweo and Lubilaie. They are less mature, and therefore have a higher probability of erosion and impact on river water quality compared to others ( Figure 2).

Morphometric Characteristics
The drainage density for all the watersheds varies between 0.5 and 0.6 km/m 2 , with a Gravellius compacity index (KG) value ranging from 1.8 to 2.5 (see Table 2), conferring an elongated shape for most of these watersheds. Only the Isalowe watershed tends toward a circular shape, with a KG value of 1.8. The watershed mean slope varies between 7.3-9.2% (Table 2), with lowest and highest values at Loweo and Isalowe watersheds, respectively.
All five watersheds are geomorphologically mature as is shown by the general aspect of the hypsometric curves ( Figure 2)-i.e. with soft plains around their main stream channels, where the variations in altitude are small despite their large surface area. This could also indicate that at basin scale, the erodibility is negligible. However, the Lusambila and Bonde watersheds are steeper than those of Isalowe, Loweo and Lubilaie. They are less mature, and therefore have a higher probability of erosion and impact on river water quality compared to others ( Figure 2).

LULC Patterns Related to Anthropization in the Different Watersheds in Yangambi
The classification of the Sentinel 2B satellite image into different LULC types identified five main thematic classes in the five watersheds selected for this study ( Figure 3). DF was the dominant thematic class for all the watersheds, while BSR represented the smallest proportion. Lusambila and Lubilaie had the largest and smallest proportion of BSR areas, respectively ( Figure 3).

Water Quality in Yangambi for the Main Domestic Water Supplies in the Different Watersheds
Most of the analyzed physical and chemical parameters of water from springs, rivers and cisterns in Yangambi meet WHO standards for drinking water. The exceptions are the acidic pH and turbidity values in surface water and unimproved water sampled from all the watersheds (Table 3). More acidic pH and small turbidity values were found in the Lubilaie and Loweo, which are the most dominated by DF cover. This could be explained by two factors; first, humic acid is released from the forest and drained to the rivers and vadose zones (discharging into the springs) through runoff; second, hypodermic water flux from the forested areas dominate in most of these watersheds. There were significant differences in pH for different water supplies throughout the study area (Table 4). Temperature and conductivity show significant differences (for all water points considered together) across the five watersheds (p-value < 0.05) ( Table 4). At the spring level, there were no significant differences for any of the parameters in the different watersheds. However, rivers showed significant differences in conductivity and turbidity ( Table 4). Values that did not meet WHO standards for drinking water for turbidity were found in Lusambila and Isalowe (Table 3), which are watersheds with a smaller proportion of DF ( Figure 4).

LULC Patterns Related to Anthropization in the Different Watersheds in Yangambi
The classification of the Sentinel 2B satellite image into different LULC types identified five main thematic classes in the five watersheds selected for this study ( Figure 3). DF was the dominant thematic class for all the watersheds, while BSR represented the smallest proportion. Lusambila and Lubilaie had the largest and smallest proportion of BSR areas, respectively ( Figure 3).

Water Quality in Yangambi for the Main Domestic Water Supplies in the Different Watersheds
Most of the analyzed physical and chemical parameters of water from springs, rivers and cisterns in Yangambi meet WHO standards for drinking water. The exceptions are the acidic pH and turbidity values in surface water and unimproved water sampled from all the watersheds (Table 3). More acidic pH and small turbidity values were found in the Lubilaie and Loweo, which are the most dominated by DF cover. This could be explained by two factors; first, humic acid is released from the forest and drained to the rivers and vadose zones (discharging into the springs) through runoff; second, hypodermic water flux from the forested areas dominate in most of these watersheds. There were significant differences in pH for different water supplies throughout the study area (Table 4). Temperature and conductivity show significant differences (for all water points considered together) across the five watersheds (p-value < 0.05) ( Table 4). At the spring level, there were no significant differences for any of the parameters in the different watersheds. However, rivers showed significant differences in conductivity and turbidity ( Table 4). Values that did not meet WHO standards for drinking water for turbidity were found in Lusambila and Isalowe (Table 3), which are watersheds with a smaller proportion of DF ( Figure 4).  Significant p-values are followed by an asterisk * for < 0.05 and *** for < 0.001.

Relationship Between the Physicochemical Quality of Stream Water Based on Indices of Landscape Configuration and Composition in Watersheds
The linearity analysis showed no significant linear relationships between the different water physicochemical parameters, except for turbidity and temperature, which showed a strong positive correlation (r = 92%). However, there were positive, nonsignificant relationships between almost all the variables of physicochemical water parameters.    Significant p-values are followed by an asterisk * for < 0.05 and *** for < 0.001.

Relationship between the Physicochemical Quality of Stream Water Based on Indices of Landscape Configuration and Composition in Watersheds
The linearity analysis showed no significant linear relationships between the different water physicochemical parameters, except for turbidity and temperature, which showed a strong positive correlation (r = 92%). However, there were positive, nonsignificant relationships between almost all the variables of physicochemical water parameters. The exceptions were for temperature and dissolved oxygen, and also for conductivity and pH, where the linearity correlation was negative (Table 5). The linear correlation test results between the different landscape indices (configuration and composition) and the different river water quality parameters revealed a significant correlation (p-value < 0.05), especially for turbidity, pH and temperature. However, almost all landscape indices were not significantly dependent (p-value > 0.05) on conductivity and dissolved oxygen, except for MPA_PF and MPA_CL, respectively. Only turbidity was negatively correlated with the proportion of DF and MPA of perturbed forests (Table 6.). The best SLR model, which explained the effect of each landscape parameter on different water physicochemical parameters, is that of pH as a function of ED-DF. Using this model, ED-DF explained up to 97% of the observed variability in pH (AIC, Akaike information criterion = −3.2; p-value > 0.05) ( Table 7).    A simple forecast of precipitation and temperature for 30 years in the future (2050) was performed using the available data ( Figures 5 and 6). The results are statistically acceptable with a symmetric mean absolute percentage errors (SMAPE) of 0.13 and 0.01 for precipitation and temperature, respectively. The predictions for the year 2050 showed a decrease in precipitation of 3% and 2% as compared to the mean precipitation of the periods 1931-2017 and 1970-2017, respectively. As for the temperature, an increase of 4% Celsius degrees is predicted for the year 2050, as compared to the initial mean temperature of 25.0 °C for the period 1970-2017. These values seem minimal as numbers, but their influence on the hydrology, water resource availability and supply within the Yangambi watershed may be huge.   A simple forecast of precipitation and temperature for 30 years in the future (2050) was performed using the available data ( Figures 5 and 6). The results are statistically acceptable with a symmetric mean absolute percentage errors (SMAPE) of 0.13 and 0.01 for precipitation and temperature, respectively. The predictions for the year 2050 showed a decrease in precipitation of 3% and 2% as compared to the mean precipitation of the periods 1931-2017 and 1970-2017, respectively. As for the temperature, an increase of 4% Celsius degrees is predicted for the year 2050, as compared to the initial mean temperature of 25.0 • C for the period 1970-2017. These values seem minimal as numbers, but their influence on the hydrology, water resource availability and supply within the Yangambi watershed may be huge.

Availability, Accessibility, Use and Management Capacity of Water Resources in Yangambi
Yangambi water supply points are located mainly around the urban-rural areas and villages. Nevertheless, several springs are located inside the forest, far from the residential areas. They were not considered in this study, as the Yangambi population barely uses them for water supply. There are only three wells, dug by the NRHS in Bangala, Likango and Ekutsu, but they are not operational. Cisterns are present exclusively in Isalowe, Lusambila and Lumumba (Isalowe, Lusambila and Loweo watersheds). These cisterns and wells are less than 1 km from settlements ( Figure 8).

Availability, Accessibility, Use and Management Capacity of Water Resources in Yangambi
Yangambi water supply points are located mainly around the urban-rural areas and villages. Nevertheless, several springs are located inside the forest, far from the residential areas. They were not considered in this study, as the Yangambi population barely uses them for water supply. There are only three wells, dug by the NRHS in Bangala, Likango and Ekutsu, but they are not operational. Cisterns are present exclusively in Isalowe, Lusambila and Lumumba (Isalowe, Lusambila and Loweo watersheds). These cisterns and wells are less than 1 km from settlements ( Figure 8). The survey results showed that the population of Yangambi in the different watersheds gets water for domestic use from: (i) springs, (ii) rivers, and (iii) community or private cisterns. Almost all respondents in the different watersheds reported using spring water, while only respondents in Isalowe and Lusambila used cisterns as springs. Water  The survey results showed that the population of Yangambi in the different watersheds gets water for domestic use from: (i) springs, (ii) rivers, and (iii) community or private cisterns. Almost all respondents in the different watersheds reported using spring water, while only respondents in Isalowe and Lusambila used cisterns as springs. Water from the Congo River was rarely used (Figure 9). The survey results showed that the population of Yangambi in the different watersheds gets water for domestic use from: (i) springs, (ii) rivers, and (iii) community or private cisterns. Almost all respondents in the different watersheds reported using spring water, while only respondents in Isalowe and Lusambila used cisterns as springs. Water from the Congo River was rarely used (Figure 9). In the different watersheds under study, on average, more than half of respondents said they lived at least 1000 m (1 km) from springs. In Lubilaie, only 25% said they lived more than 1 km from springs. In terms of distance from rivers, 63% of respondents in Lusambila said they lived more than 1 km from rivers used as water supplies. In Bonde and Lubilaie, respectively, more than 84% and 100% of respondents said they lived less than 1 km from the rivers ( Table 8). The average distance between households and water cisterns was less than 5 m (personal observation from the field). In the different watersheds under study, on average, more than half of respondents said they lived at least 1000 m (1 km) from springs. In Lubilaie, only 25% said they lived more than 1 km from springs. In terms of distance from rivers, 63% of respondents in Lusambila said they lived more than 1 km from rivers used as water supplies. In Bonde and Lubilaie, respectively, more than 84% and 100% of respondents said they lived less than 1 km from the rivers ( Table 8). The average distance between households and water cisterns was less than 5 m (personal observation from the field). The daily water consumption per person in all the watersheds studied was more than 20 liters (L) per person. The highest water consumption per person was in Isalowe (49 L per person per day) and the lowest (29 L per person per day) in Lubilaie (Table 8).
Water is mainly used for domestic purposes in Yangambi, mostly for (i) toilets, (ii) bathing, (iii) dishwashing, (iv) laundry, (v) cooking and (vi) drinking. Water from streams, cisterns and rainfall is used less for drinking and cooking, and more for toilets, bathing, dishwashing and laundry. On the other hand, water from springs and streams is used for everything, with a preference for spring water for drinking. In some settlements and villages at larger distances from springs, such as Obiloto and Yanguma (in the Bonde watershed), people sometimes use water from rivers or rainwater for drinking. Most households do not pretreat their water before using it for domestic purposes. Only a few households, especially in Isalowe and Lusambila, use clean cloths as filters, or purifying products such as chlorine, or boil drinking water before its use (Table 9). In the different settlements and villages, most people are poor and live in mud houses. In Isalowe and Lusambila, many houses are built with hard brick and have metal roofs. In Lubilaie, Loweo and Bonde, on the other hand, mud houses predominate (Table 9). Some shelters have mud walls with roofs mostly made from flowering plants (of the Marantaceae family) or leaves of oil palm trees.
The percentage of illiterate household heads was moderately low in all the watersheds, where the education level is, on average, limited to secondary school. The exceptions were in Loweo (where most household heads did not go beyond primary school), and in Lusambila and Isalowe (where several household heads attained university educations) ( Table 9). Additionally, most respondents said they had never received any training related to water. A few respondents in Isalowe, Lusambila and Loweo said they had attended some form of water-related training in the past. The survey results showed there were no organizations (governmental or non-governmental) actively involved in water management in Yangambi city, nor were there any well-known water point management committees during the study period.
The five watersheds have a Water Poverty Index (WPI) score varying between 50% and 66% ( Figure 10). The Isalowe and Lusambila watersheds have the highest scores in terms of water use and accessibility. The surveys revealed that water resource management capacity in Yangambi is generally low. The exception is the Isalowe watershed, which exceeds the 50% score ( Figure 10). The five watersheds have a Water Poverty Index (WPI) score varying between 50% and 66% ( Figure 10). The Isalowe and Lusambila watersheds have the highest scores in terms of water use and accessibility. The surveys revealed that water resource management capacity in Yangambi is generally low. The exception is the Isalowe watershed, which exceeds the 50% score ( Figure 10).

LULC Dynamics in the Yangambi Watersheds
Most of the five watersheds have DF as the highest proportion of LULC. The proportions of DF are much lower in the Isalowe and Lusambila watersheds, where they have a significant proportion of BSR, compared to Lobilo and Loweo watersheds. Analysis of the linear correlation between different landscape parameters in the five watersheds in Yangambi showed a negative relationship between the proportion of DF and most of the landscape parameters related to bare soil and residential, CL and GL. This means that as the PD of CL, ED of bare soil and residential and MPA of GL increase, the proportion of

LULC Dynamics in the Yangambi Watersheds
Most of the five watersheds have DF as the highest proportion of LULC. The proportions of DF are much lower in the Isalowe and Lusambila watersheds, where they have a significant proportion of BSR, compared to Lobilo and Loweo watersheds. Analysis of the linear correlation between different landscape parameters in the five watersheds in Yangambi showed a negative relationship between the proportion of DF and most of the landscape parameters related to bare soil and residential, CL and GL. This means that as the PD of CL, ED of bare soil and residential and MPA of GL increase, the proportion of DF patches tends to decrease in the watersheds. The CL, BSR, GL and PF landscape parameters also change positively with those showing an increase in forest landscape fragmentation (ED and PD), but also a decrease in its proportion in Yangambi.
Brandt and colleagues confirmed these relationships in 2016 [25]. They showed that the phenomenon of anthropization in the Congo Basin forests was generally accompanied by a conversion of the forest into other types of LULC. This is because the Congo Basin population depends mainly on slash-and-burn agriculture for subsistence. Hence, areas colonized by people are often characterized by intensive clearing for agriculture purposes [15]. These cropping areas then evolve either into fallow or other types of LULC dominated by herbaceous species, especially grasses (when the regeneration capacity of forest species is severely reduced by burning). This can also be explained by the location of the urban-rural area of Yangambi-largely within the YBR. The population only carries out its activities on small areas granted by the Institut National d'Etudes et de Recherches Agronomiques (INERA), and areas where there are repeated burnings on the same land. In 2016, Mikwa and colleagues [14] reported a decrease in the proportion of forest cover in Yangambi due principally to demographic pressure. The decrease could also be accentuated by management rules for preservation of the YBR, which are unclear and poorly enforced.

Water Quality in Yangambi According to WHO Standards
Apart from turbidity and pH, the physicochemical parameter analysis of the different Yangambi water points met WHO standards for drinking water [24]. Similar results in compliance with WHO standards were found in Haiti by Orelien in 2017 [26]. In the case of Haiti, of 29 physicochemical parameters analyzed, 28 complied with WHO and European Union standards. Nevertheless, possible sources of contamination influenced the quality of some water points in the YBR. These included agricultural practices, the discharge of livestock excrement into open areas, lack of demarcation of the water supply points for human consumption and animal watering, and deterioration of cisterns and roofs of houses.
Also, all the interviewees confirmed the lack of well-organized facilities for household waste management.
Comparison of the averages of the physicochemical parameters of rivers did not show any significant difference between watersheds (p-value < 0.05), except for conductivity and turbidity (p-value <0.05). The significant differences observed between the averages for turbidity and conductivity between watersheds can be explained by LULC types in Yangambi. Mwayi and Naito in 2019 found that forest massifs in Malawi could reduce sedimentation caused by watershed runoff to supply major water points [27]. This suggests the presence of the biosphere reserve in Yangambi may be the basis for the few results we found, in line with WHO standards, in the watersheds predominated by forest LULC.

Physicochemical Parameters of River Waters' Nexus Landscape Parameters in the Yangambi
The correlation analysis between different physicochemical parameters of rivers showed a positive relationship. The only exceptions were pH and conductivity-which were negatively correlated with each other-and temperature and dissolved oxygen (Table 5). This collinearity between the physicochemical parameter variables of stream water was confirmed by Xiao and Ji in 2007 [28]. Thus, only a few variables taken as indicators were needed to characterize the impact of LULC types on stream water quality. Pratt and Chang observed the negative correlation between temperature and dissolved oxygen in 2012 [29]. These authors identified biodegradable discharges in rivers as the possible cause of observed temperature rise. In fact, water quality is affected both by the biodegradable nature of these discharges and the degree of pollution by organic matter. Since microorganisms thrive on organic matter, their metabolism is largely responsible for the increase in water temperature and the decrease in dissolved oxygen in the water. Thus, the high temperatures of river waters in the highly anthropized watersheds (Lusambila and Isalowe) compared to those in the less anthropized watersheds (Bonde, Loweo and Lubilaie) could be related to the high rate of deposition of organic matter in the Lusambila and Isalowe rivers. Other factors for the rise in temperature could include atmospheric temperature at the time of water sample collection. Fortunately, our study collected samples in the morning hours during the rainy season (June-July).
The positive correlation between turbidity and pH was not significant, but could still be linked to the ferralitic soil type predominating in Yangambi, which usually undergoes deep leaching [30]. As hydrolysis is strong in equatorial environments, soluble elements such as silica and other cations are leached and transported to rivers by rainwater. The dissolution of atmospheric CO 2 and silica thus leads to the formation of bicarbonate dissolved in stream water. Bicarbonate, due to its amphoteric nature, plays a buffer role in maintaining the pH of the water around a more or less neutral value [31]. This would explain why stream water with a less acidic pH and high turbidity is found when the MPA landscape indices of bare soil increase at the watershed level (Lusambila and Isaloe) ( Figure 4).
Thus, as the proportion and MPA of DF increased, water at the watershed level became less turbid. A similar relationship with the proportion of DF was not significant, possibly due to other sources of variability not considered in this study. On the other hand, the increase in PD, MPA, ED, CL, bare soil and GL increased turbidity significantly.
The pH was positively influenced by almost all landscape parameters related to DF fragmentation, including ED_FD and PD_FD. These last two indices give an idea of the level of degradation of the DF cover. All landscape indices favoring other types of LULC (configuration and composition), other than those favoring DF, tend to make the pH of stream waters less acidic than those watersheds with high proportions of forests. This would explain why the Isalowe and Lusambila watersheds have less acidic pH and more turbid stream waters with high temperatures (Figure 4). In Loweo, Bonde and Lubilaie it is the opposite. They have more acidic pH, less turbid stream waters and lower temperatures. Therefore, the turbidity, pH and temperature of river waters in Yangambi are all influenced by the LULC types at the watershed level in Yangambi.
In a 2015 study [32] on the characterization of the Congo River waters, experts associated the acidic increase in water quality with the predominance of humic acids in forest environments. This explains the characteristic dark red coloring-due to the presence of decomposing forest organic matter-of rivers within watersheds with a large forest area. Conversely, waters in highly anthropized environments (BSR, CL, GL and PF) are less acidic with a tendency to be alkaline. The latter are referred to as white water (Figure 11). level of degradation of the DF cover. All landscape indices favoring other types of LULC (configuration and composition), other than those favoring DF, tend to make the pH of stream waters less acidic than those watersheds with high proportions of forests. This would explain why the Isalowe and Lusambila watersheds have less acidic pH and more turbid stream waters with high temperatures (Figure 4). In Loweo, Bonde and Lubilaie it is the opposite. They have more acidic pH, less turbid stream waters and lower temperatures. Therefore, the turbidity, pH and temperature of river waters in Yangambi are all influenced by the LULC types at the watershed level in Yangambi.
In a 2015 study [32] on the characterization of the Congo River waters, experts associated the acidic increase in water quality with the predominance of humic acids in forest environments. This explains the characteristic dark red coloring-due to the presence of decomposing forest organic matter-of rivers within watersheds with a large forest area. Conversely, waters in highly anthropized environments (BSR, CL, GL and PF) are less acidic with a tendency to be alkaline. The latter are referred to as white water (Figure 11).
(a) (b) Figure 11. From left to right, (a) outlets of the Lusambila stream (water with a white color tendency) and (b) the Bonde stream (water with a reddish shade). The only significant SLR models that reflected the evolution of conductivity and dissolved oxygen with the change in landscape parameters in these watersheds were those of MPA_PF and MPA_C. Unfortunately, it has been difficult to explain these two relationships. Other sources of variability not considered in this study may give more insight. Producing meaningful SLR equations for all of these variables has been difficult.
Three SLR models best explained the variability of the analyzed physicochemical parameters of river water according to landscape parameters by referring to the AIC. These models explain the variation in pH as a function of ED_PF (94%), ED_crop land (92%) and PD_DF (90%) ( Table 7). This suggests that landscape fragmentation in Yangambi has a significant and already noticeable effect on stream water quality.

Water Vulnerability and Climate Change in Yangambi
As previously discussed in the results section, it appears that different water supplies are used in Yangambi. In the old settlements constructed around the houses of INERA workers, especially in Lusambila and Isalowe, the presence of cisterns and improved water springs is noticeable. In the other watersheds (Bonde, Loweo and Lubilaie), people obtain water mostly from rivers and springs, even if these sources are not improved. Indeed, the history of Yangambi city shows that water collection infrastructures (mainly cisterns), inherited from the colonial period, are exclusively located in certain watersheds within the administrative districts of the INERA Yangambi city. These settlements are much closer to urban than rural areas. A similar situation is observed almost everywhere in sub-Saharan Africa, where lack of infrastructure is a major problem for water accessibility, especially in peri-urban and rural areas [10]. This explains the discrepancies in terms of access to quality water between rural and urban areas [9]. Also, water management policies in DRC have not paid sufficient attention to rural environments in the past. Hence, few initiatives in the framework of the establishment of hydraulic infrastructures by the NRHS or the Régie de distribution d'eau have been implemented in these areas [4]. As a result, a significant proportion of rural areas do not have water supply infrastructure. Certain characteristics of settlements in the different watersheds can explain this variability of the resource. These characteristics would explain the great diversity of water supply sources in such a small area as Yangambi.
The spatial distribution of surveyed water points using GIS shows that water sources vary in the five watersheds. Although there are several water points in Yangambi, some are in the middle of the forest (in the YRB), far away from houses ( Figure 8). People use water points closest to their homes. Thus, some watersheds are more accessible than others, depending on the available hydraulic possibilities. These hydraulic possibilities can be evaluated in terms of the number of water points close to settlements in the watersheds. This means that accessibility is good, especially in watersheds with cisterns, as they are close to houses. At the same time, in some villages, people walk more than 1 km to get water from springs or rivers. In other settlements, where the springs are far from homes and there are no cisterns, people depend exclusively on river water for all their needs, including drinking. For example: in Obiloto in the Bonde watershed, people depend mainly on river water for all their domestic needs; use in Isalowe varies between cisterns, springs and rivers. This situation affects the accessibility and the amount of water used per person per day in the different watersheds. That is why Isalowe and Lusambila are among the watersheds with the highest scores in terms of water use and accessibility.
The surveys revealed that water resource management capacity in Yangambi is generally low, except for the Isalowe watershed, which exceeds the 50% score ( Figure 10). This situation is observed in many rural areas of developing countries [33,34]. In most of these countries, as well as in Yangambi (DRC), water is available, but the capacity to manage the resource is lacking [20]. Indeed, the Congolese population in rural areas lacks access to information, education and the basic infrastructure essential to their well-being [6]. Settlements in Isalowe and Lusambila have more capacity to manage the water resource than other settlements and villages; for example, they host some INERA employees, who have a certain social status.
Similarly, in terms of the water resource use component (amount of water consumed per person per day), the Isalowe watershed has the highest score. This advantage is probably due to the water collection infrastructure-especially the cisterns, and the high density of springs and rivers close to homes ( Figure 5)-facilitating water accessibility and use. However, the same characteristics that favor Isalowe and Lusambila for water resources are at the root of their high population density. This puts strong pressure on the Yangambi forest compared to the Loweo, Lubilaie and Bonde catchment areas, whose population density is much lower. Based on different ways to assess the WPI, the vulnerability of water resources in the different watersheds varies in Yangambi. That is to say, the different socioeconomic and environmental characteristics of each village or settlement in the studied watersheds influence the different aspects of water resources in Yangambi. Other researchers have confirmed these results [21,22,34,35], showing that many parts of the world have high variability of water resources at different scales. The analysis of water resource variability based on meteorological parameters partly confirms that of Aguilar and colleagues in 2009 [36]. Based on climate data for 1955-2006, this study showed that, overall, the East Central African region had experienced a significant decline in total annual rainfall. Jemmali (2018) [37] published similar results. It identified irregularities and seasonality of rainfall as important components of the sub-Saharan climate responsible for drought and flooding in the region. The availability and accessibility of water in this region is often compromised when irregularities and seasonality are not considered. In turn, the global climate predictions in the fifth report of the Intergovernmental Panel on Climate Change and the Commission Internationale du Bassin Congo-Oubangui-Sangha [32] predict that temperatures will rise over this century, and that rainfall disruptions and prolonged periods of drought will be more frequent. Figure 12 summarizes the vulnerability impact chain of water resources in Yangambi.
responsible for drought and flooding in the region. The availability and accessibility of water in this region is often compromised when irregularities and seasonality are not considered. In turn, the global climate predictions in the fifth report of the Intergovernmental Panel on Climate Change and the Commission Internationale du Bassin Congo-Oubangui-Sangha [32] predict that temperatures will rise over this century, and that rainfall disruptions and prolonged periods of drought will be more frequent. Figure 12 summarizes the vulnerability impact chain of water resources in Yangambi.

Conclusions and Recommendations
This study aimed to analyze the vulnerability of domestic water resources in the Yangambi context. We based the assessment of the impact of LULC types on the physicochemical quality parameters of stream water and whether the population would be able to adapt to this vulnerability in the context of climate change. The physicochemical parameters of stream water within the watersheds showed that water for domestic use no longer meets WHO standards for drinking water. This situation is amplified by the conversion of DF cover to other types of land use (linked to anthropization). These practices, in turn, deteriorate physicochemical parameters such as turbidity, pH and temperature of stream water. Unfortunately, the population in Yangambi has a low capacity to manage this resource, which results in high vulnerability for climate change adaptation. This vulnerability is already apparent through climate parameters recorded at the INERA Yangambi meteorological station based in this area. Therefore, we recommend an analysis of the state of the water resource by watershed to determine the best options for water resource management in Yangambi. At the same time, an exhaustive hydrological data sampling is needed to understand the different types of LULC contributions to the variability of water resource quality in Yangambi. The installation of gauging stations at different rivers in Yangambi will be important for better monitoring of the water resource. Finally, the installation of infrastructures that can help supply water to the communities while preserving the DF within the YBR is a priority, as water quality and quantity highly depend on it. A qualitative and quantitative evaluation of the impact of LULC change on water resources' quality and availability would be interesting. It could consider detailed analyses of physicochemical (organic ions, major inorganic ions, trace elements, etc.) and biological parameters of water, as well as the variation of streamflow data within the YBR catchments.