Spatiotemporal Analysis of Water Resources in the Haridwar Region of Uttarakhand, India

: Watershed management plays a dynamic role in water resource engineering. Estimating surface runo ﬀ is an essential process of hydrology, since understanding the fundamental relationship between rainfall and runo ﬀ is useful for sustainable water resource management. To facilitate the assessment of this process, the Natural Resource Conservation Service-Curve Number (NRCS-CN) and Geographic Information Systems (GIS) were integrated. Furthermore, land use and soil maps were incorporated to estimate the temporal variability in surface runo ﬀ potential. The present study was performed on the Haridwar city, Uttarakhand, India for the years 1995, 2010 and 2018. In a context of climate change, the spatiotemporal analysis of hydro meteorological parameters is essential for estimating water availability. The study suggested that runo ﬀ increased approximately 48% from 1995 to 2010 and decreased nearly 71% from 2010 to 2018. In turn, the weighted curve number was found to be 69.24, 70.96 and 71.24 for 1995, 2010 and 2018, respectively. Additionally, a validation process with an annual water yield model was carried out to understand spatiotemporal variations and similarities. The study recommends adopting water harvesting techniques and strategies to fulﬁll regional water demands, since e ﬀ ective and sustainable approaches like these may assist in the simultaneous mitigation of disasters such as ﬂoods and droughts.


Introduction
Water assets are the most essential renewable resources required by inhabitants in all forms. Thus, the consumption of water resources needs an effective decision planning for handling both its quality and quantity by considering its spatiotemporal variations. Land use and land cover (LULC) classification describes the role of human beings in affecting the land-cover patterns with time, which ultimately reflects the amount of surface runoff that different types of surface can deal with [1].
Cities experience high surface flows when the ratio of rainfall to infiltration rate is higher [2]. The amount of runoff generated depends upon LULC, soil properties, topography, slope, vegetation cover and atmospheric conditions [3][4][5][6][7][8]. The estimation and storage of stormwater is an essential component in the hydrological cycle to maintain the equilibrium in the city. Estimating surface the spatiotemporal variations of rainfall on surface runoff due to changes in the curve number and to ascertain the applicability and effectiveness of incorporating the NRCS-CN method into geospatial tools. The study also evaluates the temporal effect of anthropogenic activities on spatial surface runoff for the study region, Haridwar, India. This will assist water planners to identify suitable zones for artificially recharge aquifers or water harvesting structures in the region.

Study Area
The study area considered is Haridwar city, Uttarakhand, India. It is one of the seven holiest places to Hindus. Hence, it is essential to estimate water yield availability and the effect of anthropogenic activities, since the area is visited by millions of devotees worldwide. Consequently, the analysis was performed to understand the temporal impact on the landscape of the study region, which is represented in Figure 1. The geophysical characterization of the study region is described in subsequent sub-sections.
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 17 number and to ascertain the applicability and effectiveness of incorporating the NRCS-CN method into geospatial tools. The study also evaluates the temporal effect of anthropogenic activities on spatial surface runoff for the study region, Haridwar, India. This will assist water planners to identify suitable zones for artificially recharge aquifers or water harvesting structures in the region.

