Spanish Soils Research Articles

Assessing pesticide residue occurrence and risks in the environment across Europe and Argentina

The widespread and extensive use of pesticides in European crop production to reduce losses from weeds, diseases, and insects may have serious consequences on the ecosystem and human health. This study aimed to identify 20 active substances of high health risk, based on their detection frequency within and across the environmental matrices (soil, crop, water, and sediment) and to identify their associated hazardous effects. A sampling campaign was conducted across 10 case study sites in Europe and 1 in Argentina and included conventional and organic farming systems. In 31% of cases, the detected substances were found at a higher concentration in the soil than in the corresponding crops, 93% of the compounds were fungicides, and the remainder were insecticides. 43% of the substances, 57% of which were insecticides, were detected only in soil. There was a clear relationship between soils and crops in terms of contamination, but not between water and sediment. Portuguese soil (wine grapes) had the highest number of substances (12) with average concentrations (AC) varying between 1 and 162 μg/kg, followed by French (11 substances in wine grapes) (1≤AC≤64 μg/kg) and Spanish soils (9 substances in vegetables) (3≤AC≤59 μg/kg). The crops corresponding to these soils contained a relatively high number of detected substances and several in high average concentrations (AC). The risk quotient was consistently higher for conventional farms than for organic farms. For the soils from conventional farms, 5 active substances (chlorpyrifos, glyphosate, boscalid, difenoconazole, lambda-cyhalothrin, and one metabolite: AMPA) were considered high risk. For water samples, 2 substances (dieldrin and terbuthylazine) found were high risk, and for sediment, there were 3 substances (metalaxyl-M, spiroxamine, and lambda-cyhalothrin). There were 6 substances detected in crops that are suspected to cause human health effects. Uncontaminated soil is a prerequisite for the adoption of sustainable alternatives to pesticides. Efforts are needed to elucidate the unknown effects of mixtures, including biocides and banned compounds in addition to the substances used in agriculture.

Soil origin corresponds with variation in growth of an invasive Centaurea, but not of non-invasive congeners.

Why only a small proportion of exotic species become invasive is an unresolved question. Escape from the negative effects of soil biota in the native range can be important for the success of many invasives, but comparative effects of soil biota on less successful exotic species are poorly understood. Studies of other mechanisms suggest that such comparisons might be fruitful. Seeds of three closely related Centaurea species with overlapping distributions in both their native range of Spain and their nonnative range of California were grown to maturity in pots to obtain an F1 generation of full sibling seeds with reduced maternal effects. Full sibling F1 seeds from both ranges were subsequently grown in pots with inoculations of soil from either the native or nonnative ranges in a fully orthogonal factorial design. We then compared plant biomass among species, regions, and soil sources. Our results indicate that escape from soil pathogens may unleash the highly invasive Centaurea solstitialis, which was suppressed by native Spanish soils but not by soils from California. In contrast, the two non-invasive Centaurea species grew the same on all soils. These results add unprecedented phylogenetically controlled insight into why some species invade and others do not.

Formulation of a Highly Effective Inoculant for Common Bean Based on an Autochthonous Elite Strain of Rhizobium leguminosarum bv. phaseoli, and Genomic-Based Insights Into Its Agronomic Performance.

