A Comparison of Streamﬂow and Baseﬂow Responses to Land-Use Change and the Variation in Climate Parameters Using SWAT

: Alteration of land use and climate change are among the main variables a ﬀ ecting watershed hydrology. Characterizing the impacts of climate variation and land use alteration on water resources is essential in managing watersheds. Thus, in this research, streamﬂow and baseﬂow responses to climate and land use variation were modeled in two watersheds, the Upper West Branch DuPage River (UWBDR) watershed in Illinois and Walzem Creek watershed in Texas. The variations in streamﬂow and baseﬂow were evaluated using the Soil and Water Assessment Tool (SWAT) hydrological model. The alteration in land use between 1992 and 2011 was evaluated using transition matrix analysis. The non-parametric Mann–Kendall test was adopted to investigate changes in meteorological data for 1980–2017. Our results indicate that the baseﬂow accounted for almost 55.3% and 33.3% of the annual streamﬂow in the UWBDR and Walzem Creek watersheds, respectively. The contribution of both land use alteration and climate variability on the ﬂow variation is higher in the UWBDR watershed. In Walzem Creek, the alteration in streamﬂow and baseﬂow appears to be driven by the e ﬀ ect of urbanization more than that of climate variability. The results reported herein are compared with results reported in recent work by the authors in order to provide necessary information for water resources management planning, as well as soil and water conservation, and to broaden the current understanding of hydrological components variation in di ﬀ erent climate regions.


Introduction
Ecosystems and humans are fundamentally dependent on different water resources. Thus, for the general development of any country, the quality and the quantity of these water resources flowing through rivers are of vital importance to socio-economic development [1]. Issues related to changes in water resources are commonly evaluated around the globe [2][3][4]. In the United States, evaluation of streamflow and baseflow has been documented [5][6][7]. However, the quantitative change in streamflow and baseflow has yet to be evaluated across different climatic conditions. Climate alterations and human actions both act as stressors to place severe pressure on water resources [8,9]. The variations in climate and land use directly impact total streamflow, interflow, surface runoff, and baseflow, causing events of droughts and floods that impact the sustainability of these resources and the social ecosystem [10]. Several studies have examined alterations in streamflow due to changes in temperature and precipitation [11][12][13], urbanization [14], and land use change [2,15]. Baseflow is the portion of streamflow sustained in a river by delayed pathways. Baseflow is often assumed to be equal to groundwater recharge [16]. It provides a relatively high water quality with a

Walzem Creek
Walzem Creek is located in Bexar County in the state of Texas and in the San Antonio East USGS quad ( Figure 1b). Currently, except for the lower most portion of the watershed, the majority of Water 2020, 12, 191 5 of 29 Walzem Creek is characterized by dense, urban development. The lower portion of the watershed is characterized by a mix of vegetation and forests and normally dry except during rain events. The Walzem Creek is a part of Upper San Antonio Watershed and covers approximately 109 km 2 (USGS Gauge 08178800), with a mean annual precipitation ranging from 320 to 1200 mm and average annual temperatures ranging from 19.3 to 22.3 • C. Mean elevation in the area is 204 m, with a minimum and maximum of 152 and 266 m above sea level, respectively. Similar to the UWBDR watershed, recently, most of Walzem Creek is covered with developed areas (84.5%); however, other land uses include wetlands (7.2%), shrublands (4.2%), and forests (2.5%). However, it was characterized by only 64.4% of residential areas, 17.5% of planted cover, and 8.7% and 7.0% of forest and shrublands covers, respectively, at the beginning of the 1990s. This area is a large portion of the Upper San Antonio Watershed; hence, it contributes a large amount of total streamflow. According to the main Koppen-Geiger climate classes for US counties, the San Antonio, Bexar County area lies at the border between warm, humid, equatorial zone and fully hot arid and steppe zone [41]. Therefore, this watershed is representative of semi-arid regions.

Materials and Methods
The data described herein include spatial topography, Digital Elevation Model (DEM), land use and soil data, and hydro-meteorological data. Data analysis procedures and methods used are detailed extensively in the work by Aboelnour et al. [8]. A flow chart depicting procedures used in this study is summarized in Figure 2.
Water 2020, 12, x FOR PEER REVIEW 6 of 30

Materials and Methods
The data described herein include spatial topography, Digital Elevation Model (DEM), land use and soil data, and hydro-meteorological data. Data analysis procedures and methods used are detailed extensively in the work by Aboelnour et al. [8]. A flow chart depicting procedures used in this study is summarized in Figure 2.

Spatial Data
Two raster land use maps for the years 1992 and 2011 were obtained from the National Map Viewer (NMV). Digital Elevation Model (DEM) data were acquired from the Geospatial Data Gateway (GDG) with a resolution of 10 m. Soil Survey Geographic Data (SSURGO) data were used in this research with a resolution ranging from 1:12,000 to 1:63,630. Land use, soil type, and slope were then used to divide the delineated sub-basins into a small series of uniform HRUs that represent the smallest representative units within the watershed [42].

Hydro-Meteorological Data
The required datasets used include daily observed streamflow data at gauged USGS stations for the period 1980 to 2017. The streamflow data were used to separate baseflow from surface runoff,

Spatial Data
Two raster land use maps for the years 1992 and 2011 were obtained from the National Map Viewer (NMV). Digital Elevation Model (DEM) data were acquired from the Geospatial Data Gateway (GDG) with a resolution of 10 m. Soil Survey Geographic Data (SSURGO) data were used in this research with a resolution ranging from 1:12,000 to 1:63,630. Land use, soil type, and slope were then Water 2020, 12, 191 6 of 29 used to divide the delineated sub-basins into a small series of uniform HRUs that represent the smallest representative units within the watershed [42].

Hydro-Meteorological Data
The required datasets used include daily observed streamflow data at gauged USGS stations for the period 1980 to 2017. The streamflow data were used to separate baseflow from surface runoff, and for the SWAT model calibration and validation. In addition, long-term daily meteorological datasets for the same period  were collected from the National Climate Data Center (NCDC). The meteorological weather stations were 12 km and 0.8 km away from the borders of the UWBDR and Walzem Creek watersheds, respectively.

Baseflow Separation
Baseflow measurements were separated from daily streamflow data acquired from USGS gauged stations using the automatic baseflow digital filter method (BFlow). The BFlow filter separates streamflow data into baseflow and surface runoff by passing the observed streamflow through the filtering equation three times [43,44]: where BF is the baseflow, α is the filter parameter (0.925), Q is the total streamflow, and t is the time step. Equation (1) is applied only when BF ≤ Q t [45]. BFlow is a conservative filter that enables the user to filter streamflow data to calculate the baseflow, and also to generate a tabular dataset or graphical hydrograph interface from USGS gaging stations. Herein, BFlow filtered streamflow data three times (Equation (1)), and it is commonly observed that the 1-pass baseflow is consistent with manually estimated baseflow and thus was subsequently used in this study [45].

