Monitoring of surface movements in Cardona has mainly been performed trough GNSS/GPS surveys and DInSAR data, complemented with levelling surveys.

The GNSS/GPS monitoring started in 1997, being repeated in a yearly basis. The current situation of the local network of Cardona is that shown in Fig. 2. Two local vertexes have been used as base points, the so-called CK29 and CK30, which belong to the CK Network covering the whole Potash Basin (Gili et al., 2010). The other 41 points are control points. As the network aims to detect general movements in the whole area (big units) and not sudden/localized movements (small sinkholes or local slides, for instance), the density of control points is relatively low (inter-point distances ranging between 200 m and 1 km).

A GNSS/GPS precision system has been used, specifically the Fast-Static Method (also called “GPS with post-processing”). The recording rate (Δt) was 15 s in the first campaigns, and 5 s since October 2010. The elevation angle threshold is 10^∘ (registering the signal from the “low” satellites tends to improve the vertical precision of the final elevations, our main goal). Duration of the logging sessions are typically 30 min per control point. With this methodology, the a priori precision (σ) of the results is 1 cm in plan and 2 cm in elevation.

The DInSAR monitoring was performed using images from the ERS, ENVISAT, COSMO-SkyMed, and SENTINEL satellite missions. ERS and ENVISAT data (1992–2010) have a spatial pixel resolution of 40×40 m and a revisit time of 35 and 30 d. COSMO-SkyMed, data (2010–2015) has a spatial pixel resolution of $\sim \mathrm{15}\times \mathrm{15}$ m and an irregular revisiting time from 4 to 16 d, depending on the available satellites. Finally SENTINEL data (2015–2019) has a spatial pixel resolution of 20×20 m and a constant revisit time of 6 d, which allow a better monitoring of the processes occurring in Cardona. During the last decades, DInSAR techniques have been widely improved, developed and applied (Blanco-Sànchez et al., 2008; Ferretti et al., 2000), proving to be an excellent tool for monitoring regional land surface movements of wide areas with sub-centimetric precision. In this study, an advanced DInSAR processor, known as Persistent Scatterer Interferometry (PSI), was applied. This process discriminate the different components of the interferometric signal, obtaining the final ground motion measurement with a higher precision (Mora et al., 2007b). One of the main drawbacks of this processing technique is that it requires a low variation of the radar signal response during the analysis period in order to obtain high quality measurements (interferometry coherence). As a consequence, soil covers experimenting great variations along the year return poor coherence. This limitation is the reason why this technique allows the measurements of surface movements mainly in urban areas or isolated buildings, but is limited in areas with important surface changes in time such as vegetated areas or areas with human-induced changes, such as landfill/excavation activities. Land use change limitation is notorious on the western section of the Cardona area, reducing significantly the data available in this area. To corroborate the DInSAR-detected movements, a high-precision levelling network was implemented in Cardona urban area. The levelling works began on January 2006 and finished on October 2012 after thirteen campaigns, and were performed using an automatic electronic (digital) level, measuring elevation increments along some traverses (Fig. 2).

It should be kept in mind that each monitoring network has a different resolution, extension, advantages and limitations, showing some differences on the measured points. Therefore, it is necessary to determine the results of each technique individually and later discuss the comparison between different techniques, in order to obtain unique coherent interpretation of the monitored movements.

GNSS/GPS measurements show two areas with different behaviours along the Vall Salada lineation.

Figure 3Terrain movement velocities obtained from the 20 year movement dataset of GNSS/GPS measures.

The area to the west, next to the Muntanya de Sal, presents the greatest subsidence (negative values) of the area, with rates up to −12 mm yr⁻¹. On the contrary, on the east section of the Vall Salada, a heave (positive) sustained velocity near to 20 mm yr⁻¹ is detected in five points between the old (Terrera Vella) and the new (Terrera Nova) mine waste areas (Fig. 3).

Figure 4Levelling based and DInSAR based terrain movement velocities determined from movement measures gathered in its respective datasets.

DInSAR measures also deliver a differential behaviour. Greater subsidence is detected at the south-west of the urban centre of Cardona, with values ranging between −10 and −15 mm yr⁻¹. On the contrary, heave values between 10 and 20 mm yr⁻¹ are measured on the north eastern section of Vall Salada (Fig. 4).

Levelling campaign measures indicate subsidence of the urban area of Cardona with values ranging from ∼0 to 7.5 mm yr⁻¹, which are consistent with DInSAR measures in this area (Fig. 4).

