Machine learning algorithms for spatial data. Case studies: Environmental pollution, natural hazards, renewable resources

This paper presents general problems and approaches for the spatial data analysis using machine learning algorithms. Machine learning is a very powerful approach to adaptive data analysis, modelling and visualisation. The key feature of the machine learning algorithms is that they learn from empirical data and can be used in cases when the modelled environmental phenomena are hidden, nonlinear, noisy and highly variable in space and in time. Most of the machines learning algorithms are universal and adaptive modelling tools developed to solve basic problems of learning from data: classification/pattern recognition, regression/mapping and probability density modelling. In the present report some of the widely used machine learning algorithms, namely artificial neural networks (ANN) of different architectures and Support Vector Machines (SVM), are adapted to the problems of the analysis and modelling of geo-spatial data. Machine learning algorithms have an important advantage over traditional models of spatial statistics when problems are considered in a high dimensional geo-feature spaces, when the dimension of space exceeds 5. Such features are usually generated, for example, from digital elevation models, remote sensing images, etc. An important extension of models concerns considering of real space constrains like geomorphology, networks, and other natural structures. Recent developments in semi-supervised learning can improve modelling of environmental phenomena taking into account on geo-manifolds. An important part of the study deals with the analysis of relevant variables and models' inputs. This problem is approached by using different feature selection/feature extraction nonlinear tools. To demonstrate the application of machine learning algorithms several interesting case studies are considered: digital soil mapping using SVM, automatic mapping of soil and water system pollution using ANN; natural hazards risk analysis (avalanches, landslides), assessments of renewable resources (wind fields) with SVM and ANN models, etc. The dimensionality of spaces considered varies from 2 to more than 30. Figures 1, 2, 3 demonstrate some results of the studies and their outputs. Finally, the results of environmental mapping are discussed and compared with traditional models of geostatistics.

Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.

Debris flow susceptibility mapping at a regional scale along the National Road N7, Argentina

A first assessment of debris flow susceptibility at a large scale was performed along the National Road N7, Argentina. Numerous catchments are prone to debris flows and likely to endanger the road-users. A 1:50,000 susceptibility map was created. The use of a DEM (grid 30 m) associated to three complementary criteria (slope, contributing area, curvature) allowed the identification of potential source areas. The debris flow spreading was estimated using a process- and GISbased model (Flow-R) based on basic probabilistic and energy calculations. The best-fit values for the coefficient of friction and the mass-to-drag ratio of the PCM model were found to be ? = 0.02 and M/D = 180 and the resulting propagation on one of the calibration site was validated using the Coulomb friction model. The results are realistic and will be useful to determine which areas need to be prioritized for detailed studies.

Mapping of Environmental Data Using Kernel-Based Methods

Recently, kernel-based Machine Learning methods have gained great popularity in many data analysis and data mining fields: pattern recognition, biocomputing, speech and vision, engineering, remote sensing etc. The paper describes the use of kernel methods to approach the processing of large datasets from environmental monitoring networks. Several typical problems of the environmental sciences and their solutions provided by kernel-based methods are considered: classification of categorical data (soil type classification), mapping of environmental and pollution continuous information (pollution of soil by radionuclides), mapping with auxiliary information (climatic data from Aral Sea region). The promising developments, such as automatic emergency hot spot detection and monitoring network optimization are discussed as well.

Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis.

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

Predictive analysis and mapping of indoor radon concentrations in Switzerland

