Exploration of alpine orographic precipitation patterns with radar image processing and clustering techniques

A methodology of exploratory data analysis investigating the phenomenon of orographic precipitation enhancement is proposed. The precipitation observations obtained from three Swiss Doppler weather radars are analysed for the major precipitation event of August 2005 in the Alps. Image processing techniques are used to detect significant precipitation cells/pixels from radar images while filtering out spurious effects due to ground clutter. The contribution of topography to precipitation patterns is described by an extensive set of topographical descriptors computed from the digital elevation model at multiple spatial scales. Additionally, the motion vector field is derived from subsequent radar images and integrated into a set of topographic features to highlight the slopes exposed to main flows. Following the exploratory data analysis with a recent algorithm of spectral clustering, it is shown that orographic precipitation cells are generated under specific flow and topographic conditions. Repeatability of precipitation patterns in particular spatial locations is found to be linked to specific local terrain shapes, e.g. at the top of hills and on the upwind side of the mountains. This methodology and our empirical findings for the Alpine region provide a basis for building computational data-driven models of orographic enhancement and triggering of precipitation. Copyright (C) 2011 Royal Meteorological Society .

Histopathological examination of nerve samples from pure neural leprosy patients: obtaining maximum information to improve diagnostic efficiency

Nerve biopsy examination is an important auxiliary procedure for diagnosing pure neural leprosy (PNL). When acid-fast bacilli (AFB) are not detected in the nerve sample, the value of other nonspecific histological alterations should be considered along with pertinent clinical, electroneuromyographical and laboratory data (the detection of Mycobacterium leprae DNA with polymerase chain reaction and the detection of serum anti-phenolic glycolipid 1 antibodies) to support a possible or probable PNL diagnosis. Three hundred forty nerve samples [144 from PNL patients and 196 from patients with non-leprosy peripheral neuropathies (NLN)] were examined. Both AFB-negative and AFB-positive PNL samples had more frequent histopathological alterations (epithelioid granulomas, mononuclear infiltrates, fibrosis, perineurial and subperineurial oedema and decreased numbers of myelinated fibres) than the NLN group. Multivariate analysis revealed that independently, mononuclear infiltrate and perineurial fibrosis were more common in the PNL group and were able to correctly classify AFB-negative PNL samples. These results indicate that even in the absence of AFB, these histopathological nerve alterations may justify a PNL diagnosis when observed in conjunction with pertinent clinical, epidemiological and laboratory data.

Seasonal variation and natural infection of Lutzomyia antunesi (Diptera: Psychodidae: Phlebotominae), an endemic species in the Orinoquia region of Colombia

Lutzomyia antunesi has been commonly reported in outbreaks of cutaneous leishmaniasis (CL) in the Orinoquia region of Colombia. The bionomics of this species were studied in the municipality of Villavicencio (Meta, Colombia). Sandflies were captured over the course of one week per month for one year in intradomiciliary, peridomiciliary and extradomiciliary housing areas. The captures were performed from 06:00 pm-06:00 am using CDC light traps and the females were processed for polymerase chain reaction (PCR) to detect Leishmania spp. A total of 22,097 specimens and 19 species were captured of which Lu. antunesi (89%) and Lutzomyia walkeri (5%) were the most abundant. Other species recognised as anthropophilic (Lutzomyia panamensis, Lutzomyia gomezi, Lutzomyia flaviscutellata and Lutzomyia fairtigi) were present in very low abundance (< 2%). Natural infection with Leishmania spp was detected using PCR in Lu. antunesi, Lu. panamensis and Lu. flavicutellata, showing infection rates of 1%, 4.8% and 7.5%, respectively. The present paper provides information on various ecological aspects of Lu. antunesi. An analysis of seasonality shows that this species increases in abundance in the hottest months (December, January and February), directly correlating with the maximum temperature and inversely correlating with precipitation. The natural infection rate is associated with the peaks of highest abundance.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.

Predicting the geographic distribution of Lutzomyia longipalpis (Diptera: Psychodidae) and visceral leishmaniasis in the state of Mato Grosso do Sul, Brazil

