Survey of land subsidence – case study: The land subsidence formation in artificial recharge ponds at South Hamadan Power Plant,northwest of Iran

Journal of Earth System Science - The artificial recharge is a technique of aquifer conservation for land subsidence. But in this article, the phenomenon of land subsidence and the resulting cracks...     

Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)

The North China Plain (NCP) has been suffering from groundwater storage (GWS) depletion and land subsidence for a long period. This paper collects data on GWS changes and land subsidence from in situ groundwater-level measurements, literature, and satellite observations to provide an overview of the evolution of the aquifer system during 1971–2015 with a focus on the sub-regional variations. It is found that the GWS showed a prolonged declining rate of ?17.8?±?0.1 mm/yr during 1971–2015, with a negative correlation to groundwater abstraction before year ~2000 and a positive correlation after ~2000. Statistical correlations between subsidence rate and the GWS anomaly (GWSA), groundwater abstraction, and annual precipitation show that the land subsidence in three sub-regions (Beijing, Tianjin, and Hebei) represents different temporal variations due to varying driver factors. Continuous drought caused intensive GWS depletion (?76.1?±?6.5 mm/yr) and land subsidence in Beijing during 1999–2012. Negative correlations between total groundwater abstraction and land subsidence exhibited after the 1980s indicate that it may be questionable to infer subsidence from regional abstraction data. Instead, the GWSA generally provides a reliable correlation with subsidence. This study highlights the spatio-temporal variabilities of GWS depletion and land subsidence in the NCP under natural and anthropogenic impacts, and the importance of GWS changes for understanding land subsidence development.     

The Lagoon of Venice： geological setting,evolution and land subsidence

The paper deals with the geological setting, history and subsidence of the Venetian Plain. Major attention is paid to the Pleistocene-Holocene stratigraphic sequence in the Lagoon of Venice, in relation to its origin that datesback to 6-7 kyr BP. Geological land subsidence, which played an important role in the origin and the evolution of the lagoon, and anthropogenic subsidence, that has recently assumed a major importance for the Venetian environment, are discussed. Considering also the sealevel rise, 23 cm loss in land elevation has occurred in the last century, leading to increased flooding events and environmental problems that require protective works.     

Review: Water–rock interactions and related eco-environmental effects in typical land subsidence zones of China

Land subsidence is common in some regions of China. Various eco-environmental problems have arisen due to changes in water–rock interactions in these subsided areas, for which a comprehensive understanding of the hydrogeological setting is needed. This paper presents the general status of land subsidence in three typical subsided areas of China through the compilation of relevant data, and reviews some typical changes in the water–rock interactions in subsided areas along with related eco-environmental issues. It is found that the subsidence development and distribution are controlled by the groundwater-withdrawal intensity externally, and by the thickness and compressibility of unconsolidated sediments internally. The physical changes and related effects of water–rock interactions in subsided areas include: (1) the decreased ground elevation that caused floods, waterlogged farmland, etc.; (2) the differential subsidence that caused ground fissures; and (3) the change of seepage field that caused substantial reduction of the water resource. Chemically, the changes and related effects of water–rock interactions include: (1) the change to the chemical environment or processes due to the hydrogeologic structure alteration, which caused groundwater pollution; and (2) hydrologic mixing (seawater intrusion, artificial recharge; exchange with adjacent aquifers or aquitards), which degraded the groundwater quality. Further research on the subsided areas in China is suggested to reveal the mechanisms regarding biological and gaseous (meteorological) changes from the perspective of interacting systems among water, rocks, biological agents and gases.     

Helium isotopic fingerprints of the heavy land subsidence left in groundwater of the Saga Plain,near Beppu–Shimabara graben,Kyushu, Japan

The Saga Plain is near Beppu–Shimabara graben, a region of potential active volcanism. In the graben, mantle He, which has a high ³He/⁴He ratio of 1.1 × 10⁻⁵, escapes easily from the underlying subduction zone. In groundwater of the Saga Plain, except in the Shiroishi district, this ratio gradually increased as the dissolved He content increased, to a maximum of 5 × 10⁻⁶. In central Shiroishi, however, the ratio reached a minimum of 8.7 × 10⁻⁷ with increasing dissolved He content, suggesting that groundwater in central Shiroishi has selectively accumulated radiogenic He, which has a very low ratio of 1 × 10⁻⁸, rather than reflecting the regional He, which is rich in mantle He. This can be explained if groundwater in Shiroishi has become mixed with fossil pore water drawn from impermeable marine clay aquitard layers. The withdrawal of pore water has also caused severe land subsidence in central Shiroishi.     

Monitoring of land subsidence and ground fissures in Xian,China 2005–2006: mapped by SAR interferometry