Common bean is a poor symbiotic N-fixer, with a low response to inoculation owing to its promiscuous nodulation with competitive but inefficient resident rhizobia. Consequently, farmers prefer to fertilize them rather than rely on their capacity for Biological Nitrogen Fixation (BNF). However, when rhizobial inoculants are based on autochthonous strains, they often have superior BNF performance in the field due to their genetic adaptations to the local environment. Nevertheless, there is scant information at the genomic level explaining their superiority or on how their genomes may influence the inoculant performance. This information is especially important in technologically advanced agri-systems like Europe, where environmental concerns and increasingly stringent fertilizer regulations are encouraging a return to the use of rhizobial inoculants, but based upon strains that have been thoroughly characterized in terms of their symbiotic performance and their genetics. The aim of this study was to design an inoculant formulation based on a superior autochthonous strain, Rhizobium leguminosarum bv. phaseoli LCS0306, to assess its performance in the field, and to determine the genomic features contributing to the high effectiveness of its symbiosis with common bean. Plants inoculated with the autochthonous strain LCS0306 fixed significantly more nitrogen than those with the allochthonous strains R. phaseoli ATCC 14482T and R. etli CFN42T, and had grain yield similar to the nitrogen-fertilized controls. Inoculation with LCS0306 was particularly efficacious when formulated with a carrier based upon a mixture of perlite and biochar. Whole genome comparisons revealed no differences in the classical symbiotic genes of strain LCS0306 within the symbiovar phaseoli. However, its symbiotic superior performance might be due to its genomic versatility, as it harbors a large assortment of genes contributing to fitness and competitiveness. It is concluded that inoculation with elite rhizobia formulated with perlite-biochar carriers might constitute a step-change in the sustainable cultivation of common bean in Spanish soils.

Chemometric assessment of soil organic matter storage and quality from humic acid infrared spectra

The knowledge of biogeochemical mechanisms involved in soil organic carbon (SOC) storage is crucial to control its release to the atmosphere. In particular, the chemical composition of soil organic matter (SOM) plays an important role in the performance of the C storage and resilience in soils. The structural information provided by infrared spectroscopy (IR) of soil humic acid (HA) was used in the assessment of the C storage potential of 35 Spanish soils. Partial least squares (PLS) regression using the intensities of the points of the IR spectra of the HAs (4000–400 cm−1) as descriptors shows that a relationship exists between IR spectral pattern and the SOC content. This was also the case for E4 (humification index based on HA optical density at 465 nm). In addition, the chemical characteristics of the HAs correlated with the SOC levels were identified from digital data treatments of the IR spectra. Additional application of principal component analysis (PCA) and multidimensional scaling (MDS) suggested that bands assigned to carboxyl and amide structures were characteristic in HAs from soils with low C content, whereas HA spectra from soils with high C levels showed a conspicuous band pattern suggesting structural units of lignin from slightly transformed plant residues. The spectral profiles were analyzed in detail by an approach based on digital subtraction of IR spectra obtained by averaging those from HAs extracted from soils in the upper and lower quartiles of the SOC distribution. The results showed that significant relationships exist between the molecular composition of HAs and SOC levels and E4 values in a way in which aromatic, carboxyl and amide groups were predominant in HAs from soils with low SOC content, whereas lignin-derived structures were more characteristic of HAs from soils with high SOC content.

Adsorción-desorción de diurón y ametrina en suelos de Colombia y España

Aunque en la actualidad existe una gran preocupación por la contaminación y los posibles efectos sobre la salud que genera el uso frecuente de herbicidas, los estudios ambientales en los países tropicales son escasos, especialmente en Colombia. Por tanto, en el presente trabajo se evaluó la adsorción-desorción de diurón y ametrina en un suelo de Colombia (SC) y de España (SE). Los ensayos de adsorción-desorción se realizaron con moléculas de los herbicidas marcados con 14C, empleando la técnica batch equilibrium. En ambos suelos y para ambos herbicidas hubo un buen ajuste a las isotermas de Freundlich con R2 mayores a 0,99. Para diurón, en SE (Kfa = 12,82 ± 1,23) hubo mayor poder de adsorción frente a SC (4,56 ± 0,26). Por el contrario, para la ametrina se encontró un valor de Kfa = 6,47 ± 0,25 en SC y de 3,05 ± 0,08 en SE, mostrando mayor interacción con SC en comparación con SE. Para diurón se observó un grado de histéresis mayor en SC y para ametrina en SE, mostrando que la interacción de los dos herbicidas con los componentes de cada suelo es diferente dependiendo de su constitución mineralógica y su valor de pH.

