Forma Tarraconis: GIS aplicado a la arqueología urbana

Desde aquel incipiente interés renacentista del siglo XVI hasta la profesionalización de la arqueología en la actualidad, la ingente actividad realizada en Tarragona ha necesitado de formas de registrarla y organizarla, a fin de ser estudiada en el contexto diacrónico general de la ciudad. Esta tarea se ha visto favorecida en el último decenio del siglo pasado por la aparición y adopción de una herramienta informática que no fue diseñada específicamente para la Arqueología. El GIS o SIG, se ha convertido en un instrumento básico e indispensable para cualquier ciudad que tenga un mínimo bagaje patrimonial. Prueba de ello son las publicaciones de Cartas Arqueológicas obtenidas a partir de este software, o sin ir mucho más lejos, los trabajos presentados en este simposio. Presentamos a continuación la experiencia de Tarragona en la creación de un programa, bajo soporte GIS, para la gestión de la actividad arqueológica. Previamente contextualizaremos esta experiencia ya que es justa heredera de iniciativas individuales y de equipos y proyectos previos que trataron, con mucha menor suerte que no éxito, de organizar y poner orden en la maltrecha documentación arqueológica de nuestra ciudad. Describiremos posteriormente las funcionalidades básicas de este programa, su arquitectura y las aplicaciones realizadas en el desarrollo de sendos proyectos de estudio de la ciudad. En este último caso nos centraremos, principalmente, en el proyecto Planimetria Arqueològica de Tàrraco, impulsado por el Museu d’Història del Ayuntamiento de Tarragona y el Institut Català d’Arqueologia Clàssica, y que ha contado con la colaboración de la Generalitat de Catalunya para su edición. Seguidamente veremos como el uso de la información en capas registrada en el GIS combinada con la inclusión de la cartografía histórica de una ciudad permiten el análisis de problemas arqueológicos irresolubles sólo con los datos estrictamente arqueológicos. En este caso se aplicará sobre la validación de la morfología urbana de Tarraco propuesta por Macias (2000a). Por último se hará una reflexión sobre cuales podrían ser las mejoras aplicables sobre el programa, en previsión de futuras revisiones del mismo.

Nutrient cycling and field-based partial nutrient balances in two mountain oases of Oman

Little is known about nutrient fluxes as a criterion to assess the sustainability of traditional irrigation agriculture in eastern Arabia. In this study GIS-based field research on terraced cropland and groves of date palm (Phoenix dactylifera L.) was conducted over 2 years in two mountain oases of northern Oman to determine their role as hypothesized sinks for nitrogen (N), phosphorus (P) and potassium (K). At Balad Seet 55% of the 385 fields received annual inputs of 100–500 kg N ha^-1 and 26% received 500–1400 kg N ha^-1. No N was applied to 19% of the fields which were under fallow. Phosphorus was applied annually at 1–90 kg ha^-1 on 46% of the fields, whereas 27% received 90–210 kg ha^-1. No K was applied to 27% of the fields, 32% received 1–300 kg K ha^-1, and the remaining fields received up to 1400 kg ha^-1. At Maqta N-inputs were 61–277 kg ha^-1 in palm groves and 112–225 kg ha^-1 in wheat (Triticum spp.) fields, respective P inputs were 9–40 and 14–29 kg ha^-1, and K inputs were 98–421 and 113–227 kg ha^-1. For cropland, partial oasis balances (comprising inputs of manure, mineral fertilizers, N2-fixation and irrigation water, and outputs of harvested products) were similar for both oases, with per hectare surpluses of 131 kg N, 37 kg P, and 84 kg K at Balad Seet and of 136 kg N, 16 kg P and 66 kg K at Maqta. This was despite the fact that N2-fixation by alfalfa (Medicago sativa L.), estimated at up to 480 kg ha^-1 yr^-1 with an average total dry matter of 22 t ha^-1, contributed to the cropland N-balance only at the former site. Respective palm grove surpluses, in contrast were with 303 kg N, 38 kg P, and 173 kg K ha^-1 much higher at Balad Seet than with 84 kg N, 14 kg P, and 91 kg K ha^-1 at Maqta. The data show that both oases presently are large sinks for nutrients. Potential gaseous and leaching losses could at least partly be controlled by a decrease in nutrient input intensity and careful incorporation of manure.

Integration of remote sensing and GIS for modelling flash floods in Wadi Hudain catchment, Egypt