Study Area
The study area considered is Haridwar city, Uttarakhand, India. It is one of the seven holiest places to Hindus. Hence, it is essential to estimate water yield availability and the effect of anthropogenic activities, since the area is visited by millions of devotees worldwide. Consequently, the analysis was performed to understand the temporal impact on the landscape of the study region, which is represented in Figure 1. The geophysical characterization of the study region is described in subsequent sub-sections. The Haridwar district is extended up to an area of about 2400 sq. km and located in the southwest part of the Uttarakhand state, India. Haridwar is at the height of 316 m above the mean sea level and lies within the North and Northeast of the Shivalik hills and the Ganges River in the South. The whole area is divided into four parts: mountainous region, upper bajada, lower bajada and alluvial plain. Major land-use classes in the study area are urban, agricultural, forest and barren land. Most of the study area (35% approximately) is forest. Moderate growth is observed, especially in roads and urban classes, as observed in 2018.
The study of climate data indicated that this area corresponds to a moderate subtropical humid region with a temperature drop in March (29.1 °C) that begins to rise until reaching its maximum in May (39.2 °C). In mid-June, when the monsoon starts, the temperature starts falling. It is within the range from 10.5 °C to 6.1 °C from November to February during the winter season. Annual rainfall is about 1200 mm. About 84% of this figure falls within the monsoon season and the rest are in nonmonsoon periods. Maximum precipitation occurs in the Himalayan foothills and decreases gradually while shifting towards the south. In the summer, the mean monthly wind speed is highest in May and June, corresponding to 7.4 and 7.2 km/h, and it decreases to its minimum in October (2.6 km/h). The Haridwar district is extended up to an area of about 2400 sq. km and located in the southwest part of the Uttarakhand state, India. Haridwar is at the height of 316 m above the mean sea level and lies within the North and Northeast of the Shivalik hills and the Ganges River in the South. The whole area is divided into four parts: mountainous region, upper bajada, lower bajada and alluvial plain. Major land-use classes in the study area are urban, agricultural, forest and barren land. Most of the study area (35% approximately) is forest. Moderate growth is observed, especially in roads and urban classes, as observed in 2018.
The study of climate data indicated that this area corresponds to a moderate subtropical humid region with a temperature drop in March (29.1 • C) that begins to rise until reaching its maximum in May (39.2 • C). In mid-June, when the monsoon starts, the temperature starts falling. It is within the range from 10.5 • C to 6.1 • C from November to February during the winter season. Annual rainfall is about 1200 mm. About 84% of this figure falls within the monsoon season and the rest are in non-monsoon periods. Maximum precipitation occurs in the Himalayan foothills and decreases gradually while shifting towards the south. In the summer, the mean monthly wind speed is highest in May and June, corresponding to 7.4 and 7.2 km/h, and it decreases to its minimum in October (2.6 km/h).

Materials and Methods
The flowchart corresponding to the methodology to estimate spatial surface runoff is illustrated in Figure 2. Different inputs to the models were prepared in raster format to analyze spatial and temporal zonal effects. The study has considered temporal variation since 1981, by performing analysis for three years i.e., 1995 (1981-1995), 2010 (1996-2010), and 2018 . The analysis for the years 1995 and 2010 have considered and averaged the past 15 years of datasets whereas for the 2018 year, 8 years of past record is considered for evaluation. Thus, the preparation of input variables to estimate the spatial variation in runoff in different years and understand their correlation is described in subsequent subsections.
2018 year, 8 years of past record is considered for evaluation. Thus, the preparation of input variables to estimate the spatial variation in runoff in different years and understand their correlation is described in subsequent subsections.

Land Use Land Cover
Landsat images were acquired from the United States Geological Survey (USGS) and clipped according to the extent of the study area for the years 1995, 2010, and 2018. Preprocessing analyses were performed on these images using the ERDAS-Imagine software. The LULC classification was achieved by adopting the nearest neighbor classifier and an object-based image analysis approach. As a result, the study area was classified into seven classes, i.e., water, forest, sand, eroded land, urban, rangeland and agriculture with a spatial resolution of 30 m. The classification was validated using ground truth points and Google Earth images.

Soil Depth and Texture
The soil map was acquired from the National Bureau of soil survey and land-use planning (NBSSLUP) at a scale of 1:250,000. In grid format, the resolution was resampled from 1200 m to 30 m to achieve the same spatial resolution of LULC using ArcGIS. Different information about the soil is available in the map such as soil depth, slope, percentage of carbon content, texture, temperature, erosion, drainage, and mineralogy. The required data, viz. soil depth and texture, was transformed into the raster format to perform as an input in the hydrological model.

Land Use Land Cover
Landsat images were acquired from the United States Geological Survey (USGS) and clipped according to the extent of the study area for the years 1995, 2010, and 2018. Preprocessing analyses were performed on these images using the ERDAS-Imagine software. The LULC classification was achieved by adopting the nearest neighbor classifier and an object-based image analysis approach. As a result, the study area was classified into seven classes, i.e., water, forest, sand, eroded land, urban, rangeland and agriculture with a spatial resolution of 30 m. The classification was validated using ground truth points and Google Earth images.

