Archaeological Sites in Small Towns—A Sustainability Assessment of Northumberland County

: Background: Land-use change, resulting from urban expansion, has led to the unprecedented importance of rethinking our interactions with the environment in general. The consequences of anthropogenic activity on natural land, as well as sustainable land use types, is a signiﬁcant concern to the detriment of small towns. In this article, I assess the changes in the urban footprint morphology of Northumberland County while evaluating the consequences of urban fragmentation on archaeological sites. Northumberland has a total population of 85,598 divided into seven municipalities and the Alderville First Nation reserve. Northumberland County is a commuting distance of one hour from Toronto, one of North America’s most important economic hubs. It is of the highest importance to measure and understand the current characteristics of land-use change and its impacts on the sustainability of small towns and local communities. Methods: A total of archaeological sites are used to test the impacts of urbanization and land-use change throughout the perimeter of Northumberland county. Spatial analysis is conducted employing spatial autocorrelation techniques (global and local) to assess the existence of hotspots of archaeological sites and the potential endangerment of these sites due to urban pressure. A hexagonal topology is used to derive a regional understanding of impacted regions within the county, allowing the richness of archaeological heritage sites as well as the possible concerns of these spatial clusters to be examined. Results: It is shown that a signiﬁcant number of archaeological sites are present in Northumberland county, with a clear indication of a hotspot around Rice Lake. Future urbanization as well as land-use fragmentation may have an impact on the sustainability of these archaeological sites. It is shown that hexagonal geometries may establish relevant spatial surfaces to compare urban and archaeological sites as criteria to understand the endangerment of archaeological sites. The availability of the Global Urban Footprint (GUF) creates a particularly rigorous framework for the detection of impact on archaeological sites. Conclusion: The ﬁndings suggest that archaeological sites in small towns can beneﬁt from high-resolution urban land-use footprints to detect impacts on archaeological sites and the relations of land-use change and intrasite management. This is of utmost importance in times of signiﬁcant anthropogenic change brought by increasing urban pressure on small and medium-sized towns.


GIS, Archaeology and Small Towns
Archaeological material evidence has led to the collection of spatial data related to anthropogenic and archaeological artefacts with a significant level of geographical detail [1]. Archaeological material evidence, as is the case for archaeological sites with specific spatial characteristics, may offer specific spatial characteristics [2]. The combination of anthropogenic sites at a spatial level supports not only the understanding of our territorial identity but incites the possibility of envisioning history as a product for regional development [3]. Archaeological sites foster a vital role in the historical knowledge of local communities and offer differentiation between small towns. Canadian archaeology brings a well-documented richness of anthropological interactions since prehistory [4]. The available remains from pre and post-settlement sites paint a complex picture of archaeological wealth that should be a part of our shared heritage and identity [5]. The aboriginal identity, in particular its richness and historical importance for the Canadian identity, is a driving force of the very core of Canadian values, sustainability and shared living heritage [6].
The region of southern Ontario has been a cradle for different societies, and inter-site analysis allows for an understanding of the topology of human behavior since the Paleo-Indian occupation [7]. The combination of archaeological evidence, together with site location, allows understanding of regional patterns of spatial topology merged with environmental characteristics [8][9][10]. By employing spatial analysis, archaeological evidence may bring the opportunity to assess past territorial behavior, as well as document our cultural complexity [11].
Through a systematic analysis of archaeological evidence, territorial uniqueness may further enable understanding of how to gear the service sector to provide the advantage of using cultural heritage in relevant economies such as tourism, and assure that the archaeological sites are protected and remain an intangible value for generations to come [12]. Within the context of southern Ontario, the richness of local and regional history may be explored efficiently by understanding the typology and distribution of archaeological sites. This is furthered by the potential of sensitizing society to the existing heritage, and by offering solutions to meso-economic issues, sustainable and more integrated directions of regional development [13]. The effects of tourism, however, must be carefully analyzed as to not overburden the carrying capacity of landscape, environment and ecosystems.
Excessive tourism may lead unavoidably to land-use change, pollution and the vulnerability of ecosystems [14][15][16]. The understanding of cultural legacy at the regional level starts with the application of spatial technologies in the context of detecting patterns of the urban footprint concerning the territorial dimension of land use planning [17,18]. Although at times left unconsidered in Canadian planning, archaeology does have an important role and has long relied on survey data, remote sensing [19], spatial analysis [20] and, recently, agent based models [21] and other techniques such as machine learning and artificial intelligence [22] to facilitate the historical understanding of past landscapes. Preservation of archaeological sites allow sustainable landscapes to be supported, particularly within small towns that rely on (i) product differentiation, (ii) seasonal tourism, (iii) endogenous economies and (iv) conservation of their small-town landscape. Archaeological sites can have an utmost important role and support local communities. By protecting archaeological sites, small towns have the chance to renew their position of unique hubs for historical and cultural landscapes.