The City of Xian, China, has been experiencing significant land subsidence and ground fissure activities since 1960s, which have brought various severe geohazards including damages to buildings, bridges and other facilities. Monitoring of land subsidence and ground fissure activities can provide useful information for assessing the extent of, and mitigating such geohazards. In order to achieve robust Synthetic Aperture Radar Interferometry (InSAR) results, six interferometric pairs of Envisat ASAR data covering 2005–2006 are collected to analyze the InSAR processing errors firstly, such as temporal and spatial decorrelation error, external DEM error, atmospheric error and unwrapping error. Then the annual subsidence rate during 2005–2006 is calculated by weighted averaging two pairs of D-InSAR results with similar time spanning. Lastly, GPS measurements are applied to calibrate the InSAR results and centimeter precision is achieved. As for the ground fissure monitoring, five InSAR cross-sections are designed to demonstrate the relative subsidence difference across ground fissures. In conclusion, the final InSAR subsidence map during 2005–2006 shows four large subsidence zones in Xian hi-tech zones in western, eastern and southern suburbs of Xian City, among which two subsidence cones are newly detected and two ground fissures are deduced to be extended westward in Yuhuazhai subsidence cone. This study shows that the land subsidence and ground fissures are highly correlated spatially and temporally and both are correlated with hi-tech zone construction in Xian during the year of 2005–2006.     

Hydrogeochemistry and quality of groundwater of coastal unconfined aquifer in Amol–Ghaemshahr plain,Mazandaran Province,Northern Iran

Groundwater of the unconfined aquifer (1,100 sq. km) of a two-tier coastal aquifer located in the Amol–Ghaemshahr plain, Mazandaran Province, Northern Iran, is classified into fresh and brackish water types. Fresh groundwater (FGW) samples (n = 36) are characterized by Ca²⁺ > Na⁺ > Mg²⁺ > K⁺ and HCO₃ ^? > Cl^? > SO₄ ^2? > NO₃ ^?. Spearman’s rank correlation coefficient matrices, factor analysis data, values of the C-ratio (av. = 0.89) and CAI and values of the molar ratios of Ca²⁺/HCO₃ ^?, Ca²⁺/SO₄ ^2?, Mg²⁺/HCO₃ ^? and Mg²⁺/SO₄ ^2? indicate that the ionic load in the FGW is derived essentially from carbonic acid-aided weathering of carbonates and aluminosilicates, saline/sea water trapped in the aquifer sediments (now admixed with the groundwater) and ion exchange reactions. Values of the CAI and Na⁺/Cl^? molar ratio suggest that the part of the Ca²⁺ (±Mg²⁺) content in 23 FGW samples is derived from clay minerals of the aquifer matrix, and part of the Na⁺ content in 20, 12, and 3 FGW samples is derived, respectively, from alkali feldspar weathering, clay minerals of the aquifer matrix and rain water and/or halite. Brackish groundwater (BGW) samples (n = 4) contain Cl^? as the dominant anion and their average total ionic concentration (38.65 meq/L) is 1.79 times higher than that of the FGW samples (21.50 meq/L). BGW pockets were generated by non-conservative mixing of FGW with the upconed saline water from the underlying saline groundwater zone of the semi-confined aquifer along bore wells involved in excessive extraction of groundwater from the unconfined aquifer. Groundwater belongs essentially to “high salinity, low sodium” irrigation water class.     

Application of “panel-data” modeling to predict groundwater levels in the Neishaboor Plain, Iran

The aim of this research was to predict groundwater levels in the Neishaboor plain, Iran, using a ??panel-data?? model. Panel-data analysis endows regression analysis with both spatial and temporal dimensions. The spatial dimension pertains to a set of cross-sectional units of observation. The temporal dimension pertains to periodic observations of a set of variables characterizing these cross-sectional units over a particular time span. Firstly, the available observation wells in the Neishaboor plain were clustered according to their fluctuation behavior using the ??Ward?? method, which resulted in six areal zones. Then, for each cluster, an observation well was selected as its representative, and for each zone, values of monthly precipitation and temperature, as independent variables, were estimated by the inverse-distance method. Finally, the performance of different panel-data regression models such as fixed-effects and random-effects models were investigated. The results showed that the two-way fixed-effects model was superior. The performance indicators for this model (R ²?=?0.97, RMSE?=?0.05?m and ME?=?0.81?m) reveal the effectiveness of the method. In addition, the results were compared with the results of an artificial-neural-network (ANN) model, which demonstrated the superiority of the panel-data model over the ANN model.     

Energetic and spatial characterization of seismicity in the Algeria–Morocco region

