Distribution and risk factors for Plasmodium and helminth co-infections: a cross-sectional survey among children in Bagamoyo district, coastal region of Tanzania.

BACKGROUND: Plasmodium and soil transmitted helminth infections (STH) are a major public health problem, particularly among children. There are conflicting findings on potential association between these two parasites. This study investigated the Plasmodium and helminth co-infections among children aged 2 months to 9 years living in Bagamoyo district, coastal region of Tanzania. METHODS: A community-based cross-sectional survey was conducted among 1033 children. Stool, urine and blood samples were examined using a broad set of quality controlled diagnostic methods for common STH (Ascaris lumbricoides, hookworm, Strongyloides stercoralis, Enterobius vermicularis, Trichuris trichura), schistosoma species and Wuchereria bancrofti. Blood slides and malaria rapid diagnostic tests (mRDTs) were utilized for Plasmodium diagnosis. RESULTS: Out of 992 children analyzed, the prevalence of Plasmodium infection was 13% (130/992), helminth 28.5% (283/992); 5% (50/992) had co-infection with Plasmodium and helminth. The prevalence rate of Plasmodium, specific STH and co-infections increased significantly with age (p < 0.001), with older children mostly affected except for S. stercoralis monoinfection and co-infections. Spatial variations of co-infection prevalence were observed between and within villages. There was a trend for STH infections to be associated with Plasmodium infection [OR adjusted for age group 1.4, 95% CI (1.0-2.1)], which was more marked for S. stercoralis (OR = 2.2, 95% CI (1.1-4.3). Age and not schooling were risk factors for Plasmodium and STH co-infection. CONCLUSION: The findings suggest that STH and Plasmodium infections tend to occur in the same children, with increasing prevalence of co-infection with age. This calls for an integrated approach such as using mass chemotherapy with dual effect (e.g., ivermectin) coupled with improved housing, sanitation and hygiene for the control of both parasitic infections.

Impact of model complexity on cross-temporal transferability in Maxent species distribution models: An assessment using paleobotanical data

Maximum entropy modeling (Maxent) is a widely used algorithm for predicting species distributions across space and time. Properly assessing the uncertainty in such predictions is non-trivial and requires validation with independent datasets. Notably, model complexity (number of model parameters) remains a major concern in relation to overfitting and, hence, transferability of Maxent models. An emerging approach is to validate the cross-temporal transferability of model predictions using paleoecological data. In this study, we assess the effect of model complexity on the performance of Maxent projections across time using two European plant species (Alnus giutinosa (L.) Gaertn. and Corylus avellana L) with an extensive late Quaternary fossil record in Spain as a study case. We fit 110 models with different levels of complexity under present time and tested model performance using AUC (area under the receiver operating characteristic curve) and AlCc (corrected Akaike Information Criterion) through the standard procedure of randomly partitioning current occurrence data. We then compared these results to an independent validation by projecting the models to mid-Holocene (6000 years before present) climatic conditions in Spain to assess their ability to predict fossil pollen presence-absence and abundance. We find that calibrating Maxent models with default settings result in the generation of overly complex models. While model performance increased with model complexity when predicting current distributions, it was higher with intermediate complexity when predicting mid-Holocene distributions. Hence, models of intermediate complexity resulted in the best trade-off to predict species distributions across time. Reliable temporal model transferability is especially relevant for forecasting species distributions under future climate change. Consequently, species-specific model tuning should be used to find the best modeling settings to control for complexity, notably with paleoecological data to independently validate model projections. For cross-temporal projections of species distributions for which paleoecological data is not available, models of intermediate complexity should be selected.

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau.

BACKGROUND: The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section. RESULTS: We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns. CONCLUSION: Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

Evolutionary patchwork of an insecticidal toxin shared between plant-associated pseudomonads and the insect pathogens Photorhabdus and Xenorhabdus.

BACKGROUND: Root-colonizing fluorescent pseudomonads are known for their excellent abilities to protect plants against soil-borne fungal pathogens. Some of these bacteria produce an insecticidal toxin (Fit) suggesting that they may exploit insect hosts as a secondary niche. However, the ecological relevance of insect toxicity and the mechanisms driving the evolution of toxin production remain puzzling. RESULTS: Screening a large collection of plant-associated pseudomonads for insecticidal activity and presence of the Fit toxin revealed that Fit is highly indicative of insecticidal activity and predicts that Pseudomonas protegens and P. chlororaphis are exclusive Fit producers. A comparative evolutionary analysis of Fit toxin-producing Pseudomonas including the insect-pathogenic bacteria Photorhabdus and Xenorhadus, which produce the Fit related Mcf toxin, showed that fit genes are part of a dynamic genomic region with substantial presence/absence polymorphism and local variation in GC base composition. The patchy distribution and phylogenetic incongruence of fit genes indicate that the Fit cluster evolved via horizontal transfer, followed by functional integration of vertically transmitted genes, generating a unique Pseudomonas-specific insect toxin cluster. CONCLUSIONS: Our findings suggest that multiple independent evolutionary events led to formation of at least three versions of the Mcf/Fit toxin highlighting the dynamic nature of insect toxin evolution.

Ligand binding to heparan sulfate proteoglycans induces their aggregation and distribution along actin cytoskeleton

Cell surface heparan sulfate proteoglycans (HSPGs) participate in molecular events that regulate cell adhesion, migration, and proliferation. The present study demonstrates that soluble heparin-binding proteins or cross-linking antibodies induce the aggregation of cell surface HSPGs and their distribution along underlying actin filaments. Immunofluorescence and confocal microscopy and immunogold and electron microscopy indicate that, in the absence of ligands, HSPGs are irregularly distributed on the fibroblast cell surface, without any apparent codistribution with the actin cytoskeleton. In the presence of ligand (lipoprotein lipase) or antibodies against heparan sulfate, HSPGs aggregate and colocalize with the actin cytoskeleton. Triton X-100 extraction and immunoelectron microscopy have demonstrated that in this condition HSPGs were clustered and associated with the actin filaments. Crosslinking experiments that use biotinylated lipoprotein lipase have revealed three major proteoglycans as binding sites at the fibroblast cell surface. These cross-linked proteoglycans appeared in the Triton X-100 insoluble fraction. Platinum/carbon replicas of the fibroblast surface incubated either with lipoprotein lipase or antiheparan sulfate showed large aggregates of HSPGs regularly distributed along cytoplasmic fibers. Quantification of the spacing between HSPGs by confocal microscopy confirmed that the nonrandom distribution of HSPG aggregates along the actin cytoskeleton was induced by ligand binding. When cells were incubated either with lipoprotein lipase or antibodies against heparan sulfate, the distance between immunofluorescence spots was uniform. In contrast, the spacing between HSPGs on fixed cells not incubated with ligand was more variable. This highly organized spatial relationship between actin and proteoglycans suggests that cortical actin filaments could organize the molecular machinery involved in signal transduction and molecular movements on the cell surface that are triggered by heparin-binding proteins.

Soft tissue artifact distribution on lower limbs during treadmill gait: Influence of skin markers' location on cluster design.

Segment poses and joint kinematics estimated from skin markers are highly affected by soft tissue artifact (STA) and its rigid motion component (STARM). While four marker-clusters could decrease the STA non-rigid motion during gait activity, other data, such as marker location or STARM patterns, would be crucial to compensate for STA in clinical gait analysis. The present study proposed 1) to devise a comprehensive average map illustrating the spatial distribution of STA for the lower limb during treadmill gait and 2) to analyze STARM from four marker-clusters assigned to areas extracted from spatial distribution. All experiments were realized using a stereophotogrammetric system to track the skin markers and a bi-plane fluoroscopic system to track the knee prosthesis. Computation of the spatial distribution of STA was realized on 19 subjects using 80 markers apposed on the lower limb. Three different areas were extracted from the distribution map of the thigh. The marker displacement reached a maximum of 24.9mm and 15.3mm in the proximal areas of thigh and shank, respectively. STARM was larger on thigh than the shank with RMS error in cluster orientations between 1.2° and 8.1°. The translation RMS errors were also large (3.0mm to 16.2mm). No marker-cluster correctly compensated for STARM. However, the coefficient of multiple correlations exhibited excellent scores between skin and bone kinematics, as well as for STARM between subjects. These correlations highlight dependencies between STARM and the kinematic components. This study provides new insights for modeling STARM for gait activity.

