Age-dependent effects of carotenoids on sexual ornaments and reproductive performance of a long-lived seabird

Understanding the mechanisms mediating age-dependent reproductive effort in wild animals is a major challenge in evolutionary ecology. The acquisition of molecules able to deter somatic damage could potentially improve reproductive performance, especially in physiologically compromised individuals. Carotenoids are antioxidants and immunoenhancers that animals can only obtain through diet. We tested in male blue-footed boobies whether carotenoid supplementation during incubation improves condition and reproductive performance, particularly of older males. Old birds showed lower levels of peroxidative damage than middle-aged males; however, changes of males' antioxidant capacity and peroxidative damage were not affected by carotenoid supplementation irrespective of age. Experimental old males displayed more colourful feet than control old males after 13 days of first carotenoid supplementation. The effect of carotenoids on male reproductive performance was dependent upon age: young males receiving extra-carotenoids produced chicks with smaller growth rates than control young males, whereas old males receiving extra-carotenoids produced chicks with higher growth rates than control old males. Our results highlight that carotenoids influence the relationship between age and reproductive performance and suggest a positive effect of carotenoids on reproduction of old males.     

Weather effects on avian breeding performance and implications of climate change

The influence of recent climate change on the world's biota has manifested broadly, resulting in latitudinal range shifts, advancing dates of arrival of migrants and onset of breeding, and altered community relationships. Climate change elevates conservation concerns worldwide because it will likely exacerbate a broad range of identified threats to animal populations. In the past few decades, grassland birds have declined faster than other North American avifauna, largely due to habitat threats such as the intensification of agriculture. We examine the effects of local climatic variations on the breeding performance of a bird endemic to the shortgrass prairie, the Lark Bunting (Calamospiza melanocorys) and discuss the implications of our findings relative to future climate predictions. Clutch size, nest survival, and productivity all positively covaried with seasonal precipitation; yet relatively intense daily precipitation events temporarily depressed daily survival of nests. Nest survival was positively related to average temperatures during the breeding season. Declining summer precipitation may reduce the likelihood that Lark Buntings can maintain stable breeding populations in eastern Colorado although average temperature increases of up to 3 degrees C (within the range of this study) may ameliorate declines in survival expected with drier conditions. Historic climate variability in the Great Plains selects for a degree of vagility and opportunism rather than strong site fidelity and specific adaptation to local environments. These traits may lead to northerly shifts in distribution if climatic and habitat conditions become less favorable in the drying southern regions of the Great Plains. Distributional shifts in Lark Buntings could be constrained by future changes in land use, agricultural practices, or vegetative communities that result in further loss of shortgrass prairie habitats.     

Extrinsic and intrinsic determinants of winter foraging and breeding phenology in a temperate seabird

In temperate regions, winter presents animals with a number of challenges including depressed food abundance, increased daily energy requirements, higher frequency of extreme weather events and shortened day length. Overcoming these constraints is critical for overwintering survival and scheduling of future breeding of long-lived species and is likely to be state dependent, associated with intrinsic abilities such as food acquisition rates. We examined the relationship between environmental and intrinsic factors on overwintering foraging and subsequent breeding phenology of the European shag Phalacrocorax aristotelis, a diurnal marine predator. We tested a range of hypotheses relating to overwintering foraging time and location. We found that individuals greatly increased their foraging time in winter to a peak of more than 90% of available daylight at the winter solstice. The seasonal patterns of foraging time appear to be driven by a combination of light levels and weather conditions and may be linked to the availability of the shag's principal prey, the lesser sandeel Ammodytes marinus. There was no evidence that shags dispersed south in winter to increase potential foraging time. Foraging time decreased after the winter solstice and, crucially, was correlated with subsequent breeding phenology, such that individuals that spent less time foraging in February bred earlier. The relationship was much stronger in females than males, in line with their more direct control of timing of breeding. Our results demonstrate that pre-breeding intrinsic foraging ability is critical in determining breeding phenology.     

Inter-annual variability in the breeding performance of six tropical seabird species: influence of life-history traits and relationship with oceanographic parameters

