Earth Observation and Cloud Computing in Support of Two Sustainable Development Goals for the River Nile Watershed Countries

: In September 2015, the members of United Nations adopted the 2030 Agenda for Sustainable Development with universal applicability of 17 Sustainable Development Goals (SDGs) and 169 targets. The SDGs are consequential for the development of the countries in the Nile watershed, which are a ﬀ ected by water scarcity and experiencing rapid urbanization associated with population growth. Earth Observation (EO) has become an important tool to monitor the progress and implementation of speciﬁc SDG targets through its wide accessibility and global coverage. In addition, the advancement of algorithms and tools deployed in cloud computing platforms provide an equal opportunity to use EO for developing countries with limited technological capacity. This study applies EO and cloud computing in support of the SDG 6 “clean water and sanitation” and SDG 11 “sustainable cities and communities” in the seven Nile watershed countries through investigations of EO data related to indicators of water stress (Indicator 6.4.2) and urbanization and living conditions (Indicators 11.3.1 and 11.1.1), respectively. Multiple approaches including harmonic, time series and correlational analysis are used to assess and evaluate these indicators. In addition, a contemporary deep-learning classiﬁer, fully convolution neural networks (FCNN), was trained to classify the percentage of impervious surface areas. The results show the spatial and temporal water recharge pattern among di ﬀ erent regions in the Nile watershed, as well as the urbanization in selected cities of the region. It is noted that the classiﬁer trained from the developed countries (i.e., the United States) is e ﬀ ective in identifying modern communities yet limited in monitoring rural and slum regions. seven Nile watershed countries. Multiple SDG indictors are studied including: the water stress (6.4.2), land use and urban

