Analysis of Deformation Dynamics in Guatemala City Metropolitan Area Using Persistent Scatterer Interferometry

: The analysis of deformation dynamics in Guatemala city and its surrounding region presented in this paper holds signiﬁcant relevance due to the high vulnerability of this area to natural disasters, combined with its rapid urbanization, similar to most Central American cities, contrasting with a lack of InSAR and deformation studies in the region. A total of 226 SAR images from Sentinel-1 A and B satellites in both ascending and descending geometries were processed with the Persistent Scatterer Interferometry (PSI) technique employing the SNAP-StaMPS integrated processing chain. The study area encompasses the Metropolitan Region of Guatemala, which is characterized by a diverse and active geological framework, with a historical record of earthquakes, intense groundwater extraction, and local subsidence phenomena, causing ﬁssures and sinkholes. Four active areas were identiﬁed in the study area, each covering more than 50 hectares, with subsidence velocities greater than 10 mm/yr. This study provides valuable insights into fostering the sustainable development of this region by identifying deformation patterns, characterizing main active areas, and evaluating associated risks for disaster management and prevention. The results can also aid informed decision-making processes and guide urban planning and resource management strategies in other Central American countries. The application of InSAR studies is crucial for improving safety and sustainability in urban environments and natural resource management in vulnerable regions.