Beyond contact-based transmission networks: the role of spatial coincidence.

Animal societies rely on interactions between group members to effectively communicate and coordinate their actions. To date, the transmission properties of interaction networks formed by direct physical contacts have been extensively studied for many animal societies and in all cases found to inhibit spreading. Such direct interactions do not, however, represent the only viable pathways. When spreading agents can persist in the environment, indirect transmission via 'same-place, different-time' spatial coincidences becomes possible. Previous studies have neglected these indirect pathways and their role in transmission. Here, we use rock ant colonies, a model social species whose flat nest geometry, coupled with individually tagged workers, allowed us to build temporally and spatially explicit interaction networks in which edges represent either direct physical contacts or indirect spatial coincidences. We show how the addition of indirect pathways allows the network to enhance or inhibit the spreading of different types of agent. This dual-functionality arises from an interplay between the interaction-strength distribution generated by the ants' movement and environmental decay characteristics of the spreading agent. These findings offer a general mechanism for understanding how interaction patterns might be tuned in animal societies to control the simultaneous transmission of harmful and beneficial agents.

Spatial pattern of landslides in Swiss Rhone valley

The present study analyses the spatial pattern of quaternary gravitational slope deformations (GSD) and historical/present-day instabilities (HPI) inventoried in the Swiss Rhone Valley. The main objective is to test if these events are clustered (spatial attraction) or randomly distributed (spatial independency). Moreover, analogies with the cluster behaviour of earthquakes inventoried in the same area were examined. The Ripley's K-function was applied to measure and test for randomness. This indicator allows describing the spatial pattern of a point process at increasing distance values. To account for the non-constant intensity of the geological phenomena, a modification of the K-function for inhomogeneous point processes was adopted. The specific goal is to explore the spatial attraction (i.e. cluster behaviour) among landslide events and between gravitational slope deformations and earthquakes. To discover if the two classes of instabilities (GSD and HPI) are spatially independently distributed, the cross K-function was computed. The results show that all the geological events under study are spatially clustered at a well-defined distance range. GSD and HPI show a similar pattern distribution with clusters in the range 0.75?9 km. The cross K-function reveals an attraction between the two classes of instabilities in the range 0?4 km confirming that HPI are more prone to occur within large-scale slope deformations. The K-function computed for GSD and earthquakes indicates that both present a cluster tendency in the range 0?10 km, suggesting that earthquakes could represent a potential predisposing factor which could influence the GSD distribution.

