A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.

Propagation dans l'environnement des microorganismes et autres particules organiques issus des élevages intensifs et conséquences pour la santé des riverains

Les élevages d'animaux de rente hébergent de plus en plus de bêtes. Cette situation génère une accumulation de poussière organique, constituée de particules inertes et de microorganismes, issus de la nourriture, de la litière, des matières fécales, des pellicules de la peau, des poils, etc. L'activité des animaux et l'activité professionnelle favorisent une remise en suspension de cette poussière, qui peut se propager à l'extérieur. Ces émissions de particules organiques dans l'environnement soulèvent des inquiétudes pour la santé des riverains. Ces craintes sont légitimes, puisque les problèmes respiratoires, allergiques ou toxiques sont bien connus chez les travailleurs agricoles exposés à de fortes doses de poussières organiques. Un autre risque sanitaire lié aux élevages intensifs d'animaux est la dissémination de bactéries résistantes aux antibiotiques dans l'environnement avec, pour éventuelle conséquence, une transmission de ces souches aux personnes résidant à proximité. Cette problématique est bien connue dans les élevages de porcs fréquemment colonisés par des SARM (Staphylococcus aureus résistant à la méticilline), qui sont transmis aux éleveurs. Les deux études analysées ci-dessous ont investigué cette problématique de dissémination des particules organiques dans l'environnement et les conséquences sur la santé des riverains. La première a étudié le lien entre le fait de résider à proximité de fermes d'élevage d'animaux et la prévalence de maladies respiratoires. La deuxième a étudié le risque de colonisation nasale par des SARM dans une population de vétérans vivant à proximité d'élevages intensifs de porcs.

Do seismic waves sense fracture connectivity?

A defining characteristic of fractured rocks is their very high level of seismic attenuation, which so far has been assumed to be mainly due to wave-induced fluid flow (WIFF) between the fractures and the pore space of the embedding matrix. Using oscillatory compressibility simulations based on the quasi-static poroelastic equations, we show that another important, and as of yet undocumented, manifestation of WIFF is at play in the presence of fracture connectivity. This additional energy loss is predominantly due to fluid flow within the connected fractures and is sensitive to their lengths, permeabilities, and intersection angles. Correspondingly, it contains key information on the governing hydraulic properties of fractured rock masses and hence should be accounted for whenever realistic seismic models of such media are needed.

Roles of astrocytic sodium dynamics in the neurometabolic coupling : a study by cellular imaging