The Importance of Aboriginal Archaeological Sites in Ontario
The complexity and controversy of the first settlement of the Americas is predominantly a result of historical, cultural, archaeological and genetic evidence [23,24], where present technology is used to gain a better understanding of this century-old debate. The controversies of the first settlement revolve around three distinct paradigms relating to the early humans settled in the New World [25]. These questions are best answered by the assessment of distinct hypotheses that build up the framework of archaeological and cultural evidence among which the following are the most relevant: (a) the Clovis hypothesis [26], (b) the coastal entry hypothesis [27] and (c) the Solutrean Hypothesis [28]. The Clovis hypothesis sets out as the most commonly accepted hypothesis that the first settlement of the Americas started during the late Ice Age when hunter-gatherers were entering North America through the Bering Land Bridge. The geomorphological characteristics of the Bering Land Bridge join the regions of Siberia and Alaska. With the variation of sea level, during the Wisconsinan glaciation (between 75,000 to 11,000 BCE), this became a possible migration path from Asia. With the colonization of Beringia potentially sometime between 20,000 and 10,000 BCE, a path of possible entry between the glaciers allowed entry to North America [29]. This allowed the Clovis people, as big-game hunters, to follow along with the corridor into the Great Plains, although with limitations concerning the geographical extent due to ice sheets further north. The apparent extinction of large animals in the Americas [30] corroborates with this settlement of humans in North America sometime around 12,000 [31].
With growing evidence, the coastal entry hypothesis suggests arrival by the first humans in America by boats [32]. Pre-Clovis cultures in South and North America seem remnant of an earlier date than the recent ice age [33]. Material evidence brought from circa 35,000 BCE identifies the possibility of navigating inland by boats through the coasts and rivers of Central and South America. The geographical choice of settlement in the warmer regions leads up to the potential of pre-Clovis people expanding northwards during the ice recession around 11,000 years ago.
For more than a century have similarities between Siberians and North American Indians been pointed out. It is, however, an analysis performed by [34] in regards to dental characteristics that disputes the concept of sinodonty, a morphological similarity between northern Asia and the Americas, with potential links to over 40,000 BCE. Very recent findings strongly connect, through archaeology again, the pre-Clovis timeline through findings of hammerstones and stone anvils dating to 130,000 years ago [35]. Nevertheless, a response publication by Haynes [36], debates the stratigraphic characteristics of these findings, shedding doubt on the archaeological evidence itself found in the Cerutti Mastodon site. This clearly shows the juxtaposing positions on pre-Clovis and Clovis humans and is a very recent example of the complex decade-old debate of the first settlement in the Americas. Ontario may have a leading role in regard to the spatial morphology through archaeological heritage to better understand prehistorical site distribution. It is through aboriginal history and archaeology that we may have a better idea of our Ontario's prehistorical origins. The exploration of these subjects brings a picture of intangible heritage and an opportunity for cherishing our past in the remote and rural landscape of small towns.