Disjunct populations of European vascular plant species keep the same climatic niches

Aim Previous research on how climatic niches vary across species ranges has focused on a limited number of species, mostly invasive, and has not, to date, been very conclusive. Here we assess the degree of niche conservatism between distant populations of native alpine plant species that have been separated for thousands of years. Location European Alps and Fennoscandia. Methods Of the studied pool of 888 terrestrial vascular plant species occurring in both the Alps and Fennoscandia, we used two complementary approaches to test and quantify climatic-niche shifts for 31 species having strictly disjunct populations and 358 species having either a contiguous or a patchy distribution with distant populations. First, we used species distribution modelling to test for a region effect on each species' climatic niche. Second, we quantified niche overlap and shifts in niche width (i.e. ecological amplitude) and position (i.e. ecological optimum) within a bi-dimensional climatic space. Results Only one species (3%) of the 31 species with strictly disjunct populations and 58 species (16%) of the 358 species with distant populations showed a region effect on their climatic niche. Niche overlap was higher for species with strictly disjunct populations than for species with distant populations and highest for arctic-alpine species. Climatic niches were, on average, wider and located towards warmer and wetter conditions in the Alps. Main conclusion Climatic niches seem to be generally conserved between populations that are separated between the Alps and Fennoscandia and have probably been so for 10,000-15,000 years. Therefore, the basic assumption of species distribution models that a species' climatic niche is constant in space and time - at least on time scales 104 years or less - seems to be largely valid for arctic-alpine plants.

Emergent geomorphic-vegetation interactions on a subalpine alluvial fan