With an estimated length ranging between 5499 km and 7088 km, it is the longest river in Africa and widely regarded as the longest river in the world [1]. The Nile serves as the primary water source to downstream nations namely, Sudan and Egypt [2]. However, water stress/scarcity are expected because of limited water to meet the various irrigation needs of the Nile Basin nations, as well as hydraulic engineering projects among the countries such as hydropower dams being built in Sudan (Merowe dam), Ethiopia (Grand Renaissance dam) and Uganda (Karuma dam) and other stressors driven by climate change. It is noted that these projects will not reduce the tensions between countries but possibly introduce more challenges, bringing the uncertainty of the future economic development and regional stability. Additionally, the local water measurements of the Nile are outdated and limited for scientific studies. For example, the latest records of stations in the Global Runoff Data Centre (https://www.bafg.de/GRDC/EN/01_GRDC/grdc_node.html) dataset date back to 1984.
Countries in the Nile are undergoing rapid population growth and urbanization. Seven Nile Basin countries are expected to double their populations between 1995 and 2025 [3]. Most of the countries have agriculturally-based economies, and rapid population growth will further pressure natural resources and often contribute to environmental stress. The limited financial resources and technological abilities of the majority of these nations result in them being unable to properly take mitigation actions. Additionally, rapid urbanization is anticipated in the Nile countries due to the increasing movement of people from rural to urban areas associated with population growth and poverty. On the one hand, the rapid urbanization is driven by the labor migration and it possibly coincided with process of industrialization [4]. On the other hand, the unplanned urban regions often lack the infrastructure and institutions needed to supply adequate water and sanitation, as well as reliable housing and public transportation, in order to protect human and environmental health. The migrated populations form or join the low-income communities (slum urbanizations) that are especially impacted by inadequate environmental services and by water and air pollution. The situation is typical in the Nile countries, since the ongoing South Sudanese civil war has sent at least over 2 million people migrating to neighboring countries as refugees [5]. The uncontrolled growth (from 245,000 in 1956 to more than 7 million in 2012) of Sudan's capital Khartoum has resulted in 50% of residents are living in informal and squatter settlements on the outskirts of the city [6].
On the historical 70th anniversary of the United Nations (UN) in September 2015, the heads of state and delegates gathered and adopted the 2030 Agenda for Sustainable Development [7]. This comprehensive sustainability framework was built on the historical experiences of human society and a shared expectation for the next 15 years. The 17 Sustainable Development Goals (SDGs) incorporate various social, economic, environmental, and developmental aspects, which consist of 169 targets and 230 indicators [8], and have been endorsed by all countries with respective national implementation plans for their different and very diverse development context. It can be concluded that the SDGs constitute a vast system that is complicated, diverse, dynamic, and interconnected, where success in achieving one goal is often linked to solving issues associated with other goals. However, according to SDG reports from the UN [9] and Chinese Academy of Sciences [10], the lack of accurate and timely data, as well as the lag in developing proper methodologies to process the data, has become a bottleneck that the full implementation of the 2030 SDGs will be hampered if these issues are not effectively resolved. Therefore, in order to achieve the SDGs and effectively assess their progress with the full strength of science, technology, and innovation (STI), the UN has established the "Technology Facilitation Mechanism" (TFM) [11]. Presently, a pressing priority of TFM is to make breakthroughs towards Tier II indicators (methods established but with poor data) and Tier III indicators (methods under development with either poor data or no data).
To support and track the progress in fulfilling SDGs, the UN suggested that Earth Observation (EO) and geo-spatial data should be utilized to mitigate the shortage of in-situ data availability [12]. EO provides reliable, updated and cost-effective data resources [13], which have been applied to monitor the SDG-related parameters from regional to global scales on a consistent and regular basis [14][15][16][17][18]. Nevertheless, the EO datasets are also characterized by massive quantity, multiple sources, multi-temporal, multi-dimensional heterogeneous structure, and high complexity. The study of EO data in support of SDGs, therefore, requires powerful computing capabilities and advanced analytical methods.
Recently, the advancement of cloud computing, including parallel computing, offer near unlimited computer processing capabilities, as well as free access to a large volume of EO data already stored on the remote cloud drives. There has been an emerging increase in high-performance cloud computing platforms, such as the National Aeronautics and Space Administration (NASA) Earth Exchange (NEX), Amazon Web Service (AWS) and Google Earth Engine (GEE) [19]. The GEE directly accesses the Google's cloud resources to undertake all the processing and analysis, which does not require large processing powers of computers or expensive software. It is characterized as an easily accessible and user-friendly front-end platform providing a convenient environment for interactive data representation and algorithm development. For instance, Donchyts et al. [20] employed the GEE for mapping global surface water changes over the past 30 years with a high spatial resolution (30 m), an effort that would not have been possible without the powerful processing and analytical capabilities from cloud computing. GEE is also widely used for environmental research in the Nile watershed and adjacent regions [14,[21][22][23][24][25]. In particular, it is a promising tool that allows researchers from developing countries (e.g., Nile Basin countries) to have the same ability to undertake state-of-art studies as those in the most advanced nations. Considering the urgent situation of water data availability and urbanization stress in the Nile Basin countries, the EO and cloud computing evidently appear to be one of the priorities to help achieve the SDGs.
This paper demonstrates the intersection between EO, cloud computing and related technologies and demonstrates how they contribute to the assessment of two SDGs in the Nile watershed countries, namely SDG 6 "clean water and sanitation" and SDG 11 "sustainable cities and communities", with focus on the water stress situation and urbanization process. This is carried out through studying the water resources using multiple hydrological variables, as well as the change detection of the land cover in five major Nile Basin cities using advanced algorithms implemented on a cloud-computing platform.

Study Area
In this study, we investigate seven countries in the River Nile watershed, namely, Egypt, Sudan, South Sudan, Ethiopia, Uganda, Kenya and Tanzania, which totally occupy up to 97.91% of the areas of the watershed calculated from the shapefile accessed from the Nile Basin Initiative (https://www.nilebasin.org) ( Table 1). Sudan has the biggest size (33.92%) compared to other countries, while both South Sudan and Uganda have their entire territories within the watershed (left of Figure 1). In Section 2.2, we will study the water stress situation in the countries mentioned above. Furthermore, in Section 2.3, we will investigate the changes of impervious surface in five specific capital cities, over the recent years, in some of these countries. The five capital cities, in this region, targeted here are namely: Cairo, Khartoum, Juda, Addis Ababa and Kampala.

Data and Methods to Study Sustainable Development Goal (SDG) 6 Indicator on Water Stress
Target 6.4 of SDG 6 is to "substantially increase water-use efficiency across all sectors and ensure sustainable withdrawals and supply of freshwater to address water scarcity and substantially reduce the number of people suffering from water scarcity." Its Indicator 6.4.2, "Level of water stress: freshwater withdrawal as a proportion of available freshwater resources", is explicitly used for monitoring water stress situations at the country level. This indicator (water stress) has been defined as the ratio percentage (%) between total freshwater withdrawn (TFWW) by all major sectors and total renewable freshwater resources (TRWR), after having taken into account environmental water requirements (EFR). It is calculated with the following equation: A low level of water stress shows a situation where the combined withdrawal is marginal in proportion to the total renewable freshwater resources and has, therefore, little potential impact on the sustainability of the resources or on the potential competition between users. A high level of water stress indicates an opposite situation where the combined withdrawal represents a substantial share of the resources, with potentially larger impacts on the sustainability of the resources and potential

Data and Methods to Study Sustainable Development Goal (SDG) 6 Indicator on Water Stress
Target 6.4 of SDG 6 is to "substantially increase water-use efficiency across all sectors and ensure sustainable withdrawals and supply of freshwater to address water scarcity and substantially reduce the number of people suffering from water scarcity." Its Indicator 6.4.2, "Level of water stress: freshwater withdrawal as a proportion of available freshwater resources", is explicitly used for monitoring water stress situations at the country level. This indicator (water stress) has been defined as the ratio percentage (%) between total freshwater withdrawn (TFWW) by all major sectors and total renewable freshwater resources (TRWR), after having taken into account environmental water requirements (EFR). It is calculated with the following equation: A low level of water stress shows a situation where the combined withdrawal is marginal in proportion to the total renewable freshwater resources and has, therefore, little potential impact on the sustainability of the resources or on the potential competition between users. A high level of water stress indicates an opposite situation where the combined withdrawal represents a substantial share of the resources, with potentially larger impacts on the sustainability of the resources and potential situations of conflicts and competition between users, resulting in negative effects on economic development and food security. Since 1994, the Food and Agriculture Organization of the United Nations (FAO) has been monitoring the parameters of water stress through its global water information system, AQUASTAT [26]. Table 2 demonstrates the water stress indicator and related parameters of the Nile countries between 2007 and 2017 from AQUASTAT database (http://www.fao.org/nr/water/aquastat/main/index.stm).  TRWR is calculated as being the sum of total internal and external renewable water resources, also equal to total renewable groundwater and surface water: Total internal renewable water resources (IRWR) is calculated as the sum of internally produced surface water and groundwater with deducting their overlap: Total external renewable water resources (ERWR) is calculated as being the sum of external renewable surface water and groundwater: Remote Sens. 2020, 12, 1391 6 of 25 TFWW is usually calculated as being: the sum of total water withdrawal with deducting the direct use of wastewater, direct use of agricultural drainage water and use of desalinated water. TFWW can be also calculated as the sum of total withdrawal of fresh surface water and groundwater: Groundwater usage and recharge. All the countries, except Egypt, have their groundwater recharge internally. Egypt has total renewable groundwater (TRGW) of 1.5 × 10 9 m 3 /year yet with higher FGWW of 7.5 and 6.5 × 10 9 m 3 /year for the year 2012 and 2017, respectively. The decreased FGWW contributes to a lower TFWW that mitigates the water stress value from 124% to 117.3%.

•
Overgeneralized spatial and temporal data. The country level data at a five year period overlooks the distribution of the water stress within the country and may hide specific situations that only exist for part of the year.
Therefore, the current data-collection process of Indicator 6.4.2 is not comprehensive and can be assisted by the EO data. In this study, the EO hydrological data is used to show the spatial distribution of water resources among the Nile Basin countries in a higher temporal scale. The tasks are performed to: (1) show the changing trend of water resources related to water stress situations (e.g., precipitation, evapotranspiration, groundwater, etc.); (2) present the connection between groundwater and hydrological parameters related to surface water; (3) compare the groundwater withdrawal and renewable groundwater from natural resources; and (4) explore the spatial and temporal recharge patterns of groundwater from internal and external resources. The hydrological datasets is introduce in the Section 2.2.1.

Hydrological Datasets
• NASA-United States Department of Agriculture (USDA) Global Soil Moisture Dataset The NASA-USDA Global soil moisture dataset [27][28][29][30] includes information of surface and subsurface soil moisture, soil moisture profile (%), as well as surface and subsurface soil moisture anomalies represented as unitless and standardized values (negative/positive either refers to drier/wetter than normal conditions, respectively) computed using a 31-days moving window. The dataset has spatial resolution of 0.25 • × 0.25 • at global scale and temporal resolution of 3 days starting from 1 January 2010. This dataset is generated through the data assimilation approach, which employs a modified two-layer Palmer model using a 1-D ensemble Kalman filter (EnKF) to integrate the satellite-derived Soil Moisture Active Passive (SMAP) Level 3 soil moisture observations. This approach turns out to be helpful to improve the model-based soil moisture predictions, particularly for the Nile watershed region, which lacks good quality hydrological measurements. The unit of soil moisture is in mm.

• Climate Hazards Group Infrared Precipitation with Station Dataset
The precipitation dataset is obtained from the Climate Hazards Group Infrared Precipitation with Station (CHIRPS) dataset (version 2.0) [31]. CHIRPS is a 30+ year quasi-global rainfall dataset, which is developed by the United States Geological Survey (USGS) in collaboration with the Earth Resource Remote Sens. 2020, 12, 1391 7 of 25 Observation and Science Center (https://www.usgs.gov/centers/eros). The dataset incorporates 0.05 • resolution satellite imagery with in situ station data with temporal resolution as pentad (i.e., the first 5 pentads in a month have 5 days and the last pentad contains all the days from the 26th to the end of the month) to create gridded precipitation time series for trend/correlational analysis and drought/flood monitoring.

• Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS) Dataset
The evapotranspiration (ET) and runoff data is obtained from the FEWS NET FLDAS [32,33]. The FLDAS is designed particularly to assist the food security investigations for the developing countries with sparse data accessibility. In this research, the version of FLDAS dataset used in this study is the Noah version 3.6.1 surface model with downscaled CHIRPS-6 hourly precipitation inputs using the NASA Land Surface Data Toolkit (https://lis.gsfc.nasa.gov/software/ldt). The model uses a combination of the new version of Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data and CHIRPS as the forcing data. The dataset provides monthly average value with spatial resolution of 0.1 arc degrees. Besides, the total runoff is computed from the sum of surface storm runoff (Qs_tavg) and baseflow-groundwater runoff (Qsb_tavg). The unit of the FLDAS data is translated from kg m −2 s −1 to mm s −1 (1 kg m −2 = 1 mm H 2 O).
• Gravity Recovery and Climate Experiment (GRACE) Mission Monthly Gravitational Anomalies GRACE Tellus Monthly Mass Grids [34][35][36] provides monthly gravitational anomalies relative to a 2004-2010 time-mean baseline. The dataset uses the units of "equivalent liquid water thickness" (LWE) which represent the deviations of mass in terms of vertical extent of water in centimeters. The GRACE mission operated from April, 2002 to June, 2017. The GRACE satellites measure the spatial and temporal variability in the Earth's gravity field and the monthly gridded terrestrial water storage (TWS) products are available from the NASA Jet Propulsion Laboratory (JPL) Tellus website (https://grace.jpl.nasa.gov/data/get-data/monthly-mass-grids-land/) on 1 horizontal resolution grids and global coverage. The unit of the GRACE data is also translated from centimeters to mm.

Multi-Assessment Methodologies of Water Stress Addressing SDG 6
Harmonic analysis is a method involving the representation of functions or signals as a superposition of elementary waves. In this study, to estimate the variation of multiple hydrological factors (e.g., precipitation), we build the harmonic model H(t) with elements of a constant band (β 0 ), a linear term of slope (β 1 ) and harmonic terms of amplitudes (β 2 , β 3 , β 4 and β 5 ). The term β 1 , associated with the linear part of the factor, represents the increasing/decreasing trend, whereas a constant band β 0 represents the extent of consistency of the time series. Moreover, f represents the fundamental frequency. The β 1 of the harmonic analysis can show linear trend of a series regardless of the seasonal variations. The t is the record for each parameter including monthly LWE, total monthly precipitation, total monthly water lost from evapotranspiration and runoff, as well as averaged monthly surface soil moisture. All the data used are accessed from the Google Earth Engine platform through sampling of random points (200 to 300 points according to the country size) within the overlapped region of both Nile watershed and countries' territories. For example, the sample points used in the Egypt region are presented in Figure 1. In this research, the value β 1 shows the monthly trend in the Nile watershed regardless seasonal variations.. Positive value presents increasing trend, and negative value decreasing.
The precipitation, ET, surface soil moisture, runoff and LWE products are used to generate the changing trend map of the Nile watershed region during May 1st, 2013 till October 20th, 2019 using month values: (1) precipitation: total precipitations (in millimeters) within a month; (2) ET and runoff: the monthly averaged value (in mm s −1 ) multiplies by the total seconds of 30 days (2,592,000 seconds); (3) soil moisture: the monthly averaged value (in millimeters); and (4) LWE groundwater: the monthly records (in millimeters). In addition, the relationship between the parameters and groundwater recharge pattern is studied by the lag correlation analysis using their anomalies value X a calculated as Equation (7).
with X as the monthly value and X as the monthly mean value.
In the previous comparative global study of water availability in the major river basins [37], the basin averaged water balance was estimated through Equation (8) and compared to the GRACE water LWE anomaly through Equation (9).
where precipitation (P), ET, runoff (R), and change in surface and subsurface water storage ∆S are monthly and basin-averaged values for each of the basins. The left-hand variables of Equation (9) can be estimated from the satellite observations or model outputs, while the right-hand ∆S is estimated thought the GRACE observations (∆S GRACE ), as well as other factors related to human activities (∆S human ), such as of exploitation of groundwater in the basin for agricultural or industrial purposes.

Data and Methods to Study SDG 11 Indicator on Urbanization Process
Target 11.3 of SDG 11 is to "enhance inclusive and sustainable urbanization and capacity for participatory, integrated and sustainable human settlement planning and management in all countries." Its Indicator 11.3.1, "Ratio of land consumption rate to population growth rate" (LCRPGR) shown in Equation (10), is defined as the ratio of land consumption rate (LCR) to population growth rate (PGR) and is used to describe the urbanization process and population growth.

LCRPGR =
land consumption rate annual population growth rate = LCR PGR (10) In particular, it is discovered that urban impervious surfaces retrieved from the remote-sensing imagery can effectively reflect urban surface and land use information [10]. Some countries (e.g., US [38], China [39]) have developed regional high-resolution urban land mapping datasets. However, many developing countries do not have the technical and financial resources to monitor their urban development. Therefore, the goal of this study is to develop a transfer learning method using freely assessable EO and land cover datasets, as well as computing resources to monitor the impervious surface changes of urban cities in the Nile Basin countries. This study will enable the developing countries to investigate the relationship between land use and population growth rate in urban environment. In addition, we also evaluate the effectiveness of this method to target slums urbanization situation, in order to monitor Indicator 11.1.1 "proportion of urban population living in slums, informal settlements or inadequate housing".

Modeling of Impervious Surfaces Addressing SDG 11
The modeling approach in this study is built on the infrastructure provided by Google's full-stack geospatial modeling and analysis platform: (1) Google Earth Engine as the geospatial analysis platform; (2) Earth Engine Data Catalog as the geospatial data archive that can be accessed and analyzed in the cloud computing platform; (3) TensorFlow as the machine learning platform to implement model design; (4) artificial intelligence (AI) Platform as to build and run the models and (5) Colab as the Jupyter notebook server for the workflow development. Firstly, the Earth Engine requests the training data from the Earth Engine Data Catalog stored in the cloud storage. Secondly, the cloud storage finds the training data and sends it to the AI Platform, while the AI Platform receives the model designed by the TensorFlow and its application programming interfaces (APIs) programmed in the Colab server. When the model is built completely, GEE can access the model to make predictions through its new functionality API of ee.Model.fromAIPlatformPredictor and output the results. This update of GEE's API to have new integration with Google Cloud's AI Platform expands its modeling scalability and enhances its analytical capabilities. The traditional APIs in the GEE include the machine-learning algorithms such as decision trees, K-means and support vector machine. However, these methods have limited capability to extract information from high-dimensional (number of bands >>3) imagery about not just the spectral values from one pixel, but spatial patterns of many pixels (i.e., the spatial context) represented as big patches of pixels (e.g., 256 × 256). Therefore, the fully convolutional neural networks (FCNNs) [41] are designed to handle these difficulties. FCNN contains one or more convolutional layers, in which inputs are neighborhoods of pixels. Therefore, it learns from patches of data, instead of single pixels, resulting in a network that is not fully connected, but is suited to identifying spatial patterns. In this study, the FCNN model is trained from training data from NLCD dataset and predicts a continuous range from 0 to 1 output in each pixel from 256 × 256 neighborhoods of pixels to show the percentage of the impervious surface levels. The example provided in the Colab document (https://colab.research.google.com/github/google/earthengine-api/blob/master/ python/examples/ipyn-b/UNET_regression_demo.ipynb) shows the details of how to export patches of data to train the network and how to overtile image patches for inference, to eliminate tile boundary artifacts. The Landsat 8 surface reflectance data has been rescaled to range from 0 to 1 for both seven reflectance bands and the two thermal temperature bands were between 0 • C to 100 • C. Therefore, the processed Landsat 8 composite image in 2016 is ready to be trained with the NLCD impervious surface as an objective field. The training patches are obtained from cities across the continental US shown in Figure 2a (red as training regions and blue as evaluation regions). Total 20,000 and 8000 patches are used for training and evaluation, respectively, including a total 50 epochs (one epoch is a full processing by the learning algorithm over the entire training set). The model uses the Stochastic gradient descent (SGD) and momentum method as the optimizer, the mean of squares of errors as the loss function, and the root mean squared error as the metrics. The model is optimized by minimizing the loss and monitoring the metrics values at the end of each epoch. The modelling process took 1 day 6 hour and 57 mins, and the process is displayed in the Figure 2b.         Figure 6 shows the correlation maps of correlation coefficient (Corr) and p-value between P-ET-R and LWE up to 3 months lags apart. The Corr value for the regions with p-value higher than 0.05 (red color) is not considered plausible. In the case of lag = 0, the positive relationship is found near the south of Egypt. South Sudan, some parts of Ethiopia and Uganda demonstrate a negative relationship, which may be due to the coincidence between high LWE and low precipitation seasons. In the case of lag = 1, only some regions of Egypt, Sudan and Ethiopia show a weaker positive relationship between P-ET-R and LWE. However, for lags of 2-3 months, the strong positive connections between P-ET-R and LWE are illustrated in the most of the regions in South Sudan, Ethiopia and southern part of Sudan, indicating the recharge process to the groundwater. This result also matches the observed relationship between the peaks of P-ET-R and LWE in Sudan, South Sudan and Ethiopia in Figure 5.  Figure 6 shows the correlation maps of correlation coefficient (Corr) and p-value between P-ET-R and LWE up to 3 months lags apart. The Corr value for the regions with p-value higher than 0.05 (red color) is not considered plausible. In the case of lag = 0, the positive relationship is found near the south of Egypt. South Sudan, some parts of Ethiopia and Uganda demonstrate a negative relationship, which may be due to the coincidence between high LWE and low precipitation seasons. In the case of lag = 1, only some regions of Egypt, Sudan and Ethiopia show a weaker positive relationship between P-ET-R and LWE. However, for lags of 2-3 months, the strong positive connections between P-ET-R and LWE are illustrated in the most of the regions in South Sudan, Ethiopia and southern part of Sudan, indicating the recharge process to the groundwater. This result also matches the observed relationship between the peaks of P-ET-R and LWE in Sudan, South Sudan and Ethiopia in Figure 5.     Table 3 shows the ratio of highly impervious surface (value >50%) areas for the years of 2013 and 2019 and the changes ∆. Three cites (Khartoum, Addis Ababa and Kampala) illustrate an expansion in the percentage of impervious surface to some extent, with additional red or yellow colored regions. For Cairo, the ratio unexpectedly decreased since the additional impervious surface areas are seen at the eastern region at the New Cairo City. This may be due to the classifier, which cannot identify some regions of human habitats as it is observed that massive  To further explore the cause of misclassification of the impervious surface in the Cairo region, Figure 8 is created to show the surface reflectance signatures of the samples (selected locations shown on the left of Figure 8) in four types of land categories, namely, human habitats not marked as high impervious surface (Unmarked Habitat), human habitats marked as high impervious surface (Marked Habitat), barren areas (Barren) and vegetation areas such as croplands or grasslands (Vegetation). The barren surface shows a highest surface reflectance over B2-B7, while vegetation surface has high B5 (Near Infrared) in contrast to B4 (Red). However, the unmarked and marked habitats only have slight differences of surface reflectance values in the bands B6 (Shortwave Infrared 1) and B7 (Shortwave Infrared 2). Meanwhile, the values of marked habitats have wider confidence intervals than the unmarked in each band.  To further explore the cause of misclassification of the impervious surface in the Cairo region, Figure 8

Discussion of the Results of the SDG 6 Study
This research investigated hydrological parameters related to water stress situations in the Nile Basin countries, especially highlighting the usage of EO to solve the spatial and temporal data shortage and inconsistency related to the renewable freshwater from both surface water and groundwater resources. The changing trends of multiple hydrological parameters related to the surface water resources and their relationships with groundwater has been demonstrated (Figures 3  and 4). The groundwater recharge from the surface water also has been shown in Figure 5 and Figure  6 in the time-series and cross-correlation analysis. The impact on groundwater recharge under various climate change scenarios has been investigated by previous studies [42][43][44], which commonly concluded that, globally, the amount of groundwater is decreasing, especially in the Northern Africa and Middle East regions. In Figure 3, the decreasing trend of LWE in Egypt and part of Sudan confirms what was found in the previous study that the basin sinks with the negative water balance in the downstream Nile countries [45]. For the Nile Basin, our results of ΔS as P-ET-R show a slightly different pattern compared with what was presented in the previous study. This may be due to the choices of different hydrological datasets (e.g. Tropical Rainfall Measuring Mission (TRMM) vs. CHIRPS datasets for precipitation, Global Land Data Assimilation System (GLDAS) vs. FLDAS for runoff, Moderate Resolution Imaging Spectroradiometer (MODIS) vs. FLDAS for ET). As the main consumption of the freshwater, the agriculture sector consumes over 80% of the total water resources in Egypt [46]. Meanwhile, the abstraction from groundwater grows rapidly with the expansion of irrigation activities, industrialization and urbanization [46]. In Figure 3, the precipitation, ET and runoff are stable for the northern part of the Nile watershed but with decreasing LWE. In Figure 5 Table 2. Therefore, the loss of the groundwater is very likely due to overexploitation. Consequently, the decreasing LWE may correspond to the impacts of abstraction of groundwater to meet increasing agricultural and other needs. Moreover, this study also directly compared the P-ET-R monthly anomalies values and the LWE values through correlational analysis ( Figure 6). The approach used both time series analysis,

Discussion of the Results of the SDG 6 Study
This research investigated hydrological parameters related to water stress situations in the Nile Basin countries, especially highlighting the usage of EO to solve the spatial and temporal data shortage and inconsistency related to the renewable freshwater from both surface water and groundwater resources. The changing trends of multiple hydrological parameters related to the surface water resources and their relationships with groundwater has been demonstrated (Figures 3 and 4). The groundwater recharge from the surface water also has been shown in Figures 5 and 6 in the time-series and cross-correlation analysis. The impact on groundwater recharge under various climate change scenarios has been investigated by previous studies [42][43][44], which commonly concluded that, globally, the amount of groundwater is decreasing, especially in the Northern Africa and Middle East regions. In Figure 3, the decreasing trend of LWE in Egypt and part of Sudan confirms what was found in the previous study that the basin sinks with the negative water balance in the downstream Nile countries [45]. For the Nile Basin, our results of ∆S as P-ET-R show a slightly different pattern compared with what was presented in the previous study. This may be due to the choices of different hydrological datasets (e.g., Tropical Rainfall Measuring Mission (TRMM) vs. CHIRPS datasets for precipitation, Global Land Data Assimilation System (GLDAS) vs. FLDAS for runoff, Moderate Resolution Imaging Spectroradiometer (MODIS) vs. FLDAS for ET). As the main consumption of the freshwater, the agriculture sector consumes over 80% of the total water resources in Egypt [46]. Meanwhile, the abstraction from groundwater grows rapidly with the expansion of irrigation activities, industrialization and urbanization [46]. In Figure 3, the precipitation, ET and runoff are stable for the northern part of the Nile watershed but with decreasing LWE. In Figure 5 Table 2. Therefore, the loss of the groundwater is very likely due to overexploitation. Consequently, the decreasing LWE may correspond to the impacts of abstraction of groundwater to meet increasing agricultural and other needs. Moreover, this study also directly compared the P-ET-R monthly anomalies values and the LWE values through correlational analysis ( Figure 6). The approach used both time series analysis, correlation maps (up to 3 months lags) to show the spatial water transport, and recharge pattern among different regions in the Nile watershed.
In Figure 6, southern Egypt has shown a positive relationship between P-ET-R and LWE, indicating possible external groundwater recharge from the upstream water resources to Egypt since it is listed as the only country receiving external groundwater in the Nile Basin countries (Table 2). It is also noted that the magnitude of the P-ET-R is much larger for water-rich basins such as Amazon and Mekong than Nile [37]. The study validates the use of P-ET-R to show the spatial and temporal variability as it reflects the increase or decrease in the water storage in coordination with the GRACE dataset, especially for some countries such as Ethiopia, Sudan and South Sudan ( Figure 5). Moreover, the accurate estimation of water stress situations is also significantly related to the freshwater withdrawn (TFWW) by all major sectors, as well as environmental water requirements (EFR). Table 2 shows that the missing TFWW data causes water stress unable to be calculated. Meanwhile, the unchanged EFR data indicates consistent environmental flow requirements between 2007 to 2017, which is questionable for the recent changes of their natural and social conditions (dams development, population growth, etc.). Therefore, the specific studies on collecting high-quality TFWW and EFR data are suggested to be undertaken for the actual report of water stress situations. The study of water stress also helps water resources management cooperation between Nile Basin countries in support of SDG Target 6.5 "implement integrated water resources management at all levels, including through transboundary cooperation as appropriate".

Discussion of the Results of the SDG 11 Study
For Indicator 11.3.1, the impervious surface maps generated from the research provides updated information of the land cover that can be used to monitor the population distribution and the quality of their settlement as a result of urban expansion of the selected countries ( Figure 7). However, the model trained from the US data sources cannot identify some parts of the human habitats in Cairo and Khartoum as impervious surfaces. The optical spectral signatures in Figure 8 only show slightly differences between the marked and unmarked regions in terms of the confidence levels of the averaged surface reflectances, which cannot clearly discriminate different surface types. It is known that the model does not simply use the spectral differences as the parameters but also includes the imagery context or texture. Therefore, it is possible that the training data from US cities lack certain texture or types of impervious surface, which are common in Cairo or other cities. It is suggested that both unmarked regions in Cairo and Khartoum are the towns with similar constructions and homogeneous texture and materials. It is also found that for the entire Nile Delta region, many towns or villages are also classified as non-impervious surfaces (Figure 9 regions A-E), where the percentage value of impervious surfaces is even lower than surrounding crop fields (also seen in Figure 7). This indicates the some non-impervious surface regions are still human settlements, while the constructions are built by the permeable materials and occasionally shaded from higher vegetation (e.g., palm trees).
Additionally, for Indicator 11.1.1, the difference between the low-income community (rural and slum areas) (characterized with homogeneous houses and imagery textures) and high-income community (chartered with high urban greening and diverse buildings) can alter the model outputs.
For instances, Figure 10 shows the classified impervious surface of Karachi, Pakistan, and its slum Orangi Town (red square region), which is regarded as the largest slum in the world with an estimated population from 1.5 to 2.4 million [47]. It is noted that the FCNN model shows more details of the distribution of the land cover than the MODIS product. However, the MODIS product classifies Orangi Town as impervious surface yet the model does not mark Orangi Town but marks the adjacent region with rich heterogeneous textures (right bottom corner cross the red line) as impervious surface. Orangi Town has houses of similar structure and homogeneous texture. It is reasonable to speculate that such typical texture of impervious surface, as well as the building materials of the houses and villages (such as sand, woods or mud) in the arid and semi-arid regions of Egypt and Pakistan are not typical and included in the US training datasets. However, such high-resolution impervious surface datasets are only limited in the US cities. This situation reflects a data inequality between developed and developing countries and that even the results of a modeling approach to help investigate the  In the model training experience of this study, the Google Cloud Platform (GCP) charged $25.99 to generate the model within 31 hours. It is a relatively inexpensive for researchers with limited financial capacity, considering GCP also provides $300 credit for new users which covers all the cost.  In the model training experience of this study, the Google Cloud Platform (GCP) charged $25.99 to generate the model within 31 hours. It is a relatively inexpensive for researchers with limited financial capacity, considering GCP also provides $300 credit for new users which covers all the cost.