“…However, the sub-watersheds were still too large (total surface area of 1 st , 2 nd , and 3 rd sub-watersheds were 1,413.27 ha, 936.54 ha, and 2,948.4 ha, respectively) to be properly run with the GeoWEPP program as our attempts failed several times. Therefore, as applied in previous studies [11,[20][21], in order to both run the program easily and obtain detailed data on soil loss and sediment yield in the GCW, all three sub-watersheds were further sub-divided into smaller hydrological units (SHUs) using the hydrologic In this subdivision approach, we believe that running each SHU with GeoWEPP resulted in better and more detailed soil loss and sediment yield data [22] since we could enter detailed soil, land cover, and management data for each SHU specifically compared to entering one large data set (e.g., climate) for the whole study area. In addition, this method allowed us to better evaluate the areas experiencing soil loss and sediment amounts as we can easily see some basic characteristics (land use, soil properties, slope, elevation, and bedrock type, etc.)…”