Following perturbation, an ecosystem (flora, fauna, soil) should evolve as a function of time at a rate conditioned by external variables (relief, climate, geology). More recently, biogeomorphologists have focused upon the notion of co-development of geomorphic processes with ecosystems over very short through to very long (evolutionary) timescales. Alpine environments have been a particular focus of this co-development. However, work in this field has tended to adopt a simplified view of the relationship between perturbation and succession, including: how the landform and ecosystem itself conditions the impact of a perturbation to create a complex spatial response impact; and how perturbations are not simply ecosystem destroyers but can be a significant source of ecosystem resources. What this means is that at the within landform scale, there may well be a complex and dynamic topographic and sedimentological template that co-develops with soil, flora and fauna. Here, we present and test a conceptual model of this template for a subalpine alluvial fan. We combine detailed floristic inventory with soil inventory, determination of edaphic variables and analysis of historical aerial imagery. Spatial variation in the probability of perturbation of sites on the fan surface was associated with down fan variability in the across-fan distribution of fan ages, fan surface channel characteristics and fan surface sedimentology. Floristic survey confirmed that these edaphic factors distinguished site floristic richness and plant communities up until the point that the soil-vegetation system was sufficiently developed to sustain plant communities regardless of edaphic conditions. Thus, the primary explanatory variable was the estimated age of each site, which could be tied back into perturbation history and its spatial expression due to the geometry of the fan: distinct plant communities were emergent both across fan and down fan, a distribution maintained by the way in which the fan dissipates potentially perturbing events.

Will climate change increase the risk of plant invasions into mountains?

Mountain ecosystems have been less adversely affected by invasions of non-native plants than most other ecosystems, partially because most invasive plants in the lowlands are limited by climate and cannot grow under harsher high-elevation conditions. However, with ongoing climate change, invasive species may rapidly move upwards and threaten mid- then high-elevation mountain ecosystems. We evaluated this threat by predicting current and future potential distributions of 48 invasive plant species distributed in Switzerland (CH) and New South Wales (NSW), two areas where climate interacts differently with the elevation gradient. Using a species distribution modeling approach combining two scales, which builds on high-resolution data (< 250 m) but accounts for the global climatic niche of species, we found that different environmental drivers limit the elevation range of invasive species in the two regions, leading to region-specific species responses to climate change. Whereas the optimal suitability for plant invaders is predicted to markedly shift from the lowland to the montane or subalpine zone in CH, such an upward shift is far less pronounced in NSW where montane and subalpine elevations are currently already suitable. Non-native species able to invade the upper reaches of mountains in a future climate will be cold-tolerant in the Swiss Alps but preferring wet soils in the Australian Alps. Other plant traits were only marginally associated with elevation limits. These results demonstrate that a more systematic consideration of future distributions of invasive species is required in conservation plans of not yet invaded mountainous ecosystems.

Evolutionary Ecology of Mountain Birch in Subarctic Stress Gradients: Interplay of Biotic and Abiotic Factors in Plant-Plant Interactions and Evolutionary Processes

