Search

Abstract In coal mining areas, surface subsidence poses significant risks to human life and property. Fortunately, surface subsidence caused by coal mining can be monitored and predicted by using various methods, e.g., probability integral method and deep learning (DL) methods. Although DL methods show promise in predicting subsidence, they often lack accuracy due to insufficient consideration of spatial correlation and temporal nonlinearity. Considering this issue, we propose a novel DL-based approach for predicting mining surface subsidence. Our method employs K-means clustering to partition spatial data, allowing the application of a gate recurrent unit (GRU) model to capture nonlinear relationships in subsidence time series within each partition. Optimization using snake optimization (SO) further enhances model accuracy globally. Validation shows our method outperforms traditional Long Short-Term Memory (LSTM) and GRU models, achieving 99.1% of sample pixels with less than 8 mm absolute error.

Hydrological changes in groundwater coupled with earthquakes had been documented in previous studies by global researchers. Although few reports investigate multiple geochemical parameters that respond to earthquakes, trace elements received less attention, whereas they were suggested to be more sensitive to small earthquakes than the commonly used geochemical parameters. Beijing is located in the Zhangjiakou-Bohai (Zhang-Bo) seismic belt of North China, and although the occurrence of small earthquakes is frequent, the great historic earthquake in the Sanhe-Pinggu area M8 in 1679 in the adjoining southeast of Beijing gained widespread public attention. To find effective precursors that are significant for operational earthquake forecasting of the Beijing area, we carried out a one year test research project through weekly collection of groundwater samples during June 2021 to June 2022 from the seismic monitoring well of Wuliying in northwest Beijing. The 41 trace elements chemical compositions were analyzed for each sample. During the project ongoing period, the biggest earthquake with a magnitude of ML3.3 occurred in the Chaoyang District of Beijing on 3 February 2022. The content changes in these trace elements were systematically monitored before and after the earthquake. Through retrospective research, it was found that a few sensitive trace elements were anomalous to be coupled to the earthquake, including Li, Sc, Rb, Mo, Cs, Ba, W, U, Sr, Mn, Ni, and Zn. In addition to trace elements, we examined stable isotopes of hydrogen and oxygen and the existing hydrological data on groundwater level, temperature, major ions, and gases to assess the validity of geochemistry as a monitoring and predictive tool. We only found that F- (fluorine) ions and He (helium) gas had apparent shifts related to the earthquakes, while no shifts in the groundwater level were observed. Such characteristics of multiple geochemical parameters indicate that trace elements are likely to be more sensitive to crustal strain than the groundwater level and major ions. We assumed a most likely mechanism of the combination of mixing and water–rock interactions to explain the phenomenon. The probable scenario was that minor stresses caused by the earthquakes might create micro-cracks in bedrocks, thereby leading to a small volume of chemically distinct water mixing with the original water of the aquifer, and finally, the earthquake-induced rock fractures enhance the water–rock interactions, resulting in the post-seismic recovery of trace elements and δ18O value migration to the GWML. More testing works to find other sensitive sites to investigate multiple geochemical characteristics aiming at long-term to short-term earthquake prediction in the Beijing area and Zhang-Bo seismic belt are in progress.

AbstractTo systematically and thoroughly analyze the sensitivity and reliability of the MSBAS regularization for the estimation of surface deformation over a mine in combination with an application example, this study processed 101 Sentinel-1A/B SAR images, constructed and solved the 2D deformation models using SVD and Tikhonov regularization methods with different orders and parameters, and estimated the vertical and east-west surface deformation time series in a mine of China. Then, this study collected the leveling-monitoring vertical surface deformation data on three leveling points, and compared and analyzed the sensitivity and reliability of the MSBAS regularization methods for estimating vertical surface deformation. The results indicate that different regularization orders and parameters can lead to thousands of times differences in condition numbers and significant differences in ill-posed degree of the deformation models. The zero-order Tikhonov regularized deformation model with regularization parameter of 0.1 has the minimum condition number and the equation is not ill-posed. The first-order Tikhonov regularized deformation model with regularization parameter of 0.001 has the maximum condition number and the equation is seriously ill-posed. As a result, the estimates of 2D surface deformation models with different parameters and orders are also different in terms of numerical values and accuracy. Compared with the leveling-monitoring data, the first- and second-order MSBAS regularization methods with parameter 0.1 have the minimum fluctuation and the maximum correlation coefficients between the estimated values and the leveling-monitoring values, and are also closest to the leveling-monitoring results with the highest accuracy.

