Estelle Villanneau

Which persistent organic pollutants can we map in soil using a large spacing systematic soil monitoring design? A case study in Northern France.

Persistent organic pollutants (POPs) impact upon human and animal health and the wider environment. It is important to determine where POPs are found and the spatial pattern of POP variation. The concentrations of 90 molecules which are members of four families of POPs and two families of herbicides were measured within a region of Northern France as part of the French National Soil Monitoring Network (RMQS: Réseau de Mesures de la Qualité des Sols). We also gather information on five covariates (elevation, soil organic carbon content, road density, land cover and population density) which might influence POP concentrations. The study region contains 105 RMQS observation sites arranged on a regular square grid with spacing of 16 km. The observations include hot-spots at sites of POP application, smaller concentrations where POPs have been dispersed and observations less than the limit of quantification (LOQ) where the soil has not been impacted by POPs. Fifty nine of the molecules were detected at less than 50 sites and hence the data were unsuitable for spatial analyses. We represent the variation of the remaining 31 molecules by various linear mixed models which can include fixed effects (i.e. linear relationships between the molecule concentrations and covariates) and spatially correlated random effects. The best model for each molecule is selected by the Akaike Information Criterion. For nine of the molecules, spatial correlation is evident and hence they can potentially be mapped. For four of these molecules, the spatial correlation cannot be wholly explained by fixed effects. It appears that these molecules have been transported away from their application sites and are now dispersed across the study region with the largest concentrations found in a heavily populated depression. More complicated statistical models and sampling designs are required to explain the distribution of the less dispersed molecules.

Spatial distribution of lindane in topsoil of Northern France.

Lindane is a persistent organochlorine insecticide and the use of this insecticide in agriculture was banned in France in 1998. In this study we investigated the concentrations of lindane in top soil in Northern France and used robust geostatistics to map the geographical distribution of lindane. The study was based on a 16 km x 16 km grid covering an area of ca 25,000 km(2). Lindane was found in all soils, even those from non-agricultural-application areas. Very low ratios of alpha-/gamma-HCH and delta-/gamma-HCH suggested that a long time had passed since technical HCH was used in the studied area, or that emission sources of lindane were still present. A strong gradient in lindane concentration was observed, with the highest lindane concentrations in an area located in the northern region. Results suggested that some of the lindane observed in the high concentration area may have come from volatilization of old lindane applied to intensively cultivated areas, which was then transported by prevailing winds coming from the south-west and deposited in a densely inhabited depression.

First evidence of large-scale PAH trends in French soils

Polycyclic aromatic hydrocarbons (PAHs) are widespread organic pollutants. Soils are a reservoir of PAHs because some soil constituents favour PAH accumulation. Therefore, soil is a key indicator of the degree of contamination. So far, studies mapping soil PAH levels over large territories are very rare. Here, we report the first nation-wide maps of soil PAHs in France. Results were obtained within the French National Soil Monitoring Network, which is the first European network monitoring systematically soil PAHs. We used advanced geostatistics to map PAH distribution over the whole French territory. Our results show clear trends of PAH levels at the nation scale. For instance, the highest PAH levels are found in Northern and Eastern France. This high contamination is explained by the intense industrial activity of these regions during the last century. High levels of PAH are also found near some coastlines. This observation could be explained by long-range atmospheric transportation. In addition, we found that light PAHs are rarely found in French topsoils.

Spatial distribution of Lindane concentration in topsoil across France.

Lindane [γ-hexachlorocyclohexane (γ-HCH)] is an organochlorine pesticide with toxic effects on humans. It is bioaccumulative and can remain in soils for long periods, and although its use for crop spraying was banned in France in 1998, it is possible that residues from before this time remain in the soil. The RMQS soil monitoring network consists of soil samples from 2200 sites on a 16 km regular grid across France, collected between 2002 and 2009. We use 726 measurements of the Lindane concentration in these samples to (i) investigate the main explanatory factors for its spatial distribution across France, and (ii) map this distribution. Geostatistics provides an appropriate framework to analyze our spatial dataset, though two issues regarding the data are worth special consideration: first, the harmonization of two subsets of the data (which were analyzed using different measurement processes), and second, the large proportion of data from one of these subsets that fell below a limit of quantification. We deal with these issues using recent methodological developments in geostatistics. Results demonstrate the importance of land use and rainfall for explaining part of the variability of Lindane across France: land use due to the past direct input of Lindane on cropland and its subsequent persistence in the soil, and rainfall due to the re-deposition of volatilized Lindane. Maps show the concentrations to be generally largest in the north and northwest of France, areas of more intensive agricultural land. We also compare levels to some contamination thresholds taken from the literature, and present maps showing the probability of Lindane concentrations exceeding these thresholds across France. These maps could be used as guidelines for deciding which areas require further sampling before some possible remediation strategy could be applied.

Analyzing the spatial distribution of PCB concentrations in soils using below-quantification limit data

Polychlorinated biphenyls (PCBs) are highly toxic environmental pollutants that can accumulate in soils. We consider the problem of explaining and mapping the spatial distribution of PCBs using a spatial data set of 105 PCB-187 measurements from a region in the north of France. A large proportion of our data (35%) fell below a quantification limit (QL), meaning that their concentrations could not be determined to a sufficient degree of precision. Where a measurement fell below this QL, the inequality information was all that we were presented with. In this work, we demonstrate a full geostatistical analysis-bringing together the various components, including model selection, cross-validation, and mapping-using censored data to represent the uncertainty that results from below-QL observations. We implement a Monte Carlo maximum likelihood approach to estimate the geostatistical model parameters. To select the best set of explanatory variables for explaining and mapping the spatial distribution of PCB-187 concentrations, we apply the Akaike Information Criterion (AIC). The AIC provides a trade-off between the goodness-of-fit of a model and its complexity (i.e., the number of covariates). We then use the best set of explanatory variables to help interpolate the measurements via a Bayesian approach, and produce maps of the predictions. We calculate predictions of the probability of exceeding a concentration threshold, above which the land could be considered as contaminated. The work demonstrates some differences between approaches based on censored data and on imputed data (in which the below-QL data are replaced by a value of half of the QL). Cross-validation results demonstrate better predictions based on the censored data approach, and we should therefore have confidence in the information provided by predictions from this method.