In nature, variation for example in herbivory, wind exposure, moisture and pollution impact often creates variation in physiological stress and plant productivity. This variation is seldom clear-cut, but rather results in clines of decreasing growth and productivity towards the high-stress end. These clines of unidirectionally changing stress are generally known as ‘stress gradients’. Through its effect on plant performance, stress has the capacity to fundamentally alter the ecological relationships between individuals, and through variation in survival and reproduction it also causes evolutionary change, i.e. local adaptations to stress and eventually speciation. In certain conditions local adaptations to environmental stress have been documented in a matter of just a few generations. In plant-plant interactions, intensities of both negative interactions (competition) and positive ones (facilitation) are expected to vary along stress gradients. The stress-gradient hypothesis (SGH) suggests that net facilitation will be strongest in conditions of high biotic and abiotic stress, while a more recent ‘humpback’ model predicts strongest net facilitation at intermediate levels of stress. Plant interactions on stress gradients, however, are affected by a multitude of confounding factors, making studies of facilitation-related theories challenging. Among these factors are plant ontogeny, spatial scale, and local adaptation to stress. The last of these has very rarely been included in facilitation studies, despite the potential co-occurrence of local adaptations and changes in net facilitation in stress gradients. Current theory would predict both competitive effects and facilitative responses to be weakest in populations locally adapted to withstand high abiotic stress. This thesis is based on six experiments, conducted both in greenhouses and in the field in Russia, Norway and Finland, with mountain birch (Betula pubescens subsp. czerepanovii) as the model species. The aims were to study potential local adaptations in multiple stress gradients (both natural and anthropogenic), changes in plant-plant interactions under conditions of varying stress (as predicted by SGH), potential mechanisms behind intraspecific facilitation, and factors confounding plant-plant facilitation, such as spatiotemporal, ontogenetic, and genetic differences. I found rapid evolutionary adaptations (occurring within a time-span of 60 to 70 years) towards heavy-metal resistance around two copper-nickel smelters, a phenomenon that has resulted in a trade-off of decreased performance in pristine conditions. Heavy-metal-adapted individuals had lowered nickel uptake, indicating a possible mechanism behind the detected resistance. Seedlings adapted to heavy-metal toxicity were not co-resistant to others forms of abiotic stress, but showed co-resistance to biotic stress by being consumed to a lesser extent by insect herbivores. Conversely, populations from conditions of high natural stress (wind, drought etc.) showed no local adaptations, despite much longer evolutionary time scales. Due to decreasing emissions, I was unable to test SGH in the pollution gradients. In natural stress gradients, however, plant performance was in accordance with SGH, with the strongest host-seedling facilitation found at the high-stress sites in two different stress gradients. Factors confounding this pattern included (1) plant size / ontogenetic status, with seedling-seedling interactions being competition dominated and host-seedling interactions potentially switching towards competition with seedling growth, and (2) spatial distance, with competition dominating at very short planting distances, and facilitation being strongest at a distance of circa ¼ benefactor height. I found no evidence for changes in facilitation with respect to the evolutionary histories of plant populations. Despite the support for SGH, it may be that the ‘humpback’ model is more relevant when the main stressor is resource-related, while what I studied were the effects of ‘non-resource’ stressors (i.e. heavy-metal pollution and wind). The results have potential practical applications: the utilisation of locally adapted seedlings and plant facilitation may increase the success of future restoration efforts in industrial barrens as well as in other wind-exposed sites. The findings also have implications with regard to the effects of global change in subarctic environments: the documented potential by mountain birch for rapid evolutionary change, together with the general lack of evolutionary ‘dead ends’, due to not (over)specialising to current natural conditions, increase the chances of this crucial forest-forming tree persisting even under the anticipated climate change.

The attenuation of concentrations model: a new method for assessing mercury mobility in sediments

In this work we propose a new approach for the determination of the mobility of mercury in sediments based on spatial distribution of concentrations. We chose the Tainheiros Cove, located in the Todos os Santos Bay, Brazil, as the study area, for it has a history of mercury contamination due to a chloro-alkali plant that was active during 12 years. Twenty-six surface sediment samples were collected from the area and mercury concentrations were measured by cold vapour atomic absorption spectrophotometry. A contour map was constructed from the results, indicating that mercury accumulated in a "hot spot" where concentrations reach more than 1 µg g-1. The model is able to estimate mobility of mercury in the sediments based on the distances between iso-concentration contours that determines an attenuation of concentrations factor. Values of attenuation ranged between 0.0729 (East of the hot spot, indicating higher mobility) to 0.7727 (North of the hot spot, indicating lower mobility).

Active sulphide mine tailings impoundments as sources of contaminated drainage: controlling factors, methods of characterisation and geochemical constraints for mitigation