This paper presents a new approach to modelling flash floods in dryland catchments by integrating remote sensing and digital elevation model (DEM) data in a geographical information system (GIS). The spectral reflectance of channels affected by recent flash floods exhibit a marked increase, due to the deposition of fine sediments in these channels as the flood recedes. This allows the parts of a catchment that have been affected by a recent flood event to be discriminated from unaffected parts, using a time series of Landsat images. Using images of the Wadi Hudain catchment in southern Egypt, the hillslope areas contributing flow were inferred for different flood events. The SRTM3 DEM was used to derive flow direction, flow length, active channel cross-sectional areas and slope. The Manning Equation was used to estimate the channel flow velocities, and hence the time-area zones of the catchment. A channel reach that was active during a 1985 runoff event, that does not receive any tributary flow, was used to estimate a transmission loss rate of 7·5 mm h−1, given the maximum peak discharge estimate. Runoff patterns resulting from different flood events are quite variable; however the southern part of the catchment appears to have experienced more floods during the period of study (1984–2000), perhaps because the bedrock hillslopes in this area are more effective at runoff production than other parts of the catchment which are underlain by unconsolidated Quaternary sands and gravels. Due to high transmission loss, runoff generated within the upper reaches is rarely delivered to the alluvial fan and Shalateen city situated at the catchment outlet. The synthetic GIS-based time area zones, on their own, cannot be relied on to model the hydrographs reliably; physical parameters, such as rainfall intensity, distribution, and transmission loss, must also be considered.

Meteomap: Generation Of Meteorological Variables For Distributed Physically-Based Hydrological Modeling

Distributed energy and water balance models require time-series surfaces of the meteorological variables involved in hydrological processes. Most of the hydrological GIS-based models apply simple interpolation techniques to extrapolate the point scale values registered at weather stations at a watershed scale. In mountainous areas, where the monitoring network ineffectively covers the complex terrain heterogeneity, simple geostatistical methods for spatial interpolation are not always representative enough, and algorithms that explicitly or implicitly account for the features creating strong local gradients in the meteorological variables must be applied. Originally developed as a meteorological pre-processing tool for a complete hydrological model (WiMMed), MeteoMap has become an independent software. The individual interpolation algorithms used to approximate the spatial distribution of each meteorological variable were carefully selected taking into account both, the specific variable being mapped, and the common lack of input data from Mediterranean mountainous areas. They include corrections with height for both rainfall and temperature (Herrero et al., 2007), and topographic corrections for solar radiation (Aguilar et al., 2010). MeteoMap is a GIS-based freeware upon registration. Input data include weather station records and topographic data and the output consists of tables and maps of the meteorological variables at hourly, daily, predefined rainfall event duration or annual scales. It offers its own pre and post-processing tools, including video outlook, map printing and the possibility of exporting the maps to images or ASCII ArcGIS formats. This study presents the friendly user interface of the software and shows some case studies with applications to hydrological modeling.

A support system that delineates location-choice sets for firms seeking office space

Factors influencing the location decisions of offices include traffic, accessibility, employment conditions, economic prospects and land-use policies. Hence tools for supporting real-estate managers and urban planners in such multidimensional decisions may be useful. Accordingly, the objective of this study is to develop a GIS-based tool to support firms who seek office accommodation within a given regional or national study area. The tool relies on a matching approach, in which a firm's characteristics (demand) on the one hand, and environmental conditions and available office spaces (supply) on the other, are analyzed separately in a first step, after which a match is sought. That is, a suitability score is obtained for every firm and for every available office space by applying some value judgments (satisfaction, utility etc.). The latter are powered by a focus on location aspects and expert knowledge about the location decisions of firms/organizations with respect to office accommodation as acquired from a group of real-estate advisers; it is stored in decision tables, and they constitute the core of the model. Apart from the delineation of choice sets for any firm seeking a location, the tool supports two additional types of queries. Firstly, it supports the more generic problem of optimally allocating firms to a set of vacant locations. Secondly, the tool allows users to find firms which meet the characteristics of any given location. Moreover, as a GIS-based tool, its results can be visualized using GIS features which, in turn, facilitate several types of analyses.

Einfluss der Landnutzungsänderung auf die Variabilität von Grundwasserneubildungs- und Sickerwasserrate - Qualitative Grundwassergefährdungsabschätzung mit einem GIS -

