Studying a Subsiding Urbanized Area from a Multidisciplinary Perspective: The Inner Sector of the Sarno Plain (Southern Apennines, Italy)

Defining the origin of ground deformation, which can be a very challenging task, may be approached through several investigative techniques. Ground deformation can originate in response to both natural (e.g., tectonics) and anthropic (e.g., groundwater pumping) contributions. These may either act simultaneously or be somewhat correlated in space and time. For example, the location of structurally controlled basins may be the locus of enhanced human-induced subsidence. In this paper, we investigate the natural and anthropic contributions to ground deformation in the urbanized area of the inner Sarno plain, in the Southern Apennines. We used a multidisciplinary approach based on the collection and analysis of a combination of geomorphological, stratigraphical, structural, hydrogeological, GPS, and DInSAR datasets. Geomorphological, stratigraphical, and structural data suggested the occurrence of a graben-like depocenter, the Sarno basin, bounded by faults with evidence of activity in the last 39 ka. Geodetic data indicated that the Sarno basin also experienced ground deformation (mostly subsidence) in the last 30 years, with a possible anthropogenic contribution due to groundwater pumping. Hydrogeological data suggested that a significant portion of the subsidence detected by geodetic data can be ascribed to groundwater pumping from the alluvial plain aquifer, rather than to a re-activation of faults in the last 30 years. Our interpretation suggested that a positive feedback exists between fault activity and the location of area affected by human-induced subsidence. In fact, fault activity caused the accumulation of poorly consolidated deposits within the Sarno basin, which enhanced groundwater-induced subsidence. The multidisciplinary approach used here was proven to be successful within the study area and could therefore be an effective tool for investigating ground deformation in other urbanized areas worldwide.


Introduction
Ground deformation is a complex phenomenon that may derive from both natural (e.g., tectonics, landslides, and sinkholes) and anthropic (e.g., human-induced subsidence) causes. Natural-and anthropic-induced subsidence may act simultaneously [1] and, in some cases, areas affected by human-induced subsidence have been shown to be clustered within basins bounded by active faults [2,3]. This makes the analysis of the ground deformation a very challenging task, especially in urbanized areas. Here, tectonic contributions to ground deformation may be disclosed by geomorphological, stratigraphical, and structural analyses [4], whereas anthropic contribution may be detected by geodetic analysis [5].
Further difficulties in the recognition of the causes that drive ground deformation may derive from the location of the investigated urban area. In fact, natural (e.g., tectonic) contributions to ground deformation are very elusive in low seismicity areas. Difficulties  [20,34]). White box indicates location of the study area, shown in inset map B. Active faults are derived from the ITHACA database [35] and from Cinque et al. [36]. Earthquake epicentres and magnitude are derived from the Italian Earthquake Catalogue CPTI15 [37]. (B) geological map of the inner sector of Sarno plain (modified from [20]). This stratigraphic setting hosts a multi-layered alluvial aquifer, which is recharged by local effective precipitation and, subordinately, by lateral groundwater flow coming from the Sarno Mountains karst aquifer. The groundwater circulation is directed towards the Sarno River and its tributaries [38−40]. A high number of uncontrolled and unauthorized private wells (up to approximately 90-100 wells/km 2 ) use the plain groundwater for industrial, domestic, and agricultural purposes [5]. Other groundwater withdrawals are carried out at three well-fields and two gravity drainage galleries located near the Sarno Mountains karst aquifer.
The inner sector of the Sarno plain has been deeply investigated by means of DInSAR analysis [1,21], which suggested intense human-induced subsidence clustered around the Sarno village. In the same area, Valente et al. [20] suggested the occurrence of a grabenlike structure bounded by active faults, which can be inferred by geomorphological and stratigraphic analysis around the Lavorate embayment. In contrast, evidence of recent  [20,34]). White box indicates location of the study area, shown in inset map B. Active faults are derived from the ITHACA database [35] and from Cinque et al. [36]. Earthquake epicentres and magnitude are derived from the Italian Earthquake Catalogue CPTI15 [37]. (B) geological map of the inner sector of Sarno plain (modified from [20]). This stratigraphic setting hosts a multi-layered alluvial aquifer, which is recharged by local effective precipitation and, subordinately, by lateral groundwater flow coming from the Sarno Mountains karst aquifer. The groundwater circulation is directed towards the Sarno River and its tributaries [38][39][40]. A high number of uncontrolled and unauthorized private wells (up to approximately 90-100 wells/km 2 ) use the plain groundwater for industrial, domestic, and agricultural purposes [5]. Other groundwater withdrawals are carried out at three well-fields and two gravity drainage galleries located near the Sarno Mountains karst aquifer.
The inner sector of the Sarno plain has been deeply investigated by means of DInSAR analysis [1,21], which suggested intense human-induced subsidence clustered around the Sarno village. In the same area, Valente et al. [20] suggested the occurrence of a grabenlike structure bounded by active faults, which can be inferred by geomorphological and stratigraphic analysis around the Lavorate embayment. In contrast, evidence of recent vertical motions in the central and north-western sector of the Sarno Mountains pediment-Remote Sens. 2021, 13, 3323 4 of 25 alluvial plain junction area were not discussed in the previous work by Valente et al. [20] but are discussed in this paper.
Finally, the recent activity of the fault strands bounding the Sarno plain were not associated with diffuse and moderate-to-high magnitude seismicity ( Figure 1A and Supplementary Figure S1). The strongest historical seismic events affecting the Sarno plain, and its surroundings, occurred in 1499 and in 1561, both having moment magnitudes of 5.5, and whose epicentres were located~15 km to the NW and to the SE of the Sarno village urban area, respectively [36] (Figure 1).