To understand the geographic distribution of visceral leishmaniasis (VL) in the state of Mato Grosso do Sul (MS), Brazil, both the climatic niches of Lutzomyia longipalpis and VL cases were analysed. Distributional data were obtained from 55 of the 79 counties of MS between 2003-2012. Ecological niche models (ENM) of Lu. longipalpis and VL cases were produced using the maximum entropy algorithm based on eight climatic variables. Lu. longipalpis showed a wide distribution in MS. The highest climatic suitability for Lu. longipalpis was observed in southern MS. Temperature seasonality and annual mean precipitation were the variables that most influenced these models. Two areas of high climatic suitability for the occurrence of VL cases were predicted: one near Aquidauana and another encompassing several municipalities in the southeast region of MS. As expected, a large overlap between the models for Lu. longipalpis and VL cases was detected. Northern and northwestern areas of MS were suitable for the occurrence of cases, but did not show high climatic suitability for Lu. longipalpis . ENM of vectors and human cases provided a greater understanding of the geographic distribution of VL in MS, which can be applied to the development of future surveillance strategies.

Influència del clima en la fenologia de la floració de Globularia alypum: anàlisi a partir de mostres d'herbari

Globularia alypum és un arbust termòfil típic de la regió oest del Mediterrani. Aquest nanofaneròfit té un patró de floració molt ampli. Les poblacions primerenques floreixen des del Juliol fins al Desembre i les poblacions tardanes des del Gener fins a l’Abril. L’objectiu del projecte és esbrinar quins factors climàtics determinen la fenologia de la floració de l’espècie a partir de mostres d’herbari i conèixer l’abast geogràfic d’ambdós patrons de floració a la península Ibèrica i França. També l’estudi d’altres característiques fenològiques en relació amb el clima: els braquiblasts i les espigues de capítols. Els plecs d’herbari han estat georeferenciats per ser utilitzats en un entorn SIG junt amb les dades del clima i s’ha realitzat una anàlisi estadística. El resultats mostren que les poblacions primerenques creixen majoritàriament en zones costaneres mentre que les tardanes ho fan en zones continentals. Les poblacions primerenques creixen en àrees on les temperatures (mitjana, mínima, màxima i mínima del mes més fred) són més altes que les tardanes per l’efecte temperador del mar. Al analitzar totes les poblacions, el nombre de mesos des de la transició floral es correlaciona negativament amb la temperatura mínima mensual i la temperatura mitjana mensual en tots els mesos, mentre que amb la temperatura màxima mensual entre Octubre i Abril. Al estudiar sols les poblacions tardanes les correlacions són similars, però no en les primerenques on no se’n troba cap. Les correlacions entre el nombre de mesos des de la transició floral i les variables climàtiques anuals són iguals, a més de trobar-se una correlació negativa amb la temperatura mínima del mes més fred i l’índex tèrmic, i positiva amb l’índex de continentalitat. No s’ha trobat cap diferència en presència de braquiblasts segons la població, però sí en les espigues de capítols, essent més freqüents en les poblacions primerenques. Es suggereix que la temperatura és un factor més determinant per a les poblacions tardanes, mentre que es desconeix si és la precipitació en les primerenques. Tampoc es té coneixement sobre si les diferències entre els dos patrons són genètiques o fenotípiques i quin mecanisme possibilita el desenvolupament de braquiblasts i espigues de capítols.

Efficiency of Gumbel analyses for determining extreme daily precipitation in Switzerland

Gumbel analyses were carried out on rainfall time-series at 151 locations in Switzerland for 4 different periods of 30 years in order to estimate daily extreme precipitation for a return period of 100 years. Those estimations were compared with maximal daily values measured during the last 100 years (1911-2010) to test the efficiency of these analyses. This comparison shows that these analyses provide good results for 50 to 60% locations in this country from rainfall time-series 1961-1990 and 1980-2010. On the other hand, daily precipitation with a return period of 100 years is underestimated at most locations from time-series 1931-1960 and especially 1911-1940. Such underestimation results from the increase of maximal daily precipitation recorded from 1911 to 2010 at 90% locations in Switzerland.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

Reconstructing paleoelevation in eroded orogens

Hydrogen isotope and Ar-40/Ar-39 geochronological data are presented from muscovite within a crustal-scale extensional detachment of the Shuswap Metamorphic Complex, North American Cordillera. The hydrogen isotope compositions (deltaD(ms)) of precisely dated muscovite attain values as low as -156parts per thousand in the detachment mylonite, whereas footwall quartzite has a deltaD(ms) value of -81parts per thousand. The very low deltaD(ms) values in the detachment are best explained by infiltration of meteoric water, with maximum deltaD values of -135parts per thousand +/- 3parts per thousand, during extensional unroofing of the orogen at 49.0-47.9 Ma. On the basis of the empirically determined relationship between elevation and isotopic composition of precipitation, the reconstructed early Eocene paleoelevations of the orogen are 4060 +/- 250 m to 4320 +/- 250 m, at least 1000 m higher than the highest present-day peaks. We propose that the isotopic composition of surface-derived waters in extensional detachments represents a newly recognized method to estimate maximum paleoelevations attained immediately preceding extensional orogenic collapse.