Study Area
Northumberland County is an upper-tier level of municipal government situated on the north shore of Lake Ontario, east of Toronto in Central Ontario. It comprises the town of Cobourg; three municipalities, namely Port Hope, Trent Hills and Brighton; and the three townships of Hamilton, Alnwick/Haldimand and Cramahe. It includes the following towns: Cobourg, Port Hope, Brighton, Campbellford, Colborne, Hastings, Baltimore, Bewdley, Camborne, Centreton, Grafton and Warkworth. From 1802 to 1849, jointly with Durham County, it formed the Newcastle District as well as the United Counties of Northumberland and Durham between 1850 to 1973. It was on January 1st, 1974 that Northumberland County became a county, detaching itself from the Durham County. With a history of the United Empire Loyalists escaping the American Colonies in the 18th century, the crown provided plots of land to settlers. It is important to note that following the War of 1812, some of the towns within Northumberland, particularly Port Hope and Cobourg, became an important commercial center as well as a landing point for European immigrants arriving in steamers.
These unique characteristics have lent to Northumberland county and particularly to the towns of Port Hope and Cobourg, a significant historical character, which remains well documented and maintained. Overall, towns such as Cobourg, as well as the Northumberland County itself, are quite proud of their history and legacy, maintaining an active cultural agenda.

GlobeLand30
GlobeLand30 is a unique project that offers global land cover throughout the world. Created by the National Geomatics Centre of China (NGCC), it provides comprehensive land-use coverage at a resolution of 30 m. With two available time sets (2000 and 2010), image rectification consisted of a total of ten thousand data images gathered from the Landsat Thematic Mapper (TM) and the Enhanced Thematic Mapper Plus (ETM+). Utilizing a pixel-based method, this data holds an accuracy of over 80 percent at a global level.
The global nature of this data allows for a significant amount of distinct analytics that assesses, monitors and quantifies land cover change. For the province of Ontario, GlobeLand30 is a consistent mapping product that permits a thorough understanding of land-use in more remote places. Given the often-unavailable land use coverage outside of highly-dense urban cores, this land use data brings a set of novel explorative analysis in regard to land cover change, possible impacts of urbanization, loss of forest area, impacts on agricultural land and the overall environmental and regional sustainability of the dynamics of artificial land-use change.

Northumberland's Urban Footprint
The urban footprint ( Figure 1) brings a concise understanding of artificial land and allows for assessment at the global and regional levels [37]. This definition of urban land use is particularly important for the assessment of small towns given their rapid change as well as their sprawl potentialized by economic growth into agricultural land. This is aggravated by the fact that many small towns increasingly become satellites to more massive urban cores, a problem increasingly witnessed in Asia [38,39], potentialized by population growth, which may be expected in the coming decades in Ontario through immigration patterns. The integration of different anthropogenic land activities into the urban footprint allows for a spatial and temporal assessment. This brings a systemic approach to the consequences of urbanization and land-use and impacts on small towns. The urban footprint for Northumberland County was analyzed through the integration of the Global Urban Footprint (GUF) dataset from Germany's Earth Observation Centre. Input data are single look slant range complex (SSC) TerraSAR-X and TanDEM-X amplitude images acquired in Stripmap mode at 3 m ground resolution [40,41]. Enhanced Thematic Mapper Plus (ETM+). Utilizing a pixel-based method, this data holds an accuracy of over 80 percent at a global level. The global nature of this data allows for a significant amount of distinct analytics that assesses, monitors and quantifies land cover change. For the province of Ontario, GlobeLand30 is a consistent mapping product that permits a thorough understanding of land-use in more remote places. Given the often-unavailable land use coverage outside of highly-dense urban cores, this land use data brings a set of novel explorative analysis in regard to land cover change, possible impacts of urbanization, loss of forest area, impacts on agricultural land and the overall environmental and regional sustainability of the dynamics of artificial land-use change.

