In order to study the distribution of pesticide residues in beekeeping matrices, samples of live in-hive worker honey bees (Apis mellifera), fresh stored pollen and beeswax were collected during 2016–2017 from 45 apiaries located in different landscape contexts in Spain. A total of 133 samples were screened for 63 pesticides or their degradation products to estimate the pesticide exposure to honey bee health through the calculation of the hazard quotient (HQ). The influence of the surrounding environment on the content of pesticides in pollen was assessed by comparing the concentrations of pesticide residues found in apiaries from intensive farming landscapes to those found in apiaries located in mountainous, grassland and urban contexts. Beeswax revealed high levels of miticides used in beekeeping such as coumaphos, chlorfenvinphos, fluvalinate and acrinathrin, which were detected in more than 75% of samples. Pollen was predominantly contaminated by miticides but also by insecticides used in agriculture such as chlorpyrifos and acetamiprid, which showed concentrations significantly higher in apiaries located in intensive farming contexts. Pesticides residues were less frequent and at lower concentrations in live honey bees. Beeswax showed the highest average hazard scores (HQ > 5000) to honey bees. Pollen samples contained the largest number of pesticide residues and relevant hazard (HQ > 50) to bees. Acrinathrin was the most important contributor to the hazard quotient scores in wax and pollen samples. The contributions of the pesticides dimethoate and chlorpyrifos to HQ were considered relevant in samples.

In order to study the distribution of pesticide residues in beekeeping matrices, samples of live in-hive worker honey bees (Apis mellifera), fresh stored pollen and beeswax were collected during 2016–2017 from 45 apiaries located in different landscape contexts in Spain. A total of 133 samples were screened for 63 pesticides or their degradation products to estimate the pesticide exposure to honey bee health through the calculation of the hazard quotient (HQ). The influence of the surrounding environment on the content of pesticides in pollen was assessed by comparing the concentrations of pesticide residues found in apiaries from intensive farming landscapes to those found in apiaries located in mountainous, grassland and urban contexts. Beeswax revealed high levels of miticides used in beekeeping such as coumaphos, chlorfenvinphos, fluvalinate and acrinathrin, which were detected in more than 75% of samples. Pollen was predominantly contaminated by miticides but also by insecticides used in agriculture such as chlorpyrifos and acetamiprid, which showed concentrations significantly higher in apiaries located in intensive farming contexts. Pesticides residues were less frequent and at lower concentrations in live honey bees. Beeswax showed the highest average hazard scores (HQ > 5000) to honey bees. Pollen samples contained the largest number of pesticide residues and relevant hazard (HQ > 50) to bees. Acrinathrin was the most important contributor to the hazard quotient scores in wax and pollen samples. The contributions of the pesticides dimethoate and chlorpyrifos to HQ were considered relevant in samples. | This work has been suported by the Spanish Ministerio de Agricultura, Alimentación y Medio Ambiente through the project 20160020000834: “Análisis de la Pérdida de Viabilidad y Despoblamiento de las Colonias de Abejas (Apis mellifera L.) Mediante el Mapeo de Xenobióticos”. Pau Calatayud-Vernich acknowledges the Foundation “Tatiana Pérez de Guzmán el Bueno” for the financial support through PhD Environmental Fellowship Program. We thank to the Agrupación de Defensa Sanitaria Apicola (apiADS) and the particular beekeepers for their help with counseling and sampling. | Peer reviewed

Vinclozolin (Vz) is a pollutant found in aquatic environments whose antiandrogenic effects in reproduction are well known in mammals. Although its reproductive effects have been less studied in invertebrates, other effects, including genotoxicity, have been described. Therefore, in this work, we studied the genotoxic effects of Vz in the freshwater benthic invertebrate Chironomus riparius. DNA damage was evaluated with the comet assay (tail area, olive moment, tail moment and % DNA in tail), and the transcriptional levels of different genes involved in DNA repair (ATM, NLK and XRCC1) and apoptosis (DECAY) were measured by RT-PCR. Fourth instar larvae of C. riparius, were exposed to Vz for 24 h at 20 and 200 μg/L. The Vz exposures affected the DNA integrity in this organism, since a dose-response relationship occurred, with DNA strand breaks significantly increased with increased dose for tail area, olive moment and tail moment parameters. Additionally, the lower concentration of Vz produced a significant induction of the transcripts of three genes under study (ATM, NLK and XRCC1) showing the activation of the cellular repair mechanism. In contrast, the expression of these genes with the highest concentration were downregulated, indicating failure of the cellular repair mechanism, which would explain the higher DNA damage. These data report for the first time the alterations of Vz on gene transcription of an insect and confirm the potential genotoxicity of this compound on freshwater invertebrates.