We investigated inter-annual variability in the breeding performance of six tropical seabirds with different life histories. We further examined the extent to which presumed differences among years in food availability to seabirds were related to large-scale oceanic processes (El Niño Southern Oscillation and Indian Ocean Dipole) or local features around the colony (Sea Surface Temperature—SST, chlorophyll-a concentration, wind patterns, association with underwater predators). Two food shortages were recorded: (1) a short-term event operating at a local-scale, lacking any relationship with any measured oceanographic parameters and affecting only a few seabird species; (2) a generalized and protracted event, strikingly characterized by low SST, possibly induced by a La Niña and/or negative IOD event, and with detrimental impacts in the whole seabird community. Overall, species with inflexible feeding habits and fast chick growth rates were more affected by food shortage events than those with more plastic life-history traits.     

Interactive effects of male and female age on extra-pair paternity in a socially monogamous seabird

Females sometimes obtain older sires for their offspring through extra-pair interactions, but how female age influences paternity is largely unexplored and interactive effects across the age span of both sexes have not been analyzed. To test whether female choice of sire age varies with female age in the blue-footed booby (Sula nebouxii), we examined associations between ages of both partners and the probability of extra-pair paternity (EPP) in 350 broods of parents up to 22 years old in a single breeding season. Extra-pair paternity enables a female to select an alternative sire for her offspring and could function to avoid or achieve particular combinations of parental ages. A male age?×?female age interaction revealed that in young females (≤4 years), EPP decreased with increasing age of the social partner, whereas in old females (≥8 years), it increased. Moreover, sires of extra-pair (EP) chicks of young females paired to young males were on average 6.33 years older than the females’ social partners. Since female boobies control copulatory access, this pattern could imply that young females choose old sires for their proven genetic quality and that old females avoid very old males because matings with them may risk infertility or genetic defects in offspring. Taking female age into account and observing across the whole age span may be necessary for understanding female age-based mate choice.     

Understanding abundance patterns of a declining seabird: implications for monitoring.

The Kittlitz's Murrelet (Brachyramphus brevirostris) is a rare, non-colonial seabird often associated with tidewater glaciers and a recent candidate for listing under the Endangered Species Act. We estimated abundance of Kittlitz's Murrelets across space and time from at-sea surveys along the coast of Alaska (USA) and then used these data to develop spatial models to describe abundance patterns and identify environmental factors affecting abundance. Over a five-week period in the summer of 2005, we recorded 794 Kittlitz's Murrelets, 16 Marbled Murrelets (B. marmoratus), and 70 unidentified murrelets. The overall population estimate (N, mean +/- SE) during the peak period (3-9 July) was 1317 +/- 294 birds, decreasing to 68 +/- 37 by the last survey period (31 July-6 August). Density of Kittlitz's Murrelets was highest in pelagic waters of Taan Fjord (18.6 +/- 7.8 birds/km2, mean +/- SE) during 10-16 July. Spatial models identified consistent "hotspots" of Kittlitz's Murrelets, including several small areas where high densities of murrelets were found throughout the survey period. Of the explanatory variables that we evaluated, tidal current strength influenced murrelet abundance most consistently, with higher abundance associated with strong tidal currents. Simulations based on the empirically derived estimates of variation demonstrated that spatial variation strongly influenced power to detect trend, although power changed little across the threefold difference in the coefficient of variation on detection probability. We include recommendations for monitoring Kittlitz's Murrelets (or other marine species) when there is a high degree of uncertainty about factors affecting abundance, especially spatial variability.     

Differential effects of temperature variability on the transmission of a marine parasite

Temperature variability is particularly pronounced in intertidal systems. The importance of considering this variability has been increasingly recognised, especially in the context of climate change and disease dynamics. Here, we investigated the effects of temperature variability on the transmission of the intertidal trematode Maritrema novaezealandensis. The experimental treatments were 15 °C (control), 15 + 5 °C daily, 15 + 10 °C every second day, 15 + 15 °C every third day (overall equal thermal loading), and a heat wave treatment (15 + 10 °C daily). Daily 6 h incubations were carried out corresponding to daytime low tides over a 12-day period. Effects on output of transmission stages (cercariae) from infected Zeacumantus subcarinatus snail hosts and transmission success of cercariae to Paracalliope novizealandiae amphipod hosts were quantified, as well as the survival of amphipods. Results showed differential effects on output and transmission success. The number of cercariae emerging was similar for treatments with equal thermal loading, but was substantially increased in the heat wave treatment. Transmission success was highest and comparable for the treatments with regular daily temperature increases (i.e. 15 + 5 °C and heat wave), compared to other treatments. Amphipod survival was not affected by temperature treatment directly, but by the number of parasites infecting an amphipod, as well as amphipod sex. These results demonstrate that cercarial output depends mostly on total thermal loading, whereas successful infection of amphipods is determined by total time above 15 °C. Repeated exposure to ~25 °C, as expected under a heat wave scenario, therefore increases both transmission pressure and success, and hence, the risk of parasite-induced mortality in amphipods.     