Northumberland's Urban Footprint
The urban footprint ( Figure 1) brings a concise understanding of artificial land and allows for assessment at the global and regional levels [37]. This definition of urban land use is particularly important for the assessment of small towns given their rapid change as well as their sprawl potentialized by economic growth into agricultural land. This is aggravated by the fact that many small towns increasingly become satellites to more massive urban cores, a problem increasingly witnessed in Asia [38,39], potentialized by population growth, which may be expected in the coming decades in Ontario through immigration patterns. The integration of different anthropogenic land activities into the urban footprint allows for a spatial and temporal assessment. This brings a systemic approach to the consequences of urbanization and land-use and impacts on small towns. The urban footprint for Northumberland County was analyzed through the integration of the Global Urban Footprint (GUF) dataset from Germany's Earth Observation Centre. Input data are single look slant range complex (SSC) TerraSAR-X and TanDEM-X amplitude images acquired in Stripmap mode at 3 m ground resolution [40,41].

Archaeological Sites
The Archaeological sites were made available by Ontario's Ministry of Heritage, Sport, Tourism and Culture Studies. The Archaeological site collection consists of a site record for every new site

Archaeological Sites
The Archaeological sites were made available by Ontario's Ministry of Heritage, Sport, Tourism and Culture Studies. The Archaeological site collection consists of a site record for every new site discovered in the province. For the study area, a total of 294 sites were considered encompassing the following affinities (Table 1). A total of 203 archaeological sites comprise Aboriginal affinities, while 90 sites have material evidence of Euro-Canadian affinity, and 22 sites have evidence of Iroquoian affinity. It is expected that due to the stratigraphic nature of archaeological sites, several time periods would overlap at the inter-sites level. Evidence of material remains from different time periods and groups were found in 24 sites. One should note that many of these sites have not been subdued to excavations, and material evidence simply reported back from the database. The figure below ( Figure 2) shows the timeline of the discovery of archaeological sites in the region.
Sustainability 2020, 11, x FOR PEER REVIEW 5 of 13 discovered in the province. For the study area, a total of 294 sites were considered encompassing the following affinities (Table 1). A total of 203 archaeological sites comprise Aboriginal affinities, while 90 sites have material evidence of Euro-Canadian affinity, and 22 sites have evidence of Iroquoian affinity. It is expected that due to the stratigraphic nature of archaeological sites, several time periods would overlap at the inter-sites level. Evidence of material remains from different time periods and groups were found in 24 sites. One should note that many of these sites have not been subdued to excavations, and material evidence simply reported back from the database. The figure below (Figure 2) shows the timeline of the discovery of archaeological sites in the region.

Global Spatial Autocorrelation
Several steps were performed to assess the spatial characteristics of the impact of land-use change on archaeological sites in Northumberland county. Further to importing the point-data into a GIS, a first step consisted of performing a Moran's I analysis to test whether there was any evidence of spatial autocorrelation over the geographical site locations [42]. This statistic was conducted to test the null hypothesis (Ho) relating to the absence of spatial clustering of archaeological sites in Northumberland County (α = 0.05) (Equation (1) where corresponds to a binary weight matrix defined with the weight of one, given contiguity of adjacency for any value that holds = and any value without adjacency as w ij = 0. The product of the distance is defined as x i for any location in the distance to relation of its mean. This holds  30 35 1962 1966 1967 1968 1971 1973 1974 1975 1976 1977 1978 1979 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1995 1996 1998

Global Spatial Autocorrelation
Several steps were performed to assess the spatial characteristics of the impact of land-use change on archaeological sites in Northumberland county. Further to importing the point-data into a GIS, a first step consisted of performing a Moran's I analysis to test whether there was any evidence of spatial autocorrelation over the geographical site locations [42]. This statistic was conducted to test the null hypothesis (Ho) relating to the absence of spatial clustering of archaeological sites in Northumberland County (α = 0.05) (Equation (1)): where w ij corresponds to a binary weight matrix defined with the weight of one, given contiguity of adjacency for any value that holds w ij = 1 and any value without adjacency as w ij = 0. The product of the distance is defined as x i for any location i in the distance to relation of its mean. This holds as a statistic for assessing the entire spatial distribution of adjacency formed for the county. The null hypothesis was rejected, suggesting global spatial autocorrelation for archaeological sites in the county.