Vinclozolin (Vz) is a pollutant found in aquatic environments whose antiandrogenic effects in reproduction are well known in mammals. Although its reproductive effects have been less studied in invertebrates, other effects, including genotoxicity, have been described. Therefore, in this work, we studied the genotoxic effects of Vz in the freshwater benthic invertebrate Chironomus riparius. DNA damage was evaluated with the comet assay (tail area, olive moment, tail moment and % DNA in tail), and the transcriptional levels of different genes involved in DNA repair (ATM, NLK and XRCC1) and apoptosis (DECAY) were measured by RT-PCR. Fourth instar larvae of C. riparius, were exposed to Vz for 24 h at 20 and 200 μg/L. The Vz exposures affected the DNA integrity in this organism, since a dose-response relationship occurred, with DNA strand breaks significantly increased with increased dose for tail area, olive moment and tail moment parameters. Additionally, the lower concentration of Vz produced a significant induction of the transcripts of three genes under study (ATM, NLK and XRCC1) showing the activation of the cellular repair mechanism. In contrast, the expression of these genes with the highest concentration were downregulated, indicating failure of the cellular repair mechanism, which would explain the higher DNA damage. These data report for the first time the alterations of Vz on gene transcription of an insect and confirm the potential genotoxicity of this compound on freshwater invertebrates.

Natural levels of heavy metals (HM) have increased during the industrial era to the point of posing a serious threat to the environment. The use of tree species to record contamination is a well-known practice. The objective of the study was to compare HM levels under different pollution conditions: a) soil pollution due to mining waste; b) atmospheric pollution due to coal-fired power plant emissions. We report significant HM enrichment in Pinus halepensis tissues. Near a burning power plant, Pb content in a tree wood was 2.5-fold higher that in natural areas (no pollution; NP). In mining areas, Cd content was 25-fold higher than NP. The hypothesis that HM contents in tree rings should register pollution is debatable. HM uptake by pines from soil, detoxification mechanisms and resuspended local soil dust is involved in HM contents in wood and bark.

Natural levels of heavy metals (HM) have increased during the industrial era to the point of posing a serious threat to the environment. The use of tree species to record contamination is a well-known practice. The objective of the study was to compare HM levels under different pollution conditions: a) soil pollution due to mining waste; b) atmospheric pollution due to coal-fired power plant emissions. We report significant HM enrichment in Pinus halepensis tissues. Near a burning power plant, Pb content in a tree wood was 2.5-fold higher that in natural areas (no pollution; NP). In mining areas, Cd content was 25-fold higher than NP. The hypothesis that HM contents in tree rings should register pollution is debatable. HM uptake by pines from soil, detoxification mechanisms and resuspended local soil dust is involved in HM contents in wood and bark.

Mercury (Hg) is a global contaminant of environmental concern. Numerous factors influencing its bioaccumulation in marine organisms have already been described at both individual and species levels (e.g., size or age, habitat, trophic level). However, few studies have compared the trophic characteristics of ecosystems to explain underlying mechanisms of differences in Hg bioaccumulation and biomagnification among food webs and systems. The present study aimed at investigating the potential primary role of the trophic status of systems on Hg bioaccumulation and biomagnification in temperate marine food webs, as shown by their medium-to high-trophic level consumers. It used data from samples collected at the shelf-edge (i.e. offshore organisms) in two contrasted ecosystems: the Bay of Biscay in the North-East Atlantic Ocean and the Gulf of Lion in the North-West Mediterranean Sea. Seven species including crustaceans, sharks and teleost fish, previously analysed for their total mercury (T-Hg) concentrations and their stable carbon and nitrogen isotope compositions, were considered for a meta-analysis. In addition, methylated mercury forms (or methyl-mercury, Me-Hg) were analysed. Mediterranean organisms presented systematically lower sizes than Atlantic ones, and lower δ¹³C and δ¹⁵N values, the latter values especially highlighting the more oligotrophic character of Mediterranean waters. Mediterranean individuals also showed significantly higher T-Hg and Me-Hg concentrations. Conversely, Me-Hg/T-Hg ratios were higher than 85% for all species, and quite similar between systems. Finally, the biomagnification power of Hg was different between systems when considering T-Hg, but not when considering Me-Hg, and was not different between the Hg forms within a given system. Overall, the different parameters showed the crucial role of the low primary productivity and its effects rippling through the compared ecosystems in the higher Hg bioaccumulation seen in organisms from oligotrophic Mediterranean waters.

