Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

What is "clinical data"? Why and how can they be collected during field surveys on medicinal plants?

ETHNOPHARMACOLOGICAL RELEVANCE: "Reverse pharmacology", also called "bedside-to-bench" or "field to pharmacy" approach, is a research process starting with documentation of clinical outcome as observed by patients with different therapeutic regimens. The treatment most significantly associated with cure is selected for future studies: first, clinical safety and efficacy; then in vivo and vitro studies. Some clinical data, i.e. details on patient status and progress, can be collected during ethnobotanical surveys; they will help clinical researchers and, once effectiveness and safety are established, will also help users of traditional medicine make safer and more effective choices. To gather clinical data successfully, ethnopharmacologists need to be backed by an appropriate team of specialists in medicine and epidemiology. Ethnopharmacologists can also gather important data on traditional medicine safety. MATERIALS AND METHODS: The first step is to create a consensus on the meaning of "clinical data", their interest and importance. An understanding of why "a cure is not a proof of effectiveness" is a starting point to avoid faulty interpretation of the clinical observations. RESULTS: Experience showed that, with the "bedside-to-bench" approach, a treatment derived from traditional recipe can be scientifically validated (in terms of safety and effectiveness) with a cost of less than a million euros, thus providing an end-product that is affordable, available and sustainable. CONCLUSIONS: With rigorous clinical study results, medicinal plant users gain the possibility to refine heath strategies. The field surveyor may gain a better relationship with the population, once she/he is seen as bringing information useful for the quality of care in the community.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

The Exotic Palaeo-Thethys Terrane in Central Iran: New geological data from Anarak, Jandaq and Posht-e-Badam areas

