Ensemble smoothing of land subsidence measurements for reservoir geomechanical characterization

Geomechanical models are often used to predict the impact on land surface of fluid withdrawal from deep reservoirs, as well as investigating measures for mitigation. The ability to accurately simulate surface displacements, however, is often impaired by limited information on the geomechanical parameters characterizing the geological formations of interest. In this study, we employ an ensemble smoother, a data assimilation algorithm, to provide improved estimates of reservoir parameters through assimilation of measurements of both horizontal and vertical surface displacement into geomechanical model results. The method leverages the demonstrated potential of remote sensing techniques developed in the last decade to provide accurate displacement data for large areas of the land surface. For evaluation purposes, the methodology is applied to the case of a disk‐shaped reservoir embedded in a homogeneous, isotropic, and linearly elastic half space, subject to a uniform change in fluid pressure. Multiple sources of uncertainty are investigated, including the radius, R, the thickness, h, and the depth, c, of the reservoir; the pore pressure change, Δp; porous medium's vertical uniaxial compressibility, c_M, and Poisson's ratio, ν, and the ratio, s, between the compressibilities of the medium during loading and unloading cycles. Results from all simulations show that the ensemble smoother has the capability to effectively reduce the uncertainty associated with those parameters to which the variability and the spatial distribution of land surface displacements are most sensitive, namely, R, c, c_M, and s. These analyses demonstrate that the estimation of these parameters values depends on the number of measurements assimilated and the error assigned to the measurement values. Copyright © 2014 John Wiley & Sons, Ltd.     

Measure of the exponential splitting of the homoclinic tangle in four-dimensional symplectic mappings

Using four-dimensional symplectic maps as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to the single resonances. We measure an exponential dependence of the size of the lobes of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by the Nekhoroshev theorem. The variation of the size of the lobes turns out to be correlated to the diffusion coefficient.     

A numerical study of the hyperbolic manifolds in a priori unstable systems. A comparison with Melnikov approximations

Using numerical methods we study the hyperbolic manifolds in a model of a priori unstable dynamical system. We compare the numerically computed manifolds with their analytic expression obtained with the Melnikov approximation. We find that, at small values of the perturbing parameter, the topology of the numerically computed stable and unstable manifolds is the same as in their Melnikov approximation. Increasing the value of the perturbing parameter, we find that the stable and unstable manifolds have a peculiar topological transition. We find that this transition occurs near those values of the perturbing parameter for which the error terms of Melnikov approximations have a sharp increment. The transition value is also correlated with a change in the behaviour of dynamical quantities, such as the largest Lyapunov exponent and the diffusion coefficient.     

Low-T thermochronometric evidence for post-thrusting (< 11 Ma) exhumation in the Western Outer Carpathians,Poland

Thermal indicators record exhumation of sedimentary units from depths in excess of 6 km over most of the Outer Carpathian fold and thrust belt in Poland. Apatite fission track data, showing cooling ages ranging between 32.1 ± 4.8 and 7.0 ± 0.8 Ma, indicate that exhumation was partially coeval with shortening. However, new thermochronometric information obtained as part of this study allowed us to unravel a previously undetected, post-thrusting exhumation stage. The integration of new field data and structural analysis with low-T thermochronometry suggests that termination of thrusting – at ca. 11 Ma in the area of the present study – was followed by gravity disequilibria within the orogenic wedge. The related extension and denudation phenomena appear to have played a primary role in the recent (< 10 Ma) tectonic evolution of the Western Outer Carpathians, exerting a major control on exhumation processes in this key area of the Alpine–Carpathian mountain system.     

The role of preconditioning in the solution to FE coupled consolidation equations by Krylov subspace methods

The repeated solution in time of the linear system arising from the finite element integration of coupled consolidation equations is a major computational effort. This system can be written in either a symmetric or an unsymmetric form, thus calling for the implementation of different preconditioners and Krylov subspace solvers. The present paper aims at investigating when either a symmetric or an unsymmetric approach should be better used. The results from a number of representative numerical experiments indicate that a major role in selecting either form is played by the preconditioner rather than by the Krylov subspace method itself. Two other important issues addressed are the size of the time integration step and the possible lumping of the flow capacity matrix. It appears that ad hoc block constrained preconditioners provide the most robust algorithm independently of the time step size, lumping, and symmetry. Copyright © 2008 John Wiley & Sons, Ltd.     

Allozyme genotypes and tolerance to copper stress in Hediste diversicolor (Polychaeta: Nereididae)

