Hydrological Modeling with Respect to Impact of Land-Use and LandCover Change on the Runoff Dynamics in Godavari River Basin Using the HEC-HMS Model

Hydrological modeling and the hydrological response to land-use/land-cover changes induced by human activities have gained enormous research interest over the last few decades. The study presented here analyzes the spatial and qualitative changes in the rainfall–runoff that have resulted from the land-cover changes between 1985–2014 in the Godavari River Basin using the Hydrologic Engineering Centre-Hydrologic Modeling System(HEC-HMS) model and remote sensing—GIS (geographic information system) techniques. The purpose of this paper is to analyze the dynamics of land-use/land-cover (LULC) changes for the years 1985, 1995, 2005, and 2014 for the Godavari Basin. The findings reveal an increase of 0.64% of built-up land, a decrease of 0.92% in shrubland, and an increase of 0.56% in waterbodies between 1985–2014. The LULC change detection results between the years 1985–2014 indicated a drastic change in the cropland, forest, built-up land, and water bodies among all of the other classes. The urbanization and agricultural activities are the major reasons for the increase of cropland, built-up land, and water bodies, at the expense of decreases in shrubland and forest. The study had an overall classification accuracy of 92% and an overall Kappa coefficient of 0.9. The HEC-HMS model is used to simulate the hydrology of the Godavari Basin. The analyses carried out were mainly focussed on the impact of LULC changes on the streamflow pattern. The surface runoff was simulated for the year 2014 to quantify the changes that have taken place due to changes in LULC. The observed and the simulated peak streamflow was found to be the same i.e., 56,780 m3/s on 9 September 2014. In the validation part, the linear regression method was used to correlate the observed and simulated streamflow data at the prominent gauge station of the Badrachalam outlet for the Godavari River Basin and give a correlation coefficient value of 0.83. It was found that the HEC-HMS model is compatible and works better for the rainfall–runoff modeling, as it takes into account the various parameters that are influencing the process. The hydrological modeling that was carried out using the HEC-HMS model has brought out the significant impact of LULCC on rainfall–runoff at the Pranhita sub-basinscale, indicating the model’s ability to successfully accommodate all of the environmental and landscape variables. The study indicates that deforestation at the cost of urbanization and cropland expansions leads to decreases in the overall evapotranspiration (ET) and infiltration, with an increase in runoff. The results of the study show that the integration of remote sensing, GIS, and the hydrological model (HEC-HMS) can solve hydrological problems in a river basin.


Introduction
Water Resources Management is an important and integrated approach that includes all of the hydrological components and their linkages with one another.It is very important to understand and quantify the hydrological components for the efficient planning and management of the water resources.Human modifications such as land cover change, irrigation, etc., bring extreme changes to the hydrological variations.The water and energy fluxes of the earth system are integrated through the hydrological cycle, and in turn, it affects the fluxes of the system.Hence, the land-use/land-cover (LULC) changes have a great role in influencing the hydrological cycle.Over the past several years, remote sensing data have had an active part in mapping the LULC changes in various parts of India [1][2][3][4][5][6][7][8][9][10][11][12].Among them, Pooja et al. [7] have quantified the land use/cover of the Gagas watershed in the Almora district using a survey of India topographic sheet from the year 1965 and LISS III satellite data for the year 2008 over a period of 43 years.Rawat et al. [9,10] carried out a study on the land use/land cover of five major towns of Kumaun Himalaya in Uttarakhand (India), and found that the built-up area has sharply increased due to the construction of new buildings in agricultural and vegetation lands.Amin et al. [8] studied the land-use/land-cover mapping of Srinagar city in Kashmir Valley, and observed that significant changes had taken place in the city between 1990-2007.
The LULC changes in and around the river basins bring drastic changes to the environment, both locally and globally.Studies on the consequences of LULC changes on the hydrological cycles have been carried out by [13][14][15].It is very important to assess and understand the hydrological components of a river basin and the impact of LULC for its efficient management of water resources.Hence, mathematical hydrological models-namely the lumped, distributed, and semi-distributed models-play a key role.Among all, the semi-distributed models have their parameters partially varying in space, as the basin is divided into a number of smaller sub-basins.For e.g., the Hydrologic Engineering Centre-Hydrologic Modeling System(HEC-HMS), Soil Water Assessment Tool (SWAT), Hydrologiska Byråns Vattenbalansavdelning model (HBV), etc., hydrological modeling using various hydrological models have been successfully applied for various parts of the river basins in India [16][17][18][19][20][21][22].Dadhwal et al. [16] attempted to model the hydrology of the Mahanadi basin using the VIC hydrological model and assessed the land-cover change impacts on the streamflow at different locations along the river.Pakorn Petchprayoon et al. emphasized the impact of land use/land cover in alleviating the flood problem.Durga Rao et al. [17] studied the hydrological modeling of large river catchments on spatial and non-spatial data.Leong et al. [23] used remote sensing, geographic information system (GIS), and hydrological model for rainfall-runoff modeling.Arash Asadi et al. [18] used HEC-HMS in combination with HEC-GeoHMS for rainfall-runoff modeling of a small sub-basin, Delibajak, and the Kabkian basin by adopting the SCS-CN method.Meenu et al. [19] studied the hydrological impacts of climate change in Tunga-Bhadra river basin of India with HEC-HMS and SDSM.Halwatura et al. [24] worked and found that the HEC-HMS model works good for runoff simulation in a tropical catchment.Roy et al. [20] worked on the calibration and validation of the HEC-HMS model for a river basin in Eastern India, and demonstrated that the use of semi-annual parameters improves the hydrological model.Shinde et al. [22] used the HEC-HMS model to simulate rainfall-runoff modeling for the catchment of Wan reservoir, which is a tributary of Purna River.Nirav Kumar et al. [25] studied the rainfall-runoff modeling for the upstream of Hadaf dam, which is located in the middle Gujarat region, and observed that SCS CN loss model produces slightly better results than the initial and constant loss model.Hence, a proper understanding of the rainfall-runoff relation at different scales for the Godavari river basin is important.
Despite this background, an attempt of hydrological modeling for the Godavari Basin using the HEC-HMS model has not been carried out so far.The present study aims to assess the impact of LULC on the availability of water resources in this region.The main objectives of this study are (i) to assess the past and potential land-cover changes between 1985-2005, and prepare the land-use and land-cover map of the year 2014; (ii) detect the land-use and land-cover changes over the past 30 years; (iii) to assess the water resources in Godavari Basin using the HEC-HMS model and compare the output obtained by the model with the observed ground data such as discharge at the Badrachalam outlet of the basin for the year 2014; and (iv) to study the impact of the LULC change on the rainfall-runoff in the Pranhita watershed, Godavari River Basin.

