Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Improved statistical analysis of low abundance phenomena in bimodal bacterial populations.

Accurate detection of subpopulation size determinations in bimodal populations remains problematic yet it represents a powerful way by which cellular heterogeneity under different environmental conditions can be compared. So far, most studies have relied on qualitative descriptions of population distribution patterns, on population-independent descriptors, or on arbitrary placement of thresholds distinguishing biological ON from OFF states. We found that all these methods fall short of accurately describing small population sizes in bimodal populations. Here we propose a simple, statistics-based method for the analysis of small subpopulation sizes for use in the free software environment R and test this method on real as well as simulated data. Four so-called population splitting methods were designed with different algorithms that can estimate subpopulation sizes from bimodal populations. All four methods proved more precise than previously used methods when analyzing subpopulation sizes of transfer competent cells arising in populations of the bacterium Pseudomonas knackmussii B13. The methods' resolving powers were further explored by bootstrapping and simulations. Two of the methods were not severely limited by the proportions of subpopulations they could estimate correctly, but the two others only allowed accurate subpopulation quantification when this amounted to less than 25% of the total population. In contrast, only one method was still sufficiently accurate with subpopulations smaller than 1% of the total population. This study proposes a number of rational approximations to quantifying small subpopulations and offers an easy-to-use protocol for their implementation in the open source statistical software environment R.

Simulations of action of DNA topoisomerases to investigate boundaries and shapes of spaces of knots.

The configuration space available to randomly cyclized polymers is divided into subspaces accessible to individual knot types. A phantom chain utilized in numerical simulations of polymers can explore all subspaces, whereas a real closed chain forming a figure-of-eight knot, for example, is confined to a subspace corresponding to this knot type only. One can conceptually compare the assembly of configuration spaces of various knot types to a complex foam where individual cells delimit the configuration space available to a given knot type. Neighboring cells in the foam harbor knots that can be converted into each other by just one intersegmental passage. Such a segment-segment passage occurring at the level of knotted configurations corresponds to a passage through the interface between neighboring cells in the foamy knot space. Using a DNA topoisomerase-inspired simulation approach we characterize here the effective interface area between neighboring knot spaces as well as the surface-to-volume ratio of individual knot spaces. These results provide a reference system required for better understanding mechanisms of action of various DNA topoisomerases.

CT-guided sacroiliac joint injection: Easy or difficult? : 53

Purpose: Fluoroscopy-guided sacroiliac joint (SIJ) injections are technically difficult to perform because of the complex anatomy with helicoidal conformation of the joint. Our study describes the procedure of CT-guided SIJ injection, its feasibility and its rate of success. Methods and materials: Retrospective study included 46 consecutive patients. The procedure was performed by 3 MSK radiologists and consisted in a puncture with a posterior approach in the inferior articular part of SIJ, then in an injection of iodinated contrast agent (1ml) with CT control of SIJ space opacification and finally in an injection of slowacting corticosteroid. The SIJ approach was noticed as correct if there was an inferior articular puncture and if the needle was in the articular space, and as impossible if there was ankylosis or osteophytosis. The study was divided in two successive periods: period 1 (4 first months) and period 2 (12 last months). Results: SIJ opacification was successful in 57% (26/46). We observed a learning curve: opacification was succeeded in 66% (23/35) and there was incorrect approach in 9% (3/35) during period 2 versus respectively 27% (3/11) and 45% (5/11) during period 1. Causes of failure were incorrect approach in 40% (6/20 too low, 2/20 too high), impossible approach in 30% (6/20) and unexplained in 30% (6/20). Mean duration of procedure was about 28 minutes. No complication occurred. Conclusion: CT guided SIJ injection is safe and successful in 66% after a training period. The success depends on SIJ correct approach and also on anatomical lesions.

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.

Simulations of electrophoretic collisions of DNA knots with gel obstacles

Gel electrophoresis can be used to separate nicked circular DNA molecules of equal length but forming different knot types. At low electric fields, complex knots drift faster than simpler knots. However, at high electric field the opposite is the case and simpler knots migrate faster than more complex knots. Using Monte Carlo simulations we investigate the reasons of this reversal of relative order of electrophoretic mobility of DNA molecules forming different knot types. We observe that at high electric fields the simulated knotted molecules tend to hang over the gel fibres and require passing over a substantial energy barrier to slip over the impeding gel fibre. At low electric field the interactions of drifting molecules with the gel fibres are weak and there are no significant energy barriers that oppose the detachment of knotted molecules from transverse gel fibres.

Simulations in Evolution. III. Randomness as a Generator of Opportunities.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous reports presented a histogram model to simulate the evolution of populations of individuals classified into bins according to an unspecified, quantifiable phenotypic character, and whose number in each bin changed generation after generation under the influence of fitness, while the total population was maintained constant. The histogram model also allowed Shannon entropy (SE) to be monitored continuously as the information content of the total population decreased or increased. Here, a simple Perl (Practical Extraction and Reporting Language) application was developed to carry out these computations, with the critical feature of an added random factor in the percent of individuals whose offspring moved to a vicinal bin. The results of the simulations demonstrate that the random factor mimicking variation increased considerably the range of values covered by Shannon entropy, especially when the percentage of changed offspring was high. This increase in information content is interpreted as facilitated adaptability of the population.

A guide for easy- and difficult-to-treat hypertension.

Although the blood pressure (BP) of many patients can be controlled using standard combinations, treatment of hypertension frequently represents a clinical challenge to the primary care physician. This article will review best practices for managing patients with easy- and difficult-to-treat hypertension, including preferred antihypertensive combinations, optimizing adherence and persistence, recognizing white-coat hypertension, and intensifying therapy for treatment-resistant patients. Each physician must decide based on his or her own level of experience at what point a patient becomes too challenging and would benefit from referral to a hypertension specialist for more intensive management and to complete the exclusion of secondary forms of arterial hypertension. With intensive pharmacotherapy, many patients with difficult-to-treat hypertension can achieve BP control. If it fails, interventional strategies (e.g., renal denervation) are a valid option to get BP controlled.

Predicting the deleterious effects of mutation load in fragmented populations.

Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure (F(ST)) and within-deme heterozygosity (H(S)). A similar part of variance in median time to extinction was explained by a combination of local population size (N) and heterozygosity (H(S)). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.

Quantitative integration of high-resolution hydrogeophysical data through Monte-Carlo-type conditional simulations

Geophysical techniques can help to bridge the inherent gap with regard to spatial resolution and the range of coverage that plagues classical hydrological methods. This has lead to the emergence of the new and rapidly growing field of hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters and their inherent trade-off between resolution and range the fundamental usefulness of multi-method hydrogeophysical surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative integration of the resulting vast and diverse database in order to obtain a unified model of the probed subsurface region that is internally consistent with all available data. To address this problem, we have developed a strategy towards hydrogeophysical data integration based on Monte-Carlo-type conditional stochastic simulation that we consider to be particularly suitable for local-scale studies characterized by high-resolution and high-quality datasets. Monte-Carlo-based optimization techniques are flexible and versatile, allow for accounting for a wide variety of data and constraints of differing resolution and hardness and thus have the potential of providing, in a geostatistical sense, highly detailed and realistic models of the pertinent target parameter distributions. Compared to more conventional approaches of this kind, our approach provides significant advancements in the way that the larger-scale deterministic information resolved by the hydrogeophysical data can be accounted for, which represents an inherently problematic, and as of yet unresolved, aspect of Monte-Carlo-type conditional simulation techniques. We present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to corresponding field data collected at the Boise Hydrogeophysical Research Site near Boise, Idaho, USA.