We estimate the energetic and spatial characteristics of seismicity in the Algeria–Morocco region using a variety of seismic and statistical parameters, as a first step in a detailed investigation of regional seismic hazard. We divide the region into five seismotectonic regions, comprising the most important tectonic domains in the studied area: the Moroccan Meseta, the Rif, the Tell, the High Plateau, and the Atlas. Characteristic seismic hazard parameters, including the Gutenberg–Richter b-value, mean seismic activity rate, and maximum possible earthquake magnitude, were computed using an extension of the Aki–Utsu procedure for incomplete earthquake catalogs for each domain, based on recent earthquake catalogs compiled for northern Morocco and northern Algeria. Gutenberg–Richter b-values for each zone were initially estimated using the approach of Weichert (Bull Seismol Soc Am 70:1337–1346, 1980): the estimated b-values are 1.04 ± 0.04, 0.93 ± 0.10, 0.72 ± 0.03, 0.87 ± 0.02, and 0.77 ± 0.02 for the Atlas, Meseta, High Plateau, Rif, and Tell seismogenic zones, respectively. The fractal dimension D ₂ was also estimated for each zone. From the ratio D ₂/b, it appears that the Tell and Rif zones, with ratios of 2.09 and 2.12, respectively, have the highest potential earthquake hazard in the region. The Gutenberg–Richter relationship analysis allows us to derive that in the Tell and Rif, the number of earthquake with magnitude above Mw 4.0, since 1925 normalized to decade and to square cell with 100-km sides is equal to 2.6 and 1.91, respectively. This study provides the first detailed information about the potential seismicity of these large domains, including maximum regional magnitudes, characteristics of spatial clustering, and distribution of seismic energy release.

     

Energetic and spatial characterization of seismicity in the Algeria–Morocco region

We estimate the energetic and spatial characteristics of seismicity in the Algeria–Morocco region using a variety of seismic and statistical parameters, as a first step in a detailed investigation of regional seismic hazard. We divide the region into five seismotectonic regions, comprising the most important tectonic domains in the studied area: the Moroccan Meseta, the Rif, the Tell, the High Plateau, and the Atlas. Characteristic seismic hazard parameters, including the Gutenberg–Richter b-value, mean seismic activity rate, and maximum possible earthquake magnitude, were computed using an extension of the Aki–Utsu procedure for incomplete earthquake catalogs for each domain, based on recent earthquake catalogs compiled for northern Morocco and northern Algeria. Gutenberg–Richter b-values for each zone were initially estimated using the approach of Weichert (Bull Seismol Soc Am 70:1337–1346, 1980): the estimated b-values are 1.04 ± 0.04, 0.93 ± 0.10, 0.72 ± 0.03, 0.87 ± 0.02, and 0.77 ± 0.02 for the Atlas, Meseta, High Plateau, Rif, and Tell seismogenic zones, respectively. The fractal dimension D ₂ was also estimated for each zone. From the ratio D ₂/b, it appears that the Tell and Rif zones, with ratios of 2.09 and 2.12, respectively, have the highest potential earthquake hazard in the region. The Gutenberg–Richter relationship analysis allows us to derive that in the Tell and Rif, the number of earthquake with magnitude above Mw 4.0, since 1925 normalized to decade and to square cell with 100-km sides is equal to 2.6 and 1.91, respectively. This study provides the first detailed information about the potential seismicity of these large domains, including maximum regional magnitudes, characteristics of spatial clustering, and distribution of seismic energy release.     

Spatiotemporal distribution of flood disasters in Asia and influencing factors in 1980–2019

Natural Hazards - Floods are common natural disasters in Asia. Flood datasets from 48 countries in Asia were collected to investigate the spatiotemporal distribution and influencing factors, using...     

Delineation of groundwater potential zones using remote sensing,GIS, and AHP technique in Tehran–Karaj plain,Iran

Evaluation of groundwater resources in dry areas without enough data is a challenging task in many parts of the world, including Tehran–Karaj plain in Iran, which includes Tehran, the capital city of Iran and Karaj, one of Iran’s biggest cities. Water demand due to increasing agricultural and industrial activities caused many problems in the field of water resources management. In this study, the potential of groundwater resources was evaluated using remote sensing, geographic information system (GIS), and analytic hierarchy process (AHP) for the first time. Digital Elevation Model from Shuttle Radar Topography Mission was used to generate a slope map and drainage density map. Three Landsat-8 satellite images were utilized to provide lineament density and land cover/land use maps. Geological and soil type maps were provided from the Geological Survey and Mineral Explorations of Iran (GSI). Tropical Rainfall Measuring Mission data were used to prepare average annual precipitation map. Discharge values from 102 pumping wells in the time period of 2002–2014 were used to evaluate the results. Seven data layers were prepared, and the geodatabase was made in GIS. The layers and their classes were assigned weights using AHP method. Finally, the layers were overlaid based on their weights, and the potential map of groundwater resources was generated. The area was classified into five zones with very high, high, moderate, low, and very low potentials. The zones covered 5.95, 32.90, 22.70, 10.20, and 28.25% of the study area, respectively. The results showed good agreement with the field data obtained from discharge wells.     

