A limited area model intercomparison on the "Montserrat-2000" flash-flood event using statistical and deterministic methods

In the scope of the European project Hydroptimet, INTERREG IIIB-MEDOCC programme, limited area model (LAM) intercomparison of intense events that produced many damages to people and territory is performed. As the comparison is limited to single case studies, the work is not meant to provide a measure of the different models' skill, but to identify the key model factors useful to give a good forecast on such a kind of meteorological phenomena. This work focuses on the Spanish flash-flood event, also known as "Montserrat-2000" event. The study is performed using forecast data from seven operational LAMs, placed at partners' disposal via the Hydroptimet ftp site, and observed data from Catalonia rain gauge network. To improve the event analysis, satellite rainfall estimates have been also considered. For statistical evaluation of quantitative precipitation forecasts (QPFs), several non-parametric skill scores based on contingency tables have been used. Furthermore, for each model run it has been possible to identify Catalonia regions affected by misses and false alarms using contingency table elements. Moreover, the standard "eyeball" analysis of forecast and observed precipitation fields has been supported by the use of a state-of-the-art diagnostic method, the contiguous rain area (CRA) analysis. This method allows to quantify the spatial shift forecast error and to identify the error sources that affected each model forecasts. High-resolution modelling and domain size seem to have a key role for providing a skillful forecast. Further work is needed to support this statement, including verification using a wider observational data set.

Introduction of sensor spectral response into image fusion methods. Application to wavelet-based methods

Usual image fusion methods inject features from a high spatial resolution panchromatic sensor into every low spatial resolution multispectral band trying to preserve spectral signatures and improve spatial resolution to that of the panchromatic sensor. The objective is to obtain the image that would be observed by a sensor with the same spectral response (i.e., spectral sensitivity and quantum efficiency) as the multispectral sensors and the spatial resolution of the panchromatic sensor. But in these methods, features from electromagnetic spectrum regions not covered by multispectral sensors are injected into them, and physical spectral responses of the sensors are not considered during this process. This produces some undesirable effects, such as resolution overinjection images and slightly modified spectral signatures in some features. The authors present a technique which takes into account the physical electromagnetic spectrum responses of sensors during the fusion process, which produces images closer to the image obtained by the ideal sensor than those obtained by usual wavelet-based image fusion methods. This technique is used to define a new wavelet-based image fusion method.

Measures of health and disease in Africa: are current methods giving us useful information about trends in cardiovascular diseases?

An enormous burst of interest in the public health burden from chronic disease in Africa has emerged as a consequence of efforts to estimate global population health. Detailed estimates are now published for Africa as a whole and each country on the continent. These data have formed the basis for warnings about sharp increases in cardiovascular disease (CVD) in the coming decades. In this essay we briefly examine the trajectory of social development on the continent and its consequences for the epidemiology of CVD and potential control strategies. Since full vital registration has only been implemented in segments of South Africa and the island nations of Seychelles and Mauritius - formally part of WHO-AFRO - mortality data are extremely limited. Numerous sample surveys have been conducted but they often lack standardization or objective measures of health status. Trend data are even less informative. However, using the best quality data available, age-standardized trends in CVD are downward, and in the case of stroke, sharply so. While acknowledging that the extremely limited available data cannot be used as the basis for inference to the continent, we raise the concern that general estimates based on imputation to fill in the missing mortality tables may be even more misleading. No immediate remedies to this problem can be identified, however bilateral collaborative efforts to strength local educational institutions and governmental agencies rank as the highest priority for near term development.

Novel methods improve prediction of species' distributions from occurrence data

