Analysis of the influence of groundwater on land subsidence in Beijing based on the geographical weighted regression (GWR) model

Geo-spatial epidemiology of gallbladder cancer in Bihar, India

Recently, a substantial increase in gallbladder cancer (GBC) cases has been reported in Bihar, India. The region's groundwater can naturally contain harmful concentrations of arsenic, which appears to be epidemiologically linked to the unusually high incidence. However, the root causes remain largely unexplored. Recent findings of uranium in the state's groundwater may also have associations. This study investigates the geo-spatial epidemiology of GBC in Bihar, India-with a focus on the correlation between environmental carcinogens, particularly arsenic and uranium in groundwater, and the incidence of GBC. Utilizing data from 8460 GBC patients' registration records over an 11-year period at a single health center, the research employs Semi-parametric Geographically Weighted Poisson Regression (S-GWPR) to account for non-stationarity associations and explores significant factors contributing to GBC prevalence at a subdistrict level. The S-GWPR model outperformed the standard Poisson regression model. The estimates suggest that arsenic and uranium concentrations in groundwater did not present significant associations; however, this could be due to the lower resolution of this data at the district level, necessitating higher resolution data for accurate estimates. Other socio-environmental factors included demonstrated significant regional heterogeneity in their association with GBC prevalence. Notably, each 1 % increase in the coverage of well- and canal-irrigated areas is associated with a maximum of 3.0 % and 5.2 % rise in the GBC incidence rate, respectively, likely attributable to carcinogen exposure from irrigation water. Moreover, distance to the health center and domestic electricity connections appear to influence the number of reported GBC cases. The latter suggests that access to electricity might have facilitated the use of groundwater pumps-increasing exposure to carcinogens. The results underscore the necessity for targeted health policies and interventions based on fine-resolution spatial analysis, as well as ongoing environmental monitoring and research to better understand the multifaceted risk factors contributing to GBC.

Unveiling the driving factors of urban land subsidence in Beijing, China

Land subsidence, an insidious and gradual geological phenomenon, presents a latent threat to future urban development and socio-economic progress. Beijing City, renowned for its high population density, has encountered significant challenges associated with land subsidence. In this study, we leverage time-series interferometric synthetic aperture radar (time-series InSAR) method to analyze the spatio-temporal patterns of land subsidence in Beijing. Furthermore, we quantify the contributions of natural and anthropogenic factors to land subsidence. Our findings reveal that land subsidence primarily occurs in the plain area of Beijing, exhibiting an average rate of -5.6 mm/year (Positive values indicate uplift, while negative values indicate subsidence.). Notably, several large-scale subsidence centers are identified, with the maximum subsidence rate reaching an alarming -232.7 mm/year. The assessments indicate that geological factors, specifically fault activity, account for 33 % of the observed land subsidence, while human activities contribute to the remaining 67 %, with groundwater overexploitation playing a prominent role. The insights gained from this study provide a foundation for understanding the causative factors behind urban land subsidence and can aid in the formulation of effective intervention policies targeting this critical issue.

Water, Water Everywhere, but Every Drop Unique: Challenges in the Science to Understand the Role of Contaminants of Emerging Concern in the Management of Drinking Water Supplies

The protection and management of water resources continues to be challenged by multiple and ongoing factors such as shifts in demographic, social, economic, and public health requirements. Physical limitations placed on access to potable supplies include natural and human-caused factors such as aquifer depletion, aging infrastructure, saltwater intrusion, floods, and drought. These factors, although varying in magnitude, spatial extent, and timing, can exacerbate the potential for contaminants of concern (CECs) to be present in sources of drinking water, infrastructure, premise plumbing and associated tap water. This monograph examines how current and emerging scientific efforts and technologies increase our understanding of the range of CECs and drinking water issues facing current and future populations. It is not intended to be read in one sitting, but is instead a starting point for scientists wanting to learn more about the issues surrounding CECs. This text discusses the topical evolution CECs over time (Section 1), improvements in measuring chemical and microbial CECs, through both analysis of concentration and toxicity (Section 2) and modeling CEC exposure and fate (Section 3), forms of treatment effective at removing chemical and microbial CECs (Section 4), and potential for human health impacts from exposure to CECs (Section 5). The paper concludes with how changes to water quantity, both scarcity and surpluses, could affect water quality (Section 6). Taken together, these sections document the past 25 years of CEC research and the regulatory response to these contaminants, the current work to identify and monitor CECs and mitigate exposure, and the challenges facing the future.

