# Geoelectrical Surveys for Characterization of the Coastal Saltwater Intrusion in Metapontum Forest Reserve (Southern Italy)

Authors: Antonio Satriani; Antonio Loperte; Vito Imbrenda; Vincenzo Lapenna
Journal: International Journal of Geophysics
 (2012)

---

## Abstract

A geoelectrical survey was carried out in the Metapontum Forest Reserve located along the Ionian coast of the Basilicata region (Southern Italy). In this work we used the method of two-dimensional electrical resistivity tomography for obtaining high-resolution electrical images in the investigated site. In particular, three electrical resistivity tomography, all orthogonal to the coastline, in the investigated area were carried out. To complete and integrate the geophysical data, soil and groundwater samplings, seventeen and five, respectively, were analyzed using chemical physical techniques. Geoelectrical survey, supported by laboratory analysis of soil and water samples have revealed the presence of a process of saltwater in coastal Forest Reserve of Metapontum, which have caused the decline of the existing pine forest with the consequent erosion and desertification problems. The results have disclosed the way to identify and discriminate large areas affected by intensive soil salinization and high resolution electrical images of the subsurface electrical resistivity plays a key role in delineating the saltwater intrusion front in coastal areas. Furthermore, our integrated study represents a contribution to the future programs for the protection, planning, and management of the terrestrial and marine resources in this coastal area.


---

## Body

## 1. Introduction

Coastal areas are of great environmental, economic, social, and cultural relevance. Therefore, the implementation of suitable monitoring and protection actions is fundamental for their preservation and for assuring future use of this resource. Such actions have to be based on an ecosystem perspective for preserving coastal environment integrity and functioning and for planning sustainable resource management of both marine and terrestrial components (EU Recommendation on Integrated Coastal Zone Management—ICZM—European Commission EC/413/2002).Post [1] has defined coastal aquifers as the subsurface equivalents of coastal areas where continental fresh groundwater and seawater meet. Coastal plains are often contaminated by salt waters and the process associated to the marine water entering an aquifer is generally called seawater intrusion. Different approaches have been adopted to estimate seawater intrusion. For example, many authors [2–5] have employed geochemical methods based on measures of electrical conductivity, of chloride concentration, and other cation and anion concentrations, parameters which generally highlight seawater contamination.Nowadays great attention is focused on innovative geophysical methods for solving hydrogeological and environmental problems [6–10]. The geophysical methods for the subsurface exploration can contribute to the different phases of coastal environment characterization. Balia et al. [11] have applied a combination of geophysical techniques at the Muravera plain in southern Sardinia, Italy, which has been significantly affected by water salinization; Batayneh [12] used electrical resistivity methods for detecting subsurface fresh and saline water in the eastern Dead Sea coastal aquifers in Jordan; Bauer et al. [13] applied electrical resistivity tomography (ERT) to map the salinity distribution in the subsurface in the Okavango Delta, a large inland delta in Botswana; while Koukadaki et al. [14] characterized the hydrogeology of a coastal aquifer using geophysical measurements and, in particular, electrical tomography. ERT, moreover, was applied to monitor seawater intrusion by Kruse et al. [15], Nowroozi et al. [16], Abdul Nassir et al. [17], and Sherif et al. [18]. Cimino et al. [19] carried out vertical electrical sounding in the coastal plain of Acquedolci (Northern Sicily- Italy) to assess seawater intrusion, while Khalil [20] applied a direct current resistivity geoelectric technique to delineate salt water intrusion from the Gulf of Suez. Moreover, Urish and Frohlich [21], Ebraheem et al. [22], and Choudhury and Saha [23] have performed geoelectrical surveys in different coastal area around the world.In this paper we present the results of geoelectrical surveys, integrated with geochemical measurements, performed for the characterization of a coastal area along the Ionian side of the Basilicata region (Figure1). This study area, located at the mouth of Bradano river (Metapontum, 40°22′ N, 16°51′ E), is mainly interested by a narrow shore (10–30 m) of fine sandy formations and has been affected by saltwater intrusion phenomena. This has caused the decline of the existing pine forest entailing significant ecological and environmental problems to the Metapontum Forest Reserve. Primary and secondary salinity of soil are linked to natural factors such as the climate, natural drainage patterns, topographic features, geological structure, distance to the sea, and by the indiscriminate exploitation of the soil and groundwater resources [1, 24, 25].Figure 1
Location of the study area. The yellow ellipse highlights the natural reserve of Metapontum.In order to manage rightly this kind of problem, it is necessary to have a thorough knowledge of the phenomenon, and the first objective of this study is to identify, through the use of different techniques and methodology, salt intake as the factor concurring or even determining the serious phenomenon of total decline and the deterioration of vegetation in the pine forest.The geoelectrical investigations and chemical-physical analysis were carried out, to evaluate the seawater intrusion effect. In particular, geoelectrical surveys allowed to assess the extent and depth of the saltwater intrusion process, while the chemical-physical characterization was performed by considering some soil properties such as pH values, salt contents measured as the electric conductivity (ECe), exchangeable sodium, and cation exchange capacity. Moreover, this study aims to better delineate the saltwater intrusion front and to support the environmental protection and regeneration programs carried out by local and national authorities for restoring the continuity of the woodland belt along the coast and consequently to obtain a reclassification of the landscape and protection of areas behind.