Soil and Water Assessment Tool (SWAT) Model Calibration and Validation
The monotonic trends in the historical meteorological data were evaluated using the modified Mann-Kendall (MK) test developed by Hamed and Rao [46]. Based on the abrupt change in the trends in precipitation and temperature using the MK test, the study period from 1980 to 2017 was split into two time spans, 1980-1998 and 1999-2017, with a breakpoint in 1998. The period 1980-1998 was assigned as a baseline for model calibration and validation. The model simulation time was segmented into a warm up period (1980)(1981)(1982)(1983), calibration period (1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993), and validation period (1994)(1995)(1996)(1997)(1998)). The SWAT model calibration and validation were performed using the land use map of 1992 and streamflow data from 1980 to 1998 for each of the selected watersheds. Model optimization, sensitivity analysis, calibration, validation, and uncertainty analysis of parameters were carried out using the Sequential Uncertainty Fitting program algorithm (SUFI-2) approach within the SWAT-CUP interface developed by Abbaspour et al. [47]. Based on Aboelnour et al. [8], the twenty hydrologic parameters listed in Table 1 were used in this study for the UWBDR and Walzem Creek watersheds calibration of streamflow and baseflow. However, sensitivity analysis using the SUFI-2 global sensitivity analysis was carried out in the first stage due to the presence of many parameters within the SWAT model [44]. Only parameters sensitive for the watersheds were then used in the calibration process and optimized based on monthly values [48]. Both automatic and manual calibration were carried out to allow qualitative and quantitative comparisons of the values, to fine tune the values of the auto-calibrated parameters, and to decrease the differences between the observed and simulated outputs [49].

Model Sensitivity Analysis
The global sensitivity analysis procedures showed that the sensitive parameters obtained from the LEC by Aboelnour et al. [8] were critical in the case of the UWBDR watershed, but with a different rank order. It was also found that these rankings were impacted by the selected objective function in the model. For example, curve number (CN2), soil evaporation compensation factor (ESCO), snowfall temperature (SFTMP), melt factor for snow (SMFMN), baseflow recession constant (ALPHA_BF), and deep aquifer percolation fraction (RCHRG_DP) were the most critical parameters in UWBDR when the Kling-Gupta Efficiency (KGE) was selected to be the objective function incorporated into the model (Table 2). These parameters characterize surface runoff, soil properties, and groundwater. In the case of Walzem Creek, the parameters in Table 2 were consistent with other SWAT parameter sensitivity analyses completed for semi-arid regions. The SWAT model is highly sensitive to surface runoff and basin parameters when the watershed is characterized by inconsistent rainfall events [50,51]. ALPHA_BF followed by CN2 were the most sensitive parameters in Walzem Creek. In contrast to the other watersheds, snowfall and snow melt parameters were not sensitive in Walzem Creek since there was no persistent snowpack. The high ALPHA_BF constant in Walzem Creek indicated a rapid response to groundwater recharge. However, the lower baseflow recession constant in the UWBDR indicated large storage discharge and slow drainage in the shallow aquifer, which might be attributed to the complex geological structure of the watershed such as the presence of folds and faults [36]. The high deep aquifer percolation parameter (RCHRG_DP) in Walzem Creek indicated the increase of water movement to the deep aquifer. SOL_AWC represented the soil moisture content and hence played a role in surface runoff and was considered to be directly proportional to the soil's ability to hold water, affecting streamflow.

Statistical Criteria and Model Evaluation Performance
The performance of the SWAT model can be computed using statistical indices and graphical comparisons [52]. For the simulated streamflow and baseflow, the coefficient of determination (R 2 ), Nash-Sutcliffe model efficiency (ENS), PBIAS, and modified KGE were adopted to evaluate the model performance [53,54]. The monthly statistical streamflow and baseflow values for the calibrated models were adopted to evaluate the model performance. The performance of the SWAT model is considered good on a monthly basis when R 2 > 0.75; ENS and KGE > 0.7; and PBIAS ≤ 15 according to Moriasi et al. [55] and Thirel et al. [56].

Scenarios Separating the Impact of Land Use Change and Climate Change
In this research, the "change-fix" approach used by Aboelnour et al. [8] was applied to evaluate the streamflow and baseflow as a response to separate and combined impacts of urbanization and climate alteration. Land use maps of 1992 and 2011 were used to represent the two time periods. The land use map of 1992 was adopted to represent the patterns in the first period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998), herein called TS1. On the other hand, the 2011 land use map was used to represent the second time span (1999-2017), herein called TS2.
A combination of four simulations were developed to evaluate the natural and human impacts on hydrology: (1) 1992 land use and TS1 climate data of 1980-1998 (X1); (2) 2011 land use and TS1 climate data of 1980-1998 (X2); (3) 1992 land use and TS2 climate data of 1999-2017 (X3); and (4) 2011 land use and TS2 climate data of 1999-2017 (X4). The well-calibrated SWAT model, using the land use data of 1992 and first climate period, was used to run the other four scenarios (X1-X4). The simulated output values obtained from these scenarios were compared to the corresponding baseline model. X1 represents the baseline scenario with the corresponding circumstances, while the difference between X4 and X1 simulation describes the combined effects of land use change and climate variation. The comparison between X1 and X2 attempts to depict the separate impact of land use change. Finally, the differences between X3 and X1 outputs emphasize the individual impact of climate alteration.

Trends in Hydrologic Components
Statistical analyses were performed on climatological variables using the modified non-parametric Mann-Kendall (MK) test [46] to evaluate possible transition points, trends, and their significance in the time series from 1980 to 2017. The modified MK test statistic, τ, is standardized and can be used in comparing variables that experience differences in their magnitude [57]. A positive slope magnitude indicates an upward trend and vice versa [58]. As shown in Table 3, the slope and the τ -statistics for annual streamflow and baseflow were all positive, except for the baseflow trend in Walzem Creek, Water 2020, 12, 191 9 of 29 which showed a significant decrease in monotonic trend. However, the null hypothesis was accepted in the case of annual streamflow, as it showed an insignificant increasing trend ( Figure 3). Results also show that the annual baseflow increased at a significance level greater than 0.1 for the UWBDR watershed, which indicates a slightly increasing trend. However, a significant increasing trend in streamflow during 1980-2017 was detected for the UWBDR watershed ( Figure 3).
In the case of UWBDR, the trend of τ-test statistics and the slope of precipitation and temperature were positive and are provided in Table 4. The results show a difference in the monotonic trends of annual temperature and precipitation. For the time series from 1980 to 2017, the annual air temperature increased at a significance level greater than 0.001, which indicates that the long-term trend of temperature is statistically significant. The annual precipitation increased only at a significance level greater than 0.1, indicating a minor increase of precipitation over time and that the trend is statistically insignificant at the 95% confidence level. On the other hand, the average annual precipitation after the change point in Walzem Creek exhibited a slight decrease from the average before the change point. However, the temperature at Walzem Creek showed an increasing trend at the 0.001 significance level, which indicates that the climate at Walzem Creek became warmer and The increase in average precipitation played an important role in the increasing trend of streamflow for the UWBDR watershed, while the slight increase in streamflow at Walzem Creek was accompanied by decreased precipitation and an increase in temperature as well. Moreover, human activity, such as construction of urban areas on agricultural areas, played a vital role in the amount of streamflow and baseflow. A combination of temperature increase and either a reduction or increase in rainfall are likely the main reasons for climate variation affecting the global water balance. In other words, the magnitude and the directions of these changes will affect any particular change in streamflow and baseflow [20].
The relationship between baseflow and human impacts and climate change varied. The reduction of annual baseflow in the Walzem Creek watershed may be attributed to the reduction of cultivated area and implementation of imperviousness, which in turn has a negative impact on the infiltration rate by increasing the surface runoff, specifically during the wet season of the year [49]. On the other hand, the increasing trend in the annual baseflow in the UWBDR was similar to the trend of the Little Eagle Creek (LEC) watershed mentioned by Aboelnour et al. [8]. This increase is likely caused by several factors, including the influx of water from outside the watershed during the process of urban development and infrastructure and leakage from water supply pipes. Lerner (2002) reported that urbanized catchments are usually associated with leakage rates of 20-50% in sewer systems and septic tanks, causing large amounts of groundwater discharge [59]. Wastewater from the West Chicago Moraine may also provide a significant amount of water, which likely originates outside the catchment. In addition, detention basins play essential roles in increasing baseflow in urban catchments, as water is retained at the surface due to an increasing amount of surface runoff, and then slowly released into the stream as a form of baseflow [8]. Lastly, physiological features may also contribute to this observed trend, including features such as the topography, geology and soil types that result from glacial melting with high porous media, which can play a significant role in increasing infiltration [60].

