An Investigation of Volcanic Ground Deformation Using InSAR Observations at Tendürek Volcano (Turkey)

: Although approximately one-tenth of the world’s population lives near volcanoes, most of the 1500 active volcanoes are not monitored by ground-based instruments because of the cost and difﬁculty of access. Since the development of interferometric synthetic aperture radar (InSAR) in the 1990s, recent advances have allowed the near-real-time detection of surface deformations, one of the earliest markers of volcanic activity. According to the Global Volcanism Program, Turkey’s most recent eruption (involving gas and ash) occurred in the Tendürek volcano in 1885. An explosion in the Tendürek volcano, which continues to actively output gas and steam, would be a critical issue for the life and property of the people living nearby. In this context, we processed the Sentinel-1 data collected by the European Space Agency using the Stanford Method for Persistent Scatterers, and the surface deformations of the Tendürek volcano were investigated. In addition, we applied two different atmospheric correction approaches (linear phase-based tropospheric correction and the Generic Atmospheric Correction Online Service for InSAR) to reduce atmospheric effects and found that the linear phase-based tropospheric correction model produced lower standard deviation values. Subse-quently, the mean deformation velocity maps, displacement time series, and deformation components in the line-of-sight direction were calculated. The results showed that the most severe subsidence was − 11 mm/yr on the upper slopes of the Tendürek volcano. Although the lower slopes of the subsidence region have a lower settlement rate, the subsidence has a peak-caldera-centered location.