Geomorphological Setting
The Sarno plain is an almost flat alluvial plain that is limited by the Gulf of Naples to the SW, and the carbonate ridges of the Sarno Mountains and the Lattari Mountains to the NE and to the SE, respectively ( Figure 1). The size of the plain is about 15 × 8 km, with elevation ranging from 0 to 30 m a.s.l. near the pediment of the abovementioned mountain fronts. The flat morphology of the plain is interrupted by a slight topographic high, up to 10 m high and NW-SE trending, named the San Valentino Torio high (SVT) [20], which is bounded by rectilinear scarps towards the NE. The SVT high is dissected by the Sarno river in the area between the towns of Striano and San Valentino Torio ( Figure 1B). The Sarno river is the main river flowing within the plain and is formed by the confluence of the Acqua della Foce River, the Acqua del Palazzo River, and the Acqua di San Marino River ( Figure 1B). The courses of these three rivers were locally rectified by human-made interventions, which mainly affected the Acqua della Foce and Acqua del Palazzo river courses. However, the Acqua di San Marino River is less rectified, and its course is quite articulated, with a river bend that follows the NW-SE trend of the rectilinear scarp bounding the SVT towards the NE [20].
The SVT limits, to the SW, a very low-topography area that passes, in the NE, to the pediment of the Sarno Mountains. This low-topography area, hereinafter referred to as the Sarno basin, is NW-SE oriented, with a size of 2 × 10 km, and is bounded, towards the NE, by the alluvial fans that developed in the pediment of the Sarno Mountains [20]. To the east of the Sarno basin, the Sarno Mountains are roughly NW-SE trending and interrupted by two large embayments, the Episcopio and Lavorate embayments. Both embayments are characterised by the presence of two generations of coalescent alluvial fans that pass to a debris slope towards the mountain front [20]. In the Lavorate embayment, the alluvial fans are locally interrupted by the Fossa San Vito sinkhole. This sinkhole is located in the footwall block of the fault that was proposed to bound the Sarno Mountains mountain front [19].

Materials and Methods
We adopted a multidisciplinary approach to investigate ground deformation in the inner sector of the Sarno plain. Geomorphological, stratigraphic, and structural analyses of the mountain front-pediment-alluvial plain system were coupled with GNSS and DInSAR data analysis. Moreover, hydrogeological and piezometric data of the Sarno River alluvial plain aquifer are re-interpreted.
In the following paragraphs, each methodological step will be described separately and in detail.

Geomorphological Analysis
Geomorphological analysis was carried out by combining the analysis of topographic maps at detail scale, GIS-aided LIDAR data (1 m horizontal resolution) [41], and field surveys. Large-scale topographic maps (Technical Map of the Campania Region at scale 1:5000) were investigated to identify structural controlled landforms along the southern slope of the Sarno Mountains and in the adjoining pediment and alluvial plain. Large-scale topography of the Sarno River alluvial plain was also investigated using LIDAR data in GIS software (Arcgis 10.8 ©), which allowed for the recognition of landscape features Remote Sens. 2021, 13, 3323 5 of 25 not detectable by 1:5000 topographic maps. Finally, field surveys allowed for both the verification and integration of fault-controlled landforms recognised in the previous steps.

Stratigraphic Analysis
Stratigraphic analysis consisted of the collection and re-interpretation of both published and unpublished data from 71 boreholes (Supplementary Figure S2). Published data comprised of 14 boreholes from the ISPRA database [42] and 34 boreholes from scientific papers [43][44][45], whereas unpublished data comprised of 23 boreholes from the Piano Urbanistico Comunale of the Sarno municipality. Borehole depths ranged from 15 to 315 m. The borehole distribution allowed three geological NE-SW oriented cross-sections to be constructed, which were representative of the subsurface stratigraphy in the Episcopio and Lavorate embayments and in the Sarno village urban area.

Structural Geology Analysis
To better detail the geology of the Sarno Mountains to the east of Sarno village, we carried out geological structural mapping. A published geological map of the area [46] provided the base map for determining the lithological units, and mapping was mostly aimed at the collection of structural data. Data, such as bedding dips, fault dips, kinematic information, and fault displacements, were all recorded in the field at individual points of measurement. Using Stereonet software, these data were then used as inputs to plot planes and poles to assist in the inspection of the spatial arrangement of the data. Finally, the data allowed for the structure of the area of interest to be characterised.