## 2. Study Area

The study area is located along the Ionian coastal plain, known as Metapontum plain. The area is made up of the regressive filling with marine sediments of a tectonic graben (Bradanic graben) during Pliocene to Pleistocene age. The area is bordered in the NE by the limestones of the Murgian foreland and in the SW by the clastic sediments of the Apennine chain. In the middle of the graben the full sedimentary sequence starts with sands, blue clays, covered by sands and gravels. The alluvial deposits in the estuarine environment are constituted of stratified sands to silty clays sediments. In this case, they are grey and very similar for aspect, but not for consistency, to the Pleistocene blue clays [4, 26, 27]. In particular, the Forest Reserve with an extent of about 240 hectares, situated in the north of the Ionian coastline between the mouth of rivers Bradano and Basento, is characterized by a superficial sandy substrate, and by clays and silts of river origin in the back dune zone (Figure 2).Figure 2
Geological map of the study area. (1) Current alluvional deposits; (2) dunes and beach sands; (3) late alluvional deposits sandy and conglomerates; (4) marine terraces; (5) gray marly clays.The hydrogeology of the study area (Figure3) is related to lithological and structural features of the same area, and they determine two types of aquifers.Schematic lithological section [4]. Legend: (1) soil; (2) clays or silty clays (yellow, brown, grey; (a)) locally sandy clays (b); (3) pebbles (a) in a sandy (b) or clayey matrix (c); pebbles locally cemented (d); (4) grey sands with clayey strata; (5) grey or yellow sands and silty sands (a), locally clayey sands (b) or with gravels (c), locally sandstone strata (d); (6) piezometric level.
(a)(b)An aquifer is constituted by marine terraces deposits and alluvial river valleys deposits, where the marine terraces aquifers show hydraulic conductivity from medium to high, equal or higher than the second type of aquifer constituted by coastal plain deposits.The bottom of the this aquifer is a silty-clayey bed which lies under the sea level permitting the seawater intrusion, this phenomenon involves the coastal plain for a width of 1–1.5 km and it is less evident moving inwards where the altitude of the clayey bottom of the aquifer becomes progressively higher than the sea level [4].The Metapontum coast is characterized by a strong erosive process. Marine erosion has really removed the dune behind the beach moving in the inland and causing the decline of the historical pinewood. Dunes and their vegetation are a natural defence against erosion by wind and sea and the variation of the line of the seashore is in close connection with the development of coastal dunes and their preservation [28]. Moreover, reduction of the coastal dune system in the study area facilitates the entry of salt water; dunes in fact, for their elevation and good infiltration capacity, give a sufficient freshwater recharge and a pressure above sea level, allowing hydrostatic control of saline intrusion [29].Metapontum wooded strip is an artificial formation, planted in the first decade of the 50s in order to preserve the coast and the inland cultivated areas from salt sea wind. The forest is composed of Aleppo pine trees (Pinus halepensis) and domestic pine trees (Pinus pinea) prevalently. Secondary species are acacia saligna (Acacia cianophylla) and eucalyptus (Eucalyptus globulus, E. camaldulensis), besides other species that are typical of Mediterranean scrub as Pistacia lentiscus, Juniperus macrocarpa and Juncos acutus, often mixed with the type of vegetation which prefers sandy soil (called “psammofila”). The rearrangement of soil was done with the realization of “baulature,” that is, elevation of soil in order to preserve tree roots from seawater intrusion.On the basis of meteorological data, the climate of the area is semiarid with hot and dry summers; however, the evolution of the climate of the last period shows a trend to a reduction of annual precipitation and an increase in average temperatures [30].