Prediction of species' distributions is central to diverse applications in ecology, evolution and conservation science. There is increasing electronic access to vast sets of occurrence records in museums and herbaria, yet little effective guidance on how best to use this information in the context of numerous approaches for modelling distributions. To meet this need, we compared 16 modelling methods over 226 species from 6 regions of the world, creating the most comprehensive set of model comparisons to date. We used presence-only data to fit models, and independent presence-absence data to evaluate the predictions. Along with well-established modelling methods such as generalised additive models and GARP and BIOCLIM, we explored methods that either have been developed recently or have rarely been applied to modelling species' distributions. These include machine-learning methods and community models, both of which have features that may make them particularly well suited to noisy or sparse information, as is typical of species' occurrence data. Presence-only data were effective for modelling species' distributions for many species and regions. The novel methods consistently outperformed more established methods. The results of our analysis are promising for the use of data from museums and herbaria, especially as methods suited to the noise inherent in such data improve.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Crescent and star shapes of members of the Chlamydiales order: impact of fixative methods.

Members of the Chlamydiales order all share a biphasic lifecycle alternating between small infectious particles, the elementary bodies (EBs) and larger intracellular forms able to replicate, the reticulate bodies. Whereas the classical Chlamydia usually harbours round-shaped EBs, some members of the Chlamydia-related families display crescent and star-shaped morphologies by electron microscopy. To determine the impact of fixative methods on the shape of the bacterial cells, different buffer and fixative combinations were tested on purified EBs of Criblamydia sequanensis, Estrella lausannensis, Parachlamydia acanthamoebae, and Waddlia chondrophila. A linear discriminant analysis was performed on particle metrics extracted from electron microscopy images to recognize crescent, round, star and intermediary forms. Depending on the buffer and fixatives used, a mixture of alternative shapes were observed in varying proportions with stars and crescents being more frequent in C. sequanensis and P. acanthamoebae, respectively. No tested buffer and chemical fixative preserved ideally the round shape of a majority of bacteria and other methods such as deep-freezing and cryofixation should be applied. Although crescent and star shapes could represent a fixation artifact, they certainly point towards a diverse composition and organization of membrane proteins or intracellular structures rather than being a distinct developmental stage.

Bacterial contamination of platelet concentrates: pathogen detection and inactivation methods

Whereas the reduction of transfusion related viral transmission has been a priority during the last decade, bacterial infection transmitted by transfusion still remains associated to a high morbidity and mortality, and constitutes the most frequent infectious risk of transfusion. This problem especially concerns platelet concentrates because of their favorable bacterial growth conditions. This review gives an overview of platelet transfusion-related bacterial contamination as well as on the different strategies to reduce this problem by using either bacterial detection or inactivation methods.

Effect of collection and preprocessing methods on neutrophil elastase plasma concentrations.

OBJECTIVES: Elevated plasma levels of the elastase alpha 1-proteinase inhibitor complex (E-alpha 1 PI) have been proposed as a marker of bacterial infection and neutrophil activation. Liberation of elastase from neutrophils after collection of blood may cause falsely elevated results. Collection methods have not been validated for critically ill neonates and children. We evaluated the influence of preanalytical methods on E-alpha 1 PI results including the recommended collection into EDTA tubes. DESIGN AND METHODS: First, we compared varying acceleration speeds and centrifugation times. Centrifugation at 1550 g for 3 min resulted in reliable preparation of leukocyte free plasma. Second, we evaluated all collection tubes under consideration for absorption of E-alpha 1 PI. Finally, 12 sets of samples from healthy adults and 42 sets obtained from critically ill neonates and children were distributed into the various sampling tubes. Samples were centrifuged within 15 min of collection and analyzed with a new turbidimetric assay adapted to routine laboratory analyzers. RESULTS: One of the two tubes containing a plasma-cell separation gel absorbed 22.1% of the E-alpha 1 PI content. In the remaining tubes without absorption of E-alpha 1 PI no differences were observed for samples from healthy adult patients. However, in samples from critically ill neonates or children, significantly higher results were obtained for plain Li-heparin tubes (mean = 183 micrograms/L), EDTA tubes (mean = 93 micrograms/L), and citrate tubes (mean = 88.5 micrograms/L) than for the Li-hep tube with cell-plasma separation gel and no absorption of E-alpha 1 PI (mean = 62.4 micrograms/L, p < 0.01). CONCLUSION: Contrary to healthy adults, E-alpha 1 PI results in plasma samples from critically ill neonates and children depend on the type of collection tube.