Introduction
Natural disasters have negative socioeconomic effects on societies, causing severe threats to people, other living things, and the natural and built environment [1]. Although the effects of volcanic eruptions have been relatively minor compared to other natural disasters in the last two centuries, the eruption of the Pelee volcano (France) in 1902 caused 30,000 deaths, the eruption of the Nevado del Ruiz volcano (Colombia) in 1985 caused 25,000 deaths, and the eruption of the Eyjafjallajokull volcano (Iceland) in 2010 caused the closure of European airspace and losses of more than $1.7 billion in revenue [2][3][4][5]. Furthermore, it is estimated that more than 29 million people worldwide live within 10 km of approximately 1500 Holocene volcanoes, which have been active within the last 10,000 years, and the rapidly increasing world population and expansion of urban centers also increase the exposure of societies to volcanic eruptions [6][7][8]. Although improvements in volcanic activity forecasting and monitoring have shown that the knock-on effects of volcanic eruptions can be minimized and/or avoided, a significant proportion of Holocene volcanoes are not regularly monitored by ground-based methods [8][9][10]. Monitoring surface deformation, which is one of the earliest signs of volcanic activity and can be detected by ground-based methods (e.g., tiltmeter and GPS), requires installation in remote and inaccessible areas, and the temporal and spatial resolution of the measurements is limited [11]. Moreover, recent advances in interferometric synthetic aperture radar (InSAR) have enabled the creation of deformation maps with high temporal and spatial resolution, supporting the detection of unrest in many volcanoes [12,13]. In addition to its frequent use for investigating volcanoes [14][15][16], InSAR has also proven to be a powerful geodetic tool for monitoring ground deformations caused by earthquakes [17][18][19], landslides [20][21][22], and land subsidence [23][24][25][26], owing to its ability to cover large areas, work in all weather conditions, and operate all day. However, the use of InSAR is greatly limited by temporal decorrelation, atmospheric effects, and orbital errors [27]. To overcome these limitations, researchers have developed time series analysis, namely multi-temporal InSAR techniques, e.g., [28][29][30]. Permanent scatterer interferometry (PSI), which is one of these techniques and was first proposed by [28], defines scatterers as exhibiting stable scattering behavior over time; however, this technique is limited in nonurban areas because of the low PS density [31,32]. Unlike other PSI techniques, the Stanford method for persistent scatterers (StaMPS) can be used for all terrains because it uses the spatial correlation of the interferometric phase to detect surface deformations occurring in non-urban areas and does not require a temporal deformation model describing how the deformation changes over time [31,33]. Considering these advantages, StaMPS is frequently used for the long-term monitoring of ground deformations in volcanic areas where other PSI approaches are less successful. In previous studies, StaMPS has been applied worldwide for monitoring surface deformations caused by volcanoes, including the Alcedo volcano in Ecuador [33], Piton de la Fournaise in France [34], the Changbaishan Tianchi volcano in China [35], the Santorini volcano in Greece [36], the Shinmoe-dake volcano in Japan [37], the Medicine Lake volcano in the United States [38], and Mount Etna in Italy [39].
In Turkey, which ranks 14th among 95 countries in terms of volcanic threat and has 10 volcanic systems that have been active within the last 10,000 years according to the Global Volcanism Program, there are limited studies in which the surface deformations caused by volcanoes are monitored using InSAR. In one of these studies, ref. [40] used ERS satellite images and single and double-pass interferometric techniques to detect surface deformations in the Tendürek volcano. As a result of the study carried out between 1993 and 2000, a total of 10-20 cm subsidence was detected at the summit of the Tendürek volcano. Another study performed in the same area, reported by [41], aimed to detect the surface deformations of the Tendürek volcano using 41 Envisat images obtained from descending and ascending orbits between 2003 and 2020. In their study, using the small baseline subset technique, an 11 mm/yr subsidence was reported as a result of the narrowing of a sill located 4.5 km below the volcano peak. The most comprehensive study on surface deformations caused by volcanoes in Turkey was carried out by [42]. They aimed to determine the surface deformations that occur in 10 volcanic systems in Turkey using approximately 5000 Sentinel-1 images obtained from the ascending and descending orbits between 2015 and 2020. The study mainly focused on the seasonal snow cover that reduces the interferometric consistency encountered when monitoring the volcanoes with the InSAR technique and the atmospheric problems that cause misinterpretation of deformation. In the study, in which the multi-time InSAR technique was used, it was reported that no deformation on the cm scale was detected in the 10 volcanic systems in Turkey, which shows that there are significant differences among the results from the limited number of studies on the Tendürek volcano. Considering the studies above, the last detailed InSAR study of the Tendürek volcano was performed with data from 13 years ago [41]. In another subsequent research, it was examined on a national scale, and a detailed study of the Tendürek volcano was not carried out [42]. Therefore, the present study offers two novelties in terms of the literature. The first of these is the use of up-to-date data and a detailed study specific to the Tendürek volcano. The second innovative aspect of the present study is that the permanent scatterer interferometry technique is used for the first time in this volcano. In addition, the current surface deformation maps and findings obtained from the Tendürek volcano, as a result of this study, will contribute significantly to the literature.
In this context, the main purpose of the present study is to determine the current surface deformations of the Tendürek volcano using StaMPS and to reveal the evolution of surface deformations over 5 years using the SAR time series. A total of 150 Sentinel-1 images were obtained from the ascending and descending acquisition geometries between 2018 and 2022, and surface deformations of the Tendürek volcano were determined using the SentiNel application platform (SNAP), snap2stamps, StaMPS, and atmospheric InSAR noise-reduction toolbox (TRAIN). In addition, we applied two different atmospheric correction approaches (linear phase-based tropospheric correction and the General Atmospheric Correction Online Service (GACOS) for InSAR) to minimize the errors caused by atmospheric interactions. Finally, the current deformation velocity maps in the line-of-sight direction (LOS), the displacement time series, and the components of the deformation (horizontal vertical) in the LOS direction were calculated for the Tendürek volcano and its vicinity.