This study analysed the occurrence of genotypic shifts in laboratory populations of Hediste diversicolor (Polychaeta: Nereididae) exposed to copper stress. Specimens of H. diversicolor were collected at three sites, up to 10 km apart, in the estuarine area of the Pialassa lagoons (North Adriatic Sea, Italy) and were used in acute toxicity tests. Specimens were assigned to copper exposure (0.34 mg/l Cu2+) or control conditions. Each combination of Treatment and Site was replicated in two tanks containing 35 specimens of H. diversicolor. The genotypic structure of both dead and survived specimens was analysed by allozyme electrophoresis at six loci (ALD, FH, HBDH LDH, PGI, SDH). Under copper exposure, specimens with the genotypes ALD100/100 and PGI102/102 had significantly lower mortalities than other genotypes Results were consistent across the three sites, suggesting that, under laboratory conditions, effects of copper stress on H. diversicolor is related to individual genotypes at ALD and PGI loci.     

Volcanic geomorphology and tectonics of the Aeolian archipelago (Southern Italy) based on integrated DEM data

The geomorphological and morphometric analysis of the sea floor topography surrounding the Aeolian Islands, South Tyrrhenian Sea, Italy, provides insights into the relationships between the volcanological evolution of the islands and their tectonic features. We constructed geomorphological maps of the submarine portions of the seven large edifices constituting the islands on the basis of a DEM with a 5 m resolution step. These maps include constructional and destructional landforms such as submarine volcanic vents located west of Lipari and north of Alicudi, and hummocky surfaces recognised north of Lipari and Salina. The latter landforms, together with the occurrence of large scars affecting the main edifices on land, suggest that sector collapses affected some islands. Geomorphological data indicate that the location of subaerial and submarine vents is strongly controlled by local tectonic structures striking WNW-ESE (Alicudi-Filicudi sector), NNW-SSE (Salina-Lipari-Vulcano sector) and NE-SW (Panarea-Stromboli sector). The islands can be divided into two groups on the basis of some morphometric parameters: a first group with a pancake-like shape, Dp/D (abrasion platform diameter/basal diameter) higher than 0.40 and H/D (total height/basal diameter) lower than 0.13, and a second group with a conical shape, characterised by Dp/D lower than 0.34 and H/D higher than 0.14. These ratios and other morphometric parameters reflect the different volcanological and structural evolution of the Aeolian Islands. The pancake-like shaped complexes have been created, in addition to their submarine stage, by extrusive and highly explosive activity, whereas the cone-shaped edifices have been characterised by effusive or moderate explosive activity.Editorial responsibility: C Kilburn     

Thermodynamics of aluminosilicate garnet solid solution: new experimental data,an optimized model,and thermometric applications

 We have experimentally determined the displacement of the equilibrium Grossular + 2 Kyanite + Quartz ⇆ 3 Anorthite (GASP) as a function of garnet composition in the systems Mg-Ca-Mn, Fe-Mg-Ca and Fe-Mg-Ca-Mn at 1000°C. The results were treated along with selected experimental and observational data available in the literature as well as binary parameters from other workers to obtain a set of mutually compatible binary mixing parameters of the quaternary (Fe,Mg,Ca,Mn)- aluminosilicate garnet solid solution. Attempts to determine equilibrium garnet composition in the GASP equilibrium in the Ca-Mg binary were unsuccessful due to the formation of pyroxene. Calculations of binary and ternary miscibility gaps show that the P,T,X combination required for unmixing of garnet solid solution is not realized by natural samples. The solution model was applied to account for compositional effects on Fe-Mg exchange between garnet and ortho- or clino-pyroxene. Applications of the revised thermometric formulations to selected natural assemblages yield P-T conditions which are much less sensitive to compositional effects compared to the other available formulations, and are consistent with independent constraints. Received: 21 April 1995 / Accepted: 17 May 1996     

Large roots dominate the contribution of trees to slope stability

Tree roots provide surface erosion protection and improve slope stability through highly complex interactions with the soil due to the nature of root systems. Root reinforcement estimation is usually performed by in situ pullout tests, in which roots are pulled out of the soil to reliably estimate the root strength of compact soils. However, this test is not suitable for the scenario where a soil progressively fails in a series of slump blocks – for example, in unsupported soils near streambanks and road cuts where the soil has no compressive resistance at the base of the hillslope. The scenario where a soil is unsupported on its downslope extent and progressively deforms at a slow strain rate has received little attention, and we are unaware of any study on root reinforcement that estimates the additional strength provided by roots in this situation. We therefore designed two complementary laboratory experiments to compare the force required to pull the root out. The results indicate that the force required to pull out roots is reduced by up to 50% when the soil fails as slump blocks compared to pullout tests. We also found that, for slump block failure, roots had a higher tendency to slip than to break, showing the importance of active earth pressure on root reinforcement behaviour, which contributes to reduced friction between soil and roots. These results were then scaled up to a full tree and tree stand using the root bundle and field-measured spatial distributions of root density. Although effects on the force mobilized in small roots can be relevant, small roots have virtually no effect on root reinforcement at the tree or stand scale on hillslopes. When root distribution has a wide range of diameters, the root reinforcement results are controlled by large roots, which hold much more force than small roots. © 2019 John Wiley & Sons, Ltd.