Based on the data from piezometers, well logs, geophysical surveys and the interferometric synthetic aperture radar (InSAR) technique, this study investigates the main causes of land subsidence in Damaneh City, Iran. The size, openings and direction of fissures were measured by micrometer and compass. The locations of fissures and wall cracks were determined by GPS. The geoelectrical data were used to determine the composition, thickness, depth and shape of lower parts of the aquifer. Groundwater fluctuations were evaluated by available piezometers. The InSAR technique was used to measure land deformation from space and to map the dense changes of surface displacements. The results indicate that the main cause of ground subsidence is the decline of groundwater heads and changes in composition and thickness of compressible lacustrine sediments. The subsidence map obtained from the radar data of ASAR sensor of ENVISAT satellite shows that the subsidence zone is mainly in northern city that is underlain by very thick fine sediments. The subsidence rates from March to December 2005 and from July 2011 to January 2012 are 6.7 and 7.0 cm/a, respectively. The results also show good correlations among the formation of earth features, decline in groundwater head and thickness of fine-grained sediments. We recommend that the groundwater withdrawal for agricultural and industrial sectors should be restricted and urban expansion in the northern part of the city should be constrained.
Land subsidence is either a gradual settling or a sudden sinking of earth surface accompanied with some horizontal displacements (Bates and Jackson, 1980). The natural occurrence of this phenomenon is usually resulted from such geological processes as solubility, thawing of ice, vibration and densification of deposits. A major man-made occurrence of subsidence in the world is the withdrawal of fluids from subsurface (Bell, 1999; Adiyaman, 2012). Ground subsidence unfortunately can result in irreparable damages to water resources, infrastructures, agricultures, industries, and residences. Ground subsidence can also cause changes in the direction and velocity of groundwater flow (Holzer and Galloway, 2005).
According to Galloway et al. (1999), groundwater withdrawal is the most common cause of subsidence. Groundwater withdrawal from aquifers for a long period of time can decrease pore water pressure and increase effective stress (Burbey, 2002). As effective stress increases, permeability and thickness of sediments decrease, and density and weight of sediments increase, leading to permanent changes in soil grain configuration and volume. This process continues even after groundwater withdrawal stops (Adiyaman, 2012).
Subsidence phenomenon is reported to have affected many major cities such as Mexico City of Mexico (Ortega-Guerrero et al., 1999), Taipei of China (Hou et al., 2005), Bangkok of Thailand (Phien-Wej et al., 2006), Osaka and Tokyo of Japan, Arizona of USA (Wei, 2006), Venice of Italy, Jakarta of Indonesia, Kolkata of India (Chatterjee et al., 2006), Kerman of Iran (Sharifikia, 2009; Esmaili and Motagh, 2009), Shanghai and Beijing of China (Wang et al., 2009) and Kashan of Iran (Ghazifard et al., 2016). The ground subsidence phenomenon in Iran is mostly due to excessive groundwater withdrawal, high sediment thickness, low bearing capacity of sediments and solution of subsurface formations (Ghazifard et al., 2016).
It is highly desired to create a suitable system to monitor the subsidence phenomenon. Methods such as GPS (Sneed and Galloway, 2000; Bell et al., 2002), leveling (Bell and Price, 1991; Bell et al., 2002), optical remote sensing, air photo, extensometer and boreholes (Pope and Burbay, 2004; Carruth et al., 2007) have been used to measure the amount of deformation and displacement of ground surface. Advanced GPS and radar detectors are usually used to infer the information associated with ground subsidence (Galloway et al., 2000). In recent years, synthetic aperture radar interferometry has been used to measure vertical displacement with time (Hoffmann et al., 2001; Bell et al., 2002; Rosen et al., 2004; Galloway and Burbey, 2011). In this method, a series of microwave pulses are sent from satellite toward the earth surface and the satellite receives the reflected waves. The obtained interferogram of phase differences can be used to estimate the vertical displacements of earth surface (Massonnet and Feigl, 1998; Galloway and Hoffmann, 2007). This radar technique can be used in any different climate conditions and verified by GPS and leveling data.
This study deals with the land subsidence in the east and northeast parts of the Damaneh City of Iran where excessive groundwater withdrawal and declines of groundwater head have been undergone for years due to increasing population and intensifying agricultural activities, causing ground subsidence. The subsidence-generated earth fissures and rupture of houses and roads and uplifting of well casings were observed for the first time in 2013.
Damaneh City (33° 00′ 20′ ′ -33° 01′ 32′ ′ N, 50° 28′ 43′ ′ -50° 29′ 58′ ′ E) is located in a semi-arid area on the west side of Isfahan Province in Central Iran. The city covers an area of 2 km2 and is surrounded by Gahiz swamp, Darehbid and Ashjerd farmlands (Fig. 1). The altitude of Damaneh City is between 2330-2370 m a.s.l. and the average slope is about 5% toward north. Gahiz swamp, on the east side of the city, was filled with water about 50 years ago and has been used as agricultural land in recent decades.
Most of the study area is underlain by the Quaternary sediments that is underlain by the Jurassic shale. In the north part of the city, the Quaternary sediments can be over 100 m thick while the sediments in the south part is less than 10 m.
Damaneh aquifer is a part of a larger unconfined aquifer with an area of 220 km2. It is composed of Quaternary alluvium which was transported from the southern and northern highlands and deposited on this vast plain (Fig. 2). The aquifer storage coefficient is about 2.8%. The mean annual precipitation is about 343.5 mm and annual mean temperature is about 9.5° C. The annual water withdrawal through the 800 wells is about 200× 106 m3. About 97% of water withdrawal is used for agricultural irrigations.
According to the data from 16 piezometers in the study area, the average decline of groundwater head in recent years is over 1 m/a and the decline can reach over 3 m/a in dry years. The total decline of groundwater head is about 30 m from 1984 to 2012.
For this research, data of location, trends and opening of fissures were obtained by field investigations. In order to determine the composition of the aquifer materials, 5 boreholes were drilled (40 m deep) within the city limit. Data of boreholes were obtained from the Department of Roads and Urban Development of Iran. The geophysical method was used to determine the composition, thickness, depth and shape of lower parts of the aquifer and bedrock topography. Data about fluctuations of groundwater table with 16 available piezometers during a 30-year period were obtained from Isfahan Regional Water Company. The air photo was used to compare land use changes during the period 1955-2015). The radar interferometry technique was used to determine the rate and zone of subsidence. In this study, 8 ASAR images of 1.5 levels (unprocessed data) were used. The radar data from 457 and 250 tracks were evaluated. The ASAR-Envisat data were provided by the European Survey Agency (ESA) through CAT-1 C1P.10465.
The sizes and openings of fissures were measured by micrometer, the directions were measured by compass and the locations of about 45 fissures and wall cracks were determined by GPS. The borehole exploration permitted by the Department of Roads and Urban Development of Iran was to determine the type and thickness of sediments and the geoelectrical resistivity test was used to prepare the contour maps of sediment thicknesses and bedrock topography. With geoelectrical method, different sediment layers can be distinguished according to the specific electrical resistivity.
Piezometers data were spatially interpolated in GIS environment to study water table fluctuations. InSAR can derive information from the interferograms that are formed by phase differences between two high resolution SAR images for the same area (Saracin et al., 2014). The differential interferometric synthetic aperture radar (D-InSAR) is formed and the effect of topographic phase is calculated and removed (Raucoules et al., 2007; Liu et al., 2014) using digital elevation model and transforming elevation (motion) into phase. Radar data were processed by SARSCAPE in ENVI environment. In this method, interferometry with short time and spatial baseline were used (Lubis et al., 2011). The ASTER DEM (ASTER, Advanced Space borne Thermal Emission and Reflection Radiometer; DEM: Digital Elevation Model) with 30 m spatial resolution was used to remove topography phase. Goldstein filter was also used to remove noises and to strengthen radar signal. Minimum Cost Flow (MCF) was used in unwrapping phase. Interferograms were prepared from five radar data of 457 tracks (three interferograms) and three radar data of 250 tracks (three interferograms).
The locations of fissures, cracks and subsidence zones were designated and correlated with thickness of different sediment layers, bedrock topography and decline in groundwater head.
The results show that the locations of all fissures and cracks are in the northern city with east-west trends (Fig. 3a). Although the size of fissures opening varies, the biggest openings are observed in the northern and northeastern parts of the city at the positions with the lowest elevations (Fig. 3b).
Except uplifted well casings, no fissure was detected at eastern Damaneh City. But damages to houses, power lines, gas pipes and water mains are rather visible in northern city (Fig. 4).
For comparing the trend of land use changes, four air photos taken in 1955 were obtained and geo-referenced. The photos (Fig. 5a) were compared with Google Earth pictures taken in 2015 (Fig. 5b) using ArcGIS. As shown in Figure 5a, the residential areas and agricultural lands were very limited and the dark tone of Gahiz swamp shows that water table was rather high in the area. But, in Figure 5b, extensively developed residential areas and agricultural lands are radially observable. The city area expanded from 0.31 to 2.00 km2 during the period 1955-2015.
Figure 6a shows that Damaneh City sits on highly compressible lacustrine sediments in 2015 and Figure 6b shows that the situation of fissures in northern city is near the boundary of two different types of sediments in 1955.
In this study, the declines of groundwater head from March to December 2005 (Fig. 7a) and from January to July 2012 (Fig. 7b) were interpolated based on the data from 16 piezometers in the city. The highest declines of groundwater head in 2005 and 2012 were 2.13 and 7.14 m (Figs. 7a and b), respectively. This decline increases toward the western city where the area was dominated by coarser alluvium deposits. The sharp decline in groundwater head from 2005 to 2012 is mostly due to the increase in water withdrawal for agricultural irrigation.
Subsidence features were observed in 5 boreholes drilled in the northern city (Fig. 8). Boreholes 1, 4 and 5 were 30 m deep. Borehole 2 and 3 were 27 and 40 m deep, respectively. As shown in Figure 8, the sediments of boreholes 1, 2 and 3 were finer than those in boreholes 4 and 5.
According to the prepared profiles of the subsidence zones, the biggest subsidence occurred in the positions where finer sediments are present (boreholes 1, 2 and 3).
The geoelectrical resistivity tests were conducted at 57 points along 8 profiles. The positions of the geophysical survey are shown in Figure 9. Two NW-SE profiles, one was from V1 to DC2 and and another one was from V1 to V14, were prepared, types, thicknesses and stratigraphy of the sediments along with the log data of these two profiles are shown in Figures 10a and b. Most of fissures are observed on the thickest layers with fine and compressible sediments.
The results of the interpolated geophysical data along with the changes of sediments thickness and the depth of bedrocks by ArcGIS are shown in Figures 11a and b.
In this study, interferometry technique was also used to determine the zone and rate of subsidence. In this regard the band C of radar data from ASAR sensor of ENVISAT satellite related to 457 and 250 tracks were used. Table 1 shows the radar data used for differential interferometry and six interferograms used in this study and the parameters of created interferograms.
The subsidence zone and the rate of subsidence for different period of time are shown in Figure 12. The subsidence zone is merely detectable in the north part of the city. The biggest subsidence rates were 6.7 and 7 cm/a (Fig. 12f) for the period spanning from 25 March to 30 December, 2005 (Fig. 12c) and from 12 July, 2011 to 8 January, 2012, respectively. These two figures show that the locations of fissures are within the subsidence zone.
The subsidence map obtained from interferogram studies and contour map of decline in groundwater head obtained from 16 piezometers were presented in Figure 13. The radar pictures from 25 March to 30 December, 2005 (Fig. 13a) and 12 July, 2011 to 8 January, 2012 (Fig. 13b) are well correspondent with the contour maps of decline in groundwater head.
A subsidence map from July 2011 to January 2012 was prepared and correlated with the contour maps of sediment thickness (Fig. 14a) and bedrock topography (Fig. 14b). Figure 14a shows that the thickness of fine sediments increased toward north where fissures obtained from radar processing and field measurements are intensively present. Map of bedrock topography also shows that the depth of bedrock increases toward north that is underlain by highly compressible fine sediments.
Ground subsidence and earth fissures have caused extensive damages to houses, infrastructures and roads in Damaneh City. The fissures in north part of the city are located on the sediments that used to be a swamp with soft clayey sediments and high water table. The results of this study show that groundwater withdrawal and associated decline in groundwater head have major effects on ground subsidence. The decline in groundwater head can result in reduction of pore water pressure and cause the effective stress to be increased and soil grains become more compacted (Burbey, 2002). Bonì et al. (2015) also found that the thickness of alluvial deposit has a direct relationship with the amount of subsidence in Guadalentin Basin of Spain. According to Helm (1994) and Bouwer (1977), various thicknesses of compressible soils can cause differential settlement and fissuring of ground surface. In this regard, differential settlement causes compaction at deeper layers and tension at shallower layers. As a result of differential settlement, the sediments tilt toward the areas with higher compaction. Adiyaman (2012) reported that the amount of subsidence is controlled by sediment type, thickness, stiffness, compressibility of soil and loading history. Adiyaman explained that the subsidence even at uniform decline rate of groundwater head will not be uniform due to the changes of soil type. The geophysical studies indicated that all of the fissures are observed in the north part of the city where the thickest layers of fine and compressible sediments with higher consolidation potential are present. The geophysical data also indicated that the presence of coarse-grain sediments in the southern part of city did not allow any consolidation and subsidence occurrences.
In recent years, interferometry is widely used to measure the deformation and displacement of the earth surface (Raucoules et al., 2007; Lubis et al., 2011; Bonì et al., 2015); and it is more attractive for displacement detection due to the high resolution of C-band SAR data (Lubis et al., 2011). In this study interferometry technique (ASAR images) was used to determine the zone and rate of subsidence. Bonì et al. (2015) reported that there is a clear relation between decline in groundwater head and the obtained subsidence map from interferometry images. The results of present study show that the subsidence zone is merely within northern part of the city and the highest subsidence rate ranges from 6.7 to 7.0 cm/a.
Comparison between the map of subsidence rate obtained from interferogram studies and the rate of decline in groundwater head obtained from 16 piezometers indicated that the maximum subsidence in the north part of city were correlated with the high decline in groundwater head. The locations of fissures can be detected along the border lines with thick and fine compressible sediments by correlating the contour maps of the sediment thickness with bedrock topography. This was also verified by radar results. These findings also agree with the studies conducted by Freyre et al. (2003) and Sheng et al. (2003). According to the studies conducted by Bonì et al. (2015), after developing a soil thickness map and estimating the subsidence rate from different available satellite data, the relationship between soil thickness and subsidence rates can be established.
Based on field investigations, piezometer data, geophysical survey, borehole studies, land use changes in the period 1955-2015 and radar interferometry technique, we concluded that the subsidence zone is located in the north part of Damaneh City and the major geological factor is the presence of very thick fine and compressible sediments. The amount of subsidence were 6.7 and 7 cm for the periods of March to December 2005 and July 2011 to January 2012, respectively. It can then be concluded that this phenomenon was caused by an increase in number of deep wells in the area (about 800 wells in and around the city), intense water withdrawal and decline in groundwater head. Therefore, the main cause of subsidence is not only due to the excessive groundwater pumping for intense agricultural activities but it is also due to housing development on thick and fine compressible sediments in the north part of the city. We suggested that groundwater withdrawal from the thick fine sediments should be restricted and the urban expansion in the northern direction should be constrained to prevent further damages to the city.
This research was conducted in department of geology at the University of Isfahan. The authors would like to thank the European Space Agency (ESA) for providing SAR images (ENVISAT data).
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|