Ziel der Arbeit ist die Entwicklung eines allgemeingültigen, im mittleren Maßstabsbereich durchgeführten, GIS-Bewertungsverfahrens zur Ermittlung der Grundwasserneubildungsrate unter Berücksichtigung einer langzeitlichen Landnutzungsänderung sowie die Darstellung eines Szenarios zur Grundwassergefährdungsabschätzung für ausgewählte rüstungsspezifi-sche organische Verbindungen. Das Untersuchungsgebiet befindet sich in einem während der beiden Weltkriege stark militär- und rüstungstechnisch genutzten Raum. Die Beurtei-lungsszenarien der Grundwassergefährdung erfolgen rein qualitativ.Die Abhängigkeit der Grundwasserneubildungsrate von der Landnutzungsänderung ist deut-lich zu erkennen. Mit zunehmender Flächenversiegelung, einher gehend mit veränderten Sickerwasserraten und des oberflächennahen Abflusses, reduziert sich die Grundwasser-neubildungsrate deutlich. Die Abschätzung der Grundwassergefährdung für die sprengstofftypischen Verbindungen erfolgt im ersten Schritt emissionsbezogen für die ungesättigte Zone unter Berücksichtigung der Standort- und Bodeneigenschaften und der physikalisch-chemischen Stoffeigenschaften. Im zweiten Schritt erfolgt eine immissionsbezogene Betrachtung der potenziellen Belastun-gen von Trinkwasserbrunnen.Trotz starker Generalisierungseffekte bietet dieser stoffspezifische Bewertungsansatz in Verbindung mit physikalischen Methoden für organische Verbindungen einen guten ersten Überblick über eine potenzielle Gefährdungen des Grundwassers. Die Datenbestände der entwickelten GIS-Anwendung sind leicht erweiterbar, so dass eine Abschätzung der potenzi-ellen Gefährdung auch für andere organische Stoffgruppen sowie auch eine Übertragung auf andere Untersuchungsgebiete möglich ist.

Development of an optimization model for biofuel facility size and location and a simulation model for design of a biofuel supply chain

To mitigate greenhouse gas (GHG) emissions and reduce U.S. dependence on imported oil, the United States (U.S.) is pursuing several options to create biofuels from renewable woody biomass (hereafter referred to as “biomass”). Because of the distributed nature of biomass feedstock, the cost and complexity of biomass recovery operations has significant challenges that hinder increased biomass utilization for energy production. To facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization and tapping unused forest residues, it is proposed to develop biofuel supply chain models based on optimization and simulation approaches. The biofuel supply chain is structured around four components: biofuel facility locations and sizes, biomass harvesting/forwarding, transportation, and storage. A Geographic Information System (GIS) based approach is proposed as a first step for selecting potential facility locations for biofuel production from forest biomass based on a set of evaluation criteria, such as accessibility to biomass, railway/road transportation network, water body and workforce. The development of optimization and simulation models is also proposed. The results of the models will be used to determine (1) the number, location, and size of the biofuel facilities, and (2) the amounts of biomass to be transported between the harvesting areas and the biofuel facilities over a 20-year timeframe. The multi-criteria objective is to minimize the weighted sum of the delivered feedstock cost, energy consumption, and GHG emissions simultaneously. Finally, a series of sensitivity analyses will be conducted to identify the sensitivity of the decisions, such as the optimal site selected for the biofuel facility, to changes in influential parameters, such as biomass availability and transportation fuel price. Intellectual Merit The proposed research will facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization in the renewable biofuel industry. The GIS-based facility location analysis considers a series of factors which have not been considered simultaneously in previous research. Location analysis is critical to the financial success of producing biofuel. The modeling of woody biomass supply chains using both optimization and simulation, combing with the GIS-based approach as a precursor, have not been done to date. The optimization and simulation models can help to ensure the economic and environmental viability and sustainability of the entire biofuel supply chain at both the strategic design level and the operational planning level. Broader Impacts The proposed models for biorefineries can be applied to other types of manufacturing or processing operations using biomass. This is because the biomass feedstock supply chain is similar, if not the same, for biorefineries, biomass fired or co-fired power plants, or torrefaction/pelletization operations. Additionally, the research results of this research will continue to be disseminated internationally through publications in journals, such as Biomass and Bioenergy, and Renewable Energy, and presentations at conferences, such as the 2011 Industrial Engineering Research Conference. For example, part of the research work related to biofuel facility identification has been published: Zhang, Johnson and Sutherland [2011] (see Appendix A). There will also be opportunities for the Michigan Tech campus community to learn about the research through the Sustainable Future Institute.

