Long-distance wolf recolonization of France and Switzerland inferred from non-invasive genetic sampling over a period of 10 years

In the early 1900s, the wolf (Canis lupus) was extirpated from France and Switzerland. There is growing evidence that the species is presently recolonizing these countries in the western Alps. By sequencing the mitochondrial DNA (mtDNA) control region of various samples mainly collected in the field (scats, hairs, regurgitates, blood or tissue; n = 292), we could (1) develop a non-invasive method enabling the unambiguous attribution of these samples to wolf, fox (Vulpes vulpes) or dog (Canis familiaris), among others; (2) demonstrate that Italian, French and Swiss wolves share the same mtDNA haplotype, a haplotype that has never been found in any other wolf population world-wide. Combined together, field and genetic data collected over 10 years corroborate the scenario of a natural expansion of wolves from the Italian source population. Furthermore, such a genetic approach is of conservation significance, since it has important consequences for management decisions. This first long-term report using non-invasive sampling demonstrates that long-distance dispersers are common, supporting the hypothesis that individuals may often attempt to colonize far from their native pack, even in the absence of suitable corridors across habitats characterized by intense human activities.

Cost-effectiveness of surgery in low- and middle-income countries: a systematic review.

BACKGROUND: There is increasing interest in provision of essential surgical care as part of public health policy in low- and middle-income countries (LMIC). Relatively simple interventions have been shown to prevent death and disability. We reviewed the published literature to examine the cost-effectiveness of simple surgical interventions which could be made available at any district hospital, and compared these to standard public health interventions. METHODS: PubMed and EMBASE were searched using single and combinations of the search terms "disability adjusted life year" (DALY), "quality adjusted life year," "cost-effectiveness," and "surgery." Articles were included if they detailed the cost-effectiveness of a surgical intervention of relevance to a LMIC, which could be made available at any district hospital. Suitable articles with both cost and effectiveness data were identified and, where possible, data were extrapolated to enable comparison across studies. RESULTS: Twenty-seven articles met our inclusion criteria, representing 64 LMIC over 16 years of study. Interventions that were found to be cost-effective included cataract surgery (cost/DALY averted range US$5.06-$106.00), elective inguinal hernia repair (cost/DALY averted range US$12.88-$78.18), male circumcision (cost/DALY averted range US$7.38-$319.29), emergency cesarean section (cost/DALY averted range US$18-$3,462.00), and cleft lip and palate repair (cost/DALY averted range US$15.44-$96.04). A small district hospital with basic surgical services was also found to be highly cost-effective (cost/DALY averted 1 US$0.93), as were larger hospitals offering emergency and trauma surgery (cost/DALY averted US$32.78-$223.00). This compares favorably with other standard public health interventions, such as oral rehydration therapy (US$1,062.00), vitamin A supplementation (US$6.00-$12.00), breast feeding promotion (US$930.00), and highly active anti-retroviral therapy for HIV (US$922.00). CONCLUSIONS: Simple surgical interventions that are life-saving and disability-preventing should be considered as part of public health policy in LMIC. We recommend an investment in surgical care and its integration with other public health measures at the district hospital level, rather than investment in single disease strategies.

Twofold cost of reproduction: an increase in parental effort leads to higher malarial parasitaemia and to a decrease in resistance to oxidative stress.

Parental effort is usually associated with high metabolism that could lead to an increase in the production of reactive oxidative species giving rise to oxidative stress. Since many antioxidants involved in the resistance to oxidative stress can also enhance immune function, an increase in parental effort may diminish the level of antioxidants otherwise involved in parasite resistance. In the present study, we performed brood size manipulation in a population of great tits (Parus major) to create different levels of parental effort. We measured resistance to oxidative stress and used a newly developed quantitative PCR assay to quantify malarial parasitaemia. We found that males with an enlarged brood had significantly higher level of malarial parasites and lower red blood cell resistance to free radicals than males rearing control and reduced broods. Brood size manipulation did not affect female parasitaemia, although females with an enlarged brood had lower red blood cell resistance than females with control and reduced broods. However, for both sexes, there was no relationship between the level of parasitaemia and resistance to oxidative stress, suggesting a twofold cost of reproduction. Our results thus suggest the presence of two proximate and independent mechanisms for the well-documented trade-off between current reproductive effort and parental survival.

The nine-year sustained cost-containment impact of swiss pilot physicians-pharmacists quality circles.