Backgrounds and aimsImmunotherapies have formed an entirely new treatment paradigm for hepatocellular carcinoma (HCC). Tertiary lymphoid structure (TLS) has been associated with good response to immunotherapy in most solid tumors. Nonetheless, the role of TLS in human HCC remains controversial, and recent studies suggest that their functional heterogeneity may relate to different locations within the tumor. Exploring factors that influence the formation of TLS in HCC may provide more useful insights. However, factors affecting the presence of TLSs are still unclear. The human gut microbiota can regulate the host immune system and is associated with the efficacy of immunotherapy but, in HCC, whether the gut microbiota is related to the presence of TLS still lacks sufficient evidence.MethodsWe performed pathological examinations of tumor and para-tumor tissue sections. Based on the location of TLS in tissues, all patients were divided into intratumoral TLS (It-TLS) group and desertic TLS (De-TLS) group. According to the grouping results, we statistically analyzed the clinical, biological, and pathological features; preoperative gut microbiota data; and postoperative pathological features of patients.ResultsIn a retrospective study cohort of 60 cases from a single center, differential microbiota analysis showed that compared with the De-TLS group, the abundance of Lachnoclostridium, Hungatella, Blautia, Fusobacterium, and Clostridium was increased in the It-TLS group. Among them, the enrichment of Lachnoclostridium was the most significant and was unrelated to the clinical, biological, and pathological features of the patients. It can be seen that the difference in abundance levels of microbiota is related to the presence of TLS.ConclusionOur findings prove the enrichment of Lachnoclostridium-dominated gut microbiota is associated with the presence of It-TLS in HCC patients.

The energy inside the Earth can not only be released outward through earthquakes and volcanoes but also can be used by humans in the form of geothermal energy. Is there a correlation between different forms of energy release? In this contribution, we perform detailed seismic and geothermal research in the Beijing area. The results show that the geothermal resources in Beijing belong to typical medium-low temperature geothermal resources of the sedimentary basin, and some areas are controlled by deep fault activities (e.g., Xiji geothermal well (No. 17)). The heat sources are upper mantle heat, radioactive heat in granite, and residual heat from magma cooling. The high overlap of earthquakes and geothermal field locations and the positive correlation between the injection water and earthquakes indicate that the exploitation and injection water will promote the release of the earth’s energy. The energy releases are partitioned into multiple microearthquakes, avoiding damaging earthquakes (ML ≥ 5) due to excessive energy accumulation. Therefore, the exploitation of geothermal resources may be one way to reduce destructive earthquakes. Furthermore, the use of geothermal resources can also reduce the burning of fossil energy, which is of great significance in dealing with global warming.

To quantitatively investigate the relationship between earthquakes and ionospheric anomalies, this paper presents a statistical study of pre-earthquake vertical total electron content (VTEC) variations. A total of 1522 shallow (≤60 km) strong (Mw≥6.0) earthquakes in the global area during 2000-2020 are selected, and classified according to different magnitudes, latitudes and focal depths. A quartile-based process with different lengths of sliding windows, equaling 10 days, 15 days and 27 days, respectively, has been utilized to detect VTEC anomalies. The abnormal level is first defined, and then VTEC anomalies occurrence probabilities (Po) and occurrence rates (PE) within 1-10 days before 1522 earthquakes have been calculated. Besides, VTEC anomalies occurrence rates of the background days (PN) are also calculated. The results show that the significant correlation between Po and epicentral latitudinal locations could be observed within 1-10 days before earthquakes. The values of Po increase with larger magnitudes in the equatorial and low-latitude regions, but decrease with greater magnitudes in the mid- and high-latitude regions to some degree. Within 1-5 days before earthquakes, the overall trend of PE shows an increase with larger magnitudes, but the correlation between the values of PE and magnitudes is relatively weak in the southern mid- and high-latitude regions. There is no evident causality between PN and the magnitude, and most of the values of PE/PN are larger than 1, indicating that VTEC anomalies within a few days before earthquakes are probably related with the forthcoming earthquakes. Moreover, when the abnormal level exceeds 60%, different sliding window lengths have a significant impact on the values of Po and PE in the mid- and high-latitude regions. In particular, there are obvious systematic deviations between the values of Po obtained from different sliding windows in the southern mid- and high-latitude regions. However, the selection of the optimal sliding window needs to be further studied.

