ERT and GPR Prospecting Applied to Unsaturated and Subwater Analogue Archaeological Site in a Full Scale Laboratory

: Geophysical techniques are widely applied in the archaeological ﬁeld to highlight variations of the physical behaviour of the subsoil due to the presence of ancient and buried remains., Considerable efforts are required to understand the complexity of the relationship between archaeological features and their geophysical response where saturated conditions occur. In the case of lacustrine and wetland scenarios, geophysical contrasts or electromagnetic signal attenuation effects drastically reduce the capabilities of the geophysical methodologies for the detection of structures in such conditions. To identify the capability of the electrical and electromagnetic methods in different water-saturated scenarios, an experimental activity was performed at the Hydrogeosite CNR laboratory. The test allowed us to analyze the limits and potentialities of an innovative approach based on the combined use of the ground-penetrating radar and 2D and 3D electrical resistivity tomographies. Results showed the effectiveness of the ground-penetrating radar for detecting archaeological remains also in quasi-saturated and underwater scenarios despite the em signal attenuation phenomena; whilst the results obtained involving the resistivity tomographies offered a new perspective for the archaeological purposes due to the use of the loop–loop shaped array. Moreover, the radar signal attenuation, resolution and depth of investigation do not allow to fully characterize the archaeological site as in the case of the scenarios with a limited geophysical contrast (i.e., water-saturated and arid scenarios). The experimental tests show that these limits can be only partially mitigated through the integration of the geophysical methodologies and further efforts are necessary for improving the results obtainable with an integrated use of the adopted geophysical methodologies.


Introduction
Among various active geophysical techniques, electric and electromagnetic (em) methods are strongly effective for the detection of archaeological features located in the subsoil at different depths and scenarios [1][2][3][4][5][6][7]. Localization of anthropic structures placed in the soil is possible due to the contrast of the em physical properties between the materials constituting the buried objects and the subsoil where they are "preserved". The presence of the "geophysical contrast" between the analyzed structures and background soil represents the necessary conditions for successful research in the archaeogeophysical field. However, the presence of subsoil with a high-water content (close to the saturation) could be a limit for the research due to em signal attenuation problems and "homogenization effect" of paper defines different GPR and DC conventional and non-conventional acquisition procedures, in order to highlight the best setting in water saturated/quasi-saturated subsoil and flooded conditions, and a first level of integration is shown. In particular, we applied GPR on a boat and ERT on the water to demonstrate the potentiality of the methods applied in tandem for shallow underwater research. Further, the effectiveness of the non-conventional loop-loop array for ERT analyses was analyzed and compared with classical setting. With the loop-shaped acquisition scheme, we have tried to overcome the limitations analysing the effective capacity of the DC method. The qualitative integration of the GPR and ERT data has been tried in order to reduce the uncertainties for the archaeological features detection and preliminary quantitative analyses of the geophysical data regarding the real sizes and positions of the buried objects, which are discussed in this paper.
The outline of this paper is as follows. After a brief overview of GPR and ERT is applied to the archaeological field, a description of the archaeological test site is shown. Then, the results obtained in two main phases of the experiments are presented and discussed. Finally, conclusions and future perspectives are debated in the last paragraph.

Theoretical Notes
Geophysical techniques have become of crucial importance in the preliminary phase of archaeological site detection and mapping; indeed, they are able to provide quick and inexpensive high-quality information on the presence and distribution of remains in many different archaeological scenarios. The probability of a successful application rapidly increases if multi-methodological approaches are adopted, according to a logic of objective complementarity of information and global convergence toward a high quality multi-parametric imaging of the buried structures [41].
In detail, GPR represents an excellent tool to support the archaeologists to identify and reconstruct the real geometry of the objects located in the soil. The GPR lies in electromagnetic (em) theory where the full electromagnetic field is mathematically described by Maxwell's equations. GPR allows for the study of the scattering phenomena of em waves caused by variations of some physical properties of the investigated medium as dielectric permittivity (ε, F/m), electrical conductivity (σ, S/m), and magnetic permeability (µ, H/m). These physical variations generate reflections related to contrasts in the em impedance; further, the em properties of the soil barely influence the velocity of propagation of em waves and the attenuation of the energy introduced in the subsoil. In low-loss conditions, velocity and the attenuation of em waves can be approximated as follows [42]: where c = 1 µ 0 ε 0 is the velocity of light in a vacuum (m/s), µ 0 and ε 0 are the magnetic permeability (4π × 10 −7 H/m) and dielectric permittivity (8.854 × 10 −12 F/m) in the free space and k the dielectric constant k = ε ε 0 (dimensionless) with magnetic permeability considered negligible and ε the dielectric permittivity of the medium (F/m).
The above two expressions show that the dielectric permittivity (ε) controls em wave velocity, while the electrical conductivity (σ) has a large effect on attenuation. For this reason, GPR works well where the soil conductive is low [42][43][44][45].
Direct current (DC) electrical methods are among the oldest and most popular techniques for the non-invasive geophysical investigations. Among them, ERT allows to investigate the horizontal and vertical electrical resistivity variations of the subsurface materials potentially induced by the presence of anomalous bodies. , 12, 1126 DC method is based on the measurement of an electric field artificially created in the ground with suitable electronic devices; it normally consists of two pairs of electrodes fixed in the ground, of which: a pair constitutes the current injection circuit, the other the measuring circuit of the potential difference (dV) generated in the ground by the passage of the current itself. The fundament of this technique is Ohm's Law, which indicates that the potential difference dV (V) at the ends of a conductor, at a given temperature T, is proportional to the electric current I (A) passing through it by means of a quantity constant and typical of the conductor, said resistance R (Ω): where ρ (Ωm) is the electrical resistivity, L (m) is the conductor length, and A (m 2 ) is the conductor area. Generally, a switched square wave is the current waveform used [64]. The data acquired are expressed in form of apparent resistivity (ρ a ) ρ a = ∆V/IK (5) where ∆V (V) is the measured potential, I (A) the transmitted current, and K (m) the geometric factor, which depend on the position of electrodes. It is possible to define different configuration: Wenner, Shlumberger, dipole-dipole, gradient arrays, pole-pole, pole-dipole. Moreover, due to new acquisition instruments with a high number of channels, it is possible to make a personal disposal of the electrodes with different regular or irregular geometry [65]. The apparent resistivity is then interpreted in terms of real resistivity and depth by means of inversion software. The aim of the inversion procedure is to compute the 'best' set of resistivity values, which satisfies both the measured dataset and some a priori constraints, in order to stabilize the inversion and constrain the final image [66,67]. However, the method is affected by limitations mainly due to the resolution and time to perform the investigations. In the past, 2D ERTs were successfully used for archaeological prospecting [68][69][70][71][72][73][74][75][76][77][78]. Furthermore, several authors have effectively performed 3D ERT applied in the archaeological field to detect the excavation of buried structures [1,[78][79][80][81][82][83]. Resistivity methods, including ERT, are less used then GPR in archaeological surveys for the limited resolution, but in the presence of conductive context as in the case of soils with high content of water, salt water or clay), they could be a good alternative to the GPR that suffers attenuation problems. Furthermore, an integration and a comparison of the results obtained via the two techniques is a preferable approach to support archaeologists and enhance the quality of the interpretation of the geophysical data as showed by the existing literature [1,[84][85][86][87][88].