Study Area: The Godavari River Basin
The Godavari River is the largest of all of the peninsular rivers, with a flow length of 1465 km.It has its origin in Sahyadris near Triambakeswar, at an elevation of 1067m in the Nasik district of Maharashtra, and falls into the Bay of Bengal at Antarvedi, East Godavari.The basin occupies 73 • 26' E to 83 • 07' E longitudes and 16 • 16' N to 23 • 43' N latitudes in the Deccan plateau, with a catchment area of 3,014,503 km 2 (Figure 1).The basin covers the states of Maharashtra (48.6%),Telangana (18.8%),Madhya Pradesh (10%), Chhattisgarh (10.9%),Orissa (5.7%), Andhra Pradesh (4.5%), and Karnataka (1.4%).The majority of the basin is covered with agricultural land (59.57%), while 3.6% is covered with water bodies.The Godavari Basin receives almost 84% of the annual rainfall on an average during the southwest monsoon, which sets in mid-June and ends by mid-October.The major problem in the basin is the frequent flooding in its lower reaches, as the area is on the coastal zone, which is cyclone-prone.Badrachalam, Kunavaram, and the deltaic portions of the river are prone to floods frequently.The delta areas face drainage congestion due to flat topography.The other major concern faced by the basin is the frequent drying up of the Godavari River in the drier months.The river water is highly alkaline, as the coal-fired power stations nearby dump their huge reserves of ash into the river.Badrachalam outlet of the basin for the year 2014; and (iv) to study the impact of the LULC change on the rainfall-runoff in the Pranhita watershed, Godavari River Basin.

Study Area: The Godavari River Basin
The Godavari River is the largest of all of the peninsular rivers, with a flow length of 1465 km.It has its origin in Sahyadris near Triambakeswar, at an elevation of 1067m in the Nasik district of Maharashtra, and falls into the Bay of Bengal at Antarvedi, East Godavari.The basin occupies 73°26′ E to 83°07′ E longitudes and 16°16′ N to 23°43′ N latitudes in the Deccan plateau, with a catchment area of 3,014,503 km 2 (Figure 1).The basin covers the states of Maharashtra (48.6%),Telangana (18.8%),Madhya Pradesh (10%), Chhattisgarh (10.9%),Orissa (5.7%), Andhra Pradesh (4.5%), and Karnataka (1.4%).The majority of the basin is covered with agricultural land (59.57%), while 3.6% is covered with water bodies.The Godavari Basin receives almost 84% of the annual rainfall on an average during the southwest monsoon, which sets in mid-June and ends by mid-October.The major problem in the basin is the frequent flooding in its lower reaches, as the area is on the coastal zone, which is cyclone-prone.Badrachalam, Kunavaram, and the deltaic portions of the river are prone to floods frequently.The delta areas face drainage congestion due to flat topography.The other major concern faced by the basin is the frequent drying up of the Godavari River in the drier months.The river water is highly alkaline, as the coal-fired power stations nearby dump their huge reserves of ash into the river.

