The implementation and application of Geographic Information Systems (GIS) and spatial analyses have become standard practice in many archaeological projects. In this study, we demonstrate how GIS can play a crucial role in the study of taphonomy, i.e., understanding the processes that underpinned the creation of archaeological deposits, in this case the distribution of artifacts across an archeological site. The Stelida Naxos Archeological Project (SNAP) is focused on the exploration of a Paleolithic-Mesolithic stone tool quarry site located on the island of Naxos, Greece. An extensive pedestrian survey was conducted during the 2013 and 2014 archeological field seasons. An abundance of lithic material was collected across the surface, with some diagnostic pieces dating to more than 250 Kya. Spatial statistical analysis (Empirical Bayesian Kriging) was conducted on the survey data to generate predictive distribution maps for the site. This study then determined the contextual integrity of the surface artifact distributions through a study of geomorphic processes. A digital surface model (DSM) of the site was produced using Unmanned Aerial Vehicle (UAV) aerial photography and Structure from Motion (SfM) terrain modeling. The DSM employed to develop a Revised Universal Soil Loss Equation (RUSLE) model and hydrological flow models. The model results provide important insights into the site geomorphological processes and allow categorization of the diagnostic surface material locations based on their contextual integrity. The GIS analysis demonstrates that the surface artifact distribution has been significantly altered by post-depositional geomorphic processes, resulting in an overall low contextual integrity of surface artifacts. Conversely, the study identified a few areas with high contextual integrity, loci that represent prime locations for excavation. The results from this study will not only be used to inform and guide further development of the archeological project (as well as representing significant new data in its own right), but also contributes to current debates in survey archaeology, and in mapping and prospection more generally. This project demonstrates the benefit of using spatial analysis as a tool for planning of pedestrian surveys and for predictive mapping of artifact distributions prior to archaeological excavations. / Thesis / Master of Science (MSc)

[Truncated abstract] Landforms rehabilitated after bauxite mining can be vulnerable to soil loss by water erosion processes. On most rehabilitated sites, management controls such as deep ripping, contour mounding and landscaped sub-catchments limit erosion. Despite these measures, severe gully erosion that is anecdotally associated with steep slopes can damage rehabilitated areas and affect downstream drinking water resources. A review of erosion dynamics reveals that gullies develop episodically and in a non-linear manner. They often initiated as a near surface process and are influenced by natural climatic drivers. Despite this, local site characteristics including soil and landform can predispose an area to gully erosion. Moreover, erosion models, becoming more-widely utilized within the mining industry, may provide useful tools with which to measure, analyse, and manage gully erosion. One of these models, SIBERIA was tested to determine its suitability for application a tool to help manage erosion risk. We first surveyed 26 eroding and erosion-prone rehabilitated hillslopes to determine the common form and setting for gully erosion on these rehabilitated bauxite mines. A conceptual model was developed to include and explore the interplay between the common causes of the gullies surveyed. The conceptual model accounts for slope steepness but suggests that additionally, certain triggers and threshold effects operating under different site conditions are as influential (or even more influential) than slope steepness as determinants of gully erosion occurrence and severity. ... Soil properties and soil erodibility had some subtle influence on landform stability and erosion risk. The most-erodible media occurred where either: mine floor material was mixed with topsoil/ overburden; and/or the topsoil/overburden layer was thin or its coverage is patchy resulting in slaking subsoil, hardsetting soil and surface crusts. When erodible surface media were combined with steeper (>8[degrees]) or longer (>50 m) slopes or with any major erosion trigger, rill and gully development was greatly intensified. The SIBERIA simulation model was calibrated and its simulated outputs were compared to known locations of gully erosion on a steep, rehabilitated pit from the Willowdale mine. At a resolution of one metre, SIBERIA was able to simulate the approximate dimensions of gullies. However, SIBERIA could not simulate the exact location of individual gully headcuts. Additionally, SIBERA was able to simulate the effect of different microtopographic surface treatments but this was only achieved by increasing the grid resolution to 25 cm and reducing the size of the area simulated due to model constraints. Locations of gully headcuts were overlain onto a grid-based, Digital Elevation Model (DEM) using a Geographic Information System (GIS). The spatial distribution of gully headcut locations was compared to DEM derivatives such as slope and flow accumulation. Positive, and predictive relationships allow between the steepness of the slope of the pre-mining landform and the cell count of the area contributing to flow (catchment), as determined by GIS, may allow a mine scale indication of erosion risk using simple GIS desktop analysis.