Local Spatial Autocorrelation
The Getis-Ord statistic was calculated by first determining a hexagon grid for archaeological site locations. This hexagon grid establishes the count of total archaeological sites falling within a 1-kilometer perimeter. While several approaches allow for spatial density estimation, I considered that the unavailability of small census tracts within the county should allow for grid analytics to create smaller spatial boundaries. This enabled an assessment of significantly greater spatial detail that otherwise would not have been possible. The usage of hexagon geometry provided the ideal setting for further analysis of local spatial autocorrelation. This approach led to the creation of a density surface of archaeological sites that defined an aggregation of hotspots and coldspots. The Getis-Ord algorithm is performed as follows: where w ij is the spatial weight matrix following a 1 km distance (d), and w ij (d) is assumed as 1. The maps show densities of injury patient residences as hot spots and cold spots, with red representing the highest concentrations of injury and blue the lowest.

The Spatial Geomorphology of Small Towns
Archaeological sites showed a positive spatial autocorrelation. The rejection of the null hypothesis confirmed a positive association between the location of archaeological sites and geographical space with a Global Moran's I of 0.307 (Figure 3).
Sustainability 2020, 11, x FOR PEER REVIEW 6 of 13 as a statistic for assessing the entire spatial distribution of adjacency formed for the county. The null hypothesis was rejected, suggesting global spatial autocorrelation for archaeological sites in the county.

Local Spatial Autocorrelation
The Getis-Ord statistic was calculated by first determining a hexagon grid for archaeological site locations. This hexagon grid establishes the count of total archaeological sites falling within a 1kilometer perimeter. While several approaches allow for spatial density estimation, I considered that the unavailability of small census tracts within the county should allow for grid analytics to create smaller spatial boundaries. This enabled an assessment of significantly greater spatial detail that otherwise would not have been possible. The usage of hexagon geometry provided the ideal setting for further analysis of local spatial autocorrelation. This approach led to the creation of a density surface of archaeological sites that defined an aggregation of hotspots and coldspots. The Getis-Ord algorithm is performed as follows: where is the spatial weight matrix following a 1 km distance (d), and ( ) is assumed as 1. The maps show densities of injury patient residences as hot spots and cold spots, with red representing the highest concentrations of injury and blue the lowest.

The Spatial Geomorphology of Small Towns
Archaeological sites showed a positive spatial autocorrelation. The rejection of the null hypothesis confirmed a positive association between the location of archaeological sites and geographical space with a Global Moran's I of 0.307 (Figure 3). Historical path-tendency of resettlement suggests that human settlements are located in similar areas where past settlements existed, albeit subject to human bias. As further explored in the conclusion section, archaeological sites are also often discovered due to anthropic intervention due to economic activity. This is a fundamental aspect of exploring the sustainability of archaeological preservation and raises the importance of using GIS efficiently for heritage management, as well as the monitoring that should be given during urbanization. The stratigraphic analysis confirms that Historical path-tendency of resettlement suggests that human settlements are located in similar areas where past settlements existed, albeit subject to human bias. As further explored in the conclusion section, archaeological sites are also often discovered due to anthropic intervention due to economic activity. This is a fundamental aspect of exploring the sustainability of archaeological preservation and raises the importance of using GIS efficiently for heritage management, as well as the monitoring that should be given during urbanization. The stratigraphic analysis confirms that several sites have inter-archaeological material evidence from different epochs. This furthers the need to articulate careful planning of cities to protect archaeological sites in a multi-dimensional framework of material evidence. The confirmed spatial hotspots brought a robust understanding of the exact location of regions within small towns that need to focus on continuous monitoring and assessment (Figure 4a,b).
Sustainability 2020, 11, x FOR PEER REVIEW 7 of 13 several sites have inter-archaeological material evidence from different epochs. This furthers the need to articulate careful planning of cities to protect archaeological sites in a multi-dimensional framework of material evidence. The confirmed spatial hotspots brought a robust understanding of the exact location of regions within small towns that need to focus on continuous monitoring and assessment (Figure 4a,b).  Of the identified hotspots and coldspots, the hotspots are particularly noteworthy of further investigation. Of the 17 identified hotspots, five of the hexagons had a p = 0.001, that is, a robust spatial clustering of archaeological sites. The usage of a hexagon mesh clearly allowed a topological division of a terrain where otherwise spatial information would be nonexistent, due to the coarseness of census tract information and the size of forward sortation areas (FSA). It is of extreme interest to understand the local and regional interactions of these hexagons within the highly clustered profile of archaeological sites. This was conducted by comparing the data resulting from the GlobeLand30 initiatives as well as the Global Urban Footprint.