Combining demographic and count-based approaches to identify source-sink dynamics of a threatened seabird.

Identifying source-sink dynamics is of fundamental importance for conservation but is often limited by an inability to determine how immigration and emigration influence population processes. We demonstrate two ways to assess the role of immigration on population processes without directly observing individuals dispersing from one population to another and apply these methods to a population of Marbled Murrelets (Brachyramphus marmoratus) in California (USA). In the first method, the rate of immigration (i) is estimated by subtracting local recruitment (recruitment from within the population due to reproduction) estimated with demographic data from total recruitment (f; recruitment from within the population plus recruitment from other populations) estimated using temporal symmetry mark-recapture models developed by R. Pradel. The second method compares population growth rates estimated with temporal symmetry models (lambdaTS) and/or population growth rates estimated from counts of individuals over multiple sampling periods (lambdaC) with growth estimates from a stage-structured projection matrix model (lambdaM). Both lambdaTS and lambdaC incorporate all demographic processes affecting population change (birth, death, immigration, and emigration), whereas matrix models are usually constructed without incorporating immigration. Thus, if lambdaTS and lambdaC are > or = 1 and lambdaM < 1, the population is sustained by immigration and is considered to be a sink. Using the first method, recruitment estimated with temporal symmetry models was high (f= 0.182, SE = 0.058), the mean adult birth rate, as estimated using the ratio of juveniles to > or = 1 year old individuals (observed during ship-based surveys) was low (bA = 0.039, SE = 0.014), and immigration was 0.160 (SE = 0.057). Using the second method, murrelet numbers in central California were stable (lambdaC = 1.058, SE = 0.047; lambdaTS = 1.064, SE = 0.033), but were projected to decline 9.5% annually in the absence of immigration (lambdaM = 0.905, SE = 0.053). Our results suggest that Marbled Murrelets in central California represent a sink population that is stable but would decline in the absence of immigration from larger populations to the north. However, the extent to which modeled immigration is due to permanent recruitment or temporarily dispersing individuals that simply mask population declines is uncertain.     

Multiple plant exploiters on a shared host: testing for nonadditive effects on plant performance.

The combined impact of multiple plant parasites on plant performance can either be additive (the total damage equals the sum of the individual effects) or nonadditive (synergistic or antagonistic damage). Two statistical models are available for testing the independent (=additive) effects of two factors. Here we suggest that the natural history of the plant-parasite system should motivate the choice of a statistical model to test for additivity. Using in-field, manipulative experiments, we examined the interactions between the herbivorous mite Calacarus flagelliseta Fletchmann, De Moraes and Barbosa (Acari: Eriophyidae), the fungal pathogen Oidium caricae F. Noack (a powdery mildew), and their host plant Carica papaya L. in Hawaii. First, we found that herbivorous mites had a moderate negative effect on powdery mildew: when mites were absent, powdery mildew colonies were larger and more numerous. Second, we showed that each plant parasite, when evaluated alone, significantly reduced several measures of plant performance. Third, we found that the combined impact of mites and mildew on plant performance is mostly additive and, for a few variables, less than additive. Finally, we explored compensatory responses and found no evidence for nonlinearities in the relationship between plant performance and cumulative parasite impact. Plants are almost universally subject to attack by multiple herbivores and pathogens; thus a deeper understanding of how multiple plant parasites shape each other's population dynamics and plant performance is essential to understanding plant-parasite interactions.     

A physically based statistical model of sand abrasion effects on periphyton biomass

