Development of GIS Interface Tool for GAMES Model and Its Application to an Agricultural Watershed in Southern Ontario

Soil erosion is an important economic and environmental concern throughout the world. In order to assess soil erosion risk and conserve soil and water resources, soil erosion modeling at the watershed scale is imperative. The Guelph model for evaluating effects of Agricultural Management System on Erosion and Sedimentation (GAMES) is tailor-made for such applications; it, however, requires a significant amount of spatial information which needs to be pre-processed using a Geographic Information System (GIS). The GAMES model currently lacks any such automated tools. As such, the GAMES was loosely coupled to a GIS interface to manage the large spatial input data and to produce efficient cartographic representations of model output results. The developed interface tool was tested to simulate the Kettle Creek paired watershed in Southern Ontario, Canada. Result demonstrated that the GIS-assisted procedure increased the ability of the GAMES model in simulating such a spatially varied watershed and made the process more efficient and user-friendly. Furthermore, the quality of reporting and displaying resultant spatial output was also significantly improved. The developed GAMES interface could be applied to any watershed, and the enhancement could be used to assess soil erosion risk and conserve soil and water resources in an effective way.


Introduction
Prior to the use of the computer-based hydrologic model and Geographic Information System (GIS) in hydrology and hydraulic projects, practitioners assembled a number of maps, imageries, stream network, and other relevant data from field surveys to perform manual preliminary hydrologic analyses such as catchment delineation [1] [2].These manual methods were very laborious, time-consuming and inefficient [3].The advancement in GIS technology, digital spatial data, and efficient computers has precluded the need for many of these labor-intensive techniques [4].Developments in the field of GIS helped users to organize all forms of geographically referenced data as layers and various GIS software/tools such as the ArcGIS [5], enabled hydrologists to pre-process, view and analyze project-specific digital elevation model (DEM), soil, land-use, and hydro-meteorological data spatially [6].Besides, several hydrologic functions are in-built to perform some of the preliminary investigations which have reduced the level of approximation and subjectivity involved [7].Use of GIS technology in any hydrology project significantly improved the quality, accuracy, and timeliness of the final product [8].Furthermore, GIS greatly simplifies the model set up and improves model performance [9].Additionally, the GIS advancements enabled users to handle a large amount of data in a short time frame.
Another critical component of GIS is an enhancement of the visual understanding of the spatial inputs and output results by managers, decision-makers and the public-at-large [10].Changing the inputs of land use or introducing best management practices (BMPs) to compare the subsequent output can be very helpful to identify the problem and ease the decision-making process.Similarly, specialized interfaces when linked to a GIS application have increased the popularity and usability of spatially distributed hydrological models and led them to greater use and wide acceptance [11].
Developments of GIS interface for the hydraulic/hydrologic model started in 1975 when the HEC-1 model was integrated using a grid-based model.The resulted model was named as HEC-SAM.In this approach, the GIS software was simply used as a database to feed input data to the model [12].There are several different approaches to integrating GIS with simulation models such as; the embedding method, loose coupling and tight coupling methods [13] [14] [15] [16].The loose coupling approach is the simplest as the GIS software and the hydrologic model exchange files, read some of its input data files from geo-database and produce resultant output in a format that allows processing and displays the results within the GIS software.This is a standard approach since it requires little or no modifications to model software [17] [18].From the successes of the past two decades, hydrological modelers have recognized the benefits of integration of GIS software and hydrological models [19] [20].The advancement in GIS technology and improved techniques for the integration of GIS with hydrological modeling has made hydrological models more efficient and user-friendly [21] [22] [23] [33].Open Journal of Civil Engineering The Guelph model for evaluating Effects of Agricultural Management Systems on Erosion and Sedimentation (GAMES) [24], developed at University of Guelph, Guelph, Ontario, Canada, is seasonal event-based soil erosion and sediment yield simulation model for analyzing watershed management practices and effects on nonpoint source pollution changes [25].However, it was not coupled with any GIS tool.As such, model inputs at user-specified spatial units were not automated which made the model difficult and inefficient to use in areas with high-resolution spatial maps.Furthermore, analyses of model outputs which could greatly be enhanced with the use of any GIS platform was hampered.The GAMES model thus needed a much-needed upgrade in this respect.
In this research, a GIS interface modeling tool has been developed which involved two distinct phases.This is the first attempt of this kind specifically to the GAMES model and hence is a novel work.In the first phase, the GAMES model was loosely coupled to a GIS interface by developing a GAMES-GIS model in a GIS system.Then, a GAMES input data generator model was developed in Microsoft Excel using the Visual Basic code in order to convert the processed GIS data in a format compatible with the input format of the GAMES model.
In the second phase, the previous two steps process with many manual-based sub-processes were automated by creating a user-friendly computer interface using C# language.During these processes, no change was made in the source codes of the GAMES model.Hence, the main objective of this study is to present the development and evaluation of the GIS interface to the GAMES model in an agricultural watershed in Southern Ontario, Canada.This has resulted in enhancing the capabilities of GAMES model and making it more efficient and user-friendly.

