Search

Aseisimic earth fissures are complex consequences of groundwater withdrawal and natural hydrogeologic conditions. This paper aims to improve the understanding of the mechanism of earth fissuring and investigate the relative importance of various factors to fissure activity, including bedrock geometry, piezometric depletion, compressibility and thickness of the exploited aquifer. For these purposes, a test case characterized by an impermeable and incompressible rock ridge in a subsiding basin is developed, where stress/displacement analyses and fissure state are predicted using an interface-finite element model. Three different methods for global sensitivity analysis are used to quantify the extent of the fissure opening to the aforementioned factors. The conventional sampling based Sobol' sensitivity analysis is compared to two surrogate based methods, the general polynomial chaos expansion based Sobol' analysis and a feature importance evaluation of a gradient boosting decision tree model. Numerical results indicate that earth fissure is forming in response to tensile stress accumulation above the ridge associated to pore-pressure depletion, inducing the fissure opening at land surface with further downward propagation. Sensitivity analysis highlights that the geometry of bedrock ridge is the most influential feature. Specifically, the fissure grows more when the ridge is steeper and closer to the land surface. Pore pressure depletion is a secondary feature and required to reach a certain threshold to activate the fissure. As for this specific application, the gradient boosting tree is the most suitable method for its better performance in capturing fissure characteristics. [ABSTRACT FROM AUTHOR]

Anthropogenic land subsidence can be evaluated and predicted by numerical models, which are often built over deterministic analyses. However, uncertainties and approximations are present, as in any other modeling activity of real-world phenomena. This study aims at combining data assimilation techniques with a physically-based numerical model of anthropogenic land subsidence in a novel and comprehensive workflow, to overcome the main limitations concerning the way traditional deterministic analyses use the available measurements. The proposed methodology allows to reduce uncertainties affecting the model, identify the most appropriate rock constitutive behavior and characterize the most significant governing geomechanical parameters. The proposed methodological approach has been applied in a synthetic test case representative of the Upper Adriatic basin, Italy. The integration of data assimilation techniques into geomechanical modeling appears to be a useful and effective tool for a more reliable study of anthropogenic land subsidence. [ABSTRACT FROM AUTHOR]

The reactivation of faults and the generation of fractures can be caused by stress changes due to injection and/or production of fluids into and/or from the subsurface. The simulation of these processes, which could be associated with (micro-)seismicity, is affected by a high uncertainty. The aim of this work is at developing a mathematical framework to quantify and possibly reduce the prior modeling uncertainties by assimilation of seismic data. The mechanics of fault (re-)activation is simulated by a Finite Element (FE) numerical model where the discontinuous displacements between the fault surfaces are suitably considered using appropriate Interface Elements (IEs). The study is carried out by using a stochastic approach, with a global sensitivity analysis (gSA) based on Sobol' indices initially performed to estimate the influence of the input parameters on the model solution. Then, a Markov Chain Monte Carlo (MCMC) sampling technique based on the generalized Polynomial Chaos expansion (gPC) surrogate solution is used to update the prior information conditioned on seismic observations. The methodology is tested on a 3D synthetic test case. The uncertain input is the natural stress regime and the Mohr-Coulomb parameters characterizing the fault activation criterion. A good reduction of the prior uncertainty is obtained, showing that the assimilation of seismic data can have a promising potential for improving the subsurface characterization. [ABSTRACT FROM AUTHOR]

A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10079-6 [ABSTRACT FROM AUTHOR]

The Vietnamese Mekong Delta was formed by rapid transgression during the second half of the Holocene by deposition of mainly unconsolidated, fine-grained (clayey) sediments undergoing high compaction rates. The natural subsidence can seriously impact the already vulnerable delta plain as its low elevation exposes the delta to global sea level rise, flooding, salinization. Human activities such as groundwater pumping, infrastructural loading, sand mining and dam construction have exacerbated the effects of natural consolidation. Here we present a novel modeling study that has allowed to reproduce the formation and evolution of the Mekong delta over the past 4000 years. Using an adaptive finite-element mesh, the model properly simulates accretion and natural consolidation characterizing the delta evolution. Large soil grain motion and the delayed dissipation of pore-water overpressure are accounted for. We find that natural compaction of Holocene deposits following delta evolution exceeds predicted values of absolute sea level rise. The unprecedented high rates (up to ~20 mm/yr) threaten the lower delta plain with permanent inundation and inevitably reduce the designed service life of flood defense structures along the coast. Total subsidence and sediment delivery to the delta plain will determine its future elevation and vulnerability to relative sea level rise. [ABSTRACT FROM AUTHOR]