RESUME LARGE PUBLIC Le système nerveux central est principalement composé de deux types de cellules :les neurones et les cellules gliales. Ces dernières, bien que l'emportant en nombre sur les neurones, ont longtemps été considérées comme des cellules sans intérêts par les neuroscientifiques. Hors, les connaissances modernes à leurs sujets indiquent qu'elles participent à la plupart des tâches physiologiques du cerveau. Plus particulièrement, elles prennent part aux processus énergétiques cérébraux. Ceux-ci, en plus d'être vitaux, sont particulièrement intrigants puisque le cerveau représente seulement 2 % de la masse corporelle mais consomme environ 25 % du glucose (substrat énergétique) corporel. Les astrocytes, un type de cellules gliales, jouent un rôle primordial dans cette formidable utilisation de glucose par le cerveau. En effet, l'activité neuronale (transmission de l'influx nerveux) est accompagnée d'une augmentation de la capture de glucose, issu de la circulation sanguine, par les astrocytes. Ce phénomène est appelé le «couplage neurométabolique » entre neurones et astrocytes. L'ion sodium fait partie des mécanismes cellulaires entrant en fonction lors de ces processus. Ainsi, dans le cadre de cette thèse, les aspects dynamiques de la régulation du sodium astrocytaire et leurs implications dans le couplage neurométabolique ont été étudiés par des techniques d'imagerie cellulaires. Ces études ont démontré que les mitochondries, machineries cellulaires convertissant l'énergie contenue dans le glucose, participent à la régulation du sodium astrocytaire. De plus, ce travail de thèse a permis de découvrir que les astrocytes sont capables de se transmettre, sous forme de vagues de sodium se propageant de cellules en cellules, un message donnant l'ordre d'accroître leur consommation d'énergie. Cette voie de signalisation leur permettrait de fournir de l'énergie aux neurones suite à leur activation. RESUME Le glutamate libéré dans la fente synaptique pendant l'activité neuronale, est éliminé par les astrocytes environnants. Le glutamate est co-transporté avec des ions sodiques, induisant une augmentation intracellulaire de sodium (Na+i) dans les astrocytes. Cette élévation de Na+i déclenche une cascade de mécanismes moléculaires qui aboutissent à la production de substrats énergétiques pouvant être utilisés par les neurones. Durant cette thèse, la mesure simultanée du sodium mitochondrial (Na+mit) et cytosolique par des techniques d'imagerie utilisant des sondes fluorescentes spécifiques, a indiqué que les variations de Na+i induites par le transport du glutamate sont transmises aux mitochondries. De plus, les voies d'entrée et de sortie du sodium mitochondrial ont été identifiées. L'échangeur de Na+ et de Ca2+ mitochondrial semble jouer un rôle primordial dans l'influx de Na+mit, alors que l'efflux de Na+mit est pris en charge par l'échangeur de Na+ et de H+ mitochondrial. L'étude du Na+mit a nécessité l'utilisation d'un système de photoactivation. Les sources de lumière ultraviolette (UV) classiques utilisées à cet effet (lasers, lampes à flash) ayant plusieurs désavantages, une alternative efficace et peu coûteuse a été développée. Il s'agit d'un système compact utilisant une diode électroluminescente (LED) à haute puissance et de longueur d'onde de 365nm. En plus de leurs rôles dans le couplage neurométabolique, les astrocytes participent à la signalisation multicellulaire en transmettant des vagues intercellulaires de calcium. Ce travail de thèse démontre également que des vagues intercellulaires de sodium peuvent être évoquées en parallèle à ces vagues calciques. Le glutamate, suite à sa libération par un mécanisme dépendent du calcium, est réabsorbé par les transporteurs au glutamate. Ce mécanisme a pour conséquence la génération de vagues sodiques se propageant de cellules en cellules. De plus, ces vagues sodiques sont corrélées spatialement avec une consommation accrue de glucose par les astrocytes. En conclusion, ce travail de thèse a permis de montrer que le signal sodique astrocytaire, déclenché en réponse au glutamate, se propage à la fois de façon intracellulaire aux mitochondries et de façon intercellulaire. Ces résultats suggèrent que les astrocytes fonctionnent comme un réseau de cellules nécessaire au couplage énergétique concerté entre neurones et astrocytes et que le sodium est un élément clé dans les mécanismes de signalisations cellulaires sous-jacents. SUMMARY Glutamate, released in the synaptic cleft during neuronal activity, is removed by surrounding astrocytes. Glutamate is taken-up with Na+ ions by specific transporters, inducing an intracellular Na+ (Na+i) elevation in astrocytes which triggers a cascade of molecular mechanisms that provides metabolic substrates to neurons. Thus, astrocytic Na+i homeostasis represents a key component of the so-called neurometabolic coupling. In this context, the first part of this thesis work was aimed at investigating whether cytosolic Na+ changes are transmitted to mitochondria, which could therefore influence their function and contribute to the overall intracellular Na+ regulation. Simultaneous monitoring of both mitochondrial Na+ (Na+mit) and cytosolic Na+ changes with fluorescent dyes revealed that glutamate-evoked cytosolic Na+ elevations are indeed transmitted to mitochondria. The mitochondrial Na+/Ca2+ exchangers have a prominent role in the regulation of Na+mit influx pathway, and Na+mit extrusion appears to be mediated by Na+/H+ exchangers. To demonstrate the implication of Na+/Ca2+ exchangers, this study has required the technical development of an UV-flash photolysis system. Because light sources for flash photolysis have to be powerful and in the near UV range, the use of UV lasers or flash lamps is usually required. As an alternative to these UV sources that have several drawbaks, we developped a compact, efficient and lowcost flash photolysis system which employs a high power 365nm light emitting diode. In addition to their role in neurometabolic coupling, astrocytes participate in multicellular signaling by transmitting intercellular Ca2+ waves. The third part of this thesis show that intercellular Na+ waves can be evoked in parallel to Ca2+ waves. Glutamate released by a Ca2+ wave-dependent mechanism is taken up by glutamate transporters, resulting in a regenerative propagation of cytosolic Na+ increases. Na+ waves in turn lead to a spatially correlated increase in glucose uptake. In conclusion, the present thesis demonstrates that glutamate-induced Na+ changes occurring in the cytosol of astrocytes propagate to both the mitochondrial matrix and the astrocytic network. These results furthermore support the view that astrocytic Na+ is a signal coupled to the brain energy metabolism.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

Real-time, in vivo intracellular recordings of caterpillar-induced depolarization waves in sieve elements using aphid electrodes.

Plants propagate electrical signals in response to artificial wounding. However, little is known about the electrophysiological responses of the phloem to wounding, and whether natural damaging stimuli induce propagating electrical signals in this tissue. Here, we used living aphids and the direct current (DC) version of the electrical penetration graph (EPG) to detect changes in the membrane potential of Arabidopsis sieve elements (SEs) during caterpillar wounding. Feeding wounds in the lamina induced fast depolarization waves in the affected leaf, rising to maximum amplitude (c. 60 mV) within 2 s. Major damage to the midvein induced fast and slow depolarization waves in unwounded neighbor leaves, but only slow depolarization waves in non-neighbor leaves. The slow depolarization waves rose to maximum amplitude (c. 30 mV) within 14 s. Expression of a jasmonate-responsive gene was detected in leaves in which SEs displayed fast depolarization waves. No electrical signals were detected in SEs of unwounded neighbor leaves of plants with suppressed expression of GLR3.3 and GLR3.6. EPG applied as a novel approach to plant electrophysiology allows cell-specific, robust, real-time monitoring of early electrophysiological responses in plant cells to damage, and is potentially applicable to a broad range of plant-herbivore interactions.

Childhood epilepsy with neuropsychological regression and continuous spike waves during sleep: epilepsy surgery in a young adult.

We describe the case of a man with a history of complex partial seizures and severe language, cognitive and behavioural regression during early childhood (3.5 years), who underwent epilepsy surgery at the age of 25 years. His early epilepsy had clinical and electroencephalogram features of the syndromes of epilepsy with continuous spike waves during sleep and acquired epileptic aphasia (Landau-Kleffner syndrome), which we considered initially to be of idiopathic origin. Seizures recurred at 19 years and presurgical investigations at 25 years showed a lateral frontal epileptic focus with spread to Broca's area and the frontal orbital regions. Histopathology revealed a focal cortical dysplasia, not visible on magnetic resonance imaging. The prolonged but reversible early regression and the residual neuropsychological disorders during adulthood were probably the result of an active left frontal epilepsy, which interfered with language and behaviour during development. Our findings raise the question of the role of focal cortical dysplasia as an aetiology in the syndromes of epilepsy with continuous spike waves during sleep and acquired epileptic aphasia.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Female sterility associated with increased clonal propagation suggests a unique combination of androdioecy and asexual reproduction in populations of Cardamine amara (Brassicaceae).

BACKGROUND AND AIMS: The coexistence of hermaphrodites and female-sterile individuals, or androdioecy, has been documented in only a handful of plants and animals. This study reports its existence in the plant species Cardamine amara (Brassicaceae), in which female-sterile individuals have shorter pistils than seed-producing hermaphrodites. METHODS: Morphological analysis, in situ manual pollination, microsatellite genotyping and differential gene expression analysis using Arabidopsis microarrays were used to delimit variation between female-sterile individuals and hermaphrodites. KEY RESULTS: Female sterility in C. amara appears to be caused by disrupted ovule development. It was associated with a 2.4- to 2.9-fold increase in clonal propagation. This made the pollen number of female-sterile genets more than double that of hermaphrodite genets, which fulfils a condition of co-existence predicted by simple androdioecy theories. When female-sterile individuals were observed in wild androdioecious populations, their ramet frequencies ranged from 5 to 54 %; however, their genet frequencies ranged from 11 to 29 %, which is consistent with the theoretically predicted upper limit of 50 %. CONCLUSIONS: The results suggest that a combination of sexual reproduction and increased asexual proliferation by female-sterile individuals probably explains the invasion and maintenance of female sterility in otherwise hermaphroditic populations. To our knowledge, this is the first report of the coexistence of female sterility and hermaphrodites in the Brassicaceae.