Organic material decomposition and nutrient dynamics in a mulch system enriched with leguminous trees in the Amazon

The new techniques proposed for agriculture in the Amazon region include rotational fallow systems enriched with leguminous trees and the replacement of biomass burning by mulching. Decomposition and nutrient release from mulch were studied using fine-mesh litterbags with five different leguminous species and the natural fallow vegetation as control. Samples from each treatment were analyzed for total C, N, P, K, Ca, Mg, lignin, cellulose content and soluble polyphenol at different sampling times over the course of one year. The decomposition rate constant varied with species and time. Weight loss from the decomposed litter bag material after 96 days was 30.1 % for Acacia angustissima, 32.7 % for Sclerolobium paniculatum, 33.9 % for Iinga edulis and the Fallow vegetation, 45.2 % for Acacia mangium and 63.6 % for Clitoria racemosa. Immobilization of N and P was observed in all studied treatments. Nitrogen mineralization was negatively correlated with phenol, C-to-N ratio, lignin + phenol/N ratio, and phenol/phosphorus ratios and with N content in the litterbag material. After 362 days of field incubation, an average (of all treatments), 3.3 % K, 32.2 % Ca and 22.4 % Mg remained in the mulch. Results confirm that low quality and high amount of organic C as mulch application are limiting for the quantity of energy available for microorganisms and increase the nutrient immobilization for biomass decomposition, which results in competition for nutrients with the crop plants.

Comparison of estimation methods of soil strength in five soils

In agriculture, the soil strength is used to describe the susceptibility to deformation by pressure caused by agricultural machine. The purpose of this study was to compare different methods for estimating the inherent soil strength and to identify their suitability for the evaluation of load support capacity, compaction susceptibility and root growth. The physical, chemical, mineralogical and intrinsic strength properties of seven soil samples, collected from five sampling pits at different locations in Brazil, were measured. Four clay (CS) and three sandy clay loam (SCL) soils were used. The clay soils were collected on a farm in Santo Ângelo, RS (28 º 16 ' 16 '' S; 54 º 13 ' 11 '' W 290 m); A and B horizons at the Universidade Federal de Lavras, Lavras, MG (21 º 13 ' 47 '' S; 44 º 58 ' 6'' W; 918 m) and on the farm Sygenta, in Uberlandia, MG (18 º 58 ' 37 '' S; 48 º 12 ' 05 '' W 866 m). The sandy clay loam soils were collected in Aracruz, ES (19 º 47 ' 10 '' S; 40 º 16 ' 29 '' W 81 m), and on the farm Xavier, Lavras, MG (21 º 13 ' 24 '' S; 45 º 05 ' 00 '' W; 844 m). Soil strength was estimated based on measurements of: (a) a pneumatic consolidometer, (b) manual pocket (non-rotating) penetrometer; and (c) automatic (rotating) penetrometer. The results of soil strength properties were similar by the three methods. The soil structure had a significant influence on soil strength. Results of measurements with both the manual pocket and the electric penetrometer were similar, emphasizing the influence of soil texture. The data showed that, to enhance the reliability of predictions of preconsolidation pressure by penetrometers, it is better to separate the soils into the different classes, rather than analyze them jointly. It can be concluded that the consolidometer method, although expensive, is the best when evaluations of load support capacity and compaction susceptibility of soil samples are desired.

Root system distribution of sugar cane as related to nitrogen fertilization, evaluated by two methods: monolith and probes