Résumé Le « terrane » d'Anarak-Jandak occupe une position géologique clé au nord-ouest du Microcontinent Centre-East Iranien (CE1M), connecté avec le Bloc du Grand Kavir et la ceinture métamorphique de Sanandaj-Sirjan. Nous discutons ici l'origine de ces différentes unités, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, pour conclure finalement de leur affinité paléotéthysienne. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur-Dévonien inférieur, pour se terminer au Trias par la collision des blocs Cimmériens dérivé du Gondwana avec le Bloc du Turan d'affinité asiatique (événement Eocimmérien). La plus importante unité métamorphique affleurant au sud-ouest de la région de Jandak-Anarak-Kaboudan est une épaisse séquence silicoclastique à grains fins contenant des blocs ophiolitiques (marginal-sea-type), et des associations basalte-gabbro à signatures géochimiques de type supra-subduction. Dans la région de Nakhlak, nous avons daté ces gabbros par la méthode U-Pb à 387f0.11 Ma ; les roches métamorphiques pélitiques ont donné des âges de refroidissement Ar-Ar pour la muscovite de 320 à 333 Ma. Ce complexe d'accrétion "varisque" a été métamorphisé dans le faciès schiste vert-amphibolite au cours de l'accrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma par la méthode U/Pb), qui affleure aujourd'hui à l'extrémité nord-ouest du terrane d'Anarak-Jandak . La subduction vers le nord de l'océan Paléotéthys depuis le Paléazoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de guyots (Anarak, Kaboudan, et Meraji Seamounts) et de hauts sous-marins, entrés en collision oblique avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries {âges Ar-Ar de 280 à 230 Ma). De plus, le magmatisme bimodal de Chah Gorbeh est caractérisé d'une part par des roches de type trondjémite-gabbros (262 Ma), d'autre part par des laves en coussin de type basaltes alcalins-rhyolites; ces roches magmatiques ont recoupé l'ophiolite d'Anarak lors de la mise en place de cette dernière dans la fosse interne de subduction. Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, i1 a été accrété le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge probable Triasique. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé dans les dépôts infra-arc Dévonien supérieur-Carbonifère de Godar-e-Siah, ainsi que dans la succession d'avant-arc de Nakhlak. Pendant l'intervalle Paléozoïque supérieur-Trias, la région de Jandak a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, elle-même comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites d'arc à collisionnel datés à 215±15 Ma. Dans la région de Yazd, témoin de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur; il en a été de même pour les dépôts de plate-forme Paléozoïque supérieur. L'érosion, qui dans ce dernier cas a atteint le Permien, pourrait être liée au bombement flexural de la marge passive. La collision finale n'a pas induit de déformations trop importantes, et se caractérise par la mise en place de nappes sur la marge passive. Cet événement est scellé par des dépôts molassique du Lias. D'un point de vue régional, la zone s'étendant actuellement de la Mer Noire au Pamir a été soumise à six épisodes d'extension-compression du Jurassique inférieur (début du l'ouverture en position arrière-arc de la Néotéthys) à l'Eocène moyen. Par exemple, le terrane d'AnarakJandak, probablement situé entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de sa patrie d'origine au début du Crétacé supérieur. Des preuves de cet événement se retrouvent dans les séries de plate-forme de Khur (préservation de séries syn-rift puis de marge passive). Les ophiolites de Nain et de Sabzevar sont de plus interprétée comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation par la plaque indienne de l'Eurasie a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent Iranien et de la formation du CEIM. Cette rotation est responsable du transport du terrane d'Anarak-Jandak vers sa position actuelle en Iran Central, et de la dislocation de Terranes de moindre importance, comme le bloc de Posht-e Badam. Depuis le Miocène supérieur, et à la suite de la collision entre l'Arabie et l'Iran, le ternane d'Anarak-Jandak a subi des déformations liées à l'activité d'une zone de cisaillement dextre parallèle à la suture du Zagros, à l'arrière de l'arc magmatique d'Uromieh-Dokhtar. Résumé large public Le Microcontinent Centre-Est Iranien occupe une position géologique clé au centre de l'Iran. Les différentes unités qui le composent, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, sont maintenant rajeunies et liés à la fermeture de l'océean Paléotéthys. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur à Dévonien inférieur, pour se terminer au Trias par la collision des- blocs Cimmériens, dérivés du Gondwana, avec le Bloc du Turan d'affinité asiatique. Dans la marge active asiatique de la Paléotéthys, nous avons daté les restes d'un océan marginal à 387±0.11 Ma. Ce complexe d'accrétion a été métamorphisé au cours de la réaccrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma), qui affleure aujourd'hui à l'extrémité nord-ouest du « terrane » d'Anarak-Jandak correspondant à la plus grande partie de la région étudiée. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé et daté du Dévonien supérieur-Carbonifère. Pendant l'intervalle Paléozoïque supérieur-Trias, la région a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites datés à 215±15 Ma. La subduction vers le nord de l'océan Paléotéthys depuis le Paléozoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de volcans sous-marins, entrés en collision avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries (280 à 230 Ma). Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, il a été mis en place le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge Triasique. Dans la région de Yazd, on trouve les témoins de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur, marqué par la flexuration de la marge passive lorsqu'elle rentra en collision avec la marge active asiatique. Cet événement est scellé par des dépôts molassique à charbon du Lias. Le «terrane» d'Anarak-Jandak, probablement situé à l'origine entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de cette région au début du Crétacé supérieur lors de l'ouverture d'un bassin d'arrière-arc, engendré, cette fois, par la subduction de l'océan Néotéthys situé au sud des blocs cimmériens. Des preuves de cet événement se retrouvent dans les séries syn-rift, puis de marge passive de Khour. Les ophiolites de Nain et de Sabzevar sont interprétées comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation de l'Eurasie par la plaque indienne a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent centre-Iranien. Cette rotation de près de 90° est responsable du transport du « terrane » d'Anarak-Jandak vers sa position actuelle. Abstract The Anarak-Jandaq terrane occupies a strategic geological situation at the north-western part of the Central-East Iranian Microcontinent (CEIM) and in connection with the Great Kavir Block and Sanandaj-Sirjan metamorphic belt. Our recent findings redefine the origin of these mentioned areas so far attributed to the Precambrian-Early Palaeozoic orogenic episodes, to be now directly related to the tectonic evolution of the Palaeo-Tethys Ocean, commenced by Late Ordovician-Early Devonian rifting events and terminated in the Triassic by the Eocimmerian tectonic event due to the collision of the Cimmerian blocks with the Asiatic Turan block. The most distributed metamorphic unit that is exposed from the south-west of Jandaq to the Anarak and Kaboudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea-basin ophiolitic blocks including basalt-gabbro association with supra-subduction-geochemical signature. These gabbros in the Nakhlak area were dated by U/Pb method at 387.6 ± 0.11 Ma and the metamorphic pelitic rocks yielded a range of 320 to 333 Ma muscovite-cooling ages based on 40Ar/39 Ar method. This "Variscan" accretionary complex was metamorphosed in greenschist-amphibolite facies during accretion to the Lower Cambrian Airekan granitic belt (549 ± 15 Ma by U/Pb method) that crops out at the northwestern edge of the Anarak-Jandaq terrane. Continued northward subduction of the Palaeo-Tethys Ocean during the entire Late Palaeozoic-Middle Triassic brought huge amount of oceanic material to the subduction zone. One chain of Carboniferous-Triassic oceanic rises and seamounts (the Anarak, Kaboudan, and Meraji Seamounts) obliquely collided with the accretionary wedge and created a mild HP metamorphic event (280-230 Ma based on 40Ar/39Ar results). Bimodal magmatism of the Chah Gorbeh area is characterized by a 262 Ma trondjemite-gabbro as well as pillow alkalibasalts-rhyolites which intruded the Anarak ophiolite when it was being emplaced within the inner-wall trench. The mainly Late Permian-Triassic Doshakh wedge was accreted along the continent and metamorphosed under lower greenschist facies and the probable Triassic Bayazeh flysch filled the foreland basin during the final closure. The Palaeo-Tethys magmatic arc products have been well preserved in the Late Devonian-Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. During the Late Palaeozoic-Triassic times, the Jandaq area has been affected by back-arc extension and probably the Arusan ophiolitic belt is the remnant of this narrow basin comparable to the Aqdarband ophiolitic remnant in north-east Iran. This metamorphic belt was intruded by 215 ± 15 Ma arc to collisional granites. In the passive margin of the Cimmerian block, on the Yazd region, the Silurian-Early Devonian syn-rift succession as well as the nearly continuous Upper Palaeozoic platform-type deposition was interrupted during the Middle to Late Triassic time, local erosion down to Devonian levels may be related to flexural bulge erosion. The collision event was not so strong to generate intensive deformation but was accompanied by some nappe thrusting onto the passive margin. It is finally unconformably covered by Liassic continental molassic deposits. Related to the onset of Neo-Tethyan back-arc opening in Early Jurassic to Mid-Eocene times, six periods of extensional-compressional events have differently influenced an elongated area, extending from the West Black Sea to Pamir. The Anarak-Jandaq terrane which was situated somewhere in this affected area, probably between the Kopeh Dagh and North Afghan platform, was completely detached from its source at the beginning of the Late Cretaceous

