Decision Support System Development to Groundwater Management and Aquifer Vulnerability Assessment: Hydrogeological Information System of Monastir (HISM) ()
1. Introduction
The Monastir region, characterized by water scarcity, is experiencing a significant demographic and industrial growth. The problem of groundwater degradation is exacerbated by the climatic conditions and water requirements for different anthropogenic activities and by pollution sources [1]. This degradation is generally expressed by the groundwater level depletion and the decline in the quality of these waters. Reliability and validity of groundwater analysis strongly depend on the availability of large volumes of high-quality data. Putting all data into a coherent and logical structure supported by a computing environment helps to ensure validity and availability and provides a powerful tool for hydrogeological studies [2].
In this situation, a DSS for the management and the planning of these resources seems essential. As a matter of fact, the Geographic Information Systems (GIS) offer modern tools to create maps, integrate information, visualize scenarios, solve complex problems and implement effective solutions based on the data collected from various sources and different formats [3] [4] [5]. State-of-the-art studies had disclosed the higher impact of applying GIS systems on groundwater analysis and hydrogeological decision support. Many studies [6] [7] had contributed to the development of approaches and methods integrated within geographical information systems in order to help users to understand current states of groundwater resources in their regions. In such projects [8] [9] [10] [11] , GIS systems had been deployed to improve the management of groundwater resources supply and urban allocation. On the other side, geographic GIS systems are just deployed as geodatabases [12]. In this paper, we propose to develop a GIS system that integrates at the same time already mentioned features (database and spatial analysis approach).
From a technical perspective, the design of most information system is based on the MERISE approach [13] [14] [15]. In our study, the design and implementation of the DSS are based on Object-Oriented approach. The object approach combines the data and treatments ensuring, thus a certain level of coherence. Indeed, the UML approach allows model complex systems and express in a way integrated the structure and dynamics with the same concept more close to the reality of things.
In this context, the objectives were: 1) to develop a hydrogeological database linked to an Information System enabling the storing, editing, and displaying of the geographic referenced data, 2) to process and analyze the spatially distributed data, 3) to properly support aquifer-vulnerability assessments, and 4) to introduce a conceptual model of groundwater functioning. A GIS-managed hydrogeological database has been developed in order to provide data used in vulnerability assessment techniques and analyses.
The paper is structured as follows. This introductory section is followed by a general overview of GIS and object approach. The third section offers a detailed description of the methodology employed from design to software implementation and development. The application of the DSS to the superficial coastal aquifer of Moknine is discussed in the fourth section. The main conclusions are underlined in the last section.
2. A General Overview of GIS and Object Approach
2.1. GIS Systems
GIS is now well recognized as a science and can be considered as a relatively new and emerging area [16]. It is defined as a system for input, storage, manipulation, and output of geographically referenced data [17]. The GIS technology is used for managing spatial and non-spatial data, linking models, and providing interfaces between the models and their users [18].
Geographic data can be represented by GIS systems as either objects or fields. The objects represent entities are characterized by geometry, topology, and non-spatial attribute values. The fields represent the real data attributes and components [19]. In fact, this distinction between objects and fields in GIS is often associated with vector and raster data models. On one hand, the vector models represent spatial phenomena through differences in the distribution of properties of points, lines, and areas. In such system, each layer is an adapted combination of one or more classes of geometrical features. On the other hand, a raster model consists of a rectangular array of cells with values being assigned to each cell. In such model, each cell is usually restricted to a single value. Consequently, the spatial representation of distributed parameters and variables requires an architecture based on multiple layers and relationship models where different data needs to be well manipulated, analyzed and combined [2] [20].
2.2. Object-Oriented Development
Object-oriented development is a methodology that syntheses approach of both conceptual modeling and object-oriented programming [21]. Object-Oriented Analysis (OOA) addresses conceptual expressions of real-world processes or systems within the object-oriented paradigm [22]. Object-Oriented Design (OOD) can be used to decompose a problem into manageable parts and depict logical aspects of a system using consistent notations [23].
These methodologies are based on the concepts of objects and classes with relations of abstraction, encapsulation and inheritance which are powerful features of object orientation. Encapsulation is a mechanism for improving the flexibility and comprehensibility of a system allowing the development time saving, the integrity of the analyzed data and the correct writing of services [21] [24]. Polymorphism allows different instantiations of objects where each one can be presented in an appropriate way according to the others which allow the system extension [25]. Inheritance is a mechanism by which a class (subclass) refines the behavior and properties of some other class [21].
All these concepts can be performed using Unified Modeling Language (UML). In fact, researchers have adopted the UML as graphical notation to support object-oriented modeling and application development [25]. UML can be used to develop conceptual models of processes and facilitate implementation of databases and systems depending on the computational environment [25]. UML is a modeling language using text and graphical notations that assures specification, analysis, design and implementation of object-oriented system development process. It can be used with all processes throughout the development life cycle and across different implementation technologies [21]. Each object is identified using a unique value known as an object identifier (OID). This unique identifier enables an object to distinguish itself from other objects when such behavior is specified [23]. Similar objects are grouped into classes. A class can be defined as a collection of physical or logical things uniformly described using the same attributes, operations, relationships, and semantics [26]. Classes are the primary organizational unit of object-oriented systems, providing a means for describing the structure and behavior of individual system objects [25]. Classes define properties and behavior of sets of objects. An attribute represents some property of the class [21]. Behavior of objects is realized via operations. An operation is a function or transformation that may be applied on objects in a class. The execution of a method may change the state of an object [21]. A method is an abstraction of something that an object can do, or have done to it [27]. Relationships provide a pathway for communication between objects [28]. There are connections between things and appear on UML class diagrams as lines connecting classes [22]. Generalization and associations are the most frequently occurring type of relationship. The application of generalization enables description of systems and processes using a hierarchy of classes [29]. Associations represent static relationships between classes.
3. Methodological Approach
A practical methodology in this study is presented as flowchart in Figure 1 [30]. The first phase involves searching for, retrieving and collecting data related to the study area. This is the longest and most expensive phase [31]. In the second
Figure 1. Flowchart of methodology for data acquisition and processing [30].
phase, the collected data should be modeling and introduced in digital form. The DSS is developed in Visual Basic (VB.net) using the Visual Studio Integrated Development Environment (IDE) with the geographic information system library “DotSpatial”. The reporting techniques as vector layers, cross analysis and queries were exploited and combined to develop new thematic and representative maps in order to introduce a conceptual model of groundwater functioning and study their vulnerabilities to pollution.
3.1. Meta-Data Exploration and Processing
The purpose of developing the hydrogeological database was to store the main data and integrate useful information by the hydrogeologist. Data are collected from existing databases and maps as well as through new field measurements. The sources of information can be of various origins. The necessary data (topographic, hydrographic, geologic, soil, piezometric ...) are gathered from three different sources:
· Field surveys for the inventory and exploration of the region;
· Consultation of old research and historical data on the region;
· The vectorization of the different existing maps.
Existing and required data types were examined in order to design the database scheme. Parameters and information were reclassified and regrouped at several times. Many hydrogeological parameters and relationships were analyzed in order to be stored in the database. Maximum information, minimum data redundancy, reduction of storage capacity, and optimum retrievability of data for analysis were the constraints that defined the final scheme [2]. This data management tool allows updating and checking of data in real time, automating data storage, management procedures as well as time-series analysis and the production of detailed hydrogeological maps [28].
3.2. Data Modeling
The conceptual design of database is a process to abstract the user needs to information model [32]. The cartographic modeling is the first step in the geographical information system realization. This function is used to model the real data according to the needs and the desired objectives. Table 1 illustrates the mapping model of the objects for the DSS. Topography is represented by contour lines. Point information describes the groundwater and rainfall layers. The groundwater level measurement, hydrochemistry and permeability layers derived from processed attributes data.
The system design is a creative and an essential phase in the development process. Applying of conceptual modeling aims at considering more adequate requirements in the development of the Information Systems according to their operating environments.
For the DSS conceptual modeling, we chose the UML. In fact, it provides a standard way to visualize a design using diagrams. UML includes a set of graphic notation techniques to create visual models of software-intensive systems [26]. There are two major kinds of UML diagram: static and dynamic diagrams. For the static diagrams, we present the use-case and classes diagrams. For the dynamic diagrams, we present the sequence diagrams.
3.2.1. Use-Case Diagram
The use-case diagram describes a set of actions that the DSS should perform in collaboration with the users (actors) of the system to provide valuable results. Figure 2 illustrates the use-case diagram which is a basic element of object-oriented modeling. The functionality of the system is determined by examining the needs of each actor that will be stated as an interaction family. Moreover, the user of the application will implement the data (water well, sample, analysis, rainfall stations, water tables and maps), select the appropriate parameter,
Table 1. Layers of hydrogeological database.
update the selected parameter (add, modify, delete) and edit the reports (Quality monitoring, piezometric monitoring, water tables and Rainfall). The system displays the selected parameter and saves the information.
3.2.2. Class Diagram
The class diagram is one of the most widely used diagrams from the UML specification [22] Describing the static view of a system using UML is a process that relies on class diagrams [25].
The class diagram captures details about entities of the system and static relationships between them. Figure 3 describes the structure of the class’ diagram of our DSS. In fact, it makes it possible to define the components of the final system. At this phase, every class has a unique name, attributes (label, size, type ...) and all the applicable operations to this class. The data dictionary presents the following classes: Water well, Sample, Analysis type, Analysis, Nature, Lithology, Type of measurement, Rainfall stations, Rainfall, Water table, Historic water table, Map, Governorate, Delegation, Pollution sources and Actors.
3.2.3. Sequence Diagram
The sequence diagram describes how processes operate (order and interactions). It also shows object interactions arranged in time sequence. For each task of the system, the sequence diagram designs the objects and classes scenario and the sequence of messages exchanged between the objects. Figure 4 presents the sequence diagram of a shapefile.
3.3. Data Acquisition (DAQ)
This descriptive component includes the quantitative and qualitative data of the water wells inventoried by the regional office of agriculture development of
Monastir. It includes the field surveys carried out in this framework as well as the climatic, piezometric, geochemical and the pollution sources data. The spatial component is a cartographic database that stores data in which the objects are represented by geometric entities (point, arc and polygon) such as pedology, land use, topographic and geological data. Acquisition of these later is based mainly, on three processes: Scanning, georeferencing and digitization of the topographic and geologic maps.
3.4. Software Implementation and Development
In order to achieve this process, we propose to use both the object-oriented development approach and the georelational models attached to GIS databases. The DSS is developed in Visual Basic (VB.net) using the Visual Studio Integrated Development Environment (IDE) that implements the geographic information system library “DotSpatial”. We had applied the conceptual model in the Database Management System (DBMS). This latter is based on the relational schema which describes the different relationships between the data tables. The georelational model consists of linking a relational database to geometrical features. The modeled entities are organized into categories sharing common characteristics. A table represents each category. The different attributes occur as columns of the table, and the rows assure the data registration. Relationships “one-to-one” or “one-to-many” can be established between tables [33].
The design of semantic data begins firstly with the creation of an SQL database (MS Access) using data entry forms. The Application Menu represents the main menu of the developed application that implements the Windows Ribbon framework (Figure 5). It is formed by three main interfaces namely: “Hydrogeological Information System of Monastir (HISM)”, “Reports” and “GIS Data”. This latter allows the user to view, add, modify, and delete the data of water well, sample, rainfall, water tables and maps. Figure 6 and Figure 7 describes input form for rainfall records and importing and/or exporting digitized layers respectively. The conversion of water points and rainfall stations from geographic into metric coordinates is executed within the DSS (Figure 8). For each piezometer and water drilling, the stratigraphy and lithology are seized for calculated the equivalent permeability vertical for vadose zone and the equivalent permeability horizontal for aquifer (Figure 9).
Figure 6. Input form for rainfall records.
Figure 7. Input form for importing and/or exporting digitized layers.
The “Reports” component grants viewing and updating groundwater level measurement, quality and rainfall over time. Figure 10 shows the report of variation of salinity as a function of time. The interface for the “Water Table” shows water well equivalent permeability of the vadose zone and the aquifer at water drilling with an object identifier “19536/4” (Figure 11).
Figure 12 describes the structure of the dynamic interface of the Hydrogeological Information System (HISM) of Monastir composed of four components. The menu bar enables to add and remove themes, change the layout and print. A toolbar is used to perform spatial analysis such as identification, zoom, move,
Figure 8. Input form for geographic coordinates.
Figure 9. Input form for stratigraphy and lithology at water point.
Figure 10. Report of variation of salinity as a function of time.
Figure 11. Equivalent permeability of the vadose zone and the aquifer at water point.
Figure 12. Dynamic interface of the Hydrogeological Information System of Monastir (HISM).
select, remove selection, attribute table and measure distances. The table of contents helps to add the geographic layers, viewing themes area with their planimetric coordinates.
4. Application of DSS: Case of Moknine Superficial Coastal Aquifer
Three main distinct applications of HISM are recognized: 1) the management of hydrogeological data, 2) hydrogeological maps elaboration, 3) vulnerability assessment based on overlay and index methods.
This DSS was applied to the superficial coastal aquifer system of Moknine. This region is located in the north-eastern part of Tunisia, between the latitudes 3936 and 3951 N and the longitudes 663 and 685 E (Figure 13). The Moknine region is influenced by the dry and hot climate from the south and the Mediterranean humid and temperate climate from the north. The annual average temperature is about 19˚C. Rainfall constitutes the main hydrological factor in the region, with an annual average of 350 mm [34].
4.1. Hydrodynamic Functioning Model
The Moknine superficial coastal aquifer system consists of sandy levels, sandstones and calcareous encrustations separated by sandy-clay levels. There is considerable variability in its hydrodynamic characteristics where the Teboulba region
has high hydraulic gradients. The mapping of piezometry revealed the existence of perched groundwater within the groundwater system of the Moknine region. At the level of Bekalta, flow is from a piezometric dome to the sebkha on one side and the sea on the other. Monitoring of the piezometry history at the Teboulba region highlighted the impact of artificial recharge and the behavior of perched groundwater within the phreatic system.
The advanced groundwater analysis in the Moknine region has shown that they are mostly brackish. The salty waters characterize the coastal area of Bekalta and the northwest of the sebkha of Moknine. Freshwater was identified (0.6 to 2 g/l) in the Teboulba area and east of the Moknine sebkha. In the upstream part, and south of Moknine, it is essentially the dissolution/precipitation of minerals that dominate. The mineralization of the waters in the regions of Ksibet elMediouni-Lamta-bouhjar-KsarHellal-Moknine and south of the city of Teboulba results from mixtures, mainly by the contribution of domestic sewage and seawater intrusion. The mixing and ion exchange process are responsible for the mineralization of the waters at the western edge of the Moknine Sebkha. North and east of Bekalta, there is a mix with marine waters.
Figure 14 illustrates Hydrodynamic functioning model of Moknine superficial coastal aquifer system. This model shows three compartments of different characteristics. At the Ksibet ElMediouni-Lamta-Bouhjar-Sayada and Fadhlin-Al Bhira compartments, the reservoir has a reduced thickness and sodium chloride facies. In the north-west, local recharge leads to a dilution of waters where the salinity varies from 1.6 to 3 g/l. To the east, in the region of Al Fadlin-Al Bhira, salinity is high (6 g/l). In the Teboulba region, the model reveals the presence of good quality perched aquifers and reservoir thickening with mixed Na-Ca-Cl and Na-HCO3 facies.
Figure 14. Hydrodynamic functioning model of Moknine superficial coastal aquifer system.
4.2. Groundwater Vulnerability Mapping
4.2.1. GALDIT Method
The GALDIT indexing method was developed by Chachadi and Lobo-Ferreira (2001) [35] for vulnerability assessment to seawater intrusion. This method uses a set of six parameters influencing the seawater intrusion in the aquifer: Groundwater occurrence, Aquifer hydraulic conductivity, height of groundwater Level above sea level, Distance from the shore, Impact of existing status of seawater intrusion and Thickness of aquifer being mapped. In this model, spatial datasets on depth of groundwater, aquifer thickness, aquifer type, hydrochemical analysis and the aquifer hydraulic conductivity are combined.
The GALDIT model was applied to evaluate the groundwater vulnerability of the Moknine phreatic aquifer system. Calculation of the GALDIT index and can be expressed as follows [36] :
The subscripts Ri and Wi denote the rating and the weight for each parameter, respectively.
4.2.2. Development of Thematic Maps
The aquifer characteristics of the study area were carried out by examining exploratory data of 9 wells. The aquifer is mainly unconfined and the rating of aquifer type corresponds to 7.5.
The estimate of the aquifer hydraulic conductivity has been made using calculations of equivalent permeability [37]. Based on the reference permeability coefficients assigned for each lithology [38] and the formula given by Castany (1967), the horizontal equivalent permeability is calculated for each piezometer. The hydraulic conductivity was classified into three classes and was rated 7.5 - 10 (Table 2).
The height of groundwater above mean sea level (L) is calculated from groundwater well measurements during 2005. The values ranged from −12 m to 23 m.
According to the classification of Chachadi and Labo-Ferraira (2001) [25] , for the distance from the shore, the rating attributed varies from 2.5 to 10.
Table 2. GALDIT parameter weights and ranges for Moknine superficial coastal aquifer system.
Geochemistry prospecting techniques and hydrogeochimical analysis are required to determine impact of existing status of seawater intrusion which can be computed using the ratio of (Cl/HCO + CO). Three classes are recorded for this seawater intrusion indicator (Table 2).
All the parameter layers were then superposed in HISM environment to produce the GALDIT map which identifies areas more likely to be susceptible to seawater intrusion.
GALDIT vulnerability map of the Moknine superficial coastal aquifer (Figure 15) shows three degrees of vulnerability; low, moderate and high with 39.9%, 43.8% and 16.3% of the area, respectively. Except for the Lamta-Sayada-Bouhjar region, the coastal zone is the most susceptible and the most exposed to seawater intrusion. This high vulnerability can be explained by the limited distance from the coast and the low piezometric levels.
5. Conclusions
The present research emphasizes the importance of developing a Decision Support System (DSS) as a helpful and excellent decision-making tool for groundwater management.
The design and implementation of this DSS called Hydrogeological Information system for Monastir region (HISM) are a new concept based on the Object Oriented approach (OO) whose:
· The modeling language is the Unified Modeling Language (UML);
· The programming language is Visual Basic (VB.net) in an environment of Visual Studio integrated development;
· The “Dotspatial” library is used for the management of geographical information layers;
· The database management system (DBMS) is relational: Microsoft Access.
A relational hydrogeological database coupled to a Hydrogeological Information System (HISM) allows the rapid and efficient processing of large volumes of
Figure 15. Vulnerability to seawater intrusion degree for Moknine superficial coastal aquifer system.
spatial data from multiple sources. The interactive interface (HISM) is the link between the two components supported by the application. It ensures the visualization of the various geographical layers, the storage and update of the data, the change of scale and the restitution of thematic maps.
Through spatial and attribute requests, the DSS permits the display and analysis of geographic information layers to define the hydrodynamic functioning model of aquifer systems and mapping their vulnerability to seawater intrusion.
Indeed, the compilation and the synthesis of geological, topographical, structural and piezometric data allowed the construction of the hydrodynamic conceptual model of the groundwater system in the study area. It is characterized by a lenticular geometry where permeable levels are inserted in a semi-permeable mass with locally perched hydrogeological entities. The hydrochemical approach enabled the development of thematic maps as well as the origins and processes responsible for groundwater quality. The mineralization of the groundwater system is controlled by the dissolution of evaporating and carbonates minerals, mixing reactions, ion exchange and seawater intrusion. The contribution of marine waters to the groundwater quality is manifested in the northeast of Bekalta where the salinity is greater than 6 g/l.
The groundwater vulnerability of the Moknine is evaluated using the GALDIT empirical model. This method yields six hydrogeologic parameters, which are Groundwater occurrence, Aquifer hydraulic conductivity, height of groundwater Level above sea level, Distance from the shore, Impact of existing status of seawater intrusion and Thickness of aquifer being mapped. The northeastern part of the study area is the most vulnerable area.
In order to protect the water resource of the study area some recommendations can be proposed:
· A systematic inventory along the northern and eastern coastal zone to monitor the spatial and vertical salinity evolution;
· Setting up of a surveillance and warning system in strong vulnerable areas and at high risk of contamination, mainly the Bekalta region.