Mercury (Hg) is a global contaminant of environmental concern. Numerous factors influencing its bioaccumulation in marine organisms have already been described at both individual and species levels (e.g., size or age, habitat, trophic level). However, few studies have compared the trophic characteristics of ecosystems to explain underlying mechanisms of differences in Hg bioaccumulation and biomagnification among food webs and systems. The present study aimed at investigating the potential primary role of the trophic status of systems on Hg bioaccumulation and biomagnification in temperate marine food webs, as shown by their medium-to high-trophic level consumers. It used data from samples collected at the shelf-edge (i.e. offshore organisms) in two contrasted ecosystems: the Bay of Biscay in the North-East Atlantic Ocean and the Gulf of Lion in the North-West Mediterranean Sea. Seven species including crustaceans, sharks and teleost fish, previously analysed for their total mercury (T-Hg) concentrations and their stable carbon and nitrogen isotope compositions, were considered for a meta-analysis. In addition, methylated mercury forms (or methyl-mercury, Me-Hg) were analysed. Mediterranean organisms presented systematically lower sizes than Atlantic ones, and lower δ13C and δ15N values, the latter values especially highlighting the more oligotrophic character of Mediterranean waters. Mediterranean individuals also showed significantly higher T-Hg and Me-Hg concentrations. Conversely, Me-Hg/T-Hg ratios were higher than 85% for all species, and quite similar between systems. Finally, the biomagnification power of Hg was different between systems when considering T-Hg, but not when considering Me-Hg, and was not different between the Hg forms within a given system. Overall, the different parameters showed the crucial role of the low primary productivity and its effects rippling through the compared ecosystems in the higher Hg bioaccumulation seen in organisms from oligotrophic Mediterranean waters.

Erythromycin is an antibiotic employed in the treatment of infections caused by Gram positive microorganisms and the increasing use has made it a contaminant of emerging concern in aqueous ecosystems. Cerium oxide nanoparticles (CeO₂ NPs), which are known to have catalytic and antioxidant properties, have also become contaminants of emerging concern. Due to the high reactivity of CeO₂ NPs, they can interact with erythromycin magnifying their effects or on the other hand, considering the redox potential of CeO₂ NPs, it can alleviate the toxicity of erythromycin. The present study was carried out to assess the toxicity of both single compounds as well as mixed on Chlamydomonas reinhardtii and Phaeodactylum tricornutum (freshwater and marine microalgae respectively) employed as target species in ecotoxicological tests. Mechanisms of oxidative damage and those harmful to the photosynthetic apparatus were studied in order to know the toxic mechanisms of erythromycin and the joint effects with CeO₂ NPs. Results showed that erythromycin inhibited the microalgae population growth and effective quantum yield of PSII (E.Q.Y.) in both microalgae. However, the freshwater microalgae Chlamydomonas reinhardtii was more sensitive than the marine diatom Phaeodactylum tricornutum. Responses related to the photosynthetic apparatus such as E.Q.Y. was affected by the exposure to erythromycin of both microalgae, as chloroplasts are target organelle for this antibiotic.Mixed experiments (CeO₂ NPs + erythromycin) showed the protective role of CeO₂ NPs in both microalgae preventing erythromycin toxicity in toxicological responses such as the growth of the microalgae population and E.Q.Y.