Soil Depth and Texture
The soil map was acquired from the National Bureau of soil survey and land-use planning (NBSSLUP) at a scale of 1:250,000. In grid format, the resolution was resampled from 1200 m to 30 m to achieve the same spatial resolution of LULC using ArcGIS. Different information about the soil is available in the map such as soil depth, slope, percentage of carbon content, texture, temperature, erosion, drainage, and mineralogy. The required data, viz. soil depth and texture, was transformed into the raster format to perform as an input in the hydrological model.
According to its characteristics, soil is classified into four hydrologic soil groups (HSG), i.e., A, B, C and D. The groups are assigned according to the runoff potential and infiltration capacity of the soil. Groups A to D have high to low infiltration capacity and low to high runoff potential. This parameter assists in computing the composite CN and estimating the fraction of rainfall that flows as surface runoff within a city [72].

Precipitation and Temperature
The daily series of temperature and precipitation were collected in gridded format with a spatial resolution of 1 • and 0.25 • , respectively, from the Indian Meteorological Department (IMD). The dataset used in the study is from 1981 to 2018 [73][74][75][76].

Surface Runoff
A combination of the NRCS-CN method with GIS tools was applied to estimate spatial runoff. The NRCS-CN method considers some readily available tables and curves and is based upon a simple empirical formula to obtain surface runoff. The curve number (CN) is an important parameter that influences the estimation of surface runoff. High values of CN represent high amounts of surface runoff and low infiltration capacity of the soil, and vice versa. The CN, which is a function of land use, soil property and slope, was obtained using LULC and soil maps [77].
where Q is direct runoff, P represents total rainfall, S is the potential maximum retention, and I a is the initial abstraction. I a is a function of S that can be expressed as Equation (1) becomes, Here, S can be estimated from the P-Q data for a constant value of I a (0.2 S), and represented in terms of CN using Equation (4).
The curve number is obtained from three variables: land use land cover, hydrologic soil group (HSG), and antecedent moisture condition (AMC) [72]. Hence, the weighted curve number for the study region can be determined by the following formula.
where CN w represents the weighted curve number, A is the total area, CN i is the curve number of sub-catchment i, where sub-catchments (or sub-areas) ranges from 1 to any number N, A i is the area of sub-catchment i, and N is the total number of sub-catchments in a watershed.

Annual Water Yield
The InVEST model has the potential to capture the alteration in the flows due to changes in the ecosystem [60]. This model entirely works according to an empirical function derived from the Budyko framework, which provides the ratio of potential evapotranspiration (PET) to actual evapotranspiration (AET) [78]. The InVEST model captures the spatial variability in precipitation, vegetation, PET, and soil depth to estimate spatial changes in the LULC. The model works on grid format and estimates the heterogeneity in the parameters (LULC, precipitation, temperature, HSG, etc.) affecting the water flows. The spatial water yield is estimated annually for each land use class as follows: where, P(x) represents the spatial annual precipitation and AET(x) refers to the actual annual evapotranspiration at each pixel x. Precipitation and PET determines the mean AET of the watershed. Some of the watershed characteristics, i.e., topography, soil, etc., plays a secondary role. The index of dryness refers to the ratio of annual PET to precipitation to estimate the annual AET.
where ω(x) is a non-physical parameter that represents the natural climatic soil properties and PET at each pixel. PET is further estimated by the following expression.
where Kc(x) represents the vegetation evapotranspiration coefficient, which is a function of the LULC characteristics [79]. Furthermore, ET 0 is the annual reference evapotranspiration, which is determined on the basis of alfalfa grass grown as Equation (11). In addition, ω(x) is an empirical parameter determined by Donohue et al. [80] as follows.
where z is the seasonality factor, the value of which varies from 1 to 30. Further, AWC is a function of the volumetric plant available water content (mm) and is determined through Equation (10).
The maximum depth in the soil up to which the roots can penetrate is referred to as root restricting layer depth. Additionally, the depth up to which 95% of the root biomass occurs is defined as root depth. PAWC is the plant available water content, which generally stands for the difference between field capacity and wilting point. ET o is estimated by the modified Hargreaves method, which is a function of average daily temperature (T avg ), temperature range (TD), ET o , and extraterrestrial radiation (RA). T avg ( • C), is the mean of the mean daily maximum and mean daily minimum temperatures, and TD ( • C) represents the temperature range obtained as the difference between the mean daily maximum and daily minimum temperatures.
To estimate RA, the expression is given as follows: where d r is the inverse relative distance Earth-Sun, RA is extraterrestrial radiation [MJm −2 d −1 ], G sc is a solar constant equal to 0.0820 MJm −2 min −1 , δ is solar declination (rad), w s is sunset hour angle (rad) and ϕ is latitude (rad).
In the present study, this approach was utilized to validate the CN based water yield estimates. The InVEST model was applied to understand the spatiotemporal trend of annual water yield with runoff. Moreover, it enables validating both the range and the spatial variation of temporal runoff for the selected years. Further, it provides a second base to deploy decision planning based on spatial water availability within the study area.