Focused time-lapse inversion of radio and audio magnetotelluric data

Geoelectrical techniques are widely used to monitor groundwater processes, while surprisingly few studies have considered audio (AMT) and radio (RMT) magnetotellurics for such purposes. In this numerical investigation, we analyze to what extent inversion results based on AMT and RMT monitoring data can be improved by (1) time-lapse difference inversion; (2) incorporation of statistical information about the expected model update (i.e., the model regularization is based on a geostatistical model); (3) using alternative model norms to quantify temporal changes (i.e., approximations of l(1) and Cauchy norms using iteratively reweighted least-squares), (4) constraining model updates to predefined ranges (i.e., using Lagrange Multipliers to only allow either increases or decreases of electrical resistivity with respect to background conditions). To do so, we consider a simple illustrative model and a more realistic test case related to seawater intrusion. The results are encouraging and show significant improvements when using time-lapse difference inversion with non l(2) model norms. Artifacts that may arise when imposing compactness of regions with temporal changes can be suppressed through inequality constraints to yield models without oscillations outside the true region of temporal changes. Based on these results, we recommend approximate l(1)-norm solutions as they can resolve both sharp and smooth interfaces within the same model. (C) 2012 Elsevier B.V. All rights reserved.

Pliensbachian magnetostratigraphy: new data from Paris Basin (France)