Spatial-temporal differences and influencing factors of coupling coordination between urban quality and technology innovation in the Guangdong-Hong Kong-Macao Greater Bay Area

The coordinated development of urban quality and technology innovation is an important element of China's technology innovation development strategy in the new era. Based on entropy TOPSIS, coupling coordination models, the gravity center and standard deviation ellipse method, the geographic probe, the GWR, and other methods, we explore the spatial variation and influencing factors of the coupling coordination relationship between urban quality and technology innovation in the Guangdong-Hong Kong-Macao Greater Bay Area from 2011 to 2020. It is found that: (1) the spatial distribution of the coupling coordination shows a characteristic of "high in the middle and low in the surroundings," and (2) the level of benign interaction in the central region is becoming more prominent. The center of gravity of coupling coordination moves toward the northeast, and the standard deviation ellipse shows a contraction trend away from the southwest. (3) Agglomeration capacity, human capital, cultural development, and infrastructure can significantly drive the improvement of the coupling coordination of urban quality and technology innovation, and the two-factor influence is significantly increased after the interaction. (4) The feedback effects of the coupling and coordination states of different cities on each factor have significant spatial differences and show the characteristics of hierarchical band distribution.

A Novel Groundwater Burial Depth Prediction Model Based on Two-Stage Modal Decomposition and Deep Learning

The variability of groundwater burial depths is critical to regional water management. In order to reduce the impact of high-frequency eigenmodal functions (IMF) generated by complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) on the prediction results, variational modal decomposition (VMD) is performed on the high frequency IMF components after the primary modal decomposition. A convolutional neural network-gated recurrent unit prediction model (CNN-GRU) is proposed to address the shortcomings of traditional machine learning which cannot handle correlation information and temporal correlation between time series. The CNN-GRU model can extract the implicit features of the coupling relationship between groundwater burial depth and time series and further predict the groundwater burial depth time series. By comparing the prediction results with GRU, CEEMDAN-GRU, and CEEMDAN-CNN-GRU models, we found that the CEEMDAN-VMD-CNN-GRU prediction model outperformed the other prediction models, with a prediction accuracy of 94.29%, good prediction results, and high model confidence.

Developing spatio-temporal approach to predict economic dynamics based on online news

Economic forecasting is a scientific decision-making tool, and it is one of the important basis for the government to formulate economic plans, predict the implementation of the plan, and guide the implementation of the plan. Current knowledge about the use of online news in the prediction of economic patterns in China is limited, especially considering the spatio-temporal dynamics over time. This study explored the spatio-temporal patterns of economic output values in Yinzhou, Ningbo, China between 2018 and 2021, and proposed generalized linear model (GLM) and Geographically weighted regression (GWR) model to predict the dynamics using online news data. The results indicated that there were spatio-temporal variations in the economic dynamics in the study area. The online news showed a great potential to predict economic dynamics, with better performance in the GWR model. The findings suggested online news combining with spatio-temporal approach can better forecast economic dynamics, which can be seen as a pre-requisite for developing an online news-based surveillance system The advanced spatio-temporal analysis enables governments to garner insights about the patterns of economic dynamics over time, which may enhance the ability of government to formulate economic plans and to predict the implementation of the plan. The proposed model may be extended to greater geographic area to validate such approach.

Effects of Human Activities on Urban Vegetation: Explorative Analysis of Spatial Characteristics and Potential Impact Factors

Since the 21st century, large cities around the world have experienced the transition from economically destructive development to a harmonious eco-environment. Understanding the dynamic relationships between human activities and urban eco-environment in this transition is a challenging and essential topic. The normalized difference vegetation index (NDVI) can reflect the urban vegetation cover status well. Socio-economic indexes can present the intensity and spatiality of human activities quantitatively. This work aims to use traditional regression models and machine learning algorithms to analyze the impact of socio-economic factors on NDVI accurately. Random forest regression (RFR) was performed to initially assess the contributions of all factors on NDVI, which was the numerical basis for feature selection. Subsequently, detailed dynamic relationship simulations were implemented using geographically weighted regression. In the case of Wuhan in China, the results showed that the goodness-of-fit of NDVI with socio-economic factors generally exceeded 50%. The influence coefficients changed from negative to positive, and 2010 was the turning point, indicating that human activities gradually played a favorable role in protecting vegetation during this transition period. The urban–rural interface, which was located between urban centers and marginal urban suburbs, was the area where human activities contributed most to vegetation. Thus, policy makers should focus on planning and managing housing construction and vegetation planting in urban–rural interface to relieve the population burden of the central area and improve the environmental conditions of the urban eco-environment subconsciously.