It is estimated that around 230 people die each year due to radon (222Rn) exposure in Switzerland. 222Rn occurs mainly in closed environments like buildings and originates primarily from the subjacent ground. Therefore it depends strongly on geology and shows substantial regional variations. Correct identification of these regional variations would lead to substantial reduction of 222Rn exposure of the population based on appropriate construction of new and mitigation of already existing buildings. Prediction of indoor 222Rn concentrations (IRC) and identification of 222Rn prone areas is however difficult since IRC depend on a variety of different variables like building characteristics, meteorology, geology and anthropogenic factors. The present work aims at the development of predictive models and the understanding of IRC in Switzerland, taking into account a maximum of information in order to minimize the prediction uncertainty. The predictive maps will be used as a decision-support tool for 222Rn risk management. The construction of these models is based on different data-driven statistical methods, in combination with geographical information systems (GIS). In a first phase we performed univariate analysis of IRC for different variables, namely the detector type, building category, foundation, year of construction, the average outdoor temperature during measurement, altitude and lithology. All variables showed significant associations to IRC. Buildings constructed after 1900 showed significantly lower IRC compared to earlier constructions. We observed a further drop of IRC after 1970. In addition to that, we found an association of IRC with altitude. With regard to lithology, we observed the lowest IRC in sedimentary rocks (excluding carbonates) and sediments and the highest IRC in the Jura carbonates and igneous rock. The IRC data was systematically analyzed for potential bias due to spatially unbalanced sampling of measurements. In order to facilitate the modeling and the interpretation of the influence of geology on IRC, we developed an algorithm based on k-medoids clustering which permits to define coherent geological classes in terms of IRC. We performed a soil gas 222Rn concentration (SRC) measurement campaign in order to determine the predictive power of SRC with respect to IRC. We found that the use of SRC is limited for IRC prediction. The second part of the project was dedicated to predictive mapping of IRC using models which take into account the multidimensionality of the process of 222Rn entry into buildings. We used kernel regression and ensemble regression tree for this purpose. We could explain up to 33% of the variance of the log transformed IRC all over Switzerland. This is a good performance compared to former attempts of IRC modeling in Switzerland. As predictor variables we considered geographical coordinates, altitude, outdoor temperature, building type, foundation, year of construction and detector type. Ensemble regression trees like random forests allow to determine the role of each IRC predictor in a multidimensional setting. We found spatial information like geology, altitude and coordinates to have stronger influences on IRC than building related variables like foundation type, building type and year of construction. Based on kernel estimation we developed an approach to determine the local probability of IRC to exceed 300 Bq/m3. In addition to that we developed a confidence index in order to provide an estimate of uncertainty of the map. All methods allow an easy creation of tailor-made maps for different building characteristics. Our work is an essential step towards a 222Rn risk assessment which accounts at the same time for different architectural situations as well as geological and geographical conditions. For the communication of 222Rn hazard to the population we recommend to make use of the probability map based on kernel estimation. The communication of 222Rn hazard could for example be implemented via a web interface where the users specify the characteristics and coordinates of their home in order to obtain the probability to be above a given IRC with a corresponding index of confidence. Taking into account the health effects of 222Rn, our results have the potential to substantially improve the estimation of the effective dose from 222Rn delivered to the Swiss population.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

Getting the most sulfate from soil: Regulation of sulfate uptake transporters in Arabidopsis.

Sulfur (S) is an essential macronutrient for all living organisms. Plants require large amounts of sulfate for growth and development, and this serves as a major entry point of sulfate into the food web. Plants acquire S in its ionic form from the soil; they have evolved tightly controlled mechanisms for the regulation of sulfate uptake in response to its external and internal availability. In the model plant Arabidopsis thaliana, the first key step in sulfate uptake is presumed to be carried out exclusively by only two high-affinity sulfate transporters: SULTR1;1 and SULTR1;2. A better understanding of the mode of regulation for these two transporters is crucial because they constitute the first determinative step in balancing sulfate in respect to its supply and demand. Here, we review the recent progress achieved in our comprehension of (i) mechanisms that regulate these two high-affinity sulfate transporters at the transcriptional and post-transcriptional levels, and (ii) their structure-function relationship. Such progress is important to enable biotechnological and agronomic strategies aimed at enhancing sulfate uptake and improving crop yield in S-deficient soils.

High-resolution spatial mapping of changes in the neurochemical profile after focal ischemia in mice.