Spatial relationship between the atmospheric circulation and the precipitation measured in the western Swiss Alps by means of the analogue method

An adaptation technique based on the synoptic atmospheric circulation to forecast local precipitation, namely the analogue method, has been implemented for the western Swiss Alps. During the calibration procedure, relevance maps were established for the geopotential height data. These maps highlight the locations were the synoptic circulation was found of interest for the precipitation forecasting at two rain gauge stations (Binn and Les Marécottes) that are located both in the alpine Rhône catchment, at a distance of about 100 km from each other. These two stations are sensitive to different atmospheric circulations. We have observed that the most relevant data for the analogue method can be found where specific atmospheric circulation patterns appear concomitantly with heavy precipitation events. Those skilled regions are coherent with the atmospheric flows illustrated, for example, by means of the back trajectories of air masses. Indeed, the circulation recurrently diverges from the climatology during days with strong precipitation on the southern part of the alpine Rhône catchment. We have found that for over 152 days with precipitation amount above 50 mm at the Binn station, only 3 did not show a trajectory of a southerly flow, meaning that such a circulation was present for 98% of the events. Time evolution of the relevance maps confirms that the atmospheric circulation variables have significantly better forecasting skills close to the precipitation period, and that it seems pointless for the analogue method to consider circulation information days before a precipitation event as a primary predictor. Even though the occurrence of some critical circulation patterns leading to heavy precipitation events can be detected by precursors at remote locations and 1 week ahead (Grazzini, 2007; Martius et al., 2008), time extrapolation by the analogue method seems to be rather poor. This would suggest, in accordance with previous studies (Obled et al., 2002; Bontron and Obled, 2005), that time extrapolation should be done by the Global Circulation Model, which can process atmospheric variables that can be used by the adaptation method.

Cell disposition of raltegravir and newer antiretrovirals in HIV-infected patients: high inter-individual variability in raltegravir cellular penetration.

Objectives The site of pharmacological activity of raltegravir is intracellular. Our aim was to determine the extent of raltegravir cellular penetration and whether raltegravir total plasma concentration (C(tot)) predicts cellular concentration (C(cell)). Methods Open-label, prospective, pharmacokinetic study on HIV-infected patients on a stable raltegravir-containing regimen. Plasma and peripheral blood mononuclear cells were simultaneously collected during a 12 h dosing interval after drug intake. C(tot) and C(cell) of raltegravir, darunavir, etravirine, maraviroc and ritonavir were measured by liquid chromatography coupled to tandem mass spectrometry after protein precipitation. Longitudinal mixed effects analysis was applied to the C(cell)/C(tot) ratio. Results Ten HIV-infected patients were included. The geometric mean (GM) raltegravir total plasma maximum concentration (C(max)), minimum concentration (C(min)) and area under the time-concentration curve from 0-12 h (AUC(0-12)) were 1068 ng/mL, 51.1 ng/mL and 4171 ng·h/mL, respectively. GM raltegravir cellular C(max), C(min) and AUC(0-12) were 27.5 ng/mL, 2.9 ng/mL and 165 ng·h/mL, respectively. Raltegravir C(cell) corresponded to 5.3% of C(tot) measured simultaneously. Both concentrations fluctuate in parallel, with C(cell)/C(tot) ratios remaining fairly constant for each patient without a significant time-related trend over the dosing interval. The AUC(cell)/AUC(tot) GM ratios for raltegravir, darunavir and etravirine were 0.039, 0.14 and 1.55, respectively. Conclusions Raltegravir C(cell) correlated with C(tot) (r = 0.86). Raltegravir penetration into cells is low overall (∼5% of plasma levels), with distinct raltegravir cellular penetration varying by as much as 15-fold between patients. The importance of this finding in the context of development of resistance to integrase inhibitors needs to be further investigated.

Are the maximum hardness and minimum polarizability principles always obeyed in nontotally symmetric vibrations?

An overview is given on a study which showed that not only in chemical reactions but also in the favorable case of nontotally symmetric vibrations where the chemical and external potentials keep approximately constant, the generalized maximum hardness principle (GMHP) and generalized minimum polarizability principle (GMPP) may not be obeyed. A method that allows an accurate determination of the nontotally symmetric molecular distortions with more marked GMPP or anti-GMPP character through diagonalization of the polarizability Hessian matrix is introduced