Microstructural Variation of Clay during Land Subsidence and the Correlation between Macroscopic and Microscopic Parameters

The nonlinear deformation, visco-elasto-plasticity and other macroscopic properties of soil are the concentrated manifestations of its microstructural state. In order to study the microstructural characteristics and variations of the clay under the action of additional stress caused by groundwater exploitation, borehole sampling was carried out on the clay layers at different depths in a typical land subsidence area. Consolidation tests, freeze-drying, ion sputtering, and scanning electron microscope (SEM) were conducted in order to scan and analyze the microstructure of the test samples at different scales. The Particles and Cracks Analysis System (PCAS) was used to quantify the microscopic parameters, the variations of the microstructural parameters with consolidation loads at different sizes were revealed, and the correlation between the macroscopic and microscopic parameters were discussed. The results show that: (1) the microstructural characteristics of soils with different buried depths have directivity, to a certain extent; (2) as the consolidation load increases, the average unit area and average form factor of the soil microstructure generally decrease, the structural arrangement of the unit gradually tends to be orderly, and the average pore area, apparent void ratio and the number of pores generally show a decreasing trend; (3) under the action of a consolidation load, when the microstructure at a relatively large scale is basically stable, the microstructure at a smaller scale will continue to adjust; (4) the apparent void ratio has a good linear regression relationship with the conventional void ratio, and the apparent void ratio has a good exponential growth relationship with the compressibility.

Evolution Assessment of Mining Subsidence Characteristics Using SBAS and PS Interferometry in Sanshandao Gold Mine, China

Ground subsidence is a common geological phenomenon occurring in mining areas. As an important Chinese gold mine, Sanshandao Gold Mine has a mining history of 25 years, with remarkable ground subsidence deformation. Mining development, life security, property security and ecological protection all require comprehension of the ground subsidence characteristics and evolution in the mining area. In this study, the mining subsidence phenomenon of the Sanshandao Gold Mine was investigated and analyzed based on Persistent Scatterer Interferometry (PSI) and small baseline subset (SBAS). The SAR (synthetic aperture radar) images covering the study area were acquired by the Sentinel-1A satellite between 2018 and 2021; 54 images (between 22 February 2018 and 25 May 2021) were processed using the PSI technique and 24 images (between 11 April 2018 and 12 July 2021) were processed using the SBAS technique. In addition, GACOS (generic atmospheric correction online service) data were adopted to eliminate the atmospheric error in both kinds of data processing. The interferometric synthetic aperture radar (InSAR) results showed a basically consistent subsidence area and a similar subsidence pattern. Both InSAR results indicated that the maximum LOS (line of sight) subsidence velocity is about 49 mm/year. The main subsidence zone is situated in the main mining area, extending in the northwest and southeast directions. According to the subsidence displacement of several representative sites in the mining area, we found that the PSI result has a higher subsidence displacement value compared to the SBAS result. Mining activities were accompanied by ground subsidence in the mining area: the ground subsidence phenomenon is exacerbated by the increasing mining quantity. Temporally, the mining subsidence lags behind the increase in mining quantity by about three months. In summary, the mining area has varying degrees of ground subsidence, monitored by two reliable time-series InSAR techniques. Further study of the subsidence mechanism is necessary to forecast ground subsidence and instruct mining activities.

Investigating Ecosystem Service Trade-Offs/Synergies and Their Influencing Factors in the Yangtze River Delta Region, China

A comprehensive understanding of the ecosystem services (ESs) trade-off/synergy relationships has become increasingly important for ecological management and sustainable development. This study employed the Yangtze River Delta (YRD) region in China as the study area and investigated the spatiotemporal changes in three ESs, namely, carbon storage (CS), water purification (WP), and habitat quality (HQ). A trade-off/synergy degree (TSD) indicator was developed that allowed for the quantification of the trade-off/synergy intensity, and the spatial pattern of the TSD between ESs in the YRD region to be analyzed. Furthermore, a geographically weighted regression (GWR) model was used to analyze the relationship between the influencing factors and trade-offs/synergies. The results revealed that CS, WP, and HQ decreased by 0.28%, 2.49%, and 3.38%, respectively, from 2005 to 2015. The TSD indicator showed that the trade-off/synergy relationships and their magnitudes were spatially heterogeneous throughout the YRD region. The coefficients of the natural and socioeconomic factors obtained from the GWR indicated that their impacts on the trade-offs/synergies vary spatiotemporally. The impact factors had both positive and negative effects on the trade-offs/synergies. The findings of this study could improve the understanding of the spatiotemporal dynamics of trade-offs/synergies and their spatially heterogeneous correlations with related factors.