This study investigated the role of the number of differentialinterferograms and coherent threshold values on the accuracy of SBAS InSAR (small baseline subset interferometric synthetic aperture radar) results for specific applications in Jiyang Coal. Fifty-eight imageswere utilized to form fourdifferential interferogram timeseries with different numbers of interferograms and coherence thresholds. The four SBAS InSAR-monitored results, mining information of 15 working faces, and levelling-monitored results of 260 levelling points werecompared.The greater number of differential interferograms and lower coherent threshold values could better reflect the spatial distribution and variation trend of mining subsidence and demonstrate the advantages of SBAS InSAR. However, an excessive number of differential interferograms and excessively low coherent threshold values would increase the data processingdifficulty and the differences between SBAS InSAR- and levelling-monitoredresults. Location, spatial distributionandscopeofSBAS InSAR-monitored mining subsidence were consistent with the mining progress of the working faces.The accuracy ofSBASInSAR-monitored subsidence values followed a certain spatio-temporal variation law. The variation trend of absolute differences between SBAS InSAR-and levelling-monitored subsidence values was similar to the shape of the levelling-monitoredsubsidence basin, which was helpful to study the correction method of SBAS InSAR-monitored results.

Groundwater overexploitation is a critical issue in the North China Plain (NCP), resulting in groundwater level decline and surface subsidence for the last half-century. This problem, however, has been greatly alleviated by the South-to-North Water Diversion (SNWD) Project since 2015. Monitoring of this process has been steadily improved in recent years using water level and geodetic observations. Here, we characterize the water storage change at the Huairou groundwater reserve site (HGRS) in Beijing due to the SNWD by combining Interferometric Synthetic Aperture Radar (InSAR) data of the Sentinel-1 satellites, continuous Global Positioning System (GPS) data, and well water level data observed during the same time. InSAR observations revealed subsidence up to ~400 mm in the Beijing plain but uplift at ~40 mm in the HGRS during 2015–2019, and more than 70% of the uplift occurred from October 2018 to January 2019. By integrating the most significant uplift deformation during October 2018 to January 2019 with water level observations at the same time, we estimated the storativity of the confined aquifer system at HGRS as 1.68–7.82×10−3, weighing in the correction for effective stress and surface deformation for various situations. Based on the estimated aquifer storativity and the observed water level change in the unconfined and confined aquifer, the recharged water storage for the confined and unconfined aquifers was estimated as 1.20–1.39×107 m3 and ~2.86×108 m3 from 6 October 2018 to 22 January 2019, respectively, which is about 4% and 91% of the surface water recharge through river channels in the same period due to the SNWD Project. Our study demonstrates that integration of geodetic and hydrological data can provide crucial information for the assessment of groundwater circulation and assistance of groundwater management.

Three-dimensional surface deformation is resolved using deformations obtained from multiple line of sight differential interferometry synthetic aperture radar (multi-LOS D-InSAR). However, the random model is susceptible to both accidental errors (such as errors due to noises) and system errors (such as orbital errors and topographic errors). To solve this problem, we present a weight-determining method based on high-precision global positioning system (GPS) observations. Specifically, the precision of multi-LOS D-InSAR deformation observations is evaluated based on high-precision GPS deformation observations from the same period. Then, by considering the influence of system errors and accidental errors on the deformation of a pixel observed by multi-LOS D-InSAR, a combined fitting model associated with the positions and elevations of a pixel is built. Finally, the precision of deformations from multi-LOS D-InSAR is evaluated pixel-by-pixel by using the built model. As a result, the weights of observations in the model for resolving three-dimensional surface deformation using multi-LOS D-InSAR can be determined. We verify the feasibility and precision of the weight-determining method through a simulation experiment and a real data experiment, respectively. The verification results show that the proposed weight-determining method based on high-precision GPS observations provides higher reliability for the resolved three-dimensional surface deformation.