Games Model
The GAMES model [24] provides estimates of soil loss by water erosion and subsequent delivery of sediment from a field to stream in agricultural water- ( ) where, A is computed site soil loss (t ha −1 yr −1 ), R is a seasonal or event rainfall erosivity factor (MJ mm ha −1 h −1 year −1 ), K is a soil erodibility factor (t h MJ −1 mm −1 ), L is a slope length factor (unitless) that is often combined with S is a slope steep-Open Journal of Civil Engineering ness factor, to yield a unitless terrain factor (LS), C is a unitless crop and land-use management factor, and P is a dimensionless erosion control practice factor.The

Development of the GIS Interface Tool
A GIS interface for the GAMES model was developed in two phases.
In the first phase, an interface was developed using two different platforms (the ArcGIS and Microsoft EXCEL).The interface developed in the ArcGIS platform is termed as the GAMES-GIS model and the one developed in the Visual Basic EXCEL platform is termed as the GAMES Input Data Generator (GAMES-IDG) model [31].
The GAMES-GIS model was developed by stringing the geo-processes to run through specially designed dialogue boxes and all the input data layers; such as study area, soil properties, land use, drainage network, Digital Elevation Model (DEM) through the ArcGIS dialogues (Figure 1).
The GAMES-GIS model uses the capabilities of ArcGIS to extract, process and organize necessary information from the input layers to create the reference ta-  In the second phase, an interactive window interface (Figure 2) was built to integrate all the processes established in the first phase.This interface application was developed using C# language.The user has to input the project title, and description in the interface text fields and the information regarding annual or seasonal or event rainfall erosion index values need to be supplied in the provided textboxes of the window.A drop-down menu is also provided to select the simulation season.The interface divides the whole simulation process into four steps.Firstly, it creates an input file for GAMES model.Secondly, it runs the GAMES model.Thirdly, it processes the GAMES output file.Finally, it converts the tabular output file into spatial layers by using the ArcGIS platform.Appropriate buttons are provided to carry out these steps in the interface.
As can be seen in Figure 2, the "Create Games Input" activates the process of generating an input file for GAMES model.In this step, the input data is loaded and previously developed GAMES_GIS model, and GAMES-IDG model performs Open Journal of Civil Engineering

Application
The developed interface was applied to model the erosion and sediment yield dynamics of the Kettle Creek paired watershed of southwestern Ontario, Canada.

