Stabilization Methodology in Foundation Soils by ERT-3D Application in Estepona, South Spain

Featured Application: A non-invasive solution developed to preserve damaged buildings using an injection of cement grout in the soils to stabilize it. Abstract: The paper proposes a novel methodology for the stabilization of shallow foundations, with a simplified model combined with 3D electrical resistivity tomography (ERT-3D and consolidation injections. To determine its usefulness, the method has been applied in a case located in Estepona (southern Spain). The chosen tomography model is the dipole–dipole configuration, with an optimized distance between electrodes of 0.80 m for a better visualization of the foundation subsoil; with this parameterization, a total of 72 electrodes were installed in the analyzed case. In this work, the depth of the anomaly in the building’s supporting subsoil was detected ranging from 2.00 m to 3.90 m deep. The study also delineates areas of high resistivity variations (50–1000 Ω m) in the middle and eastern end of the field. These data have been validated and corroborated with a field campaign. The results of the ERT-3D monitoring are presented, once the investment data has been processed with the RES3DINV software, from the beginning to the end of the stabilization intervention. The novelty occurs with the interaction between the tomography and the foundation consolidation injections, until the final stabilization. This is a very useful methodology in case of emergency consolidation, where there is a need to minimize damage to the building. Thus, people using this combined system will be able to practically solve the initial anomalies of the subsoil that caused the damages, in a non-invasive way, considerably lowering the value of the resistivities.

is common to find many stress-induced cracks and other defect-related issues due to the various structures and foundations of historic, public, and private buildings, and some of these problems are disasters in the making [12][13][14]. This necessitates proposing a non-invasive intervention methodology for the subsoil of these buildings to avoid collapse.
One of the most used non-invasive techniques for subsoil exploration is geophysics and, more specifically, electrical resistivity tomography, in both two and three dimensions [15][16][17][18][19][20][21], and metrological perspectives of tomography in civil engineering [22]. In recent years, this method has proven to be an efficient tool, not only to monitor wall degeneration [23], but also to detect other types of problems in building foundations, such as subsoil degradation [14,24,25]. Spatial variability in load-bearing rock/soil beneath a building puts stress on poorly supported structures, causing failure [10,26]. Therefore, a geophysical investigation offers a faster, non-invasive means of obtaining detailed, credible information about the subsurface under a building. Electric tomography can also image the ground's distribution and structural deformation, both of which offer credible information regarding the strength the rock/soil is likely to offer a building [24,27,28]. Electrical resistivity imaging has, for decades, been very effective in illuminating the subsurface, and apt at providing information about the soil's physical properties for economic, environmental, and geological engineering. As mentioned, this technique is a useful and non-invasive method to diagnose subsoil problems in shallow foundations of buildings [3,[29][30][31][32], which represent 25% of the claims reported in Spain, according to recent data provided by Aseguradora Mutua de Arquitectos Superiores, a public limited company (ASEMAS S.A.). However, the three-dimensional geophysical research approach ERT 3D [33][34][35] is better able than the two-dimensional approach ERT 2D to characterize the subsurface and determine heterogeneity in measured rock properties along the vertical (Z) and orthogonal horizontal (X, Y) axes [36,37].Thus, the determination of the variation in the soil properties along the three orthogonal directions would allow the evaluation of the spatial variation in the strength of the foundations, imposed by the heterogeneous properties of the soil [38].
Electrical resistivity tomography is an appropriate procedure for detecting and controlling underground consolidation and stabilization, specifically the ERT-3D tomography technique [39][40][41][42]. As mentioned above, the method consists of setting out parallel lines of observation, which cover the study area, to obtain underground data; however, the key point lies in processing the data. Fortunately, important advances have been made in 3D resistive imaging and its inversion processing through applying the powerful resistance inversion software, RES3DINV [43].This is possibly the most widespread, accurate, simple, and affordable software for data inversion in electrical tomography, and, thus, it has been used in this work. Several factors have been reported to impact the variation of electrical resistivity in the subsurface. These include variation in rock type, rock fabrics, rock deformation, water saturation, different degrees of weathering, etc. [44,45]. These factors, which are known to impact spatial variation on ground electrical resistivity, are also capable of impacting variation in other physical properties of rocks [46,47].With different electrode distances, ERT-3D offers the possibility of locating these gaps/cavities and possible holes under the foundation of a building, which can cause severe deformation and settlement. Both consolidating and stabilizing soil through injections have been proven effective [48][49][50][51] and are used to mitigate or even solve differential settling. The injection material generally depends on the lithological morphology of the soil in question. Synthetic resin is frequently used as an injection fluid [52], but this paper proposes injections of cement grout as a more economical and versatile solution.
For a soil consolidation project to be successful, electrical resistivity data must be available in advance to provide initial information on the subsurface structure. In our methodology for consolidating and stabilizing buildings with shallow foundations, we will use electrical tomography as a tool that detects empty cavities, which are typical examples of anthropic (resistive) fillings with low bearing capacity in subsoils [3]. Thus, this work proposes ERT -3D to monitor the subsoil in real time, from the beginning to the final consolidation through the different injection passes into the subsoil, until its stabilization. As an example of a case of application of our combined methodology, we have chosen a historic building with serious stability problems in its foundations [5], located in Estepona, on the Spanish south coast. The subsoil was consolidated with injections of cement grout applying the electrical methodology, resulting in the preservation of the building, and avoiding its possible collapse.

The Model Development
The ERT method discharges an electric current into the ground and measures the potential difference at two determined points on the surface. The suitability of this method lies in the fact that irregularities in the subsoil beneath a building can be identified as contrasts or anomalies in the subsoil's electrical properties.
This method is based on Ohm's law: where: Ρa = apparent resistivity [36]; k = geometric constant that depends only on the reciprocal positions of the current and potential electrodes; ΔV = potential difference; l = intensity of the injected current.
The apparent values of resistivity depend on the real resistivity distribution in the tested area. The true resistivity distribution can be estimated through a reverse procedure, based on minimizing an adequate function [15,39,44]. The solution to this problem is not unique. For the same set of data, a wide range of models can calculate the same apparent resistivity values. A preliminary lithological analysis is usually carried out on the subsoil's nature, to reduce the range of possible models, which can be incorporated into the reverse subroutine.
The solution method used minimizes the difference between the apparent resistivities measured and those calculated by the RES3DINV software, which uses the limited softness inversion formulation, restricting the model's change in resistivity values [15,36,43,53]. This study has used RES3DINVx64, which implements a smooth routine based on least squares and is practically the only commercial software of its type available [36,54,55]. The inversion routine used by the program is based on the smoothness-constrained least-squares method [56]. The basic smoothness-constrained least-squares method is based on the following equation.
(J T J + λF) Δqk=J T g−λ Fqk-1 (2) where: J T = Transpose of J. J = Jacobian matrix of partial derivatives. λ = Damping factor. q = Disturbance vector. k = Iteration number g = Data mismatch vector. αx, αy, αz = weights for roughness filters Cx, Cy = horizontal roughness filters Cz = vertical roughness filter This method is advantageous because of its versatility since the damping and roughness factor filters adjust to the different types of data. This program uses the Gauss-Newton method, which recalculates the Jacobian matrix, after each interaction [57]. The interpretations of electrical tomography profiles are made using RES3DINVx64 [34,36] for resistivity and induced polarization. As mentioned above, this calculation software is based on the least-squares method with forced smoothing, modified with the Quasi-Newton optimization technique. The inversion method designs a subsoil model using rectangular prisms and determines the resistivity values for each of them, minimizing the difference between the apparent resistivity values observed [15,58]. The results of our model show an uncertainty in its acceptable magnitude. They are in a critical state with respect to their structural stability, load capacity and, under these conditions, the safety of the building's habitability. This is at an acceptable level of uncertainty associated with these increasing risks of instability [5].

Application
Some results reported in the literature and obtained in the last eleven years regard the application of the ERT-3D method to study geotechnical anomalies in the subsoils on slopes after landslides located in different geographical contexts.Such reports make us consider the ERT method as a tool and methodology very suitable to investigate these geotechnical anomalies of the subsoil during the pre-event and post-event phases of a disaster cycle [30,31,[59][60][61] or simply a poor foundation support subsoil (anthropic fill).
In fact, during the pre-event phase, the resistivity contrasts that characterize ERT-3D allow defining the geological environment of the subsoil. They allow identifying areas of high-water content that could be responsible for reactivation events. In the post-event phase, ERT-3Denablesus to reconstruct the damaged or altered subsoil body by also providing information on the volume of the removed or altered material. This information can help to better plan future mitigation activities. Our original application not only detects these subsoil alterations in real time, but we also propose an original tool for subsequent mitigation, in case of affecting existing buildings or infrastructure in a non-invasive way and thus preventing their collapse.
One of the biggest drawbacks of 3D tomography for the investigation of shallow and disturbances in the subsurface was the fact that it did not provide continuous acquisitions over time, which made it unsuitable for the study of the dynamic nature of shallow and disturbances.
Fortunately, the development of systems for the continuous acquisition of electrical resistivity time and software for data inversion [61] are paving the way to test this method during the emergency phase, such as the application and methodology developed in the present work, where in real time we are applying the geoelectric investment data and carrying out the consolidation of the subsurface area, increasing its resistivity in the emergency case, and thus being able to recover existing constructions in the altered subsoil.
The possibility of using ERT-3D to monitor geotechnical changes and alterations in the first layers of a ground clearance-settlement area will add important information during the emergency phase. At this time, the preliminary results obtained in our methodology when applying the ERT-3D for this purpose are very encouraging. The proof of this is that we were able to recover the building that was based on the geotechnically altered substrate, where the grout cement, filler material proposed in our methodology for consolidating the subsoil, satisfactorily filled the gaps that water and air occupied in the altered substrate. Cement grout has advantages over other fillers, such as its economy, easy means of performing on site, ease of dosing, and versatility, among others.
In applying our methodology in these cases of emergencies where there is a removal of surface soil after a shallow of the subsurface that affects buildings and infrastructures [62,63], we have observed that a low resistivity zone lies in the upper area from the surface to a depth up to 4.00 m, and the zone shows a higher resistivity of greater than 400 Ω m. According to the data obtained in this work, the upper area of the soil unit consists of mostly silt soil with granule from the erosion of the bedrock. Therefore, we thought that the anomaly in the upper area is attributed to the silt soil and granule with low resistiv-ity. This zone of low resistivity coincides with the zone of removal and presence of water, which we consolidated with our grout cement filling methodology and in real time the monitoring was carried out to verify the satisfactory level of the landfill.
The procedure followed by the authors clearly defines the altered surface material with cavities and gaps of the rocky matrix (phyllite) by lower resistivities in the bedrock (50-150 Ω m). This methodology characterizes the subsoil altered by shallow bedrock material mainly composed of (by) clay material (colluvial) with high pore-hole content and high resistivities. The modelling of the subsoil analyzed with our system was monitored with the ERT of the grout process using the 3D RESINV software [63]. Another novelty of our tool is the geological characterization of this type of shallow phenomenon, and its affectation to existing buildings on the slope that in the South of Spain, in La Cordillera Bética, occurs with recurrence [62][63][64]. To have a tool like this proposal to characterize and mitigate damage to existing buildings and infrastructures in a non-invasive way in this type of geological formation is an advantage.
The results after the application of the methodology in the case study showed the effectiveness of the diagnostic and intervention methodology for mitigating the serious damage suffered by the building, preventing the collapse and destruction of the building while preserving the safety of its inhabitants.

Phases: Developed Methodology
Anew methodology has been developed consisting of different phases and applied to the case study: 1. It is necessary to interpret the characteristics of the subsurface soil/rock in an electrical resistivity tomography methodology. The ERT-3D tomograms will show how subsurface materials affect the resistivity of the subsurface. 2. Using a geophysical method to boreholes will provide a detailed knowledge of the studied site, with no limits over analysis and misinterpretation. 3. This will delineate depths and thicknesses of subsurface layers and identify the stiffness and distribution of subsurface materials in the Baetic Mountain Range. 4. The electrical resistivity of metamorphic rocks in the Baetic Mountain Range (phyllite) is dependent on the degree of fracturing and the resistivity of the interstitial and pore water in the rock and its water content in such facture zones. Thus, the electrical resistivity ranges from 3000 of Ω m to a fraction of 1.00 Ω m. 5. Finally, a non-invasive system of filling of holes through the injection of cement grout in the subsoil will be applied, monitored in real time with the ERT and thus increase the resistivity of the subsoil, and with it, its consolidation.

Case Study
This research shows the results of applying the proposed methodology to a specific case of differential settlements in a building in anat-risk area after soil removal [5]-see Figure 1. The land movements occurred after heavy rains [64] in the 2009-2010 hydrological years, but measurements presented in this paper were taken in 2012.
. The applied methodology is based on ERT-3D and consists of placing electrodes along profiles separated from each other according to the resolution, depth, and objectives to be covered. The lower the separation, the greater the resolution; the greater the separation, the greater the depth. An optimal distance of 0.80 m. has been determined in our methodology between electrodes to obtain a balance between resolution tomographic profiles and depth in accordance with the problem in question.
Prior to performing the tomographic profiles, it is recommended that preliminary geological research be carried out in the area where the buildings are located [64]. In this case, the affected buildings were located at Paraje del Arroyo, La Cala, Estepona, Malaga (Spain), which is in the south of the Baetic Mountain Range (South Iberian Peninsula). The damage occurred at coordinates 36.461094, −5.160498 ( Figure 1). In the first phase, it is proposed that a granulometric analysis and a Standard Penetration Test (SPT) be executed in the geotechnical surveys. In the present case, soil consistency increased with depth (Table 1). Change in the geotechnical response of the soil occurred at an approximate depth between one and four meters, according to the SPT hits along the analyzed sample (Table 1). Note that the numbers in bold correspond to values that are too low. The grey background switching to white marks a discontinuity and an important change of capacity and bearing resistance in the subsoil, coincident with level I of Table 2.
Based on geological and geotechnical studies, the materials extracted in the area depict the lithological levels listed in Table 2. The analyzed soils were mainly rocky and clayey, of varying thickness, which could be verified in field work and through laboratory geotechnical tests.

Damaged Building Analysis
The most common deformations in a building are related to foundation differential movements. These alterations result in structural deformations in the building and angular distortions in its foundation [3,65]. This provokes stress throughout the construction, and when this stress limit is exceeded, cracking or breakages occur.
The Spanish Building Technical Code (2006) defines the maximum angular distortion values for a building's ultimate limit of state and service, conditioning this value to the type of structure. In the service limit state, the allowable angular distortion value is L/500 for reticulated structures, and L/300 for isostatic structures and for load-bearing walls (the studied case). For L, the length in a straight line between the axes of the footings is analyzed. However, the structure's ability to assume these deformations will depend, among other factors, on the stiffness of each element comprising it, so any fixed distortion value could be conservative or give rise to non-tolerable deformations [66][67][68][69].
According to current urban regulations, the rural house examined in this study would be difficult to demolish and rebuild. Therefore, a major challenge to this research and methodology was finding a solution to recover the building without having to demolish it. The main building in this study was a rectangular two-story building, 17.  This was a residential building and was used as a rural house (country holiday house); hence, guest and owner safety was paramount. The distribution comprised two building components: two floors, with an open area on the ground, and an upper, noble floor ( Figure 2). Both the main building and the entrance foyer were designed in the late 1980s and completed in the early 2000s.The building's structure consisted of load-bearing walls and solid brick pilasters on the façade's main porch (vertical structure); these rested on a shallow foundation base of a reinforced concrete simple wall footing, formed by three parallel simple wall footings (A, B, C) ( Figure 2). The horizontal structures of the building were24 cm thick unidirectional slabs formed of semi-resistant joists, supported in three spaces by load-bearing walls and pilasters resting on a shallow foundation composed of three longitudinal foundations in a trench.
From the previous geotechnical and structural analysis of the building illustrated in Figures 2 and 3, a problem of support in the foundation for differential seats was deduced. Of the A, B, and C foundations, B and C had become unstable (sunk downward) with respect to A (Figure4). This phenomenon caused the building to swing forward, which caused serious damage to its structure, as can be seen in Figure 4. The foundation of wall A also suffered a small movement in its seat (S1) (Figure 4). . Cross-section scheme, presenting the likely behavior of the building's foundation and structure. There was a differential settlement of the S1 and S2 foundations due to a possible settlement-removal of the soil and the lack of bearing capacity of the resistant subsoil. Displacement occurred through the variable Φ = 95 mm, which is the depth B and C have descended, with respect to A. The seat of the foundation for wall A (S1) has lowered a depth of 25 mm, remaining within the admissible limit according to the Spanish Building Technical Code (CTE).
Due to the type of movement and structural typology, displacement occurred as shown in Figure 4. When load-bearing walls deform, they collapse, and destruction of the building becomes a major threat.
To determine the building deformation as a function of the angular distortion, it was calculated as the differential settlement, defined as the settlement difference between foundations A and C, which were the extreme foundations affected by the settlements, applying the following equation. Knowing that in this case, it was the same height that has descended B and C with respect to A: where δAC = foundation vertical differential displacement between points A and C; SA = suffered differential settlement at point A; SC = suffered differential settlement at point C.

The equation of angular distortion is also applied [3,70-72]:
= where:β = angular distortion; LAC = distance between foundation points A and C. The Spanish Building Technical Code (2006), also known as CTE, establishes limitations on β movements, according to building structure typology. Studies regarding β began in the 1940s [73] and continue today [74][75][76]. The CTE limitation for the β value, for building damages, must be >L/300 for load-bearing walls. In this case, the differential settlement between foundations A and C was δAC = 95 mm and β = 0.0086, which are above the admissible limit established by the CTE, indicating imminent building collapse.
Authors also link angular distortion β values with damages suffered in enclosures (shell and internal partitions) and structural elements. The present case exceeded the limit values established by these authors, where the dimensions of cracks in enclosures must be >25 mm to be considered a serious structural risk, according to a programmed damage classification [77]. As shown in Figure 4, this building had serious damage, due to cracks in enclosures and structural elements (load-bearing walls), which were measured at greater than 50 mm. Therefore, according to this classification, the building was close to collapse. Due to the type of structure, foundation, and damages suffered, the methodology proposed in this research seems appropriate for the case. It would allow detection of zones of soil weakness where the differential settlement has been produced and consolidation of the subsoil from outside the construction, thus solving the building's issues non-invasively.

Analysis and Stabilisation Methodology
According to the data-regulations, the state of the building, the undeveloped land surrounding the building, and the necessity of keeping the building stable (i.e., avoiding imminent collapse)-a simplified model was proposed, in which ERT-3D would be used, with subsequent processing of the reversal data via the RES2DINV software [43], as seen in Figure 5. This analysis will consider the geotechnical data given in the previous sections. The developed methodology allowed the researchers to see the reality of the soil, as well as any possible cavities in the subsoil supporting the building's foundation. These cavities were filled with grout to consolidate the soil, and then, via contrast campaigns (ERT-3D), the researchers checked to determine if the filling worked, if it moved, etc. Thus, the proposed methodology provided a complete, three-dimensional analysis of the subsoil supporting the damaged foundation, where the differential settlements occurred. Later, along with this original geophysical technique, the researchers monitored, through successive contrast campaigns, the consolidation and underpinning of the subsoil with different passes of consolidation injections. The injection material depends on the lithological morphology of the soil to be consolidated. Due to the type of subsoil in this case, as well as the material's low cost and ease of use and control, this paper proposed injections of controlled cement grout [78,79]. Synthetic resin is frequently used as an injection fluid, though it does not allow for high solicitation [52]. Three beneficial effects of injecting cement grout into the subsoil have been observed: (1) filling existing gaps, (2) soil compaction, and (3) interstitial water reduction and/or elimination.

Approach Methodology: Application
Had the proposed methodology not been applied in this case, the building could have been damaged to a critical point, making its repairs too expensive to pursue. Other, more invasive proceedings, such as underpinning through micropyles [5,[80][81][82], were not feasible due to the state of the building and the structural typology of the load-bearing walls. However, making consolidation injections into the subsoil allowed recovery, to a certain degree, of the differential settlement suffered and rendered the building useable again.
The first necessary action involved using electrical tomography to analyze the subsoil beneath the affected foundations. In an initial campaign, this study accomplished various tests with spacings of 0.80, 1.00, and 1.50 m to identify the best resolution for checking the subsoil's state. A spacing of 0.80 m was chosen, based on another recent research [11].
The Syscal Switch 48 (V114 ++) from IRIS Instruments was used as electrical tomography equipment. It is a multi-electrode piece of equipment with an integrated computer capable of managing up to 900 electrodes, with resolution characteristics of resolution/accuracy:1 µV/0.2% [43]. The equipment power source is 250 W and 2.5 A, which generates 880 Vp-p pulses, and the manufacturer has incorporated a transmitter and receiver into the system. Among the equipment's features are a time injection adjustment, an apparent resistivity and chargeability automatic processor, a 3D real-time resistivity control, a voltage and current injection curve control, an integrated PC, and a commutation processor. The locations of the two topographies' parallel profiles (E1-E2) were shown via coordinates given by a GPS from Garmin Etrex Ventura [14], with a WG84 coordinate system. The results are listed in Table 3. Table 3. UTM coordinates of the tomography profiles; location of the profiles according to the distribution of Figure 6.

Start UTM Coordinates of Tomography
Profile E1. In the developed methodology, two parallel profiles (E1 and E2) were selected, covering an area of approximately 240 m 2 , which included the main affected building area in accordance with Figure 6 and Table 3. The electrodes were placed in parallel. ERT -3D (TE1-3D) delimited tomography surface, covering the main building's foundation structure, within the parallel profiles proposed: E1 and E2. Circles refer to the E1 profile, while squares refer to E2; orange refers to the initial UTM coordinates, while green refers to the final ones (Table 3).
RES3DINV-3D data inversion software for electrical imaging and induced polarization (IP) [83]-was used to automatically invert the acquired data from apparent resistivity and produce a 3D resistivity model. The two parallel profiles were executed using a dipole-dipole matrix. The dipole-dipole matrix was chosen due to its sensitivity to lateral variations of resistivity [41,84]; thus, structures such as gaps and cavities that were damaging the building's subsoil could be detected. This dipole-dipole matrix provided us with a sharper horizontal resolution of the whole, which is an important advantage in this type of building pathology. As discussed previously, an inter-electrode separation of 0.80 m was estimated as suitable for this methodology. The two profiles extended along 28 m, with 72 electrodes distributed according to the coordinates listed in Table 3 and Figure 6. Figure 7 shows a 3D electrical resistivity model from the studied terrain and the depths reached, in this case, up to 5.80 m depth. Both profiles show resistivity values represented in a range of colors for better and easier observation of the variations in subsoil vertical and horizontal resistivities. Zones with resistivity values between 0-55 Ω m corresponded to levels of clay or colluvial material, represented in green; zones between 55-1000 Ω m, shown in orange, corresponded to fillings or colluvium removed with low compactness; and resistivity >1000 Ω m, in red, indicated anomalies due to the presence of interstitial gaps. Figure 7. Advanced inversion algorithm of 3D synthetic data (RES3DINV software). Electric Ground Section ERT-3D: cross-section diagram of the ground, with the probable behavior of the soil and foundation. Cross-section scheme, with the likely behavior of the soil and foundation. Figure 7 shows resistivity anomalies presented in slices of subsoil supporting the foundation of the damaged building. The red color shows a higher resistivity value, up to 3000Ωm, which suggests that there are cavities causing the resistivity to rise. From the obtained resistivities, combined with the geotechnical data (SPT standards) [85,86], the physicochemical properties of geological materials can be determined.
The survey depth was established as 5.80 m, considering the previous geotechnical results of Table 1. From that depth, level III began to appear, determined by a rocky matrix of phyllite. The main anomalies were located on the main load-bearing walls, specifically A and B.
The structure of the building had been seriously damaged by the seat of its foundation and the turning of the building, specifically in the main facade between pilasters and the main bearing wall (wall B). There were also apparent interstitial anomalies in the subsoil up to 3.90 m deep in the dielectric profile E2 (wall B), and 2.00 m in the dielectric profile E1 (wall A). The anomalies and settlement in the foundation of wall B (S2) were of greater importance than the settlement at the base of wall A (S1), with a characteristic angular distortion of β = 0.0086, hence the forward inclination of the building. In sum-mary, according to the algorithm developed in Figure 6, the first phase of the tomographical campaign was carried out, by analyzing the resistivities of the subsoil supporting the building and considering previous data, including the initial geotechnical data. Anomalies and areas with low compaction and/or gaps were identified and isolated in ERT-3D for the building's main load-bearing walls (Figure 8), so this knowledge might be acted on in the future. Given the powerful information about the subsoil under the building's main load-bearing walls, the second phase of this methodology was proposed: consolidation injections into the subsoil for those detected anomalies (flow diagram, Figure 5). Figure 9 isolates the areas to be treated with injections via the RES3DINV investment program. The damages observed in the analyzed building and the initial geotechnical tests confirmed the diagnosis from the ERT-3D tomography methodology. The initial ERT-3D results-that is, the application of the first phase of the flow chart ( Figure 5)-showed high coincidence in the geological composition of the supporting subsoil beneath the damaged building: altered clay edges and colluvial sand from the meteorisation of the rock matrix (phyllite).
This innovative intervention technique allowed the researchers to evaluate the induced effects in the field and, in light of them, carry out possible modifications in the distribution of injection points to consolidate the subsoil. In the procedure's second phase, after the injections in the anomalies were made, the ERT-3D contrast was repeated to check the consolidation results.
The second phase of the proposed procedure ( Figure 5) intended to obtain uniformity of chemical and physical features in the stabilized subsoil [65,70,86] until the desired goal was reached. As shown in Figure 10, the contrast tomography gave a favorable result after the cement grout injections were applied. The resistivities up to 3000 Ω m that caused distortions (cavities), affecting the stability of the building in Figure 8, were generally reduced drastically, even at values below 55 Ω m. This means the anomalies and stability issues of the building were resolved. In this way, the affected building could be recovered, and demolition avoided. The proposed intervention methodology is, therefore, highly effective in buildings damaged by settlements in shallow foundations, due to a deficient stress capacity of the subsoil, in this case by a removed anthropic fill. In this case study, a building on the verge of collapse, was recovered. Figure 10, showing the tomography performed after consolidation injections, indicates that up to 3000 Ω m resistivity anomaly areas have disappeared and left the building's subsoil consolidated and stabilized.
Four beneficial effects of our combined interactive methodology (tomography-injection) have been observed on subsoil where a building foundation rests: (1) Fills existing gaps improving its bearing capacity.
(3) Reduces and/or eliminates interstitial water in subsoil gaps. (4) Recovery of severely damaged buildings with hardly any collateral effects.

Results and Discussion
The interpretation and adaptation of this paper's proposed methodology were satisfactory and successful. The building, on the verge of collapse, recovered its stability and was ready for use again. This procedure is promising for buildings that have reached their ultimate limit state but must be maintained either because they could not be rebuilt due to urban regulations (as in this case), or because they are listed as historic and cultural heritage buildings to be preserved at all costs.
The advantage of this proposed method is that ERT-3D allows researchers to determine if the foundation settlement and subsidence of the subsoil has been permanently eliminated after consolidation injections. Consolidating foundations through grout injections is inherently problematic because researchers do not know how much injection grout or how many injections are necessary, and they are unable to control or physically see the result, as it is underground. Therefore, an intervention methodology, such as the one presented here uniting ERT-3D monitoring with the consolidation injection process, is essential for checking the results. The advantage of this ERT-3D methodology, as well as consolidation grout injection compared to other reclining systems (e.g., micropyles), is that it is a non-invasive method for an already damaged building. This makes the system useful in extreme cases, like the one analyzed here.
In the application of our methodology, two electrical profiles were made, E1 and E2 ( Figure 6), an upper part of approximately 3.90 m thick was defined for E2 (wall B) and 2.00 m for E1 (wall A) with high resistivity, ranging from 500 to up to 3000 Ω m. This section corresponded to the colluvial material (porous) that covers the metamorphic rock (phyllites). Under this lithology was the rocky matrix, formed by altered phyllites with less resistivity, of the order of 55 Ω m. The upper colluvial material, which was in contact with the underside of the building's footings, showed a significant degree of alteration, which means a high resistivity (from 1000 to up to 3000 Ω m). However, the resistivity values were quite variable within this mass of colluvial material, and this is linked to two causes:(1) lack of compaction and (2) presence of gaps and cavities in the subsoil where the foundation of the building supported, which coincided with the levels I and II altered. These data confirm the instability results of the shallow foundation of the building. The building's final state, after the application of this methodology, is shown in Figure  11. Figure 11. Current state of the main building. Image taken after the intervention and application of the methodology. Had this methodology not been used, the only solution would have been demolition of the building.

Conclusions and Follow-Up
The ERT-3D technique proposed in this research proved to be versatile, fast, and cost-effective for detecting subsurface anomalies. In the present case study, the researchers were able to analyze the subsurface to a depth of 5.80 m, deep below the severely damaged building. Nevertheless, anomalies only existed at a depth of 3.90 m, below the foundation of load-bearing wall B. When the subsoil has holes, those spaces are filled with air or with water which, when evaporated, are dielectric. This means that the terrain presents a strong gradient anomaly and very high resistivity values, which causes serious damage to the foundation and structure of the building. Our methodology solves this problem, by reducing the resistivity of these initial anomalies (cavities) by means of the combined grout injection in two passes, with the control of resistivities through 3D tomography, and returning the stabilization of the foundation of the building in a non-invasive way. It has been applied to a case study that has been totally effective, where initially we had resistivities up to 3000 Ω m and it was effectively reduced with our methodology to 55 Ω m.
In the present case study, to implement the methodology, the subsoil cavities and weaknesses were first identified and located beneath the building's foundations. Secondly, a consolidation grout injection campaign was carried out, enhancing the subsoil conditions to an optimum result. The consolidation fluid injections eliminated the high value of initial subsoil resistivities (up to 3000 Ω m) and moved those resistivities to values between 55 and 1000 Ω m, filling in those cavities. Through this process, ERT-3D showed itself to be a fundamental, minimally invasive tool for research. It was useful both during the project phase and in the follow-up, underpinning work when the consolidation injections were placed under the foundations.
The limitations of our methodology increase with the depth of the altered subsoil, from 5.80 m. As has been shown in our research, it is difficult to have a reliable electrical section of resistivities to intervene later with consolidation injections. Although it would be necessary to use another geophysical methodology and another type of deep consolidation intervention using a pile or micropyle, our methodology is appropriate for emergency surface consolidations in disturbed porous subsoils. As they work safely outside the affected building, this methodology offers researchers an accurate lithological model and correctly highlights subsurface anomalies. Hence, ERT-3D allows them to determine the causes of the instability of the foundation and building and how it can be later consolidated and stabilized, in addition to monitoring the chosen solution.