The Impact of Reforestation Induced Land Cover Change (1990–2017) on Flood Peak Discharge Using HEC-HMS Hydrological Model and Satellite Observations: A Study in Two Mountain Basins, China

: Understanding the e ﬀ ect of land use and land cover (LULC) type change on watershed hydrological response is essential for adopting applicable measures to control ﬂoods. In China, the Grain to Green Program (GTGP) and the Natural Forest Conservation Program (NFCP) have had a substantial impact on LULC. We investigate the e ﬀ ect of these conservation e ﬀ orts on ﬂood peak discharge in two mountainous catchments. We used a series of Landsat images ranging from 1990 to 2016 / 2017 to evaluate the LULC changes. Further to this, the hydrological responses at the basin and sub-basin scale were generated by the Hydrologic Modeling System (HEC-HMS) under four LULC scenarios. Between 1990 and 2016 / 2017, both catchments experienced an increase in forest and urban land by 18% and 2% in Yanhe and by 16% and 8% in Guangyuan, respectively. In contrast, the agricultural land decreased by approximately 30% in Yanhe and 24% in Guangyuan, respectively. The changes in land cover resulted in decrease in ﬂood peak discharge ranging from 14% in Yanhe to 6% in Guangyuan. These ﬁndings provide a better understanding on the impact of reforestation induced LULC change on spatial patterns of typical hydrological responses of mountainous catchment and could help to mitigate ﬂash ﬂood hazards in other mountainous regions. Author Contributions: Conceptualization, R.L. and C.K.; Methodology, C.K.; Software, Y.Z.; Validation, Y.W.; Formal Analysis, M.M.; Investigation, J.G.; Resources, R.L. and J.G.; Data Curation, Z.G.; Writing—Original Draft Preparation, C.K. and R.L.; Visualization, D.E.R.R.; Supervision, R.L.; Project Administration, R.L.; Funding


Introduction
Globally, flash floods are the most frequent and devastating natural disaster that pose a significant threat to human life and cause damage to the economy [1][2][3][4]. Modrick and Georgakakos [5] defined flash Water 2020, 12 floods as an event where streamflow increases and decreases fast without warning as a result of severe rainfall. LULC and climate changes can affect catchments' hydrologic processes and cause changes in flash flood frequencies [6]. For example, changes in LULC as a result of various anthropogenic activities alter the catchment's characteristics and aggravate the fast hydrological responses to result in destructive floods in low-lying areas [7]. Moreover, flash floods are likely to occur in mountainous catchments of a few hundred square kilometers or less, where intense rainfall is converted into runoff from upper stream areas [8,9]. In China, floods have caused loss of life and infrastructure damages in mountainous and hilly areas [10]. Studies [11,12] showed that between 2006 and 2015, flood disasters caused the death of 6641 people and affected half a billion other people in China. Additionally, they provoked more than $87.5 billion USD in economic loss. For example, the flash flood that hit the capital city of China, Beijing, in 2012, caused the death of 77 people, with an estimated economic loss of US $1.6 billion. Some other studies [13] showed that 641 Chinese regions and cities are likely to be affected by floods. Furthermore, Miao et al. [14] observed that most of the flash floods in China occur regularly in mountainous catchments with poor vegetation cover.
The Chinese government introduced policies and other initiatives to reduce the impact of flooding. For example, the Natural Forest Conservation Program (NFCP) and Grain to Green Program (GTGP), which are the two main wide programs of reforestation and afforestation, were introduced with the aim of reducing flash floods and erosion in many regions [15,16]. Other measures such as forecasting and warning systems and sustainable land use (e.g., green spaces) for flood mitigation were introduced in various regions [17]. For instance, vegetation and green spaces have various functions on catchment hydrology, including reducing flash floods, as they serve as barriers to surface runoff [18]. Several studies [19,20] demonstrated that the GTGP and NFCP conservation policies influenced the land use/cover changes in various regions in China. However, few studies investigated the effect of these two policies on catchment hydrological response in mountainous regions of China.
In the past decades, some authors showed that the decrease in forest covers was the main cause of the increase in surface runoff as well as flooding [21][22][23]. Other studies [24,25] attributed the increase in runoff and peak flows to both climate and land use changes. However, Liu et al. [26] and Chu et al. [27] argued that the change in land use only causes a slight change in runoff and streamflow. Alkama et al. [28] also found no strong correlation between changes of land use types and runoff generation. Wilk et al. [29] demonstrated that the alternations in land use/cover are insignificant to the hydrological response in larger watershed because surface runoff takes place at a small part of the watershed area, which causes a slight change in peak discharge. Moreover, Cavalcante et al. [30] pointed out that it may be difficult to understand the hydrological response of land use changes in large-scale basins that have a mixture of land use classes in different phases of protection and regeneration. The aforementioned studies showed that the relationships between changes in land cover and catchment hydrological response remain a controversial topic and will be constantly researched in the future.
A limited amount of research has assessed the effectiveness of the NFCP and GTGP in reducing flash flood and erosion in mountainous regions in China [31,32]. Yang et al. [33] and Zhao et al. [34] assessed the impact of GTGP on vegetation cover change and soil erosion in the Loess Plateau and noted a decrease in soil erosion. Zhang et al. [35] observed a reduction in soil erosion in most southern regions of China as result of land cover change induced by GTGP and NFCP. On the contrary, Xie et al. [36] found that the change in land cover caused by afforestation activities in China's Three-North region had insignificant effects on hydrological cycle. Most of these studies are inconsistent concerning the impact of land cover changes on catchment hydrological response in various regions in China.
More important, some of these studies were focused on the relationship between vegetation cover changes and soil erosion in non-mountainous regions and mostly using statistical methods and on short timescales. Therefore, there is a need to extend previous studies to mountain area and conduct quantitative analysis based on long historical data and hydrological modeling. The goal of the present study is to investigate the effect of alternation in land cover induced by the NFCP and GTGP on flash flood peak discharge in two medium-sized (500-8000 km 2 ) mountainous watersheds in China. This study does not consider the impact of climate change. In summary, it attempts to address the following specific objectives: (1) to quantify the changes in land use land cover induced by NFCP and GTGP in the two catchments; (2) to investigate the temporal evolution of flash flood peak discharge with diverse historical land cover changes using the HEC-HMS hydrological model; (3) to analyze the spatiotemporal trends of peak discharge and vegetation index at the sub-basin scale. The results of this study will improve our understanding of how changes in land cover affect catchment hydrological response and serve as a reference formulating mitigation strategies to prevent future flash floods disaster in mountainous regions.