Study Area
The Kettle Creek paired watershed drains two distinct watersheds of Madter (380 ha) and Holtby (355 ha) into their respective municipal drain which eventually flows to the Kettle Creek river, further flowing south through Elgin County into Lake Erie (Figure 3).Soils (Figure 4

Input Data
The

Model Simulation Results
As with any deterministic model, the accuracy of the results depends on the precise   5 shows, for the illustration purposes, some of the output maps such as total erosion and sediment yield rates for spring, summer and fall seasons.
The summary tables (Table 1) and GIS output maps (Figure 5) clearly illustrate the potential soil erosion and sediment yield of both the Matder and Holtby watersheds.It should be noted that the GAMES model has not been quantitatively validated on the field to field basis, as it is not the scope of the paper.
However, qualitative field observations confirmed that the high soil loss and sediment yields areas simulated by the model exhibited significant soil losses and sediment yields.Similarly, areas predicted to have minimal soil losses, and sediment yield does exhibit these characteristics.
Hence, it has been illustrated that the developed GIS interface of the GAMES model has enhanced spatial data pre-processing and input, computational capability, and most importantly the output.Availability of GIS maps such as illustrated in Figure 5, would be crucial to identify the critical source areas of erosion and sediment yield, and that targeted BMPs could be placed over these areas in order to reduce the erosion as erosion is the source-related problem [33].

Conclusion
The study presented a novel development of a GIS-based interface tool for the GAMES model to manage the large spatial input data and to produce efficient cartographic representations of model output results.The developed interface tool was tested to simulate the Kettle Creek paired watershed in Southern Ontario, Canada.Based on the application, it can be concluded that the developed interface greatly enhanced the GAMES model as it allowed capturing, storing, sheds.The model operates on the discretization of the watershed into field sized hydrologic response units (HRUs) with identical characteristics of land-use, soil type, and class of slope.This discretization process allows for the modeling of the basin with homogeneous irregularly shaped HRUs of the watershed.The model also accommodates spatial changes/variations in soil and crop management conditions.The model comprises of two main components, soil loss component, and sediment delivery component.The potential soil erosion component is based on the Universal Soil Loss Equation (USLE) [26]:

K
and C could vary from rainfall event to event and season to season.The input parameters required for the model include land slope, channel slope, soil type, soil erodibility, cropping factors and location and efficiency of sediment detention structures.The sediment delivery from each HRU to downstream HRU and from HRU to the watershed's stream channels is determined by a delivery ratio, calculated based on physical characteristics of the HRU such as roughness and flow length assuming homogeneous characteristics of eroded sediments.The model allows the estimation of potential soil loss and its delivery to streams during various seasons or events, thereby providing useful insights into the annual or seasonal or event variations and their distribution within the watershed.The model, therefore, provides a powerful tool for watershed soil erosion and sediment control planning.It facilitates implementation of efficient and economical ways to apply remedial control strategies and selection of BMPs for different areas of the watershed[27] [28].The model has been evaluated in many Canadian provinces and other countries and has been adopted as an approach to control the nonpoint source of pollution by many governmental and non-governmental agencies[24] [29][30].
bles.Firstly, the Clip function of the ArcGIS is used to clip the soil, land use and drainage network feature files and the Mask function is used for masking the DEM input file for a given specified study area.The second step starts with filling the sinks of the raw DEM file and creates the slope, flow direction, and flow accumulation files.The flow paths are delineated from the flow accumulation raster using the "Greater than Equal" function and the input threshold flow accumulation value.The flow paths are fragmented into sections with a unique ID number which is later used to determine the HRUs.Then, flow direction file is used to define stream links using the Stream Link function and flow paths.Once Open Journal of Civil Engineering

Figure 2 .
Figure 2. A screenshot of the GIS integrated GAMES model interface application.

Figure 3 .
Figure 3. Two watersheds (Madter and Holtby) of the Kettle Creek paired watershed, Ontario, Canada with stream network.Also shown in the inset, is the location of the paired watershed with respect of Lake Huron and Lake Erie.
basic GIS data requirements of the developed interface include a sub-watershed boundary layer (polygon), digital elevation model (raster, Figure 4(a)), soil layer (polygon, Figure 4(b)), land-use layer (polygon, Figure 4(c)), and stream network layer.The annual rainfall erosion index (R) value of 100 MJ mm ha −1 h −1 yr −1 , for the nearest climatic station (London, Ontario, Canada) was obtained from the study by Wall et al. 1983 [32].The seasonal (R) value for spring, summer and fall were calculated to be 28, 48 and 17, respectively, using a procedure suggested by Wall et al. 1983 [32].The GAMES-GIS model generates drainage lines (Figure 4(d)) and divided the watershed into land HRUs (Figure 4(e)) and created different layers of start points, endpoints and drainage lines (Figure 4(f)).

Figure 4 .
Figure 4. (a) DEM, (b) soil and (c) land use, (d) drainage lines (e) delineated HRUs and (f) start and end point of flow paths at both Madter and Holtby sub-watersheds of the Kettle Creek paired watershed.

Figure 5 .
Figure 5. Erosion rate and sediment yield rate for spring, summer and fall seasons.

Table 1 .
Average erosion rate and sediment yield at different seasons of the year at two watersheds of Kettle Creek paired watershed.dataandselection of the parameter values of the model.The GAMES model is based both on parameters which can be physically measured and parameters which must be estimated with the help of literature and experience.The input data values such as cell area, cell slope, length of cell flow path, length of the flow path to the watershed outlet, and the drainage network are all physically measurable, and determination of these parameters can be accomplished with a reasonable degree of accuracy.The parameter values of USLE such as R, K, C and delivery parameters such as Hc and n were selected using the literature and experience.The GAMES model was run using the developed interface.The simulations for the Madter and Holtby sub-watersheds were performed separately.The model divided the Madter and Haltby watersheds into 352 and 360 HRU's respectively.The GAMES model created a tabular resultant output text file that included potential soil erosion, erosion rate, delivery ratio to the adjacent HRU, delivery ratio to the stream, sediment load, and sediment yield rate for every HRU.The developed GIS interface tool converted the output text file into excel format and loaded it to the ArcGIS.The excel file was then joined with the attribute table of the HRUs polygon file for displaying the results in spatial map format.Figure input