The Full-Scale Archaeologeophysical Test
The experiment was realized at the Hydrogeosite laboratory (CNR-IMAA Marsico Nuovo -PZ), where a concrete pool of 250 m 3 (12 × 7 × 3 m) is located (Figure 1). The pool is filled with silica sand (95% SiO 2 ) characterized by an average diameter equal to 0.09 mm (very fine sand), a porosity of about 45-50%, and a hydraulic conductivity of about 10 −5 m/s (see Table 1).
The pool is equipped with 17 wells and a hydraulic system (draining ring system) for the reproduction of a phreatic aquifer that allows the variation of groundwater level (Figure 1a).
In this context, an archaeological site was reconstructed and buried by sand and different water content conditions were reproduced in order to simulate wet and lacustrine conditions typical of wetlands even assuming extreme events (for example, submarine landslide or fast coast natural erosion). For this experiment, a limited area of 6 m × 4 m  The pool is equipped with 17 wells and a hydraulic system (draining ring system) for the reproduction of a phreatic aquifer that allows the variation of groundwater level (Figure 1a). In this context, an archaeological site was reconstructed and buried by sand and different water content conditions were reproduced in order to simulate wet and lacustrine conditions typical of wetlands even assuming extreme events (for example, submarine landslide or fast coast natural erosion). For this experiment, a limited area of 6 m × 4 m (Figure 1b) inside the large pool was selected, and the archaeological features were buried until a depth of about 1.50 m from the original ground surface.
After the excavation activities, the area was filled with moist sand to obtain a quite homogeneous filling. We have proceeded layer by layer regularly compacting the sand by several cycles of tap water charge and discharge, from the bottom to the surface. In order to make the reconstruction as faithful as possible to real case, a study on the historical and archaeological area was carried out, taking in account the material used, the construction types, and techniques and methods of materials. The archaeological context was obtained following the dictates reported by ancient sources as those set forth by Vitruvio (De Architectura) and several studies published over time [89][90][91][92][93][94][95][96][97]. We adopted stylistic and structural elements belonging to different historical periods, places, and ethnos that find their coexistence in the middle Republican age and in the late-Imperial phase but even employed in more recent times, such as the capuchin tombs and 'box-grave' composed of large clay tegulae (tiles) used in the IV-III B. C. century but also in the Imperial period, or the opus coementicium (i.e., a hydraulic mortar laid in alternate courses with aggregate) continuously employed from its invention to the present time. In detail, it was involved in the construction of two different areas for the living context and funerary world, as shown in Figure 2b. A marble column, a capuchin burial, and a tiles burial were placed,  After the excavation activities, the area was filled with moist sand to obtain a quite homogeneous filling. We have proceeded layer by layer regularly compacting the sand by several cycles of tap water charge and discharge, from the bottom to the surface. In order to make the reconstruction as faithful as possible to real case, a study on the historical and archaeological area was carried out, taking in account the material used, the construction types, and techniques and methods of materials. The archaeological context was obtained following the dictates reported by ancient sources as those set forth by Vitruvio (De Architectura) and several studies published over time [89][90][91][92][93][94][95][96][97]. We adopted stylistic and structural elements belonging to different historical periods, places, and ethnos that find their coexistence in the middle Republican age and in the late-Imperial phase but even employed in more recent times, such as the capuchin tombs and 'box-grave' composed of large clay tegulae (tiles) used in the IV-III B. C. century but also in the Imperial period, or the opus coementicium (i.e., a hydraulic mortar laid in alternate courses with aggregate) continuously employed from its invention to the present time. In detail, it was involved in the construction of two different areas for the living context and funerary world, as shown in Figure 2b. A marble column, a capuchin burial, and a tiles burial were placed, as shown in Figure 2c along the S1 section (from left to right). Instead, along the S2 section as shown in Figure 2d (from left to right), part of a building and a paved road fragment were constructed. The building, defined by two dividing walls converging to form a corner, had incorporated a floor covered with mosaic; in this case, a small collapse near the structure   WP2: all the site was covered by water to reproduce an underwater archaeological site ( Figure 3b).
1. WP1: top of the shallower structures placed at 0.30 m from the original surface (Fig-ure 3a); 2. WP2: all the site was covered by water to reproduce an underwater archaeological site (Figure 3b). In WP1, the water level (WL) was placed at three different heights equal to 0.20 m (WL1), 1.00 m (WL2), and 1.50 m (WL3) from the surface of investigation (Figure 3a), while in WP2 the surveys were performed in the presence of a water column of 0.3 m above the sand surface (Figure 3b). The entire experiment was realized with the use of tap water characterized by the electrical resistivity value of about 30 Ωm.

Geophysical Surveys
Geophysical surveys were performed with the use of the SIR-3000 TerraSurveyor (GSSI System) ground penetrating radar ( Figure 4a) and Syscal Pro Switch 96 (Iris Instruments) georesistivimeter ( Figure 4b). Data acquisition and processing were set considering the two different scenarios related to WP1 (Figure 4a,b) and WP2 (Figure 4c,d).
In detail, GPR surveys were carried out in WP1 using a reference grid where the distance from each line was equal to 0.20 m and the investigations were made in both the main directions (red arrows in Figure 5a). In WP2, the presence of the water did not allow the 3D acquisition and only the two sections S1 and S2 were investigated (Figure 5c). In this case, the antenna was placed on a small boat dragged with pulleys fixed at the concrete wall of the pool to allow for the collection of the data without the use of an odometer. A marker every 0.075 m was imposed for assigning the real position to the acquired traces. All the data were acquired in continuous and reflection mode with the two-time window of 70 ns at the frequency of 400 MHz. The scan samples were set at 512 with a resolution of 16 bits and a transmission rate of 100 KHz.

