Multi-scale aquifer characterization and groundwater flow model parameterization using direct push technologies     

Bart Rogiers  Thomas Vienken  Matej Gedeon  Okke Batelaan  Dirk Mallants  Marijke Huysmans  Alain Dassargues 《Environmental Earth Sciences》2014,72(5):1303-1324

Direct push (DP) technologies are typically used for cost-effective geotechnical characterization of unconsolidated soils and sediments. In more recent developments, DP technologies have been used for efficient hydraulic conductivity (K) characterization along vertical profiles with sampling resolutions of up to a few centimetres. Until date, however, only a limited number of studies document high-resolution in situ DP data for three-dimensional conceptual hydrogeological model development and groundwater flow model parameterization. This study demonstrates how DP technologies improve building of a conceptual hydrogeological model. We further evaluate the degree to which the DP-derived hydrogeological parameter K, measured across different spatial scales, improves performance of a regional groundwater flow model. The study area covers an area of ~60 km² with two overlying, mainly unconsolidated sand aquifers separated by a 5–7 m thick highly heterogeneous clay layer (in north-eastern Belgium). The hydrostratigraphy was obtained from an analysis of cored boreholes and about 265 cone penetration tests (CPTs). The hydrogeological parameter K was derived from a combined analysis of core and CPT data and also from hydraulic direct push tests. A total of 50 three-dimensional realizations of K were generated using a non-stationary multivariate geostatistical approach. To preserve the measured K values in the stochastic realizations, the groundwater model K realizations were conditioned on the borehole and direct push data. Optimization was performed to select the best performing model parameterization out of the 50 realizations. This model outperformed a previously developed reference model with homogeneous K fields for all hydrogeological layers. Comparison of particle tracking simulations, based either on the optimal heterogeneous or reference homogeneous groundwater model flow fields, demonstrate the impact DP-derived subsurface heterogeneity in K can have on groundwater flow and solute transport. We demonstrated that DP technologies, especially when calibrated with site-specific data, provide high-resolution 3D subsurface data for building more reliable conceptual models and increasing groundwater flow model performance.  相似文献   

Using Multiple-Point Geostatistics for Tracer Test Modeling in a Clay-Drape Environment with Spatially Variable Conductivity and Sorption Coefficient   总被引：1，自引：0，他引：1  

Marijke Huysmans  Philippe Orban  Elke Cochet  Mathias Possemiers  Benedicta Ronchi  Katherine Lauriks  Okke Batelaan  Alain Dassargues 《Mathematical Geosciences》2014,46(5):519-537

This study investigates the effect of fine-scale clay drapes on tracer transport. A tracer test was performed in a sandbar deposit consisting of cross-bedded sandy units intercalated with many fine-scale clay drapes. The heterogeneous spatial distribution of the clay drapes causes a spatially variable hydraulic conductivity and sorption coefficient. A fluorescent tracer (sodium naphthionate) was injected in two injection wells and ground water was sampled and analyzed from five pumping wells. To determine (1) whether the fine-scale clay drapes have a significant effect on the measured concentrations and (2) whether application of multiple-point geostatistics can improve interpretation of tracer tests in media with complex geological heterogeneity, this tracer test is analyzed with a local three-dimensional ground-water flow and transport model in which fine-scale sedimentary heterogeneity is modeled using multiple-point geostatistics. To reduce memory needs and calculation time for the multiple-point geostatistical simulation step, this study uses the technique of direct multiple-point geostatistical simulation of edge properties. Instead of simulating pixel values, model cell edge properties indicating the presence of irregularly shaped surfaces are simulated using multiple-point geostatistical simulations. Results of a sensitivity analysis show under which conditions clay drapes have a significant effect on the concentration distribution. Calibration of the model against measured concentrations from the tracer tests reduces the uncertainty on the clay-drape parameters. The calibrated model shows which features of the breakthrough curves can be attributed to the geological heterogeneity of the aquifer and which features are caused by other processes.  相似文献   

Using geodetic GPS receivers to measure vegetation water content   总被引：1，自引：0，他引：1  

Wei Wan  Kristine M. Larson  Eric E. Small  Clara C. Chew  John J. Braun 《GPS Solutions》2015,19(2):237-248

  相似文献   

A Full-Sun Magnetic Index from Helioseismology Inferences     

I. González Hernández  M. Díaz Alfaro  K. Jain  W. K. Tobiska  D. C. Braun  F. Hill  F. Pérez Hernández 《Solar physics》2014,289(2):503-514

Solar magnetic indices are used to model the solar irradiance and ultimately to forecast it. However, the observation of such indices is generally limited to the Earth-facing hemisphere of the Sun. Seismic maps of the far side of the Sun have proven their capability to locate and track medium–large active regions at the non-visible hemisphere. We present here the possibility of using the average signal from these seismic far-side maps, combined with similarly calculated near-side maps, as a proxy to the full-Sun magnetic activity.  相似文献   

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)     

Jens Orm  Maurits Lindstrm  Alain Lepinette  Jesús Martinez‐Frias  Enrique Diaz‐Martinez 《Meteoritics & planetary science》2006,41(10):1605-1612

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.  相似文献   

Limit analysis for saturated porous media without fluid flow calculation     

Alain Corfdir 《国际地质力学数值与分析法杂志》2004,28(6):543-561