BACKGROUND: Six pioneer physicians-pharmacists quality circles (PPQCs) located in the Swiss canton of Fribourg (administratively corresponding to a state in the US) were under the responsibility of 6 trained community pharmacists moderating the prescribing process of 24 general practitioners (GPs). PPQCs are based on a multifaceted collaborative process mediated by community pharmacists for improving compliance with clinical guidelines within GPs' prescribing practices. OBJECTIVE: To assess, over a 9-year period (1999-2007), the cost-containment impact of the PPQCs. METHODS: The key elements of PPQCs are a structured continuous quality improvement and education process; local networking; feedback of comparative and detailed data regarding costs, drug choice, and frequency of prescribed drugs; and structured independent literature review for interdisciplinary continuing education. The data are issued from the community pharmacy invoices to the health insurance companies. The study analyzed the cost-containment impact of the PPQCs in comparison with GPs working in similar conditions of care without particular collaboration with pharmacists, the percentage of generic prescriptions for specific cardiovascular drug classes, and the percentage of drug costs or units prescribed for specific cardiovascular drugs. RESULTS: For the 9-year period, there was a 42% decrease in the drug costs in the PPQC group as compared to the control group, representing a $225,000 (USD) savings per GP only in 2007. These results are explained by better compliance with clinical and pharmacovigilance guidelines, larger distribution of generic drugs, a more balanced attitude toward marketing strategies, and interdisciplinary continuing education on the rational use of drugs. CONCLUSIONS: The PPQC work process has yielded sustainable results, such as significant cost savings, higher penetration of generics and reflection on patient safety, and the place of "new" drugs in therapy. The PPQCs may also constitute a solid basis for implementing more comprehensive collaborative programs, such as medication reviews, adherence-enhancing interventions, or disease management approaches.

Cost evaluation and comparison of three decision strategies to revascularize : results of the suspected CAD protocol of the european CMR registry

Background: Cardiac magnetic resonance (CMR) is accepted as a method to assess suspected coronary artery disease (CAD). Nonetheless, invasive coronary angiography (CXA) combined or not with fractional flow reserve (FFR) remains the main diagnostic test to evaluate CAD. Little data exist on the economic impact of the use of these procedures in a population with a low to intermediate pre-test probability. Objective: To compare the costs of 3 decision strategies to revascularize a patient with suspected CAD: 1) strategy guided by CMR 2) hypothetical strategy guided by CXA-FFR, 3) hypothetical strategy guided by CXA alone.

Strategy selection: An introduction to the modeling challenge

Modeling the mechanisms that determine how humans and other agents choose among different behavioral and cognitive processes-be they strategies, routines, actions, or operators-represents a paramount theoretical stumbling block across disciplines, ranging from the cognitive and decision sciences to economics, biology, and machine learning. By using the cognitive and decision sciences as a case study, we provide an introduction to what is also known as the strategy selection problem. First, we explain why many researchers assume humans and other animals to come equipped with a repertoire of behavioral and cognitive processes. Second, we expose three descriptive, predictive, and prescriptive challenges that are common to all disciplines which aim to model the choice among these processes. Third, we give an overview of different approaches to strategy selection. These include cost‐benefit, ecological, learning, memory, unified, connectionist, sequential sampling, and maximization approaches. We conclude by pointing to opportunities for future research and by stressing that the selection problem is far from being resolved.

Minimally invasive fixation versus conservative treatment of undisplaced scaphoid fractures: a cost-effectiveness study.

This study compares the direct and indirect costs of conservative and minimally invasive treatment for undisplaced scaphoid fractures. Costs data concerning groups of non-operated and operated patients were analysed. Direct costs were higher in operated patients. Although highly variable, indirect costs were significantly smaller in operated patients and the total costs were higher in non-operated patients. In conclusion, operative treatment of scaphoid fractures is initially more expensive than conservative treatment but markedly decreases the work compensation costs.

Assessment of diesel exhaust particulate exposure and surface characteristics in association with levels of oxidative stress biomarkers

