On some numerical aspects of an active strain model in cardiac mechanics

We are interested in the development, implementation and testing of an orthotropic model for cardiac contraction based on an active strain decomposition. Our model addresses the coupling of a transversely isotropic mechanical description at the cell level, with an orthotropic constitutive law for incompressible tissue at the macroscopic level. The main differences with the active stress model are addressed in detail, and a finite element discretization using Taylor-Hood and MINI elements is proposed and illustrated with numerical examples.

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.

Ultrastructural and morphometrical analyses of the brown adipocyte nucleus in a hibernating dormouse.

The fine morphology, size, and perichromatin granule frequency were analysed in brown adipocyte nuclei from hibernating, arousing, and euthermic dormice, Muscardinus avellanarius. Unusual nuclear structural constituents such as nuclear amorphous bodies, coiled body-like constituents and bundles of nucleoplasmic filaments were described as typical of hibernating nuclei. Morphometrical findings showed significant difference in total nuclear and nucleolar size in the three physiological conditions investigated as well as decreasing frequency of perichromatin granules in nuclei of hibernating to arousing to euthermic animals. A possible involvement of these granules in the intranuclear transport or storage of pre-mRNA is discussed in the context of other experimental evidence.

Systematic review of population pharmacokinetic analyses of imatinib and relationships with treatment outcomes

Several population pharmacokinetic (PPK) analyses of the anticancer drug imatinib have been performed to investigate different patient populations and covariate effects. The present analysis offers a systematic qualitative and quantitative summary and comparison of those. Its primary objective was to provide useful information for evaluating the expectedness of imatinib plasma concentration measurements in the frame of therapeutic drug monitoring. The secondary objective was to review clinically important concentration-effect relationships to provide help in evaluating the potential suitability of plasma concentration values. Nine PPK models describing total imatinib plasma concentration were identified. Parameter estimates were standardized to common covariate values whenever possible. Predicted median exposure (Cmin) was derived by simulations and ranged between models from 555 to 1388 ng/mL (grand median: 870 ng/mL and interquartile "reference" range: 520-1390 ng/mL). Covariates of potential clinical importance (up to 30% change in pharmacokinetic predicted by at least 1 model) included body weight, albumin, α1 acid glycoprotein, and white blood cell count. Various other covariates were included but were statistically not significant or seemed clinically less important or physiologically controversial. Concentration-response relationships had more importance below the average reference range and concentration-toxicity relationships above. Therapeutic drug monitoring-guided dosage adjustment seems justified for imatinib, but a formal predictive therapeutic range remains difficult to propose in the absence of prospective target concentration intervention trials. To evaluate the expectedness of a drug concentration measurement in practice, this review allows comparison of the measurement either to the average reference range or to a specific range accounting for individual patient characteristics. For future research, external PPK model validation or meta-model development should be considered.

Numerical analysis of the impact of charcoal production on soil hydrological behavior, runoff response and erosion susceptibility

The impact of charcoal production on soil hydraulic properties, runoff response and erosion susceptibility were studied in both field and simulation experiments. Core and composite samples, from 12 randomly selected sites within the catchment of Kotokosu were taken from the 0-10 cm layer of a charcoal site soil (CSS) and adjacent field soils (AFS). These samples were used to determine saturated hydraulic conductivity (Ksat), bulk density, total porosity, soil texture and color. Infiltration, surface albedo and soil surface temperature were also measured in both CSS and AFS. Measured properties were used as entries in a rainfall runoff simulation experiment on a smooth (5 % slope) plot of 25 x 25 m grids with 10 cm resolutions. Typical rainfall intensities of the study watershed (high, moderate and low) were applied to five different combinations of Ks distributions that could be expected in this landscape. The results showed significantly (p < 0.01) higher flow characteristics of the soil under charcoal kilns (increase of 88 %). Infiltration was enhanced and runoff volume reduced significantly. The results showed runoff reduction of about 37 and 18 %, and runoff coefficient ranging from 0.47-0.75 and 0.04-0.39 or simulation based on high (200 mm h-1) and moderate (100 mm h-1) rainfall events over the CSS and AFS areas, respectively. Other potential impacts of charcoal production on watershed hydrology were described. The results presented, together with watershed measurements, when available, are expected to enhance understanding of the hydrological responses of ecosystems to indiscriminate charcoal production and related activities in this region.

Expression analyses of three members of the AtPHO1 family reveal differential interactions between signaling pathways involved in phosphate deficiency and the responses to auxin, cytokinin, and abscisic acid.