Erythromycin is an antibiotic employed in the treatment of infections caused by Gram positive microorganisms and the increasing use has made it a contaminant of emerging concern in aqueous ecosystems. Cerium oxide nanoparticles (CeO2 NPs), which are known to have catalytic and antioxidant properties, have also become contaminants of emerging concern. Due to the high reactivity of CeO2 NPs, they can interact with erythromycin magnifying their effects or on the other hand, considering the redox potential of CeO2 NPs, it can alleviate the toxicity of erythromycin. The present study was carried out to assess the toxicity of both single compounds as well as mixed on Chlamydomonas reinhardtii and Phaeodactylum tricornutum (freshwater and marine microalgae respectively) employed as target species in ecotoxicological tests. Mechanisms of oxidative damage and those harmful to the photosynthetic apparatus were studied in order to know the toxic mechanisms of erythromycin and the joint effects with CeO2 NPs. Results showed that erythromycin inhibited the microalgae population growth and effective quantum yield of PSII (E.Q.Y.) in both microalgae. However, the freshwater microalgae Chlamydomonas reinhardtii was more sensitive than the marine diatom Phaeodactylum tricornutum. Responses related to the photosynthetic apparatus such as E.Q.Y. was affected by the exposure to erythromycin of both microalgae, as chloroplasts are target organelle for this antibiotic. Mixed experiments (CeO2 NPs + erythromycin) showed the protective role of CeO2 NPs in both microalgae preventing erythromycin toxicity in toxicological responses such as the growth of the microalgae population and E.Q.Y. | This research has been funded by the Junta de Andalucía (PE2011-RNM-7812 project and FQM-110 group) and the Spanish National Research Plan (CTM2016-75908-R and Project MINECO/FEDER MAT2013-40823-R). We are also grateful to Pilar Yeste and José Manuel Gatica from SCICYT of Cadiz University (UCA) for the use of its Electron Microscopy division facilities and characterization of NPs. Furthermore, we are grateful to Catalina Fernández from IFAPA for her support in the characterization of NPs through DLS. CVM Araújo is also grateful to the Juan de la Cierva-Incorporación contract (IJCI-2014-19318) provided by the Spanish Ministry of Economy and Competitiveness. | Peer reviewed

The Mediterranean Sea has been described as one of the most affected areas by marine litter in the world. Although effects on organisms from marine plastic litter ingestion have been investigated in several oceanic areas, there is still a lack of information from the Mediterranean Sea. The main objectives of this paper are to review current knowledge on the impact of marine litter on Mediterranean biodiversity, to define selection criteria for choosing marine organisms suitable for use as bioindicator species, and to propose a methodological approach to assessing the harm related to marine litter ingestion in several Mediterranean habitats and sub-regions. A new integrated monitoring tool that would provide the information necessary to design and implement future mitigation actions in the Mediterranean basin is proposed.According to bibliographic research and statistical analysis on current knowledge of marine litter ingestion, the area of the Mediterranean most studied, in terms of number of species and papers in the Mediterranean Sea is the western sub-area as well as demersal (32.9%) and pelagic (27.7%) amongst habitats.Applying ecological and biological criteria to the most threatened species obtained by statistical analysis, bioindicator species for different habitats and monitoring scale were selected. A threefold approach, simultaneously measuring the presence and effects of plastic, can provide the actual harm and sub-lethal effects to organisms caused by marine litter ingestion. The research revealed gaps in knowledge, and this paper suggests measures to close the gap. This and the selection of appropriate bioindicator species would represent a step forward for marine litter risk assessment, and the implementation of future actions and mitigation measures for specific Mediterranean areas, habitats and species affected by marine litter ingestion.

The Mediterranean Sea has been described as one of the most affected areas by marine litter in the world. Although effects on organisms from marine plastic litter ingestion have been investigated in several oceanic areas, there is still a lack of information from the Mediterranean Sea. The main objectives of this paper are to review current knowledge on the impact of marine litter on Mediterranean biodiversity, to define selection criteria for choosing marine organisms suitable for use as bioindicator species, and to propose a methodological approach to assessing the harm related to marine litter ingestion in several Mediterranean habitats and sub-regions. A new integrated monitoring tool that would provide the information necessary to design and implement future mitigation actions in the Mediterranean basin is proposed. According to bibliographic research and statistical analysis on current knowledge of marine litter ingestion, the area of the Mediterranean most studied, in terms of number of species and papers in the Mediterranean Sea is the western sub-area as well as demersal (32.9%) and pelagic (27.7%) amongst habitats. Applying ecological and biological criteria to the most threatened species obtained by statistical analysis, bioindicator species for different habitats and monitoring scale were selected. A threefold approach, simultaneously measuring the presence and effects of plastic, can provide the actual harm and sub-lethal effects to organisms caused by marine litter ingestion. The research revealed gaps in knowledge, and this paper suggests measures to close the gap. This and the selection of appropriate bioindicator species would represent a step forward for marine litter risk assessment, and the implementation of future actions and mitigation measures for specific Mediterranean areas, habitats and species affected by marine litter ingestion | This work was supported by the Deutsche Gesellschaft für Internationale Zusammenarbeit (GIZ), Grant number: 13.6257.3-002.00–81197130, within the Plastic Busters project-SDSN Solution Initiative. The authors would like to thank the Union for the Mediterranean and the UN Sustainable Development and Solution Network-Mediterranean that labelled and endorsed the project and the resulting paper. | Peer reviewed