Study Area
The present study was carried out in two mountainous catchments with different size and climate conditions, located in China. The two catchments have undergone changes in land cover at different degrees in the last 25 years, and hence, serve as ideal sites to evaluate the impact of reforestation induced land use land cover change on flash flood peak discharge. The first study area is the Yanhe catchment, and is located in the middle regions of Loess Plateau, in the northern Shaanxi province in China (36 • 21 -37 • 19 N; 108 • 38 -110 • 29 E) ( Figure 1c). The catchment covers an area of 7724 km 2 with an altitude range from 484 to 1797 m above sea level. It is characterized by a steep gullies region and has a warm temperate continental semi-arid monsoon climate. The Yanhe catchment has an average annual precipitation of 495 mm with over 70% falls in the summer (June-September) and average annual temperature varies from 8.8 to 10.2 • C. The main land use types are grassland, forest and agricultural land. This study catchment is affected by soil erosion and floods, and is among the regions where the Chinese government introduced many land use conservation programs [12].
The second study area is the Guangyuan catchment, one of the Sichuan basin tributaries and is located in the Southwest of China (32 • 16'-32 • 35' N; 105 • 50'-106 • 11' E) ( Figure 1d). This study catchment has a drainage area of 578 km 2 . The elevation of the catchment ranges from 800 to 1400 m above sea level and the slopes range from 0% to 44%, respectively. The catchment has a subtropical monsoon climate with an annual rainfall that ranges between 1200 to 1600 mm. Land use is dominated by forest and agriculture. Urban areas are found mainly in low laying areas in proximity to the main river. The catchment is located in the region that was selected for forest cover expansion under various conservation programs such as GTGP and NFCP. Moreover, it has experienced various floods in the past that have resulted in various types of disasters such as loss of human lives and economic losses [37]. Hence, the Guangyuan catchment serves as an ideal site to evaluate the impact of reforestation induced land use land cover change on flood peak discharge.