After ischemic stroke, the ischemic damage to brain tissue evolves over time and with an uneven spatial distribution. Early irreversible changes occur in the ischemic core, whereas, in the penumbra, which receives more collateral blood flow, the damage is more mild and delayed. A better characterization of the penumbra, irreversibly damaged and healthy tissues is needed to understand the mechanisms involved in tissue death. MRSI is a powerful tool for this task if the scan time can be decreased whilst maintaining high sensitivity. Therefore, we made improvements to a (1) H MRSI protocol to study middle cerebral artery occlusion in mice. The spatial distribution of changes in the neurochemical profile was investigated, with an effective spatial resolution of 1.4 μL, applying the protocol on a 14.1-T magnet. The acquired maps included the difficult-to-separate glutamate and glutamine resonances and, to our knowledge, the first mapping of metabolites γ-aminobutyric acid and glutathione in vivo, within a metabolite measurement time of 45 min. The maps were in excellent agreement with findings from single-voxel spectroscopy and offer spatial information at a scan time acceptable for most animal models. The metabolites measured differed with respect to the temporal evolution of their concentrations and the localization of these changes. Specifically, lactate and N-acetylaspartate concentration changes largely overlapped with the T(2) -hyperintense region visualized with MRI, whereas changes in cholines and glutathione affected the entire middle cerebral artery territory. Glutamine maps showed elevated levels in the ischemic striatum until 8 h after reperfusion, and until 24 h in cortical tissue, indicating differences in excitotoxic effects and secondary energy failure in these tissue types. Copyright © 2011 John Wiley & Sons, Ltd.

Fungi, bacteria and soil pH: the oxalate-carbonate pathway as a model for metabolic interaction

The oxalatecarbonate pathway involves the oxidation of calcium oxalate to low-magnesium calcite and represents a potential long-term terrestrial sink for atmospheric CO2. In this pathway, bacterial oxalate degradation is associated with a strong local alkalinization and subsequent carbonate precipitation. In order to test whether this process occurs in soil, the role of bacteria, fungi and calcium oxalate amendments was studied using microcosms. In a model system with sterile soil amended with laboratory cultures of oxalotrophic bacteria and fungi, the addition of calcium oxalate induced a distinct pH shift and led to the final precipitation of calcite. However, the simultaneous presence of bacteria and fungi was essential to drive this pH shift. Growth of both oxalotrophic bacteria and fungi was confirmed by qPCR on the frc (oxalotrophic bacteria) and 16S rRNA genes, and the quantification of ergosterol (active fungal biomass) respectively. The experiment was replicated in microcosms with non-sterilized soil. In this case, the bacterial and fungal contribution to oxalate degradation was evaluated by treatments with specific biocides (cycloheximide and bronopol). Results showed that the autochthonous microflora oxidized calcium oxalate and induced a significant soil alkalinization. Moreover, data confirmed the results from the model soil showing that bacteria are essentially responsible for the pH shift, but require the presence of fungi for their oxalotrophic activity. The combined results highlight that the interaction between bacteria and fungi is essential to drive metabolic processes in complex environments such as soil.

Plasmodium falciparum merozoite surface protein 2: epitope mapping and biological activity analysis of specific antibodies

Background: Plasmodium falciparum(P. falciparum) merozoite surfaceprotein 2 (MSP-2) is one of bloodstage proteins that are associated withprotection from malaria. MSP-2 consistsof a highly polymorphic centralrepeat region flanked by a dimorphicregion that defines the two allelicfamilies, 3D7 and FC27; N- and Cterminalregions are conserved domains.Long synthetic peptides (LSP)representing the two allelic familiesof MSP-2 and constant regions arerecognized by sera from donors livingin endemic areas; and specific antibodies(Abs) are associated with protectionand active in antibody dependentcellular inhibition (ADCI) in vitro.However, the fine specificity ofAb response to the two allelic familiesof MSP-2 is unknown. Methods: Peptidesrepresenting dimorphic regionof 3D7 and FC27 families and theirC-terminal (common fragment to thetwo families) termed 3D7-D (88 aa),FC27-D (48 aa) and C (40 aa) respectivelywere synthesized. Overlapping20 mer peptides covering dimorphicand constant regions of two familieswere also synthesized for epitopemapping. Human sera were obtainedfrom donors living in malaria endemicareas. SpecificDand CregionsAbs were purified from single or poolhuman sera. Sera from mice were obtainedafter immunization with thetwo families LSP mixture in three differentadjuvants: alhydrogel (Alum),Glucopyranosyl Lipid Adjuvant-Stableoil-in-water Emulsion (GLA-SE)and Virosome. For ADCI, P. falciparum(strain 3D7) parasite wasmaintained in culture at 0.5% parasitemiaand 4% hematocrit in air tightbox at love oxygen (2%) and 37 ºC.Results: We identified several epitopesfrom the dimorphic and constantregions of both families of MSP-2, inmice and humans (adults and children).In human, most recognizedepitopes were the same in differentendemic regions for each domain ofthe two families of MSP-2. In mice,the differential recognition of epitopewas depending on the strain of mouseand interestingly on the adjuvantused. GLA-SE and alum as adjuvantswere more often associated with therecognition of multiple epitopes thanvirosomes. Epitope-specific Abs recognizednative merozoites of P.falciparum and were active in ADCIto block development of parasite.Conclusion: The delineation of a limitednumber of epitopes could be exploitedto develop MSP-2 vaccinesactive on both allelic families ofMSP-2.