Hydrogeological Analysis
To analyse hydrogeological trends within the study area, we integrated published piezometric data collected between 1978 and 2003 [19,38,[47][48][49] with new data collected from September 2019 to March 2020.
Furthermore, monthly rainfall data, recorded by the rain gauge station at the Sarno village ( Figure 1B) and collected from 1992 to 2020, were also considered ( Table 1).  [50] All these data were analysed in QGIS ® and allowed for the reconstruction of three multi-temporal maps (for the years 1978, 1995, 2003) of groundwater levels (m a.s.l.; data not shown), and one map of the spatial-temporal variation of the piezometric heads. The latter map was constructed by overlapping piezometric data for the years 1992 and 2003 and by interpolating data points with the inverse distance weighted technique. We were also able to estimate the lowering groundwater level velocity (m/y) by dividing the piezometric head variation, observed between 1992 and 2003, by the time lapsed during the same period.

Geodetic Analysis
GNSS and DInSAR are two techniques of satellite geodesy, which have been widely used in recent years to measure crustal deformation associated with a number of natural phenomena, such as earthquakes, volcanism, landslides, and subsidence (e.g., [51][52][53][54][55][56][57][58]). GNSS can provide millimetre-level accuracy in horizontal directions and centimetre-level accuracy in the vertical direction [59] with unprecedented time resolution, but only at selected points, even if in dense networks. On the contrary, DInSAR techniques can be used to measure areal displacements in the direction of radar line-of-sight (LoS) with an accuracy of up to a portion of the microwave wavelength, namely an accuracy of several millimetres to centimetres [60]. For realistic crustal deformation monitoring, these two techniques perform better when used jointly. Due to the scarce availability of permanent GNSS stations in the studied area, hereafter we used GNSS only for validation purposes and cross-checking of the DInSAR data. Details of geodetic data used in this study are presented in Table 2. The available GNSS data were collected at the permanent station PACA, located in Palma Campania ( Figure 1A). The whole data set encompasses 17.3 years, between 2003 and 2020 ( Table 2). The GNSS data are a time series of 3D (north, east, up) site positions, produced by the Nevada Geodetic Laboratory [61] as daily solutions obtained using GIPSY-OASIS II software in Precise Point Positioning (PPP) mode. The solid earth and pole tide contributions were reduced according to the "IERS 2010 Conventions" [62], whereas ocean tide loading accounted for the FES2004 model [63]. The daily GNSS solutions referred to the IGS14 realisation of the International Terrestrial Reference Frame 2014 [64]. North, east, and up CGNSS components were analysed. Since the horizontal components displayed negligible rates after removal of the Eurasia plate movement, for this study we focused only on the vertical (up) component of the ground deformation.

DInSAR 1993-2020
DInSAR data were obtained from four different satellite constellations that have covered the investigated area since 1993 (Table 2). In detail, data was obtained from the ERS1/2 and ENVISAT constellations of the European Space Agency (ESA) for the period 1993-2010; the Cosmo-SkyMed satellite constellation of the Italian Space Agency (ASI) for the period 2012-2013, using and interpreting data obtained as part of the "II Not-ordinary Plan of Environmental Remote Sensing" [65,66]; and, finally, ESA s SENTINEL-1 data, for the period 2016-2020, were processed and interpreted using SUBSIDENCE software, which implements the coherent pixel technique-temporal phase coherence (CPT-TPC) approach [67,68], developed at the Remote Sensing Laboratory (RSLab) of the Universitat Politecnica de Catalunya of Barcelona.
DInSAR data produce relative LoS motion measurements with respect to a reference point. To refer the GNSS and SAR motions to the same frame, the 3D GNSS vectors defined in the local reference system have to be projected into the SAR LoS using the relationship: GNSS LoS = sin (θ) cos (φ)E + sin (θ) sin (φ)N + cos (θ) UP [69] Where, (θ) and (φ) are the incidence angle and ground track, respectively, and E, N, UP satellite are the components (east, north, and up) of the GNSS displacement.
When just the vertical components of the GNSS vectors are of interest, the projection into the SAR LoS can be performed using the simplified relationship from [69]: GNSS LoS = cosθ * UP

Geomorphological Analysis
Results of geomorphological analysis are reported in the geomorphological map of Figure 2. Here, we combined fault-controlled landforms detected in this study with landforms recognised by Valente et al. [20].  Fault scarps characterise the southern slope of the Sarno Mountains mountain front and are NW-SE trending and SW facing. Their height ranges from 250 m, to the NW of the Acqua della Foce River, to 400 m, to the NE of the Sarno urban area ( Figure 2). Furthermore, NE-SW trending, SE facing fault scarps were also recognised in the NW sector of both the Episcopio and the Lavorate embayments ( Figure 3A). All the fault scarps exhibit a planar form, being poorly dissected by the river network, with a slope also exceeding 40 • . In the Episcopio embayment, both NE-SW trending, SE facing, and NW-SE trending, SW facing triangular facets were recognised. Remote Sens. 2021, 13, x FOR PEER REVIEW 9 of 28

Tectonostratigraphy of the Sarno Plain
The detailed stratigraphic and structural setting of the surface and subsurface Quaternary filling of the Sarno plain and the Episcopio and Lavorate embayments was described based on field surveys and borehole data ( Figure 4 and Supplementary Figure 2). Field surveys confirmed the geological setting described in Section 2.1. Pivotal to our study are deep boreholes, which allowed for the detailing of the subsurface stratigraphy of the alluvial plain deposits and their structural configuration. Rectilinear scarps were recognised at the pediment-alluvial plain junction. These included a NW-SE trending, SW facing, 3 to 5 m high scarp that cut the second-generation alluvial fan in the Lavorate embayment, already recognised by Valente et al. [20]. Furthermore, field work allowed recognition of another rectilinear scarp, which occurs to the north of the previous one. It is NW-SE trending, SW facing, 3-5 m high and has been recognised in the urban area of Sarno village ( Figure 3B,C).
The abundance of fault-controlled landforms suggests the presence of a complex fault pattern, thus enhancing the role of tectonics in shaping the mountain front. Further detail on the fault systems that affect the Sarno mountains in this area is provided in Section 4.3.

Tectonostratigraphy of the Sarno Plain
The detailed stratigraphic and structural setting of the surface and subsurface Quaternary filling of the Sarno plain and the Episcopio and Lavorate embayments was described based on field surveys and borehole data ( Figure 4 and Supplementary Figure S2). Field surveys confirmed the geological setting described in Section 2.1. Pivotal to our study are deep boreholes, which allowed for the detailing of the subsurface stratigraphy of the alluvial plain deposits and their structural configuration. Borehole data in the Lavorate embayment (cross section A−A' in Figure 5) indicated that the embayment is filled with ~100 m of deposits. They consist of ~25 m of secondgeneration alluvial fan deposits passing downward to the 25−30 m thick CI, to the 25−30 m thick pre-CI deposits, and to the 15−30 m thick marine to transitional deposits. The above-described sequence lies on the carbonate substratum that is located at −50 m a.s.l. In the adjacent alluvial plain, the sequences described above are thicker, with thickness exceeding 200 m, and the carbonate substratum is not found up to −250 m a.s.l. These data suggest the presence of a SW-dipping fault, named F1 in Figures 4 and 5, that initially downthrows the carbonate substratum, with a minimum throw of ca. 200 m. Fault F1 was also recognised by the geomorphological analysis (Section 4.1). Considerable change in the thickness of the marine to transitional deposits across fault F1 ( Figure 5) indicated synkinematic accumulation of this sequence. This is suggestive of an evolutionary motion model in which fault F1 initially downthrows the top of the carbonate substratum with a synchronous, syn-faulting deposition of the marine to transitional deposits. The fault was then sealed by the pre-CI deposits and the CI and was then reactivated in more recent times and displaced them with a much smaller throw (ca. 15 m). Furthermore, in the hanging wall block of fault F1, a 5 m thick peat layer is present in the uppermost portion of the alluvial plain deposits. Both pieces of evidence suggest a recent motion of fault F1. The cross-section continues to the SW up to the SVT. In this area, subsurface data indicated offsets of both the CI, the pre-CI, and the marine to transitional deposits. These data al- The subsurface stratigraphic and structural setting of the study area was synthesized in three geological cross-sections, labelled A-A , B-B and C-C (Figures 4 and 5). Stratigraphic units were grouped as follows (from the base to the top): carbonate substratum; marine to transitional deposits; pre-Campania Ignimbrite (pre-CI) alluvial deposits; Campania Ignimbrite (CI); alluvial plain deposits and post-Campania Ignimbrite (post-CI) pyroclastics; alluvial fan deposits; and lacustrine, peat, and travertine layers interlayered within the post-CI deposits.
Borehole data in the Lavorate embayment (cross section A-A in Figure 5) indicated that the embayment is filled with~100 m of deposits. They consist of~25 m of second-generation alluvial fan deposits passing downward to the 25-30 m thick CI, to the 25-30 m thick pre-CI deposits, and to the 15-30 m thick marine to transitional deposits. The above-described sequence lies on the carbonate substratum that is located at −50 m a.s.l. In the adjacent alluvial plain, the sequences described above are thicker, with thickness exceeding 200 m, and the carbonate substratum is not found up to −250 m a.s.l. These data suggest the presence of a SW-dipping fault, named F1 in Figures 4 and 5, that initially downthrows the carbonate substratum, with a minimum throw of ca. 200 m. Fault F1 was also recognised by the geomorphological analysis (Section 4.1). Considerable change in the thickness of the marine to transitional deposits across fault F1 ( Figure 5) indicated syn-kinematic accumulation of this sequence. This is suggestive of an evolutionary motion model in which fault F1 initially downthrows the top of the carbonate substratum with a synchronous, syn-faulting deposition of the marine to transitional deposits. The fault was then sealed by the pre-CI deposits and the CI and was then reactivated in more recent times and displaced them with a much smaller throw (ca. 15 m). Furthermore, in the hanging wall block of fault F1, a 5 m thick peat layer is present in the uppermost portion of the alluvial plain deposits. Both pieces of evidence suggest a recent motion of fault F1. The cross-section continues to the SW up to the SVT. In this area, subsurface data indicated offsets of both the CI, the pre-CI, and the marine to transitional deposits. These data allowed us to hypothesise the presence of another fault, named F2 in Figures 4 and 5, associated with downward increasing throw values from a few metres near the surface to up to ca. 50 m at the top of the marine to transitional deposits. This fault may be also envisaged by the rectilinear scarp that bounds the SVT in this area and by anomalies in the river network ( Figure 2).

Structural Geology of the Sarno Mountains
The carbonate substratum rocks, cropping out to the east of Sarno village, are wellbedded with bedding planes that are, overall, gently to moderately north-dipping ( Figure  6c). The bedding is cross-cut by at least three fault systems displaying different orientations.
The most dominant of these is associated with hundred-metre-long fault segments that are, overall, oriented NW-SE (Figure 6d), roughly parallel to the above-described faults F1 and F2 (see Section 4.2). Although it was not possible to find kinematic indicators associated with this fault system, we interpreted it to be associated with dip-slip extensional kinematics on the basis of literature data from adjacent areas [31,71].
Another fault system that was identified in the field is oriented WNW-ESE and is made up of sub-vertical faults, which, in some cases, steeply dipped toward the NNE (Figure 6d). Field evidence suggested that this fault system could be associated with mostly strike-slip kinematics. The best example of a ca. 1 km long fault segment, belonging to this system, is cropping out immediately east of Sarno village ( Figure 6). Finally, a  Figure 4. Cross-section A-A is centred in the Lavorate embayment and it is modified from Valente et al. [20]. Cross-section B-B is centred in the Episcopio embayment; cross-section C-C is centred in the Sarno urban area.
Cross-section B-B in the Episcopio embayment ( Figure 5) enhanced a similar subsurface setting as that of the Lavorate embayment. Unfortunately, the poor detail of deep boreholes, which describe the presence of a~100 m thick sequence of undistinguished gravelly-sandy layers with abundant pyroclastics and tuff units, do not allow the subsurface stratigraphic setting to be detailed. Only the alluvial fan deposits may be envisaged by geomorphological analysis, even if it is not possible to define the thickness of this unit. Furthermore, it is possible to constrain the depth of the carbonate substratum in the Episcopio embayment that is located at ca. 0 m a.s.l. Similar to section A-A , the SW-ward deepening of the top of the carbonate substratum could be ascribed to the presence of the SW-dipping fault F1.
Cross-section C-C ( Figure 5) allowed for the subsurface setting of the Sarno urban area to be detailed. Here, borehole data indicated a more complex fault configuration in which the carbonate substratum was downthrown across SW-dipping fault segments that partition the bulk displacement in the footwall of fault F1. In detail, one fault could be identified near the NE end of cross section C-C , which bifurcates upward in two fault segments. These fault segments are, however, unlinked along the strike, as shown in Figure 4 (see Camanni et al. [70] for similar fault geometries); each of these downthrow the top of the carbonate substratum of ca. 20 to 30 m. Furthermore, in the footwall of fault F1, the marine to transitional deposits are not present, thus testifying that this area was a topographic high during the deposition of this unit. Peat layers and travertine accumulated both in the footwall and in the hanging wall of fault F1. A thin layer of lacustrine deposits also accumulated at the base of the post-CI deposits, in the hanging wall of fault F1. Furthermore, the reconstructed CI top surface is NE-dipping, both in the hanging wall of fault F1 (elevation ranging from 0 m a.s.l. in borehole B24 to −5 m a.s.l. in borehole B27), and between fault F1 and its NE splay (elevation ranging from 8 m a.s.l. in borehole B51 to −1 m a.s.l. in borehole B28). These data testify to fault activity in the last 39 ka. Activity, very recently, of both fault F1 and its NE splay was also indicated by the north-eastward thickening of the alluvial plain deposits and post-CI pyroclastics, which can be interpreted as the effect of syn-tectonic deposition of this unit. Furthermore, very recent activity of the NE splay of fault F1 was also manifested by the offset of 3-5 m thick peat layers placed at the top of the post-CI deposits. This small splay is located in the area where geomorphological analysis allowed for recognition of a NW-SE trending, SW facing, 5 m high rectilinear scarp that cleaves the urban area of Sarno village (Figures 2 and 3B,C). These data suggest that the rectilinear scarp may be interpreted as the surface expression of this fault activity.

Structural Geology of the Sarno Mountains
The carbonate substratum rocks, cropping out to the east of Sarno village, are wellbedded with bedding planes that are, overall, gently to moderately north-dipping ( Figure 6C). The bedding is cross-cut by at least three fault systems displaying different orientations.
The most dominant of these is associated with hundred-metre-long fault segments that are, overall, oriented NW-SE ( Figure 6D), roughly parallel to the above-described faults F1 and F2 (see Section 4.2). Although it was not possible to find kinematic indicators associated with this fault system, we interpreted it to be associated with dip-slip extensional kinematics on the basis of literature data from adjacent areas [31,71].
Another fault system that was identified in the field is oriented WNW-ESE and is made up of sub-vertical faults, which, in some cases, steeply dipped toward the NNE ( Figure 6D). Field evidence suggested that this fault system could be associated with mostly strike-slip kinematics. The best example of a ca. 1 km long fault segment, belonging to this system, is cropping out immediately east of Sarno village ( Figure 6). Finally, a roughly E-W fault system was also recognised ( Figure 6D), whose kinematics are unknown. For each of the aforementioned fault systems, an accurate definition of the displacements of individual fault segments was not straightforward, largely due to the lack of markers within the carbonate substratum rocks.
Mapping also allowed us to collect bedding and fault data within the Quaternary deposits that flank the SW side of the Sarno Mountains in at least two locations around the Lavorate embayment ( Figure 6). The bedding is overall E-W oriented and mostly shallowly dipping towards the south, as it was mapped within the Quaternary slope deposited on the south facing mountain slope ( Figure 6A). However, one bedding data point displayed a northward dip. This is associated with a very recent fault that developed within these Quaternary deposits east of Sarno village and folded the bedding, determining this anomalous bedding dip. This fault belongs to a system comprised of faults oriented mostly NE-SW ( Figure 6B), which were interpreted to be the most recent within the mountain front because it deforms Quaternary deposits. Offsets of beds across segments belonging to this fault system indicated primarily extensional displacements, with values spanning from a few centimetres to a few tens of metres.

Hydrogeological Analysis
The groundwater level evolution in the study area during the period of 1992-2003 is depicted in Figure 7. Due to strong groundwater withdrawal from the plain aquifer, the water table suffered a general piezometric lowering. According to [72], the maximum values (up to −5m) of the lowering of the water table level were observed in the north-western and south-western sectors of the study area and the urban area of Sarno village. Finally, cross-cutting map evidence seemed to suggest that NW-SE faults are cut by E-W faults in the NW of Episcopio, and that the recent NE-SW faults in the Quaternary deposits are cutting across the WNW-ESE sub-vertical fault to the east of Sarno village ( Figure 6). However, the field data collected for this study did not allow for the inference of any definitive indication of the relative chronology between the recognised fault systems, an aim which was far beyond the scope of this article.

Hydrogeological Analysis
The groundwater level evolution in the study area during the period of 1992-2003 is depicted in Figure 7. Due to strong groundwater withdrawal from the plain aquifer, the water table suffered a general piezometric lowering. According to [72] The groundwater lowering rate varied between 0.45 and 0.09 m/yr. The highest rates, occurring between 1992 and 2003, were concentrated in the central portion of the study area, around the town of Sarno and in the sector between the F1 and F2 faults.
As the study period was characterised by relatively stable amounts of monthly rainfall (Figure 8), the groundwater lowering likely resulted from other causes. The most probable of these could be intense groundwater pumping from the plain aquifer, which is equivalent to about 30 x 10 6 m 3 /yr [72], for agricultural, domestic, and industrial purposes. As the study period was characterised by relatively stable amounts of monthly rainfall (Figure 8), the groundwater lowering likely resulted from other causes. The most probable of these could be intense groundwater pumping from the plain aquifer, which is equivalent to about 30 × 10 6 m 3 /yr [72], for agricultural, domestic, and industrial purposes. As the study period was characterised by relatively stable amounts of month fall (Figure 8), the groundwater lowering likely resulted from other causes. Th probable of these could be intense groundwater pumping from the plain aquifer, is equivalent to about 30 x 10 6 m 3 /yr [72], for agricultural, domestic, and industri poses.

Geodetic Data (GNSS and DInSAR) Analysis
The next two sections are devoted to the overview of the main results obtained from the analysis of the geodetic data sets.

GNSS Analysis for the Period 2003-2020
The up component of the GNSS dataset (2003-2020) clearly showed that the PACA permanent station was relatively stable, with little vertical deformation ( Figure 9a) and a low rate of subsidence experienced (about 1 mm/yr), which can be at the limit of significance for the accepted accuracies for this methodology. For a more accurate determination of the vertical ground motion, a first order linear trend was estimated from the CGNSS data. The results clearly depicted a periodic trend (Figure 9a), which is a well-known feature of the continuous GNSS time series [73]. Indeed, long-period noise in GNSS position time-series data is ubiquitous and often the most limiting factor preventing proper detection of subtle ground deformation in the geodetic signals. The continental hydrology loading was proven to be the largest contributor to surface deformation at seasonal and even longer timescales [74]. We estimated the elastic ground displacement, due to the hydrological loading over the Sarno plain, using the MERRA2 (modern-era retrospective analysis for research and applications) [75] and GLDAS (global land data assimilation system) [76] hydrological models. Computation of the ground displacement at the selected sites used in this study was made available at the EOST loading service [77]. The output grid of these models has, for MERRA2, a spatial resolution of 0.5 • × 0.625 • latitude by longitude and an hourly temporal resolution, while GLDAS has a spatial resolution of 0.25 • × 0.25 • and a three-hour temporal resolution. As the Sarno plain covers an area of 15 × 8 km, a single grid cell of the MERRA2 or GLDAS models covers the entire area.
The annual phase and the amplitude of 5-7 mm of elastic displacements, provided by both MERRA2 and GLDAS models, match the de-trended geodetic time series of the PACA CGNSS station. The first-order linear fit preliminarily applied to the CGNSS time series indicated a mean geodetic velocity of -1.28 ± 0.03 mm/yr (Figure 9a). After periodic trend removal with the hydrology model, a more reliable mean subsidence rate was obtained (Figure 9b). Regardless, even after the removal of the periodic trend, GNSS data provided evidence of stability or very low subsidence (<1.5 mm/yr) for Palma Campania.

DInSAR Analysis for the Period 1993-2020
The DInSAR data used for this study covered the period of 1993-2020. In particular, the ERS1/2 and ENVISAT datasets ( Figure 10A,D) were first analysed in both acquisition geometries (ascending and descending). As can be seen in Figure 10, the coverage of targets in the study area were well distributed, characterised by a density of about 60 targets/km 2 for the ERS1/2 dataset and about 250 targets/km 2 for the period of 2002-2010 (ENVISAT). In particular, for the periods of both 1993-2000 and 2002-2010, displacement rates in some sectors even higher than 10 mm/yr were recorded. Subsequently, the analysis was extended to the period of 2012-2013 due to the availability of high-resolution data acquired by the Cosmo-SkyMed satellite constellation ( Figure 10E). Notably, the density of targets was much higher than previously, with about 2600 targets/km 2 .
Regarding displacement rates, again, values in some sectors were found to be greater than 10 mm/yr, particularly near the Sarno village urban area.
Moreover, the observation window for displacement rates was extended to 2020 using an additional 138 SAR images from SENTINEL-1, collected in the period of 2016-2020. Figure 10F shows the mean displacement rate for SENTINEL-1 images.
As with the other periods analyzed, it can be seen that the largest rates of displacement were, once again, concentrated within the faults F1 and F2. Although the SENTINEL-1 images were also in C-band and characterised by medium resolution, it was nevertheless possible to identify a high density of targets, i.e., about 1500 per km 2 . Figure 9. GNSS daily solutions collected at PACA permanent station located in Palma Campania for the up component. Linear trend, which quantifies the subsidence rate, is plotted (red line). A tentative modelling of the hydrological contribution of the observed ground displacement based on MERRA2 and GLDAS models is plotted as well (a). GNSS daily solutions, after correction of the hydrological contribution and final assessment of the subsidence rate. The black brackets mark the time range for which the DInSAR dataset is available for comparison (b).

DInSAR Analysis for the Period 1993-2020
The DInSAR data used for this study covered the period of 1993-2020. In particular, the ERS1/2 and ENVISAT datasets (Figure 10a-d) were first analysed in both acquisition geometries (ascending and descending). As can be seen in Figure 10, the coverage of targets in the study area were well distributed, characterised by a density of about 60 targets/km 2 for the ERS1/2 dataset and about 250 targets/km 2 for the period of 2002-2010 (ENVISAT). In particular, for the periods of both 1993-2000 and 2002-2010, displacement rates in some sectors even higher than 10 mm/yr were recorded. Subsequently, the analysis was extended to the period of 2012-2013 due to the availability of high-resolution data acquired by the Cosmo-SkyMed satellite constellation (Figure 10e). Notably, the density of targets was much higher than previously, with about 2600 targets/km 2 . Linear trend, which quantifies the subsidence rate, is plotted (red line). A tentative modelling of the hydrological contribution of the observed ground displacement based on MERRA2 and GLDAS models is plotted as well (a). GNSS daily solutions, after correction of the hydrological contribution and final assessment of the subsidence rate. The black brackets mark the time range for which the DInSAR dataset is available for comparison (b).

Combination of GNSS and DInSAR datasets
To validate the interferometric results, because no GNSS data were available in the Sarno plain, we used the GNSS data collected at the PACA permanent station, which is located a few kilometres from the study area ( Figure 11). The comparison between ENVISAT and GNSS time-series data is displayed in Figure 12. It is worth noting that, to obtain the vertical component identified by the interferometric results, the two acquisition geometries (ascending and descending) were decomposed according to procedures presented in the literature [78,79]. As a result, both the DInSAR and GNSS detected a ground displacement rate of the same amplitude (0.4-0.5 mm/yr) even if this was at the limit of detectability.  Figure 4). located a few kilometres from the study area ( Figure 11). The comparison between ENVI-SAT and GNSS time-series data is displayed in Figure 12. It is worth noting that, to obtain the vertical component identified by the interferometric results, the two acquisition geometries (ascending and descending) were decomposed according to procedures presented in the literature [78,79]. As a result, both the DInSAR and GNSS detected a ground displacement rate of the same amplitude (0.4-0.5 mm/yr) even if this was at the limit of detectability.   Figure 13 shows the time series extracted from the DInSAR images for three selected points, which were deemed more suitable to characterise the kinematics of the study area. The DInSAR data were up components, projected according to the suitable satellite ground tracks. From the comparison between the GNSS and DInSAR time series, we derived a confidence limit for the DInSAR data. For the period when both GNSS and DIn-SAR data are available, namely, 2003-2010 (Figure 9b), the two techniques provided evidence on the stability of the PACA site located on the Apennine foothill, or a low-rate subsidence, whereas a different behaviour was evident from the DInSAR data for the Sarno plain.
Therefore, after validating the interferometric results in the study area, it was possible to reconstruct the whole time series of displacements, by calculating the vertical components, for the period of 1993-2020 ( Figure 13). It is evident that the area was subject to subsidence phenomena, as was well known in the literature [1,21]. Some targets, belonging to specific sectors within the Sarno basin, were selected (SRN1, SRN2, and SRN3 in Figure 11). From the time series, it is possible to observe an almost linear trend ( Figure  13). The most relevant displacements were observed in the central-northern sector of the  Figure 13 shows the time series extracted from the DInSAR images for three selected points, which were deemed more suitable to characterise the kinematics of the study area. The DInSAR data were up components, projected according to the suitable satellite ground tracks. From the comparison between the GNSS and DInSAR time series, we derived a confidence limit for the DInSAR data. For the period when both GNSS and DInSAR data are available, namely, 2003-2010 (Figure 9b), the two techniques provided evidence on the stability of the PACA site located on the Apennine foothill, or a low-rate subsidence, whereas a different behaviour was evident from the DInSAR data for the Sarno plain. Remote Sens. 2021, 13, x FOR PEER REVIEW 21 of 28 Figure 13. Time series of the DInSAR data (up component) used to characterise the vertical displacements within the Sarno basin (a-c). Refer to Figure 11 to locate the selected points. Therefore, after validating the interferometric results in the study area, it was possible to reconstruct the whole time series of displacements, by calculating the vertical components, for the period of 1993-2020 ( Figure 13). It is evident that the area was subject to subsidence phenomena, as was well known in the literature [1,21]. Some targets, belonging to specific sectors within the Sarno basin, were selected (SRN1, SRN2, and SRN3 in Figure 11). From the time series, it is possible to observe an almost linear trend ( Figure 13). The most relevant displacements were observed in the central-northern sector of the plain (SRN1 and SRN2), where the highest rates of piezometric depressions were recorded during the period of 1992-2003 (Figure 7). The central sector (SRN2) was found to be the most affected by subsidence phenomena, with a cumulative displacement of about -180 mm and a rate of nearly 6 mm/yr (Figure 13b), whereas the area less affected by subsidence, namely the south-eastern sector (SRN3), had subsidence rates of the order of 3 mm/yr (Figure 13c).

Discussion
Overall, the data confirmed the presence of a graben-like structure, the Sarno basin, in the inner sector of the Sarno plain, which was first suggested by Valente et al. [20]. This area was the locus of relevant subsidence, both in the long-term (as suggested by the accumulation of around 20 m of lacustrine, travertine, and peat layers in the last 39 ka), and in the short-term (as suggested by geodetic data for the 1992-2020 time span).
The comparison of all the datasets allowed for the distinguishing of the natural (e.g., tectonics) and larger anthropogenic (e.g., groundwater pumping) contributions to subsidence within the Sarno basin.
The tectonic contribution to subsidence was indicated by geomorphological and stratigraphic datasets (Sections 4.1 and 4.2). Evidence of recent tectonic activity include the rectilinear scarps that ruptured the surface. These occurred both to the NE (e.g., in the surroundings of the Lavorate embayment and in the Sarno village) and to the SW (e.g., along the northern sector of the SVT) of the Sarno basin. Associated with these rectilinear scarps are, at depth, offset CI and post-CI deposits. As the CI is chronologically constrained at 39 ka [80], the fault slip rate during the last 39 ka may be derived. Valente et al. [20] estimated the slip rates for both faults F1 and F2, which were named, in that paper, faults F3 and F4, respectively. Slip rates in the last 39 ka were found to be 0.2 mm/yr for fault F1 and 0.1 mm/yr for fault F2. Notably, the slip rate for the fault F1 was estimated near the Lavorate embayment. In the Sarno urban area, fault F1 caused a~10 m vertical offset of the top surface of the CI. This would result in an average slip rate of~0.25 mm/yr during the last 39 ka, slightly larger than that seen in the adjoining Lavorate basin. It must also be considered that the top surface of the CI may have experienced erosion after its deposition, thus the value of the slip rate may be slightly different. Unfortunately, borehole data in the hanging wall block near the Sarno village urban area does not reach the base of the CI. As a consequence, we could not evaluate the offset of the CI base surface, which could have provided a more robust estimation of slip rate.
The slip rate variation along fault F1 was consistent with the often-observed variation in offset values along a normal fault [81], which are maxima in the central portion of the fault (e.g., a slip rate of 0.25 mm/yr in the Sarno urban area) and decrease towards its tips (e.g., a slip rate of 0.1-0.2 mm/yr in the Lavorate embayment). Furthermore, slip rate values around 0.1-0.25 mm/yr are consistent with the slip rate values detected in the inner sector of the Southern Apennines [40,82].
These data in combination suggest recent activity of the fault strands bounding the Sarno basin towards the NE (fault F1 in Figures 4 and 5) and the SW (fault F2 in Figures 4 and 5). Indeed, these faults are small and incapable of generating large earthquakes and breaking the surface in one single seismic event. The rectilinear scarps, however, can be interpreted as resulting from multiple low to moderate energy slip events. Furthermore, the orientation of the F1 and F2 faults is coherent with that of the more pervasive fault system cropping out in the adjacent Sarno Mountains (Section 4.3). Anthropogenic contribution to subsidence was indicated by the geodetic datasets and the hydrogeological data. DInSAR datasets also suggested the occurrence of a cluster of subsidence within the Sarno basin (Section 4.5). Hydrogeological data showed that piezometric variation was spread across the Sarno plain, with variable rates within the Sarno basin (Section 4.4). Furthermore, high values of lowering of the water table levels were observed both within the Sarno basin and outside of it (e.g., the area with values of −3/−5 m in Figure 7). Worthy of note is the fact that subsidence clustered only within the Sarno basin, thus enforcing the idea that subsidence occurred where unconsolidated deposits were present in the subsurface of the Sarno basin.
Furthermore, Figure 14 depicts the groundwater withdrawal rates at SF, MP, and SM well-fields (refer to Figure 7 for locations), located within Sarno Mountains karst aquifer, and monthly rainfall during the 2003-2020 period. The temporal variations of both hydrogeological parameters did not show significant anomalies but, rather, smooth trends. Consequently, we reasoned that the observed vertical ground deformation was not triggered by a decrease in the natural groundwater recharge of the plain aquifer because the precipitation regime did not display any substantial change during the same interval of time.
by a decrease in the natural groundwater recharge of the plain aquifer because the precipitation regime did not display any substantial change during the same interval of time.
This field evidence presents new results that were not reported by or discussed in previous studies [21,49,83]. These findings allow for a better understanding of the anthropogenic cause of the ground deformation mechanisms observed in the plain aquifer.
Indeed, even in the absence of updated monitoring data on the groundwater table and pumping rates from the alluvial aquifer, these results indirectly confirmed that the land subsidence observed in the period of 1993-2020 ( Figure 13) was mainly attributable to the widespread and uncontrolled groundwater abstraction, for industrial and irrigation use, within the plain aquifer. In a retrospective analysis, this hydrogeological interpretation is reasonably consistent with a strong link between the trend of ground deformation and groundwater table depletion observed in the period of 1992-2003.
Our results suggested that a positive feedback exists between fault activity and the distribution of subsidence. In fact, fault activity caused the formation of the Sarno basin and the consequent accumulation of unconsolidated deposits. The presence of these deposits amplified the deformation effect induced by groundwater pumping during the past 30 years, thus resulting in the significant subsidence rate detected by geodetic data, which is clustered within the Sarno basin.  This field evidence presents new results that were not reported by or discussed in previous studies [21,49,83]. These findings allow for a better understanding of the anthropogenic cause of the ground deformation mechanisms observed in the plain aquifer.

Conclusion
Indeed, even in the absence of updated monitoring data on the groundwater table and pumping rates from the alluvial aquifer, these results indirectly confirmed that the land subsidence observed in the period of 1993-2020 ( Figure 13) was mainly attributable to the widespread and uncontrolled groundwater abstraction, for industrial and irrigation use, within the plain aquifer. In a retrospective analysis, this hydrogeological interpretation is reasonably consistent with a strong link between the trend of ground deformation and groundwater table depletion observed in the period of 1992-2003.
Our results suggested that a positive feedback exists between fault activity and the distribution of subsidence. In fact, fault activity caused the formation of the Sarno basin and the consequent accumulation of unconsolidated deposits. The presence of these deposits amplified the deformation effect induced by groundwater pumping during the past 30 years, thus resulting in the significant subsidence rate detected by geodetic data, which is clustered within the Sarno basin.

Conclusions
We carried out a multidisciplinary analysis of the inner sector of the Sarno plain, in the Southern Apennines of Italy, with the aim of characterizing the subsidence experienced in both the short-and long-term in this area of low to moderate seismicity. We approached this problem by combining the geomorphological, stratigraphic, structural, hydrogeological, GNSS and DInSAR datasets. Comparison of all datasets allowed us to discriminate between the natural (i.e., tectonics) and the anthropogenic (i.e., groundwater pumping) contributions to subsidence in the Sarno basin.
We identified a sector of the inner Sarno plain, named the Sarno basin, that experienced enhanced subsidence both in the long-and short-term. Geomorphological, stratigraphic, and structural data suggest that the Sarno basin is bounded by two faults (i.e., faults F1 and F2 in Figures 4 and 5) that show evidence of activity during the last 39 ka, and slip rates ranging from 0.1 to 0.25 mm/yr. We named these faults the Sarno Fault (fault F1) and the San Valentino Torio Fault (fault F2). Both are active faults, with the Sarno Fault also showing evidence of surface rupture, thus allowing for the interpretation that it is an active and capable fault. Notably, neither the Sarno fault nor the San Valentino Torio fault were included in any official national catalogues of active tectonic structures, e.g., the ITHACA [39] and DISS [84] databases.
The geodetic data indicated that the Sarno basin is also currently experiencing subsidence, which has occurred for at least the past 30 years. This subsidence has a rate of 3-5 mm/yr and is localized in the area between the bounding faults of the basin. Hydrogeological data suggested that a significant part of this subsidence can be ascribed to groundwater pumping from the alluvial plain aquifer.
We also determined that the fault activity and observed subsidence rate are related. In fact, fault activity caused the formation of the Sarno basin and the consequent accumulation of unconsolidated deposits, which are currently acting as amplifiers of subsidence that is mainly induced by groundwater exploitation.
Finally, our multidisciplinary approach was proven to be reliable for the study of ground deformation in urbanized areas. In principle, it could be applicable globally. In the inner Sarno plain, further investigation might provide more data about active tectonics. These data might include new boreholes across both the Sarno and the San Valentino Torio faults, chronological analysis on datable stratigraphic layers (e.g., 14 C analysis of peat layers) and paleoseismological analysis.