The magnitude and frequency of discharge and fine sediment delivery to rivers can influence riverine food webs through the frequency of scour of algae from the streambed. Models that simulate changes in algal biomass are not very accurate for long periods with frequent low-magnitude flow events. During these periods, sand is mobilized over a stable gravel bed and periphyton losses are patchy at the reach scale. At the patch scale, we examine if an established threshold for rapid sand transport is also a periphyton perturbation threshold. We also develop and validate a statistical rock scale periphyton saltation abrasion model (PSAM) to simulate the abrasive effects of sand, transported by a hopping motion called saltation, on post-flow event biomass. Data were collected from 15 riffles of a Canadian Atlantic salmon river. The threshold clearly divided bed patches with high biomass and low transport rates, from those patches with low biomass and high transport rates. A dimensionally balanced PSAM regression model including explained 57% of the variance in post-flow event biomass. The validated model indicates that periphyton biomass decreases with increasing sand transport rates . Biomass was higher if the microscopic algae were protected from abrasion by growing either above the near-bed layer of saltating sand or within a mat containing more resistant macroalgae (e.g. Nostoc). The use of in our models facilitates testing of our findings in other hydro-sedimentary environments because W^* is a dimensionless scaling parameter that is well established in sediment transport literature. New insight is provided regarding modelling local heterogeneity in post-flow event biomass. These developments are essential to enable more accurate assessments of how periphyton biomass will change with the increase in the recurrence frequency of small flow events (and sand supply) associated with urbanization and climate change.     

Consequences of dominance: a review of evenness effects on local and regional ecosystem processes

The composition of communities is strongly altered by anthropogenic manipulations of biogeochemical cycles, abiotic conditions, and trophic structure in all major ecosystems. Whereas the effects of species loss on ecosystem processes have received broad attention, the consequences of altered species dominance for emergent properties of communities and ecosystems are poorly investigated. Here we propose a framework guiding our understanding of how dominance affects species interactions within communities, processes within ecosystems, and dynamics on regional scales. Dominance (or the complementary term, evenness) reflects the distribution of traits in a community, which in turn affects the strength and sign of both intraspecifc and interspecific interactions. Consequently, dominance also mediates the effect of such interactions on species coexistence. We review the evidence for the fact that dominance directly affects ecosystem functions such as process rates via species identity (the dominant trait) and evenness (the frequency distribution of traits), and indirectly alters the relationship between process rates and species richness. Dominance also influences the temporal and spatial variability of aggregate community properties and compositional stability (invasibility). Finally, we propose that dominance affects regional species coexistence by altering metacommunity dynamics. Local dominance leads to high beta diversity, and rare species can persist because of source-sink dynamics, but anthropogenically induced environmental changes result in regional dominance and low beta diversity, reducing regional coexistence. Given the rapid anthropogenic alterations of dominance in many ecosystems and the strong implications of these changes, dominance should be considered explicitly in the analysis of consequences of altered biodiversity.     

Direct and interactive effects of enemies and mutualists on plant performance: a meta-analysis

Plants engage in multiple, simultaneous interactions with other species; some (enemies) reduce and others (mutualists) enhance plant performance. Moreover, effects of different species may not be independent of one another; for example, enemies may compete, reducing their negative impact on a plant. The magnitudes of positive and negative effects, as well as the frequency of interactive effects and whether they tend to enhance or depress plant performance, have never been comprehensively assessed across the many published studies on plant-enemy and plant-mutualist interactions. We performed a meta-analysis of experiments in which two enemies, two mutualists, or an enemy and a mutualist were manipulated factorially. Specifically, we performed a factorial meta-analysis using the log response ratio. We found that the magnitude of (negative) enemy effects was greater than that of (positive) mutualist effects in isolation, but in the presence of other species, the two effects were of comparable magnitude. Hence studies evaluating single-species effects of mutualists may underestimate the true effects found in natural settings, where multiple interactions are the norm and indirect effects are possible. Enemies did not on average influence the effects on plant performance of other enemies, nor did mutualists influence the effects of mutualists. However, these averages mask significant and large, but positive or negative, interactions in individual studies. In contrast, mutualists ameliorated the negative effects of enemies in a manner that benefited plants; this overall effect was driven by interactions between pathogens and belowground mutualists (bacteria and mycorrhizal fungi). The high frequency of significant interactive effects suggests a widespread potential for diffuse rather than pairwise coevolutionary interactions between plants and their enemies and mutualists. Pollinators and mycorrhizal fungi enhanced plant performance more than did bacterial mutualists. In the greenhouse (but not the field), pathogens reduced plant performance more than did herbivores, pathogens were more damaging to herbaceous than to woody plants, and herbivores were more damaging to crop than to non-crop plants (suggesting evolutionary change in plants or herbivores following crop domestication). We discuss how observed differences in effect size might be confounded with methodological differences among studies.     