Supervised Classification
The LULC classification for the study period was obtained using supervised classification, as represented in Figure 3. The study area was divided into seven classes, i.e., forest, water, sand, eroded land, urban, rangeland and agriculture.

Supervised Classification
The LULC classification for the study period was obtained using supervised classification, as represented in Figure 3. The study area was divided into seven classes, i.e., forest, water, sand, eroded land, urban, rangeland and agriculture. Table 1    From Table 1, it is understood that a large part of the study area is covered under forest, which is approximately 30% of the total area (254 sq. km). The upper Siwalik region is covered with forest and some built-up areas that include urban forests. Agriculture and rangeland extends to approximately 40%, as agriculture is a major occupation in this region. Water bodies account for approximately 1.4% of the total area since 1995. Most of this water is transmitted through the canal  From Table 1, it is understood that a large part of the study area is covered under forest, which is approximately 30% of the total area (254 sq. km). The upper Siwalik region is covered with forest and some built-up areas that include urban forests. Agriculture and rangeland extends to approximately 40%, as agriculture is a major occupation in this region. Water bodies account for approximately 1.4% of the total area since 1995. Most of this water is transmitted through the canal in the study area, with some small ponds and lakes situated near the residential area. The flashy stream contains water only in the monsoon season. Eroded land is formed due to the erosion of the surface caused by high discharge. The slopes are steep (more than 12%) and streams contain large flows in the monsoon season. To verify the accuracy of this supervised classification, an accuracy assessment was performed as tabulated in Table 2. User accuracy corresponds to the error of commission whereas producer accuracy defines the error of omission in the LULC classification. In other words, user accuracy is the ratio of the number of pixels correctly identified to the number of pixels claimed to be in the respective maps class. The producer accuracy is the ratio of the number of pixels correctly identified in reference plots to the number actually in that reference class. Thus, user and producer accuracy were computed for each class of the study area as described in Table 2.

Trend Analysis of Land Cover Classes
The trend analysis chart is shown in Figure 4. It shows that the forest area reduced from 29.57% (250.46 sq. km) in 1995 to 28.12% (238.176 sq. km) in 2018. The change in forest area is not significant as compared to the increased urban area because of the policies strategized by the government to protect forest. Water class has remained constant since 1995, because most of the water demand is met by the canal. Very little change is observed in the water class as the water flows in the flashy streams. Built-up area increased from 5.5%, (46.58 sq. km) in 1995 to 8.71%, (73.77 sq. km) in 2018. Most of the urbanization is in the lower Siwalik and the plain region. Total cultivable land, which was approximately 40% in 1995, also decreased and became 38% in the year 2018, representing urbanization as directly proportional to cultivation. The eroded land increased by approximately 3% of the area directing the influence of flashy streams and anthropogenic activities in the area. Moreover, the region experienced high sediment loads, which makes it more vulnerable to flood risk. Thus, suitable low impact development (LID) techniques should be implemented at specified zones to mitigate potential negative impacts and prevent socio-economic loss.
was approximately 40% in 1995, also decreased and became 38% in the year 2018, representing urbanization as directly proportional to cultivation. The eroded land increased by approximately 3% of the area directing the influence of flashy streams and anthropogenic activities in the area. Moreover, the region experienced high sediment loads, which makes it more vulnerable to flood risk. Thus, suitable low impact development (LID) techniques should be implemented at specified zones to mitigate potential negative impacts and prevent socio-economic loss.