Introduction
Synthetic Aperture Radar Interferometry (InSAR) is a cost-effective and highly accurate technology for studying the deformation dynamics of surfaces [1][2][3]. InSAR has been extensively used to analyze land deformation resulting from geologic and seismological hazards [4][5][6], as well as landslide risk studies [7][8][9][10] and mining activities [11,12], among others. Moreover, InSAR has proved to be a valuable tool for urban development strategies aimed at improving sustainability and safety in urban environments [13][14][15][16][17][18][19][20][21]. Additionally, InSAR has also been used to analyze urban deformation caused by the exploitation of groundwater [1,[22][23][24][25]. For instance, Mexico City is one of the fastest sinking metropolises in the world, with a displacement rate of up to 2.5 cm/month. Many urban areas have experienced these phenomena in varying degrees; thus, subsidence is becoming a major factor in urban development planning [1]. Jakarta has sunk by about 4 m, prompting Indonesia to build a new capital. This subsidence began affecting permanent constructions as early as 1978, causing such phenomena as cracks, expansion of flooding areas, lowering of groundwater level, and increased inland seawater intrusion [26]. Another example is the Semarang case [17], where excessive groundwater extraction and natural consolidation of alluvial soil led to land subsidence, with displacement velocities reaching up to 16.5 cm/yr. In summary, InSAR has proven to be a versatile tool for studying and monitoring various types of surface deformations, including those caused by natural and anthropogenic factors, making it an essential tool for urban development strategies.
Persistent Scatterer Interferometry (PSI) [2,27] is a technique that relies on the identification of point-like features on the Earth's surface, called Persistent Scatterers (PSs), which have a strong and stable radar backscatter signal. PSI focuses on stable and coherent scatterers such as buildings and infrastructure, providing high-precision measurements of surface deformation.
PS techniques are particularly abundant in built-up areas, allowing this technique to achieve high spatial resolutions of deformation dynamics in urban areas. By quantifying the phase difference between the radar signals from different images, PSI can detect small changes in the distance between the satellite and the persistent scatterers, allowing for the detection and monitoring of ground subsidence, landslides, and other surface deformations with high accuracy. PSI is then particularly useful for monitoring urban deformation processes that can have a significant impact on infrastructure, property, and public safety. PSI can also provide valuable information for urban planning and management, allowing authorities to identify areas of potential risk and take appropriate action. Overall, the use of PSI in urban areas offers a cost-effective and efficient way to monitor surface deformation and improve the safety and sustainability of cities. Despite the extensive use of InSAR in several parts of the world, such as Europe, there is a lack of such studies in Central American countries. Only a few studies have utilized InSAR to study surface deformations in Guatemala, specifically, in volcanology [28,29], tectonics [30], and urban environments [17,18,20,23,29,31].
The lack of InSAR studies in Central American countries is particularly concerning given the region's high vulnerability to natural disasters and anthropogenic activities [32]. For instance, several cities in the region are located in areas prone to earthquakes, landslides, and volcanic activity, increasing the likelihood of surface deformation events. Furthermore, rapid urbanization and population growth in the region have led to an increase in infrastructure development and extraction of natural resources, which can contribute to land subsidence and other forms of surface deformation. Therefore, the use of InSAR in Central American countries can provide valuable insights into surface deformation dynamics, inform risk assessments, and guide decision-making processes for sustainable urban development and natural resource management.
This study analyzes the deformation dynamics, with PS InSAR, of part of the Metropolitan Region of Guatemala, known officially as Región Metropolitana de Guatemala (RMG) (Figures 1 and 2). This area is integrated by Guatemala city (formed by 22 zones) and four adjacent cities: Santa Catarina de Pinula, Petapa, Villanueva, and Mixco. The RMG is influenced by natural factors, such as geological features, extreme climate with intense precipitation, and socio-demographic factors, which makes it an interesting case study. The purpose of this paper is to emphasize the significance of using InSAR for an accurate and comprehensive analysis of deformation in urban areas to ensure sustainable development. Examining and evaluating risks is an initial step towards reducing their impact, as acknowledged by numerous international reports on Disaster Risk Reduction (DRR) and sustainable development [33][34][35].
DRR plays a crucial role in promoting social and economic growth and ensuring a sustainable future. Central American urban areas are geologically young volcanic areas. The physiographic architecture of Central America is primarily determined by the northwestern course of the Central American trench and the Central American volcanic front, which were formed by Cenozoic subduction [36]. Soils in the region, combined with heavy rainfall, topographic complexity, and a seismically active area, can be prone to great instability during earthquakes and heavy precipitation [37]. The Population Census XII and the Housing Census VII reveal that 14,901,286 people reside in Guatemala, with a yearly population growth rate of 1.8% between 2002 and 2018. The Guatemala Department, which includes Guatemala City and 16 other municipalities, houses the most significant percentage of the country's population, 20.2% of the total [38]. This population growth has resulted in increased infrastructure development and exploitation of natural resources, such as underground water. Between 2017 and 2019, 2-3 million square meters of construction was authorized, with the municipality of Guatemala accounting for 45% of these new construction areas. Housing projects make up a large proportion of these projects, with an increase of 1,363,777 buildings between 2002 and 2018 [39]. Construction and urban growth must be conducted within the framework provided by risk assessments, with a focus on disaster risk reduction. Land use sensitive to disaster risks, urban planning, safe construction, and sound infrastructure will not only contribute to water management, but will also protect lives and properties, and thus benefit entire cities [39]. Studies have demonstrated that water pumping for urban and agricultural use is one of the primary factors that trigger land subsidence [10,16,[40][41][42][43][44]. InSAR is an accurate tool that can be applied to groundwater management and used to monitor changes in urban areas and predict disasters to avoid economic, environmental, and human losses. The conservation and efficient use of water resources is a vital concern of the RMG government. Factors such as climate change, environmental degradation, and pollution have contributed to water scarcity in Guatemala [45,46]. Due to population growth, there has been an increase in groundwater exploitation, resulting in lower water levels in some wellfields [45][46][47]. In addition, deforestation of hydrographic basins has led to an increase in surface runoff, decreasing infiltration into the soil. The sealing of land through construction and urbanization of areas, overuse of aquifers by public and private entities, lack of regulations and control on the use of water resources, and other factors, threaten the availability of this vital resource [48]. Guatemala City is located in an area that has been affected by highly destructive earthquakes throughout history. Specifically it has experienced 21 major earthquakes, with an intensity on the Mercalli scale greater than VIII [49]. The study area contains important faults and volcanic elements, with Mixco and Pinula fault structures being the main faults located within the study area, and Santiaguito, Fuego, and Pacaya being the main volcanoes [50].
Several sinkholes and some other local events have been identified in the study area. A relevant sinkhole cited in the literature [51] occurred in the year 2007, measuring 30 m in diameter and 60 m in depth. In the short period between June and October 2022, nine sinkhole events [52], with associated collapsed buildings, were recorded in the study. As mentioned above, this study focuses on diagnosing the patterns of deformation in the main urbanized zones of RMG, and aims to locate and quantify the severity of the deformation phenomenon in the area.
The use of InSAR technology is being considered as a potentially economically sustainable approach to detect areas of instability. Given the reported cases in the literature of environmental characteristics and population growth in the area, attending to this case study is deemed important. Understanding the conditions of change is crucial in meeting the needs of the affected areas, and identifying the characteristics of the most affected areas can aid in finding the causes and developing solutions to avoid future problems.

Study Area
The study area encompasses the municipalities of Guatemala City, Mixco, Villa Nueva, San Miguel de Petapa, and Santa Catarina de Pinula. They integrate the RMG (Región Metropolitana de Guatemala), spanning a total area of 47,218 hectares.
The Guatemala City valley lies at the intersection of three tectonic plates: North America, Caribbean, and Cocos. The Cocos plate subducts under the Caribbean plate in the Mesoamerican trench, giving rise to the Central American mountain range of volcanoes. This range has attained its maximum growth in Guatemala, and is characterized by calderic cones and depression outcrops, with acid volcanism represented by rhyolites, ignimbrites, and pyroclasts. The Guatemala City valley can be defined as a graben or trench bounded on the east by the Pinula and on the west by the Mixco faults [50,53]. The geological formations in the area consist of granitic intrusive rocks, rhyolitic lava flows, sediments, and tuffs, which constitute the Tertiary volcanic group. Additionally, there are Quaternary volcanic deposits in the form of tephras (pyroclasts) that culminate in alluvial deposits. The geology was reclassified into four classes: alluvium, sedimentary rock, tephra, and other volcanic rocks ( Figure 1) [50,[54][55][56][57].
The formation process of the Guatemala City valley is intricately linked with the Motagua and Jalpatagua faults, and the zone of distension that arose at their intersection. During the Middle Miocene, the Motagua and Pinula faults became increasingly active [58], resulting in the creation of a weaker zone that facilitated the rise of magma to the volcanic structures of the Pinula and El Naranjo volcanoes. While this zone initially served as a channel for magma to rise, continuous movements of the faults in the distension zone eventually led to the formation of the Mixco fault. The shape of the Mixco fault closely follows the contour of the relict structure of the Pinula volcano, and its movements induced the subsidence of the valley and its volcanic structures. The boundary between the depression and the volcanic structure is delineated by the El Trébol fault within the Guatemala City valley [50]. Guatemala City was established as the capital of the country in 1773, after the former capital of Santiago de los Caballeros was devastated by the Santa Marta earthquake [59,60] .
Analyzing isophreatic lines from 1978, 1995, and 2012 [46], a decline in the water table over the past three decades of exploitation was observed. In sedimentary rock deposits, the water table has decreased by −3 to −5.5 m/year, while volcanic aquifers have decreased by around −3.4 m/yr. Due to uncontrolled exploitation of aquifers without legal regulation, excessive groundwater extraction is taking place in the area, risking premature depletion of the reserves [46].

MT-InSAR Processing and Vectorial Decomposition
A total of 226 SAR images were processed, including both ascending and descending geometries from Sentinel-1 constellation satellites A and B ( Figure 2). All images are Single Look Complex (SLC) TOPSAR data acquired in Interferometric Wide (IW) swath mode with VV polarization, for the time frame of January 2017 to September 2021 (Table 1). The SeNtinel Application Platform (SNAP) version 9.0.0 developed by the European Space Agency (ESA), snap2stamps, and the Stanford Method for Persistent Scatterers (StaMPS) software were used following SNAP-StaMPs integrated processing for Sentinel-1 PSI [61].
To achieve satisfactory outcomes, it is crucial to appropriately configure the input parameters for StaMPS processing to prevent the occurrence of non-natural artifacts. StaMPS parameters significantly influence the final results. The characteristics of the analyzed movements impact the processing results, emphasizing the importance of appropriate settings. In our case, default parameter values led to an irregular stepped trend in several PS time series. To address this, parameters related to atmospheric filtering, phase unwrapping, and estimation of the Spatially Correlated Look Angle (SCLA) error have been set to values explained by Balbi et al. 2021 [62]. Specifically, these modified parameters introduced in the final steps of StaMPS that allowed us to obtain a more natural trend are as follows: "unwrap_time_win" set to 24 days (default 730 days) to smooth the phase in time by estimating noise for each pair of neighboring pixels; "unwrap_grid_size" set to 10 m (default 200 m) for spacing of the resampling grid; "unwrap_gold_n_win" set to 8 (default 32) for the size of the window used in the Goldstein filter; "scla_deramp" set to "yes" (default "no") to estimate the phase ramp for each interferogram; and "scn_time_win" set to 50 days (default 365) for the window size of the low-pass temporal filter.
Both geometries were referenced to a compact area with highly temporal coherence (>0.8) near the GNSS station GUAT (14.590, −90.520) [63]. The reference point is located in a flat area (slopes close to zero degrees) over the tephra geology class. Line-of-Sight (LOS) ascending and descending series were combined in a vectorial decomposition to obtain the vertical and east-west displacements. Each point in the ascending geometry is paired with the nearest neighbor point in the descending geometry, not further than a certain predefined threshold. This distance threshold has been defined according to the accuracy of the position of the PS. Under the assumption that there is no movement in the north-south planimetry direction, it is possible to decompose LOS movements on the vertical and east-west planimetry directions. For each ascending/descending pair of points, we compute decomposed movements following Equations (1) and (2). This decomposition approach has been successfully applied in previous studies [64,65]: where d asc LOS and d desc LOS are the measured displacements in the Line-Of-Sight (LOS) direction for the ascending and descending datasets, θ asc and θ desc are the looking angles for the ascending and descending LOS directions, φ asc and φ desc are the horizontal angles of the LOS directions measured from the north in a clockwise direction, and d up and d east are the decomposed mean deformations in the vertical and east directions (Figure 3). The combination of ascending and descending PS was performed in the vector domain [64,66], combining PSs based on spatial proximity. Any PS located in one geometry that lacked neighboring PSs in the opposite geometry within a designated search radius of 10 m were omitted.

Deformation Analysis and Identification of Critical Areas of Subsidence
With the aim of highlighting the existing patterns over the vertical velocity deformation maps, a Jenks classification with six classes was performed with QGIS software [67,68]. Class ranges were rounded to closer multiples of 10, taking into account the range and variability of vertical deformations in the area of study. Velocity profiles have been drawn over the areas with the largest subsidences to perform a qualitative analysis on the trends and patterns of the vertical and east-west deformation average velocities. Kriging techniques [69] have been applied to delineate the critical areas of subsidence. In the generated rasters, regions displaying subsidence rates exceeding 10 mm/yr were carefully identified and outlined as areas with the most significant subsidence patterns. The administrative zones and the Global Human Settlement Built-up area product [70] were employed to accurately quantify the extent and severity of the impact on each of these administrative zones.
In order to delve deeper into the analysis and identify consistent patterns of deformation, a comprehensive examination was conducted within the study area. Specifically, relevant Areas of Interest (AOIs) were delineated and characterized, aiming to pinpoint distinct zones exhibiting significant subsidence. These AOIs were defined as areas where vertical deformations exceeded 10 mm/yr and encompassed an area larger than 50 hectares.
The deformation patterns of each AOI were examined independently and in conjunction with other variables. Deformation maps were analyzed jointly with slope (%), elevation (m.a.s.l.), aspect, and geology class. For this purpose, we used a DTM from Alos Palsar-1 [71] and the geology map classification developed by Intituto Geográfico Nacional in Guatemala [54][55][56][57]. This map was reclassified on the four main geology classes that can be found in the area: sedimentary rocks, alluvium, tephra, and other volcanic rocks.
In each of the AOIs, areas representing a subsidence trend have been selected. Subsequently, the time series of deformation for these areas were analyzed to assess variations in subsidence trends and accelerations during the specified period.

Results
A total of 269,205 PSs were obtained after the decomposition into vertical and east-west directions. A point density of 189 PS/km 2 was obtained for the urban surface inside the study area (Figure 4). According to the Global Human Settlement Built-up area product [70], 48% of the analysis area is urbanized, while 51% is vegetated.
The majority of PSs are localized in constructed zones, as dense vegetation hinders the obtention of PSs derived from C band data. The data also show that 3.9% (891 ha) of the urbanized area included in the study area is affected by deformations greater than 10 mm/yr ( Table 2).
The most affected areas by deformation in terms of surface are Petapa and Mixco, located in the south and northwest of the study area, respectively, with surfaces that reach 355 and 185 ha of urban area. Santa Catarina de Pinula, Villanueva, and Zone 5 (Guatemala City) also reach subsidences of over 10 mm/yr in areas larger than 50 ha (Table 2).   Figure 5 presents the distribution of the PS vertical velocities for the main geology classes. This distribution shows that the largest subsidence is located in the alluvium class. The presence of numerous and widely scattered outliers for sedimentary, tephra, and volcanic rock indicates a higher degree of variability and potentially complex deformation patterns in these geology classes.
The majority of the alluvium class is situated in the municipality of Petapa, which coincides with the area where the highest deformation magnitudes have been observed. Non-meaningful results were obtained in the comparison of deformation patterns with the rest of the analyzed variables, such as slope, elevation, and aspect.

Areas of Interest: Critical Areas of Subsidence
Four different AOIs have been identified and characterized, all of them with subsidence velocities greater than 10 mm/yr on areas exceeding 50 ha ( Figure 6). These areas are located within the following zones: (i) Mixco and Zone 19 (AOI 1), (ii) historical area (AOI 2), (iii) Petapa and Villanueva (AOI 3), and (iv) Santa Catarina de Pinula (AOI 4). In all these areas subsidence velocities exceed 20 mm/yr, reaching up to 6 cm/yr in Petapa/Villanueva. The largest areas affected by subsidences over 10 mm/yr are located in AOIs 1 (Mixco/Zone 19) and 3 (Petapa/Villanueva), which reach 216 and 440 ha, respectively, representing 1% and 2%, respectively, of the total urbanized surface in the RMG.   Figure 9 presents the corresponding profiles of vertical and horizontal deformation for the selected cross sections. Subsidence bowllike patterns, or depression cones, with a radial displacement in planimetry towards the center of maximum subsidence, were detected in AOI 1 (Mixco and Zone 19), AOI 2 (historical area), and AOI 4 (Santa Catarina de Pínula). In AOI 1, three subsidence zones were identified; the one at Zone 19 exhibits a subsidence bowl pattern, while in the Mixco subsidence area, western PSs were not obtained due to dense vegetation, hindering the confirmation of a subsidence bowl pattern. However, we observed westward movements in the eastern part of that region. The third one, located south from Zone 19, shows net subsidence.   In AOI 2, the subsidence and east-west deformations expand over the entire AOI, indicating either a wide subsidence bowl or a combination of diverse subsidence processes, as shown in the profiles in Figure 9. AOI 4 also shows a pattern compatible with a subsidence bowl or depression cone pattern. In these AOIs (1, 2, and 4), the subsidence velocities reached 20 mm/yr, while the planimetry magnitudes remained below 13 mm/yr.
The results for AOI 3 show that the deformation points are mainly located within the boundaries of the Amatitlán caldera, with a significant concentration of PSs displaying an eastward movement within the Amatitlán caldera and along the southern part of the Mixco fault system [60]. In this AOI, the planimetry magnitudes remained below 13 mm/yr, as in the other AOIs, but subsidence velocities reached up to 60 mm/yr.

Discussion
The obtained PS density is in agreement with other studies carried out in urban areas using PSI InSAR [72]. As expected, information gaps or low-PS density areas correspond to regions with vegetation or complex topography. In situations where subsidence is the predominant process, it can be postulated that the horizontal displacement is directly proportional to the first spatial derivative of the vertical deformation [66,73]. By applying this hypothesis, the horizontal displacement reaches zero in the maximum subsidence point and in the edges of the subsidence bowl-like patterns. Conversely, it attains maximum values where the first spatial derivative of vertical deformation is at its maximum, which aligns with observations made in other cases worldwide [23,44,66,74,75].
In AOIs 1, 2, and 4, we identified a pattern coherent with these subsidence bowl deformation processes, or depression cones, associated with groundwater extraction. This pattern is characterized by convergent horizontal deformations (eastward movements in the western margin and westward movements in the eastern margin) towards the areas of maximum subsidence. North-south displacements cannot be estimated due to the geometry of the image acquisition and further decomposition; therefore, the movements along the north-south axis could be in the same order of magnitude as the east-west movements, fitting the theoretical 3D cone shape. These conical strains have maximum deformation velocities in the center of the areas that taper to the edges. This phenomenon was documented by Hu (2009) [76] in the Yang-tze delta, revealing a discernible correlation between groundwater extraction and a temporal and spatial pattern of land subsidence. A previous study in Guatemala [23] confirms the patterns appreciated in the present study, compatible with depression cones caused by aquifer extraction. Zhu et al. (2015) [44] identified an area in northern Beijing with deformations of 55 mm/yr in the vertical component, partially matching the theoretical distribution of a groundwater subsidence cone. This study attributes the irregularities to the unequal thickness of the deposits. However, groundwater volume analysis using the InSAR technique is not fully feasible without accurate knowledge of the lithology and stratigraphy of the study area [16]. Calibration data are also required for the validation of groundwater management models. The need for consistent information is one of the major challenges in hydrogeology [77]. Kim Figure 10) . There are also two subsidence areas located in administrative Zones 10 and 14 (south of AOI2), consistent with [23]. However, our results do not show any meaningful subsidence in central administrative Zone 16 (eastern AOI 2) due to a lack of PSs.
Deformation time series for AOI 1 ( Figure 10) shows a sudden increase in subsidence velocity starting from October 2018 to May 2019. This area is affected by the Mixco fault system [50,53,60]. The change in the time series coincides in time with two registered earthquakes of over 5.4 Mw reported on October 12th and 18th 2018 [78,79].
The deformation times series of AOIs 2 and 3 show trend changes after October 2018: time series of AOI 2 decelerate, meanwhile, time series in AOI 3 accelerate. It seems considerably important that one must analyze the possible effects of seismic events as a trigger of deformation dynamics. The velocity profile of AOI 2 ( Figure 9) shows a subsidence bowl pattern in the first part of the profile. The distribution of the vertical deformations in the profile shows different overlapping areas of subsidence, which suggests a combination of different driving processes (water extraction, geological activity, anthropogenic works, etc.).
In AOI 3, the alluvial sector in the south of the study area presents higher values of subsidence, coinciding with an area of water extraction. The deformation pattern is intense and more uniform than in the other classes of geology. This behavior can be related to three phenomena: first, water withdrawal in Ojo de Agua, which is the group of water wells that extracts around 72% of the groundwater from the aquifers of the central valley [80,81]; second, alluvium areas happen to be prone to seismic hazards [82][83][84][85]; third, its location at the boundaries of the Amatitlán caldera and the presence of associated faults [60].
The characterization of the deformation patterns in the RMG opens up new avenues for investigating their origins and enhances urban and infrastructure management practices. Further research in this area would enable the identification and anticipation of local events, such as landslides and sinkholes, as well as the monitoring of infrastructure during construction and operation. Developing a historical database of such incidents would facilitate the application of machine learning algorithms to predict future events and plan appropriate measures to mitigate potential harm and losses.
Moreover, the integration of data from other SAR sensors (such as L-band and Xband) and combination with other InSAR algorithms [86] based on distributed scatterers would prove advantageous, particularly in regions with challenging topography or dense vegetation. L-band SAR data could improve the monitoring in vegetated areas, whereas distributed scatterer processing would provide results in the less coherent areas. Conducting small-scale investigations that incorporate physical and environmental data, along with detailed GNSS information, would further support future research efforts and enhance the validity of our findings.
As an illustrative example, we refer to a recent news report by García and Kestller in September 2022 [87], highlighting how several buildings in the municipality of Petapa have been affected. Petapa municipality reached maximum subsidence velocities exceeding 6 cm/yr in our study. This demonstrates the real-world implications and significance of studying deformation patterns for the betterment of local communities and infrastructure management.

Conclusions
The results of this study indicate that the application of InSAR in the Central American region can provide valuable insights into surface deformation dynamics and associated risks. InSAR also provides valuable information for urban planning and management, allowing authorities to identify areas of potential disaster and take appropriate action and/or planning measures.
This research confirms that multi-temporal InSAR PSI processing is a competent method to characterize deformation patterns in the Metropolitan Region of Guatemala (RMG). Nonetheless, the intricate topography and dense vegetation in the area posed challenges in accurately characterizing the deformation patterns. Further research is necessary to combine data from different SAR bands, such as using the vegetation penetration capabilities of L-band and combining PSI and distributed scatterer algorithms, to complete the characterization of the identified deformation patterns. Critical areas in the RMG were identified for the monitored period. The delineation was based on qualitative analysis and Jenks classifications. Future studies will focus on clustering and machine learning methods considering time series of relevant variables and time-precise variations based on registered events (such as earthquake records).
Additional research is needed to characterize some of the observed deformations, such as the subsidence bowl patterns, which are currently depicted in the literature as depression cones induced by water withdrawal processes. Furthermore, these patterns should be verified with the water extraction and precipitation data registers.
It is also recommended that we analyze the influence of tectonics and geological factors in the observed deformation patterns, excluding other driving forces, such as anthropogenic activities. The study area is traversed by dense geological structures (faults, calderas, etc.) and complex river networks, making it difficult to characterize the relationships between deformation dynamics and causative factors. In the case of tectonics, a reliable database of events and affected areas and structures is required to conduct a proper analysis. The collection and consolidation of these data is a challenge due to the lack of resources or governance structure in some countries.