A comparison of landslide inventories along the shoreline of Steamboat Island Peninsula, Thurston County, WA: office based protocol vs. field based protocol

The mountain ranges and coastlines of Washington State have steep slopes, and they are susceptible to landslides triggered by intense rainstorms, rapid snow melts, earthquakes, and rivers and waves removing slope stability. Over a 30-year timespan (1984-2014 and includes State Route (SR) 530), a total of 28 deep-seated landslides caused 300 million dollars of damage and 45 deaths (DGER, 2015). During that same timeframe, ten storm events triggered shallow landslides and debris flows across the state, resulting in nine deaths (DGER, 2015). The loss of 43 people, due to the SR 530 complex reactivating and moving at a rate and distance unexpected to residents, highlighted the need for an inventory of the stateís landslides. With only 13% of the state mapped (Lombardo et al., 2015), the intention of this statewide inventory is to communicate hazards to citizens and decision makers. In order to compile an accurate and consistent landslide inventory, Washington needs to adopt a graphic information system (GIS) based mapping protocol. A mapping protocol provides consistency for measuring and recording information about landslides, including such information as the type of landslide, the material involved, and the size of the movement. The state of Oregon shares similar landslide problems as Washington, and it created a GIS-based mapping protocol designed to inform its residents, while also saving money and reducing costly hours in the field (Burns and Madin, 2009). In order to determine if the Oregon Department of Geology and Mineral Industries (DOGAMI) protocol, developed by Burns and Madin (2009), could serve as the basis for establishing Washingtonís protocol, I used the office-based DOGAMI protocol to map landslides along a 40-50 km (25-30 mile) shoreline in Thurston County, Washington. I then compared my results to the field-based landslide inventory created in 2009 by the Washington Division of Geology and Earth Resources (DGER) along this same shoreline. If the landslide area I mapped reasonably equaled the area of the DGER (2009) inventory, I would consider the DOGAMI protocol useful for Washington, too. Utilizing 1m resolution lidar flown for Thurston County in 2011 and a GIS platform, I mapped 36 landslide deposits and scarp flanks, covering a total area of 879,530 m2 (9,467,160 ft2). I also found 48 recent events within these deposits. With an exception of two slides, all of the movements occurred within the last fifty years. Along this same coastline, the DGER (2009) recorded 159 individual landslides and complexes, for a total area of 3,256,570 m2 (35,053,400 ft2). At a first glance it appears the DGER (2009) effort found a larger total number and total area of landslides. However, in addition to their field inventory, they digitized landslides previously mapped by other researchers, and they did not field confirm these landslides, which cover a total area of 2,093,860 m2 (22,538,150 ft2) (DGER, 2009). With this questionable landslide area removed and the toes and underwater landslides accounted for because I did not have a bathymetry dataset, my results are within 6,580 m2 (70,840 ft2) of the DGERís results. This similarity shows that the DOGAMI protocol provides a consistent and accurate approach to creating a landslide inventory. With a few additional modifications, I recommend that Washington State adopts the DOGAMI protocol. Acquiring additional 1m lidar and adopting a modified DOGAMI protocol poises the DGER to map the remaining 87% of the state, with an ultimate goal of informing citizens and decision makers of the locations and frequencies of landslide hazards on a user-friendly GIS platform.

Development and evaluation of a geographic information system based method to estimate flooding susceptibility in an area of Broward County, Florida

The objective of this study was to develop a GIS-based multi-class index overlay model to determine areas susceptible to inland flooding during extreme precipitation events in Broward County, Florida. Data layers used in the method include Airborne Laser Terrain Mapper (ALTM) elevation data, excess precipitation depth determined through performing a Soil Conservation Service (SCS) Curve Number (CN) analysis, and the slope of the terrain. The method includes a calibration procedure that uses "weights and scores" criteria obtained from Hurricane Irene (1999) records, a reported 100-year precipitation event, Doppler radar data and documented flooding locations. Results are displayed in maps of Eastern Broward County depicting types of flooding scenarios for a 100-year, 24-hour storm based on the soil saturation conditions. As expected the results of the multi-class index overlay analysis showed that an increase for the potential of inland flooding could be expected when a higher antecedent moisture condition is experienced. The proposed method proves to have some potential as a predictive tool for flooding susceptibility based on a relatively simple approach.