Curve Number Analysis
Once the raster map was prepared for the soil and LULC, the maps were integrated into the ArcGIS environment to create the curve number map. Hence, the curve number maps obtained for the study area are represented in Figure 5. The higher value of curve number defines the higher runoff for that particular zone and vice versa. The results portray a clear picture of the zones that are experiencing high surface runoff. Thus, to handle drought and flooding situations, water harvesting schemes should be implemented in these specific areas having an abundant amount of surface water availability. In this scenario, the zones having higher CN value represents potential suitable sites for water harvesting.

Curve Number Analysis
Once the raster map was prepared for the soil and LULC, the maps were integrated into the ArcGIS environment to create the curve number map. Hence, the curve number maps obtained for the study area are represented in Figure 5. The higher value of curve number defines the higher runoff for that particular zone and vice versa. The results portray a clear picture of the zones that are experiencing high surface runoff. Thus, to handle drought and flooding situations, water harvesting schemes should be implemented in these specific areas having an abundant amount of surface water availability. In this scenario, the zones having higher CN value represents potential suitable sites for water harvesting.

Variation in Area Ratio
The curve number is a function of the LULC and soil property of the area. Figure 6 shows the spatial distribution of the CN during the years 1995, 2010 and 2018. Area ratio represents the ratio of a particular land covered with a specific number of CN values to the total area. It is understood from

Variation in Area Ratio
The curve number is a function of the LULC and soil property of the area. Figure 6 shows the spatial distribution of the CN during the years 1995, 2010 and 2018. Area ratio represents the ratio of a particular land covered with a specific number of CN values to the total area. It is understood from Figure 6 that the CN is increasing from 1995 to 2018. The initial range of CN  belongs to the forest area, which is decreasing from 1995 to 2018. The agriculture and crop areas are also reduced, and the maximum deviation in the CN is in the range of 83-87, which belongs to wasteland and built up area. A decrease in forest and cultivated land and an increase in wasteland and built-up area are observed in the city due to the effect of urban sprawl.

Runoff Analysis
Monthly runoff was estimated for the study area by providing rainfall data as an input to the NRCS-CN method, as represented in Figure 7. For the year 1995, maximum rainfall was observed in August, which corresponds to AMC III condition and a maximum amount of runoff for this period. About 85% of the yearly rainfall is observed in July to August, contributing to maximum runoff in these months. AMC I condition exists for the rest of the months, as they received less rainfall and thus resulted in a lower amount of runoff. Thus, an adjustment table was referred to select the curve number corresponding to different soil moisture conditions [81].
For the year 2010, the area receives a lower amount of rainfall in the non-monsoon period, whilst approximately 87% of total rainfall occurs during the monsoon period. The AMC III condition is applicable for the monsoon region, thereby causing a higher amount of runoff in this period. Some of the rainfall is observed in the winter season due to the western disturbances, but runoff is nearly zero because of the non-saturation condition of the soil. As for 2018, rainfall is very low as compared to previous years. Because of the lesser rainfall in the area, the AMC II condition is applicable for June, July and August. Hence, the estimated runoff is also very reduced in relation to previous years.
The observed mean annual rainfall for the Haridwar district is approximately 1174.3 mm, which is also 85% of the rainfall observed for the monsoon period, i.e., from July to September. Out of the three years under study, 2010 was found to receive the maximum amount of rainfall and runoff, as shown in Figure 7. As the data of the Ganges basin is classified, it is challenging to present the complete results along with their validation in the paper.
The computation was performed to analyze the impact of climate change and anthropogenic activities on the regional surface runoff potential. The study infers that runoff increased approximately 48% from 1995 to 2010 and 71% from 2010 to 2018 (Figure 8). The values of weighted