Sources of Data
The dataset used in this study includes the satellite images data from the Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) for the years 1990, 2000, 2002, 2008, 2010, 2016 and 2017. Both Landsat sensors have the same spatial resolution of 30 m and cover each of our study area. The free available geometrically corrected Landsat L1T data were downloaded (http://www.earthexplorer.usgs.gov) and used in the present study. We selected clear satellite images with less than 10% cloud and cloud shadows. The Landsat images from the same season of the year and with close time of acquisition by the sensors were selected. This facilitated in comparing Landsat images of different periods. We also used land use maps of China (http://www.resdc.cn) for validation of the Landsat images. The soil data of the two study catchments (soil type and texture) were downloaded from the world soil database (http://www.fao.org/data/en/) [38]. The hourly rainfall data (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011) were provided by the China Meteorological Administration and were collected from various meteorological stations located inside each study catchment and within a maximum distance of 10 km from the basin perimeter ( Figure 1). The daily discharge data (2000-2012) of the Yanhe catchment were acquired from the Yellow River Conservancy Commission (YRCC) and used for model calibration and validation. Lastly, the Digital Elevation Model (DEM) with a spatial resolution of 30 m from the Shuttle Radar Topography Mission (SRTM) was accessed from the United States Geological Survey website.
Water 2020, 12, x FOR PEER REVIEW 4 of 23 Figure 1. The location of the two study catchments used in this study in China (a), location of study area catchments in Shaanxi and Sichuan provinces (b), the Yanhe (c) and Guangyuan (d) study area boundary; elevation with the stream network is shown with blue solid lines and red triangles represent meteorological stations.
The second study area is the Guangyuan catchment, one of the Sichuan basin tributaries and is located in the Southwest of China (32°16'-32°35' N; 105°50'-106°11' E) (Figure 1d). This study catchment has a drainage area of 578 km 2 . The elevation of the catchment ranges from 800 to 1400 m above sea level and the slopes range from 0% to 44%, respectively. The catchment has a subtropical monsoon climate with an annual rainfall that ranges between 1200 to 1600 mm. Land use is dominated by forest and agriculture. Urban areas are found mainly in low laying areas in proximity to the main river. The catchment is located in the region that was selected for forest cover expansion under various conservation programs such as GTGP and NFCP. Moreover, it has experienced various floods in the past that have resulted in various types of disasters such as loss of human lives and economic losses [37]. Hence, the Guangyuan catchment serves as an ideal site to evaluate the impact of reforestation induced land use land cover change on flood peak discharge.

Sources of Data
The dataset used in this study includes the satellite images data from the

Satellite Image Processing and Land Use Land Cover Classification
The atmospheric effect removed Landsat reflectance images were used in this study. The supervised image classification method was used [39] to perform land use land cover classification. High-resolution imagery from Google Earth pro served as the reference map for training during Landsat image classification. Additionally, other information from China land use maps (2000,2005,2010 and 2015 accessible at http://www.resdc.cn) related to the land cover pattern of the two study areas were taken into consideration during on-screen training sample selection. Google Earth pro provides a tool to access historical high-resolution imagery of the area of interest with dates similar or close to that of the Landsat images [40]. In our case, all the imagery covering both study catchments were available in Google Earth pro. Other studies [41,42] also demonstrated that Google Earth is a reliable source for selecting training samples in satellite image classification. We selected randomly 20 training samples per land cover class [43] and used the Maximum Likelihood Classifier (MLC) algorithm of supervised classification for land cover classification of each Landsat image. The MLC algorithm is a robust and optimal classifier and is the most used algorithm in the assessment of land cover changes that produces high accurate results [44]. Therefore, we used the training samples for selecting representative classes for each category with the aim of defining various spectral signatures. From the spectral signature, each pixel in the Landsat satellite image was assigned to the appropriate class [45]. Finally, all LULC types were identified and classified into five land cover classes (agriculture, grassland, forest, urban and water bodies), named based on the nomenclature developed by the Chinese Academy of Science (CAS) [46]. Then, based on the Normalized Difference Vegetation Index (NDVI), we reassessed the misclassified pixels from the initial classification and reclassified accordingly. The isolated pixels in each classified image were removed by applying the majority spatial filters [47]. The classified maps showed satisfactory results for the overall classification accuracies of 87% and 91% and Kappa coefficients of 0.86 and 0.89, (1 represents the best classification [48]) for the Yanhe and Guangyuan catchments, respectively. The accuracy assessment results were satisfactory following the standard suggested by Anderson [45] for Landsat image classification.

Land Use Land Cover Change and Vegetation Analysis
According to Lu et al. [49], assessing the LULC change requires the analysis of satellite images representing the same area taken at different times. Therefore, to allow comparison and detect any changes in our study areas across selected years, the classified Landsat images were analyzed. We quantified all classified LULC maps by using the post classification change extraction method [50,51]. Therefore, the statistical information and area (in km 2 ) of each LULC class from the classified maps were computed to generate a change matrix. The similar classes structure used during the classification of Landsat images simplified the computation. Furthermore, the spatial analysis tools in ArcGIS 10.3 were used to better visualize the observed trends in land cover change in the two study catchments.
Vegetation plays a crucial role in flash flood control and regulation because areas with dense vegetation and forest diminish surface runoff [52][53][54][55]. In this study, the vegetation distribution and status in the study catchments in various historical periods was analyzed using the normalized difference vegetation index (NDVI). We calculated the NDVI values from the surface reflectance Landsat images as shown in Equation (1) [56]. The NDVI maps were compared to assess vegetation changes in each year. Moreover, the NDVI maps were used to validate the classified LULC maps.
where ρNIR is the near infrared band and ρRED is the red band.

Hydrological Modeling Using HEC-GeoHMS and HEC-HMS
The Hydrologic Engineering Center Hydrologic Modeling System (HEC-HMS) is an event-based hydrological model and was developed to simulate the hydrological processes of various types of catchments in different areas of research, including watershed runoff studies, flood damage studies through the analysis of hydrography, floodplain management etc., [57]. The model consists of four modules for event simulation including basin and meteorological modules, time series data, and control specifications. The basin model includes processes such as basin and river routing while the meteorological model includes rainfall data. The control specifications are used for timing processes of the model. After integration of all the required data, the HEC-HMS applies all the models to simulate the hydrological process. The outcome of the HEC-HMS model is the hydrograph at the catchment outlet. More details on the model functionality of the HEC-HMS model are provided in the user and technical references [57,58]. Sanyal et al. [59] suggested the use of the HEC-HMS model in a data scarce situation as the main inputs data to run the model can be obtained through remote Water 2020, 12, 1347 6 of 22 sensing means. Considering that there are no sufficient ground data in our two study catchments, we selected HEC-HMS version 4.2.1 to quantify the impact of land cover changes on flash flood peak discharge, based on the application of suitable loss, transform and routing methods. All the main input data required to run the HEC-HMS model are prepared in the HEC-GeoHMS model. The HEC-GeoHMS model is an extension of Arc GIS software. The units of the HEC-HMS model are provided in detail below.

The Basin Model
The main input data required to run the HEC-HMS Model were the Digital Elevation Model (DEM), land cover maps, soil map and rainfall data [58] (Figure 2). The basin model incorporates the data relevant to the physical attributes of the HEC-HMS model. We used HEC-GeoHMS, an extension in the ArcGIS environment, to delineate the study catchments from DEM. Through a number of processing steps, known as the terrain processing to the DEM, we derived the streams network and subdivide our study catchments into interconnected sub-basins. In addition, other catchment characteristics such as elevation, slope, stream length and the area of each sub-catchment were obtained. The hydrologic elements of the basin model for the two catchments are presented in Figure 3.

The HEC-HMS Model Control Specifications
The model control specifications incorporate data related to the timing of the HEC-HMS model such as the time steps, start and end date, and simulation time adopted in the model. In the present study, the time steps were set to 1 h and were selected based on the computed shortest lag times of the small sub-basin in the HEC-HMS model [58].

Curve Number Loss Method
To calculate the rainfall losses occurring in the catchment, we selected the soil conservation

Routing Method
The routing method accounts for the route the runoff through the channel reaches to the outlet of the catchment and generates a flow hydrograph [64]. In this study, we used the Muskingum-Cunge routing method, which uses two parameters K and X. The travel time K simulates a delay in runoff as it passes through the channel whereas the degree of storage values X (assumed to be 0.2) is used to regulate the influence of inflow and outflow hydrograph [58].

Meteorological Model
The rainfall data are an essential input during simulation of flash flood event using hydrological models, as high intensity rainfall falling for a few hours may cause flooding [60]. In the HEC-HMS model, the rainfall data are required to set up the meteorological model and are one of the main inputs to run the model. In this study, we examined rainfall data at hourly temporal resolution during 2002 to 2011 and selected the largest rainstorm events that caused flash floods in each study catchment. For the Yanhe catchment, we selected intense rainfall events from the 18-22 August 2010 for the model simulation. Similarly, for the Guangyuan catchment, the rainfall events that took place on 22-25 August 2010 were given as input in the model. The distribution of these rainfall events in the catchment was calculated using the kriging interpolation method [61] based on the rainfall values collected at rain gauging stations located in each studied catchment.

The HEC-HMS Model Control Specifications
The model control specifications incorporate data related to the timing of the HEC-HMS model such as the time steps, start and end date, and simulation time adopted in the model. In the present study, the time steps were set to 1 h and were selected based on the computed shortest lag times of the small sub-basin in the HEC-HMS model [58].

Curve Number Loss Method
To calculate the rainfall losses occurring in the catchment, we selected the soil conservation service (SCS) curve number (CN) loss method [62]. The SCS-CN method is extensively used in ungauged watershed because of minimal data requirements [63]. It uses the CN parameter, which reflects land features such as LULC and soil type characteristics of the two study catchments. It has been successful in simulating hydrological processes in many regions [59,64]. The SCS-CN method uses CN values for each sub-basin and the initial abstraction values. The low CN value signifies higher infiltration losses while the higher CN correspond to less infiltration losses. Thus, the ground surface with a higher CN value is more likely to generate higher runoff yield and cause floods. Equation (2) presents the calculation of loss using SCS-CN method [62].
where Q is the accumulated direct runoff depth, P is accumulated storm rainfall, Ia is the initial abstraction (initial loss) and is given by 0.2S, and S is the potential maximum retention and is a measure of the capability of a catchment to absorb and retain storm rainfall. The potential maximum retention (S) is related to the curve number (CN) by the following equation (Equation (3)): where CN is the soil conservation service (SCS) curve number (CN), the SCS-CN is an index that illustrates the combination of soil hydrologic group, LULC classes, and antecedent soil moisture conditions. We generated the parameter CN by combining the hydrologic soil groups (HSG) and land use/cover maps of each studied catchment in different periods. The HSG data were prepared using the soil data of each studied that contained all the information related to soil texture and soil types of each of the study catchment. The HSG categorizes soil classes into four groups (A, B, C and D), where soils in group A has high infiltration potential and soils in group D has low infiltration rate [58]. We referred to the CN table defined by the USDA [58] to generate the CN maps of each catchment. The cover type was determined by comparing cover type as named in the CN table to the LULC type in each studied catchment ( Table 1). The obtained CN maps show the physical parameters in each catchment and were used as input during model simulation by HEC-HMS. In this study, we used the SCS unit hydrograph method to convert the excess rainfall into a flow hydrograph at the catchment outlet [65]. The basin lag time and time of concentration (Tc) are the main input parameters. The basin lag time was calculated for each sub-basin using the curve number values and sub-basin slope length, whereas the Tc was estimated based on sub-basin characteristics (e.g., topography) and reach length in the HEC-GeoHMS software environment.

Routing Method
The routing method accounts for the route the runoff through the channel reaches to the outlet of the catchment and generates a flow hydrograph [64]. In this study, we used the Muskingum-Cunge routing method, which uses two parameters K and X. The travel time K simulates a delay in runoff as it passes through the channel whereas the degree of storage values X (assumed to be 0.2) is used to regulate the influence of inflow and outflow hydrograph [58].

Simulated Land Use Land Cover Change Scenarios
To evaluate and understand the effect of the two conservation policies (NFCP and GTGP) on flood peak discharge, we set up four scenarios of LULC for each studied catchment. These LULC scenarios represent the land cover conditions from 1990 to 2016 for the Yanhe catchment and from 1990 to 2017 for the Guangyuan catchment, respectively. For example, the first scenario was performed using the 1990 LULC map and was considered as the baseline scenario (S1). Similarly, the second (S2), the third (S3) and the fourth (S4) scenarios were established based on other historical LULC maps [66]. The land cover change rate and the change of land cover area (in km 2 ) under these four scenarios in the two catchments are presented in Figure 4. For each catchment, we conducted four hydrological simulations (S1, S2, S3 and S4) using the HEC-HMS hydrological model while keeping constant other model input parameters that describe the catchment characteristics and boundary conditions, hence, to observe and compare the changes in flood peak discharge in each historical time period due to the isolated impact of LULC changes alone [67,68].
flood peak discharge, we set up four scenarios of LULC for each studied catchment. These LULC scenarios represent the land cover conditions from 1990 to 2016 for the Yanhe catchment and from 1990 to 2017 for the Guangyuan catchment, respectively. For example, the first scenario was performed using the 1990 LULC map and was considered as the baseline scenario (S1). Similarly, the second (S2), the third (S3) and the fourth (S4) scenarios were established based on other historical LULC maps [66]. The land cover change rate and the change of land cover area (in km 2 ) under these four scenarios in the two catchments are presented in Figure 4. For each catchment, we conducted four hydrological simulations (S1, S2, S3 and S4) using the HEC-HMS hydrological model while keeping constant other model input parameters that describe the catchment characteristics and boundary conditions, hence, to observe and compare the changes in flood peak discharge in each historical time period due to the isolated impact of LULC changes alone [67,68]. Following the argument of Wilk and Hughes [69], we assessed the impact of changed vegetation cover on flood peak discharge at the sub-basin scale in each study catchment. This was done by comparing the model outputs of the baseline scenario S1 to that of scenario S4 in the four sub-basins that have undergone significant changes in land cover and examine any change in peak discharge during our study period. Moreover, we analyzed the spatial NDVI changes in each sub-basin and compared them to changes in peak discharge.

Land Use Land Cover Classification and Change Detection Analysis
The LULC classified maps for both the Yanhe and Guangyuan catchments are presented in Figures 5 and 6, and the corresponding detailed changes in land cover types (in %) were analyzed and quantified as compiled in Table 2 and Table 3, respectively. Both studied sites have undergone extensive land cover changes due to NFCP and GTGP conservation programs. Following the argument of Wilk and Hughes [69], we assessed the impact of changed vegetation cover on flood peak discharge at the sub-basin scale in each study catchment. This was done by comparing the model outputs of the baseline scenario S1 to that of scenario S4 in the four sub-basins that have undergone significant changes in land cover and examine any change in peak discharge during our study period. Moreover, we analyzed the spatial NDVI changes in each sub-basin and compared them to changes in peak discharge.

Land Use Land Cover Classification and Change Detection Analysis
The LULC classified maps for both the Yanhe and Guangyuan catchments are presented in Figures 5 and 6, and the corresponding detailed changes in land cover types (in %) were analyzed and quantified as compiled in Tables 2 and 3, respectively. Both studied sites have undergone extensive land cover changes due to NFCP and GTGP conservation programs.
The LULC changes in the Yanhe catchment were investigated from 1990 to 2016 (Table 2 and Figure 5). The results showed that agricultural land was the main type of land use and occupied proportions of 58.39% and 27.68% in 1990 and 2016, respectively. We observe a continuous decrease in agricultural land since the year 1990. Grassland was the second main land cover type in the Yanhe catchment and increases continuously from 31.90% in 1990 to 43.91% in 2016. During the same period, the forest area exhibited an obvious increase from 9.37% to 26.48%. The forest cover is found in the middle area of the catchment and expanded toward downstream areas. This demonstrated that the Yanhe catchment experienced a fast reforestation program over the last 26 years. The urban land has the least coverage and increased from 0.20% to 1.80% during the entire study period. Similarly, the area of water bodies is very small (0.14%) and exhibited a slight changes from 1990 to 2016. The spatial analysis reveals that during 1990-2016, the diminishing area of agricultural land (2343.49 km 2 ) is mainly attributed to the increase in grassland and forest cover in this study catchment.  The LULC changes in the Yanhe catchment were investigated from 1990 to 2016 (Table 2 and Figure 5). The results showed that agricultural land was the main type of land use and occupied proportions of 58.39% and 27.68% in 1990 and 2016, respectively. We observe a continuous decrease in agricultural land since the year 1990. Grassland was the second main land cover type in the Yanhe catchment and increases continuously from 31.90% in 1990 to 43.91% in 2016. During the same period, the forest area exhibited an obvious increase from 9.37% to 26.48%. The forest cover is found in the middle area of the catchment and expanded toward downstream areas. This demonstrated that the Yanhe catchment experienced a fast reforestation program over the last 26 years. The urban land has the least coverage and increased from 0.20% to 1.80% during the entire study period. Similarly, the area of water bodies is very small (0.14%) and exhibited a slight changes from 1990 to 2016. The spatial  For the Guangyuan catchment, in the year 1990, agricultural land was the most dominant LULC type, occupying 56.55% of the total area, followed by forest cover type 41.98%, urban area 1.06% and water bodies 0.41% (Table 3 and Figure 6). During the whole study period 1990-2017, the agricultural land cover type declined the most from 56.55% to 32.81% while the other classes such as forest cover and urban area increased. There is a steady increase in forest cover (16.49%) in the first period of our For the Guangyuan catchment, in the year 1990, agricultural land was the most dominant LULC type, occupying 56.55% of the total area, followed by forest cover type 41.98%, urban area 1.06% and water bodies 0.41% (Table 3 and Figure 6). During the whole study period 1990-2017, the agricultural land cover type declined the most from 56.55 to 32.81% while the other classes such as forest cover and urban area increased. There is a steady increase in forest cover (16.49%) in the first period of our analysis (1990-2000) and a slight decrease (−0.8) between 2000 and 2017, probably due to the rapid urbanization (5.61%) in the later period. In general, from 1990 to 2017, the forest cover exhibited an overall increase of 15.7%. Meanwhile, the urban area expanded eastward of the catchment along the floodplain from 1.06 to 9.16%. The spatial patterns of LULC indicated that mainly agricultural land in the floodplain regions was converted to built-up areas in the surroundings of the 1990 urban center boundary of this study catchment. Further analysis of land cover and slope maps revealed that agricultural land was mainly transformed to forest on steep slopes areas (>25%).
Some other studies [12,36] observed the same trend related to the conversion of agricultural land into forest cover and grassland in many mountainous regions of China, including the location of the two studied catchments. The increase in forest coverage is likely attributed to the GTGP and NFCP conservation policies implemented in the late 1990s as a way of reducing soil erosion and floods in many regions in China [15,70].

Vegetation Analysis
The NDVI classification maps show spatiotemporal changes in vegetation cover for the two study catchments (Figures 7 and 8). The analysis of the NDVI classification maps indicates that the vegetation cover area experienced an increasing trend in both catchments. In the Yanhe catchment, from 1990 to 2016, the NDVI increased by 0.14 (from 0.78 to 0.92). Spatially, the lowest NDVI values were observed in the upstream areas whereas the increasing trends in NDVI were observed in the midstream towards downstream areas of the catchment (Figure 7). These regions with higher NDVI values have been the most affected by increases in forest and grassland cover.
The NDVI classification maps show spatiotemporal changes in vegetation cover for the two study catchments (Figures 7 and 8). The analysis of the NDVI classification maps indicates that the vegetation cover area experienced an increasing trend in both catchments. In the Yanhe catchment, from 1990 to 2016, the NDVI increased by 0.14 (from 0.78 to 0.92). Spatially, the lowest NDVI values were observed in the upstream areas whereas the increasing trends in NDVI were observed in the midstream towards downstream areas of the catchment (Figure 7). These regions with higher NDVI values have been the most affected by increases in forest and grassland cover. For the Guangyuan catchment, the mean NDVI value increased from 0.60 in 1990 to 0.85 in 2017. The highest increment in the NDVI value was observed in upstream hilly areas of the study catchment ( Figure 8). The spatial analysis of LULC maps shows the increase in forest coverage on The increasing trend in NDVI indicates that both catchments were affected by various reforestation and afforestation programs (e.g., NFCP and GTGP) introduced since the 1990s. These conservation projects reflect the effort made to restore vegetation cover in various mountainous regions in China, including the two study catchments. Overall, our results are consistent with other research findings [20] conducted in these regions where the increasing trends in NDVI were noted after the implementation of the Grain to Green Program.

Impact of Reforestation Driven LULC Changes on Flood Peak Discharge
The HEC-HMS model was applied in each catchment for four different land cover scenarios (S1, S2, S3 and S4) using the same meteorological input. The simulated flash flood hydrographs at the catchment outlet decreased between the baseline scenario S1 and fourth scenario S4 (Table 4 and For the Guangyuan catchment, the mean NDVI value increased from 0.60 in 1990 to 0.85 in 2017. The highest increment in the NDVI value was observed in upstream hilly areas of the study catchment ( Figure 8). The spatial analysis of LULC maps shows the increase in forest coverage on steep slopes (>25%). On the other hand, we observed a sharp decline in NDVI values in the low elevation areas near the catchment outlet and alongside the river floodplain. The analysis of the LULC maps revealed the highest expansion in urban land (4.33%) in these regions from 1990 to 2017.
The increasing trend in NDVI indicates that both catchments were affected by various reforestation and afforestation programs (e.g., NFCP and GTGP) introduced since the 1990s. These conservation projects reflect the effort made to restore vegetation cover in various mountainous regions in China, including the two study catchments. Overall, our results are consistent with other research findings [20] conducted in these regions where the increasing trends in NDVI were noted after the implementation of the Grain to Green Program.

Impact of Reforestation Driven LULC Changes on Flood Peak Discharge
The HEC-HMS model was applied in each catchment for four different land cover scenarios (S1, S2, S3 and S4) using the same meteorological input. The simulated flash flood hydrographs at the catchment outlet decreased between the baseline scenario S1 and fourth scenario S4 (Table 4 and Figure 9). The magnitude of changing flood peak discharge was different in each catchment. For the Yanhe catchment, we observed a larger decreasing trend in flash flood peak discharge from 1990 to 2016 (Figure 9a). A maximum decrease in flood peak discharge of 14.35% between scenario S4 (2016 LULC) and baseline scenario S1 (1990 LULC) was recorded. Similar decreasing trends in peak discharge were noted in scenario S2 (2002 LULC) and scenarios S3 (2010 LULC) as shown in Table 4. These findings show that reforestation induced land cover changes probably led to a continuous decrease in flood peak discharge since the 1990s. From 1990 to 2016, the Yanhe catchment had undergone an expansion in forest and grassland cover and a decrease in agricultural land (see Section 3.1). Additionally, this change in land cover was reflected in changing NDVI values, which increased by 0.14 during the same period. As discussed by other authors [16,71], the vegetation cover intercepts rainfall and increases surface roughness (Manning's coefficients), which result in higher soil infiltration capacity and a decrease in surface runoff and streamflow.
In the Guangyuan catchment, the HEC-HMS simulated results showed that between 1990 to 2017, the flood peak discharge decreases up to −6.2% in scenario S4 (2017 LULC) compared to the baseline scenario S1 (1990 LULC) at the catchment outlet ( Figure 9b). However, scenario S2 (2000 LULC) was the most effective in reducing flood peak discharge (up to −7.5%). In general, all three scenarios (S2, S3, and S4) showed a decreasing tendency in peak discharge compared to the baseline scenario (S1) as presented in Table 4. Many possible reasons explain the observed unsteady changes in flash flood peak discharge noted from 1990 to 2016. First, for the baseline scenario S1 (LULC1990), the analysis of the land cover map showed that in the year 1990, agricultural land was the main LULC type, occupying more than 56.55% of the total area of the study catchment. Previous studies [72] argued that the expansion of agricultural land led to the increase in streamflow, which might explain the cause of the highest flood peak discharge observed during the baseline scenario (S1). Second, the conversion of agricultural land to forest and grasslands between 1990 and 2000 had more effect on the larger percentage reduction in flood peak discharge registered in scenario S2 (2000LULC). Third, we observed a high rise in built-up areas (4.3%) from 2000 to 2008 on the fringe of the 1990 city boundary mainly along the floodplain. This probably increased the surface runoff, which led to little increase in peak discharge in scenario S3 (2008LULC). the cause of the highest flood peak discharge observed during the baseline scenario (S1). Second, the conversion of agricultural land to forest and grasslands between 1990 and 2000 had more effect on the larger percentage reduction in flood peak discharge registered in scenario S2 (2000LULC). Third, we observed a high rise in built-up areas (4.3%) from 2000 to 2008 on the fringe of the 1990 city boundary mainly along the floodplain. This probably increased the surface runoff, which led to little increase in peak discharge in scenario S3 (2008LULC). Overall, our results demonstrated that the changes in land cover driven by afforestation activities (e.g., GTGP and NFCP) across the two mountainous catchments were probably the main cause of decreasing flood peak discharge during the 1990-2016/2017. Moreover, the difference in the catchment area and location of changed land cover might have caused the observed variation in peak discharge at the outlet of each catchment. As discussed by Jansson et al. [22], the location of vegetation cover in the watershed often contributes to the amount of changes in streamflow at the outlet. Our findings agree with other previous studies [73,74] conducted in various watersheds in China, where the increase in woodland resulted in a lower runoff and water discharge. Elsewhere, authors reported the reduction in forest cover driven by human activities as the main cause of increases in surface runoff and flood peak discharge [75][76][77]. Table 4. Modeling comparison of different LULC changes scenarios for the Yanhe and Guangyuan watersheds. Overall, our results demonstrated that the changes in land cover driven by afforestation activities (e.g., GTGP and NFCP) across the two mountainous catchments were probably the main cause of decreasing flood peak discharge during the 1990-2016/2017. Moreover, the difference in the catchment area and location of changed land cover might have caused the observed variation in peak discharge at the outlet of each catchment. As discussed by Jansson et al. [22], the location of vegetation cover in the watershed often contributes to the amount of changes in streamflow at the outlet. Our findings agree with other previous studies [73,74] conducted in various watersheds in China, where the increase in woodland resulted in a lower runoff and water discharge. Elsewhere, authors reported the reduction in forest cover driven by human activities as the main cause of increases in surface runoff and flood peak discharge [75][76][77].

Spatial Variability of Flood Peak Discharge at the Sub-Basin Scale
The spatial distributions of changes in peak discharge at the sub-basin scale in the Yanhe and Guangyuan catchments between S1 and S4 are shown in Figures 10 and 11, respectively. These figures show that the extent of peak discharge decrease differed among sub-basins. For the Yanhe catchment, from 1990 to 2016, the highest decreases in peak discharge were observed in the midstream and downstream sub-basins with a maximum percentage change of approximately 17.7% ( Figure 10).
The spatial analysis of land cover changes showed that these sub-basins were associated with the replacement of agricultural land by forest and grassland cover. Moreover, these sub-basins had higher NDVI, with a maximum value of 0.95 in 2016, which confirms the presence of higher vegetation coverage. On the other hand, the lowest decrease in peak discharge, were noted in the northern portion of sub-basins with a minimum percentage change of 7%. The spatial distribution pattern of land cover type shows that these sub-basins largely match with agriculture activities and experienced a little increase in grassland and forest coverage during our study period. Furthermore, we selected individual sub-basins with significant changes in various land cover type and assess their hydrological response during 1990-2016. For example, the expansion of urban land in sub-basin W850 resulted in a slight decrease in peak discharge (4.9%). In contrast, the increase in forest and grassland in sub-basin W960, resulted in the highest decrease in peak discharge up to 14% during the same study period (Figure 12a). Overall, all sub-basins showed a decreasing trend at various degrees in peak discharge from 1990 to 2016 in this study catchment. The observed variation in peak discharge among sub-basins might have been caused by the higher rate conversion of agricultural land to forest and grassland cover between the year 1990 and 2017 in these sub-basins. In the Guangyuan catchment, the highest decrease in peak discharge (up −19.06%) occurred in the upstream (southern and eastern) sub-basins during 1990-2017 ( Figure 11). Conversely, the downstream sub-basins showed an increasing tendency in peak discharge (up to 39.27%). The spatial patterns of land cover and the analysis of NDVI maps indicate a significant increase in urban land in downstream sub-basins and an expansion of forest cover in upstream sub-basins. Therefore, the urbanized sub-basins were probably the main source of surface runoff that increased the peak discharge in downstream areas while the forest coverage reduced the peak discharge in upper areas. Moreover, the analysis of some individual sub-basins provided further understanding on the effect of sub-basins land cover change on hydrological responses. For example, from 1990 to 2017 in the urbanized sub-basin W900, the streamflow changed from 73.60 to 102.5 m 3 /s. On the contrary, in the forested sub-basin W1180, the streamflow decreased from 77.4 to 64.8 m 3 /s during the same period (Figure 12b). These findings indicate that the changes in land cover driven by GTGP and NFCP certainly affected more the sub-basins hydrological response behavior than the whole catchment. In the Guangyuan catchment, the highest decrease in peak discharge (up −19.06%) occurred in the upstream (southern and eastern) sub-basins during 1990-2017 ( Figure 11). Conversely, the downstream sub-basins showed an increasing tendency in peak discharge (up to 39.27%). The spatial patterns of land cover and the analysis of NDVI maps indicate a significant increase in urban land in downstream sub-basins and an expansion of forest cover in upstream sub-basins. Therefore, the urbanized sub-basins were probably the main source of surface runoff that increased the peak discharge in downstream areas while the forest coverage reduced the peak discharge in upper areas. Moreover, the analysis of some individual sub-basins provided further understanding on the effect of sub-basins land cover change on hydrological responses. For example, from 1990 to 2017 in the urbanized sub-basin W900, the streamflow changed from 73.60 to 102.5 m 3 /s. On the contrary, in the forested sub-basin W1180, the streamflow decreased from 77.4 to 64.8 m 3 /s during the same period (Figure 12b). These findings indicate that the changes in land cover driven by GTGP and NFCP certainly affected more the sub-basins hydrological response behavior than the whole catchment. The analysis of hydrological behavior at the sub-basin scale in both catchments was consistent with the magnitude of change of NDVI values (see Section 3.2). We observed no significant correlation (p > 0.10 for both catchments) between changes in NDVI and flood peak discharge at the sub-basin scale. For example, the increase in NDVI by 19% led to a decrease in peak discharge by 8% in one sub-basin of the Guangyuan catchment. In general, sub-basins that showed higher decreasing trends in flood peak discharge were characterized by larger increases in NDVI values than sub-basins that exhibited low or non-decrease in peak discharge. Moreover, the sensitivity analysis exposed a positive correlation (r 2 = 0.99) between changes in land cover (increase/decrease in CN) and simulated peak discharge. For example, the increase in CN value by 5% in one sub-basin of the Yanhe catchment has resulted in an increase in peak discharge by 4%. These results indicate that the risk of flash floods in the catchment increases with changing land cover types towards a more urbanized environment. Furthermore, these analyses helped us to better understand the impact of land cover dynamics on catchment hydrological response. Our findings are in line with the study by Liu and Shi [11] who attributed the increase in surface runoff in various urban areas in China to the conversion of floodplains into constructed lands. Moreover, other numerous studies [78,79] discussed the benefit of increases in forest cover in the small catchment and concluded that an increase in vegetation coverage results in a reduction in streamflow. The analysis of hydrological behavior at the sub-basin scale in both catchments was consistent with the magnitude of change of NDVI values (see Section 3.2). We observed no significant correlation (p > 0.10 for both catchments) between changes in NDVI and flood peak discharge at the sub-basin scale. For example, the increase in NDVI by 19% led to a decrease in peak discharge by 8% in one sub-basin of the Guangyuan catchment. In general, sub-basins that showed higher decreasing trends in flood peak discharge were characterized by larger increases in NDVI values than sub-basins that exhibited low or non-decrease in peak discharge. Moreover, the sensitivity analysis exposed a positive correlation (r 2 = 0.99) between changes in land cover (increase/decrease in CN) and simulated peak discharge. For example, the increase in CN value by 5% in one sub-basin of the Yanhe catchment has resulted in an increase in peak discharge by 4%. These results indicate that the risk of flash floods in the catchment increases with changing land cover types towards a more urbanized environment. Furthermore, these analyses helped us to better understand the impact of land cover dynamics on catchment hydrological response. Our findings are in line with the study by Liu and Shi [11] who attributed the increase in surface runoff in various urban areas in China to the conversion of floodplains into constructed lands. Moreover, other numerous studies [78,79] discussed the benefit of increases in forest cover in the small catchment and concluded that an increase in vegetation coverage results in a reduction in streamflow.

Limitations
This study was conducted in mountainous catchments with limited ground discharge data, and due to their inaccessible locations, detailed field surveys for estimation of hydrological parameters were not conducted. Additionally, we used freely available data with medium spatial resolution (e.g., soil data). This may have introduced some uncertainties in the modeling. Nevertheless, by following standard guidance on model input data preparation and processing, these errors have been minimized. Hence, our study provides a reference for better understanding the importance of two main conservation policies introduced in some regions of China to mitigate flash flood and how such program could be useful in other regions. The present study did not consider the impact of climate change due to insufficient long historical climate data, thus, care should be taken while interpreting our results. Some other studies [4,[80][81][82] also reported reliable results while using un-calibrated models.

Summary and Conclusion
The present study investigated the effect of GTGP and NFCP induced land cover change on flood peak discharge in two mountainous catchments in China since the 1990s. Our results indicate the increase in forest coverage by about 18% in Yanhe and 16% in Guangyuan at the expense of agricultural land, reflected in changes in NDVI value and LULC. Similarly, we observed the increase in urban land in some sub-basins in both catchments. These changes in LULC resulted in different hydrological responses in the two catchments. The HEC-HMS simulation results under four different LULC change scenarios showed a decreasing trend in flash flood peak discharge by ~14% in Yanhe and ~6% in Guangyuan during the whole study period (1990-2016/2017). Further to this, the spatial analysis of catchment hydrological response to land cover changes at the sub-basin scale was much more significant. The forested sub-basins demonstrated a high decreasing rate in peak discharge ranging from 17% in Yanhe to 19% in Guangyuan while few urbanized sub-basins showed an increasing trend in peak discharge up to 39%. It appears that the expansion of forest and grassland caused by the GTGP and NFCP policies certainly had great potential in reducing the flood peak discharge in the two studied catchments. Despite providing reference for better understanding the importance of implementing sustainable land cover policies (e.g., reforestation) in mountainous regions, these findings provide useful information for decision makers to establish proper flood control measures. Considering that, the present study only explored the influence of conservation programs (e.g., GTGP) on mountainous catchment hydrological responses, future research shall explore the effect of other factors such as climate conditions and socioeconomic variables.

Limitations
This study was conducted in mountainous catchments with limited ground discharge data, and due to their inaccessible locations, detailed field surveys for estimation of hydrological parameters were not conducted. Additionally, we used freely available data with medium spatial resolution (e.g., soil data). This may have introduced some uncertainties in the modeling. Nevertheless, by following standard guidance on model input data preparation and processing, these errors have been minimized. Hence, our study provides a reference for better understanding the importance of two main conservation policies introduced in some regions of China to mitigate flash flood and how such program could be useful in other regions. The present study did not consider the impact of climate change due to insufficient long historical climate data, thus, care should be taken while interpreting our results. Some other studies [4,[80][81][82] also reported reliable results while using un-calibrated models.

Summary and Conclusions
The present study investigated the effect of GTGP and NFCP induced land cover change on flood peak discharge in two mountainous catchments in China since the 1990s. Our results indicate the increase in forest coverage by about 18% in Yanhe and 16% in Guangyuan at the expense of agricultural land, reflected in changes in NDVI value and LULC. Similarly, we observed the increase in urban land in some sub-basins in both catchments. These changes in LULC resulted in different hydrological responses in the two catchments. The HEC-HMS simulation results under four different LULC change scenarios showed a decreasing trend in flash flood peak discharge by~14% in Yanhe and~6% in Guangyuan during the whole study period (1990-2016/2017). Further to this, the spatial analysis of catchment hydrological response to land cover changes at the sub-basin scale was much more significant. The forested sub-basins demonstrated a high decreasing rate in peak discharge ranging from 17% in Yanhe to 19% in Guangyuan while few urbanized sub-basins showed an increasing trend in peak discharge up to 39%. It appears that the expansion of forest and grassland caused by the GTGP and NFCP policies certainly had great potential in reducing the flood peak discharge in the two studied catchments. Despite providing reference for better understanding the importance of implementing sustainable land cover policies (e.g., reforestation) in mountainous regions, these findings provide useful information for decision makers to establish proper flood control measures. Considering that, the present study only explored the influence of conservation programs (e.g., GTGP) on mountainous catchment hydrological responses, future research shall explore the effect of other factors such as climate conditions and socioeconomic variables.