Data Used
To find the hydrologic response due to LULC change, the decadal vectorized LULC maps (1985, 1995, and 2005) were obtained from the Indian Space Research Organisation (ISRO) Geosphere Biosphere Programme, soil maps were obtained from NBSS and LUP (Indian National Bureau of Soil Survey and Land Use Planning), and daily rainfall data from the IMD (India Meteorological Department) were used for the study area.The LULC map for the year 2014 that was used in the study was prepared using Landsat 8 Operational Land Imager (OLI) satellite data, which was accessed from Earth Explorer data portal (https://earthexplorer.usgs.gov/).For calibration and validation, the daily discharge data was obtained from the Central Water Commission (CWC) for the gauging sites in the basin.ERDAS Imagine and ArcGIS tools were used for preparing the required input data sets in the present study.The soil textural map of the study area at 1:250,000 scales were obtained from National Bureau of Soil Survey and Land Use Planning, India (NBSS and LUP) in

Data Used
To find the hydrologic response due to LULC change, the decadal vectorized LULC maps (1985, 1995, and 2005) were obtained from the Indian Space Research Organisation (ISRO) Geosphere Biosphere Programme, soil maps were obtained from NBSS and LUP (Indian National Bureau of Soil Survey and Land Use Planning), and daily rainfall data from the IMD (India Meteorological Department) were used for the study area.The LULC map for the year 2014 that was used in the study was prepared using Landsat 8 Operational Land Imager (OLI) satellite data, which was accessed from Earth Explorer data portal (https://earthexplorer.usgs.gov/).For calibration and validation, the daily discharge data was obtained from the Central Water Commission (CWC) for the gauging sites in the basin.ERDAS Imagine and ArcGIS tools were used for preparing the required input data sets in the present study.The soil textural map of the study area at 1:250,000 scales were obtained from National Bureau of Soil Survey and Land Use Planning, India (NBSS and LUP) in Nagpur.Shuttle Radar Topographic Mission (SRTM) Digital Elevation Model (DEM) data of 30-m resolution was used for extracting the topographic characteristics of the basin.The daily rainfall data was obtained from the Indian Meteorological Department (IMD) for the period of 1985 to 2014 (Figure 2).The basin receives almost 84% of its annual rainfall on average during the southwest monsoon, which sets in mid-June and ends by mid-October.It was found that the rainfall varies temporally and spatially across the basin.The data shows the maximum rainfall (1475 mm during 2013-2014) and minimum (855 mm during 1991-1992) during the last 30 years (1985-2014).HEC-HMS v. 4.2 and HEC-GeoHMS were the tools used for the rainfall-runoff modeling, which were downloaded from the United States Army Corps of Engineers (USACE) website http:www.hec.usace.army.mil/software/hec-hms.

Land-Use and Land-Cover Classification
The satellite remote sensing data, with its repetitive nature, has proved to be quite useful in mapping the LULC patterns and their changes over time.The LULC maps for the Godavari Basin for the years 1985, 1995, and 2005 were obtained from the ISRO Geosphere Biosphere Programme (IGBP) to study the past LULC categories.The LULC classes that were originally present in the 1985, 1995, and 2005 LULC maps were merged, and a new set of LULC classes was created.The new ten (10) LULC classes were: built-up, crop land, forest, plantations, fallow land, shrubland, barren land, water bodies, grassland, and mangroves.The vectorized LULC map for the year 2014 was prepared by means of on-screen visual interpretation and digitization techniques based on interpretation elements at 1:50,000 scale by adopting a standard classification scheme [26] with the help of ERDAS imagine and ArcGIS.The visual interpretation was carried out using on-screen digitization [12].Finally, the land-cover mapping output was subjected to accuracy assessment.

Accuracy Assessment
Accuracy assessment is a necessary step to determine and quantify a LULC map.The accuracy assessment was done by comparing the classification product with the reference data, which accurately reflects the true landcover.The accuracy assessment reflects the difference between the classified data and the referenced data.The most common way to represent the classification accuracy of remotely sensed data is an error matrix [27] The user and producer accuracy are two widely used

Land-Use and Land-Cover Classification
The satellite remote sensing data, with its repetitive nature, has proved to be quite useful in mapping the LULC patterns and their changes over time.The LULC maps for the Godavari Basin for the years 1985, 1995, and 2005 were obtained from the ISRO Geosphere Biosphere Programme (IGBP) to study the past LULC categories.The LULC classes that were originally present in the 1985, 1995, and 2005 LULC maps were merged, and a new set of LULC classes was created.The new ten (10) LULC classes were: built-up, crop land, forest, plantations, fallow land, shrubland, barren land, water bodies, grassland, and mangroves.The vectorized LULC map for the year 2014 was prepared by means of on-screen visual interpretation and digitization techniques based on interpretation elements at 1:50,000 scale by adopting a standard classification scheme [26] with the help of ERDAS imagine and ArcGIS.The visual interpretation was carried out using on-screen digitization [12].Finally, the land-cover mapping output was subjected to accuracy assessment.

Accuracy Assessment
Accuracy assessment is a necessary step to determine and quantify a LULC map.The accuracy assessment was done by comparing the classification product with the reference data, which accurately reflects the true landcover.The accuracy assessment reflects the difference between the classified data and the referenced data.The most common way to represent the classification accuracy of remotely sensed data is an error matrix [27].The user and producer accuracy are two widely used measures of accuracy assessment.Producer's accuracy is the percentage of pixels in the reference data for a certain class that are correctly identified by the classifier.User's accuracy is the percentage of pixels classified as a certain class that agree with the reference data.The Kappa coefficient (K) measures the relationship between beyond chance agreement and expected disagreement.The estimate of Kappa is the proportion of agreement after chance agreement is removed from consideration.Kappa statistics/index were computed for a classified map to measure the accuracy of the results.The Kappa coefficient was calculated according to the formula given by [27].
where 'r' is the number of rows in the matrix, x ii is the number of observations in row i and column i, x i+ and x +i are the marginal totals of row i and column i, respectively, and N is the total number of observations [27].

Land-Use and Land-Cover Change Detection Analysis
Change detection is the process of identifying the differences, if any, in the LULC maps of a river basin.It describes and quantifies the differences between images of the same scene at different times.The important aspect of change detection is to identify and thereby conclude the land-use classes where changes have taken place i.e., a class that has been changed into another class over a period of time.This information is extremely useful for management decisions.

Hydrological Modeling Using the HEC-HMS Model 4.1
The HEC-HMS model was developed at the Hydrologic Engineering Center (HEC) by the US Army Corps of Engineers.It was designed to simulate the rainfall-runoff processes of dendritic drainage basins.It was also designed to be applicable in a wide range of geographic areas and suitable for simulating daily, monthly, and seasonal streamflow.The HEC-GeoHMS acts as a bridge between the HEC-HMS and ArcGIS (U.S. Army Corps of Engineers 2000Engineers , 2001Engineers , 2003)).The datasets prepared in GeoHMS are directly used in HEC-HMS to simulate the runoff of the basin.HEC-GeoHMS is used to create basin models using a terrain data digital elevation model.By defining a watershed outlet, HEC-GeoHMS automatically delineates the upstream watershed boundary and preliminary sub-basin outlines.GeoHMS creates a basin model that can be imported into HEC-HMS, and also creates a database table of parameters that were estimated from terrain and other supplementary data layers such as soil and land-use databases.
Land use/land cover is a very important parameter in hydrological modeling.The next major component in hydrological model is infiltration.Infiltration depends on soil texture, which in conjunction with the land use provides various basin parameters for the modeling.A soil textural map of the study area at 1:250,000 scale was obtained from National Bureau of Soil Survey and Land Use Planning of India (NBSS and LUP) in Nagpur, and converted into digital format (Figure 3) with the help of ArcGIS; a vector layer of soil map was thus created.The soil textural classes in the present study are clayey, loamy, clay skeletal, loamy skeletal, and outcrops (Figure 3).Among these, the clayey and loamy are the dominating classes that have the properties of low infiltration rate and more runoff.The digital elevation model (DEM) is another main input for the topographic parameter extraction.Shuttle Radar Topographic Mission (SRTM) DEM data of 30-m resolution was used to extract the topographic and hydraulic parameters of the basin.The sub-basins and drainage network are also delineated using the DEM through an automated process.The SRTM-DEM data of 30-m is used for terrain preprocessing, which includes a series of steps starting with filling sinks, flow direction, and flow accumulation, followed by stream definition, stream segmentation, catchment, and drainage line delineation in HEC-GeoHMS.A new HMS project was created by defining the Badrachalam station as the basin outlet, which includes the delineation of an upstream basin boundary.Keeping the spatial extent of the basin, the computational time, and the desired accuracy, the basin was divided into 145 sub-basins.The basin characteristics such as river length, river slope, and basin slope for each sub-basin were extracted from the DEM.
The soil conservation service curve number (SCS-CN) model was the most widely used in hydrological modeling; it was developed by United States Department of Agriculture (USDA, 1972).In HEC-HMS, for each component of the runoff process, different methods or models are available.Each method has parameters, and the values of these parameters need to be entered as inputs to the model.To compute the losses from the basin, the SCS CN loss method was used, which includes the curve number and percent of imperviousness for each basin.The inputs-the curve number grid and the impervious grid-were prepared using a LULC map, a soil map, and DEM data.The curve number depends upon the hydrological soil group (HSG), land-use type, and AMC-II (antecedent moisture condition).The SCS unit hydrograph transform model is used to compute the direct runoff from excess precipitation.The input parameter for this model is basin lag, which is calculated from the time of concentration for each sub-basin from the excel file generated by the model in the Technical Release 55 (TR55) method.The constant monthly base flow method was used to calculate the base flow, as it represents the amount of ground water that is contributing to the runoff.The base flow for every month has been calculated by means of the straight line method.The Muskingum-Cunge and kinematic wave method were used for routing rivers and channels to the basin outlet, and to estimate the runoff from each sub-basin, the bottom width and side slopes of the reaches are the parameters that are given in the model.
In order to calibrate and validate the model, the model performance is done by (i) visually inspecting and comparing the calculated and observed hydrograph; and (ii) calculating the coefficient of correlation (R 2 ).For the calibration and validation, the Badrachalam outlet of the Godavari River Basin was taken.By defining the Badrachalam station as the outlet point, the Godavari drainage basin was created, and its topographic characteristics were derived by using SRTM 30-m DEM in HEC-GeoHMS, which is an extension tool bar added to the Arc map.The stream network and sub-basins were delineated; its elevations and channel slopes were derived.The model was run for the year 2014 by using the daily rainfall data, the LULC map, a soil map, and the SRTM 30-m DEM.Runoff was computed for the Godavari River Basin.The result of the modeling process is the computation of stream flow hydrographs at the basin outlet, where the main output from the model is discharge at the outlet of the catchment.In HEC-HMS, for each component of the runoff process, different methods or models are available.Each method has parameters, and the values of these parameters need to be entered as inputs to the model.To compute the losses from the basin, the SCS CN loss method was used, which includes the curve number and percent of imperviousness for each basin.The inputs-the curve number grid and the impervious grid-were prepared using a LULC map, a soil map, and DEM data.The curve number depends upon the hydrological soil group (HSG), land-use type, and AMC-II (antecedent moisture condition).The SCS unit hydrograph transform model is used to compute the direct runoff from excess precipitation.The input parameter for this model is basin lag, which is calculated from the time of concentration for each sub-basin from the excel file generated by the model in the Technical Release 55 (TR55) method.The constant monthly base flow method was used to calculate the base flow, as it represents the amount of ground water that is contributing to the runoff.The base flow for every month has been calculated by means of the straight line method.The Muskingum-Cunge and kinematic wave method were used for routing rivers and channels to the basin outlet, and to estimate the runoff from each sub-basin, the bottom width and side slopes of the reaches are the parameters that are given in the model.
In order to calibrate and validate the model, the model performance is done by (i) visually inspecting and comparing the calculated and observed hydrograph; and (ii) calculating the coefficient of correlation (R 2 ).For the calibration and validation, the Badrachalam outlet of the Godavari River Basin was taken.By defining the Badrachalam station as the outlet point, the Godavari drainage basin was created, and its topographic characteristics were derived by using SRTM 30-m DEM in HEC-GeoHMS, which is an extension tool bar added to the Arc map.The stream network and sub-basins were delineated; its elevations and channel slopes were derived.The model was run for the year 2014 by using the daily rainfall data, the LULC map, a soil map, and the SRTM 30-m DEM.Runoff was computed for the Godavari River Basin.The result of the modeling process is the computation of stream flow hydrographs at the basin outlet, where the main output from the model is discharge at the outlet of the catchment.

Results and Discussion
The hydrological response to human activity-induced land-cover changes in major river basins have gained a lot of attention in the past few years due to their effects on water resources.

Land-Use and Land-Cover Map of Godavari Basin
A total of 10 classes were mapped in 1985, 1995, 2005, and 2014 for the basin: built-up, cropland, fallowland, plantation, forest, shrubland, grassland, mangrove, barrenland, and waterbodies (Figure 4).The area statistics are presented in the Table 1.The land-use/land-cover maps analyzed from 1985 to 2014 show that the crop land occupies the maximum area-about 58%-in all of the years from 1985 to 2014.The next predominant class found was forest cover, occupying nearly 30% of the basin area.These two land-use classes contribute maximum evapotranspiration to the basin.During 1985-2014 (Figure 4a-d, Table 1), it has been inferred that there was an increase in cropland (from 57.86% in 1985 to 58.33% in 2014), followed by a decrease in forest cover (29.52-29.33%),a decrease in shrubland (5.98-5.06%),an increase in water bodies (3.08-3.24%),an increase in fallow land (2.09-2.15%),an increase in built-up (0.51-1.15%), an increase in plantations (0.50-0.55%), and an increase in grassland (0.02-0.03%).The results indicate that the increase of built-up land and plantation from 1985 to 2014 showed the urbanization and agricultural activities of the people living in this area.The decrease of forest and shrubland was also observed, which might be utilized for agricultural activities.The results revealed that the overall LULC in the basin showed a reverse trend for the forest cover and cropland classes, i.e., decreasing forests and increasing crop land.The expansion of the built up, cropland, and water bodies may be due to the deforestation and loss of shrubland.The area under the mangroves has not undergone any major change, which may be due to the various protection legislations/ordinances formulated by the Government of India.An increase in the fallowland during 1995 was also found, which could be due to adequate rainfall (Figure 2) and the implementation of canal irrigation.The reason of the enhanced agricultural activity could be due to the construction of mini/minor irrigation tanks.Another important increase was found in the built-up land, from 0.51% in 1985 to 1.15% in 2014.A series of changes due to dam and reservoir construction led to deforestation and cropland conversions.
We have evaluated the accuracy of the LULC 2014 map using 825 random points from the total study area, covering 10 random points from all of the classes.We have used the confusion error matrix, which was created with the mapped and ground reference points, to determine the user's accuracy and Cohen's kappa accuracy (Table 2).It has been identified that out of the 825 points, 759 were classified correctly, showing an overall classification accuracy of 92%, and an overall Kappa coefficient of 0.90.

Results and Discussion
The hydrological response to human activity-induced land-cover changes in major river basins have gained a lot of attention in the past few years due to their effects on water resources.

Land-Use and Land-Cover Map of Godavari Basin
A total of 10 classes were mapped in 1985, 1995, 2005, and 2014 for the basin: built-up, cropland, fallowland, plantation, forest, shrubland, grassland, mangrove, barrenland, and waterbodies (Figure 4).The area statistics are presented in the Table 1.The land-use/land-cover maps analyzed from 1985 to 2014 show that the crop land occupies the maximum area-about 58%-in all of the years from 1985 to 2014.The next predominant class found was forest cover, occupying nearly 30% of the basin area.These two land-use classes contribute maximum evapotranspiration to the basin.During 1985-2014 (Figure 4a-d, Table 1), it has been inferred that there was an increase in cropland (from 57.86% in 1985 to 58.33% in 2014), followed by a decrease in forest cover (29.52-29.33%),a decrease in shrubland (5.98-5.06%),an increase in water bodies (3.08-3.24%),an increase in fallow land (2.09-2.15%),an increase in built-up (0.51-1.15%), an increase in plantations (0.50-0.55%), and an increase in grassland (0.02-0.03%).The results indicate that the increase of built-up land and plantation from 1985 to 2014 showed the urbanization and agricultural activities of the people living in this area.The decrease of forest and shrubland was also observed, which might be utilized for agricultural activities.The results revealed that the overall LULC in the basin showed a reverse trend for the forest cover and cropland classes, i.e., decreasing forests and increasing crop land.The expansion of the built up, cropland, and water bodies may be due to the deforestation and loss of shrubland.The area under the mangroves has not undergone any major change, which may be due to the various protection legislations/ordinances formulated by the Government of India.An increase in the fallowland during 1995 was also found, which could be due to adequate rainfall (Figure 2) and the implementation of canal irrigation.The reason of the enhanced agricultural activity could be due to the construction of mini/minor irrigation tanks.Another important increase was found in the builtup land, from 0.51% in 1985 to 1.15% in 2014.A series of changes due to dam and reservoir construction led to deforestation and cropland conversions.
We have evaluated the accuracy of the LULC 2014 map using 825 random points from the total study area, covering 10 random points from all of the classes.We have used the confusion error matrix, which was created with the mapped and ground reference points, to determine the user's accuracy and Cohen's kappa accuracy (Table 2).It has been identified that out of the 825 points, 759 were classified correctly, showing an overall classification accuracy of 92%, and an overall Kappa coefficient of 0.90.

Change Detection Analysis of the Godavari Basin
The change detection matrices for the change periods of 1985-1995; 1995-2005; and 2005-2014 were calculated, and the results are presented in Tables 3-5, respectively.The land-use/land-cover (LULC) change detection results during 1985 and 2014 also indicated a drastic change in the cropland, form of protection of natural areas and plantation resorts [28].These restoration activities are being carried out by MS Swaminathan Research Foundation (MSSRF) jointly with the Andhra Pradesh Forest Department.The increase in water bodies may be due to the construction of dams such as Bhavali, Chanai, Gautami, Kalpathari, Khadakpurna, Lal, Mukane, ShivanTakli, Upper Ghatghar, Upper Manar, and Wakad in Maharashtra; Kosarteda in Chattisghar; Ragargaon in M.P; Suddavagu in Telangana; etc.

HEC-HMS Model of the Godavari River Basin
For simulating the runoff of the basin for the year 2014, the LULC map of the year 2014 was reclassified into eight classes: built-up land, cropland, fallow land, plantation, forest, shrubland, barren land, and water bodies.The LULC map for the reclassified classes for the year 2014 was created by ArcGIS (Figure 5).The change matrix from 1985 to 2014 shows that mangrove vegetation had a positive change from 1985 to 2014, which may be due to the establishment of the Coringa Wildlife Sanctuary near Kakinada Bay.These positive changes may be due to the active involvement of the government in the form of protection of natural areas and plantation resorts [28].These restoration activities are being carried out by MS Swaminathan Research Foundation (MSSRF) jointly with the Andhra Pradesh Forest Department.The increase in water bodies may be due to the construction of dams such as Bhavali, Chanai, Gautami, Kalpathari, Khadakpurna, Lal, Mukane, ShivanTakli, Upper Ghatghar, Upper Manar, and Wakad in Maharashtra; Kosarteda in Chattisghar; Ragargaon in M.P; Suddavagu in Telangana; etc.

HEC-HMS Model of the Godavari River Basin
For simulating the runoff of the basin for the year 2014, the LULC map of the year 2014 was reclassified into eight classes: built-up land, cropland, fallow land, plantation, forest, shrubland, barren land, and water bodies.The LULC map for the reclassified classes for the year 2014 was created by ArcGIS (Figure 5).The classified LULC map of the year 2014 (Figure 5) was used for creating the CN grid and an impervious grid.The hydrological soil group is an attribute of the soil mapping unit; each soil mapping unit is assigned a particular hydrological Group: A, B, C, and D. The textural classes were grouped into hydrological soil groups [29].The study area is covered with groups B, C, and D (Figure 6).Group B soils had a moderately low runoff potential due to moderate filtration rates.Using the hydro DEM (also known as filled DEM), a union of soil and LULC maps, and the CN look-up table, the curve number grid was generated (Figure 7).Curve numbers were assigned by means of the CN look-up table (Table 6) for each land cover mixed with the type of soil [30].The classified LULC map of the year 2014 (Figure 5) was used for creating the CN grid and an impervious grid.The hydrological soil group is an attribute of the soil mapping unit; each soil mapping unit is assigned a particular hydrological Group: A, B, C, and D. The textural classes were grouped into hydrological soil groups [29].The study area is covered with groups B, C, and D (Figure 6).Group B soils had a moderately low runoff potential due to moderate filtration rates.Using the hydro DEM (also known as filled DEM), a union of soil and LULC maps, and the CN look-up table, the curve number grid was generated (Figure 7).Curve numbers were assigned by means of the CN look-up table (Table 6) for each land cover mixed with the type of soil [30].The impervious grid was generated by converting the land-use impervious percent feature to raster [31].The impervious map was generated, and is shown in Figure 8.These CN values and impervious percent values for each sub-basin were extracted (Table 7).The travel time for each sub-basin was calculated using the TR55 method.The specified hyetograph was considered for the meteorological model, where the daily rainfall data for the year 2014 was given for each sub-basin, and the runoff was simulated at 24 h time interval.The peak discharges were generated for each sub-basin and at the basin outlet.The estimated discharge at the outlet near Bhadrachalam was calibrated against the observed discharge data available from the CWC for the period 2014 (Figure 9).High discharge peaks were observed during the months of July and September.The peak discharge of 56,780 m 3 /s was simulated for 9 September 2014.The trend of the simulated and observed discharge was verified by the R 2 value, which was 0.83.
ISPRS Int.J. Geo-Inf.2018, 7, x; doi: FOR PEER REVIEW www.mdpi.com/journal/ijgimeteorological model, where the daily rainfall data for the year 2014 was given for each sub-basin, and the runoff was simulated at 24 h time interval.The peak discharges were generated for each subbasin and at the basin outlet.The estimated discharge at the outlet near Bhadrachalam was calibrated against the observed discharge data available from the CWC for the period 2014 (Figure 9).High discharge peaks were observed during the months of July and September.The peak discharge of 56,780 m 3 /s was simulated for 9 September 2014.The trend of the simulated and observed discharge was verified by the R 2 value, which was 0.83.

Impact of LULC Change on Runoff
Rainfall is the most important parameter that influences all of the water balance components such as runoff.The total rainfall in the Godavari Basin for the years 1985, 1995, 2005, and 2014 were 928 mm, mm, 908 mm, and 1474 mm, respectively.The year 1995 had the highest recorded rainfall, and these values are highly correlated with the rainfall.The runoff value depends on the rainfall.The LULC change analysis was carried out for all the watersheds for the years 1985, 1995, 2005, and 2014 and the Pranhita watershed was found to have undergone the major LULC change in terms of decrease in forest cover; decrease in shrubland; increase in built-up land, increase in cropland; increase in fallow land; increase in water bodies; and increase in barren land (Table 8).Hence, we have considered the Pranhita watershed to study the effect of LULC change on the runoff.Simulation for the runoff in the Pranhita watershed was carried out for the years 1985, 1995, 2005, and 2014.Figure 10 shows the distribution of runoff for 1985, 1995, 2005, and 2014 for the Pranhita watershed.The trend of runoff was found to be highly dependable on the amount of rainfall received in the entire basin.The yearly trend of rainfall during 1985 was found to be much less than the other years, and the runoff was also found to follow the same trend (Figures 2 and 10).Similarly, for the year 1995, the runoff was recorded as high, which was mainly due to that year also having the highest rainfall.A decrease in the runoff was observed from 1995 to 2005, which may have been due to the reverse trend in the land-cover conversions and/or human activities.In summary, a decrease in the natural cover of the forest over a period of time caused a significant rise in the surface runoff.The surface runoff increases were due to the decrease in infiltration, and hence groundwater recharge processes.The increase in built-up, which was due to urban expansion and intensive cultivation (increase in cropland), will loosen the soil, leading to soil erosion, as well as a decrease in the surface roughness, and thereby ultimately increase the surface runoff.Urbanization also tend to decrease the infiltration rates, and increase the extent of impervious surfaces, although the area over which such changes have occurred is a small fraction of the total watershed.The HEC-HMS model is particularly suitable for the land-use scenarios and their implications on hydrological processes at regional and global scales over the selected time frames.

Impact of LULC Change on Runoff
Rainfall is the most important parameter that influences all of the water balance components such as runoff.The total rainfall in the Godavari Basin for the years 1985, 1995, 2005, and 2014 were 928 mm, 1372 mm, 908 mm, and 1474 mm, respectively.The year 1995 had the highest recorded rainfall, and these values are highly correlated with the rainfall.The runoff value depends on the rainfall.The LULC change analysis was carried out for all the watersheds for the years 1985, 1995, 2005, and 2014 and the Pranhita watershed was found to have undergone the major LULC change in terms of decrease in forest cover; decrease in shrubland; increase in built-up land, increase in cropland; increase in fallow land; increase in water bodies; and increase in barren land (Table 8).Hence, we have considered the Pranhita watershed to study the effect of LULC change on the runoff.Simulation for the runoff in the Pranhita watershed was carried out for the years 1985, 1995, 2005, and 2014.Figure 10 shows the distribution of runoff for 1985, 1995, 2005, and 2014 for the Pranhita watershed.The trend of runoff was found to be highly dependable on the amount of rainfall received in the entire basin.The yearly trend of rainfall during 1985 was found to be much less than the other years, and the runoff was also found to follow the same trend (Figures 2 and 10).Similarly, for the year 1995, the runoff was recorded as high, which was mainly due to that year also having the highest rainfall.A decrease in the runoff was observed from 1995 to 2005, which may have been due to the reverse trend in the land-cover conversions and/or human activities.In summary, a decrease in the natural cover of the forest over a period of time caused a significant rise in the surface runoff.The surface runoff increases were due to the decrease in infiltration, and hence groundwater recharge processes.The increase in built-up, which was due to urban expansion and intensive cultivation (increase in cropland), will loosen the soil, leading to soil erosion, as well as a decrease in the surface roughness, and thereby ultimately increase the surface runoff.Urbanization also tend to decrease the infiltration rates, and increase the extent of impervious surfaces, although the area over which such changes have occurred is a small fraction of the total watershed.The HEC-HMS model is particularly suitable for the land-use scenarios and their implications on hydrological processes at regional and global scales over the selected time frames.all-India monsoon rainfall (11%).The daily rainfall data of the year 2014 were used to simulate the runoff of the basin.The runoff was simulated for 24 h time interval for each year.The peak discharge occurred during July and September.The graph shows that the simulated and observed discharges followed similar trends, with a correlation co-efficient of 0.83.The simulated and observed peak discharges were observed on same date: 9 September 2018.
Although the agreement between the observed and simulated discharges is good, under-estimation and over-estimation are inherent in the simulation.This is because the HEC-HMS model simulates the naturalized flows, and the observed discharge used for validation is based on the human interventions.The model performance showed good agreement at the Bhadrachalam outlet despite the presence of major reservoirs such as Sri Ram Sagar, Indravati, Jayakwadi, and Nizam Sagar, since the calibration was done at this outlet.The model had overestimates (S > O) because of the reservoir storage, due to which the observed flows were less than the simulated values.It may be concluded that the agreement between the observed and simulated values largely depends on the land-cover conditions and hydrological components in the basin, as well as their assumptions.The error may be due to the flood plain inundation in the basin.
The hydrological modeling that was carried out in this study using the HEC-HMS model clearly brought out the significant impact of LULC on the hydrological components of runoff on the Pranhita sub-basinscale, indicating the model's ability to successfully accommodate all of the environmental and landscape variables.The overall increase in annual runoff is well corroborated with the overall increase in rainfall from 1985-1995 and 2005-2014, with the conversion of forest to plantation and forest to cropland in the river basins.This study indicates that deforestation at the cost of urbanization and cropland expansions leads to the decrease in the overall ET, with an increase in runoff.This study has provided valuable insights in the perspective of the subsequent changes in hydrological components as a result of LULC for future prediction, which can be useful in developing management policies to conserve the forests in more intelligent and scientific way.Deforestation, cropland expansion, and urbanization are prominent, and will continue in the upcoming decades.This HEC-HMS modeling study can provide insights into future hydrological scenarios, which will offer the planners to take prior actions for sustainable water use.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 4 of 17 ISPRS Int.J. Geo-Inf.2018, 7, x; doi: FOR PEER REVIEW www.mdpi.com/journal/ijgiNagpur.Shuttle Radar Topographic Mission (SRTM) Digital Elevation Model (DEM) data of 30-m resolution was used for extracting the topographic characteristics of the basin.The daily rainfall data was obtained from the Indian Meteorological Department (IMD) for the period of 1985 to 2014 (Figure 2).The basin receives almost 84% of its annual rainfall on average during the southwest monsoon, which sets in mid-June and ends by mid-October.It was found that the rainfall varies temporally and spatially across the basin.The data shows the maximum rainfall (1475 mm during 2013-2014) and minimum (855 mm during 1991-1992) during the last 30 years (1985-2014).HEC HMS v. 4.2 and HEC GeoHMS were the tools used for the rainfall-runoff modeling,which were downloaded from the United .hec.usace.army.mil/software/hec-hms.

Figure 2 .
Figure 2. Annual variations in the rainfall of the Godavari Basin.

Figure 2 .
Figure 2. Annual variations in the rainfall of the Godavari Basin.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 6 of 17 ISPRS Int.J. Geo-Inf.2018, 7, x; doi: FOR PEER REVIEW www.mdpi.com/journal/ijgiBadrachalam station as the basin outlet, which includes the delineation of an upstream basin boundary.Keeping the spatial extent of the basin, the computational time, and the desired accuracy, the basin was divided into 145 sub-basins.The basin characteristics such as river length, river slope, and basin slope for each sub-basin were extracted from the DEM.The soil conservation service curve number (SCS CN) model was the most widely used in hydrological modeling; it was developed by United States Department of Agriculture (USDA, 1972).

Figure 3 .
Figure 3. Soil map of the Godavari Basin.

Figure 3 .
Figure 3. Soil map of the Godavari Basin.

Figure 5 .
Figure 5. Reclassified land-use/land-cover (LULC) map of the Godavari Basin for the year 2014 showing the various classes of built-up land, cropland, fallow land, plantation, forest, shrubland, barrenland, and waterbodies.

Figure 5 .
Figure 5. Reclassified land-use/land-cover (LULC) map of the Godavari Basin for the year 2014 showing the various classes of built-up land, cropland, fallow land, plantation, forest, shrubland, barrenland, and waterbodies.

Figure 6 .
Figure 6.Hydrological soil group map of the study area.

Figure7.
Figure7.Curve number map of the study area for the year 2014.

Figure 6 .
Figure 6.Hydrological soil group map of the study area.

Figure7.
Figure7.Curve number map of the study area for the year 2014.

Figure 7 .
Figure 7. Curve number map of the study area for the year 2014.

Figure 8 .
Figure 8. Percent imperviousness map of the study area.Figure 8. Percent imperviousness map of the study area.

Table 1 .
Table showing the various land-use/land-cover (LULC) categories of the Godavari River Basin.

Table 2 .
Error matrix or confusion matrix for the accuracy assessment of the LULC map of the Godavari River Basin.

Table 1 .
Table showing the various land-use/land-cover (LULC) categories of the Godavari River Basin.

Table 2 .
Error matrix or confusion matrix for the accuracy assessment of the LULC map of the Godavari River Basin.

Table 6 .
Curve number look-up table for each land cover type with the type of soil.

Figure 6 .
Hydrological soil group map of the study area.

Table 6 .
Curve number look-up table for each land cover type with the type of soil.

Table 6 .
Curve number look-up table for each land cover type with the type of soil.

Table 7 .
The percent imperviousness of the study area for the corresponding land-use classes.
Figure 8. Percent imperviousness map of the study area.

Table 7 .
The percent imperviousness of the study area for the corresponding land-use classes.