Assessing The Effect of Urbanization On Regional-S
Assessing The Effect of Urbanization On Regional-S
Assessing The Effect of Urbanization On Regional-S
com/scientificreports
Water supply demand and water resources pollution are consistently increasing related to population growth,
economic development, and lifestyle c hanges1. Urbanization and extensive agricultural activities emerge as
leading water resources affecting anthropogenic factors2–5. Consequently, eutrophication became the ubiquitous
problem of this e ra5–7. Usually, land cover changes increase the source of nitrate while decreasing the sink, and
the nitrate exports are elevated in suburban and urban areas compared to the natural e nvironment2,8. Nutrient
sources and the exchange across the complex surface ad subsurface aquatic boundaries of urban and agricultural
areas need to be understood to establish pollution control or correcting approaches8,9. Furthermore, the impacts
of urbanization on groundwater systems must be explored and considered in land use planning to make future
urban areas s ustainable3.
Nitrate sources in the urban area include leaky sanitary sewers, industrial spills, leaching from landfills,
stormwater runoff from roads, fertilizer use on lawns, sewage disposal, septic tanks, and infiltration from the
polluted river8,10,11. The impact of the majority of these pollution sources on water resources is a ddressed11 but
the latter, infiltration from polluted river networks, is ignored. Urban runoff is a widely known non-point type
of pollution source in urbanized river basins8,11,12. Septic tank is a frequently cited high concentration nitrate
source in urban groundwater, especially in developing countries, often with high unsewered s lums11.
Addressing non-point sources of nutrients comprehensively in agricultural areas is an ongoing c oncern8.
And several integrated regional-scale surfaces and subsurface models were developed and applied for evaluating
nutrient loading in agricultural river basins, e.g.,13–16. However, concentration differences in aquifers underlying
agricultural and urban areas are negligible and sometimes higher in urban regions because of the high number,
density, and concentration of potential s ources11. Thus, to prevent the adverse effects of urbanization on water
resources, it is essential to assess the changes in groundwater flow and contaminant transport related to pollution
1
Department of Water Resources and River Research, Korea Institute of Civil Engineering and Building Technology,
Goyang 10223, Republic of Korea. 2Civil and Environmental Engineering Department, University of Science and
Technology, Daejeon 34113, Republic of Korea. *email: [email protected]
in urban a quifers17. In river basins with hydraulically connected aquifer and surface water systems, the nutri-
ent loading is complex and needs insight into surface water-groundwater interactions (SGI) to characterize the
pollution18–20. Understanding regional-scale SGI, including the controlling factors, is also vital for conjunctive use
of surface water and groundwater resources, assessment, and control of water contamination, and sustenance of
wetlands21–23. However, human activities, especially groundwater over-abstraction and spatiotemporal recharge
change associated with Land use/Land cover (LULC) changes, continuously alter the SGI distribution3,24. An
SGI study is data extensive and often requires a multidisciplinary approach which makes it challenging for many
data-scarce regions9,19,25,26.
These interlinked challenges of surface water and groundwater management in the changing environments
require concerted efforts of hydrological, hydrogeological, and environmental scientists. Process-based coupled
surface water-groundwater models appear to be plausible for simulating complex physical processes, such as
urbanized and highly irrigated r egions12,26,27. Moreover, integrated models provide a more accurate modeling
approach, compared to separate saturated and unsaturated zones focused models, for best water resources man-
agement decisions28. Studies on nitrate loading into urban rivers by considering river-aquifer connection usually
focus on small reach or experimental catchment, e.g.,29. This study focused on the simulation of urban expan-
sion and groundwater pumping on the SGI and nutrient loading using an integrated modeling approach. This
work incorporated non-point and point pollution sources with urban expansion scenarios and evaluated using
an integrated surface-subsurface modeling approach. Of the many potential urban pollutants, we focused on
nitrate because it is a common contaminant that can cause adverse health effects in infants and animals besides
eutrophication8,11,13,16 and data scarcity in the study region.
The study region is in Africa. In Sub-Saharan Africa (SSA), where over 40% of the population lives in an urban
area, wastewaters are often discharged untreated and increasingly become a leading source of eutrophication-
causing nutrients6,30. In SSA, while population and urbanization have been increasing rapidly since 1950, most
countries could treat only less than 30% of the total wastewater p roduced6,30. While regional-scale groundwater
resources have been recently quantified, hydrogeological evidence at high resolution in many African urban
centers is often lacking30,31. The fate and nutrient transport mechanisms in soils and aquifers and the effect of
contaminants carried out from urban surfaces on nearby lakes and rivers, draining the urban areas, remain
unexplored6.
Figure 1. Description of the study area, including Major River networks, gauging stations, and reservoirs
processed using QGIS3.10.4 (https://www.qgis.org/en/site/).
intermittent precipitation. In Kiremt, the region receives an average rainfall of about 759 mm and 200 mm in
the Belge season. A warm temperate climate characterizes the area with an average temperature of 15 °C and
relative humidity varying from 40 to 73%.
SWAT model inputs and setup. The SWAT model setup needs LULC, soil properties, Digital Elevation
Model (DEM), and daily climate data to simulate the fundamental hydrological processes. SWAT divides the
model area (river basin) into subbasins and hydrologic response units (HRU). HRU represents the combina-
tion of unique LULC, soil properties, and slope classes. In this study, 30 m resolution LULC (for the years 2000,
2010, and 2020) and DEM accessed from GlobeLand30 provided by the National Geomatics Center of China47
and U.S. Geological Survey “earthexplorer” websites, respectively, were used (Fig. 2a, d). The GlobeLand30
LULC classification has ten types. However, in the study region only eight LULC classes, namely, Forest, Culti-
vated land, Grassland, Shrubland, Wetland, Water Bodies, Bareland, and Artificial surfaces (Table 1). The soil
data were acquired from Harmonized World Soil Database v1.248. Eutric Vertisols cover a large part of the area
(Fig. 2b). Daily wind speed, temperature, relative humidity, precipitation, and solar radiation from 1980 to 2013
were collected from the National Metrological Agency of Ethiopia.
Figure 2. SWAT model inputs mapped using QGIS3.10.4 (https://www.qgis.org/en/site/): (a) Digital Elevation
lasses48, (c) Slope classes, and (d) Land
Model (DEM 30 m, USGS; https://earthexplorer.usgs.gov/), (b) Soil C
47
Use/Land cover (LULC) for the year 2000 .
Table 1. Land use/Land cover classification and the percentage coverage from the total area for 2000, 2010,
and 2020 in the study region.
Figure 3. Geology and hydrogeology of Akaki River Basin, obtained from the Ministry of Water Resources,
Irrigation, and Electricity (EMWRIE), the principal boundary conditions considered in the groundwater flow
modeling, and location of contaminant point sources, processed using QGIS3.10.4 (https://www.qgis.org/en/
site/).
The study area, Akaki River Basin, was subdivided into 37 subbasins and an HRU of about 1363. The topog-
raphy varies from 3239 m around Entoto Mountain to 1813 m above sea level near the outlet of the basin. About
77% of the area has a slope of 0–15% (Fig. 2c). The HRU was generated using the zero-threshold value of LULC,
soil, and slope option.
Groundwater flow model and hydrogeological descriptions. Topographically the study area favors
the emergence of springs in the hills and valleys. The elevation descends from north to south, and more of the
central, eastern, and southeastern parts are low areas covered with thick Quaternary deposits. Details of geology
and hydrogeology data are important for building a representative model. However, like any other part of the
country, data scarcity is the principal limiting factor for groundwater-related studies in particular. And therefore,
the model developed for this study relied on prior studies41,49,50 and low resolution or short-duration data col-
lected from the Ministry of Water Resources, Irrigation, and Electricity (EMWRIE) and Addis Ababa Water and
Sewerage Authority (AAWSA).
Miocene–Pleistocene volcanic successions dominate the geology in Akaki River Basin, including acidic and
intermediate lava flows, basaltic lava flows, and pyroclastic flows forming an interlayered sequence with qua-
ternary faults39. The region has highly variable and complex aquifers. The aquifers in southern Addis Ababa
are mainly young volcanic rocks of lava flows and tectonic fractures. The thickness and hydraulic conductivity
of the unconsolidated sediments govern the subsurface infiltration, overlying weathered and fractured porous
volcanic rocks with relatively high infiltration capacity, and favoring circulation and storage of subsurface water.
Alluvial, residual, and lacustrine clay deposits dominate and vary with topography and geomorphology. The
alluvial deposits, composed of clay with increasing thickness from north to south, are the predominant geologi-
cal types along the rivers. A thin layer of residual clay soil covers the ridges and hillsides in and around the city.
Thick residual clay soils cover high plains around the basin boundary. Black cotton lacustrine deposits cover the
lower part of the region. In the area around the Akaki well-field, thick alluvial, lacustrine clay deposits overlay
the underlying volcanic rocks interbedded between the volcanic rocks such as scoria, scoriaceous basalt, and
basalts. As conceptualized by Ayenew et al.49 , the regional groundwater flow is from north to south. More details
of the geology and hydrogeology of the basin are presented by several researchers, e.g.,41,49.
The central volcanic hills and mountains, not crossed by faults, bound the river basin; the hydraulic conduc-
tivity is highly variable, 0.01–550 m/day49. The primary geologic classes in the basin are basalts, scoria, rhyolites,
trachytes, ignimbrites, trachybasalts, and tuff of varying ages (Fig. 3a). The scoria and scoriaceous basalts with
primary porosity and permeability make up highly productive aquifers in the region. The highly weathered and
fractured basalts, ignimbrite, fractured tuff, and pyroclastics are also highly productive aquifers with secondary
porosity and permeability. Fractured basalts, sparsely spaced joints and vesicles, ignimbrites, and agglomerates
form moderately productive aquifers. The transmissivity is also highly variable, increasing from the highlands
to Akaki well-field41,49.
Even though the study area is data-scarce, the hydrogeology framework is described in several previous works,
e.g.,49,51. Previously, the groundwater flow modeling was performed by using two vertical layers and 400 × 400 m
cell dimensions49. However, the thickness of the top layer is thin compared to the bottom layer, it varies from “a
few meters to 53 m”49. In this study, it was found that the data which identify the thickness of the first layer was
insufficient to build a reasonable model. And thus, a simplified convertible single-layer model was built. The
mode area was discretized into 300 × 300 m horizontal resolution, which results in 18,976 active cells. The top
and bottom boundaries were interpolated from the well-completion reports. The most important groundwater
recharge-discharge areas (boundary conditions) are groundwater pumping, springs, reservoirs, river and canal
networks, and percolation from precipitation. Even though pumping wells are distributed all over the basin,
the primary groundwater pumping area is the Akaki well-field. In the area, pumping wells with varying depths
up to 604 m below the surface and a pumping rate of 13,193 m 3/day were identified and incorporated into the
model using the Well Package. Details of the reservoir operation are unavailable, and thus, to simplify the data
requirements, it was modeled as General Head Boundary. The major rivers were imported from the SWAT model
and modeled using the River Package. Canal or small rivers in a few places were modeled using Drain Package
(Fig. 3b). The initial groundwater head was interpolated from the static water level record.
Integrated SWAT‑MODFLOW model. Both SWAT and MODFLOW models are widely used and well
known among hydrologists and hydrogeologists. And therefore, here, the general equation of SWAT and MOD-
FLOW models and a brief description of the integrated SWAT-MODFLOW code are given (Eqs. 1 and 2). The
general equation of the SWAT model is as follows42:
t
SWt = SW0 + (PCP − Qs − ET − Perc − Qgw ) (1)
i=1
where SWo and SWt are the quantities of initial and final soil water (mm/day), respectively, t is the time (days),
PCP is the precipitation (mm/day), Qs is the surface runoff (mm/day), ET is the evapotranspiration (mm/day),
Perc is the percolation (mm/day), and Qgw is the flow from the aquifer (mm/day).
MODFLOW is physically-based algorithm that solves groundwater flow through porous earth material43,52.
The model is based on the following general partial differential equation 43:
∂ ∂h ∂ ∂h ∂ ∂h ∂h
Kxx + Kyy + Kzz ± W = Ss (2)
∂x ∂x ∂y ∂y ∂z ∂z ∂t
where Kxx, Kyy, and Kzz are principal components of the hydraulic conductivity tensor in the x, y, and z special
directions, W is the source or sink, Ss is the specific storage (1/L), h is the hydraulic head (L), and t is time.
Among the advantages of using the integrated SWAT-MODFLOW is water yield computation. The basic
water yield calculation in the SWAT model is based on the following equation 42:
WY = Qs + Qgw + Qlat − Tloss (3)
3 3 3
where WY is water yield ( L /T), Qs is surface runoff ( L /T), Qgw is groundwater flow to the river ( L /T), Qlat is
lateral soil flow ( L3/T), Tloss is water loss through the riverbed ( L3/T). Tloss, in Eq. 3, is computed with an assump-
tion that when a channel does not receive groundwater, possibly it loses water through the side and bed and joins
bank storage or the deep aquifer and approximated using the following relationships:
Tloss = Kch .TQ .Pch .Lch (4)
where Kch is the effective hydraulic conductivity of the channel alluvium (L/T), TQ is the flow travel time (T), Pch
is wetted perimeter (L), and Lch is channel length (L).
SWAT-MODFLOW uses the MODFLOW River Package that is formulated based on Darcy’s law to compute
water loss from the river and groundwater discharge to the river with the following mathematical principle:
hch − hgw
Qleak = Kbed (Lch .Pch ) (5)
M
where Qleak is volumetric water flux between the river and the aquifer ( L3/T); Kbed is hydraulic conductivity of the
riverbed (L/T); M is riverbed thickness (L); hch and h gw are river stage and groundwater head (L), respectively; Pch
is wetted perimeter (L); L ch is channel length (L). Equation 5 shows that if the groundwater head is higher than
the water level in the channel, Qleak will be negative, which means water flows from the aquifer to the channel,
and Qleak will be positive if the reverse is true. In SWAT-MODFLOW water yield computation is based on Qleak
(Eq. 5) instead of Tloss (Eq. 4).
The recent SWAT-MODFLOW a lgorithm13,27,44 was used to integrate the calibrated SWAT-2012 and Newton
formulated MODFLOW-2005 MODFLOW-NWT:45. The current version of SWAT-MODFLOW44 has been applied to
study several hydrological processes, including groundwater r echarge26, regional-scale S GI44,53,54, and the effect
of natural and anthropogenic factors on water balance55–60. However, the diversity of application regions and
study objectives is limited compared to the individual SWAT and MODFLOW models.
Reactive transport in 3‑dimensions (RT3D) model. Recently, Wei et al.13 applied SWAT-MODFLOW
for the assessment of non-point source contaminant transport in agricultural areas by integrating it with Reac-
tive Transport in 3-Dimensions RT3D;61. RT3D is a contaminant and solute transport model that can simulate
dispersion, advection, and chemical reactions in a saturated groundwater s ystem61. The algorithm has the func-
tionality of solving multi-species reactive transport. The following governing equation describes the fate and
transport of aqueous and solid-phase species in multi-dimension saturate porous m edia61:
Figure 4. Conceptual diagram showing surface and subsurface processes and modeling framework, adopted
and modified from previous studies13,16,44, with copyright permission.
∂Ck ∂ Ĉk ∂ ∂Ck ∂
φ + ρb = φDij − (φvi Ck ) + qs Csk + φrc , where k = 1, 2, . . . , m (6)
∂t ∂t ∂xi ∂xj ∂xi
In Eq. 6, m is the total number of species, Ck is the concentration of the k th species (M/L3), t is time (T), Dij
in the dynamic dispersion coefficient tensor ( L2/T), xi is the distance along the respective axis coordinate (L), v i
average seepage velocity (L/T), ϕ is the porosity of the material, qs is the volumetric flux of water per unit volume
of aquifer representing sources and sinks (1/T), Csk is the concentration of the kth species in the sources or sinks
(M/L3), ρb is bulk density of porous media (M/L3), C k is the solid-phase concentration of the k th species (M/M),
rc the rate of all reactions that occur in the aqueous phase (ML3/T).
RT3D uses the groundwater hydraulic head, cell-by-cell flow data, and water sources and sinks simulated
by the MODFLOW model to establish the groundwater flow field. The detailed linkage processes with SWAT-
MODFLOW are provided in 13,16,61. SWAT-MODFLOW-RT3D combines SWAT-MODFLOW code as a base and
RT3D as a sub-routine into a single executable a lgorithm13. The integrated model uses the SWAT ‘in-stream’ algo-
rithm to rout nitrate through the stream network. In SWAT-MODFLOW, the SWAT model computed recharge
is mapped to MODFLOW grid cells, and MODFLOW simulates the groundwater flux in each cell61. RT3D uses
MODFLOW simulated groundwater flux. RT3D also receives the nitrate concentration from the SWAT model
recharge water and river network and computes the change in each grid cell and river mass loading.
In the study area, the major point sources of nitrates are industrial discharges into the river and leaching
from waste disposal sites 39,62. Sources with high concentrations such as industrial discharge and leaching from
solid waste disposal sites 63, along Little Akaki River, were selected (Fig. 3b). In these sites as identified by sev-
eral researchers, e.g., Angello et al. 63, the nitrate concentration varies from 33 mg/l to 300 gm/L. The biggest
waste disposal site called Koshe was modeled as constant concertation (300 mg/L) point source. The denitrifi-
cation is specified using single-Monod expression so that the reaction rate depends on the nitrate presence 64.
Thus, the values of the first order-rate denitrification constant (1/d) and Monod half-saturation constant (M/
L3) are required for the model setup. However, these values are not easy to measure and are rarely available for
watershed-scale studies, usually assumed based on commonly used values from the literature. In this study area,
related studies are unavailable, and therefore, we have referred to worldwide literature, e.g., 15,16,64–66, and con-
sidered 0.05/day and 10 g/m3 for the first-order denitrification and Monod half-saturation constant. The surface
and subsurface hydrological processes, models, and scenario simulations were conceptualized by modifying the
schematic diagram from previous studies 13,16,44. Different colors are used to identify the processes and the cor-
responding simulating algorithm (Fig. 4).
Model calibration and validation. The SWAT model was calibrated and validated using sequential
uncertainty fitting version 2 (SUFI-2) in the SWAT-CUP67,68 user interface. Monthly river flow at two gauging
stations, Little Akaki and Big Akaki, obtained from EMWRIE was used. The calibration and validation perfor-
mance was evaluated using the Nash–Sutcliffe Efficiency coefficient (NSE) and minimizing the Percent Bias
(PBIAS), Root Mean Square Error (RMSE), and coefficient of determination (R2).
Table 2. Initial range, calibrated value, and definition of parameters used in the calibration and validation
processes. r__default SWAT parameter is multiplied by (1 + optimized value); v__ replaced the default value.
Table 3. Detail of groundwater pumping wells based on hydrogeology map of the study area: Well depth,
pumping rate, number of wells located in each hydrogeologic zone, and pumping rate change (%).
The groundwater model was also calibrated using the Parameter ESTimation (PEST) Package. However, the
automatic calibration technique was limited to the steady-state simulation because of the transient groundwater
level data scarcity. The important calibrated parameters were hydraulic conductivity, river conductance, and drain
conductance. The zonal and pilot point techniques69 were used to calibrate the hydraulic conductivity using a
range value of 0.012 to 375 m/day. The initial groundwater head is a key input for the simulation of integrated
surface water and groundwater flow models. In this study, the initial head was based on the calibrated steady-
state MODFLOW model simulation. Nevertheless, transient simulation parameters (specific storage, specific
yield) of the aquifer calibrated after integration using a trial-and-error technique based on the groundwater level
fluctuation, groundwater discharge (or comparing the measured and simulated river flow), and the overall water
balance of the study area. In addition to these soft calibration and validation practices previous studies in the
region, e.g.,41,49,70,71, were important sources of initial ranges of parameters and were also used as a benchmark
for general the model outputs analogy.
Automatic calibration functions are unavailable in the recent version of SWAT-MODFLOW. Individually
calibrated SWAT and MODFLOW for Akaki River Basin were integrated, and moderate manual calibration was
performed (Table 2). In SWAT-MODFLOW-RT3D mainly longitudinal dispersivity and initial concentration
of nitrate in groundwater were adjusted by trial and error. We have used longitudinal dispersivity of 120–180 m
and constant porosity of 0.3.
Simulation scenarios. Two important water and nutrient mass balance affecting factors, namely LULC
changes focused on the urban expansion and groundwater pumping for urban water supply, were tested using an
integrated modeling approach. The effects were evaluated based on the baseline simulation. Therefore, the simu-
lations were baseline, the effect of groundwater pumping (scenario1), urban expansion in the period of 2000–
2010 (scenario 2), urban expansion from 2000 to 2010 and groundwater pumping (scenario 3), urban expansion
from 2010 to 2020 (scenario 4), and urban expansion between 2010 and 2020 and groundwater pumping (sce-
nario 5). The LULC data is tabulated in “Groundwater flow model and hydrogeological descriptions” Sections
(Table 1). The current groundwater pumping rate was used as a reference and increased by 15–35% to evaluate
the effect (Table 3). Related scenarios were simulated in a previous study49 to assess the aquifer response using
MODFLOW. Under all scenarios formulated and simulated in this study, the regional-scale SGI, groundwater
recharge, and nutrient loading were the focus of the analysis.
Calibration Validation
Subbasin (gauging station) R2 NSE PBIAS (%) R2 NSE PBIAS (%)
Big Akaki 0.85 0.68 − 1.12 0.84 0.72 − 0.27
Little Akaki 0.89 0.82 − 0.73 0.86 0.72 − 9.14
Table 4. SWAT-MODFLOW model calibration and validation performance indicating statistical values.
Figure 5. Calibration and validation of river flow at Little Akaki (a) and Big Akaki (b) gauging stations and the
average precipitation in the study area.
Baseline simulation results. The baseline simulation results were analyzed after the calibration and vali-
dation of the model. The statistical model evaluation using river flow shows an NSE value of 0.68–0.82, R 2
0.84–0.89, and PBIAS − 0.27 to − 9.14 (Table 4).
The model performance based on the criteria proposed by Moriasi et al.72 ranges from “good” to “very good”.
The peaks and lows of the flow match well with the corresponding basin average precipitation. However, the
performance at low and peak flows was dissimilar in the calibration and validation period. The model overesti-
mated the low flows and underestimated the peak flows. For Little Akaki, the validation performance was less
than the calibration, while for the Big Akaki, the validation was better than the calibration (Fig. 5). The SWAT-
MODFLOW based river flow showed a slight improvement at low flow compared to the SWAT model. In Little
Akaki, the flow slightly increased from 1997 onward compared to the preceding years, and rapid urbanization
could be one reason for this.
The calibrated groundwater head varied from 1800 to 2653 m (Fig. 6a). Previous research o utputs49 present a
closely related calibrated groundwater head. The simulated and observed head agreed fairly with an R 2 value of
0.98 and RMSE of 7.1(Fig. 6d). The model underestimated the groundwater level in the altitude areas.
The average groundwater discharge to the river map (Fig. 6b) shows that in the study area, the groundwater
discharge to the river occurred almost in all sections of the river. However, the seepage from the river to the
Figure 6. Map showing the calibrated starting groundwater head (a), the spatial distribution of groundwater
discharge to the river (b), the initial nitrate concentration (c) in the Akaki Aquifers prepared using the
QSWATMOD2 plugin (https://swat.tamu.edu/software/swat-modflow/), and scatter plot of measured and
simulated groundwater head (d).
aquifer was in limited areas, and the magnitude was significantly high, comparatively. The initial nitrate concen-
tration in groundwater was also evaluated using a limited number of measurements obtained from the EMWRIE
(Fig. 6c). The measured nitrate concentration varies from zero to 33 mg/L. The initial simulated concentration
varied from zero to 0.9 mg/L. The spatial distribution of high and low measured and simulated concentration
matches well with an RMSE value of 6.5, even though there is a difference at peak values. While, in a large part of
the area, the measured and simulated data show good agreement, the main reason for the peak value of measured
and simulated concentration mismatched is that the points with higher concentration have no spatial relation
and are limited in number. The spatial distribution of simulated initial nitrate concentration represented the
groundwater quality in the study area.
To show the model performance in transient groundwater flow simulation (SWAT-MODFLOW), the
groundwater discharge to the river and baseflow, analyzed from observed river flow using automated separa-
tion techniques73, were compared (Fig. 7). The baseflow and groundwater discharge matched fairly, but SWAT-
MODFLOW underestimated the discharge in the rainy period and overestimated in dry periods.
The LULC change effect was substantial on average annual basin surface runoff (surq), recharge (rch), and
lateral flow (latq). The change in lateral flow from the baseline was quite significant (~ 72%). The river seep-
age and groundwater discharge were controlled by both groundwater pumping and LULC changes (Table 5).
Figure 7. Comparison of SWAT-MODFLOW simulated groundwater discharge and baseflow in Little Akaki
and Big Aaki Rivers.
Table 5. Principal annual basin water balance components, including precipitation (prec), surface runoff
(surq), lateral flow (latq), groundwater discharge to the river (gwq), surface water discharge to the aquifer
(swgw), and recharge (rch) of the study area under all simulation scenarios (mm).
Simulation results using the 2000 LULC data shows that river seepage increased under groundwater abstraction
but decreased in all other scenarios. The seepage decreased by more than 4%. The groundwater discharge declined
with the increasing groundwater pumping while increasing under the LULC change scenarios. This result is likely
to be related to the moderate recharge increment under the LULC change scenarios.
Surface water‑groundwater interactions. The effect of the groundwater pumping and urban expan-
sion on SGI in terms of groundwater recharge, river seepage, and water yield over the simulation period was
immediately apparent. After several trial-and-error simulations, the effect of groundwater pumping, which
increased by an average of the proposed increase rate from the baseline pumping (25%), was considered. The
total groundwater abstraction in the baseline scenario was 99,000 m3/day.
The graphical comparison shows that the effect of groundwater pumping on the monthly water yield was
insignificant (Fig. 8a). However, the average value revealed that groundwater pumping caused a decrease in water
yield by up to 3%, the LULC change by 9%, and the combined LULC and groundwater pumping by up to 15%.
The urban expansion caused to increase in the groundwater recharge in Kiremt and a decrease in the Bega and
Belge seasons (Fig. 8b–d). Principally in the Belge season, the recharge showed a significant difference compared
to the baseline simulation (Fig. 8d). The recharge under scenarios 2 and 4 had an insignificant difference in all
seasons. These results further support the effect of rapid urbanization that occurred from 2010 to 2020 compared
to the preceding decade.
Maximum river seepage and groundwater discharge occurred in Belge and Kiremt, respectively (Fig. 9). The
LULC change had a decreasing effect on river seepage (Fig. 9a); the effect was significant in the Belge season.
On the other hand, the groundwater discharge increased with LULC change (Fig. 9b, scenarios 2 and 4). The
groundwater discharge difference between the seasons was trivial.
Figure 10 shows the groundwater head differences from the baseline simulation at selected wells. The fluctua-
tion was significant around Akaki well-field. The spatiotemporal groundwater level variability under all scenarios
(difference from the baseline simulation) was random and complex. However, comparatively, the groundwater
level decreased under almost all simulation scenarios, while in Big Akaki, predominantly in the upstream area,
the level showed a slight increase. The central reason is the distribution of the groundwater pumping wells. The
results revealed that the response of the aquifers to the LULC and pumping from upstream to downstream is
Figure 8. Average monthly groundwater recharge and water yield (a) and seasonal groundwater recharge
fluctuation from 1990 to 2013 (b–d) in the Akaki River Basin.
Figure 9. Seasonal surface water discharge to the aquifer (swgw) and groundwater discharge to the river (gwq)
differences under the simulation scenarios from the baseline scenario (scenarios-bassline).
highly variable. The groundwater head at the upstream area was more sensitive to LULC changes than ground-
water abstraction.
The most important flow parameter, which affects groundwater level fluctuation and other water balance
components, is groundwater recharge. In addition to the results presented in the previous section, the spatial
variability of the groundwater recharge related to the LULC change was analyzed. In the region, the groundwater
recharge reached a maximum of 162 m3/day, and it varied largely with the LULC change (Fig. 11). Surprisingly, as
Figure 10. Groundwater level difference from the baseline scenario from 1990 to 2013: positive values show
that the groundwater level decreased compared to the baseline simulation.
urban expansion increased in the period from 2000 to 2010, the area, which gets moderate to high groundwater
recharge, mainly around the boundary of the city, enlarged (compare Fig. 11a, b). The principal reason could be
the slight increase in Forest and Shrubland coverage in that period. Compared to scenario 2, the average recharge
under scenario 4, particularly in the urban area, increased slightly (Fig. 11c).
Nitrate mass seeping into the aquifer. The effect of urbanization on the average spatiotemporal dis-
tribution of nitrate mass seepage from the river to the aquifer is presented in Figs. 12 and 13. The seasonal
distribution shows that scenarios 4 and 5 increased the nitrate seepage significantly (Fig. 12a). The seepage was
high in the Belge under all simulation scenarios comparatively. These results reasonably mirror the finding on
water seepage from the river to the aquifer. However, the peaks and lows trend varied from year to year over the
simulation period (Fig. 12b–d). The potential source for this random nutrient seepage over the years could be
the groundwater recharge change corresponding to the LULC changes.
The average spatial distribution of nitrate mass seepage, evaluated on the SWAT model subbasin scale,
revealed that the seepage was moderate in a large part of the region (600–10,000 kg/day). In most of the areas,
the average nitrate seepage increased up to 17–250% from the baseline attributed to the urban expansion. How-
ever, in scenario 1, the seepage declined by about 4%. This shows that the average nitrate seepage in the region
is more sensitive to LULC change than groundwater pumping. The groundwater pumping effect was significant
around Aba Samuel lake, also known as the Akaki well-field (Fig. 13a, b). But the effect had no distinct relation-
ship with the spatial distribution of nitrate seepage in all other scenarios (Fig. 13c–f). A maximum of 116,173 kg/
day mass of nitrate seeps into the aquifer (averaged over the simulation period). Most of the area covering high
seepage was in the Little Akaki, the most urbanized subbasin of the region (Fig. 13e, f).
Figure 11. Spatial distribution of average groundwater recharge under baseline (a), scenario 2 (b), and scenario
4 (c) in the study area from SWAT-MODFLOW44 output prepared using QSWATMOD2 (https://swat.tamu.edu/
software/swat-modflow/).
over 4%, while the groundwater discharge decreased by less than 2%. The water yield decreased up to 3% under
groundwater abstraction scenarios and up to 9% in LULC change scenarios. Significant changes were found in
the Belge season. The combined urban expansion and groundwater pumping compounded the water balance
and nitrate seepage changes. The relation between water and nitrate seepage from the river to the aquifer did
not exhibit any unique relationship. High nitrate loading was found attributed to LULC changes, and the effect
of groundwater abstraction was small. On average, nitrate seepage increased by 17–250% from the baseline in
most of the area.
This overarching modeling work highlighted the effect of urbanization usually accompanied by an expansion
of the built area, groundwater over-abstraction, and surface water pollution, principally in developing countries
with poor waste management on groundwater. Even though the focus was not on providing calibrated and vali-
dated model, this work relied on limited data, and the results must be interpreted with caution. Incorporating
detailed and longtime-recorded data, different contaminant types, and simulation scenarios are suggested for
future studies. Further investigation on aquifer recovery and protection strategies is needed for sustainable water
resource management in the region.
Figure 12. Temporal distribution of average nitrate (No3) seepage to aquifer: (a) comparison of average
seasonal No3 seepage in each simulation scenario, (b–d) Seasonal No3 seepage from 1990 to 2013 under all
simulation scenarios.
Figure 13. Spatial distribution of nitrate (NO3) seepage (Kg/day) from the river to the aquifer: (a) baseline, (b)
scenario 1, (c) scenario 2, (d) scenario 3, (e) Scenario 4, and (f) Scenario 5.
Data availability
The datasets that support the findings of this study are available from the corresponding author on reasonable
request.
References
1. Liu, J. et al. Water scarcity assessments in the past, present, and future. Earth’s Future 5, 545–559 (2017).
2. Hibbs, B. J. Groundwater in urban areas. J. Contempor. Water Res. Educ. 159, 1–4 (2016).
3. Sharp, J. M. The impacts of urbanization on groundwater systems and recharge. Aquamundi 1, 51–56 (2010).
4. Fletcher, T. D., Andrieu, H. & Hamel, P. Understanding, management and modelling of urban hydrology and its consequences for
receiving waters: A state of the art. Adv. Water Resour. 51, 261–279 (2013).
5. Yu, L. et al. Urban hydrogeology: Transport routes and mixing of water and solutes in a groundwater influenced urban lowland
catchment. Sci. Total Environ. 678, 288–300 (2019).
6. Nyenje, P. M., Foppen, J. W., Uhlenbrook, S., Kulabako, R. & Muwanga, A. Eutrophication and nutrient release in urban areas of
sub-Saharan Africa — A review. Sci. Total Environ. 408, 447–455 (2010).
7. Walsh, C. J. et al. The urban stream syndrome: Current knowledge and the search for a cure. J. North Am. Benthol. Soc. 24, 706–723
(2005).
8. Cadenasso, M. L. et al. Exchanges across land-water-scape boundaries in urban systems. Ann. N. Y. Acad. Sci. 1134, 213–232
(2008).
9. Almasri, M. N. & Kaluarachchi, J. J. Modeling nitrate contamination of groundwater in agricultural watersheds. J. Hydrol. 343,
211–229 (2007).
10. Lalehzari, R., Tabatabaei, S. H. & Kholghi, M. Simulation of nitrate transport and wastewater seepage in groundwater flow system.
Int. J. Environ. Sci. Technol. 10, 1367–1376 (2013).
11. Wakida, F. T. & Lerner, D. N. Non-agricultural sources of groundwater nitrate: a review and case study. Water Res. 39, 3–16 (2005).
12. Barron, O. V., Barr, A. D. & Donn, M. J. Effect of urbanisation on the water balance of a catchment with shallow groundwater. J.
Hydrol. 485, 162–176 (2013).
13. Wei, X., Bailey, R. T., Records, R. M., Wible, T. C. & Arabi, M. Comprehensive simulation of nitrate transport in coupled surface-
subsurface hydrologic systems using the linked SWAT-MODFLOW-RT3D model. Environ. Model. Softw. 122, 104242 (2019).
14. Yang, L., Zheng, C., Andrews, C. B. & Wang, C. Applying a regional transport modeling framework to manage nitrate contamina-
tion of groundwater. Groundwater 59, 292–307 (2021).
15. Bailey, R. T., Gates, T. K. & Ahmadi, M. Simulating reactive transport of selenium coupled with nitrogen in a regional-scale irrigated
groundwater system. J. Hydrol. 515, 29–46 (2014).
16. Wei, X. & Bailey, R. T. Evaluating nitrate and phosphorus remediation in intensively irrigated stream-aquifer systems using a
coupled flow and reactive transport model. J. Hydrol. 598, 126304 (2021).
17. Onodera, S. I. et al. Effects of intensive urbanization on the intrusion of shallow groundwater into deep groundwater: Examples
from Bangkok and Jakarta. Sci. Total Environ. 404, 401–410 (2008).
18. Hare, D. K., Helton, A. M., Johnson, Z. C., Lane, J. W. & Briggs, M. A. Continental-scale analysis of shallow and deep groundwater
contributions to streams. Nat. Commun. 12(1), 1–10 (2021).
19. Conant, B., Robinson, C. E., Hinton, M. J. & Russell, H. A. J. A framework for conceptualizing groundwater-surface water interac-
tions and identifying potential impacts on water quality, water quantity, and ecosystems. J. Hydrol. 574, 609–627 (2019).
20. Sophocleous, M. Interactions between groundwater and surface water: the state of the science. Hydrogeol. J. 10, 52–67 (2002).
21. Larned, S. T., Gooseff, M. N., Packman, A. I., Rugel, K. & Wondzell, S. M. Groundwater–surface-water interactions: Current
research directions. Freshw. Sci. 34, 92–98 (2015).
22. Flores, L., Bailey, R. T. & Kraeger-Rovey, C. Analyzing the effects of groundwater pumping on an urban stream-aquifer system.
JAWRA J. Am. Water Resour. Assoc. 56, 310–322 (2020).
23. Winter, T. C., Harvey, J. W., Franke, O. L., Alley, W. M. Ground water and surface water: A single resource. U.S. Geological Survey.
Circular 1139. https://pubs.er.usgs.gov/publication/cir1139. https://doi.org/10.3133/cir1139 (1998).
24. Lerner, D. N. Identifying and quantifying urban recharge: A review. Hydrogeol. J. 10, 143–152 (2002).
25. Barthel, R. & Banzhaf, S. Groundwater and surface water interaction at the regional-scale—A review with focus on regional
integrated models. Water Resour. Manag. 30, 1–32 (2016).
26. Yifru, B. A., Chung, I. M., Kim, M. G. & Chang, S. W. Assessment of groundwater recharge in agro-urban watersheds using inte-
grated SWAT-MODFLOW model. Sustainability (Switzerland) 12, 6593 (2020).
27. Aliyari, F. et al. Coupled SWAT-MODFLOW model for large-scale mixed agro-urban river basins. Environ. Model. Softw. 115,
200–210 (2019).
28. Sarma, R. & Singh, S. K. Simulating contaminant transport in unsaturated and saturated groundwater zones. Water Environ. Res.
93, 1496–1509 (2021).
29. Beltran, J., Lautz, L. K. & Slosson, J. R. The impact of stream-groundwater exchange on seasonal nitrate loads in an urban stream.
Hydrol. Processes 35, e14324 (2021).
30. Lapworth, D. J. et al. Urban groundwater quality in sub-Saharan Africa: Current status and implications for water security and
public health. Hydrogeol. J. 25(4), 1093–1116 (2017).
31. Birhanu, B. et al. Impact of natural and anthropogenic stresses on surface and groundwater supply sources of the upper Awash
sub-basin, Central Ethiopia. Front. Earth Sci. 9, 656726 (2021).
32. Kebede, S., Charles, K., Godfrey, S., MacDonald, A. & Taylor, R. G. Regional-scale interactions between groundwater and surface
water under changing aridity: Evidence from the River Awash Basin, Ethiopia. Hydrol. Sci. J. 66, 450–463 (2021).
33. Ayenew, T., Demlie, M. & Wohnlich, S. Hydrogeological framework and occurrence of groundwater in the Ethiopian aquifers. J.
Afr. Earth Sc. 52, 97–113 (2008).
34. Aregahegn, Z. & Zerihun, M. Study on irrigation water quality in the rift valley areas of Awash river basin, Ethiopia. Appl. Environ.
Soil Sci. 2021, 1–20 (2021).
35. Mersha, A., Masih, I., de Fraiture, C., Wenninger, J. & Alamirew, T. Evaluating the impacts of IWRM policy actions on demand
satisfaction and downstream water availability in the upper Awash basin, Ethiopia. Water (Basel) 10, 892 (2018).
36. Yard, E. et al. Metals exposures of residents living near the Akaki River in Addis Ababa, ethiopia: A cross-sectional study. J. Environ.
Public Health 2015, 1–8 (2015).
37. Demlie, M. & Wohnlich, S. Soil and groundwater pollution of an urban catchment by trace metals: case study of the Addis Ababa
region, central Ethiopia. Environ. Geol. 51, 421–431 (2006).
38. Abiye, T. A., Sulieman, H., Ayalew, M. & Alemayehu, T. Use of treated wastewater for managed aquifer recharge in highly populated
urban centers: A case study in Addis Ababa, Ethiopia. Environ. Geol. 58, 55–59 (2009).
39. Yilma, M., Kiflie, Z., Windsperger, A. & Gessese, N. Application of artificial neural network in water quality index prediction: A
case study in Little Akaki river, Addis Ababa, Ethiopia. Model. Earth Syst. Environ. 4(1), 175–187 (2018).
40. Mengistu, T. D. et al. Challenges and prospects of advancing groundwater research in Ethiopian aquifers: A review. Sustainability
13, 11500 (2021).
41. Demlie, M., Wohnlich, S., Wisotzky, F. & Gizaw, B. Groundwater recharge, flow and hydrogeochemical evolution in a complex
volcanic aquifer system, central Ethiopia. Hydrogeol. J. 15, 1169–1181 (2007).
42. Arnold, J. G., Srinivasan, R., Muttiah, R. S. & Williams, J. R. Large area hydrologic modeling and assessment part I: Model develop-
ment. J Am Water Resour Assoc 34, 73–89 (1998).
43. Harbaugh, A. W. MODFLOW-2005, The U.S. Geological Survey modular ground-water model — the ground-water flow process.
Modeling techniques. US Department of the Interior, US Geological Survey, (2005).
44. Bailey, R. T., Wible, T. C., Arabi, M., Records, R. M. & Ditty, J. Assessing regional-scale spatio-temporal patterns of groundwater-
surface water interactions using a coupled SWAT-MODFLOW model. Hydrol. Process. 30, 4420–4433 (2016).
45. Niswonger, R. G., Panday, S. & Motomu, I. MODFLOW-NWT , A Newton formulation for MODFLOW-2005. In U.S. Geological
Survey Techniques and Methods 44 (2011).
46. Akoko, G., Le, T. H., Gomi, T. & Kato, T. A review of swat model application in africa. Water (Switzerland) 13, 1313 (2021).
47. Jun, C., Ban, Y. & Li, S. Open access to Earth land-cover map. Nature 514, 434–434 (2014).
48. Fao/Iiasa/Isric/Isscas/Jrc. Harmonized world soil database (version 1.2). (2012)
49. Ayenew, T., Demlie, M. & Wohnlich, S. Application of numerical modeling for groundwater flow system analysis in the Akaki
catchment, Central Ethiopia. Math. Geosci. 40, 887–906 (2008).
50. Muleta, D. & Abate, B. Groundwater hydrodynamics and sustainability of Addis Ababa city aquifer. Groundw. Sustain. Dev. 12,
100485 (2021).
51. Yitbarek, A., Razack, M., Ayenew, T., Zemedagegnehu, E. & Azagegn, T. Hydrogeological and hydrochemical framework of Upper
Awash River basin, Ethiopia: With special emphasis on inter-basins groundwater transfer between Blue Nile and Awash rivers. J.
Afr. Earth Sci. 65, 46–60 (2012).
52. McDonald, M. G. & Harbaugh, A. W. A Modular Three-Dimensional Finite-Difference Ground-Water Flow Model. https://pubs.er.
usgs.gov/publication/ofr83875. https://doi.org/10.3133/ofr83875(1988).
53. Taie Semiromi, M. & Koch, M. Analysis of spatio-temporal variability of surface–groundwater interactions in the Gharehsoo river
basin, Iran, using a coupled SWAT-MODFLOW model. Environ. Earth Sci. 78, 201 (2019).
54. Jafari, T., Kiem, A. S., Javadi, S., Nakamura, T. & Nishida, K. Fully integrated numerical simulation of surface water-groundwater
interactions using SWAT-MODFLOW with an improved calibration tool. J. Hydrol. Reg. Stud. 35, 100822 (2021).
55. Chunn, D., Faramarzi, M., Smerdon, B. & Alessi, D. Application of an integrated SWAT–MODFLOW model to evaluate potential
impacts of climate change and water withdrawals on groundwater-surface water interactions in West-Central Alberta. Water
(Basel) 11, 110 (2019).
56. Yifru, B. A., Chung, I.-M., Kim, M.-G. & Chang, S. W. Assessing the effect of land/use land cover and climate change on water
yield and groundwater recharge in east african rift valley using integrated model. J. Hydrol. Reg. Stud. 37, 100926 (2021).
57. Yifru, B. A., Chung, I.-M., Kim, M.-G. & Chang, S. W. Augmenting freshwater availability in mountain headwater streams: Assess-
ing the sustainability under baseline and future climate change scenarios. Int. Soil Water Conserv. Res. 10, 293–307 https://doi.org/
10.1016/j.iswcr.2021.08.002 (2021).
58. Guevara-Ochoa, C., Medina-Sierra, A. & Vives, L. Spatio-temporal effect of climate change on water balance and interactions
between groundwater and surface water in plains. Sci. Total Environ. 722, 137886 (2020).
59. Wei, X. & Bailey, R. T. Assessment of system responses in intensively irrigated stream-aquifer systems using SWAT-MODFLOW.
Water (Basel) 11, 1576 (2019).
60. Liu, W. et al. Assessing the impacts of groundwater abstractions on flow regime and stream biota: Combining SWAT-MODFLOW
with flow-biota empirical models. Sci. Total Environ. 706, 135702 (2020).
61. TP Clement. A Modular Computer Code for Simulating Reactive Multi-Species Transport in 3-Dimensional Groundwater Systems.
http://www.osti.gov/servlets/purl/8022-44Hkz5/webviewable/. https://doi.org/10.2172/8022 (1999).
62. Melaku, S., Wondimu, T., Dams, R. & Moens, L. Pollution status of Tinishu Akaki River and its tributaries (Ethiopia) evaluated
using physico-chemical parameters, major ions, and nutrients. Bull. Chem. Soc. Ethiop. 21, 13–22 (2007).
63. Angello, Z. A., Behailu, B. M. & Tränckner, J. Integral application of chemical mass balance and watershed model to estimate point
and nonpoint source pollutant loads in data-scarce Little Akaki river, Ethiopia. Sustainability 12, 7084 (2020).
64. Wei, X., Bailey, R. T., Records, R. M., Wible, T. C. & Arabi, M. Comprehensive simulation of nitrate transport in coupled surface-
subsurface hydrologic systems using the linked SWAT-MODFLOW-RT3D model. Environ. Model. Softw. 122, 104242 (2019).
65. Bailey, R. T., Gates, T. K. & Romero, E. C. Assessing the effectiveness of land and water management practices on nonpoint source
nitrate levels in an alluvial stream–aquifer system. J. Contam. Hydrol. 179, 102–115 (2015).
66. Sarma, R. & Singh, S. K. Simulating contaminant transport in unsaturated and saturated groundwater zones. Water Environ. Res.
93, 1496–1509 (2021).
67. Abbaspour, K. C. SWAT-CUP4: SWAT calibration and uncertainty programs-a user manual. Department of Systems Analysis,
Integrated Assessment and Modelling(SIAM), Eawag, Swiss Federal Institute of Aquatic Science and Technology, Duebendorf,
Switzerland. https://swat.tamu.edu/media/114860/usermanual_swatcup.pdf (2015)
68. Abbaspour, K. C. et al. Modelling hydrology and water quality in the pre-alpine/alpine Thur watershed using SWAT. J. Hydrol.
333, 413–430 (2007).
69. Doherty, J. Ground Water Model Calibration Using Pilot Points and Regularization. Ground Water 41, 170–177 (2003).
70. Demlie, M., Wohnlich, S., Gizaw, B. & Stichler, W. Groundwater recharge in the Akaki catchment, central Ethiopia: Evidence from
environmental isotopes (δ18O, δ2H and3H) and chloride mass balance. Hydrol. Process. 21, 807–818 (2007).
71. Demlie, M. Assessment and estimation of groundwater recharge for a catchment located in highland tropical climate in cen-
tral Ethiopia using catchment soil–water balance (SWB) and chloride mass balance (CMB) techniques. Environ. Earth Sci. 74,
1137–1150 (2015).
72. Moriasi, D. N. et al. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans ASABE
50, 885–900 (2007).
73. Arnold, J. G. & Allen, P. M. Automated methods for estimating baseflow and ground water recharge from streamflow records. J.
Am. Water Resour. Assoc. 35, 411–424 (1999).
Author contributions
B.A.Y.: Conceptualization, Methodology, Software, Original draft preparation, Reviewing and Editing. I.M.C.:
Conceptualization, Methodology, Supervision, Reviewing and Editing. M.G.K.: Data curation, Software, Visu-
alization. S.W.C.: Software, Supervision.
Funding
This research was supported by a grant from a Strategic Research Project (20220178-001) funded by the Korea
Institute of Civil Engineering and Building Technology.
Competing interests
The authors declare no competing interests.
Additional information
Correspondence and requests for materials should be addressed to I.-M.C.
Reprints and permissions information is available at www.nature.com/reprints.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
1. use such content for the purpose of providing other users with access on a regular or large scale basis or as a means to circumvent access
control;
2. use such content where to do so would be considered a criminal or statutory offence in any jurisdiction, or gives rise to civil liability, or is
otherwise unlawful;
3. falsely or misleadingly imply or suggest endorsement, approval , sponsorship, or association unless explicitly agreed to by Springer Nature in
writing;
4. use bots or other automated methods to access the content or redirect messages
5. override any security feature or exclusionary protocol; or
6. share the content in order to create substitute for Springer Nature products or services or a systematic database of Springer Nature journal
content.
In line with the restriction against commercial use, Springer Nature does not permit the creation of a product or service that creates revenue,
royalties, rent or income from our content or its inclusion as part of a paid for service or for other commercial gain. Springer Nature journal
content cannot be used for inter-library loans and librarians may not upload Springer Nature journal content on a large scale into their, or any
other, institutional repository.
These terms of use are reviewed regularly and may be amended at any time. Springer Nature is not obligated to publish any information or
content on this website and may remove it or features or functionality at our sole discretion, at any time with or without notice. Springer Nature
may revoke this licence to you at any time and remove access to any copies of the Springer Nature journal content which have been saved.
To the fullest extent permitted by law, Springer Nature makes no warranties, representations or guarantees to Users, either express or implied
with respect to the Springer nature journal content and all parties disclaim and waive any implied warranties or warranties imposed by law,
including merchantability or fitness for any particular purpose.
Please note that these rights do not automatically extend to content, data or other material published by Springer Nature that may be licensed
from third parties.
If you would like to use or distribute our Springer Nature journal content to a wider audience or on a regular basis or in any other manner not
expressly permitted by these Terms, please contact Springer Nature at