Geophysical Surveys
Geophysical surveys were performed with the use of the SIR-3000 TerraSurveyor (GSSI System) ground penetrating radar ( Figure 4a) and Syscal Pro Switch 96 (Iris Instruments) georesistivimeter (Figure 4b). Data acquisition and processing were set considering the two different scenarios related to WP1 (Figure 4a,b) and WP2 (Figure 4c,d). In detail, GPR surveys were carried out in WP1 using a reference grid where the distance from each line was equal to 0.20 m and the investigations were made in both the main directions (red arrows in Figure 5a). In WP2, the presence of the water did not allow the 3D acquisition and only the two sections S1 and S2 were investigated ( Figure 5c). In this case, the antenna was placed on a small boat dragged with pulleys fixed at the concrete wall of the pool to allow for the collection of the data without the use of an odometer. A marker every 0.075 m was imposed for assigning the real position to the acquired traces. All the data were acquired in continuous and reflection mode with the two-time window of 70 ns at the frequency of 400 MHz. The scan samples were set at 512 with a resolution of 16 bits and a transmission rate of 100 KHz.
GPR data processing consisted of the use of some basic operations, including bandpass filters, gain function removal, amplitude compensation, background removal filter, and Kirchhoff migration. To migrate the data, the em velocity evaluation was performed, studying the diffraction hyperbolas generated by the archaeological features that was constantly varied due to the raising of the water table. In order to detect and reconstruct the archaeological remains, after the examination of the individual radargrams, 3D data volumes were created and some significative time slices, converted in depth slices, were extracted to identify the buried structures. In WP2, it has been necessary to add a further GPR data processing consisted of the use of some basic operations, including bandpass filters, gain function removal, amplitude compensation, background removal filter, and Kirchhoff migration. To migrate the data, the em velocity evaluation was performed, studying the diffraction hyperbolas generated by the archaeological features that was constantly varied due to the raising of the water table. In order to detect and reconstruct the archaeological remains, after the examination of the individual radargrams, 3D data volumes were created and some significative time slices, converted in depth slices, were extracted to identify the buried structures. In WP2, it has been necessary to add a further step for the data editing, including the association of the real position to the traces acquired without the use of the odometer.
Regarding the ERT acquisitions in WP1, 3D ERTs were performed with the use of two different types of 3D arrays. In detail, a 3D acquisition based on a grid of 96 steel electrodes (pink circles in Figure 5a) distributed on an area of 7 × 4 m was adopted. First, 3D pseudo resistivity data were collected using dipole-dipole array with in-line, parallel-line and diagonal-line dipoles; the electrode spacing was 0.60 m in both the directions. Further, a 3D non-conventional array with a loop shaped (loop-loop array) was used: both current (C+ and C−) and potential (V+ and V−) electrodes were distributed along the entire site surface according to four concentric rings (Figure 5b). In the outer ring (red circles), the electrodes were placed at a reciprocal distance of 1.00 m; in the second ring (magenta circles), the mutual distance was 0.75 m; for the last two inner rings, the distances were 0.50 (green circles) and 0.25 m (blue circles), respectively.
In WP2, where underwater scenarios have been simulated, 2D ERTs were carried out with electrodes placed both on the underwater floor and on the water surface by a floating system, which was realized ad hoc. In addition, only the S1 and S2 sections were investigated ( Figure 5c). resistivity data inversion of 3D ERT was performed using ERTLab software using a quadrilateral mesh. The inversion procedure is based on a smoothness constrained leastsquared algorithm with Tikhonov model regularization, where the condition of the minimum roughness of the model is used as a stabilizing function. [98]. The 2D inversion was carried out with the ResIPy open-source software [99]. In this case, a confined rectangular mesh was used for considering the presence of the concrete walls of the pool. Moreover, in WP2, the presence of the water column has necessarily required the optimization of the starting forward model by considering, above the sand body, the presence of a 30 cm thick layer with constant electrical resistivity value equal to 30 Ωm.  ERT data were processed by means of an inverse modelling software according to an iterative process, which aims at minimizing the difference between the measured pseudosection and the calculated pseudo-section based on a starting model. In detail, apparent resistivity data inversion of 3D ERT was performed using ERTLab software using a quadrilateral mesh. The inversion procedure is based on a smoothness constrained least-squared algorithm with Tikhonov model regularization, where the condition of the minimum roughness of the model is used as a stabilizing function. [98]. The 2D inversion was carried out with the ResIPy open-source software [99]. In this case, a confined rectangular mesh was used for considering the presence of the concrete walls of the pool. Moreover, in WP2, the presence of the water column has necessarily required the optimization of the starting forward model by considering, above the sand body, the presence of a 30 cm thick layer with constant electrical resistivity value equal to 30 Ωm.  Figure 6). Regarding the size of the highly reflective bodies, there is a good agreement between the reflective areas and the presence of the structures, especially for the column and the wall. Some difficulties are related to the localization of the capuchin tomb (F) and enchytrismos (E) that are not clearly detectable for to their geometrical sizes and orientation. Further some reflections not associable to the structures, due to a low soil compaction, have caused a blurry image of the buried objects.
responding to the top of the column and stone wall. The drier scenario (WL3), with the water table deep 1.50 m, is characterized by the em velocity of 0.17 mns −1 equivalent to a εr of 2.8. The reflections imputable to the column (A), stone wall (D), rectangular tomb (B), and paved road (C) are clearly detectable ( Figure 6). Regarding the size of the highly reflective bodies, there is a good agreement between the reflective areas and the presence of the structures, especially for the column and the wall. Some difficulties are related to the localization of the capuchin tomb (F) and enchytrismos (E) that are not clearly detectable for to their geometrical sizes and orientation. Further some reflections not associable to the structures, due to a low soil compaction, have caused a blurry image of the buried objects.  Figure 7. It is interesting to note how the reflections related to the objects are weakened at the depth of one meter for the attenuation of the signal. The capuchin tomb (F) is not easily recognizable, while the enchytrismos (E) is not identifiable. The rectangular tomb (B), the paved road (C), and the wall corner (D) are detectable; moreover, the last structure is well-defined and focused, because the increase of the water content degree has induced a better compaction of the  Figure 7. It is interesting to note how the reflections related to the objects are weakened at the depth of one meter for the attenuation of the signal. The capuchin tomb (F) is not easily recognizable, while the enchytrismos (E) is not identifiable. The rectangular tomb (B), the paved road (C), and the wall corner (D) are detectable; moreover, the last structure is well-defined and focused, because the increase of the water content degree has induced a better compaction of the ground and some reflections caused by the excavation activities, visible in the previous case, disappear.
The third step was obtained with the water table level fixed at 0.20 m (WL1). This is the wetter scenario characterized by a em velocity of 0.07 mns −1 equivalent to a ε r of 18.5. The results shown in Figure 8 highlight unequivocally the attenuation imputable to the fact that the objects are all underwater. Indeed, only the column (A), the top of the wall (D), the rectangular tomb (B), and the enchytrismos (E) are clearly detectable, whilst the capuchin tomb (F) and the paved road (C) are not easily detectable.