Runoff Analysis
Monthly runoff was estimated for the study area by providing rainfall data as an input to the NRCS-CN method, as represented in Figure 7. For the year 1995, maximum rainfall was observed in August, which corresponds to AMC III condition and a maximum amount of runoff for this period. About 85% of the yearly rainfall is observed in July to August, contributing to maximum runoff in these months. AMC I condition exists for the rest of the months, as they received less rainfall and thus resulted in a lower amount of runoff. Thus, an adjustment table was referred to select the curve number corresponding to different soil moisture conditions [81].
For the year 2010, the area receives a lower amount of rainfall in the non-monsoon period, whilst approximately 87% of total rainfall occurs during the monsoon period. The AMC III condition is applicable for the monsoon region, thereby causing a higher amount of runoff in this period. Some of the rainfall is observed in the winter season due to the western disturbances, but runoff is nearly zero because of the non-saturation condition of the soil. As for 2018, rainfall is very low as compared to previous years. Because of the lesser rainfall in the area, the AMC II condition is applicable for June, July and August. Hence, the estimated runoff is also very reduced in relation to previous years. The CN values were found to have very little temporal variability for the study area (Figure 7). This suggested that LULC, soil map or variations in land-use class have low influence on surface runoff. Changes in the rainfall pattern and intensity are solely responsible for the flooding condition in the study region. The city experiences runoff mainly during July to September. In July, runoff was maximum in 2010 (approximately 100% greater than in 1995 and 245% greater than in 2018). Subsequently, flooding can be mitigated by implementing water harvesting schemes at the specified areas that experience high surface runoff. Furthermore, to validate the annual spatial trend of runoff, the InVEST model is applied to estimate the annual water yield.

Validation Using Alternative Approach
The Himalayan study region has very little data available for validation due to the political restriction and physical difficulties. The mountainous region of the Himalayas are not equipped with any gauging stations. The downstream of the catchment has measuring stations available. However, due to the data sharing constrained by government policies, it is very difficult to validate the model with the observation stations. This results in a lack of data availability for the research community resulting in little hydrological and water resources analysis in the Ganges basin. Therefore, to capture and validate the spatial trend of runoff estimated by the NRCS-CN model, the InVEST model was applied to the study region. The annual spatiotemporal variations in water yield (InVEST) and surface runoff (NRCS-CN) are represented in Figure 8. It has been observed that the range of yield and runoff for the selected years are similar; however, they vary in their spatial trend. This is a very interesting finding, whereby both approaches show spatial variability of runoff and yield despite the similarity of yield and runoff in different years. Due to the variation in the Budyko and CN approaches, different classes represent different runoff potential with the same range of surface runoff for all the selected years. The observed mean annual rainfall for the Haridwar district is approximately 1174.3 mm, which is also 85% of the rainfall observed for the monsoon period, i.e., from July to September. Out of the three years under study, 2010 was found to receive the maximum amount of rainfall and runoff, as shown in Figure 7. As the data of the Ganges basin is classified, it is challenging to present the complete results along with their validation in the paper.
The computation was performed to analyze the impact of climate change and anthropogenic activities on the regional surface runoff potential. The study infers that runoff increased approximately The CN values were found to have very little temporal variability for the study area ( Figure 7). This suggested that LULC, soil map or variations in land-use class have low influence on surface runoff. Changes in the rainfall pattern and intensity are solely responsible for the flooding condition in the study region. The city experiences runoff mainly during July to September. In July, runoff was maximum in 2010 (approximately 100% greater than in 1995 and 245% greater than in 2018). Subsequently, flooding can be mitigated by implementing water harvesting schemes at the specified areas that experience high surface runoff. Furthermore, to validate the annual spatial trend of runoff, the InVEST model is applied to estimate the annual water yield.