Spatially Varying Relationships between Land Subsidence and Urbanization: A Case Study in Wuhan, China

Land subsidence has become an increasing global concern over the past few decades due to natural and anthropogenic factors. However, although several studies have examined factors affecting land subsidence in recent years, few have focused on the spatial heterogeneity of relationships between land subsidence and urbanization. In this paper, we adopted the small baseline subset-synthetic aperture radar interferometry (SBAS-InSAR) method using Sentinel-1 radar satellite images to map land subsidence from 2015 to 2018 and characterized its spatial pattern in Wuhan. The bivariate Moran’s I index was used to test and visualize the spatial correlations between land subsidence and urbanization. A geographically weighted regression (GWR) model was employed to explore the strengths and directions of impacts of urbanization on land subsidence. Our findings showed that land subsidence was obvious and unevenly distributed in the study area, the annual deformation rate varied from −42.85 mm/year to +29.98 mm/year, and its average value was −1.0 mm/year. A clear spatial pattern for land subsidence in Wuhan was mapped, and several apparent subsidence funnels were primarily located in central urban areas. All urbanization indicators were found to be significantly spatially correlated with land subsidence at different scales. In addition, the GWR model results showed that all urbanization indicators were significantly associated with land subsidence across the whole study area in Wuhan. The results of bivariate Moran’s I and GWR results confirmed that the relationships between land subsidence and urbanization spatially varied in Wuhan at multiple spatial scales. Although scale dependence existed in both the bivariate Moran’s I and GWR models for land subsidence and urbanization indicators, a “best” spatial scale could not be confirmed because the disturbance of factors varied over different sampling scales. The results can advance the understanding of the relationships between land subsidence and urbanization, and they will provide guidance for subsidence control and sustainable urban planning.

Investigation on spatial variability and influencing factors of drinking water iodine in Xinjiang, China

Background and objectives

Xinjiang is one of the areas in China with extremely severe iodine deficiency. The health of Xinjiang residents has been endangered for a long time. In order to provide reasonable suggestions for scientific iodine supplementation and improve the health and living standards of the people in Xinjiang, it is necessary to understand the spatial distribution of iodine content in drinking water and explore the influencing factors of spatial heterogeneity of water iodine content distribution.

Methods

The data of iodine in drinking water arrived from the annual water iodine survey in Xinjiang in 2017. The distribution of iodine content in drinking water in Xinjiang is described from three perspectives: sampling points, districts/counties, and townships/streets. ArcGIS was used for spatial auto-correlation analysis, mapping the distribution of iodine content in drinking water and visualizing the distribution of Geographically Weighted Regression (GWR) model parameter. Kriging method is used to predict the iodine content in water at non-sampling points. GWR software was used to build GWR model in order to find the factors affecting the distribution of iodine content in drinking water.

Results

There are 3293 sampling points in Xinjiang. The iodine content of drinking water ranges from 0 to 128 μg/L, the median is 4.15 μg/L. The iodine content in 78.6% of total sampling points are less than 10 μg/L, and only that in the 3.4% are more than 40 μg/L. Among 1054 towns’ water samples in Xinjiang, 88.9% of the samples’ water iodine content is less than 10 μg/L. Among the 94 studied areas, the median iodine content in drinking water in 87 areas was less than 10 μg/L, those values in 7 areas were between 10–40 μg/L, and the distribution of water iodine content in Xinjiang shows clustered. The GWR model established had found that the effects of soil type and precipitation on the distribution of iodine content in drinking water were statistically significant.

Conclusions

The iodine content of drinking water in Xinjiang is generally low, but there are also some areas which their drinking water has high iodine content. Soil type and precipitation are the factors affecting the distribution of drinking water iodine content, and are statistically significant (P<0.05).