Exposure to PM10 and PM2.5 (particulate matter with aerodynamic diameter smaller than 10 μm and 2.5 μm, respectively) is associated with a range of adverse health effects, including cancer, pulmonary and cardiovascular diseases. Surface characteristics (chemical reactivity, surface area) are considered of prime importance to understand the mechanisms which lead to harmful effects. A hypothetical mechanism to explain these adverse effects is the ability of components (organics, metal ions) adsorbed on these particles to generate Reactive Oxygen Species (ROS), and thereby to cause oxidative stress in biological systems (Donaldson et al., 2003). ROS can attack almost any cellular structure, like DNA or cellular membrane, leading to the formation of a wide variety of degradation products which can be used as a biomarker of oxidative stress. The aim of the present research project is to test whether there is a correlation between the exposure to Diesel Exhaust Particulate (DEP) and the oxidative stress status. For that purpose, a survey has been conducted in real occupational situations where workers were exposed to DEP (bus depots). Different exposure variables have been considered: - particulate number, size distribution and surface area (SMPS); - particulate mass - PM2.5 and PM4 (gravimetry); - elemental and organic carbon (coulometry); - total adsorbed heavy metals - iron, copper, manganese (atomic adsorption); - surface functional groups present on aerosols (Knudsen flow reactor). Several biomarkers of oxidative stress (8-hydroxy-2'-deoxyguanosine and several aldehydes) have been determined either in urine or serum of volunteers. Results obtained during the sampling campaign in several bus depots indicated that the occupational exposure to particulates in these places was rather low (40-50 μg/m3 for PM4). Bimodal size distributions were generally observed (5 μm and <1 μm). Surface characteristics of PM4 varied strongly, depending on the bus depot. They were usually characterized by high carbonyl and low acidic sites content. Among the different biomarkers which have been analyzed within the framework of this study, mean urinary levels of 8-hydroxy-2'-deoxyguanosine increased significantly (p<0.05) during two consecutive days of exposure for non-smoker workers. On the other hand, no statistically significant differences were observed for serum levels of hexanal, nonanal and 4- hydroxy-nonenal (p>0.05). Biomarkers levels will be compared to exposure variables to gain a better understanding of the relation between the particulate characteristics and the formation of ROS by-products. This project is financed by the Swiss State Secretariat for Education and Research. It is conducted within the framework of the COST Action 633 "Particulate Matter - Properties Related to Health Effects".

Robust parametric indirect estimates of the expected cost of a hospital stay with covariates and censored data.

We consider the problem of estimating the mean hospital cost of stays of a class of patients (e.g., a diagnosis-related group) as a function of patient characteristics. The statistical analysis is complicated by the asymmetry of the cost distribution, the possibility of censoring on the cost variable, and the occurrence of outliers. These problems have often been treated separately in the literature, and a method offering a joint solution to all of them is still missing. Indirect procedures have been proposed, combining an estimate of the duration distribution with an estimate of the conditional cost for a given duration. We propose a parametric version of this approach, allowing for asymmetry and censoring in the cost distribution and providing a mean cost estimator that is robust in the presence of extreme values. In addition, the new method takes covariate information into account.

Fitness cost and impaired survival in penicillin-resistant Streptococcus gordonii isolates selected in the laboratory.

Cycling of Streptococcus gordonii (115 times) with penicillin resulted in a MIC increase of more than 100-fold, from 0.008 to 2 microg/ml. The 2-microg/ml MIC maximum was already reached after 36 passages but resulted in impaired fitness. Although the MIC did not increase further, fitness was partially recovered during the 79 additional cycles.

Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Low-cost MR-compatible moving heart phantom.

For the development and evaluation of cardiac magnetic resonance (MR) imaging sequences and methodologies, the availability of a periodically moving phantom to model respiratory and cardiac motion would be of substantial benefit. Given the specific physical boundary conditions in an MR environment, the choice of materials and power source of such phantoms is heavily restricted. Sophisticated commercial solutions are available; however, they are often relatively costly and user-specific modifications may not easily be implemented. We therefore sought to construct a low-cost MR-compatible motion phantom that could be easily reproduced and had design flexibility. A commercially available K'NEX construction set (Hyper Space Training Tower, K'NEX Industries, Inc., Hatfield, PA) was used to construct a periodically moving phantom head. The phantom head performs a translation with a superimposed rotation, driven by a motor over a 2-m rigid rod. To synchronize the MR data acquisition with phantom motion (without introducing radiofrequency-related image artifacts), a fiberoptic control unit generates periodic trigger pulses synchronized to the phantom motion. Total material costs of the phantom are US$ < 200.00, and a total of 80 man-hours were required to design and construct the original phantom. With schematics of the present solution, the phantom reproduction may be achieved in approximately 15 man-hours. The presented MR-compatible periodically moving phantom can easily be reproduced, and user-specific modifications may be implemented. Such an approach allows a detailed investigation of motion-related phenomena in MR images.