Few studies on sugar cane have evaluated the root system of the crop, in spite of its importance. This is mainly due to the difficulty of evaluation and high variability of results. The objective of this study was to develop an evaluation method of the cane root system by means of probes so as to evaluate the mass, distribution and metabolically active roots related to N fertilization at planting. For this purpose, an experiment was conducted in an Arenic Kandiustults with medium texture in Jaboticabal/SP, in a randomized block design with four replications and four treatments: control (without N) and 40, 80 and 120 kg ha-1 of N applied in the form of urea in the planting furrow of the cane variety SP81 3250. One week before harvest, a urea-15N solution was applied at the cane stalk base to detect active metabolism in the root system. Trenches of 1.5 m length and 0.6 m depth were opened between two sugar cane rows for root sampling by two methods: monoliths (0.3, 0.2 and 0.15 m wide, deep and long respectively) taken from the trench wall and by probe (internal diameter 0.055 m). For each method, 15 samples per plot were collected. The roots were separated from the soil in a sieve (2 mm mesh), oven-dried (at 65 ºC) and the dry matter was measured. Root sampling by probes resulted in root mass that did not differ from the evaluation in monoliths, indicating that this evaluation method may be used for sugar cane root mass, although neither the root distribution in the soil profile nor the rhizome mass were efficiently evaluated, due to the small sample volume. Nitrogen fertilization at planting did not result in a greater root accumulation in the sugar cane plant, but caused changes in the distribution of the root system in the soil. The absence of N fertilization led to a better root distribution in the soil profile, with 50, 34 and 16 % in the 0-0.2, 0.2-0.4 and 0.4-0.6 m layers, respectively; in the fertilized treatments the roots were concentrated in the surface layer, with on average 70, 17 and 13 % for the same layers. The metabolically active roots were concentrated in the center of the cane stool, amounting to 40 % of the total root mass, regardless of N fertilization (application of 120 kg ha-1 N or without N).

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved. The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10 X 10 mt in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.

Socio-economic inequalities in reported depression in Spain : a decomposition approach

Recent evidence questions some conventional view on the existence of income-related inequalities in depression suggesting in turn that other determinants might be in place, such as activity status and educational attainment. Evidence of socio-economic inequalities is especially relevant in countries such as Spain that have a limited coverage of mental health care and are regionally heterogeneous. This paper aims at measuring and explaining the degree of socio-economic inequality in reported depression in Spain. We employ linear probability models to estimate the concentration index and its decomposition drawing from 2003 edition of the Spanish National Health Survey, the most recent representative health survey in Spain. Our findings point towards the existence of avoidable inequalities in the prevalence of reported depression. However, besides ¿pure income effects¿ explaining 37% of inequality, economic activity status (28%), education (15%) and demographics (15%) play also a key encompassing role. Although high income implies higher resources to invest and cure (mental) illness, environmental factors influencing in peoples perceived social status act as indirect path as explaining the prevalence of depression. Finally, we find evidence of a gender effect, gender social-economic inequality in income is mainly avoidable.

Selection methods regulate evolution of cooperation in digital evolution.

A key, yet often neglected, component of digital evolution and evolutionary models is the 'selection method' which assigns fitness (number of offspring) to individuals based on their performance scores (efficiency in performing tasks). Here, we study with formal analysis and numerical experiments the evolution of cooperation under the five most common selection methods (proportionate, rank, truncation-proportionate, truncation-uniform and tournament). We consider related individuals engaging in a Prisoner's Dilemma game where individuals can either cooperate or defect. A cooperator pays a cost, whereas its partner receives a benefit, which affect their performance scores. These performance scores are translated into fitness by one of the five selection methods. We show that cooperation is positively associated with the relatedness between individuals under all selection methods. By contrast, the change in the performance benefit of cooperation affects the populations' average level of cooperation only under the proportionate methods. We also demonstrate that the truncation and tournament methods may introduce negative frequency-dependence and lead to the evolution of polymorphic populations. Using the example of the evolution of cooperation, we show that the choice of selection method, though it is often marginalized, can considerably affect the evolutionary dynamics.