Several works evaluated the toxicity of pharmaceutical drugs and climate related changes in invertebrates but few explored the combined effects of both stressors, namely considering their mode of action (MoA). Carbamazepine (CBZ) and cetirizine (CTZ) are pharmaceutical drugs detected in the environment and the toxicity derived from the combined effects of these drugs with ocean acidification (OA) is poorly explored. Thus, the present study investigated the biochemical parameters related to an oxidative stress response and the transcription of genes related to the MoA of CBZ (1.0 μg/L) and CTZ (0.6 μg/L) in the clam Ruditapes philippinarum chronically exposed (28 days) to control (7.8) and low (7.5) pH conditions. The results obtained showed that despite the clams accumulated both drugs, at low pH the clams exposed to CTZ decreased drug concentration and BCF values (CTZ uptake: 2.0 ± 0.5 ng/g fresh weight; BCF: 3.8 ± 0.9) in comparison with clams exposed to control pH (CTZ uptake: 2.9 ± 0.3 ng/g fresh weight; BCF: 5.5 ± 0.6). No oxidative stress was induced by the exposure to CBZ or CTZ at each pH level, but the transcription of several genes related with the MoA (neurotransmission, immunity and biomineralization) was altered by low pH, drug exposure and the combination of both stressors. At both pH conditions, CBZ increased the transcription of GABA receptor gene (neurotransmission) and CTZ led to a decrease of Perlucin gene (biomineralization) transcription. The transcription of MyD88 gene (immunity) decreased at low pH (7.5) combined with drug exposure (CBZ or CTZ). Thus, it was highlighted that the interaction of drug exposure and low pH conditions can change bivalves’ sensitivity to drugs or alter drugs toxicity.

Several works evaluated the toxicity of pharmaceutical drugs and climate related changes in invertebrates but few explored the combined effects of both stressors, namely considering their mode of action (MoA). Carbamazepine (CBZ) and cetirizine (CTZ) are pharmaceutical drugs detected in the environment and the toxicity derived from the combined effects of these drugs with ocean acidification (OA) is poorly explored. Thus, the present study investigated the biochemical parameters related to an oxidative stress response and the transcription of genes related to the MoA of CBZ (1.0 μg/L) and CTZ (0.6 μg/L) in the clam Ruditapes philippinarum chronically exposed (28 days) to control (7.8) and low (7.5) pH conditions. The results obtained showed that despite the clams accumulated both drugs, at low pH the clams exposed to CTZ decreased drug concentration and BCF values (CTZ uptake: 2.0 ± 0.5 ng/g fresh weight; BCF: 3.8 ± 0.9) in comparison with clams exposed to control pH (CTZ uptake: 2.9 ± 0.3 ng/g fresh weight; BCF: 5.5 ± 0.6). No oxidative stress was induced by the exposure to CBZ or CTZ at each pH level, but the transcription of several genes related with the MoA (neurotransmission, immunity and biomineralization) was altered by low pH, drug exposure and the combination of both stressors. At both pH conditions, CBZ increased the transcription of GABA receptor gene (neurotransmission) and CTZ led to a decrease of Perlucin gene (biomineralization) transcription. The transcription of MyD88 gene (immunity) decreased at low pH (7.5) combined with drug exposure (CBZ or CTZ). Thus, it was highlighted that the interaction of drug exposure and low pH conditions can change bivalves’ sensitivity to drugs or alter drugs toxicity. | published