Study Area
The Tendürek volcano is located between 39 • 22 N and 39 • 46 N latitudes and 43 • 70 E and 44 • 05 E longitudes, covering an area of approximately 600 km 2 at an altitude of 3514 m above sea level ( Figure 1a). It is located in the Eastern Anatolia region of Turkey, with Çaldıran located to the south of the volcano, Diyadin to the northwest, and Dogubeyazıt to the northeast. Approximately 2 million people live within 100 km. The last eruption (gas and ash) of the Tendürek volcano occurred in 1855. maps and findings obtained from the Tendürek volcano, as a result of this study, will contribute significantly to the literature. In this context, the main purpose of the present study is to determine the current surface deformations of the Tendürek volcano using StaMPS and to reveal the evolution of surface deformations over 5 years using the SAR time series. A total of 150 Sentinel-1 images were obtained from the ascending and descending acquisition geometries between 2018 and 2022, and surface deformations of the Tendürek volcano were determined using the Senti-Nel application platform (SNAP), snap2stamps, StaMPS, and atmospheric InSAR noise-reduction toolbox (TRAIN). In addition, we applied two different atmospheric correction approaches (linear phase-based tropospheric correction and the General Atmospheric Correction Online Service (GACOS) for InSAR) to minimize the errors caused by atmospheric interactions. Finally, the current deformation velocity maps in the line-of-sight direction (LOS), the displacement time series, and the components of the deformation (horizontal vertical) in the LOS direction were calculated for the Tendürek volcano and its vicinity.

Study Area
The Tendürek volcano is located between 39°22′ N and 39°46′ N latitudes and 43°70′ E and 44°05′ E longitudes, covering an area of approximately 600 km 2 at an altitude of 3514 m above sea level ( Figure 1a). It is located in the Eastern Anatolia region of Turkey, with Çaldıran located to the south of the volcano, Diyadin to the northwest, and Doğubeyazıt to the northeast. Approximately 2 million people live within 100 km. The last eruption (gas and ash) of the Tendürek volcano occurred in 1855.  Although Tendürek volcano has a polygenic structure, it is a shield volcano located within the strike-slip Balık Gölü fault zone ( Figure 1c). The fault passing through the south of Tendürek volcano extends from Iran to the Diyadin plain [44]. However, the stratigraphic and lithological units that crop out in the study area range from Paleozoic to Quaternary in age [45]. The oldest rock outcropping in the area is the Paleozoic-aged metamorphic units composed of mica schist, marble, and quartzite. These units are overlain by Permo-Carbonifer-aged crystallized limestone, which is located around Dogubeyazıt. The basement units are overlain unconformably by Eocene-aged flysch. These units are tectonically overlain by late Cretaceous-aged ophiolitic melange-containing Cretaceousaged limestone blocks. The basement units are unconformably overlain by Mioceneaged units represented by limestone, sandstone, and mudstone alternations. The abovementioned units are covered by Late Quaternary-aged Tendürek volcanic rocks. The Tendürek volcano continues to Dogubeyazıt in the north and reaches Çaldıran in the south and many preserved crater structures are observed in this area [46]. The volcanic activity in the region has continued for 250,000 years and has been divided into five different phases [46]. In addition, gas vents and thermal water springs observed in the region demonstrate that volcanic activity continues [44].

