Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Geophysical Investigations for Hydrogeological Characterization at the Kappelen Test Site (Switzerland)

Vertical electric soundings, 2D resistivity imaging and several logging measurements were performed at Kappelen test site to identify the various geolelectric facies that allowed determining the tabular and horizontal structure of the aquifer. The surface-based geoelectric methods allowed for a reliable characterization of the overall structure and the geometry of the aquifer, while geophysical logging methods allowed for inferring detailed hydrogeophysical characteristics, such as the electrical resistivity, total porosity, global and matrix density and hydraulic conductivity. The synoptic interpretation and integration of this broad and diverse database allows for constraining the key hydrological characteristics and hence forms the basis for the detailed hydraulic modelling of flow and transport process.

Three essays on liquidity risk

Preface In this thesis we study several questions related to transaction data measured at an individual level. The questions are addressed in three essays that will constitute this thesis. In the first essay we use tick-by-tick data to estimate non-parametrically the jump process of 37 big stocks traded on the Paris Stock Exchange, and of the CAC 40 index. We separate the total daily returns in three components (trading continuous, trading jump, and overnight), and we characterize each one of them. We estimate at the individual and index levels the contribution of each return component to the total daily variability. For the index, the contribution of jumps is smaller and it is compensated by the larger contribution of overnight returns. We test formally that individual stocks jump more frequently than the index, and that they do not respond independently to the arrive of news. Finally, we find that daily jumps are larger when their arrival rates are larger. At the contemporaneous level there is a strong negative correlation between the jump frequency and the trading activity measures. The second essay study the general properties of the trade- and volume-duration processes for two stocks traded on the Paris Stock Exchange. These two stocks correspond to a very illiquid stock and to a relatively liquid stock. We estimate a class of autoregressive gamma process with conditional distribution from the family of non-central gamma (up to a scale factor). This process was introduced by Gouriéroux and Jasiak and it is known as Autoregressive gamma process. We also evaluate the ability of the process to fit the data. For this purpose we use the Diebold, Gunther and Tay (1998) test; and the capacity of the model to reproduce the moments of the observed data, and the empirical serial correlation and the partial serial correlation functions. We establish that the model describes correctly the trade duration process of illiquid stocks, but have problems to adjust correctly the trade duration process of liquid stocks which present long-memory characteristics. When the model is adjusted to volume duration, it successfully fit the data. In the third essay we study the economic relevance of optimal liquidation strategies by calibrating a recent and realistic microstructure model with data from the Paris Stock Exchange. We distinguish the case of parameters which are constant through the day from time-varying ones. An optimization problem incorporating this realistic microstructure model is presented and solved. Our model endogenizes the number of trades required before the position is liquidated. A comparative static exercise demonstrates the realism of our model. We find that a sell decision taken in the morning will be liquidated by the early afternoon. If price impacts increase over the day, the liquidation will take place more rapidly.

Artificial Urinary Sphincter based on smart materials

Objective: There are only a few established artificial urinary sphincters for treatment of incontinence. We have developed a new device composed by three parts: the actuator, three contractile rings and a control unit. The actuator is made of Nitinol fibers, driven by microprocessor. The fibers are linked to the rings placed around the urethra. They function with alternance in their open and closed position. This concept is called piano concept. With this set-up, the constant compression on the urethra is strongly reduced. Methods: Six male sheep have been used for this study. The sphincter was open each hour for a period of 10 min., to guaranty urination. The bladder was filled with water while one cuff was closed and bladder pressure was monitored. The animals were sacrificed. Two biopsies around two cuffs of each explant and all three cuffs from each explant including urethra were analyzed. Urethra not surrounded by a cuff was taken as control. Results: The pressure exerted by the sphincter around the urethra provided continence. Simulated incontinence occurred at a pressure of 1bar measured on the bladder wall using a pressure probe. The closing force of the cuff was approx. 0·7N. No difference in tissue structure and organization of the urethra with and without artificial sphincter was observed. Conclusions: This device has several advantages compared to other urinary sphincters. It is easy to implant, has no hydraulic nature and reduces ischemic injury of the urethra by the alternance of urethral part compressed. Proof of concept in vivo has been demonstrated. Other studies are planned to determine long-term outcome.

Weight-bearing bones are more sensitive to physical exercise in boys than in girls during pre- and early puberty: a cross-sectional study.