Several works evaluated the toxicity of pharmaceutical drugs and climate related changes in invertebrates but few explored the combined effects of both stressors, namely considering their mode of action (MoA). Carbamazepine (CBZ) and cetirizine (CTZ) are pharmaceutical drugs detected in the environment and the toxicity derived from the combined effects of these drugs with ocean acidification (OA) is poorly explored. Thus, the present study investigated the biochemical parameters related to an oxidative stress response and the transcription of genes related to the MoA of CBZ (1.0 μg/L) and CTZ (0.6 μg/L) in the clam Ruditapes philippinarum chronically exposed (28 days) to control (7.8) and low (7.5) pH conditions. The results obtained showed that despite the clams accumulated both drugs, at low pH the clams exposed to CTZ decreased drug concentration and BCF values (CTZ uptake: 2.0 ± 0.5 ng/g fresh weight; BCF: 3.8 ± 0.9) in comparison with clams exposed to control pH (CTZ uptake: 2.9 ± 0.3 ng/g fresh weight; BCF: 5.5 ± 0.6). No oxidative stress was induced by the exposure to CBZ or CTZ at each pH level, but the transcription of several genes related with the MoA (neurotransmission, immunity and biomineralization) was altered by low pH, drug exposure and the combination of both stressors. At both pH conditions, CBZ increased the transcription of GABA receptor gene (neurotransmission) and CTZ led to a decrease of Perlucin gene (biomineralization) transcription. The transcription of MyD88 gene (immunity) decreased at low pH (7.5) combined with drug exposure (CBZ or CTZ). Thus, it was highlighted that the interaction of drug exposure and low pH conditions can change bivalves' sensitivity to drugs or alter drugs toxicity.

The concentrations of Cd, Cr and Pb in soil samples and As, Cd, Cr and Pb in plant specimens were analyzed in an arid area in central Iran. Plants were categorized into desert-adapted (Haloxylon ammodendron, Atraphaxis spinosa and Artemisia persica) and non-desert species. It was found that the trace element (TE) accumulating potential of the desert species (Haloxylon ammodendron and Artemisia persica) with a mean value of 0.1 mg kg⁻¹ for Cd was significantly higher than that of the majority of the non-desert species with an average of 0.05 mg kg⁻¹. Artemisia also had a high As accumulating capability with a mean level of 0.8 mg kg⁻¹ in comparison with an average of 0.2 mg kg⁻¹ for most of the other plant species. The mean values of Cr and Pb in Haloxylon ammodendron and Artemisia persica were 5 and 3 mg kg⁻¹, respectively. Among the desert-adapted plants, Atraphaxis proved to be a species with high Cr and Pb accumulating potential, as well. The geoaccumulation index and the overall pollution scores indicated that the highest environmental risk was related to Cd. Different statistical analyses were used to study the spatial patterns of soil Cd and their connections with pollution sources. The variogram was estimated using a classical approach (weighted least squares) and was compared with that of the posterior summaries that resulted from the Bayesian technique, which lay within the 95% Bayesian credible quantile intervals (BIC) of posterior parameter distributions. The prediction of cadmium values at un-sampled locations was implemented by multi-Gaussian kriging and sequential Gaussian simulation methods. The prediction maps showed that the region most contaminated by Cd was the north-eastern part of the study area, which was linked to mining activities, while agricultural influence contributed less in this respect.

The concentrations of Cd, Cr and Pb in soil samples and As, Cd, Cr and Pb in plant specimens were analyzed in an arid area in central Iran. Plants were categorized into desert-adapted (Haloxylon ammodendron, Atraphaxis spinosa and Artemisia persica) and non-desert species. It was found that the trace element (TE) accumulating potential of the desert species (Haloxylon ammodendron and Artemisia persica) with a mean value of 0.1 mg kg−1 for Cd was significantly higher than that of the majority of the non-desert species with an average of 0.05 mg kg−1. Artemisia also had a high As accumulating capability with a mean level of 0.8 mg kg−1 in comparison with an average of 0.2 mg kg−1 for most of the other plant species. The mean values of Cr and Pb in Haloxylon ammodendron and Artemisia persica were 5 and 3 mg kg−1, respectively. Among the desert-adapted plants, Atraphaxis proved to be a species with high Cr and Pb accumulating potential, as well. The geoaccumulation index and the overall pollution scores indicated that the highest environmental risk was related to Cd. Different statistical analyses were used to study the spatial patterns of soil Cd and their connections with pollution sources. The variogram was estimated using a classical approach (weighted least squares) and was compared with that of the posterior summaries that resulted from the Bayesian technique, which lay within the 95% Bayesian credible quantile intervals (BIC) of posterior parameter distributions. The prediction of cadmium values at un-sampled locations was implemented by multi-Gaussian kriging and sequential Gaussian simulation methods. The prediction maps showed that the region most contaminated by Cd was the north-eastern part of the study area, which was linked to mining activities, while agricultural influence contributed less in this respect.