Accurate performance of a rat model of schizophrenia in the water maze depends on visual cue availability and stability: a distortion in cognitive mapping abilities?

Rats were treated postnatally (PND 5-16) with BSO (l-buthionine-(S,R)-sulfoximine) in an animal model of schizophrenia based on transient glutathione deficit. The BSO treated rats were impaired in patrolling a maze or a homing table when adult, yet demonstrated preserved escape learning, place discrimination and reversal in a water maze task [37]. In the present work, BSO rats' performance in the water maze was assessed in conditions controlling for the available visual cues. First, in a completely curtained environment with two salient controlled cues, BSO rats showed little accuracy compared to control rats. Secondly, pre-trained BSO rats were impaired in reaching the familiar spatial position when curtains partially occluded different portions of the room environment in successive sessions. The apparently preserved place learning in a classical water maze task thus appears to require the stability and the richness of visual landmarks from the surrounding environment. In other words, the accuracy of BSO rats in place and reversal learning is impaired in a minimal cue condition or when the visual panorama changes between trials. However, if the panorama remains rich and stable between trials, BSO rats are equally efficient in reaching a familiar position or in learning a new one. This suggests that the BSO accurate performance in the water maze does not satisfy all the criteria for a cognitive map based navigation on the integration of polymodal cues. It supports the general hypothesis of a binding deficit in BSO rats.

3D geoelectrical monitoring to follow hydric soil behaviour during rainfall

RESUME Durant les dernières années, les méthodes électriques ont souvent été utilisées pour l'investigation des structures de subsurface. L'imagerie électrique (Electrical Resistivity Tomography, ERT) est une technique de prospection non-invasive et spatialement intégrée. La méthode ERT a subi des améliorations significatives avec le développement de nouveaux algorithmes d'inversion et le perfectionnement des techniques d'acquisition. La technologie multicanale et les ordinateurs de dernière génération permettent la collecte et le traitement de données en quelques heures. Les domaines d'application sont nombreux et divers: géologie et hydrogéologie, génie civil et géotechnique, archéologie et études environnementales. En particulier, les méthodes électriques sont souvent employées dans l'étude hydrologique de la zone vadose. Le but de ce travail est le développement d'un système de monitorage 3D automatique, non- invasif, fiable, peu coûteux, basé sur une technique multicanale et approprié pour suivre les variations de résistivité électrique dans le sous-sol lors d'événements pluvieux. En raison des limitations techniques et afin d'éviter toute perturbation physique dans la subsurface, ce dispositif de mesure emploie une installation non-conventionnelle, où toutes les électrodes de courant sont placées au bord de la zone d'étude. Le dispositif le plus approprié pour suivre les variations verticales et latérales de la résistivité électrique à partir d'une installation permanente a été choisi à l'aide de modélisations numériques. Les résultats démontrent que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle et plus apte à détecter les variations latérales et verticales de la résistivité électrique, et cela malgré la configuration non-conventionnelle des électrodes. Pour tester l'efficacité du système proposé, des données de terrain ont été collectées sur un site d'étude expérimental. La technique de monitorage utilisée permet de suivre le processus d'infiltration 3D pendant des événements pluvieux. Une bonne corrélation est observée entre les résultats de modélisation numérique et les données de terrain, confirmant par ailleurs que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle. La nouvelle technique de monitorage 3D de résistivité électrique permet de caractériser les zones d'écoulement préférentiel et de caractériser le rôle de la lithologie et de la pédologie de manière quantitative dans les processus hydrologiques responsables d'écoulement de crue. ABSTRACT During the last years, electrical methods were often used for the investigation of subsurface structures. Electrical resistivity tomography (ERT) has been reported to be a useful non-invasive and spatially integrative prospecting technique. The ERT method provides significant improvements, with the developments of new inversion algorithms, and the increasing efficiency of data collection techniques. Multichannel technology and powerful computers allow collecting and processing resistivity data within few hours. Application domains are numerous and varied: geology and hydrogeology, civil engineering and geotechnics, archaeology and environmental studies. In particular, electrical methods are commonly used in hydrological studies of the vadose zone. The aim of this study was to develop a multichannel, automatic, non-invasive, reliable and inexpensive 3D monitoring system designed to follow electrical resistivity variations in soil during rainfall. Because of technical limitations and in order to not disturb the subsurface, the proposed measurement device uses a non-conventional electrode set-up, where all the current electrodes are located near the edges of the survey grid. Using numerical modelling, the most appropriate arrays were selected to detect vertical and lateral variations of the electrical resistivity in the framework of a permanent surveying installation system. The results show that a pole-dipole array has a better resolution than a pole-pole array and can successfully follow vertical and lateral resistivity variations despite the non-conventional electrode configuration used. Field data are then collected at a test site to assess the efficiency of the proposed monitoring technique. The system allows following the 3D infiltration processes during a rainfall event. A good correlation between the results of numerical modelling and field data results can be observed since the field pole-dipole data give a better resolution image than the pole-pole data. The new device and technique makes it possible to better characterize the zones of preferential flow and to quantify the role of lithology and pedology in flood- generating hydrological processes.