Comparing Land Use for Small Towns
While a substantial amount of research has been conducted on global and regional land use/land cover change, studies of land-use change within small towns are mostly nonexistent. This is largely due to the unavailability of data sources, such as demographic data and land use/land cover, and the nonexistence of multiple timestamps for spatial data. The generation of a hexagon mesh for local analysis within smaller administrative extents may provide a useful solution when combined with global land-use change data sets. The hotspot hexagons were compared with the two global land-use cover sets ( Figure 5). Of the identified hotspots and coldspots, the hotspots are particularly noteworthy of further investigation. Of the 17 identified hotspots, five of the hexagons had a p = 0.001, that is, a robust spatial clustering of archaeological sites. The usage of a hexagon mesh clearly allowed a topological division of a terrain where otherwise spatial information would be nonexistent, due to the coarseness of census tract information and the size of forward sortation areas (FSA). It is of extreme interest to understand the local and regional interactions of these hexagons within the highly clustered profile of archaeological sites. This was conducted by comparing the data resulting from the GlobeLand30 initiatives as well as the Global Urban Footprint.

Comparing Land Use for Small Towns
While a substantial amount of research has been conducted on global and regional land use/land cover change, studies of land-use change within small towns are mostly nonexistent. This is largely due to the unavailability of data sources, such as demographic data and land use/land cover, and the nonexistence of multiple timestamps for spatial data. The generation of a hexagon mesh for local analysis within smaller administrative extents may provide a useful solution when combined with global land-use change data sets. The hotspot hexagons were compared with the two global land-use cover sets ( Figure 5).  The Global Urban Footprint (GUF) provided an excellent understanding of the existing urban footprint within the archaeological cluster. The inspection was carried out by converting the hotspot hexagons as a kml file, allowing a detailed examination of the geographical surrounding of the area to be performed on Google Earth Pro ( Figure 6). This inspection allowed us to look at the economic profile of the surrounding small towns beyond the geographical evidence found through traditional spatial analytical tools. The hotspot for archaeological sites is predominantly located to the southwest of Rice Lake, with the area encompassing smaller communities such as Bewdley, Taits Beach and Gores Landing. The proximity to the lake and the availability of cottages make these communities particularly prone to seasonal tourism, especially when considering the future expansion of cities such as Toronto [43] and Oshawa, towards transportation efficiency in the direction of Ottawa, passing through Port Hope and Cobourg. Overall, while the GlobeLand30 for 2000 and 2010 did not infer much information, the urban footprint showed a significant level that allows one to parse the spatial characteristics of land use.

The Role of GIS and Archaeology in Small Towns
Despite the rich historical legacy, past land use and the relation to archaeological sites have not been explored robustly for small towns in Ontario. However, the availability of archaeological catalogs and geo-information data concerning the environment, supported by spatial analysis and GIS, enables a better understanding of past land-use patterns [44]. As the causal relations between georeferenced sites are established, environmental and geographical characteristics can help to recover spatial dimensions of archaeological evidence together with changing urban environments [45,46]. Such socio-economic niches seem to be common throughout the entire province [47]. Understanding aboriginal archaeological evidence allows for the generation of historical landscapes that bring added value to the sustainability of local and regional development [48]. The collection of this spatial information has enabled a database to be compiled with hundreds of occurrences dated from the 1970s to the present. This becomes a rich dataset for further analysis and should be integrated into multidisciplinary frameworks relating economic growth and sustainability of small towns and communities. This inspection allowed us to look at the economic profile of the surrounding small towns beyond the geographical evidence found through traditional spatial analytical tools. The hotspot for archaeological sites is predominantly located to the southwest of Rice Lake, with the area encompassing smaller communities such as Bewdley, Taits Beach and Gores Landing. The proximity to the lake and the availability of cottages make these communities particularly prone to seasonal tourism, especially when considering the future expansion of cities such as Toronto [43] and Oshawa, towards transportation efficiency in the direction of Ottawa, passing through Port Hope and Cobourg. Overall, while the GlobeLand30 for 2000 and 2010 did not infer much information, the urban footprint showed a significant level that allows one to parse the spatial characteristics of land use.