Application of yield design to porous media usually requires a preliminary calculation of the fluid flow net. The stability analysis is then carried out with seepage forces associated with the flow. We assume here that the flow is steady and that the yield criterion is defined by a function of the effective stress tensor. The formulation that we propose here allows taking into account seepage force in the expression of the kinematic stability conditions by means of hydraulic boundary conditions without calculation of the fluid flow. One obtains a formulation of the kinematic condition similar to the case of classic, non-porous media. The method is illustrated by two examples: a cylinder subjected to fluid flow and a vertical cut. It can be adapted to various boundary conditions and to the case of a criterion defined by a function of a generalized effective stress tensor. We also give a method to derive rigorous lower bounds using approximate fluid pressure field. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Probing sunspot magnetic fields with p-mode absorption and phase shift data   总被引：1，自引：0，他引：1  

P. S. Cally †  A. D. Crouch  D. C. Braun 《Monthly notices of the Royal Astronomical Society》2003,346(2):381-389

Long-standing observations of incoming and outgoing f- and p-modes in annuli around sunspots reveal that the spots partially absorb and substantially shift the phase of waves incident upon them. The commonly favoured absorption mechanism is partial conversion to slow magneto-acoustic waves that disappear into the solar interior channelled by the magnetic field of the sunspot. However, up until now, only f-mode absorption could be accounted for quantitatively by this means. Based on vertical magnetic field models, the absorption of p-modes was insufficient. In this paper, we use the new calculations of Crouch & Cally for inclined fields, and a simplified model of the interaction between spot interior and exterior. We find excellent agreement with phase shift data assuming field angles from the vertical in excess of 30° and Alfvén/acoustic equipartition depths of around 600–800 km. The absorption of f-modes produced by such models is considerably larger than is observed, but consistent with numerical simulations. On the other hand, p-mode absorption is generally consistent with observed values, up to some moderate frequency dependent on radial order. Thereafter, it is too large, assuming absorbing regions comparable in size to the inferred phase-shifting region. The excess absorption produced by the models is in stark contrast with previous calculations based on a vertical magnetic field, and is probably due to finite mode lifetimes and excess emission in acoustic glories. The excellent agreement of phase shift predictions with observational data allows some degree of probing of subsurface field strengths, and opens up the possibility of more accurate inversions using improved models. Most importantly, though, we have confirmed that slow mode conversion is a viable, and indeed the likely, cause of the observed absorption and phase shifts.  相似文献   

Infrared interferometry with single telescopes     

Alain Chelli 《Astrophysics and Space Science》1989,160(1-2):369-371

We examine the merits of three basic interferometric techniques using two-dimensional infrared arrays: speckle interferopmetry, pupil plane interferometry, and speckle holography.Paper presented at the Symposium on the JNLT and Related Engineering Developments, Tokyo, November 29–December 2, 1988.  相似文献   

Continuous GPS and broad-scale deformation across the Rhine Graben and the Alps     

Magdala Tesauro  Christine Hollenstein  Ramon Egli  Alain Geiger  Hans-Gert Kahle 《International Journal of Earth Sciences》2005,94(4):525-537

In order to study the ongoing tectonic deformation in the Rhine Graben area, we reconstruct the local crustal velocity and the strain rate field from GPS array solutions. Following the aim of this work, we compile the velocities of permanent GPS stations belonging to various networks (EUREF, AGNES, REGAL and RGP) in central western Europe. Moreover, the strain rate field is displayed in terms of principal axes and values, while the normal and the shear components of the strain tensor are calculated perpendicular and parallel to the strike of major faults. The results are compared with the fault plane solutions of earthquakes, which have occurred in this area. A broad-scale kinematic deformation model across the Rhine Graben is provided on the basis of tectonics and velocity results of the GPS permanent stations. The area of study is divided into four rigid blocks, between which there might be relative motions. The velocity and the strain rate fields are reconstructed along their borders, by estimating a uniform rotation for each block. The tectonic behaviour is well represented by the four-block model in the Rhine Graben area, while a more detailed model will be needed for a better reconstruction of the strain field in the Alpine region.

Magdala TesauroEmail:

  相似文献   

On different types of adjustment usable to calculate the parameters of the stream power law     

Alain Demoulin  Arnaud BeckersBenoît Bovy 《Geomorphology》2012,138(1):203-208

Model parameterization through adjustment to field data is a crucial step in the modeling and the understanding of the drainage network response to tectonic or climatic perturbations. Using as a test case a data set of 18 knickpoints that materialize the migration of a 0.7-Ma-old erosion wave in the Ourthe catchment of northern Ardennes (western Europe), we explore the impact of various data fitting on the calibration of the stream power model of river incision, from which a simple knickpoint celerity equation is derived. Our results show that statistical least squares adjustments (or misfit functions) based either on the stream-wise distances between observed and modeled knickpoint positions at time t or on differences between observed and modeled time at the actual knickpoint locations yield significantly different values for the m and K parameters of the model. As there is no physical reason to prefer one of these approaches, an intermediate least-rectangles adjustment might at first glance appear as the best compromise. However, the statistics of the analysis of 200 sets of synthetic knickpoints generated in the Ourthe catchment indicate that the time-based adjustment is the most capable of getting close to the true parameter values. Moreover, this fitting method leads in all cases to an m value lower than that obtained from the classical distance adjustment (for example, 0.75 against 0.86 for the real case of the Ourthe catchment), corresponding to an increase in the non-linear character of the dependence of knickpoint celerity on discharge.  相似文献