Dataset
Sentinel-1 radar data operating in the C band (~5.6 cm wavelength) was obtained in a joint venture between the European Space Agency (ESA) and the European Commission and we used the data to detect surface deformations. Sentinel-1 data, which can be used to visualize large areas with short revisit times (6 or 12 days) and capture images in all weather conditions day and night, is available free of charge from the ESA open access portal. For this study, 76 images in descending acquisition geometries and 74 images in ascending acquisition geometries, captured between May 2018 and October 2022, were used. SAR acquisition times were 03:09 UTC for the descending orbit and 15:01 UTC for the ascending orbit. Details of the SAR data are given in Table 1. In addition, the GACOS delay products corresponding to each SAR acquisition for the GACOS correction were downloaded by sending a request to the official GACOS website (http://www.gacos.net, accessed on 13 February 2023) along with the study area location and SAR data-acquisition times.

StaMPS
StaMPS, which was developed by considering previous PSI techniques, has been frequently used in the analysis of crustal movements [33,47]. This method uses the spatial correlation of the interferometric phase to detect pixels with low phase changes in all terrain structures, with or without prominent objects such as buildings and bridges [48,49]. The process consists of three main steps: interferogram generation, PS selection, and atmospheric filtration. First, appropriate sub-swath and bursts are selected in the SAR dataset, and precise orbit files are applied to them. Then, the master image is selected considering the time interval, vertical baseline, and least atmospheric effect, and all dependent images are recorded using the master image with the help of the S-1 Back Geocoding operator. After the coregistration process, interferograms are produced and prepared for StaMPS analysis by removing the topographic phase component. The entire processing block has been performed using SNAP and the freely available snap2stamps scripts developed by [50].
The next steps of the processing block are performed in the StaMPS software, and in the first step, PS candidates are selected based on the amplitude distribution index. The amplitude distribution index is defined as the ratio of the standard deviation of each pixel to the average amplitude value and is calculated using Equation (1) [33,51]: where D A defines the amplitude dispersion value, and σ A and µ A define the standard deviation and mean of the amplitude values [33]. The threshold value used is typically between 0.4 and 0.42 [52]. After determining the candidate pixels according to the amplitude distribution, the phase noise is estimated and removed for each candidate pixel using phase analysis. The residual phase of the xth pixel and the ith interferogram is given by Equation (2) [33]: where W{·} defines the wrapping operator.
x,i represent the general terms for phase change due to movement of the pixel in the satellite LOS direction, phase change due to atmospheric delay between satellite passes, residual phase change due to satellite orbital error, phase change due to viewpoint, and phase noise, respectively [33,52]. Then, candidate PS points are filtered according to the noise features estimated and extracted in the previous step so that PS points in the time series are determined based on the amplitude and phase analysis for each pixel in each interferogram. After filtering the PS points, the wrapped phase is corrected for spatially uncorrelated look angle (SCLA) error, and then the phase unwrapping step is performed. In the last step, SCLA errors mostly caused by digital elevation model (DEM) errors are predicted and removed, after this step atmospheric filtering is performed, and the LOS direction velocities for the region of interest are obtained.
In this study, a total of 150 Sentinel-1A SAR images from the descending and ascending acquisition geometries acquired between May 2018 and October 2022 were obtained from ESA and used to detect the surface deformations caused by the Tendürek volcano. The master image selection from the SAR dataset was performed by considering temporal (<1000 days) and vertical baselines (<200 m) ( Figure 2). In addition, a one-arc-second shuttle radar topography mission (SRTM) (approximately 30 m × 30 m) was used to extract the topographic phase from the interferograms. StaMPS parameters used in the processing of SAR images are given in Table 2 [33]. Finally, we test two different approaches (linear phase-based tropospheric correction and GACOS) within the TRAIN to predict and correct errors caused by atmospheric interactions. A flowchart for the methodology that was used in this study is presented in Figure 3.

Atmospheric Filtering 2.4.1. Linear Phase-Based Tropospheric Correction
To reduce the error caused by atmospheric interactions, we used the linear phase-based model available in the TRAIN, a free module developed by [53].
where ∆∅ tropo , K∆ ∅ , h, and ∆∅ 0 refer to interferometric tropospheric phase contribution, coefficient relating phase to topography, topography, and a constant phase shift, respectively [54]. After StaMPS processing is complete, linear phase-based correction is applied as the final step. More details on these steps can be found in the references [53,55,56].

GACOS Correction
We used GACOS tropospheric delay products as a second approach to predict and correct errors caused by atmospheric interactions. Developed by Newcastle University in England in June 2017, GACOS is based on numerical weather models [57]. GACOS has become frequently used by many researchers to mitigate or eliminate the effect of atmospheric error, owing to its advantages such as near real-time usability, simplicity, and free implementation [58][59][60][61][62][63]. In the processing of GACOS delay products, SRTM DEM (90 m, S60-N60) and ASTER GDEM (90 m, N60-N83, S60-S83) data are used, which are also high-definition weather models, created by the European Centre for medium-range weather forecasts [64]. GACOS uses an iterative tropospheric decomposition model to separate stratified and turbulent components from tropospheric total delays, and high spatial resolution zenith total delay (ZTD) maps are generated. Then, the ZTD maps obtained at different dates are converted to LOS interferograms and applied after phase correction to minimize or eliminate the effect of atmospheric errors on the interferograms. Finally, the SCLA error is removed after the GACOS-based atmospheric correction and as a result deformation maps and time series in the LOS direction are produced. The procedure of GACOS-based atmospheric correction approaches is shown in Figure 3. In addition, further details about these steps can be found in the references [53,55,56].

Results
The Sentinel-1 images obtained from ascending and descending acquisition geometries between May 2018 and October 2022 and the mean deformation velocity maps and standard deviation values in the LOS direction produced as a result of linear phase-based atmospheric correction are shown in Figure 4. The red-colored areas on the average deformation rate map show that they are moving away from the satellite, whereas the blue-colored areas are approaching the satellite. This situation indicates the subsidence of a region centered on the summit caldera with a maximum of −10 mm/yr in the LOS direction, in both the ascending and descending acquisition geometries. In addition, the standard deviation values of the mean deformation maps obtained from descending and ascending acquisition geometries are 0.762 and 1.185 mm/yr, respectively (Figure 4c,d).
The mean deformation velocity maps and standard deviation values in the LOS direction produced as a result of atmospheric correction based on GACOS are shown in Figure 5. Although the GACOS and linear phase-based atmospheric correction results show a similar deformation pattern, the GACOS-based results detect greater subsidence rates. In addition, the GACOS-based standard deviation values of the mean deformation maps from descending and ascending acquisition geometries are 0.873 and 1.139 mm/yr, respectively.
Results in the LOS direction from both atmospheric correction approaches show a similar deformation pattern with a caldera-centered −10 mm/yr subsidence. However, the linear phase-based atmospheric correction approach exhibits lower standard deviation values in the areas with maximum subsidence. Therefore, we use the linear phase-based atmospheric correction results in the subsequent analyses.
To examine the deformations in the Tendürek volcano, transverse and longitudinal sections perpendicular to the volcano were selected, and profile analyses were conducted (Figure 6a,b). Profile analyses show that the deformation rate, which shows maximum subsidence at the summit of the volcano, decreases linearly with distance from the summit of the volcano (Figure 6c-f). In addition, the profiles in both the descending and ascending acquisition geometries show a bowl-shaped pattern that continues with a smaller slope from the summit of the volcano to the eastern flanks. Although the E-W directional profiles of the ascending acquisition geometries show a stronger inclination than the E-W directional profiles of the descending acquisition geometries, the results from both the descending and ascending datasets show the subsidence of a similar region centered on the summit caldera.
Time series show the same findings in more detail than the mean deformation velocity maps. Figure 7 shows the deformation time series of four certain PS pixels in ascending and descending orbits at Tendürek volcano. The results show that the Tendürek volcano subsidence gradually evolves between 2018 and 2022, and the temporal evolution of the displacement is linear.   Deformation measurements made in areas where the topography changes suddenly, such as Tendürek volcano, may be adversely affected by atmospheric effects, slope, and satellite LOS directions, and the pixel distributions in the resulting datasets may not be homogeneous. Although two different atmospheric correction approaches were applied in this study, atmospheric artifacts were completely removed, and therefore, sufficient point density could not be established in some regions, especially in the center of the caldera. For this reason, loss velocity values were created using the biharmonic spline interpolation (MATLAB ® 4 grid data method) method. In addition, the results obtained in the LOS direction from ascending and descending acquisition geometries show similar deformation patterns. This suggests that vertical deformation is dominant and greater than horizontal deformation [65,66]. However, in order to detect this situation and to determine the dominant deformation pattern in the study area to better the discussion of the analysis results, decomposition is performed using PSI measurements obtained from descending and ascending acquisition geometries. On the other hand, because the SAR satellites move in orbits close to the poles (parallel to the north-south direction) and illuminate the world with small incidence angles, the LOS direction deformation is less sensitive to north-south movements [67][68][69]. Therefore, in this study, the north-south movements were ignored, and displacements in the east-west directions were estimated using Equations (4) and (5), as shown in Figure 8.
A = cosθ asc −cosα asc sinθ asc cosθ dsc −cosα dsc sinθ dsc (5) where d LOS , d ver , d hor , α, and θ refer to the displacements along LOS direction, vertical displacement, horizontal displacement, incidence angle, and azimuth angle, respectively [70]. Results in the LOS direction from both atmospheric correction approaches show a similar deformation pattern with a caldera-centered −10 mm/yr subsidence. However, the linear phase-based atmospheric correction approach exhibits lower standard deviation values in the areas with maximum subsidence. Therefore, we use the linear phase-based  Time series show the same findings in more detail than the mean deformation velo ity maps. Figure 7 shows the deformation time series of four certain PS pixels in ascendin and descending orbits at Tendürek volcano. The results show that the Tendürek volcan subsidence gradually evolves between 2018 and 2022, and the temporal evolution of th displacement is linear. Deformation measurements made in areas where the topography changes suddenly such as Tendürek volcano, may be adversely affected by atmospheric effects, slope, and satellite LOS directions, and the pixel distributions in the resulting datasets may not b homogeneous. Although two different atmospheric correction approaches were applied in this study, atmospheric artifacts were completely removed, and therefore, sufficien  The components of the deformation show that the vertical deformation rate is between 5 and −15 mm/yr, and the horizontal deformation rate is between 5 and −5 mm/yr (Figure 8a,b). The vertical deformation component shows the subsidence of a peak caldera-centered region with a maximum of −15 mm/yr. In addition, both the northern and southern parts of the deformed region show an eastward displacement. Regarding deformations in the east-west direction (Figure 8b), the rate of displacement is much smaller, and the maximum rate is approximately −5 to 5 mm/yr. This shows that the deformation in the Tendürek volcano predominantly occurred in the vertical direction, and the horizontal deformation component is negligible.

Discussion
In this study, we aimed to determine the surface deformation occurring in the Tendürek volcano. A total of 150 Sentinel-1 images from descending and ascending orbits were used, and two different atmospheric correction approaches were applied to estimate the atmospheric phase contribution in the study area. Although both atmospheric correction approaches show a similar deformation pattern, the deformation rates obtained from the linear phase-based approach were used in the continuation of the study because the linear phase-based approach exhibited lower standard deviation values. In this context, the results obtained from both satellite geometries show that there was a maximum subsidence of −50 mm over 5 years at the Tendürek volcano. The standard deviations of the deformation rates are ∼5 mm, which is approximately 10 times less than the deformation rate in both cases.
Profile analyses of both geometries show a bowl-shaped pattern that continues from the summit of the volcano to the eastern flanks and demonstrates that the maximum subsidence rate decreases linearly with the distance from the summit of the volcano. In addition, when the deformation time series of four PS pixels selected from the ascending and descending orbits are examined, it is observed that the subsidence occurs in a linear structure and there is no increase or decrease in the deformation rate in any period. As a result, the subsidence occurred in a summit caldera-centered region of the Tendürek volcano. However, to investigate whether the horizontal movement affects the obtained deformation rates, the deformation components were calculated, and the horizontal deformation was small enough to be neglected. When the vertical deformation components were examined, a deformation of −15 mm/yr was detected in the Tendürek volcano.
As mentioned above, ref. [40] used ERS satellite images and detected a 10-20 cm sub- The components of the deformation show that the vertical deformation rate is between 5 and −15 mm/yr, and the horizontal deformation rate is between 5 and −5 mm/yr (Figure 8a,b). The vertical deformation component shows the subsidence of a peak calderacentered region with a maximum of −15 mm/yr. In addition, both the northern and southern parts of the deformed region show an eastward displacement. Regarding deformations in the east-west direction (Figure 8b), the rate of displacement is much smaller, and the maximum rate is approximately −5 to 5 mm/yr. This shows that the deformation in the Tendürek volcano predominantly occurred in the vertical direction, and the horizontal deformation component is negligible.

Discussion
In this study, we aimed to determine the surface deformation occurring in the Tendürek volcano. A total of 150 Sentinel-1 images from descending and ascending orbits were used, and two different atmospheric correction approaches were applied to estimate the atmospheric phase contribution in the study area. Although both atmospheric correction approaches show a similar deformation pattern, the deformation rates obtained from the linear phase-based approach were used in the continuation of the study because the linear phase-based approach exhibited lower standard deviation values. In this context, the results obtained from both satellite geometries show that there was a maximum subsidence of −50 mm over 5 years at the Tendürek volcano. The standard deviations of the deformation rates are ∼5 mm, which is approximately 10 times less than the deformation rate in both cases.
Profile analyses of both geometries show a bowl-shaped pattern that continues from the summit of the volcano to the eastern flanks and demonstrates that the maximum subsidence rate decreases linearly with the distance from the summit of the volcano. In addition, when the deformation time series of four PS pixels selected from the ascending and descending orbits are examined, it is observed that the subsidence occurs in a linear structure and there is no increase or decrease in the deformation rate in any period. As a result, the subsidence occurred in a summit caldera-centered region of the Tendürek volcano. However, to investigate whether the horizontal movement affects the obtained deformation rates, the deformation components were calculated, and the horizontal deformation was small enough to be neglected. When the vertical deformation components were examined, a deformation of −15 mm/yr was detected in the Tendürek volcano.
As mentioned above, ref. [40] used ERS satellite images and detected a 10-20 cm subsidence at the Tendürek volcano summit between 1993 and 2000, which is consistent with our study considering that 15-20 mm/yr corresponds to a deformation of 12-16 cm over 8 years. In the study performed by [41] using Envisat satellite images, a subsidence of 11 mm/yr was detected in an area of 10 × 14 km on the Tendürek volcano between 2003 and 2010, and the cause of this subsidence was attributed to a magmatic threshold-like source located 4.5 km below the summit. Another study conducted in the region was carried out by [42] using Sentinel-1 satellite images of 10 volcanic systems in Turkey, including the Tendürek volcano, between 2015 and 2020, and no significant deformation rate was detected in the Tendürek volcano during this period. Their observation is not consistent with the results of our study. In the study conducted by [42], a reliable interpretation cannot be made because a detailed deformation map of the Tendürek volcano or time series of deformation were not given. In addition, the fact that the two studies cover different periods, using different methods and images, may explain this difference. Overall, we obtained similar results to those obtained by [40,41], and we conclude that the deformation in the Tendürek volcano has been ongoing since 1993.

Conclusions
In conclusion, the surface deformation of Tendürek volcano was determined using 150 Sentinel-1 SAR images captured between May 2018 and October 2022 using the SNAP-StaMPS integration. In this way, the most up-to-date and comprehensive analysis of the Tendürek volcano was carried out. Since the study area is a mountainous region with a high tendency to be affected by tropospheric delay, we applied two different atmospheric correction approaches. However, in our case, the GACOS atmospheric correction method produced larger standard deviation values for the maximum deformation rate. Therefore, we used the linear phase-based atmospheric correction approach. As a result, 10.6 ± 0.76 and 10 ± 1.18 mm/yr LOS velocities were determined from descending and ascending acquisition geometries, respectively, and updated deformation maps showing annual mean velocity changes for the region were created.
This study confirms that the StaMPS-InSAR analysis on Sentinel-1 datasets has great potential for investigating surface deformations caused by volcanoes. In this context, the findings of this study may have important implications for authorities in disaster management, owing to the ability of this technique to achieve rapid results in large areas. However, this work can be improved in the future. For instance, surface deformations can be detected with higher sensitivity using a local atmospheric model and high-resolution InSAR images. It may also be useful to compare InSAR observations with ground measurements if ground-based monitoring instruments can be established in the future.  Acknowledgments: The Sentinel-1A and 1B data used in the present study were accessed through the Sentinel Scientific Data Hub (scihub.copernicus.eu, accessed on 6 January 2023) of the Copernicus Open Access Hub.

Conflicts of Interest:
The authors declare no conflict of interest.