UTILIZATION OF GIS AND SPATIAL ANALYSIS TECHNIQUES IN ARTISANAL AND SMALL SCALE MINING TO LOCATE A CENTRALIZED PROCESSING CENTRE

One of the global phenomena with threats to environmental health and safety is artisanal mining. There are ambiguities in the manner in which an ore-processing facility operates which hinders the mining capacity of these miners in Ghana. These problems are reviewed on the basis of current socio-economic, health and safety, environmental, and use of rudimentary technologies which limits fair-trade deals to miners. This research sought to use an established data-driven, geographic information (GIS)-based system employing the spatial analysis approach for locating a centralized processing facility within the Wassa Amenfi-Prestea Mining Area (WAPMA) in the Western region of Ghana. A spatial analysis technique that utilizes ModelBuilder within the ArcGIS geoprocessing environment through suitability modeling will systematically and simultaneously analyze a geographical dataset of selected criteria. The spatial overlay analysis methodology and the multi-criteria decision analysis approach were selected to identify the most preferred locations to site a processing facility. For an optimal site selection, seven major criteria including proximity to settlements, water resources, artisanal mining sites, roads, railways, tectonic zones, and slopes were considered to establish a suitable location for a processing facility. Site characterizations and environmental considerations, incorporating identified constraints such as proximity to large scale mines, forest reserves and state lands to site an appropriate position were selected. The analysis was limited to criteria that were selected and relevant to the area under investigation. Saaty’s analytical hierarchy process was utilized to derive relative importance weights of the criteria and then a weighted linear combination technique was applied to combine the factors for determination of the degree of potential site suitability. The final map output indicates estimated potential sites identified for the establishment of a facility centre. The results obtained provide intuitive areas suitable for consideration

Application of GIS in potential beekeeping assessment: Case study of Montesinho Natural Park (Portugal)

Planners require solutions that address routine work needs and seems essential to improving efficiency and productivity. There are a great number of different factors related to beekeeper activity as well the quality and productivity of different bee products. The spatial analysis is a powerful tool for overlap and relates various levels of information on a map, and consequently a very useful for beekeeping activity planning. This work proposes and applies a methodology to potential beekeeping assessment in Montesinho Natural Park, a region in the northwest of Portugal. The beekeeping potential maps were developed with the following data sources: legal standards, vegetation, land use, topography, water resources, roads, electromagnetic fields, and some honey physico-chemical analysis. The design and implementation of spatial analysis model based on Geographic Information System (GIS) to beekeeping planning activities has already been described by Anjos et al (2014). Spatial analysis techniques allows to define the potential beekeeper map supporting the beekeeper management in this region. Anjos O, Silva G, Roque N, Fernandez P, 2014. GIS based analysis to support the beekeeping planning. Book of abstracts of the International Symposium on Bee Products 3rd edition – Annual meeting of the International Honey Commission (IHC), Faculty of medicine, University of Rijeka, p:61

Small rivers in the southwestern Amazon and their role in CO(2) outgassing