Organic Phosphorus Forms in Agricultural Soils under Mediterranean Climate

Core Ideas Monoesters were the dominant OP form. Their degradation occurs under conditions of P limitation. Fe oxide type was crucial for explaining the stabilization of different organic P forms in soils. Inositol hexakisphosphate was mainly related to poorly crystalline Fe oxides. Monoesters different from inositol‐6‐phosphate were related to crystalline Fe oxides and carbonates. Future sustainable phosphorus (P) management in agriculture will require a full understanding of the effects of agricultural practices on soil organic P (OP) speciation in different regions. This work studied the effects of land use, available P status, and soil properties on OP forms in Spanish agricultural soils under Mediterranean climate. We used 31P nuclear magnetic resonance spectroscopy to characterize P in eight soils under two types of land use (field crops and olive or vine orchards), using two samples with different soil Olsen P status (“low” and “high” P samples) per soil type. Orthophosphate monoesters dominated OP in these soils, and inositol hexakisphosphates (IP6) were the most abundant identified group of monoesters. Land use affected total OP concentrations mainly through crop residue incorporation. The contribution of monoesters, including myo‐IP6, to total OP increased with higher Olsen P concentrations, but only in soils from orchards. Degradation of monoesters including IP6 occurred under conditions of P limitation, suggesting some contribution of these OP forms to the plant P supply. Total IP6 was mainly related to poorly crystalline iron (Fe) oxides, whereas other monoesters were related to crystalline Fe oxides and to active Ca carbonate equivalents. The dominant Fe oxide type, carbonates, and the low accumulation of crop residues explained the relative contents of different OP forms in these Spanish soils compared with soils from other geographical areas. The effect of oxides in the accumulation of different OP forms was attributed to the affinities of these forms for oxides and to differences in the density of binding sites depending on the oxide type.

Accumulation of De-Icing Salt and Leaching in Spanish Soils Surrounding Roadways.

The environmental implications of soil salinity caused by accumulation of de-icing salt and leaching in soils of northeastern Spain were examined. For this purpose, the concentrations of ions associated with diagnosing and managing this problem were evaluated from several measurements performed over one year along a road. This analysis demonstrated a higher concentration of soluble Na+ in the soil 3 m from a road in the northernmost part of the study area in February, which made the soil saline-sodic. Data from the rest of the study period (during the spring and summer) demonstrated that the de-icing salt moved to areas farther south by runoff water, which caused environmental impacts by modifying soil characteristics. These results suggest that leaching of Ca2+ and Mg2+ cations occurred faster in the studied systems in sodic soils. Leaching of these cations may affect plant yield, and results in environmental impacts within 3–30 m from the road. Awareness of this impact will be useful for developing future strategies for evaluating and reporting these complex relationships within Spain’s transport system and environment.

Validating proposed migration equation and parameters' values as a tool to reproduce and predict 137Cs vertical migration activity in Spanish soils

This paper shows the procedure performed to validate the migration equation and the migration parameters’ values presented in a previous paper (Legarda et al., 2011) regarding the migration of 137Cs in Spanish mainland soils.In this paper, this model validation has been carried out checking experimentally obtained activity concentration values against those predicted by the model. This experimental data come from the measured vertical activity profiles of 8 new sampling points which are located in northern Spain. Before testing predicted values of the model, the uncertainty of those values has been assessed with the appropriate uncertainty analysis. Once establishing the uncertainty of the model, both activity concentration values, experimental versus model predicted ones, have been compared. Model validation has been performed analyzing its accuracy, studying it as a whole and also at different depth intervals. As a result, this model has been validated as a tool to predict 137Cs behaviour in a Mediterranean environment.

Soil wind erosion in ecological olive trees in the Tabernas desert (southeastern Spain): a wind tunnel experiment