## 3. Field Investigations and Sampling

In Figure4 location of geoelectrical surveys and sampling points for chemical-physical analysis was indicated. The surveys have been performed during a geophysical campaigns carried out in July, and in this time of year the investigated area, classified as climatic type semiarid, is characterized by drought.Figure 4
Location map of geophysical surveys and sampling points.


### 3.1. Geoelectrical Surveys

Resistivity measurements were conducted as tomographies to determine variations with depth in soil resistivity.The surface electrical resistivity tomographies are a useful tool in determining seawater intrusion in coastal areas for its capability to discriminate the large resistivity contrast between the presence of seawater that strongly reduces the resistivity values and saturated freshwater layers. Moreover, in order to evaluate if this contrast is significant and does not depend on geological layer, a validation of the geoelectrical measurements with some chemical-physical analysis of soil samples has been performed.Three electrical resistivity tomographies (ERT) were obtained using the georesistivity meter Syscal R2 (Iris Instruments) connected to a multielectrode system. We used a Wenner-Schlumberger configuration array with electrodes spacing on the surface of 10.00 meters. The electrical tomographies were all orthogonal to the coastline.The Wenner-Schlumberger array provided horizontal and vertical resolution, and to obtain 2D resistivity models the data obtained were interpreted through the inversion algorithm, RES2DINV, proposed by Loke [31]. The measured values of apparent resistivity provide, in fact, a first preliminary image of the electrical subsurface structure denominated as the “pseudosection”. In a second step, the apparent resistivity measurements are transformed into true resistivity values using the rapid inversion algorithm of Loke and Barker [31].The inversion routine is based on a smoothness-constrained least-squares method inversion algorithm that divides the subsurface into rectangular blocks, and the resistivity of the blocks is adjusted to minimize iteratively the difference between the computed and the measured apparent resistivity values; the root mean square (RMS) error gives a measure of this difference [31, 32].


### 3.2. Soil and Water Sampling

In addition, chemical and physical analyses on soil and water samples were performed in order to validate and integrate geophysical data. Undisturbed soil samples, collected from the surface in the first 40 cm of soil using Soil Sample Ring Kits (Eijkelkamp Equipment Netherlands), were taken for chemical analyses, using standard analytical procedures (Official Gazette of Italy No 248 of 21 October 1999). Soil samples were analyzed for pH (in water and in CaCl2), for cations (Na,K,Ca, and Mg), for sulphates, chlorurs and nitrates, for organic matter and organic carbon, for cation exchange capacity, for an electrical conductivity of a saturated paste extract (ECe) or its inverse, electrical resistivity, and for texture. Groundwater samples were collected from few shallow holes across the area, since the groundwater table in the study area was high.


## 3.1. Geoelectrical Surveys

Resistivity measurements were conducted as tomographies to determine variations with depth in soil resistivity.The surface electrical resistivity tomographies are a useful tool in determining seawater intrusion in coastal areas for its capability to discriminate the large resistivity contrast between the presence of seawater that strongly reduces the resistivity values and saturated freshwater layers. Moreover, in order to evaluate if this contrast is significant and does not depend on geological layer, a validation of the geoelectrical measurements with some chemical-physical analysis of soil samples has been performed.Three electrical resistivity tomographies (ERT) were obtained using the georesistivity meter Syscal R2 (Iris Instruments) connected to a multielectrode system. We used a Wenner-Schlumberger configuration array with electrodes spacing on the surface of 10.00 meters. The electrical tomographies were all orthogonal to the coastline.The Wenner-Schlumberger array provided horizontal and vertical resolution, and to obtain 2D resistivity models the data obtained were interpreted through the inversion algorithm, RES2DINV, proposed by Loke [31]. The measured values of apparent resistivity provide, in fact, a first preliminary image of the electrical subsurface structure denominated as the “pseudosection”. In a second step, the apparent resistivity measurements are transformed into true resistivity values using the rapid inversion algorithm of Loke and Barker [31].The inversion routine is based on a smoothness-constrained least-squares method inversion algorithm that divides the subsurface into rectangular blocks, and the resistivity of the blocks is adjusted to minimize iteratively the difference between the computed and the measured apparent resistivity values; the root mean square (RMS) error gives a measure of this difference [31, 32].