Influence of the local spring warming on the breeding phenology in blackcap (Sylvia atricapilla) in Croatia

Recent papers have shown that climate change affects many species, including birds. Several papers from different areas have devoted attention to the negative correlation between the dates of breeding and year, and also negative correlation between the dates of breeding and air spring temperatures. Recent work give some firm evidence for postulating causal relationship between spring temperature and laying dates. We used 31 years (1979-2009) of data from blackcap Sylvia atricapilla in Mokrice area, Northwestern Croatia to assess whether there has been any systematic change in breeding phenology through time.Among environmental factors possibly affecting the breeding date, consideration was given to mean monthly air temperatures (April-May). Over the past three decades, the breeding date of blackcap in our study population has changed significantly: They started breeding progressively earlier (11.77 days). Two regression analysis with laying date as criterion variable showed that spring temperatures can significantly predict variation of laying date and that there are also some other unknown factors which significantly explain variation of laying date. We conclude that blackcaps across Northwestern Croatia are breeding earlier and that mean air spring temperatures is probably the most important factor causing it, among other factors.     

Ecological effects of cultivation on the machair sand dune systems of the Outer Hebrides, Scotland

The machair sand dune systems of the Outer Hebrides of Scotland are a unique habitat, which is rare within both a global and European context. Unusually, the machair habitat also represents an agricultural resource that is very important to the Hebridean people, having been subject to both grazing and cultivation throughout the historical period. Following designation as an Environmentally Sensitive Area (ESA) in 1988, the machairs of South Uist have been studied with the aim of understanding the links between agricultural practice and their plant community and ecosystem dynamics. This research focused primarily on the effects of cultivation practices and their role in maintaining plant species richness and community and ecosystem stability. Within two carefully selected areas, the successional plant communities of machair at different stages of recovery following ploughing and cultivation of cereals and potato patches or ‘lazy beds’ were identified at both a macro- and micro-level. Investigations of the vegetation recovery processes on turves taken from newly ploughed land indicated that the initial stages of recolonization are characterized primarily by rapid vegetative reproduction and growth, although re-vegetation by seeds is also an important factor. The implications of these findings for the long-term management of machair plant communities are discussed and in particular the need to maintain old cultivation practices such as shallow ploughing. The need for more detailed research into both seed banks and seed rain and into processes of vegetative reproduction is stressed. Nomenclature: Clapham et al. (1981) and Stace (1991, 1997) for vascular plants; Hubbard (1984) for grasses; Pankhurst & Mullin (1991) for the regional flora; Dobson (1992) for lichens; Watson (1981) for mosses and liverworts.     

Arkadolith土壤改良剂对内蒙沙区常见牧草生长状况的影响

通过盆栽试验,对土壤采用添加25%土壤改良剂、以及添加12.5%改良剂和12.5%腐殖质的2种处理方式,研究Arkadolith土壤改良剂对冰草（Agropyron desertorum）、沙打旺（Astragalus huangheensis）和准格尔苜蓿（Medicago sativa cv ＂Zhonggeer＂）3种沙生牧草的生长特性的影响。结果表明,3种牧草中,植株高度、生长力、发芽力以及产量都较土壤改良前有所增强,沙打旺和准格尔苜蓿的生长状况要明显好于冰草。含有多种矿物成分的Arkadolith土壤改良剂不仅可以改良贫瘠土壤,提高土壤肥力,促进牧草的增产,而且还可以减少化学肥料的施用,降低成本,做到绿色环保,因此其具有广泛的推广意义。     

Ecological effects of cultivation on the machair sand dune systems of the Outer Hebrides, Scotland

The machair sand dune systems of the Outer Hebrides of Scotland are a unique habitat, which is rare within both a global and European context. Unusually, the machair habitat also represents an agricultural resource that is very important to the Hebridean people, having been subject to both grazing and cultivation throughout the historical period. Following designation as an Environmentally Sensitive Area (ESA) in 1988, the machairs of South Uist have been studied with the aim of understanding the links between agricultural practice and their plant community and ecosystem dynamics. This research focused primarily on the effects of cultivation practices and their role in maintaining plant species richness and community and ecosystem stability. Within two carefully selected areas, the successional plant communities of machair at different stages of recovery following ploughing and cultivation of cereals and potato patches or ‘lazy beds’ were identified at both a macro- and micro-level. Investigations of the vegetation recovery processes on turves taken from newly ploughed land indicated that the initial stages of recolonization are characterized primarily by rapid vegetative reproduction and growth, although re-vegetation by seeds is also an important factor. The implications of these findings for the long-term management of machair plant communities are discussed and in particular the need to maintain old cultivation practices such as shallow ploughing. The need for more detailed research into both seed banks and seed rain and into processes of vegetative reproduction is stressed.     