Sampling of an industrial drill string from the northeastern Paris Basin (Montcornet, France) provides early Jurassic magnetostratigraphic data coupled with biochronological control. About 375 paleomagnetic samples were obtained from a 145 m thick series of Pliensbachian rocks. A composite demagnetization thermal up to 300 C and an alternating field up to 80 mT were used to separate the magnetic components. A low unblocking temperature component (<250degreesC) with an inclination of about 64 is interpreted as a present-day field overprint. The characteristic remanent component with both normal and reversed antipodal directions was isolated between 5 and 50 mT. Twenty-nine polarity intervals were recognized. Correlation of these new results from the Paris Basin with data from the Breggia Gorge section (Ticino, southern Alps, Switzerland), which is generally considered as the reference section for Pliensbachian magnetostratigraphy, reveals almost identical patterns of magnetic polarity reversals. However, the correlation implies significant paleontological age discrepancies. Revised age assignments of biostratigraphic data of Breggia as well as an objective evaluation of the uncertainties on zonal boundaries in both Breggia and Moncornet resolve the initial discrepancies between magnetostratigraphic correlations and biostratigraphic ages. Hence, the sequence of magnetic reversals is significantly strengthened and the age calibration is notably improved for the Pliensbachian, a stage for which sections combining adequate magnetic signal and biostratigraphic constraints are still very few. (C) 2002 Elsevier Science B.V. All rights reserved.

GeoKernels: Modeling of Spatial Data on GeoManifolds

This paper presents a review of methodology for semi-supervised modeling with kernel methods, when the manifold assumption is guaranteed to be satisfied. It concerns environmental data modeling on natural manifolds, such as complex topographies of the mountainous regions, where environmental processes are highly influenced by the relief. These relations, possibly regionalized and nonlinear, can be modeled from data with machine learning using the digital elevation models in semi-supervised kernel methods. The range of the tools and methodological issues discussed in the study includes feature selection and semisupervised Support Vector algorithms. The real case study devoted to data-driven modeling of meteorological fields illustrates the discussed approach.

Structural characterization of rockfall sources in Yosemite Valley from remote sensing data and field surveys

Yosemite Valley poses significant rockfall hazard and related risk due to its glacially steepened walls and approximately 4 million visitors annually. To assess rockfall hazard, it is necessary to evaluate the geologic structure that contributes to the destabilization of rockfall sources and locate the most probable future source areas. Coupling new remote sensing techniques (Terrestrial Laser Scanning, Aerial Laser Scanning) and traditional field surveys, we investigated the regional geologic and structural setting, the orientation of the primary discontinuity sets for large areas of Yosemite Valley, and the specific discontinuity sets present at active rockfall sources. This information, combined with better understanding of the geologic processes that contribute to the progressive destabilization and triggering of granitic rock slabs, contributes to a more accurate rockfall susceptibility assessment for Yosemite Valley and elsewhere.

Intégration de méthodes de data mining dans le renseignement criminel : analyse par des structures issues de la théorie des graphes dans le profilage des stupéfiants

Data mining can be defined as the extraction of previously unknown and potentially useful information from large datasets. The main principle is to devise computer programs that run through databases and automatically seek deterministic patterns. It is applied in different fields of application, e.g., remote sensing, biometry, speech recognition, but has seldom been applied to forensic case data. The intrinsic difficulty related to the use of such data lies in its heterogeneity, which comes from the many different sources of information. The aim of this study is to highlight potential uses of pattern recognition that would provide relevant results from a criminal intelligence point of view. The role of data mining within a global crime analysis methodology is to detect all types of structures in a dataset. Once filtered and interpreted, those structures can point to previously unseen criminal activities. The interpretation of patterns for intelligence purposes is the final stage of the process. It allows the researcher to validate the whole methodology and to refine each step if necessary. An application to cutting agents found in illicit drug seizures was performed. A combinatorial approach was done, using the presence and the absence of products. Methods coming from the graph theory field were used to extract patterns in data constituted by links between products and place and date of seizure. A data mining process completed using graphing techniques is called ``graph mining''. Patterns were detected that had to be interpreted and compared with preliminary knowledge to establish their relevancy. The illicit drug profiling process is actually an intelligence process that uses preliminary illicit drug classes to classify new samples. Methods proposed in this study could be used \textit{a priori} to compare structures from preliminary and post-detection patterns. This new knowledge of a repeated structure may provide valuable complementary information to profiling and become a source of intelligence.