GNSS/GPS and DInSAR results are consistent in the identification of two main areas with different terrain movement trend. Both techniques show a subsiding area located on the south and western areas of the Muntanya de Sal with rates between 10 and 15 mm yr⁻¹, and a maximum heave between 10 and 20 mm yr⁻¹ on the east section of Vall Salada. In both directions of movement the two techniques show some differences on the location of the maximum subsidence and heave points. These differences can also be detected between DInSAR and levelling data. The most plausible reason for these differences is related with the technique's resolution. GNSS/GPS and levelling techniques are point based, meaning its results are restricted to a single point of measure. On the contrary, DInSAR data represents an extension of terrain (pixel), which varies with the satellite's sensor resolution (from 40×40 to 15×15 m). Consequently, in DInSAR data higher or lower punctual displacements are masked by the general pixel movement, while GNSS/GPS and levelling control points might be in locations with higher or lower movement than its surroundings. Anyhow, results between techniques are consistent considering its resolution and precision (Gili et al., 2010; Marturià et al., 2006).

The area where GNSS/GPS and DInSAR show subsidence is located on the geological contact between the Cardona formation (salts) and the Igualada marls formation, above the old Nieves mine tunnels. This location explains the subsidence movements, as the contact presence infers ground debility and the excavated tunnels introduce a change in the ground tension's equilibrium. This change induces deformations of the surrounding ground (strain) producing a tension (stress) rearrangement to reach a new mechanical equilibrium (Zangerl et al., 2003). These changes are reflected on the surface by generating subsidence, which in absence of new processes of erosion or solution of ground materials, tend to evolve to stability as the new ground equilibrium is reached (Pipia et al., 2007).

In addition to the western areas of the Muntanya de Sal,SENTINEL data also show subsidence on the Terrera Nova area. In this location GNSS/GPS control points were not installed as this landfill deposit has been progressively removed due environmental issues. Despite this land cover change, the surface materials are similar and the DInSAR coherence has been kept. Consequently, the movements on this section detected with this technique cannot be considered as a surface movement.

Regarding the area where GNSS/GPS and DInSAR show heave, this is located between the Terrera Vella and the old Cardener river bed. In addition, DInSAR data also show a heave area on the section close to the old river bed, where no GNSS/GPS control points are installed. This area is located at the axis of the diapiric anticlinal formation showing the outcrop of the evaporitic units (salt gypsum and anhydrite). As this is an active diapiric process (Masana and Sans, 1995; Wagner et al., 1971), and the heave is stable in time, the halokinesis is considered the main cause generating heave. Anyhow, the exposure of the evaporitic lithologies together with the known issues with the water infiltrations (Lucha et al., 2008), can induce to the eventual rehydration of the gypsum and anhydrite, which can contribute to the ascending movement (Yilmaz, 2001).

The results of GNSS/GPS and DInSAR techniques indicate that the latest Earth observation missions based on DInSAR (SENTINEL), combined with advanced processing techniques, greatly improves the measure interval, spatial resolution, spatial coverage and revisiting time from previous DInSAR data, making this technology an excellent option to monitor areas with good signal response. On the other hand, data based on GNSS/GPS data deliver high quality punctual information which is extremely useful to determine local trends on any type of terrain and land cover.

The data used in this work is proprietary of the government of Catalonia through the Institut Cartogràfic I Geologic de Catalunya (ICGC) and third-party suppliers and clients. In particular cases, some of this data can be made available by ICGC under request. Queries may be addressed to ICGC mailbox https://www.icgc.cat/en/The-ICGC/Contact/Contact-mailbox (last access: 19 February 2020).

JM supervised the different ICGC projects from which the data was obtained. JAG, NL and CPP conceived and planned the location points and frequency of the GPS/GNSS measurements, carried out the GPS/GNSS field campaigns and processed the data. XRL and JM compiled and processed the DInSAR and leveling Data and contributed to the interpretation of the results. XRL took the lead in writing the manuscript. All authors provided critical feedback and helped shape the research, analysis and manuscript, and read and agreed to the published version of the work.

The authors declare that they have no conflict of interest.

This article is part of the special issue “TISOLS: the Tenth International Symposium On Land Subsidence – living with subsidence”. It is a result of the Tenth International Symposium on Land Subsidence, Delft, the Netherlands, 17–21 May 2021.