Abstract. Wind erosion is a key component of the soil degradation processes. The purpose of this study is to find out the influence of material loss from wind on soil properties for different soil types and changes in soil properties in olive groves when they are tilled. The study area is located in the north of the Tabernas Desert, in the province of Almería, southeastern Spain. It is one of the driest areas in Europe, with a semiarid thermo-Mediterranean type of climate. We used a new wind tunnel model over three different soil types (olive-cropped Calcisol, Cambisol and Luvisol) and studied micro-plot losses and deposits detected by an integrated laser scanner. We also studied the image processing possibilities for examining the particles attached to collector plates located at the end of the tunnel to determine their characteristics and whether they were applicable to the setup. Samples collected in the traps at the end of the tunnel were analyzed. We paid special attention to the influence of organic carbon, carbonate and clay contents because of their special impact on soil crusting and the wind-erodible fraction. A principal components analysis (PCA) was carried out to find any relations on generated dust properties and the intensity and behavior of those relationships. Component 1 separated data with high N and OC contents from samples high in fine silt, CO3= and available K content. Component 2 separated data with high coarse silt and clay contents from data with high fine sand content. Component 3 was an indicator of available P2O5 content. Analysis of variance (ANOVA) was carried out to analyze the effect of soil type and sampling height on different properties of trapped dust. Calculations based on tunnel data showed overestimation of erosion in soil types and calculation of the fraction of soil erodible by wind done by other authors for Spanish soils. As the highest loss was found in Cambisols, mainly due to the effect on soil crusting and the wind-erodible fraction aggregation of CaCO3, a Stevia rebaudiana cover crop was planted between the rows in this soil type and this favored retention of particles in vegetation.

Electrokinetic-Fenton technology for the remediation of hydrocarbons historically polluted sites

The feasibility of the electrokinetic-Fenton technology coupled with surfactants in the treatment of real historically hydrocarbons polluted soils has been studied. The characterisation of these soils from Spain and Romania was performed and identified as diesel and diesel-motor oil spillages, respectively. Moreover, the ageing of the spillages produced by the soil contamination was estimated showing the historical pollution of the sites (around 11 and 20 years for Romanian and Spanish soils, respectively). An ex-situ electrochemical treatment was performed to evaluate the adequacy of surfactants for the degradation of the hydrocarbons present in the soils. It was found an enhancement in the solubilisation and removal of TPHs with percentages increasing from 25.7 to 81.8% by the presence of Tween 80 for Spanish soil and from 15.1% to 71.6% for Triton X100 in Romanian soil. Therefore, the viability of coupling enhanced electrokinetic and Fenton remediation was evaluated through a simulated in-situ treatment at laboratory scale. The results demonstrated that the addition of the selected surfactants improved the solubilisation of the hydrocarbons and influenced the electroosmotic flow with a slight decrease. The efficiency of the treatment increased for both considered soil samples and a significant degradation level of the hydrocarbons compounds was observed. Buffering of pH coupled with the addition of a complexing agent showed to be important in the treatment process, facilitating the conditions for the degradation reactions that take place into the soil matrix. The results demonstrated the effectiveness of the selected techniques for remediation of the investigated soils.

Factors Controlling Magnetism of Reddish Brown Soil Profiles from Calcarenites in Southern Spain: Dust Input or In-situ Pedogenesis?

Under aerobic conditions, the A and B horizons of soils are magnetically enhanced due to neoformation of ferrimagnets through pedogenesis. This study systematically investigated soils developed on calcarenites of Neogene age in southern Spain to determine the dominant factors controlling the soil magnetism. Geochemical and clay mineral analyses indicate that aeolian dust significantly contribute to the A and B horizon material of the Spanish soil. Nevertheless, the magnetic enhancement of soils can be simply attributed to the pedogenically produced ferrimagnets in-situ. Therefore, the magnetism of Spanish soils is still linked to paleoclimatic variations regardless of the complexities of aeolian inputs from the Northwestern Africa.

Assessment of the soil organic carbon stock in Spain