Derivation of soil-specific streaming potential electrical parameters from hydrodynamic characteristics of partially saturated soils

Water movement in unsaturated soils gives rise to measurable electrical potential differences that are related to the flow direction and volumetric fluxes, as well as to the soil properties themselves. Laboratory and field data suggest that these so-called streaming potentials may be several orders of magnitudes larger than theoretical predictions that only consider the influence of the relative permeability and electrical conductivity on the self potential (SP) data. Recent work has improved predictions somewhat by considering how the volumetric excess charge in the pore space scales with the inverse of water saturation. We present a new theoretical approach that uses the flux-averaged excess charge, not the volumetric excess charge, to predict streaming potentials. We present relationships for how this effective excess charge varies with water saturation for typical soil properties using either the water retention or the relative permeability function. We find large differences between soil types and the predictions based on the relative permeability function display the best agreement with field data. The new relationships better explain laboratory data than previous work and allow us to predict the recorded magnitudes of the streaming potentials following a rainfall event in sandy loam, whereas previous models predict values that are three orders of magnitude too small. We suggest that the strong signals in unsaturated media can be used to gain information about fluxes (including very small ones related to film flow), but also to constrain the relative permeability function, the water retention curve, and the relative electrical conductivity function.

Burkholderia sartisoli sp. nov., isolated from a polycyclic aromatic hydrocarbon-contaminated soil.

A Gram-negative, rod-shaped, aerobic bacterium, designated strain RP007(T), was isolated from a polycyclic aromatic hydrocarbon-contaminated soil in New Zealand. Two additional strains were recovered from a compost heap in Belgium (LMG 18808) and from the rhizosphere of maize in the Netherlands (LMG 24204). The three strains had virtually identical 16S rRNA gene sequences and whole-cell protein profiles, and they were identified as members of the genus Burkholderia, with Burkholderia phenazinium as their closest relative. Strain RP007(T) had a DNA G+C content of 63.5 mol% and could be distinguished from B. phenazinium based on a range of biochemical characteristics. Strain RP007(T) showed levels of DNA-DNA relatedness towards the type strain of B. phenazinium and those of other recognized Burkholderia species of less than 30 %. The results of 16S rRNA gene sequence analysis, DNA-DNA hybridization experiments and physiological and biochemical tests allowed the differentiation of strain RP007(T) from all recognized species of the genus Burkholderia. Strains RP007(T), LMG 18808 and LMG 24204 are therefore considered to represent a single novel species of the genus Burkholderia, for which the name Burkholderia sartisoli sp. nov. is proposed. The type strain is RP007(T) (=LMG 24000(T) =CCUG 53604(T) =ICMP 13529(T)).