We carried out a cross-section study of the sex-specific relationship between bone mineral content and physical activity at sites with different loading in pre- and early pubertal girls and boys. There was significant sensitivity of bone mineral content of the hip to physical exercise in boys, but not in girls. BACKGROUND: Since little is known whether there are sex differences in sensitivity of bone to loading, we investigated sex differences in the cross-sectional association between measures of physical activity (PA) and bone mass and size in pre- and early pubertal children of both sexes. METHODS: We measured bone mineral content/density (BMC/BMD) and fat-free mass (FFM) in 269 6- to 13-year-old children from randomly selected schools by dual-energy X-ray absorptiometry. Physical activity (PA) was measured by accelerometers and lower extremity strength by a jump-and-reach test. RESULTS: Boys (n = 128) had higher hip and total body BMC and BMD, higher FFM, higher muscle strength and were more physically active than girls (n = 141). Total hip BMC was positively associated with time spent in total and vigorous PA in boys (r = 0.20-0.33, p < 0.01), but not in girls (r = 0.02-0.04, p = ns), even after adjusting for FFM and strength. While boys and girls in the lowest tertile of vigorous PA (22 min/day) did not differ in hip BMC (15.62 vs 15.52 g), boys in the highest tertile (72 min/day) had significantly higher values than the corresponding girls (16.84 vs 15.71 g, p < 0.05). CONCLUSIONS: Sex differences in BMC during pre- and early puberty may be related to a different sensitivity of bone to physical loading, irrespective of muscle mass.

Calibration of high-resolution geophysical data with tracer test measurements to improve hydrological predictions

Relationships between porosity and hydraulic conductivity tend to be strongly scale- and site-dependent and are thus very difficult to establish. As a result, hydraulic conductivity distributions inferred from geophysically derived porosity models must be calibrated using some measurement of aquifer response. This type of calibration is potentially very valuable as it may allow for transport predictions within the considered hydrological unit at locations where only geophysical measurements are available, thus reducing the number of well tests required and thereby the costs of management and remediation. Here, we explore this concept through a series of numerical experiments. Considering the case of porosity characterization in saturated heterogeneous aquifers using crosshole ground-penetrating radar and borehole porosity log data, we use tracer test measurements to calibrate a relationship between porosity and hydraulic conductivity that allows the best prediction of the observed hydrological behavior. To examine the validity and effectiveness of the obtained relationship, we examine its performance at alternate locations not used in the calibration procedure. Our results indicate that this methodology allows us to obtain remarkably reliable hydrological predictions throughout the considered hydrological unit based on the geophysical data only. This was also found to be the case when significant uncertainty was considered in the underlying relationship between porosity and hydraulic conductivity.

Mir211 Is Dysregulated In Conjunctival Melanocytic proliferations

Purpose Downregulation of TRPM1 mRNA, a transient receptor potential cation channel, has been identified in highly metastatic cutaneous melanoma cell lines. TRPM1 mRNA expression is inversely correlated with skin melanoma metastases. Recent evidence has demonstrated that the tumor suppressive activity of TRPM1 is due to miR211 situated in intron 6 of TRPM1. As we have previously identified a downregulation of TRPM1 mRNA expression in conjunctival melanoma, we decided to assess miR211 expression and its potential target gene IGF2R and KCNMA1 in conjunctival melanocytic proliferations. As MITF has been shown to regulate both TRPM1 and mir211 expression, we also assessed MITF expression in our series. Methods Expression of miR211 was assessed by in situ hybridization in 14 conjunctival naevi and 14 conjunctival melanoma. Integrity of miRNA in tissues was evaluated in each sample with the preservation of miR126 expression in endothelial cells. Protein expression of MITF, IGF2R and KCNMA1 was assessed by immunohistochemistry. Statistical analysis was performed with JUMP 8,0 software. In situ hybridization and immunohistochemistry were assessed independently by two observers. Results There were 7 subepithelial nevi and 7 compound nevi. There were 5 female and 9 male. The population mean age was 48.7 ± 6.4 years (SEM). miR211 was found in 11 nevi (79%). MITF was expressed in all the nevi. IGF2R was found in 13 nevi. KCNMA1 was found in 57% of the nevi.The melanoma group was composed of 9 females and 5 males with a mean age of 67 ± 4.8 years (SEM). Using the recent TNM classification, 5 tumors were belonging to the T1, 3 to theT2 and 6 to the T3 categories. miR211 was found in 5 melanoma (36%). There was a significant downregulation of miR211 in the melanoma compared to the nevi (p=0,0219). MITF was found in 13 melanoma (93%). IGF2R and KCNMA1 were respectively found in 71% and 77% of the melanoma. There was no significant differential expression of MITF, KCNMA1 and IGF2R between the nevi and the melanoma as well as no association between miR211 expression and protein expression of two potential target genes Conclusions In vivo miR211 is significantly reduced in conjunctival melanoma compared to conjunctival nevi. No correlation between mir211 expression and two potential target genes KCNMA1 and IGF2R was observed.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Characterization of lower-limbs inter-segment coordination during the take-off extension in ski jumping.