Low quality mine drainage from tailings facilities persists as one of the most significant global environmental concerns related to sulphide mining. Due to the large variation in geological and environmental conditions at mine sites, universal approaches to the management of mine drainage are not always applicable. Instead, site-specific knowledge of the geochemical behaviour of waste materials is required for the design and closure of the facilities. In this thesis, tailings-derived water contamination and factors causing the pollution were investigated in two coeval active sulphide mine sites in Finland: the Hitura Ni mine and the Luikonlahti Cu-Zn-Co-Ni mine and talc processing plant. A hydrogeochemical study was performed to characterise the tailingsderived water pollution at Hitura. Geochemical changes in the Hitura tailings were evaluated with a detailed mineralogical and geochemical investigation (solid-phase speciation, acid mine drainage potential, pore water chemistry) and using a spatial assessment to identify the mechanisms of water contamination. A similar spatial investigation, applying selective extractions, was carried out in the Luikonlahti tailings area for comparative purposes (Hitura low-sulphide tailings vs. Luikonlahti sulphide-rich tailings). At both sites, hydrogeochemistry of tailings seepage waters was further characterised to examine the net results of the processes observed within the impoundments and to identify constraints for water treatment. At Luikonlahti, annual and seasonal variation in effluent quality was evaluated based on a four-year monitoring period. Observations pertinent to future assessment and mine drainage prevention from existing and future tailings facilities were presented based on the results. A combination of hydrogeochemical approaches provided a means to delineate the tailings-derived neutral mine drainage at Hitura. Tailings effluents with elevated Ni, SO₄ ^2- and Fe content had dispersed to the surrounding aquifer through a levelled-out esker and underneath the seepage collection ditches. In future mines, this could be avoided with additional basal liners in tailings impoundments where the permeability of the underlying Quaternary deposits is inadequate, and with sufficiently deep ditches. Based on the studies, extensive sulphide oxidation with subsequent metal release may already initiate during active tailings disposal. The intensity and onset of oxidation depended on e.g. the Fe sulphide content of the tailings, water saturation level, and time of exposure of fresh sulphide grains. Continuous disposal decreased sulphide weathering in the surface of low-sulphide tailings, but oxidation initiated if they were left uncovered after disposal ceased. In the sulphide-rich tailings, delayed burial of the unsaturated tailings had resulted in thick oxidized layers, despite the continuous operation. Sulphide weathering and contaminant release occurred also in the border zones. Based on the results, the prevention of sulphide oxidation should already be considered in the planning of tailings disposal, taking into account the border zones. Moreover, even lowsulphide tailings should be covered without delay after active disposal ceases. The quality of tailings effluents showed wide variation within a single impoundment and between the two different types of tailings facilities assessed. The affecting factors included source materials, the intensity of weathering of tailings and embankment materials along the seepage flow path, inputs from the process waters, the water retention time in tailings, and climatic seasonality. In addition, modifications to the tailings impoundment may markedly change the effluent quality. The wide variation in the tailings effluent quality poses challenges for treatment design. The final decision on water management requires quantification of the spatial and seasonal fluctuation at the site, taking into account changes resulting from the eventual closure of the impoundment. Overall, comprehensive hydrogeochemical mapping was deemed essential in the identification of critical contaminants and their sources at mine sites. Mineralogical analysis, selective extractions, and pore water analysis were a good combination of methods for studying the weathering of tailings and in evaluating metal mobility from the facilities. Selective extractions with visual observations and pH measurements of tailings solids were, nevertheless, adequate in describing the spatial distribution of sulphide oxidation in tailings impoundments. Seepage water chemistry provided additional data on geochemical processes in tailings and was necessary for defining constraints for water treatment.

Distribution Patterns of Soil Properties in a Rural Mediterranean Area in Northeastern Spain

Soil properties on the Cap de Creus Peninsula, NE Spain depend primarily on scarce agricultural practices and early abandonment. In the study area, 90% of which is mainly covered by Cistus shrubs, 8 environments representing variations in land use/land cover and soil properties at different depths were identified. In each environment variously vegetated areas were selected and sampled. The soils, collected at different depths, were classified as Lithic Xerorthents according to the United States Department of Agriculture system of soil classification (USDA-NRCS 1975). Differences in soil properties were largely found according to the evolution of the plant canopy and the land use history. To identify underlying patterns in soil properties related to environmental evolution, factor analysis was performed and factor scores were used to determine how the factor patterns varied between soil variables, soil depths and selected environments. The three-factor model always accounted for 80% of the total variation in the data at the different soil depths. Organic matter was the more relevant soil property at 0–2 cm depth, whereas active minerals (silt and clay) were found to be the most relevant soil parameters controlling soil dynamics at the other depths investigated. Results showed that vineyards and olive tree soils are poorly developed and present worse conditions for mineral and organic compounds. Analysis of factor scores allowed independent assessment of soils, depth and plant cover and demonstrated that soils present the best physico-chemical characteristics under Erica arborea and meadows. In contrast, soils under Cistus monspeliensis were less nutrient rich and less well structured