Spatiotemporal modeling of land subsidence using a geographically weighted deep learning method based on PS-InSAR

The demand for water resources during urbanization forces the continuous exploitation of groundwater, resulting in dramatic piezometric drawdown and inducing regional land subsidence (LS). This has greatly threatened sustainable development in the long run. LS modeling helps understanding the factors responsible for the ongoing loss of land elevation and hence enhances the development of prevention strategies. Data-driven LS models perform well with fewer variables and faster convergence than physically-based hydrogeological models. However, the former models often cannot simultaneously reflect the temporal nonlinearity and spatial correlation (SC) characteristics of LS under complex variables. We proposed a LS spatiotemporal model which considers both nonlinear and spatial correlations between LS and groundwater level change of exploited aquifers. It is based on deep learning method and LS time series detected by permanent scatterer-interferometric synthetic aperture radar (PS-InSAR). The LS time series and hydrogeological properties are constructed as a spatiotemporal dataset for model training. The spatiotemporal LS model, geographically weighted long short-term memory (GW-LSTM), is constructed by integrating SC with LSTM. This latter is a deep recurrent neural network approach incorporating sequential data. The model is validated by a case study in the Beijing plain. The results show that the accuracy of the proposed model can be greatly improved considering the spatial correlation between LS and influencing factors. Furthermore, the comparison between the LSTM and GW-LSTM models reveals that groundwater level variation is not a unique causation of LS in the study area. The developed model deals with the spatiotemporal characteristics of LS under multiple variables and can be used to predict LS under different scenarios of groundwater level variations for the purpose of monitoring and providing evidence to support the prevention of future LS.

Comparative Study of Groundwater-Induced Subsidence for London and Delhi Using PSInSAR

Groundwater variation can cause land-surface movement, which in turn can cause significant and recurrent harm to infrastructure and the water storage capacity of aquifers. The capital cities in the England (London) and India (Delhi) are witnessing an ever-increasing population that has resulted in excess pressure on groundwater resources. Thus, monitoring groundwater-induced land movement in both these cities is very important in terms of understanding the risk posed to assets. Here, Sentinel-1 C-band radar images and the persistent scatterer interferometric synthetic aperture radar (PSInSAR) methodology are used to study land movement for London and National Capital Territory (NCT)-Delhi from October 2016 to December 2020. The land movement velocities were found to vary between −24 and +24 mm/year for London and between −18 and +30 mm/year for NCT-Delhi. This land movement was compared with observed groundwater levels, and spatio-temporal variation of groundwater and land movement was studied in conjunction. It was broadly observed that the extraction of a large quantity of groundwater leads to land subsidence, whereas groundwater recharge leads to uplift. A mathematical model was used to quantify land subsidence/uplift which occurred due to groundwater depletion/rebound. This is the first study that compares C-band PSInSAR-derived land subsidence response to observed groundwater change for London and NCT-Delhi during this time-period. The results of this study could be helpful to examine the potential implications of ground-level movement on the resource management, safety, and economics of both these cities.

Spatiotemporal analysis of COVID-19 outbreaks in Wuhan, China

Few study has revealed spatial transmission characteristics of COVID-19 in Wuhan, China. We aimed to analyze the spatiotemporal spread of COVID-19 in Wuhan and its influence factors. Information of 32,682 COVID-19 cases reported through March 18 were extracted from the national infectious disease surveillance system. Geographic information system methods were applied to analysis transmission of COVID-19 and its influence factors in different periods. We found decrease in effective reproduction number (Rt) and COVID-19 related indicators through taking a series of effective public health measures including restricting traffic, centralized quarantine and strict stay-at home policy. The distribution of COVID-19 cases number in Wuhan showed obvious global aggregation and local aggregation. In addition, the analysis at streets-level suggested population density and the number of hospitals were associated with COVID-19 cases number. The epidemic situation showed obvious global and local spatial aggregations. High population density with larger number of hospitals may account for the aggregations. The epidemic in Wuhan was under control in a short time after strong quarantine measures and restrictions on movement of residents were implanted.

Characterizing the Topographic Changes and Land Subsidence Associated with the Mountain Excavation and City Construction on the Chinese Loess Plateau