ERT Results in WP1
As shown in Figure 9, the resistivity values obtained with the use of the 3D-ERTs range between 2 and 200 Ωm. The low values of resistivity are due to the presence of a high-water content in the sand while the greater values are induced by the presence of the archaeological remains. The resistivity distribution acquired with the classical 3D ERT configuration is not associated well to the buried archaeological structures, both in the drier scenario ( Figure 9a) and the wetter scenario ( Figure 9b). On the contrary, the loop-loop array provided an enhanced information regarding their distribution and shape ( Figure 9c,d). In fact, several resistivity anomalies are well detected in the loop-loop 3D ERT image and they clearly highlight the presence of the corner wall (D), the capuchin tomb (F), and enchytrismos (E), as shown in the wetter scenario WL2 (Figure 9d). Less evident is the presence of the other structures that are located at the edges of the investigated area (column A and rectangular tomb B) or at a greater depth, as in the case of the paved road (C). Appl   The comparison of the results obtained with the two approaches shows the importance of the results obtained with the less conventional loop-loop shaped array, also highlighted by the depth slices extracted in the three different scenarios (Figures 10-12). In detail, the slices extracted at increasing depths from the resistivity model in WL3 (drier scenario) allow to detect at least five of the six structures buried in the subsoil ( Figure 10) despite the apparent noise in the shallower layer of the subsoil, induced by the excavation activities and bad compaction of the soil (Figure 10a). Increasing the depth, the presence of the structures is more easily detectable, as in the case of the corner wall (D) and the paved road (C) (Figure 10b,c). At the depth of 1.00 m, no structures are detectable for the scarce resolution of the method imputable to the few measurement points available.
The results obtained with the water table at the depth of 1.00 m (in WL2) allow to detect a lower noise in the shallowed layers due to the higher water content in the subsoil ( Figure 11), resulting in an improved geophysical contrast between the archaeological structures (more resistive) and background soil (more conductive). As in the previous case, apart from the rectangular tomb (B), the electrical anomalies can be associated for position and shape to the buried structures at the depths of 0.25 m (Figure 11a The third step was obtained with the water table level fixed at 0.20 m (WL1). This is the wetter scenario characterized by a em velocity of 0.07 mns −1 equivalent to a εr of 18.5. The results shown in Figure 8 highlight unequivocally the attenuation imputable to the fact that the objects are all underwater. Indeed, only the column (A), the top of the wall (D), the rectangular tomb (B), and the enchytrismos (E) are clearly detectable, whilst the capuchin tomb (F) and the paved road (C) are not easily detectable.  The comparison of the results obtained with the two approaches shows the importance of the results obtained with the less conventional loop-loop shaped array, also highlighted by the depth slices extracted in the three different scenarios (Figures 10-12). In detail, the slices extracted at increasing depths from the resistivity model in WL3 (drier scenario) allow to detect at least five of the six structures buried in the subsoil ( Figure 10) The results obtained with the water table at the depth of 1.00 m (in WL2) allow to detect a lower noise in the shallowed layers due to the higher water content in the subsoil (Figure 11), resulting in an improved geophysical contrast between the archaeological structures (more resistive) and background soil (more conductive). As in the previous case, apart from the rectangular tomb (B), the electrical anomalies can be associated for position and shape to the buried structures at the depths of 0.25 m (Figure 11a    In WL1 (Figure 12), the high water content does not permit the identification of electrical anomalies associated to the buried structures in any slices.

WP2-Underwater Scenario
In this phase, the underwater archaeological scenario has been analyzed. The experiment was setup by simulating a water column of 30 cm above the sand surface. In correspondence of the S1 and S2 transects, 2D GPR was performed above the water surface while ERTs were carried out above and below the water column. In WL1 (Figure 12), the high water content does not permit the identification of trical anomalies associated to the buried structures in any slices.

WP2-Underwater Scenario
In this phase, the underwater archaeological scenario has been analyzed. The ex iment was setup by simulating a water column of 30 cm above the sand surfac GPR results are shown in Figure 13. In this case, the signal attenuation is greater than the previous conditions (WP1). However, it is possible to easily detect the location and the shape of the buried archaeological structures. The em velocity was equal to 0.035 mns −1 for the water layer (ε r-w ≈ 81) and 0.065 mns −1 for the water saturated sand (ε r ≈ 25). Appl correspondence of the S1 and S2 transects, 2D GPR was performed above the water surface while ERTs were carried out above and below the water column. GPR results are shown in Figure 13. In this case, the signal attenuation is greater than the previous conditions (WP1). However, it is possible to easily detect the location and the shape of the buried archaeological structures. The em velocity was equal to 0.035 mns −1 for the water layer (εr-w ≈ 81) and 0.065 mns −1 for the water saturated sand (εr ≈ 25).
Even with the presence of water-saturated soil and the water column, the reflections imputable to the column (A), the stone wall (D), and the top of the two tombs (F and D) are clearly individuated (Figure 13a,b) while the top of the road (C) is weakly identifiable.

Discussion
Electrical DC and electromagnetic (GPR) techniques were applied to an archaeological framework reconstructed in a laboratory to test the capability of geophysics to identify some objects characterized by different materials (stone, bricks, mortar, etc.) and shapes (rectangular, circular, irregular, etc.) placed in lacustrine and very wet scenarios.
The aim of the experiments was to analyze the limits of the geophysical techniques, as GPR and ERT, when the archaeological structures are below the water level or close to the water table (high water content). The experiment was realized in two distinct phases: in the first one (WP1), ERT and GPR analyses were performed in different scenarios characterized by an increasing water content (WL1-WL2-WL3). The surveys were carried out adopting 3D GPR and ERT. The results have demonstrated the limits of the two methodologies for archaeological issues when the remains are below the water table (high saturation degree).
Regarding GPR surveys, in all the analyzed cases, the influence of the water on the propagation of the em pulses was evident. In Figure 15, the results obtained at the frequency of 400 MHz in the three different scenarios analyzed in WP1 are compared. A strong decay of the radar amplitude signals with the increasing of the groundwater level is clearly recorded. The 2D ERTs were acquired with two different approaches. The electrical resistivity distribution obtained with floating electrodes was not able to detect anomalies well that were associated with the buried objects (Figure 14a,b). It is possible to identify the shallowest features (column A, the wall D) whilst, the features at greater depth, as in the case of the paved road and tombs, are not visible. On the contrary, the electrical resistivity image obtained with electrodes placed on the underwater floor provides an improvement of the information about the position of the archaeological remains (Figure 14c,d). In fact, the shallowest (column A, the wall D) and deeper (paved road and tombs) features are highlighted well.

Discussion
Electrical DC and electromagnetic (GPR) techniques were applied to an archaeological framework reconstructed in a laboratory to test the capability of geophysics to identify some objects characterized by different materials (stone, bricks, mortar, etc.) and shapes (rectangular, circular, irregular, etc.) placed in lacustrine and very wet scenarios.
The aim of the experiments was to analyze the limits of the geophysical techniques, as GPR and ERT, when the archaeological structures are below the water level or close to the water table (high water content). The experiment was realized in two distinct phases: in the first one (WP1), ERT and GPR analyses were performed in different scenarios characterized by an increasing water content (WL1-WL2-WL3). The surveys were carried out adopting 3D GPR and ERT. The results have demonstrated the limits of the two methodologies for archaeological issues when the remains are below the water table (high saturation degree).
Regarding GPR surveys, in all the analyzed cases, the influence of the water on the propagation of the em pulses was evident. In Figure 15, the results obtained at the frequency of 400 MHz in the three different scenarios analyzed in WP1 are compared. A strong decay of the radar amplitude signals with the increasing of the groundwater level is clearly recorded. Further in the first two cases (WL2 and WL3), characterized by a lower water content in the vadose zone, the electromagnetic reflections allow to identify all the buried structures (Figure 15c-f) while in the full saturated case, the GPR anomalies are less evident (15g-h). Figure 16 shows the frequency spectrum of the mean trace for the radargram acquired above the wall and road (S2). In WL2 and WL3 (Figure 16a,b) the signal spectra are almost identical with a limited downshift for the wetter case (Figure 16b). On the contrary, the spectrum in WL1 (Figure 16c) suffers a substantial variation towards the lower frequency where the biggest variation is evident, and the archaeological structures presents are more difficulty to be identified. The frequency centroids estimated for the three cases were equals to 484 MHz (WL3), 406 MHz (WL2) and 348 (WL1) in agreement with the expected signal attenuation in a lossy medium which is affected by absorption and dispersion phenomena able to cause a decay of their amplitude [100]. Further in the first two cases (WL2 and WL3), characterized by a lower water content in the vadose zone, the electromagnetic reflections allow to identify all the buried structures (Figure 15c-f) while in the full saturated case, the GPR anomalies are less evident (15g-h). Figure 16 shows the frequency spectrum of the mean trace for the radargram acquired above the wall and road (S2). In WL2 and WL3 (Figure 16a,b) the signal spectra are almost identical with a limited downshift for the wetter case (Figure 16b). On the contrary, the spectrum in WL1 (Figure 16c) suffers a substantial variation towards the lower frequency where the biggest variation is evident, and the archaeological structures presents are more difficulty to be identified. The frequency centroids estimated for the three cases were equals to 484 MHz (WL3), 406 MHz (WL2) and 348 (WL1) in agreement with the expected signal attenuation in a lossy medium which is affected by absorption and dispersion phenomena able to cause a decay of their amplitude [100].
However, the 3D reconstructions obtained with the GPR data have demonstrated a great capability to identify all the structures in WL2 and WL3; in WL1, where the attenuation is greater, it was impossible to detect those objects placed at the greater depths or not perpendicular to the adopted acquisition grid. It is worth noting that some difficulties are also encountered with the drier scenario (WL3) where the results are noisier because the inadequate compaction of the soil and the scarce permittivity contrast did not allow to reconstruct the real shape of the archaeological structure, that are better focused in WL2 where the water content is higher; therefore, the expected geophysical contrast is higher. However, the 3D reconstructions obtained with the GPR data have demonstrated a great capability to identify all the structures in WL2 and WL3; in WL1, where the attenuation is greater, it was impossible to detect those objects placed at the greater depths or not perpendicular to the adopted acquisition grid. It is worth noting that some difficulties are also encountered with the drier scenario (WL3) where the results are noisier because the inadequate compaction of the soil and the scarce permittivity contrast did not allow to reconstruct the real shape of the archaeological structure, that are better focused in WL2 where the water content is higher; therefore, the expected geophysical contrast is higher.
In WP2 (underwater scenario), GPR showed encouraging results in presence of nonionic solutions giving optimal information about the distribution and geometrical shape of the structures despite the attenuation limits (strong variation of the spectrum frequency characterized by a centroid value of 287 MHz). Thus, GPR confirmed its usefulness for the archaeological research also in lacustrine contexts as support of the well-established sound acoustic techniques allowing to obtain a higher resolution.
Results obtained with the 3D ERTs show how important it is to define the best electrode setting; therefore, the simulated approaches in full-scale laboratory experiments help to evaluate the best solution for each kind of detection aims. In our experiments, the results highlight the capability of the loop-loop array to detect the structures buried in the subsoil is greater than the classical dipolar approach. Therefore, these results show the potentiality of this approach for the investigation in areas where obstacles can reduce the possibility to perform electrical analyses with 2D and 3D ERTs based on conventional approaches. In WP2 (underwater scenario), GPR showed encouraging results in presence of nonionic solutions giving optimal information about the distribution and geometrical shape of the structures despite the attenuation limits (strong variation of the spectrum frequency characterized by a centroid value of 287 MHz). Thus, GPR confirmed its usefulness for the archaeological research also in lacustrine contexts as support of the well-established sound acoustic techniques allowing to obtain a higher resolution.
Results obtained with the 3D ERTs show how important it is to define the best electrode setting; therefore, the simulated approaches in full-scale laboratory experiments help to evaluate the best solution for each kind of detection aims. In our experiments, the results highlight the capability of the loop-loop array to detect the structures buried in the subsoil is greater than the classical dipolar approach. Therefore, these results show the potentiality of this approach for the investigation in areas where obstacles can reduce the possibility to perform electrical analyses with 2D and 3D ERTs based on conventional approaches.
Moreover, 2D ERTs acquired with electrodes located on the water surface show that the high conductive environment and the depth of investigation of the methodology do not allow to easily identify the buried structures. On the contrary, the use of floor embedded electrodes allows to improve the depth of investigation for supporting the detection of the buried structures. Therefore, 2D ERTs can be a suitable contribution for the shallower structures, but the limitations mainly related to the problem of resolution at depth require further efforts to apply in contexts similar to the analyzed ones.
In order to overcome the limits of the single methodology used, the best approach is the integration of electrical resistivity and GPR data. In detail, through a simplified geophysical integration of the GPR and ERT 3D data, based on the model co-rendered image process [101], the integration of the methodologies provides good information for the localization of the archaeological buried remains as shown in Figure 17. 3D ERTs have provided acceptable results only in the drier conditions (WL2 and WL3). The cause is due to the low resolution of the method when the water table was shallower (WL1) and the electrical contrasts were weak. However, the strong reflections are generally correspondent to the higher electrical anomalies, and where one of the two methodology struggles, the other one can support the detection of the structure. This is the case of the capuchin tomb (F), which is not clearly detectable by the GPR data due to its orientation with respect to the acquisition direction, whereas ERT data clearly identify the presence of a strong electrical anomaly where the tomb is expected. On the contrary, the rectangular tomb (B) is well detectable with the GPR, while it is not visible for the ERT because its position too marginal within the acquisition scheme.
Moreover, GPR radargrams and ERTs acquired along S1 and S2 during WP2 (underwater conditions) are model co-rendered and the results are plotted in Figure 18. The first level of integration highlights the capability of both the methods to detect electric and em anomalies in correspondence of the buried structures despite the different resolution. Generally, the presence of strong reflections implies anomalies of the electrical behaviour of the subsoil related to the presence of archaeological remains. The results are in good agreement, especially in the upper layer where the ERTs, acquired with the electrodes fixed into the floor, have the better resolution and all the structures are distinguishable. As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, unequivocally, the importance of the comparison and integration of different methodologies for improving the localization and identification of the buried structures in humid and wetland scenarios.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 o As summarized in Table 3, GPR has permitted to quantitatively characterize, w high accuracy, the sizes of the buried objects when the strong geophysical contrast occ as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identifi tion of the structures more difficult. Regarding ERT, the higher water content of WL1 not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capab ity of ERT and GPR for archaeological purposes in humid and lacustrine scenarios ha provided interesting results. In particular, GPR measurements constituted a valid supp in humid/lacustrine scenarios, but their usefulness depends strongly on the size and p sitions of the buried structures as demonstrated in WP1. Furthermore, as shown in W GPR can work very well, both near the banks of rivers and lakes and above the wa providing an excellent resolution; however, a crucial role is played in this case by the quired investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and qu ity of the electrode contacts should be considered carefully in order to obtain results ceptable from the archaeological point of view. The best results in this work were achiev with analyses carried out with the 3D loop-loop array; whilst the use of the classical array with electrodes equispaced in two directions did not allow to easily localize the chaeological features. This result is very interesting for the applicability of the resistiv methods often not used in urban scenarios where regular disposition of the electrodes not always fully achievable for the constant presence of obstacle (trees, monuments, urb furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the me odology to identify strong anomalies in correspondence of the shallower structures sho ing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to ov come; therefore, a qualitative integration of GPR and ERT data was tried by co-render the images. In this research, the integration was effectively created by adopting 3D a 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, u quivocally, the importance of the comparison and integration of different methodolog for improving the localization and identification of the buried structures in humid a wetland scenarios.
Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, GPR has permitted high accuracy, the sizes of the buried objects when th as in the case of WL1, while the lower attenuation of tion of the structures more difficult. Regarding ERT not allow to reconstruct the real shape of the buried results are encouraging for real application in the ar

Conclusions
The geophysical tests realized at the Hydrogeos ity of ERT and GPR for archaeological purposes in provided interesting results. In particular, GPR meas in humid/lacustrine scenarios, but their usefulness d sitions of the buried structures as demonstrated in W GPR can work very well, both near the banks of ri providing an excellent resolution; however, a crucia quired investigation depth.
The applicability of ERT needs more attention; i ity of the electrode contacts should be considered c ceptable from the archaeological point of view. The b with analyses carried out with the 3D loop-loop arr array with electrodes equispaced in two directions d chaeological features. This result is very interesting methods often not used in urban scenarios where reg not always fully achievable for the constant presence furniture, etc.). In WP2, the 2D ERT profiles have hi odology to identify strong anomalies in corresponde ing the usefulness of the resistivity methods in very Despite great efforts carried out for the test, s come; therefore, a qualitative integration of GPR and the images. In this research, the integration was eff 2D data for the WP1 and WP2, respectively, and th quivocally, the importance of the comparison and in for improving the localization and identification of wetland scenarios.

Paved Road (C)
Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 of 27 As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, unequivocally, the importance of the comparison and integration of different methodologies for improving the localization and identification of the buried structures in humid and Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 o As summarized in Table 3, GPR has permitted to quantitatively characterize, w high accuracy, the sizes of the buried objects when the strong geophysical contrast occ as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identifi tion of the structures more difficult. Regarding ERT, the higher water content of WL1 not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capab ity of ERT and GPR for archaeological purposes in humid and lacustrine scenarios ha provided interesting results. In particular, GPR measurements constituted a valid supp in humid/lacustrine scenarios, but their usefulness depends strongly on the size and p sitions of the buried structures as demonstrated in WP1. Furthermore, as shown in W GPR can work very well, both near the banks of rivers and lakes and above the wa providing an excellent resolution; however, a crucial role is played in this case by the quired investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and qu ity of the electrode contacts should be considered carefully in order to obtain results ceptable from the archaeological point of view. The best results in this work were achiev with analyses carried out with the 3D loop-loop array; whilst the use of the classical array with electrodes equispaced in two directions did not allow to easily localize the chaeological features. This result is very interesting for the applicability of the resistiv methods often not used in urban scenarios where regular disposition of the electrodes not always fully achievable for the constant presence of obstacle (trees, monuments, urb furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the me odology to identify strong anomalies in correspondence of the shallower structures sho ing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to ov come; therefore, a qualitative integration of GPR and ERT data was tried by co-render the images. In this research, the integration was effectively created by adopting 3D a 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, u quivocally, the importance of the comparison and integration of different methodolog for improving the localization and identification of the buried structures in humid a

Stone Wall (D) Enchytrismos (E)
Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 of 27 As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, unequivocally, the importance of the comparison and integration of different methodologies Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 o As summarized in Table 3, GPR has permitted to quantitatively characterize, w high accuracy, the sizes of the buried objects when the strong geophysical contrast occ as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identifi tion of the structures more difficult. Regarding ERT, the higher water content of WL1 not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capab ity of ERT and GPR for archaeological purposes in humid and lacustrine scenarios ha provided interesting results. In particular, GPR measurements constituted a valid supp in humid/lacustrine scenarios, but their usefulness depends strongly on the size and p sitions of the buried structures as demonstrated in WP1. Furthermore, as shown in W GPR can work very well, both near the banks of rivers and lakes and above the wa providing an excellent resolution; however, a crucial role is played in this case by the quired investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and qu ity of the electrode contacts should be considered carefully in order to obtain results ceptable from the archaeological point of view. The best results in this work were achiev with analyses carried out with the 3D loop-loop array; whilst the use of the classical array with electrodes equispaced in two directions did not allow to easily localize the chaeological features. This result is very interesting for the applicability of the resistiv methods often not used in urban scenarios where regular disposition of the electrodes not always fully achievable for the constant presence of obstacle (trees, monuments, urb furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the me odology to identify strong anomalies in correspondence of the shallower structures sho ing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to ov come; therefore, a qualitative integration of GPR and ERT data was tried by co-render the images. In this research, the integration was effectively created by adopting 3D a 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, u quivocally, the importance of the comparison and integration of different methodolog Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, GPR has permitted to quantitatively characte high accuracy, the sizes of the buried objects when the strong geophysical cont as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the tion of the structures more difficult. Regarding ERT, the higher water content not allow to reconstruct the real shape of the buried objects, while in WL2 and results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing th ity of ERT and GPR for archaeological purposes in humid and lacustrine scen provided interesting results. In particular, GPR measurements constituted a val in humid/lacustrine scenarios, but their usefulness depends strongly on the si sitions of the buried structures as demonstrated in WP1. Furthermore, as show GPR can work very well, both near the banks of rivers and lakes and above providing an excellent resolution; however, a crucial role is played in this case quired investigation depth.
The applicability of ERT needs more attention; in particular, the resolution ity of the electrode contacts should be considered carefully in order to obtain ceptable from the archaeological point of view. The best results in this work wer with analyses carried out with the 3D loop-loop array; whilst the use of the c array with electrodes equispaced in two directions did not allow to easily loca chaeological features. This result is very interesting for the applicability of the methods often not used in urban scenarios where regular disposition of the elec not always fully achievable for the constant presence of obstacle (trees, monume furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of odology to identify strong anomalies in correspondence of the shallower structu ing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficu come; therefore, a qualitative integration of GPR and ERT data was tried by co the images. In this research, the integration was effectively created by adopti 2D data for the WP1 and WP2, respectively, and the results obtained demons quivocally, the importance of the comparison and integration of different meth Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, GPR has permitted to quantitatively high accuracy, the sizes of the buried objects when the strong geophys as in the case of WL1, while the lower attenuation of WL2 and WL3 m tion of the structures more difficult. Regarding ERT, the higher water not allow to reconstruct the real shape of the buried objects, while in results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for ity of ERT and GPR for archaeological purposes in humid and lacustr provided interesting results. In particular, GPR measurements constitu in humid/lacustrine scenarios, but their usefulness depends strongly o sitions of the buried structures as demonstrated in WP1. Furthermore, GPR can work very well, both near the banks of rivers and lakes an providing an excellent resolution; however, a crucial role is played in quired investigation depth.
The applicability of ERT needs more attention; in particular, the re ity of the electrode contacts should be considered carefully in order to ceptable from the archaeological point of view. The best results in this w with analyses carried out with the 3D loop-loop array; whilst the use array with electrodes equispaced in two directions did not allow to ea chaeological features. This result is very interesting for the applicabili methods often not used in urban scenarios where regular disposition o not always fully achievable for the constant presence of obstacle (trees, furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capa odology to identify strong anomalies in correspondence of the shallowe ing the usefulness of the resistivity methods in very conductive scenar Despite great efforts carried out for the test, some limitations ar come; therefore, a qualitative integration of GPR and ERT data was tri the images. In this research, the integration was effectively created by 2D data for the WP1 and WP2, respectively, and the results obtained quivocally, the importance of the comparison and integration of differ Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, G high accuracy, the sizes of the bur as in the case of WL1, while the lo tion of the structures more difficu not allow to reconstruct the real s results are encouraging for real ap

Conclusions
The geophysical tests realized ity of ERT and GPR for archaeolo provided interesting results. In pa in humid/lacustrine scenarios, bu sitions of the buried structures as GPR can work very well, both ne providing an excellent resolution; quired investigation depth.
The applicability of ERT need ity of the electrode contacts shoul ceptable from the archaeological p with analyses carried out with the array with electrodes equispaced chaeological features. This result methods often not used in urban s not always fully achievable for the furniture, etc.). In WP2, the 2D ER odology to identify strong anomal ing the usefulness of the resistivit Despite great efforts carried come; therefore, a qualitative inte the images. In this research, the i 2D data for the WP1 and WP2, re quivocally, the importance of the As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, une- Appl. Sci. 2022, 12, x FOR PEER REVIEW 23 of 27 As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, une-Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, GPR has permitted to quantitatively high accuracy, the sizes of the buried objects when the strong geophys as in the case of WL1, while the lower attenuation of WL2 and WL3 m tion of the structures more difficult. Regarding ERT, the higher water not allow to reconstruct the real shape of the buried objects, while in results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for ity of ERT and GPR for archaeological purposes in humid and lacustr provided interesting results. In particular, GPR measurements constitu in humid/lacustrine scenarios, but their usefulness depends strongly o sitions of the buried structures as demonstrated in WP1. Furthermore, GPR can work very well, both near the banks of rivers and lakes an providing an excellent resolution; however, a crucial role is played in quired investigation depth.
The applicability of ERT needs more attention; in particular, the re ity of the electrode contacts should be considered carefully in order to ceptable from the archaeological point of view. The best results in this w with analyses carried out with the 3D loop-loop array; whilst the use array with electrodes equispaced in two directions did not allow to ea chaeological features. This result is very interesting for the applicabili methods often not used in urban scenarios where regular disposition o not always fully achievable for the constant presence of obstacle (trees, furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capa odology to identify strong anomalies in correspondence of the shallowe ing the usefulness of the resistivity methods in very conductive scenar Despite great efforts carried out for the test, some limitations ar come; therefore, a qualitative integration of GPR and ERT data was tri the images. In this research, the integration was effectively created by 2D data for the WP1 and WP2, respectively, and the results obtained Appl. Sci. 2022, 12, x FOR PEER REVIEW As summarized in Table 3, G high accuracy, the sizes of the bur as in the case of WL1, while the lo tion of the structures more difficu not allow to reconstruct the real s results are encouraging for real ap

Conclusions
The geophysical tests realized ity of ERT and GPR for archaeolo provided interesting results. In pa in humid/lacustrine scenarios, bu sitions of the buried structures as GPR can work very well, both ne providing an excellent resolution; quired investigation depth.
The applicability of ERT need ity of the electrode contacts shoul ceptable from the archaeological p with analyses carried out with the array with electrodes equispaced chaeological features. This result methods often not used in urban s not always fully achievable for the furniture, etc.). In WP2, the 2D ER odology to identify strong anomal ing the usefulness of the resistivit Despite great efforts carried come; therefore, a qualitative inte the images. In this research, the i 2D data for the WP1 and WP2, re A further quantitative consideration can be realized comparing the geophysical detected anomalies both with the GPR and the ERTs with the real sizes and depths of the various objects buried in the subsoil. As previously shown, the used geophysical methods were able to detect the buried structures and the GPR was the best one in the different scenarios. In Figures 19 and 20, a preliminary quantitative analysis between the geophysical and expected results in the plan is shown in terms of position and size of the detected structures, which reveal the greater or lesser capability of the used methodology for the localization and characterization of the archaeological features for the different water levels adopted in WP1 where 3D acquisitions were performed. It is evident how in wetter conditions GPR is more able for detecting the real sizes of the elements. This is clearly due to the reflection focusing effects induced by the increase of the dielectric permittivity related to the background soil. GPR results further highlight a wrong positioning of the column if compared with the planned design. The best results are obtained for those elements with regular shape, as in the cases of the rectangular tomb, paved road, and stone wall. In detail, as showed in Figure 19, the anomaly of the column in WL2 is equal to the expected sizes of 1.4 m × 0.2 m; the reflections related to the rectangular tomb in WL1 highlight an anomaly really similar to the real one; similarly, the anomalies recorded in WL2 and WL1 for the stone wall with sizes of 1.6 m × 1.6 m. Further, there is a strong agreement between the reflections induced by the paved road and its planimetric development. Regarding the enchytrismos, only in WL1 were the detected sizes adequately matched with the real ones. The depth of the objects, when detected, generally complies with the real ones, and some discrepancies can be related to the heterogenous compaction degree of the sand placed above the structures. Moreover, GPR radargrams and ERTs acquired along S1 and S2 during WP2 (underwater conditions) are model co-rendered and the results are plotted in Figure 18. The first level of integration highlights the capability of both the methods to detect electric and em anomalies in correspondence of the buried structures despite the different resolution. Generally, the presence of strong reflections implies anomalies of the electrical behaviour of the subsoil related to the presence of archaeological remains. The results are in good agreement, especially in the upper layer where the ERTs, acquired with the electrodes fixed into the floor, have the better resolution and all the structures are distinguishable. Figure 18. Model co-rendered GPR and ERT data acquired along the two main sections (a,b) in WP2.
A further quantitative consideration can be realized comparing the geophysical detected anomalies both with the GPR and the ERTs with the real sizes and depths of the various objects buried in the subsoil. As previously shown, the used geophysical methods were able to detect the buried structures and the GPR was the best one in the different scenarios. In Figures 19 and 20, a preliminary quantitative analysis between the geophysical and expected results in the plan is shown in terms of position and size of the detected structures, which reveal the greater or lesser capability of the used methodology for the localization and characterization of the archaeological features for the different water levels adopted in WP1 where 3D acquisitions were performed. It is evident how in wetter conditions GPR is more able for detecting the real sizes of the elements. This is clearly due to the reflection focusing effects induced by the increase of the dielectric permittivity related to the background soil. GPR results further highlight a wrong positioning of the column if compared with the planned design. The best results are obtained for those elements with regular shape, as in the cases of the rectangular tomb, paved road, and stone wall. In detail, as showed in Figure 19, the anomaly of the column in WL2 is equal to the expected sizes of 1.4 m × 0.2 m; the reflections related to the rectangular tomb in WL1 highlight an anomaly really similar to the real one; similarly, the anomalies recorded in WL2 and WL1 for the stone wall with sizes of 1.6 m × 1.6 m. Further, there is a strong agreement between the reflections induced by the paved road and its planimetric development. Regarding the enchytrismos, only in WL1 were the detected sizes adequately matched with the real ones. The depth of the objects, when detected, generally complies with the real ones, and some discrepancies can be related to the heterogenous compaction degree of the sand placed above the structures. Despite the lower resolution, ERTs have detected some archaeological structures and in the case of the paved road and stone wall, have permitted to characterize the sizes with high detail. This is particularly true for WL2 when the water saturation degree is in intermediate conditions. Regarding the other buried objects, it is possible to define only the position; indeed, the sizes of the resulting electrical anomalies do not allow to reconstruct the expected geometries.     Despite the lower resolution, ERTs have detected some archaeological structures and in the case of the paved road and stone wall, have permitted to characterize the sizes with high detail. This is particularly true for WL2 when the water saturation degree is in intermediate conditions. Regarding the other buried objects, it is possible to define only the position; indeed, the sizes of the resulting electrical anomalies do not allow to reconstruct the expected geometries.
As summarized in Table 3, GPR has permitted to quantitatively characterize, with high accuracy, the sizes of the buried objects when the strong geophysical contrast occurs as in the case of WL1, while the lower attenuation of WL2 and WL3 makes the identification of the structures more difficult. Regarding ERT, the higher water content of WL1 do not allow to reconstruct the real shape of the buried objects, while in WL2 and WL3, the results are encouraging for real application in the archaeological field.

Conclusions
The geophysical tests realized at the Hydrogeosite laboratory for testing the capability of ERT and GPR for archaeological purposes in humid and lacustrine scenarios have provided interesting results. In particular, GPR measurements constituted a valid support in humid/lacustrine scenarios, but their usefulness depends strongly on the size and positions of the buried structures as demonstrated in WP1. Furthermore, as shown in WP2, GPR can work very well, both near the banks of rivers and lakes and above the water providing an excellent resolution; however, a crucial role is played in this case by the required investigation depth.
The applicability of ERT needs more attention; in particular, the resolution and quality of the electrode contacts should be considered carefully in order to obtain results acceptable from the archaeological point of view. The best results in this work were achieved with analyses carried out with the 3D loop-loop array; whilst the use of the classical 3D array with electrodes equispaced in two directions did not allow to easily localize the archaeological features. This result is very interesting for the applicability of the resistivity methods often not used in urban scenarios where regular disposition of the electrodes are not always fully achievable for the constant presence of obstacle (trees, monuments, urban furniture, etc.). In WP2, the 2D ERT profiles have highlighted the capability of the methodology to identify strong anomalies in correspondence of the shallower structures showing the usefulness of the resistivity methods in very conductive scenarios.
Despite great efforts carried out for the test, some limitations are difficult to overcome; therefore, a qualitative integration of GPR and ERT data was tried by co-rendering the images. In this research, the integration was effectively created by adopting 3D and 2D data for the WP1 and WP2, respectively, and the results obtained demonstrate, unequivocally, the importance of the comparison and integration of different methodologies for improving the localization and identification of the buried structures in humid and wetland scenarios.
In the future, greater efforts will be required to integrate the geophysical data aimed at obtaining quantitative information about the planimetric and volumetric structures placed in the subsoil also involving the use of advanced adopting algorithms (i.e., image data-fusion, machine learning algorithms). At the same time, a further topic of the GPR research analyses will be to identify the best strategies for the migration of the hyperbola in order to make quantitative analyses possible.