A recent estimate of CO(2) outgassing from Amazonian wetlands suggests that an order of magnitude more CO(2) leaves rivers through gas exchange with the atmosphere than is exported to the ocean as organic plus inorganic carbon. However, the contribution of smaller rivers is still poorly understood, mainly because of limitations in mapping their spatial extent. Considering that the largest extension of the Amazon River network is composed of small rivers, the authors` objective was to elucidate their role in air-water CO(2) exchange by developing a geographic information system ( GIS)- based model to calculate the surface area covered by rivers with channels less than 100 m wide, combined with estimated CO(2) outgassing rates at the Ji-Parana River basin, in the western Amazon. Estimated CO(2) outgassing was the main carbon export pathway for this river basin, totaling 289 Gg C yr(-1), about 2.4 times the amount of carbon exported as dissolved inorganic carbon ( 121 Gg C yr(-1)) and 1.6 times the dissolved organic carbon export ( 185 Gg C yr(-1)). The relationships established here between drainage area and channel width provide a new model for determining small river surface area, allowing regional extrapolations of air - water gas exchange. Applying this model to the entire Amazon River network of channels less than 100 m wide ( third to fifth order), the authors calculate that the surface area of small rivers is 0.3 +/- 0.05 million km(2), and it is potentially evading to the atmosphere 170 +/- 42 Tg C yr(-1) as CO(2). Therefore, these ecosystems play an important role in the regional carbon balance.

Assessment of groundwater recharge and water fluxes of the Guarani Aquifer System, Brazil

The groundwater recharge and water fluxes of the Guarani Aquifer System in the state of Sao Paulo in Brazil were assessed through a numeric model. The study area (6,748 km(2)) comprises Jacar,-Gua double dagger A(0) and Jacar,-Pepira River watersheds, tributaries of the Tiet River in the central region of the state. GIS based tools were used in the storage, processing and analysis of data. Main hydrologic phenomena were selected, leading to a groundwater conceptual model, taking into account the significant outcrops occurring in the study area. Six recharge zones were related to the geologic formation and structures of the semi-confined and phreatic aquifer. The model was calibrated against the baseflows and static water levels of the wells. The results emphasize the strong interaction of groundwater flows between watersheds and the groundwater inflow into the rivers. It has been concluded that lateral groundwater exchanges between basins, the deep discharges to the regional system, and well exploitation were not significant aquifer outflows when compared to the aquifer recharge. The results have shown that the inflows from the river into the aquifer are significant and have the utmost importance since the aquifer is potentially more vulnerable in these places.

The influence of experience and training in a group of novice observers: a jackknife alternative free-response receiver operating characteristic analysis

Purpose - The study evaluates the pre- and post-training lesion localisation ability of a group of novice observers. Parallels are drawn with the performance of inexperienced radiographers taking part in preliminary clinical evaluation (PCE) and ‘red-dot’ systems, operating within radiography practice. Materials and methods - Thirty-four novice observers searched 92 images for simulated lesions. Pre-training and post-training evaluations were completed following the free-response the receiver operating characteristic (FROC) method. Training consisted of observer performance methodology, the characteristics of the simulated lesions and information on lesion frequency. Jackknife alternative FROC (JAFROC) and highest rating inferred ROC analyses were performed to evaluate performance difference on lesion-based and case-based decisions. The significance level of the test was set at 0.05 to control the probability of Type I error. Results - JAFROC analysis (F(3,33) = 26.34, p < 0.0001) and highest-rating inferred ROC analysis (F(3,33) = 10.65, p = 0.0026) revealed a statistically significant difference in lesion detection performance. The JAFROC figure-of-merit was 0.563 (95% CI 0.512,0.614) pre-training and 0.677 (95% CI 0.639,0.715) post-training. Highest rating inferred ROC figure-of-merit was 0.728 (95% CI 0.701,0.755) pre-training and 0.772 (95% CI 0.750,0.793) post-training. Conclusions - This study has demonstrated that novice observer performance can improve significantly. This study design may have relevance in the assessment of inexperienced radiographers taking part in PCE or commenting scheme for trauma.

Assessment of groundwater contamination in an agricultural peri-urban area (NW Portugal): an integrated approach

The excessive use of pesticides and fertilisers in agriculture has generated a decrease in groundwater and surface water quality in many regions of the EU, constituting a hazard for human health and the environment. Besides, on-site sewage disposal is an important source of groundwater contamination in urban and peri-urban areas. The assessment of groundwater vulnerability to contamination is an important tool to fulfil the demands of EU Directives. The purpose of this study is to assess the groundwater vulnerability to contamination related mainly to agricultural activities in a peri-urban area (Vila do Conde, NW Portugal). The hydrogeological framework is characterised mainly by fissured granitic basement and sedimentary cover. Water samples were collected and analysed for temperature, pH, electrical conductivity, chloride, phosphate, nitrate and nitrite. An evaluation of groundwater vulnerability to contamination was applied (GOD-S, Pesticide DRASTIC-Fm, SINTACS and SI) and the potential nitrate contamination risk was assessed, both on a hydrogeological GIS-based mapping. A principal component analysis was performed to characterised patterns of relationship among groundwater contamination, vulnerability, and the hydrogeological setting assessed. Levels of nitrate above legislation limits were detected in 75 % of the samples analysed. Alluvia units showed the highest nitrate concentrations and also the highest vulnerability and risk. Nitrate contamination is a serious problem affecting groundwater, particularly shallow aquifers, especially due to agriculture activities, livestock and cesspools. GIS-based cartography provided an accurate way to improve knowledge on water circulation models and global functioning of local aquifer systems. Finally, this study highlights the adequacy of an integrated approach, combining hydrogeochemical data, vulnerability assessments and multivariate analysis, to understand groundwater processes in peri-urban areas.