Trends in Climatic Components
The MK test was furthermore employed to quantify the monotonic trends of precipitation and temperature in the selected watersheds. Compared to the first climate period (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998), statistical results indicated that the mean air temperature increased by 0.7 • C (from 9.7 to 10.4 • C) and 0.6 • C (from 20.7 to 21.3 • C) during TS2 at the UWBDR and Walzem Creek watersheds, respectively. Average annual precipitation increased by 9.1% (82 mm, from 890 to 972 mm) during TS2 in the UWBDR, while decreasing by 6.5% (56 mm, from 858 to 802 mm) in Walzem Creek ( Figure 4).    In the case of UWBDR, the trend of τ-test statistics and the slope of precipitation and temperature were positive and are provided in Table 4. The results show a difference in the monotonic trends of annual temperature and precipitation. For the time series from 1980 to 2017, the annual air temperature increased at a significance level greater than 0.001, which indicates that the long-term trend of temperature is statistically significant. The annual precipitation increased only at a significance level greater than 0.1, indicating a minor increase of precipitation over time and that the trend is statistically insignificant at the 95% confidence level. On the other hand, the average annual precipitation after the change point in Walzem Creek exhibited a slight decrease from the average before the change point. However, the temperature at Walzem Creek showed an increasing trend at the 0.001 significance level, which indicates that the climate at Walzem Creek became warmer and drier during the study period. While the average annual precipitation and temperatures shifted over time, these trends may not reflect the true picture as the change displayed in both may have been seasonally influenced [61]. Therefore, the MK test was further performed at a monthly scale for time series data from 1980 to 2017 ( Table 5). The results show that the monotonic trends of the monthly meteorological data for the study were different. For the UWBDR watershed, the results indicate that the monthly temperature showed increasing trends in slope in every month of the year. The monotonic increasing trends of monthly temperature were only statistically significant at a confidence level of p = 0.05 in June, September, and October. Monthly precipitation trend for November decreased significantly, while it showed an insignificant reduction in August, September, and December. The remaining months showed an insignificant increase in monthly precipitation, with the highest increment recorded in June (1.36 mm/month). On the other hand, monthly precipitation in Walzem Creek Watershed showed decreasing trends in February, May, June, October, November, and December, while increasing trends in the other months with the highest increment recorded in September (1.52 mm/month) and maximum reduction recorded in June (1.68 mm/month). Change in monthly precipitation in the Walzem Creek Watershed was insignificant. Similar to the UWBDR watershed, increases in monthly temperature trends were recorded in every month, with only February showing a significant increase in the monotonic trend at a confidence level of p = 0.05.

Changes in Land Use Characteristics
Cross tabulation analysis and post classification comparison were applied to evaluate the quantity of temporal conversions and nature of changes from one land cover category to another in land use maps of 1992 and 2011 [61][62][63]. In the UWBDR, a comparison of land use maps for the years 1992 and 2011 indicated that the most significant changes occurred in three classes: developed urban, planted, and forest ( Figure 5). In 1992, the main land use types were planted and developed areas, which occupied 76.1% of the total watershed area. However, owing to urban expansion, the proportional extent of developed areas increased from 44% to 77% from 1992 to 2011. Conversely, the proportional extent of planted and forest decreased from 35.8% to 2.7% and from 8.1% to 4.4%, respectively. The transition matrix of UWBDR land use in Table 6 explains these changes in detail. Overall, 43.6% or 39.9 km 2 of the developed area in 1992 remained unchanged, whereas 27.7 km 2 (30.2%) and 5.38 km 2 (5.9%) of the planted and forest areas, respectively, were primarily converted to developed urban areas from 1992 to 2011. In hydrological modeling, uncertainties in land use data is determined by the sensitivity of the model output to different land use data inputs. Some uncertainties might be associated with different classification algorithms used in both 1992 and 2011 NLCD land use data. Therefore, the presence of low percentages of land use changes between 1992 and 2011 is omitted. Uncertainties and accuracies in NLCD data are also dependent on the interpretation of the person(s) collecting the information and therefore may be assessed differently depending on how it was analyzed. Some uncertainties, therefore, might be applicable to the intended application, while others may have no effects [64].  uncertainties might be associated with different classification algorithms used in both 1992 and 2011 NLCD land use data. Therefore, the presence of low percentages of land use changes between 1992 and 2011 is omitted. Uncertainties and accuracies in NLCD data are also dependent on the interpretation of the person(s) collecting the information and therefore may be assessed differently depending on how it was analyzed. Some uncertainties, therefore, might be applicable to the intended application, while others may have no effects [64].  The Walzem Creek watershed also underwent some land use changes over the past few decades ( Figure 6). During the 20-year period, developed and planted areas were the two largest land use types, and they accounted for approximately 64% and 17% of the total area, respectively. The planted areas shrunk from 1992 to 2011 by 18.3 km 2 . Developed and wetland areas had the greatest increase from 64% to 92% and from approximately 0% to 7.8%, respectively. These increases were due to a large scale, continuous decrease in planted areas (17.5% to 0.8% of the watershed area) and a gradual decrease in forests (9.5% to 2.5%). The increase in wetland areas mostly occurred after 2006 (from 0.06 to 7.82 km 2 ) due to the ecological restoration program for watershed protection that enhanced the urban reaches, bringing the basin back into compliance with water resources and water quality recreation standards. On the other hand, developed areas increased to the detriment of planted and cultivated areas due to the rapid urban development and expansion in the city of San Antonio (Table  7). The Walzem Creek watershed also underwent some land use changes over the past few decades ( Figure 6). During the 20-year period, developed and planted areas were the two largest land use types, and they accounted for approximately 64% and 17% of the total area, respectively. The planted areas shrunk from 1992 to 2011 by 18.3 km 2 . Developed and wetland areas had the greatest increase from 64% to 92% and from approximately 0% to 7.8%, respectively. These increases were due to a large scale, continuous decrease in planted areas (17.5% to 0.8% of the watershed area) and a gradual decrease in forests (9.5% to 2.5%). The increase in wetland areas mostly occurred after 2006 (from 0.06 to 7.82 km 2 ) due to the ecological restoration program for watershed protection that enhanced the urban reaches, bringing the basin back into compliance with water resources and water quality recreation standards. On the other hand, developed areas increased to the detriment of planted and cultivated areas due to the rapid urban development and expansion in the city of San Antonio (Table 7).