The PHO1 protein is involved in loading inorganic phosphate (Pi) to the root xylem. Ten genes homologous to AtPHO1 are present in the Arabidopsis thaliana (L.) Heyn genome. From this gene family, transcript levels of only AtPHO1, AtPHO1;H1 and AtPHO1;H10 were increased by Pi-deficiency. While the up-regulation of AtPHO1;H1 and AtPHO1;H10 by Pi deficiency followed the same rapid kinetics and was dependent on the PHR1 transcription factor, phosphite only strongly suppressed the expression of AtPHO1;H1 and had a minor effect on AtPHO1;H10. Addition of sucrose was found to increase transcript levels of both AtPHO1 and AtPHO1;H1 in Pi-sufficient or Pi-deficient plants, but to suppress AtPHO1:H10 under the same conditions. Treatments of plants with auxin or cytokinin had contrasting effect depending on the gene and on the Pi status of the plants. Thus, while both hormones down-regulated expression of AtPHO1 independently of the plant Pi status, auxin and cytokinin up-regulated AtPHO1;H1 and AtPHO1;H10 expression in Pi-sufficient plants and down-regulated expression in Pi-deficient plants. Treatments with abscisic acid inhibited AtPHO1 and AtPHO1;H1 expression in both Pi-sufficient and Pi-deficient plants, but increased AtPHO1;H10 expression under the same conditions. The inhibition of expression by abscisic acid of AtPHO1 and AtPHO1;H1, and of the Pi-starvation responsive genes AtPHT1;1 and AtIPS1, was dependant on the ABI1 type 2C protein phosphatase. These results reveal that various levels of cross talk between the signal transduction pathways to Pi, sucrose and phytohormones are involved in the regulation of expression of the three AtPHO1 homologues.

How to assign a likelihood ratio in a footwear mark case: an analysis and discussion in the light of R v T

This article analyses and discusses issues that pertain to the choice of relevant databases for assigning values to the components of evaluative likelihood ratio procedures at source level. Although several formal likelihood ratio developments currently exist, both case practitioners and recipients of expert information (such as judiciary) may be reluctant to consider them as a framework for evaluating scientific evidence in context. The recent ruling R v T and ensuing discussions in many forums provide illustrative examples for this. In particular, it is often felt that likelihood ratio-based reasoning amounts to an application that requires extensive quantitative information along with means for dealing with technicalities related to the algebraic formulation of these approaches. With regard to this objection, this article proposes two distinct discussions. In a first part, it is argued that, from a methodological point of view, there are additional levels of qualitative evaluation that are worth considering prior to focusing on particular numerical probability assignments. Analyses will be proposed that intend to show that, under certain assumptions, relative numerical values, as opposed to absolute values, may be sufficient to characterize a likelihood ratio for practical and pragmatic purposes. The feasibility of such qualitative considerations points out that the availability of hard numerical data is not a necessary requirement for implementing a likelihood ratio approach in practice. It is further argued that, even if numerical evaluations can be made, qualitative considerations may be valuable because they can further the understanding of the logical underpinnings of an assessment. In a second part, the article will draw a parallel to R v T by concentrating on a practical footwear mark case received at the authors' institute. This case will serve the purpose of exemplifying the possible usage of data from various sources in casework and help to discuss the difficulty associated with reconciling the depth of theoretical likelihood ratio developments and limitations in the degree to which these developments can actually be applied in practice.

There and back again? Combining habitat suitability modelling and connectivity analyses to assess a potential return of the otter to Switzerland

After a steady decline in the early 20th century, several terrestrial carnivore species have recently recovered in Western Europe, either through reintroductions or natural recolonization. Because of the large space requirements of these species and potential conflicts with human activities, ensuring their recovery requires the implementation of conservation and management measures that address the environmental, landscape and social dimensions of the problem. Few examples exist of such integrated management. Taking the case of the otter (Lutra lutra) in Switzerland, we propose a multi-step approach that allows to (1) identify areas with potentially suitable habitat, (2) evaluate their connectivity, (3) verify the potentiality of the species recolonization from populations in neighbouring countries. We showed that even though suitable habitat is available for the species and the level of structural connectivity within Switzerland is satisfactory, the level of connectivity with neighbouring populations is crucial to prioritize strategies that favour the species recovery in the field. This research is the first example integrating habitat suitability and connectivity assessment at different scales with other factors in a multi-step assessment for species recovery.

Numerical algorithm for Ginzburg-Landau equations with multiplicative noise: Application to domain growth

We consider stochastic partial differential equations with multiplicative noise. We derive an algorithm for the computer simulation of these equations. The algorithm is applied to study domain growth of a model with a conserved order parameter. The numerical results corroborate previous analytical predictions obtained by linear analysis.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).