A mega project, Mountain Excavation and City Construction (MECC), was launched in the hilly and gully region of the Chinese Loess Plateau in 2012, in order to address the shortage of available land and create new flat land for urban construction. However, large-scale land creation and urban expansion significantly alters the local geological environment, leading to severe ground deformation. This study investigated the topographic changes, ground deformation, and their interactions due to the MECC project in the Yan’an New District (YND). First, new surface elevations were generated using ZiYuan-3 (ZY-3) stereo images acquired after the construction in order to map the local topographic changes and the fill thickness associated with the MECC project. Then, the interferometric synthetic aperture radar (InSAR) time series and 32 Sentinel-1A images were used to assess the spatial patterns of the ground deformation in the YND during the postconstruction period (2017–2018). By combining the InSAR-derived results and topographic change features, the relationship between the ground deformation and large-scale land creation was further analyzed. The results indicated that the MECC project in the YND has created over 22 km2 of flat land, including 10.8 km2 of filled area, with a maximum fill thickness of ~110 m. Significant uneven ground deformation was detected in the land-creation area, with a maximum subsidence rate of approximately 121 mm/year, which was consistent with the field survey. The strong correlation between the observed subsidence patterns and the land creation project suggested that this recorded uneven subsidence was primarily related to the spatial distribution of the filling works, along with the changes in the thickness and geotechnical properties of the filled loess; moreover, rapid urbanization, such as road construction, can accelerate the subsidence process. These findings can guide improvements in urban planning and the mitigation of geohazards in regions experiencing large-scale land construction.

Applications of Unmanned Aerial Systems (UASs) in Hydrology: A Review

In less than two decades, UASs (unmanned aerial systems) have revolutionized the field of hydrology, bridging the gap between traditional satellite observations and ground-based measurements and allowing the limitations of manned aircraft to be overcome. With unparalleled spatial and temporal resolutions and product-tailoring possibilities, UAS are contributing to the acquisition of large volumes of data on water bodies, submerged parameters and their interactions in different hydrological contexts and in inaccessible or hazardous locations. This paper provides a comprehensive review of 122 works on the applications of UASs in surface water and groundwater research with a purpose-oriented approach. Concretely, the review addresses: (i) the current applications of UAS in surface and groundwater studies, (ii) the type of platforms and sensors mainly used in these tasks, (iii) types of products generated from UAS-borne data, (iv) the associated advantages and limitations, and (v) knowledge gaps and future prospects of UASs application in hydrology. The first aim of this review is to serve as a reference or introductory document for all researchers and water managers who are interested in embracing this novel technology. The second aim is to unify in a single document all the possibilities, potential approaches and results obtained by different authors through the implementation of UASs.

Analysis of the Contribution Rate of the Influencing Factors to Land Subsidence in the Eastern Beijing Plain, China Based on Extremely Randomized Trees (ERT) Method

As a common geological hazard, land subsidence is widely distributed in the Eastern Beijing Plain. The pattern of evolution of this geological phenomenon is controlled by many factors, including groundwater level change in different aquifers, compressible layers of different thicknesses, and static and dynamic loads. First, based on the small baseline subset Interferometric Synthetic Aperture Radar (SBAS-InSAR) technique, we employed 47 ENVISAT ASAR images and 48 RADARSAT-2 images to acquire the ground deformation of the Beijing Plain from June 2003 to November 2015 and then validated the results using leveling benchmark monitoring. Second, we innovatively calculated additional stress to obtain static and dynamic load information. Finally, we evaluated the contribution rate of the influencing factors to land subsidence by using the Spearman’s rank correlation coefficient (SRCC) and extremely randomized trees (ERT) machine learning methods. The SBAS-InSAR outcomes revealed that the maximum deformation rate was 110.7 mm/year from 2003 to 2010 and 144.4 mm/year from 2010 to 2015. The SBAS-InSAR results agreed well with the leveling benchmark monitoring results; the correlation coefficients were 0.97 and 0.96 during the 2003–2010 and 2013–2015 periods, respectively. The contribution rate of the second confined aquifer to the cumulative land subsidence was 49.3% from 2003 to 2010, accounting for the largest proportion; however, its contribution rate decreased to 23.4% from 2010 to 2015. The contribution rate of the third confined aquifer to the cumulative land subsidence increased from 2003 to 2015. Although the contribution of additional stress engendered from static and dynamic loads to the cumulative land subsidence was slight, it had a significant effect on the uneven land subsidence, with a contribution rate of 33.8% from 2003 to 2010 and 23.1% from 2010 to 2015. These findings provide scientific support for mitigating hazards associated with land subsidence.