Soil organic carbon stock (SOCS) assessments at a national level are essential in the climate change mitigation role. Our study provides the first assessment of the SOC stored in entire Spanish surface soils obtained from one sampling protocol within a regular sampling framework in a short edaphic time. In this study, we analysed topsoil samples (0–30cm) from 4401 locations. SOC, soil bulk density and stoniness were measured to estimate SOCS. The results showed that half the Spanish areas obtained a SOC value below 1%. There is inherently wide spatial variability in SOC contents in Spain with low SOC concentration levels located in southern areas. The lowest SOC levels were associated with agricultural soils. However, no statistically significant differences were found between forestlands and grasslands. The mean SOCS in the Spanish topsoil layer was 56.57MgCha−1 and the total stored in topsoil was 2.8 Pg C. We were surprised to find that our estimates were 40% lower than those known to date for Spain. The results herein highlight the potential land use change for SOC sequestration in Spanish soils. Certain changes in agricultural practises (conservation tillage) or converting some unproductive croplands into grasslands or forestlands can increase carbon sequestration in soils.

The Olsen P/solution P relationship as affected by soil properties

AbstractLittle attention has so far been given to the relationship of Olsen P to solution P at P levels in the soil solution within the range relevant to plant growth. In this work, we investigated the relationship between these two P forms in 24 noncalcareous and 25 calcareous Spanish soils, using the concentration of P in the 0.01 m CaCl2 extract (1:10 soil:solution ratio; 30 min shaking) as a proxy for the P concentration in the soil solution. Application of the Freundlich equation (y = axb) in linearized form provided a parsimonious, largely accurate description of the Olsen P/CaCl2‐P relationship. The average value of exponent b, 0.53, suggests that Olsen P relates to the quantity of P in soil rather than to its intensity. The estimated Olsen P values at a CaCl2‐P concentration of 0.01 and 0.03 mg/L (viz., Olsen P0.01 and Olsen P0.03, respectively) varied over wide ranges (0.6–33.9 and 2.0–38.9 mg/kg, respectively). Based on the results of a regression analysis, silicate clays, carbonates and, especially, Fe oxides contribute positively to Olsen P0.01 and Olsen P0.03. The fact that, as shown here, Olsen P at a specific solution P concentration can be predicted from easily measured soil properties can facilitate decision‐making in managing fertilizer P.

Occurrence and analysis of selected pharmaceutical compounds in soil from Spanish agricultural fields

This work describes the analysis of 15 pharmaceutical compounds, belonging to different therapeutic classes (anti-inflammatory/analgesics, lipid regulators, antiepileptics, β-blockers and antidepressants) and with diverse physical-chemical properties, in Spanish soils with different farmland uses. The studied compounds were extracted from soil by ultrasound-assisted extraction (UAE) and determined, after derivatization, by gas chromatography with mass spectrometric detection (GC-MS). The limits of detection (LODs) ranged from 0.14 ng g(-1) (naproxen) to 0.65 ng g(-1) (amitriptyline). At least two compounds where detected in all samples, being ibuprofen, salicylic acid, and paracetamol, the most frequently detected compounds. The highest levels found in soil were 47 ng g(-1) for allopurinol and 37 ng g(-1) for salicylic acid. The influence of the type of crop and the sampling area on the levels of pharmaceuticals in soil, as well as their relationship with soil physical-chemical properties, was studied. The frequent and widespread detection of some of these compounds in agricultural soils show a diffuse contamination, although the low levels found do not pose a risk to the environment or the human health.

Soil carbon stocks and their variability across the forests, shrublands and grasslands of peninsular Spain