The Role of GIS and Archaeology in Small Towns
Despite the rich historical legacy, past land use and the relation to archaeological sites have not been explored robustly for small towns in Ontario. However, the availability of archaeological catalogs and geo-information data concerning the environment, supported by spatial analysis and GIS, enables a better understanding of past land-use patterns [44]. As the causal relations between georeferenced sites are established, environmental and geographical characteristics can help to recover spatial dimensions of archaeological evidence together with changing urban environments [45,46]. Such socio-economic niches seem to be common throughout the entire province [47]. Understanding aboriginal archaeological evidence allows for the generation of historical landscapes that bring added value to the sustainability of local and regional development [48]. The collection of this spatial information has enabled a database to be compiled with hundreds of occurrences dated from the 1970s to the present. This becomes a rich dataset for further analysis and should be integrated into multidisciplinary frameworks relating economic growth and sustainability of small towns and communities.
One should consider the nature, however, of the archaeological record itself. While it is certain that protecting archaeological sites in small towns is of utmost importance to address the identity and structure of cultural landscapes and the preservation of small and medium-town sustainability, it is paramount to address the spatial distribution of archaeological sites. There seems to be a relation between the evidence of the archaeological record and the expectation of finding archaeological sites due to anthropogenic activity, such as urbanization, agriculture and surveying. While future expansions of small and medium-sized towns are likely certain, the sustainability of archaeological sites is not. Intense development and urbanization will likely expose new sites in the future. Exposing these sites must be handled both with care as well as with the awareness that the archaeological record might well be enhanced, and the distribution of past material evidence better understood, exploring a more concise and consistent picture of archaeological evidence.
Since the beginning of processual archaeology, material evidence obtained from archaeological excavations has had a vital role in the interpretation of material culture [49]. Because of their inherent spatial character, past activities based on spatial archaeological data become easier to understand. Technological advances, statistical approaches, survey methods and available information have become more critical supporting methodologies for "quantifying archaeology" [50]. In the context of urban and regional planning [51], environmental determinism in archaeology and urban footprint analysis may lead to the following policy implications: (a) the possibility to observe past land use has shed light on road networks, which leads to monuments of historical interest; (b) to develop specific products available in towns, that amount to historical and archaeological preservation and knowledge dissemination; (c) the articulation between tourism and archaeology or cultural heritage leading to spatial environmental determinism already exists because of available infrastructures that complement the provision of the already available tourist industry.

Establishing a Small Town Framework for Archaeological Preservation
Over the years, archaeology has dramatically benefitted from spatial analysis and surveillance by remote sensing techniques, as well as database management and GIS in general. The ubiquity of areas in which archaeology benefits has such a broad spectrum that the correct manipulation of data and research of collected material is often non-linear. Such complexity involves many different actors with different needs and demands within the archaeological subject. Information should allow the creation of innovative scientific processes involving different actors in the archaeological frame. It is of utmost importance to position these debates in line with the development and sustainability of small communities and the growth of small towns as a driver for sustainable choices and heritage preservation [52]. Small towns have to consider the impact in the medium and long run of their policies and governance and become more resilient [53]. Sustainable growth should be a prerequisite of understanding the impact not only on archaeological sites but the consequential impact of land use/land cover change in general [54,55]. In southern Ontario, where massive urban sprawl has led to significant and often irreversible land loss, this must be carefully considered as small towns will significantly change in the coming decades.