Arsenic and humic substances of groundwater in the southwestern coastal plain of Taiwan

High concentrations of arsenic and humic substances in groundwater from the southwestern coastal plain of Taiwan were well known for their probable relationships with black-foot disease. In order to realize the relationships between the concentrations of …     

Characterization of arsenic-contaminated groundwater and human health risk assessment in Vietnam

The characteristics of arsenic-contaminated groundwater were investigated and the potential risks from groundwater were evaluated in Vietnam. Arsenic contamination in groundwater was found in Ha Nam Province in the northern part of Vietnam. Since groundwater has been used as one of the main drinking water sources in this region, groundwater （n=10） and hair （n=15） samples were collected in the Vinh Tru district, Ha Nam Province, during February of 2006. At the site, the concentrations of anions and silica were analyzed and As（Ⅲ） and As（Ⅴ） species in groundwater were separated with disposable arsenic speciation cartridge. The concentrations of arsenic in groundwater ranged from 13 to 582 μg/L （mean=366 μg/L）. According to the results of arsenic speciation, approximately 90% of arsenic in groundwater existed as As（Ⅲ） species. Average concentrations of iron and manganese were 18 and 0.7 mg/L, respectively. All samples exceeded the Vietnamese drinking water standard of 10μg/L for arsenic and 0.5 mg/L for iron. Also 70% of the groundwater samples were above the Vietnamese drinking water standard of 0.5 mg/L for manganese. However, a real intake of arsenic to human is less than analyzed arsenic concentration because sand filtration units are used to remove iron and arsenic in groundwater in this region. After treatment, arsenic concentrations decreased from 〈1 to 82 μg/L （mean=33 μg/L）. The concentrations of arsenic in hair samples and treated groundwater were compared. Especially, arsenic concentrations in female hair samples and treated groundwater had significantly positive correlation with the Spearman correlation of 0.88 and the P-value of 0.001. Based on the arsenic concentrations in treated groundwater, the assessment of human health risks was conducted.     

Constraining the timing of porphyry mineralization in northwest Iran in relation to Lesser Caucasus and Central Iran; Re–Os age data for Sungun porphyry Cu–Mo deposit

The Neo-Tethyan subduction in Iran is characterized by the Urumieh–Dokhtar magmatic arc (UDMA), formed by northeast-ward subduction of the oceanic crust beneath the central Iran. This belt coincides with the porphyry copper metallogenic belt that comprises several metallogenic zones, including Ahar–Jolfa in northwest Iran. The Ahar–Jolfa metallogenic zone encompasses two main batholiths of Qaradagh and Sheyvardagh and numerous intrusive bodies of Cenozoic, which have produced many base and precious metal deposits and prospects. The former is considered as continuation of the Meghri–Ordubad pluton in South Armenian Block (SAB), which also hosts porphyry copper deposits (PCDs). The Sungun PCD is the largest occurrence in northwest Iran. Rhenium-Osmium ages of Sungun molybdenites are early Miocene and range between 22.9 ± 0.2 and 21.7 ± 0.2 Ma. Comparison of the ages obtained here with published ages for mineralization across the region suggests the following sequence. The earliest porphyry Cu–Mo mineralization event in northwest Iran is represented by Saheb Divan PCD of late Eocene age, which is followed by the second epoch of middle Oligocene, including the Cu–Mo–Au mineralization at Qarachilar and the Haftcheshmeh PCD. Mineralization in Sungun, Masjed Daghi, Kighal and Niaz deposits corresponds to the third mineralization event in northwest Iran. The first epoch in northwest Iran postdates all Eocene mineralizations in SAB, while the second epoch is coeval with Paragachay and the first-stage of Kadjaran PCDs. Its third epoch is younger than all mineralizations in SAB, except the second stage in Kadjaran PCD. Finally, the Cu mineralization epochs in northwest Iran are older than nearly all PCDs and prospects in Central Iran (except the Bondar Hanza PCD), altogether revealing an old to young trend along the UDMA and the porphyry Cu belt towards southeast, resulted from diachronous, later closure of the Neo-Tethyan oceanic basin in central and SE Iran.     

Measurement of pesticides in surface water and sediment of Caspian Sea, Iran

Background: With development of intensive agriculture in the agricultural regions of the Northern provinces (Caspian Sea coasts), which involves frequent pesticide use in highly mechanized cash-crop cultures. To investigate and provide information on pest…