Take-off, the most important phase in ski jumping, has been primarily studied in terms of spatio-temporal parameters; little is known about its motor control aspects. This study aims to assess the inter-segment coordination of the shank-thigh and thigh-sacrum pairs using the continuous relative phase (CRP). In total 87 jumps were recorded from 33 athletes with an inertial sensor-based system. The CRP curves indicated that the thighs lead the shanks during the first part of take-off extension and that the shanks rotated faster at the take-off extension end. The thighs and sacrum first rotated synchronously, with the sacrum then taking lead, with finally the thighs rotating faster. Five characteristic features were extracted from the CRP and their relationship with jump length was tested. Three features of the shank-thigh pair and one of the thigh-sacrum pair reported a significant association with jump length. It was observed that athletes who achieved longer jumps had their thighs leading their shanks during a longer time, with these athletes also having a more symmetric movement between thighs and sacrum. This study shows that inter-segment coordination during the take-off extension is related to performance and further studies are necessary to contrast its importance with other ski jumping aspects.

Urinary Sphincter Based on Electronics and Artificial Muscles

Objectives: The AMS 800 is the current artifi cial urinary sphincter (AUS) forincontinence due to intrinsic sphincter defi ciency. Despite good clinical results,technical failures inherent to the hydraulic mechanism or urethral ischemicinjury contribute to revisions up to 60%. We are developing an electronic AUS,called ARTUS to overcome the rigors of AMS. The objective of this study wasto evaluate the technical effi cacy and tissue tolerance of the ARTUS systemin an animal model.Methods: The ARTUS is composed by three parts: thecontractile unit, a series of rings and an integrated microprocessor. The contractileunit is made of Nitinol fi bers. The rings are placed around the urethrato control the fl ow of urine by squeezing the urethra. They work in a sequentialalternative mode and are controlled by a microprocessor. In the fi rst phase athree-rings device was used while in the second phase a two-rings ARTUS wasused. The device was implanted in 14 sheep divided in two groups of six andeight animals for study purpose. The fi rst group aimed at bladder leak pointpressure (BLPP) measurement and validation of the animal model; the secondgroup aimed at verifying midterm tissue tolerance by explants at twelve weeks.General animal tolerance was also evaluated.Results: The ARTUS systemimplantation was uneventful. When the system was activated, the BLPP wasmeasured at 1.038 ± 0.044 bar (mean ± SD). Urethral tissue analysis did notshow signifi cant morphological changes. No infection and no sign of discomfortwere noted in animals at 12 weeks.Conclusions: The ARTUS proved to beeffective in continence achievement in this study. Histological results supportour idea that a sequential alternative mode can avoid urethral atrophy andischemia. Further technical developments are needed to verify long-termoutcome and permit human use.

The squeeze cell hypothesis for the activation of jasmonate synthesis in response to wounding

Jasmonates are lipid mediators that control defence gene expression in response to wounding and other environmental stresses. These small molecules can accumulate at distances up to several cm from sites of damage and this is likely to involve cell-to-cell jasmonate transport. Also, and independently of jasmonate synthesis, transport and perception, different long-distance wound signals that stimulate distal jasmonate synthesis are propagated at apparent speeds of several cmmin(-1) to tissues distal to wounds in a mechanism that involves clade 3 GLUTAMATE RECEPTOR-LIKE (GLR) genes. A search for jasmonate synthesis enzymes that might decode these signals revealed LOX6, a lipoxygenase that is necessary for much of the rapid accumulation of jasmonic acid at sites distal to wounds. Intriguingly, the LOX6 promoter is expressed in a distinct niche of cells that are adjacent to mature xylem vessels, a location that would make these contact cells sensitive to the release of xylem water column tension upon wounding. We propose a model in which rapid axial changes in xylem hydrostatic pressure caused by wounding travel through the vasculature and lead to slower, radially dispersed pressure changes that act in a clade 3 GLR-dependent mechanism to promote distal jasmonate synthesis.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Taxonomic status and origin of the shrews (Soricidae) from the Canary islands inferred from a mtDNA comparison with the European Crocidura species.

European island shrews are either relicts of the endemic Pleistocene fauna, e.g.,. Crocidura zimmermanni, or were introduced from continental source populations. In order to clarify the taxonomic status and the origin of the two shrew species from the Canary islands, a 981bp fragment of cytochrome b gene was investigated in all European Crocidura species and compared with the Canary shrew (Crocidura canariensis) and the Osorio shrew (Crocidura osorio). The first shares its karyotype with the Sicilian shrew Crocidura sicula (2N=36), the second with the Greater white-toothed shrew Crocidura russula (2N=42), suggesting possible sister species relationships. Results confirm the monophyly of taxa sharing the same karyotype. Genetic distances between C. sicula and C. canariensis suggest a separation since 5 Myr. The first was probably isolated from the North African ancestor after the Messinian desiccation; the second arrived on the Canary islands by natural jump dispersal. Within the 2N=42 cluster, a first split separated an Eastern line (Tunisia) from a western line (Morocco/Europe) of C. russula. C. osorio clusters together with C. russula from Spain, indicating conspecificy. This suggests a recent introduction from Spain by human.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.