## 3.2. Soil and Water Sampling

In addition, chemical and physical analyses on soil and water samples were performed in order to validate and integrate geophysical data. Undisturbed soil samples, collected from the surface in the first 40 cm of soil using Soil Sample Ring Kits (Eijkelkamp Equipment Netherlands), were taken for chemical analyses, using standard analytical procedures (Official Gazette of Italy No 248 of 21 October 1999). Soil samples were analyzed for pH (in water and in CaCl2), for cations (Na,K,Ca, and Mg), for sulphates, chlorurs and nitrates, for organic matter and organic carbon, for cation exchange capacity, for an electrical conductivity of a saturated paste extract (ECe) or its inverse, electrical resistivity, and for texture. Groundwater samples were collected from few shallow holes across the area, since the groundwater table in the study area was high.


## 4. Results

The 2D resistivity results obtained by the inversion algorithm performed with the software RES2DINV are shown as images in Figure5.2D electrical resistivity tomographies obtained in the investigated site. The position of each ERT is shown in Figure3.
(a)(b)(c)Electrical resistivity values are related to geological parameters of the subsurface and, in particular, resistivity values are controlled by the types of rocks and fluid. Then, the high-resolution electrical images are a powerful tool to identify conductive zones for the saltwater intrusion phenomena. Figure5 shows the three 2D electrical images obtained in the study area related with depth and true resistivity of the subsoil investigated.ERT 1 with direction NNW towards SE obtained on a rather flat regular surface, showed a first zone in the shallow part with values of 5–15 Ωm which could be associated with a layer filled with inert material, and a second one with lower resistivity values (ρ<1.5 Ωm) reflecting saturated strata.ERT 2 is nearly parallel to the direction of ERT 1, it runs from NW to SE. As shown in Figure5, the geoelectrical image presents a lateral variation in resistivity distribution, with higher electrical resistivity values (about 1–10 Ωm) towards SE in the first part of the section (between 0 and 110 m). Lower electrical resistivity values were observed towards sea side, and these values reduced are maybe due to saturated strata for a subsurface seawater flow zone.ERT 3 in direction of WNW-SE shows a significant trend of resistivity, in fact, the resistivity section is characterized by two zones with different resistivity values: the first zone, towards the inland of the Forest Reserve and throughout the thickness of soil explored, shows moderate to high resistivity values 10–220 Ωm, while the second zone, towards the sea, shows very low resistivity values.As expected, the 2D electrical resistivity tomographies show that resistivity values decrease in depth and from inland towards the Ionian Sea.We have checked, one by one, the apparent resistivity values obtained during the field survey, and the topographic correction has been introduced. Only the apparent resistivity measurements affected by a standard deviation lower than 5% have been considered. In a second step we have increased the number of iterations of the inversion software routine. After this procedure the RMS of the electrical imaging is lower than 15%. If we take into account the strong variability of the resistivity values from ~1000 Ωm for sandy layers to ~1 Ωm of fully saltwater saturated zones, this error can be considered acceptable.We used the RES2DINV, that is, a well-assessed algorithm for 2D resistivity data inversion, and there is a large scientific literature in which the RMS of the 2D electrical imaging is greater than 10%. Electrical imaging of complex geological environments are generally affected by RMS greater than 15% [10, 33] and in many cases is not possible to reduce the RMS without introducing artifacts.In literature, salt water resistivity values below 1.0 Ωm were reported, in fact, seawater has an average resistivity of 0.2 Ωm [16, 34], while the resistivity of a layer saturated by saline water and dissolved solids is in the range of 8 to 50 Ωm [13, 16, 35, 36]. Therefore, based on these values of resistivity of layers saturated by saline water and some dissolved solids reported in the literature, resistivity data obtained in this work from geoelectric investigations, highlight the presence of strata saturated with seawater.In Tables1 and 2 the results of laboratory analysis for the physical and chemical parameters of soil and groundwater samples were summarised. Laboratory analyses permit to suppose that high pH values are determined by the accumulation of sodium and magnesium on the soil exchange complex. In fact, in the study area, pH values (in CaCl2) have a range between 7.9 and 8.9 in soil samples, and between 7.3 and 7.8 in groundwater samples. A sample of sea water, instead, reported pH value of 8.2. Important differences were resulted in the conductivity (ECe) expressed in Siemens m−1 or its inverse, electrical resistivity, expressed in Ohm m. The mean ECe value on soil samples was equal to 0.613 Siemens m−1, corresponding to 1.631 Ωm with a range from 0.153 to 1.484 Siemens m−1, equal to a range of resistivity from 6.536 to 0.674 Ωm while in water samples the mean ECe value was equal to 2.720 Siemens m−1, corresponding to 0.368 Ωm, and with a range from 0.68 to 7.27 Siemens m−1, equal to a range of resistivity from 1.470 to 0.137 Ωm.Table 1
Minimum and maximum values of chemical-physical data of soil samples collected in the study area.
ParametersMinimum and Maximum ValuesNa (mg kg−1)626–3425K (mg kg−1)162–960Ca (mg kg−1)917–2830Mg (mg kg−1)203–905Nitrogen (g kg−1)0.6–1.4Organic carbon (g kg−1)4.9–21.5Organic matter (g kg−1)8.4–37.1pH (H2O)8.6–9.8pH (CaCl2)7.9–8.9Electrical conductivity (ECe) (Siemens m−1)0.153–1.484Electrical resistivity (Ωm)6.536–0.674Cation exchange capacity (meq/100 gr)7.4–17.2Chlorurs (mg kg−1)162–5894Nitrate (mg kg−1)<0.1Sulphate (mg kg−1)15.3–441Sand %33.3–85.2Silt %7.1–35.7Clay %7.5–38.2Table 2
Electrical proprieties and pH values of groundwater samples.
SamplespHElectrical conductivity (ECe)  Siemens m−1 to 25°CElectrical resistivity    ΩmG17.70.681.470G27.70.811.234G37.81.920.520G47.32.920.342G57.57.270.137Sea water8.25.660.176On the other hand, in seawater samples ECe value was equal to 5.660 Siemens m−1, corresponding to 0.177 Ωm. Furthermore, the major ions were Na+, Ca2+, and Cl−. Therefore, according to the classification NRCS [37], the area investigated presents from inland to coast, zones slightly saline, zones moderately saline, and zones close to highly saline. The high concentration of soluble salts determines high osmotic pressure that makes very difficult to absorb water and nutrients from the roots of plants, creating the limiting factor to the diffusion and survival of forest species. The salinity distribution in the subsurface of study area is reflected by the vegetation cover, in fact, the zone with soil salinity is covered by salt-tolerant grasses [30]. On the basis of the geoelectrical measurements and of the chemical-physical analysis, the presence of zones with very low electrical resistivity values, can be associated to the intrusion of salt water into the subsoil. In the coastal environment, like the study area, hydraulic connection with seawater and inland gradient that move water to inland from seawater source, are such conditions as to cause seawater intrusion [2].


## 5. Conclusions

Geoelectrical surveys reveal the presence of two principal zones characterized by different resistivity values. The former is a strip back dune, closer to the sea, with a higher clay content and located at an altitude below sea level. It is submerged during the winter by water for a poor drainage, and it is characterized by lower resistivity values. The second zone is characterized by higher resistivity values and has a quite good vegetative cover for the best environmental conditions, greater distance from the sea, and elevation above sea level. In this area the spontaneous emergence of new plant pine is frequently observed for the fact that soil conditions are favourable as there is no seawater intrusion yet.Results presented in this paper showed that the process of salt water is active on coastal Forest Reserve of Metapontum, and geoelectrical surveys have proved to be precious tools. Geoelectrical surveys have highlighted significant variations in subsurface resistivity associated with lithological characteristics of underground layers and variations in water saturation. In particular, these surveys have produced images of resistivity of the soil investigated, where the spatial distribution of brackish and saline water in the pine forest and resistivity variation with depth is very clear. These observations are also confirmed from laboratory analysis on soil and groundwater samples.



---

*Source: 238478-2012-02-28.xml*