No abstract

International audienceThis paper presents a 3D parametric fault representation for modeling the displacement field associated with faults in accordance with their geometry. The displacements are modeled in a canonical fault space where the near-field displacement is defined by a small set of parameters consisting of the maximum displacement amplitude and the profiles of attenuation in the surrounding space. The particular geometry and the orientation of the slip of each fault are then taken into account by mapping the actual fault onto its canonical representation. This mapping is obtained with the help of a curvilinear frame aligned both on the fault surface and slip direction. This formulation helps us to include more geological concepts in quantitative subsurface models during 3D structural modeling tasks. Its applicability is demonstrated in the framework of forward modeling and stochastic sequential fault simulations, and the results of our model are compared to observations of natural objects described in the literature

International audienceA recent method for modeling folds uses a fold frame with coordinates based on the structural geology of folds: fold axis direction, fold axial surface and extension direction. The fold geometry can be characterised by rotating the fold frame by the pitch of the fold axis in the axial surface and the angle between the folded foliation and the axial surface. These rotation angles can be expressed as 1D functions of the fold frame coordinates. In this contribution we present methods for extracting and automatically modeling the fold geometries from structural data. The fold rotation angles used for characterising the fold geometry can be calculated locally from structural observations. The fold rotation angles incorporate the structural geology of the fold and allow for individual structural measurements to be viewed in the context of the folded structure. To filter out the effects of later folding the fold rotation angles are plotted against the coordinates of the fold frame. Using these plots the geometry of the folds can be interpolated directly from structural data where we use a combination of radial basis function and harmonic analysis to interpolate and extrapolate the fold geometry. This contribution addresses a major limitation in existing methods where the fold geometry was not constrained from structural data. We present two case studies: a proof of concept synthetic model of a non-cylindrical fold and an outcrop of an asymmetrical fold within the Lachlan Fold belt at Cape Conran, Victoria, Australia

No abstract

International audienceThree-dimensional structural modeling is gaining importance for a broad range of quantitative geoscientific applications. However, existing approaches are still limited by the type of structural data they are able to use and by their lack of structural meaning. Most techniques heavily rely on spatial data for modeling folded layers, but are unable to completely use cleavage and lineation information for constraining the shape of modeled folds. This lack of structural control is generally compensated by expert knowledge introduced in the form of additional interpretive data such as cross-sections and maps. With this approach, folds are explicitly designed by the user instead of being derived from data. This makes the resulting structures subjective and deterministic. This paper introduces a numerical framework for modeling folds and associated foliations from typical field data. In this framework, a parametric description of fold geometry is incorporated into the interpolation algorithm. This way the folded geometry is implicitly derived from observed data, while being controlled through structural parameters such as fold wavelength, amplitude and tightness. A fold coordinate system is used to support the numerical description of fold geometry and to modify the behavior of classical structural interpolators. This fold frame is constructed from fold-related structural elements such as axial foliations, intersection lineations, and vergence. Poly-deformed terranes are progressively modeled by successively modeling each folding event going backward through time. The proposed framework introduces a new modeling paradigm, which enables the building of three-dimensional geological models of complex poly-deformed terranes. It follows a process based on the structural geologist approach and is able to produce geomodels that honor both structural data and geological knowledge

No abstract

The Weierbach experimental catchment (0.45 km 2 ) is the most instrumented and studied sub-catchment in the Alzette River basin in Luxembourg. Within the last decade, it has matured towards an interdisciplinary critical zone observatory focusing on a better understanding of hydrological and hydro-geochemical processes. The Weierbach catchment is embedded in an elevated sub-horizontal plateau, characterized by slate bedrock and representative of the Ardennes Massif. Its climate is semimarine, with precipitation being rather evenly distributed throughout the year. Base

Abstract. Within the Critical Zone, regolith plays a key role in the fundamental hydrological function of water collection, storage, mixing and release. Electrical Resistivity Tomography (ERT) is recognized as a remarkable tool for characterizing the geometry and properties of the regolith, overcoming limitations inherent to conventional borehole-based investigations. However, ERT measurements with a high vertical resolution remain restricted to shallow depths, essentially due to the requirement of small electrode spacing increments (ESI). Under these circumstances, the use of ERT measurements for large horizontal surveys remains cumbersome and time-consuming. Here we focus on the need to optimize the ESI parameter in order to adequately characterize the subsurface fabric. We use a set of synthetic three-layered soil–saprock/saprolite–bedrock models in combination with a field dataset. We demonstrate that oversized ESI can significantly affect our perception of shallow subsurface structures by missing important layers and increasing the ill-posed inverse problem effects. More precisely, we document how a thin surficial layer can influence inverted ERT results and cause a resistivity bias, both at the surface and at deeper horizons. To overcome this limitation, we propose adding interpolated levels of surficial apparent resistivity based on a limited number of ERT profiles with small ESI. We demonstrate that our protocol significantly improves the accuracy of ERT profiles based on large ESI. Our protocol is time and cost efficient – especially in the case of large-scale ERT surveys.