Abstract. Accurate estimates of C stocks and fluxes of soil organic carbon (SOC) are needed to assess the impact of climate and land use change on soil C uptake and soil C emissions to the atmosphere. Here, we present an assessment of SOC stocks in forests, shrublands and grasslands of peninsular Spain based on field measurements in more than 900 soil profiles. SOC to a depth of 1 m was modelled as a function of vegetation cover, mean annual temperature, total annual precipitation, elevation and the interaction between temperature and elevation, while latitude and longitude were used to model the correlation structure of the errors. The resulting statistical model was used to estimate SOC in the ∼8 million pixels of the Spanish Forest Map (29.3 × 106 ha). We present what we believe is the most reliable estimation of current SOC in forests, shrublands and grasslands of peninsular Spain thus far, based on the use of spatial modelling, the high number of profiles and the validity and refinement of the data layers employed. Mean concentration of SOC was 8.7 kg m−2, ranging from 2.3 kg m−2 in dry Mediterranean areas to 20.4 kg m−2 in wetter northern locations. This value corresponds to a total stock of 2.544 Tg SOC, which is four times the amount of C estimated to be stored in the biomass of Spanish forests. Climate and vegetation cover were the main variables influencing SOC, with important ecological implications for peninsular Spanish ecosystems in the face of global change. The fact that SOC was positively related to annual precipitation and negatively related to mean annual temperature suggests that future climate change predictions of increased temperature and reduced precipitation may strongly reduce the potential of Spanish soils as C sinks. However, this may be mediated by changes in vegetation cover (e.g. by favouring the development of forests associated to higher SOC values) and exacerbated by perturbations such as fire. The estimations presented here provide a baseline to estimate future changes in soil C stocks and to assess their vulnerability to key global change drivers, and should inform future actions aimed at the conservation and management of C stocks.

Two new species of Acremonium from Spanish soils

In a survey on the diversity of microfungi in Spanish soils, two new species of Acremonium were found. Both species were characterized as having more or less erect, mostly branched conidiophores bearing whorls of acicular phialides. In addition, one of these species, Acremonium asperulatum, produced abundant chlamydospores and globose rough-walled conidia. The other species, Acremonium variecolor, produced a brownish diffusible pigment and smooth-walled, subglobose conidia with apiculate base; sessile conidia inserted directly on vegetative hyphae also were present. The analysis of the sequences of the ITS region, the D1/D2 domains of the 28S rRNA gene and a fragment of the actin gene revealed relationships of both species with members of the Bionectriaceae (Hypocreales). Genetic differences were observed with morphologically similar species.

Map on predicted deposition of Cs-137 in Spanish soils from geostatistical analyses

The knowledge of the distribution of 137Cs deposition over Spanish mainland soils along with the geographical, physical and morphological terrain information enable us to know the 137Cs background content in soil. This could be useful as a tool in a hypothetical situation of an accident involving a radioactive discharge or in soil erosion studies. A Geographic Information System (GIS) would allow the gathering of all the mentioned information. In this work, gamma measurements of 137Cs on 34 Spanish mainland soils, rainfall data taken from 778 weather stations, soil types and geographical and physical terrain information were input into a GIS. Geostatistical techniques were applied to interpolate values of 137Cs activity at unsampled places, obtaining prediction maps of 137Cs deposition. Up to now, geostatistical methods have been used to model spatial continuity of data. Through semivariance and cross-covariance functions the spatial correlation of such data can be studied and described. Ordinary and simple kriging techniques were carried out to map spatial patterns of 137Cs deposition, and ordinary and simple co-kriging were used to improve the prediction map obtained through a second related variable: namely the rainfall. To choose the best prediction map of 137Cs deposition, the spatial dependence of the variable, the correlation coefficient and the prediction errors were evaluated using the different models previously mentioned. The best result for 137Cs deposition map was obtained when applying the co-kriging techniques.

Long-term amendment of Spanish soils with sewage sludge: Effects on soil functioning