The optimal selection of multi-temporal differential interferogram series is an important step to monitor ground subsidence using the Stanford method for persistent scatterers (StaMPS)-based small baseline subset (SBAS) interferometric synthetic aperture radar (InSAR). Using a deformation model and its two solution methods, least squares and singular value decomposition, we present the composing mode and optimal selection of multi-temporal differential interferogram series and show that their quality and quantity affect the accuracy of monitored deformation information of SBAS InSAR. Using 29 ENVISAT ASAR images covering urban areas of Beijing, China, a different number of optimal multi-temporal differential interferogram series are formed to monitor urban ground subsidence by the StaMPS-based SBAS method. The comparison and verification of test results indicate that the quality and quantity of multi-temporal differential interferogram series substantially impact the singularity and degree of ill condition of the deformation model, locations of the selected slowly decorrelating filtered phase (SDFP) pixels, and monitored annual mean subsidence velocities. The suitable number of multi-temporal differential interferogram series under the optimal quality is 1–2 for a month in urban ground subsidence monitoring using StaMPS-based SBAS InSAR, a higher quantity of differential interferograms of the optimal quality is not always better.

Beijing is one of the most water-stressed cities in the world. Due to over-exploitation of groundwater, the Beijing region has been suffering from land subsidence since 1935. In this study, the Small Baseline InSAR technique has been employed to process Envisat ASAR images acquired between 2003 and 2010 and TerraSAR-X stripmap images collected from 2010 to 2011 to investigate land subsidence in the Beijing region. The maximum subsidence is seen in the eastern part of Beijing with a rate greater than 100 mm/year. Comparisons between InSAR and GPS derived subsidence rates show an RMS difference of 2.94 mm/year with a mean of 2.41 ± 1.84 mm/year. In addition, a high correlation was observed between InSAR subsidence rate maps derived from two different datasets (i.e., Envisat and TerraSAR-X). These demonstrate once again that InSAR is a powerful tool for monitoring land subsidence. InSAR derived subsidence rate maps have allowed for a comprehensive spatio-temporal analysis to identify the main triggering factors of land subsidence. Some interesting relationships in terms of land subsidence were found with groundwater level, active faults, accumulated soft soil thickness and different aquifer types. Furthermore, a relationship with the distances to pumping wells was also recognized in this work.

In this paper, we tried to reveal the characteristics of the crustal deformation in both the horizontal and vertical directions in Beijing using InSAR and GNSS observations. Regarding the serious land subsidence in Beijing plain, we also analysed the mechanism of the occurrence and development of the subsidence in combination with the tectonic settings. The GNSS results reveal that the crust in Beijing shows a significant left-lateral trend movement in the horizontal direction, while the vertical direction shows a gentle rise in the mountainous region and a significant subsidence in the plain area. The INSAR results shows a detailed subsidence area and the deformation characteristics were analyzed considering the fault activity. The foundation of geological structure dominates the subsiding in the Beijing Plain. The exploitation of groundwater exacerbates the level of subsidence and has new development. The active faults controlled the development of the subsiding in present days. [ABSTRACT FROM AUTHOR]

Early detection of landslides is important for disaster prevention, which was still very hard work with traditional surveying methods. Interferometric Synthetic Aperture Radar (InSAR) technology provided us the ability to monitor displacements along the slope with wide coverage and high accuracy. In this paper, we proposed a qualitatively multi-baseline DInSAR method to early detect and map the potential landslides. Two sections of China National Highway 317 and 213 were selected as study area. With this method 10 potential landslide areas were early detected and mapped in a quick and effective way. One of them (i.e. Shidaguan landslide) collapsed on August 2017, which was coincident with our results, suggesting that this method could become an effective way to acquire the landslide early detection map to assist the future disaster prevention work. [ABSTRACT FROM AUTHOR]