Validation Using Alternative Approach
The Himalayan study region has very little data available for validation due to the political restriction and physical difficulties. The mountainous region of the Himalayas are not equipped with any gauging stations. The downstream of the catchment has measuring stations available. However, due to the data sharing constrained by government policies, it is very difficult to validate the model with the observation stations. This results in a lack of data availability for the research community resulting in little hydrological and water resources analysis in the Ganges basin. Sustainability 2020, 12, x FOR PEER REVIEW 12 of 17 However, low temporal variations are observed in the weighted curve number values between the selected years. Therefore, yield or runoff entirely depend upon the rainfall trend of the study region. To understand and deploy site specific strategic planning, it is required to analyze the monthly variation and trend of rainfall. Thus, monthly rainfall data was observed and analyzed effectively for decision-making and policy purposes. The monthly analysis provides a clear assessment of surface water runoff for the three selected years. It has been observed that the city experiences high rainfall in July, August and September which results in high surface runoff and flooding conditions. By considering the importance of the city, proper management of water resources should be strategized for the specified months based on the CN analysis, which provides a detailed assessment and location of the zones that are vulnerable to high surface runoff.
Furthermore, the study indicated that most of the months remain dry. Thus, the city will depend on natural resources or storage tanks for potable and non-potable water demands. If the storage tanks are constructed and utilized during rainfall periods, they can be used to store water for dry months and thus reduce the exploitation of natural resources. In this context, it is recommended to implement stormwater harvesting schemes and storage tanks in the desired zones. This course of action will address the increasing demand and low water resources availability, making water accessible for inhabitants and civilians. Geospatial tools can assist in obtaining more precise information in rainfallrunoff modeling about the catchment size and characteristics. Further, the analysis can be performed swiftly by integrating composite land-use classes with diverse soil types. Thus, effective planning can be deployed at the sites of interest to have an abundant availability of supply during dry periods. Additionally, awareness programs can be targeted at specified zones for water conservation and mitigating disasters such as floods and droughts. Therefore, to capture and validate the spatial trend of runoff estimated by the NRCS-CN model, the InVEST model was applied to the study region. The annual spatiotemporal variations in water yield (InVEST) and surface runoff (NRCS-CN) are represented in Figure 8. It has been observed that the range of yield and runoff for the selected years are similar; however, they vary in their spatial trend. This is a very interesting finding, whereby both approaches show spatial variability of runoff and yield despite the similarity of yield and runoff in different years. Due to the variation in the Budyko and CN approaches, different classes represent different runoff potential with the same range of surface runoff for all the selected years.
However, low temporal variations are observed in the weighted curve number values between the selected years. Therefore, yield or runoff entirely depend upon the rainfall trend of the study region. To understand and deploy site specific strategic planning, it is required to analyze the monthly variation and trend of rainfall. Thus, monthly rainfall data was observed and analyzed effectively for decision-making and policy purposes. The monthly analysis provides a clear assessment of surface water runoff for the three selected years. It has been observed that the city experiences high rainfall in July, August and September which results in high surface runoff and flooding conditions. By considering the importance of the city, proper management of water resources should be strategized for the specified months based on the CN analysis, which provides a detailed assessment and location of the zones that are vulnerable to high surface runoff.
Furthermore, the study indicated that most of the months remain dry. Thus, the city will depend on natural resources or storage tanks for potable and non-potable water demands. If the storage tanks are constructed and utilized during rainfall periods, they can be used to store water for dry months and thus reduce the exploitation of natural resources. In this context, it is recommended to implement stormwater harvesting schemes and storage tanks in the desired zones. This course of action will address the increasing demand and low water resources availability, making water accessible for inhabitants and civilians. Geospatial tools can assist in obtaining more precise information in rainfall-runoff modeling about the catchment size and characteristics. Further, the analysis can be performed swiftly by integrating composite land-use classes with diverse soil types. Thus, effective planning can be deployed at the sites of interest to have an abundant availability of supply during dry periods. Additionally, awareness programs can be targeted at specified zones for water conservation and mitigating disasters such as floods and droughts.

Conclusions
Estimation of surface runoff helps in sustainable planning and management of land use and available water content. Towards this, the NRCS-CN method was initially applied for estimating the regional composite curve number. The present methodology reduces time as well as the efforts to deploy strategic planning by focusing on spatial water availability. The method also yields comparable estimates as those obtained using the Budyko model. However, the spatial variability of water yield is observed in both approaches. It has been analyzed that the city has not experienced much temporal variation in the weighted curve number since 1995. In other words, changes in land use and topography due to the anthropogenic influences since 1995 will not have much influence on flow variability. Hence, runoff is directly proportional to rainfall, which has substantially changed due to the effect of climate change. Thus, the implementation of suitable conservation practices and structural works is suggested to reduce the pressure from freshwater resources within the city. As a fact, cities are facing natural calamities very often due to urban sprawl, such as droughts or flooding. Nevertheless, if suitable zones are identified and combined with the knowledge of water availability and demand, strategic plans can be deployed spatially with ease. Additionally, if the conservation practices are suitably placed, urban flooding can be embraced as an opportunity to treat the regions with water demands during dry periods. Henceforth, it is highly recommended to implement some low impact development techniques within the zones with high water availability, so that the dependency on natural resources is reduced and the effect of disaster could be mitigated with time.