Wastewater treatment processes generate highly biodegradable solid wastes. As their fate is an environmental issue of great concern, public administrations promote a sustainable management of urban wastes. The objective of the present study was to analyze the systematic and periodical use, for 16 years, of anaerobically digested sewage sludge as an agricultural fertilizer by assessing the effects on some soil physical–chemical, functional, and ecotoxicological properties. The results showed that the input of sludge enhances soil properties proportionally to the application doses and/or frequency. The organic amendments increased the organic matter content (and its aromaticity), the soil nitrogen, and the microbial activity, improving carbon and nitrogen mineralization processes and some enzymatic functions. However, a maximum dose was identified (40Mgha−1year−1), beyond which soil properties do not improve, and may even worsen. Regarding environmental risks, although the bioluminescent bacteria test showed no toxicity on soil extracts, potential adverse effects such as some potentially toxic elements accumulation, phytotoxicity and the likelihood of groundwater pollution by nitrates or dissolved organic matter should be taken into account. The complementarity of studying soil functioning parameters and ecotoxicological effects, together with the analysis of pollutant content, must be enhanced. This assures a more realistic assessment of long-term effects of sewage sludge-amended soils.

Distribution and efficiency of Rhizobium leguminosarum strains nodulating Phaseolus vulgaris in Northern Spanish soils: Selection of native strains that replace conventional N fertilization

Phaseolus vulgaris is a legume extensively cultivated in Spain, León province being the most important producer. This province produces selected varieties of common bean highly appreciated by their quality that warrants a Protected Geographic Indication (PGI). In this work we analysed the rhizobia present in nodules of the variety “Riñón” in several soils from León province in order to select native rhizobial strains to be used as biofertilizers. The analysis of rrs and housekeeping genes of these strains showed that they belong to two phylogenetic groups within Rhizobium leguminosarum (I and II). Although the group II strains were most abundant in nodules, very effective strains were also found in group I. Strains LCS0306 from group I and LBM1123 from group II were the best nitrogen fixers among all strains isolated and were selected for field experiments. The field research showed that the biofertilization of common bean with native and selected rhizobial strains can completely replace the fertilization with chemical N fertilizers. The biofertiliser designed in such way, was valid for the whole agroecological area, regardless the specific properties of each soil and microclimatic conditions. This conclusion can be generalised as a strategy for the development of biofertilisers in different agroecological conditions worldwide.

Selenium and heavy metals content in some Mediterranean soils

The study of metal contents in industrial, agricultural or/and polluted soils compared with natural or unpolluted soils is currently necessary to obtain reference values and to assess soil contamination. Nonetheless, very few works published appear in international journals on elements like Se, Li and Sr in Spanish soils. This study determines the total levels of Se, Li, Sr, As, Cd, Co, Cr, Cu, Ni, Pb, V, Zn, Fe, Mn and Ba in 14 natural (unpolluted) soils (Gypsisols, Leptosols, Arenosols and Acrisols), 14 agricultural soils (Anthrosols, Fluvisols and Luvisols), and 4 industrial–urban affected-surface soil horizons (Anthrosols and Fluvisols) of Eastern Spain. The geochemical baseline concentrations (GBC) and reference values (RV) have been established, and the relationships among elements and also between soil properties and elemental concentrations have been analysed. The RV obtained in this study were (mg kg −1): Se 2.68, Li 115, Sr 298, Cd 0.97, Co 35, Cr 217, Cu 46, Ni 50, Pb 137, V 120, Zn 246, Fe 124,472, Mn 2691, and Ba 743. The RV for Se and Li were used as a preliminary approach to assess soil contamination in Spanish soils. The results confirm human impact on Sr, As, Cd, Cr, Cu, Ni, Pb and Zn soil concentrations, but evidence no deviation from natural Se, Li, Co, V, Fe, Mn and Ba concentrations. The results obtained from the statistical analysis reveal significant correlations between some elements and clay and soil organic matter (SOM) contents, indicating that metal concentrations are controlled by soil composition. One particularly interesting finding is the high correlation coefficients obtained between SOM and Se, Cd, Cr, V, Fe, and Mn, and between clay and Cd, Zn, V, Fe and Mn. Once again, these facts confirm the role of SOM and clay minerals in soil functions and that soil is an ecosystem element responsible for maintaining environmental quality.