Assessing the effects of pseudo-absences on predictive distribution model performance

Modelling species distributions with presence data from atlases, museum collections and databases is challenging. In this paper, we compare seven procedures to generate pseudo-absence data, which in turn are used to generate GLM-logistic regressed models when reliable absence data are not available. We use pseudo-absences selected randomly or by means of presence-only methods (ENFA and MDE) to model the distribution of a threatened endemic Iberian moth species (Graellsia isabelae). The results show that the pseudo-absence selection method greatly influences the percentage of explained variability, the scores of the accuracy measures and, most importantly, the degree of constraint in the distribution estimated. As we extract pseudo-absences from environmental regions further from the optimum established by presence data, the models generated obtain better accuracy scores, and over-prediction increases. When variables other than environmental ones influence the distribution of the species (i.e., non-equilibrium state) and precise information on absences is non-existent, the random selection of pseudo-absences or their selection from environmental localities similar to those of species presence data generates the most constrained predictive distribution maps, because pseudo-absences can be located within environmentally suitable areas. This study shows that if we do not have reliable absence data, the method of pseudo-absence selection strongly conditions the obtained model, generating different model predictions in the gradient between potential and realized distributions.     

Spatial structure effects on the detection of patches boundaries using local operators

Landscapes exhibit various degrees of spatial heterogeneity according to the differential intensity and interactions among processes and disturbances that they are subjected to. The management of these spatially dynamical landscapes requires that we can accurately map them and monitor the evolution of their spatial arrangement through time. Such a mapping requires first the delineation of various spatial features present in the landscape such as patches and their boundaries. However, there are several environmental (spatial variability) as well as technical (spatial resolution) factors that impair our ability to accurately delineate patches and their boundaries as polygons. Here, we investigate how the spatial structure and spatial resolution of the data affect the accuracy of detecting patches and their boundaries over simulated landscapes and real data. Simulated landscapes consisted of two patches with parameterized spatial properties (patches’ level of spatial autocorrelation, mean value and variance) separated by a boundary of known location. Real data allowed the investigation of a more complex landscape where there is a known transition between two forest domains with unknown spatial properties. Boundary locations are defined using the lattice-wombling edge detector at various aggregation levels and the degree of patch homogeneity is determined using Getis-Ord’s G*. Results show that boundary detection using a local edge detector is greatly affected by the spatial conditions of the data, namely variance, abruptness of the spatial gradient between two patches and patches’ level of spatial autocorrelation. They also suggest that data aggregation is not a panacea for bringing out the ecological process creating the patches and that indicators derived from local measures of spatial association can be complementary tools for analysing spatial structures affecting boundary delineation.

Natural resource depletion and the economic performance of local districts: suggestions from a within-country analysis

SUMMARY

Studying the interactions among ecological factors and economic performance measured at the local scale is necessary to suggest policies able to mitigate natural resource depletion in complex ecosystems, like those in the Mediterranean region. The aim of this paper is to explore the relationship among natural resource depletion due to land degradation (LD) and some economic characteristics of local labour systems in Italy, a country where many areas, especially southern regions, are vulnerable. LD was estimated using an environmental sensitive area index, covering the national territory at a fine scale, which takes into account ecological factors such as climate, soil and vegetation. Economic indicators used here consider per capita income, land productivity, production and labour productivity by sector (agriculture, industry and services). A regression analysis was built-up at the LLS scale using an index change over time (1990–2000) as dependent variable and economic variables as predictors. A principal component analysis (PCA) was carried out to synthesise the outputs of the regression analysis. Results indicate a negative relationship among per capita income and LD over the whole study area. However, other variables showed a significant correlation with the dependent variable, highlighting the importance of local factors to increase land vulnerability. In order to clarify the contribution of economic factors to natural resource depletion it is necessary to drive integrate policies to combat LD in dry Mediterranean areas.