SWAT Model Calibration and Validation Results
At the UWBDR, the total observed and simulated streamflow during the calibration period were 7.52 and 7.60 m 3 /s, respectively. The resulting hydrograph from SWAT streamflow in the UWBDR also showed agreement in trends between the two (Figure 7). The best calibration achieved was an R 2 of 0.69, PBIAS of 4.86, ENS of 0.67, and KGE of 0.82. Note that KGE was used as an objective function type in the SUFI-2 calibration and validation because it could be decomposed into three terms that represented the correlation, bias, and relative variability between the measured and simulated values [65]. Hence, it allowed the simultaneous use of baseflow and streamflow in calibration and enabled comparison between different strategies. The summed observed and simulated streamflow during the validation period were 9.03 and 8.27 m 3 /s, respectively. Streamflow validation showed a higher performance than the calibration with an R 2 of 0.84, PBIAS of 23.1, ENS of 0.68, and KGE of 0.67 (Table 8).  [52]. However, the model underestimated the simulated streamflow for the validation period at a monthly time step during low streamflow, which indicates that there may be uncertainty in the results of SWAT simulations for urban watersheds. The lower performance of the SWAT model in the UWBDR may be attributed to the fact that the climate data obtained from the main weather station were located outside the basin, and the distribution of the climate stations with a complete record was sparse. In addition, the overestimation of some peaks in baseflow could be related to the existence of the West Chicago Moraine outwash plain, creating circumstances that promote fast groundwater movement from the moraine through the outwash. Ratios of baseflow to the total annual streamflow were 55.3% and 60.8% for both measured and simulated streamflow, respectively. This discrepancy is acceptable because all of the separation methods of baseflow using different filters are subject to uncertainties [36].  Jan-84 Jan-85 Jan-86 Jan-87 Jan-88 Jan-89 Jan-90 Jan-91 Jan-92 Jan-93 Jan-9 4 J a n -9 5 J a n -9 6 J a n -9 7 J a n - 9 8 Precipitation (  On the other hand, the total annual baseflow during the calibration and validation periods for both measured and simulated data were 8.02 and 7.82 m 3 /s, respectively. Goodness-of-fit measures were evaluated to test the performance of baseflow predictions. The R 2 for the calibration period was 0.67, with a PBIAS of −1.08, ENS of 0.60, and KGE of 0.80 (Table 8). Figure 8 shows the results of model calibration and validation of baseflow at the UWBDR. Overall, there was reasonably good agreement between computed and simulated baseflow. Further, the model performance was validated using data for the subsequent time period. It was observed that the computed baseflow from the USGS streamflow values were reasonably close to the simulated ones. The evaluation indices R 2 , PBIAS, ENS, and KGE were 0.79, 8.43, 0.58, and 0.79 for the baseflow of the validation period, respectively. Unlike the UWBDR and the LEC watersheds reported by Aboelnour et al. [8], the baseflow proportion of the observed and simulated streamflow at the Walzem Creek watershed were 33.3% and 26.8%, respectively, which indicated that surface runoff was a major supply component for the stream. Figure 9 shows the comparison between the simulated and observed monthly streamflow for the calibration and validation periods. USGS records show that the total monthly streamflow for Walzem Creek was 18.7 m 3 /s, while the simulated one was 19.5 m 3 /s. However, streamflow was overestimated for most of the light rainfall events (dry climate periods) and showed very good agreement with the large rainfall events (wet periods). Previous studies have shown that SWAT performed better under more humid climatic conditions [65,66]. In addition, SWAT has some problems with precisely accounting for water loss through infiltration and evapotranspiration, especially during dry climate seasons, and evaluating the soil moisture storage [67][68][69].
During the streamflow calibration period, the R 2 , ENS, PBIAS, and KGE were 0.87, 0.87, −4.31, and 0.91, respectively, while they were 0.83, −3.83, 0.70, and 0.52 during the validation period ( Table  8). The SWAT performance for the monthly streamflow during both the calibration and validation periods was very good [52]. Moreover, the high values of R 2 and ENS in the calibration and validation periods indicated that, with calibrated parameters, the SWAT model was useful to simulate streamflow in semi-arid regions and to further quantify the hydrological impacts of climate variation and land use change over water balance components. Although the SWAT performance for the streamflow validation period was not as good as the calibration period, the results show that its performance was still good, implying that SWAT is applicable to Walzem Creek. The reason that SWAT validation performance was less than the calibration performance is most likely due to the occurrence of an extreme flooding event in October 1998, in which a strong flood killed at least 25 people and caused hundreds of millions of dollars in damages across counties in the southern and eastern regions of San Antonio. The SWAT model poorly matched the peak flow of this large event.
The results also indicate that the simulated values of baseflow were slightly lower than those of the computed ones from observed USGS records. The computed monthly baseflow from USGS records and the simulated one were 6.23 and 5.22 m 3 /s, respectively, during the whole calibration and validation periods. Figure 10  In general, the results suggested that the SWAT model performed satisfactorily in the UWBDR watershed according to the criteria set by Moriasi et al. [52]. However, the model underestimated the simulated streamflow for the validation period at a monthly time step during low streamflow, which indicates that there may be uncertainty in the results of SWAT simulations for urban watersheds. The lower performance of the SWAT model in the UWBDR may be attributed to the fact that the climate data obtained from the main weather station were located outside the basin, and the distribution of the climate stations with a complete record was sparse. In addition, the overestimation of some peaks in baseflow could be related to the existence of the West Chicago Moraine outwash plain, creating circumstances that promote fast groundwater movement from the moraine through the outwash. Ratios of baseflow to the total annual streamflow were 55.3% and 60.8% for both measured and simulated streamflow, respectively. This discrepancy is acceptable because all of the separation methods of baseflow using different filters are subject to uncertainties [36].
Unlike the UWBDR and the LEC watersheds reported by Aboelnour et al. [8], the baseflow proportion of the observed and simulated streamflow at the Walzem Creek watershed were 33.3% and 26.8%, respectively, which indicated that surface runoff was a major supply component for the stream. Figure 9 shows the comparison between the simulated and observed monthly streamflow for the calibration and validation periods. USGS records show that the total monthly streamflow for Walzem Creek was 18.7 m 3 /s, while the simulated one was 19.5 m 3 /s. However, streamflow was overestimated for most of the light rainfall events (dry climate periods) and showed very good agreement with the large rainfall events (wet periods). Previous studies have shown that SWAT performed better under more humid climatic conditions [65,66]. In addition, SWAT has some problems with precisely accounting for water loss through infiltration and evapotranspiration, especially during dry climate seasons, and evaluating the soil moisture storage [67][68][69].
semi-arid region of Walzem Creek, and confirmed that the optimized and calibrated model can be applied to evaluate the responses of the basin's hydrology to land use and climate change.
However, considering the study area was in a semi-arid region and only one meteorological station within the catchment was used, it was difficult to detect whether the climatic conditions in the entire watershed were precisely captured. In addition, the design of the SWAT model may not fully capture the groundwater flow characteristics. However, the outputs are expected to be accurate and reliable since the model was calibrated and validated using observed streamflow.

Impacts of Land Use Change
One of the vital parameters assessed for sustainable management of water resources is water yield. Total water yield is the aggregate amount of water entering the main channel after leaving the HRUs during a time step and can be computed using the following equation [29]: where is the total water yield (mm); is the surface runoff (mm); and are the contributions of lateral flow and groundwater to streamflow (mm), respectively; and is the transmission loss through the bed from the tributary channels in the HRU (mm).
The SWAT simulation suggested that the conversion of the existing planted land cover to urban areas in the UWBDR watershed caused a minor increase in the annual mean water yield by 0.5% During the streamflow calibration period, the R 2 , ENS, PBIAS, and KGE were 0.87, 0.87, −4.31, and 0.91, respectively, while they were 0.83, −3.83, 0.70, and 0.52 during the validation period ( Table 8). The SWAT performance for the monthly streamflow during both the calibration and validation periods was very good [52]. Moreover, the high values of R 2 and ENS in the calibration and validation periods indicated that, with calibrated parameters, the SWAT model was useful to simulate streamflow in semi-arid regions and to further quantify the hydrological impacts of climate variation and land use change over water balance components. Although the SWAT performance for the streamflow validation period was not as good as the calibration period, the results show that its performance was still good, implying that SWAT is applicable to Walzem Creek. The reason that SWAT validation performance was less than the calibration performance is most likely due to the occurrence of an extreme flooding event in October 1998, in which a strong flood killed at least 25 people and caused hundreds of millions of dollars in damages across counties in the southern and eastern regions of San Antonio. The SWAT model poorly matched the peak flow of this large event.
The results also indicate that the simulated values of baseflow were slightly lower than those of the computed ones from observed USGS records. The computed monthly baseflow from USGS records and the simulated one were 6.23 and 5.22 m 3 /s, respectively, during the whole calibration and validation periods. Figure 10 shows the comparison between the simulated and the computed monthly baseflow values at the Walzem Creek watershed in the calibration and validation periods. In the calibration period, the baseflow of the computed and simulated results had a similar trend. Meanwhile, the values of R 2 , ENS, PBIAS, and KGE were 0.85, 21.65, 0.76, and 0.70, respectively, with a P-factor of 70% and R-factor of 0.62. In the validation period, these measures were 0.70, −5.12, 0.68, and 0.79, respectively. The statistical measure results indicate a "very good" to "good" match between the simulated baseflow in the calibration and validation periods and the computed records [52]. However, SWAT overestimated the computed baseflow during the validation period, which was exemplified in the negative values of PBIAS. The statistical indicator and the similar trend between the computed and simulated results showed that the SWAT model was adequate in the semi-arid region of Walzem Creek, and confirmed that the optimized and calibrated model can be applied to evaluate the responses of the basin's hydrology to land use and climate change.
However, considering the study area was in a semi-arid region and only one meteorological station within the catchment was used, it was difficult to detect whether the climatic conditions in the entire watershed were precisely captured. In addition, the design of the SWAT model may not fully capture the groundwater flow characteristics. However, the outputs are expected to be accurate and reliable since the model was calibrated and validated using observed streamflow.

Impacts of Land Use Change
One of the vital parameters assessed for sustainable management of water resources is water yield. Total water yield is the aggregate amount of water entering the main channel after leaving the HRUs during a time step and can be computed using the following equation [29]: where is the total water yield (mm); is the surface runoff (mm); and are the contributions of lateral flow and groundwater to streamflow (mm), respectively; and is the transmission loss through the bed from the tributary channels in the HRU (mm).
The SWAT simulation suggested that the conversion of the existing planted land cover to urban areas in the UWBDR watershed caused a minor increase in the annual mean water yield by 0.5%

Impacts of Land Use Change
One of the vital parameters assessed for sustainable management of water resources is water yield. Total water yield is the aggregate amount of water entering the main channel after leaving the HRUs during a time step and can be computed using the following equation [29]: where W YLD is the total water yield (mm); SUR Q is the surface runoff (mm); LAT Q and GW Q are the contributions of lateral flow and groundwater to streamflow (mm), respectively; and T LOSS is the transmission loss through the bed from the tributary channels in the HRU (mm). The SWAT simulation suggested that the conversion of the existing planted land cover to urban areas in the UWBDR watershed caused a minor increase in the annual mean water yield by 0.5% ( Table 9). The variation could be explained by the reduction in the extent of forests and planted areas and implementation of imperviousness, leading to the reduction of evapotranspiration and infiltration, and increase in surface runoff. However, the reduction of evapotranspiration and the increase in surface runoff were considered not significant in the UWBDR watershed. This could explain the minor increase in the annual mean water yield at the area. Other than the total water yield, the SWAT simulation also suggested a considerable change in baseflow due to the effect of urbanization. It was observed that baseflow increased by 2 mm (accounting for 3.0%) when only the effect of land use dynamics between the two different periods was considered. Figures 11 and 12 show the distribution of the monthly average water yield and baseflow simulated by SWAT, respectively, for the four scenarios for the UWBDR. We observed that the average monthly water yield was concentrated in the late fall/spring seasons and accounted for 29% in the land use change scenario (X2). The change in water yield tended to be positive under the X2 scenario except for the winter season. On the other hand, land use change had minimal effect on baseflow, with no obvious change between X1 and X2. Baseflow variation showed increasing trends in warm months from May to September, and then decreased from October to April. Such increase may be attributed to leakage from an outwash plain at the base of West Chicago Moraine and the increased precipitation during the wet season.
average monthly water yield was concentrated in the late fall/spring seasons and accounted for 29% in the land use change scenario (X2). The change in water yield tended to be positive under the X2 scenario except for the winter season. On the other hand, land use change had minimal effect on baseflow, with no obvious change between X1 and X2. Baseflow variation showed increasing trends in warm months from May to September, and then decreased from October to April. Such increase may be attributed to leakage from an outwash plain at the base of West Chicago Moraine and the increased precipitation during the wet season.   average monthly water yield was concentrated in the late fall/spring seasons and accounted for 29% in the land use change scenario (X2). The change in water yield tended to be positive under the X2 scenario except for the winter season. On the other hand, land use change had minimal effect on baseflow, with no obvious change between X1 and X2. Baseflow variation showed increasing trends in warm months from May to September, and then decreased from October to April. Such increase may be attributed to leakage from an outwash plain at the base of West Chicago Moraine and the increased precipitation during the wet season.   The results in Table 10 show that the average annual water yield increased by 8.0% due to the urbanization effect in Walzem Creek (X2-X1). Meanwhile, urbanization caused the baseflow to experience a reduction by 26.0%. Based on the proposed approach, the average annual evapotranspiration and surface runoff variability during the three scenarios were further analyzed to provide deeper insight into how climate and land use dynamics interacted with hydrologic systems in Walzem Creek watershed. In semi-arid regions, hydrologic systems could be very sensitive to climate variability. Evapotranspiration was an important component of the hydrologic process, often nearly equaling precipitation in the catchment water balance, and, under given climate conditions, it was mainly affected by vegetation cover [70]. Under the same precipitation conditions, decreased evapotranspiration brought an increase in baseflow and streamflow, while increased evapotranspiration led to the reduction in both [71]. This is illustrated in our findings shown in Table 10, in which evapotranspiration experienced a minor reduction by 1.5% due to the land use alteration. However, the reduction of groundwater discharge reported at Walzem Creek Watershed is mainly due to urbanization, agriculture loss, and deforestation. Of note, the evapotranspiration rate is higher in the UWBDR Watershed than in the Walzem Creek Watershed, despite having a higher potential evapotranspiration (PET), as a result of several key differences. The PET recorded in Walzem Creek Watershed was higher than that of the UWBDR Watershed for all scenarios. The PET ranged from 1016 to 1091 mm for the UWBDR Watershed, but was between 1800 and 2084 mm for the Walzem Creek Watershed. However, the amount of precipitation in the Walzem Creek Watershed is lower than that of the UWBDR Watershed, resulting in dryer soils that limit ET. Moreover, the areas of vegetation and forest cover in the UWBDR Watershed are higher for both 1992 and 2011 (Tables 6 and 7). In addition, the UWBDR watershed experienced higher increases in temperature in TS2 compared to TS1 (0.7 • C). Finally, the UWBDR Watershed average daily wind speed was 4.5 and 4.3 ms −1 in TS1 and TS2, respectively. In comparison, the average for the Walzem Creek watershed was 3.9 and 3.8 ms −1 in TS1 and TS2, respectively. Figure 13 illustrates the monthly impacts of land use change, climate, and their joint effect on Walzem Creek's water yield. Land use change had a more pronounced effect for all months in conjunction with a higher monthly average of rainfall in the first period of time (TS1). For example, the monthly average precipitation in June was 120.2 mm in TS1, and decreased to 75.0 mm in TS2. The contributions of land use impacts on monthly water yield were the highest in May, June, October, and December. Conversely, deforestation and urban expansion resulted in a reduction in monthly baseflow in all months from with the highest reduction recorded in the summer season from May to July with a total of 33.9% ( Figure 14). sensitive to climate variability. Evapotranspiration was an important component of the hydrologic process, often nearly equaling precipitation in the catchment water balance, and, under given climate conditions, it was mainly affected by vegetation cover [70]. Under the same precipitation conditions, decreased evapotranspiration brought an increase in baseflow and streamflow, while increased evapotranspiration led to the reduction in both [71]. This is illustrated in our findings shown in Table  10, in which evapotranspiration experienced a minor reduction by 1.5% due to the land use alteration. However, the reduction of groundwater discharge reported at Walzem Creek Watershed is mainly due to urbanization, agriculture loss, and deforestation. Of note, the evapotranspiration rate is higher in the UWBDR Watershed than in the Walzem Creek Watershed, despite having a higher potential evapotranspiration (PET), as a result of several key differences. The PET recorded in Walzem Creek Watershed was higher than that of the UWBDR Watershed for all scenarios. The PET ranged from 1016 to 1091 mm for the UWBDR Watershed, but was between 1800 and 2084 mm for the Walzem Creek Watershed. However, the amount of precipitation in the Walzem Creek Watershed is lower than that of the UWBDR Watershed, resulting in dryer soils that limit ET. Moreover, the areas of vegetation and forest cover in the UWBDR Watershed are higher for both 1992 and 2011 (Tables 6  and 7). In addition, the UWBDR watershed experienced higher increases in temperature in TS2 compared to TS1 (0.7 °C). Finally, the UWBDR Watershed average daily wind speed was 4.5 and 4.3 ms −1 in TS1 and TS2, respectively. In comparison, the average for the Walzem Creek watershed was 3.9 and 3.8 ms −1 in TS1 and TS2, respectively. Figure 13 illustrates the monthly impacts of land use change, climate, and their joint effect on Walzem Creek's water yield. Land use change had a more pronounced effect for all months in conjunction with a higher monthly average of rainfall in the first period of time (TS1). For example, the monthly average precipitation in June was 120.2 mm in TS1, and decreased to 75.0 mm in TS2. The contributions of land use impacts on monthly water yield were the highest in May, June, October, and December. Conversely, deforestation and urban expansion resulted in a reduction in monthly baseflow in all months from with the highest reduction recorded in the summer season from May to July with a total of 33.9% ( Figure 14).  The streamflow changes in the UWBDR watershed appeared to occur in the same manner as changes for the LEC watershed discussed by Aboelnour et al. [8], with some minor differences ( Figure  15). For instance, under X2, streamflow was reduced in winter months by 5.4% to 35.8%. In addition, the average annual water year experienced an increase of 6.7% in the LEC watershed, while it was  The streamflow changes in the UWBDR watershed appeared to occur in the same manner as changes for the LEC watershed discussed by Aboelnour et al. [8], with some minor differences ( Figure 15). For instance, under X2, streamflow was reduced in winter months by 5.4% to 35.8%. In addition, the average annual water year experienced an increase of 6.7% in the LEC watershed, while it was simulated to be only 0.5% and 8.0% in the UWBDR watershed and Walzem Creek watershed, respectively. In contrast to the LEC watershed, which showed a reduction in average annual baseflow as a result of reducing infiltration rate due to urbanization by 28.8%, the UWBDR watershed experienced an average increase in baseflow regardless of the urbanization trend. The reduction in baseflow in the LEC watershed could be caused by over-exploitation and excessive pumping of groundwater used in industry and production [7], while the minor increase in average annual baseflow in the UWBDR might be attributed to flooding of underground structure and the leakage of the groundwater into wastewater systems. The significant decrease in average annual baseflow at Walzem Creek might be due to clearing vegetation, deforestation, and increasing imperviousness, which in turn led to the reduction of evapotranspiration and groundwater discharge while increasing surface runoff. Urbanization is usually associated with measures that play a vital role in accelerating the removal of water from the catchment and stream system, especially during heavy rainfall events. Compacted soil, channelization, and imperviousness allow water to flow rapidly as a result of lower hydraulic resistance of channels and land surfaces of urbanized catchments, which might be an explanation for decreasing baseflow in the Walzem Creek Watershed. The streamflow changes in the UWBDR watershed appeared to occur in the same manner as changes for the LEC watershed discussed by Aboelnour et al. [8], with some minor differences ( Figure  15). For instance, under X2, streamflow was reduced in winter months by 5.4% to 35.8%. In addition, the average annual water year experienced an increase of 6.7% in the LEC watershed, while it was simulated to be only 0.5% and 8.0% in the UWBDR watershed and Walzem Creek watershed, respectively. In contrast to the LEC watershed, which showed a reduction in average annual baseflow as a result of reducing infiltration rate due to urbanization by 28.8%, the UWBDR watershed experienced an average increase in baseflow regardless of the urbanization trend. The reduction in baseflow in the LEC watershed could be caused by over-exploitation and excessive pumping of groundwater used in industry and production [7], while the minor increase in average annual baseflow in the UWBDR might be attributed to flooding of underground structure and the leakage of the groundwater into wastewater systems. The significant decrease in average annual baseflow at Walzem Creek might be due to clearing vegetation, deforestation, and increasing imperviousness, which in turn led to the reduction of evapotranspiration and groundwater discharge while increasing surface runoff. Urbanization is usually associated with measures that play a vital role in accelerating the removal of water from the catchment and stream system, especially during heavy rainfall events. Compacted soil, channelization, and imperviousness allow water to flow rapidly as a result of lower hydraulic resistance of channels and land surfaces of urbanized catchments, which might be an explanation for decreasing baseflow in the Walzem Creek Watershed.

Impacts of Climate Variation
In comparison to the land use change scenario, the climate variation scenario caused the average annual water yield to increase by 9.9% as a result of a prominent increase in precipitation at the UWBDR watershed. Baseflow also showed an increase when only climate variation was considered (X3); however, this increase was much more pronounced than the change in water yield, with an amount of 53.8 mm (81.6%) ( Table 9). These results indicate that both land use change and climate variability played a role in increasing baseflow. However, climate change played a more pronounced role than land use change in impacting the hydrologic regime of the UWBDR during the recent past, due mainly to the increase in precipitation. This can also be seen in Table 9, in which the surface runoff decreased by 9.1% and evapotranspiration increased by 7.0%. Together, these results indicate that the climate alteration contributes more substantially to the effects observed on hydrological components compared to urbanization.
Similar to the land use change scenario, the average monthly water yield was predominantly observed in late fall/spring at the UWBDR watershed. Of note, the highest change in monthly water yield is observed in July (39%) due to the X3 (climate change) scenario. The change in water yield tended to be positive in months that experienced a significant increase in precipitation in the second period (TS2) compared to the first one (TS1) (Figure 11). On the other hand, the results show an increase in average monthly baseflow under the effect of climate change only, impacts of X3 in all months, although the highest growth was detected in the warmest months of the year (May to September) ( Figure 12).
The climate change scenario had the maximum impact on the average annual water yield, causing it to decrease by 11.9% for Walzem Creek Watershed, while it caused the average annual baseflow to decrease by 42.9% (51.7 mm) compared to the baseline scenario (X1) ( Table 10). This may be attributed to the significant reduction in the precipitation pattern and the increase in temperature in TS2 as compared to TS1, where the climate became warmer and drier. Therefore, these likely played an important role in the contribution to the total streamflow for Walzem Creek. Climatic variables, specifically precipitation, largely determined the runoff hydrograph. Precipitation reduction in the second climatic period (TS2) resulted in the significant decline of surface runoff by 20.8 mm (8.2%), and a reduction in evapotranspiration by 2.5%, within the X3 scenario (Table 10). These results indicate that impact of climate variability on baseflow and evapotranspiration was larger than the land use alteration scenario; however, both scenarios had opposite impacts on average annual water yield and surface runoff. Overall, the impacts of climate variation were greater than those of land use change.
On a monthly basis, the highest negative impacts of climate change over the monthly average water yield were detected in June, October, and December, with amounts of 19.2, 12.5, and 12.3 mm, respectively, where the average monthly precipitation was significantly higher in TS2 as compared to TS1. On the other hand, monthly water yield increased at the end of summer and beginning of fall seasons, especially in September; in which it increased by 17.5 mm. It could be inferred that climate variation had a lasting negative effect on water yield ( Figure 13). Similarly, the climate change scenario (X3) caused a reduction in monthly baseflow in all months except August, September, October, and December, with the highest difference recorded in September (3.1 mm). The increase of baseflow in these months was mainly due to changes in precipitation and temperature patterns from TS1 to TS2. For example, TS2 experienced less precipitation as compared to TS1, while the temperature was higher in TS2 compared to TS1. Hence, baseflow played a role in water contribution to total streamflow when the weather got warmer and drier in the semi-arid watershed ( Figure 14). The probable climate alteration impacting most of the globe is mainly determined by the combination of temperature increase and either decrease or increase in rainfall intensity, and any particular baseflow response will depend on the direction and magnitude of both precipitation and temperature. For instance, the change in average monthly precipitation was positive in September, while it was negative in October. Meanwhile, average monthly baseflow showed an increase in both September and October, which might also be explained by the changes in monthly temperature between the two scenarios as it was higher in October than in September.
The streamflow changes under the X3 scenario, which was considered the climate change scenario, were remarkably similar at the UWBDR watersheds to those at the LEC watershed, in which all months exhibited an increase of 12.2-34.5% (Figure 15). In addition, the relative change in streamflow percentage in the UWBDR watershed was higher than the change in the LEC watershed, suggesting that streamflow change is more sensitive to climate change than to land use dynamics. Climate change had a similar impact on the average annual water yield in both the UWBDR and the LEC watersheds, in which it increased by 17.9% and 9.9% in the LEC and UWBDR watersheds, respectively. However, negative impacts due to the X3 scenario occurred in Walzem Creek, indicating that urbanization and climate change had opposite impacts in the semi-arid region of the Walzem Creek Watershed. On the other hand, the climate change caused the average annual baseflow to increase by 15.2% and 81.6% at the LEC and UWBDR watersheds, respectively. However, it declined by 42.9% at the Walzem Creek Watershed. In addition, the average annual surface runoff exhibited an increase in the LEC Watershed due to the impact of climate change by an amount of 22.7%, but decreased in Walzem Creek by 8.2%.
These findings imply that the runoff hydrographs of a catchment are largely impacted by climatic variables, especially precipitation, which in turn affects the percolation of soil water to the groundwater.

Combined Impacts of Both Land Use Change and Climate Variations
To evaluate the combined impacts of land use and climate change, the simulated results under the X4 scenario were compared to the calibrated baseline scenario. The annual mean water yield increased by 11.1% as a response to the X4 scenario at the UWBDR watershed (Table 9). These changes, compared to X2 and X3 scenarios, emphasize that the joint effects of land use change and climate variability led to consistent growth in water yield in the UWBDR watershed. Furthermore, the effect of climate variation was larger than that of the land use dynamic on the total water yield. This can be clearly seen by the X3 and X4 scenarios, in which the mean annual precipitation showed an increase of 81.4 mm, resulting in an increase in the mean annual water yield. These changes are similar to the changes reported in the LEC watershed discussed by Aboelnour et al. [8], resulting from the combined impacts of land use and climate change. In contrast to the LEC watershed, where baseflow decreased, the X4 scenario for the UWBDR watershed led to an increase in the average annual baseflow to 7.1 mm (10.8%) ( Table 9). This difference might be attributed to the prevalence of negative urbanization impacts for the LEC watershed, in contrast to the significant positive effects of climate variation for the UWBDR watershed.
Similar to the climate scenario, we observed that the average monthly water yield was concentrated in the late fall and early spring in the X4 scenario, totaling 39% of the annual yield. In general, positive changes were detected in all months under different scenarios except for November, January, and February. However, the variation due to the joint effects tended to be higher in all months with a higher precipitation pattern in the second period of time (TS2) than in TS1 ( Figure 11). For instance, the effect of land use scenario (X2) was higher than those of X3 and X4 in August, as the average monthly rainfall was 127.2 mm in TS1, while it was only 99.8 mm in TS2. Meanwhile, the combined effect of land use change and climate variability and the sole effect of climate change had greater impacts on water yield in July, as the average monthly precipitation was 83.7 in TS1, increasing to 114.8 mm in TS2. Baseflow variations showed increasing trends in warm months from May to September, then decreased from October to April in conjunction with the joint effect of climate variation and land use change ( Figure 12). The increase in baseflow may be mostly due to an increase in rainfall, and could be explained by fluctuations in both precipitation and temperature between TS1 and TS2. The freeze-thaw processes of the active layer could have changed the soil infiltration capacity and the volume of subsurface water storage, thus impacting baseflow as well [72].
Results from the X4 scenario in the Walzem Creek Watershed indicate that the average annual water yield decreased by 25.7%, while the average annual baseflow showed a consistent reduction by 67.9% (Table 10). Additionally, the annual evapotranspiration was negatively impacted by the joint effect of climate variation and land use change, decreasing by 4.3%. The decline in evapotranspiration is mainly caused by reduction in green cover (Table 7). Compared to X1, the combined effects of land use change and climate variability under X4 decreased surface runoff by 4.2 mm (1.7%). Therefore, with the concurrent reduction in evapotranspiration, average annual water yield, and baseflow had a significant decrease under the X4 scenario These findings indicate that changes of average annual water yield and surface runoff under the joint effects of climate variation and urbanization were smaller than that the changes due to the impacts of the sole impact of climate variation. In other words, the climate alteration had a dominant role, while the land use variation experienced a counteractive role affecting water yield and runoff. The land use change reduced the negative impacts of climate variation by 3.0% and 6.5% for annual water yield and surface runoff, respectively. However, the joint effect of climate variability and land use change on baseflow was higher than the sole impact of the land use change scenario and climate variation scenario. X4 had the greatest negative impact on evapotranspiration. Thus, when the impacts of individual land use change scenario and sole climate alteration happened in the same direction (increase/decrease), the impacts will be intensified when both changes occur at the same time. Of note, the joint effect of climate variation and the urbanization scenario are not a simple summation of each of the individual impacts; however, it represents the interaction of both climate and land use change represented by the SWAT model.
On a monthly basis, the contribution of the joint effects of both climate variability and land use change tended to be similar to the contribution of climate change impacts but with a smaller magnitude for the monthly water yield, with the maximum difference recorded in September (19.0 mm) due to a significant increase in rainfall. Monthly water yield exhibited the highest reduction in June (18.6), due to the notable decline in precipitation in this month. It should be noted that the impact of land use change played a counteractive role for water yield from the Walzem Creek Watershed (Figure 13). Similarly, the highest average monthly change for baseflow was recorded in September (2.5 mm) as a result of the combined impacts of urbanization and climate variation ( Figure 14). Generally, the behavior of average monthly baseflow under the combined impact of land use and the climate change scenario was consistent with the changes under the individual impact of climate variation. At a monthly timescale, the streamflow for Walzem Creek Watershed increased only in August, September, and November considering the climate change scenario and the combined scenario of land use and climate variation, with the highest increase in September by 63.7% (X3) and 71.3% (X4), while it exhibited a reduction in all other months with the highest decline recorded in June by 35.0% (X3) and (33.8) (X4). Note that streamflow in Walzem Creek showed a minor increase in all months when considering the impacts of land use change scenario (X2), except for July and August that showed a minor reduction of 0.7% and 1.1%, respectively ( Figure 16). Moreover, the streamflow rate tended to decrease when considering scenarios X3 and X4, except in August, September, and November, due to the increase in precipitation during these months. The impact of the combined effect of land use change and climate variability showed the same behaviors as the sole impact of climate variability in the Walzem Creek watershed. This situation is well demonstrated by the monthly streamflow variation in the watershed (Figure 16), with the greatest streamflow increase estimated in September at 71.3%. Meanwhile, the highest reduction in monthly streamflow when considering the individual impacts of climate change was estimated in June at 35.0%. These changes were mainly the result of incremental, dynamic precipitation patterns between the two periods, TS1 and TS2. For instance, September experienced the highest increase in rainfall with an amount of 40.4 mm, while June showed the highest reduction in monthly precipitation with an amount of 45.2 mm in TS2 as compared to TS1.

Summary and Conclusions
Urbanization and climate change play an important role in altering the spatiotemporal distribution of water resources and hydrologic components. Streamflow and baseflow are two critically important components of hydrology that are essential to sustain water demands by various sectors, such as agriculture and industry, and are vulnerable to these changes. Therefore, it is of vital significance to understand the behaviors of these components under the separate and combined impacts of climate variation and land use dynamics in different climate regions. In this research, we followed the methodology discussed by Aboelnour et al. [8] for computing streamflow and baseflow for diverse watersheds.
Findings of this research indicate that the climate became warmer and wetter for both the Compared to the LEC watershed, in which the urbanization had the prevailing negative effect on baseflow while climate changes caused increases in both flows, in Walzem Creek watershed, both land use change and climate change had an impact on streamflow and baseflow. However, our study showed that the climate change impacts played a more important role than land use dynamics and urban expansion on streamflow and baseflow in semi-arid regions.

Summary and Conclusions
Urbanization and climate change play an important role in altering the spatiotemporal distribution of water resources and hydrologic components. Streamflow and baseflow are two critically important components of hydrology that are essential to sustain water demands by various sectors, such as agriculture and industry, and are vulnerable to these changes. Therefore, it is of vital significance to understand the behaviors of these components under the separate and combined impacts of climate variation and land use dynamics in different climate regions. In this research, we followed the methodology discussed by Aboelnour et al. [8] for computing streamflow and baseflow for diverse watersheds.
Findings of this research indicate that the climate became warmer and wetter for both the UWBDR and LEC watersheds evaluated by Aboelnour et al. [8] but warmer and drier at the Walzem Creek watershed. The combined effect of these changes showed nonlinear responses to the water balance component. Changes at the UWBDR watershed were remarkably similar to those for the LEC watershed, with the exception that the climate variation was shown to have the greater impact on streamflow, surface runoff, and baseflow, while land use change exerted a relatively small influence on the flow. In other words, in the UWBDR watershed, when the direction of the changes caused by urbanization and climate variation occur in the same direction, the changes of the combined impacts will be intensified. Of note, increasing surface runoff was considered a negative impact as it further strengthened environmental stress by generating more surface erosion and sedimentation. On the other hand, urbanization influenced streamflow positively, while it affected baseflow negatively in the semi-arid Walzem Creek Watershed. However, the climate change had negative impacts on all water components in the area. This might be attributed to the change in rainfall pattern between the two climate periods. The small reduction in mean annual precipitation in the TS2 produced a considerable reduction in runoff. Therefore, the impact of the combined scenario will be amplified when the individual impacts of land use alteration and climate variation are in the same direction (positive/negative). These findings indicate the necessity of evaluating the influences of urbanization and climate alteration separately when assessing the hydrologic effects in urban catchments.
Generally, with the variation in spatiotemporal properties of precipitation, and increasing hazardous events associated with water, such as droughts and floods, stress on water resources will increase and will further encourage the development of mitigation approaches. Based on this research, findings will provide practical suggestions for policy makers on how to sustain water resources more efficiently in relation to its variability as a response to urbanization, land use, and climate change. These changes can be problematic and incur great cost to establish new infrastructure, especially in undeveloped nations. Therefore, policy makers need to develop policies to address these types of changes, taking into account the individual influences of human activities and climate variation, for instance, improving infrastructure to be more resilient to human activities, constructing dams following proper regulations on water resources, and limiting the amount of deforestation, which threatens some hydrological components. In addition, outcomes of this study can be used in quantifying the potential impacts of future projected climate change and land use change. Nevertheless, it might be found that the driving factors interact to impact streamflow and baseflow through chain effects, in which one factor is trying to increase/decrease the magnitude of the other. Hence, more studies are crucial to evaluate this potential future impact on the hydrological system, with the emphasis on the interactive effect of environmental change drivers when predicting future change.
While this research showed the separate and combined impacts of human activity and climate alteration using the SWAT model, modelers should be aware that other types of uncertainties associated with the model exist that may result from observed data, the parameterization process, or from the conceptual model itself. One of the potential shortcomings of this study is that the urbanization processes is an integrated part of the watershed, along with climate alteration. Therefore, it is difficult to discern whether the separate effects of human action and climate change were able to be truly simulated and this issue might therefore create a biased condition. Thus, a suggestion to avoid this limitation in future research is to hypothesize an extreme land use/land cover change that is sensitive to the change instead of a natural system simulated by the model.