Stable carbon isotope composition of benthic foraminifera from sediments of the Skagerrak, North Sea

The sediment cores 225514 and 225510 were recovered from 420 and 285 m water depth, respectively. They were investigated for their benthic foraminiferal delta13C during the last 500 years. Both cores were recovered from the southern flank of the Skagerrak. The delta13C values of Uvigerina mediterranea and other shallow infaunal species in both cores indicate that organic matter rain rates to the seafloor varied around a mean value until approximately AD 1950 after which they increased. This increase might result from changes in the North Atlantic Current System and a co-occurring persistently high North Atlantic Oscillation index state in the 1980s to 1990s, rather than from anthropogenic eutrophication. Using delta13C mean values of multiple species, we reconstruct delta13C gradients of dissolved inorganic carbon (DIC) within pore waters for the time periods AD 1500 to 1950 and AD 1950 to 2000. The calculated delta13CDIC ranges, interpreted as indicating total organic matter remineralization due to respiration, are generally bigger in Core 225514 than in Core 225510. Since mean delta13C values of U. mediterranea suggest that organic matter rain rates were similar at both locations, differences in total organic matter remineralization are attributed to differing oxygen availability. However, oxygen concentrations in the overlying bottom water masses are not likely to have differed significantly. Thus, we suggest that organic matter remineralization was controlled by oxygen availability within the sediments, reflecting strong differences in sedimentation rates at the two investigated core sites. Based on the assumptions that tests of benthic foraminiferal species inhabiting the same microhabitat depth should show equal delta13C values unless they are affected by vital effects and that Globobulimina turgida records pore water delta13CDIC, we estimate microhabitat-corrected vital effects for several species with respect to G. turgida: >0.7 per mil for Cassidulina laevigata, >1.3 per mil for Hyalinea balthica, and >0.7 per mil for Melonis barleeanus. Melonis zaandami seems to closely record pore water delta13CDIC.

Radiocarbon ages and 230Th data of five sediment cores from the Björn and Gardar Drift, North Atlantic

In locations of rapid sediment accumulation receiving substantial amounts of laterally transported material the timescales of transport and accurate quantification of the transported material are at the focus of intense research. Here we present radiocarbon data obtained on co-occurring planktic foraminifera, marine haptophyte biomarkers (alkenones) and total organic carbon (TOC) coupled with excess Thorium-230 (230Thxs) measurements on four sediment cores retrieved in 1649-2879 m water depth from two such high accumulation drift deposits in the Northeast Atlantic, Björn and Gardar Drifts. While 230Thxs inventories imply strong sediment focussing, no age offsets are observed between planktic foraminifera and alkenones, suggesting that redistribution of sediments is rapid and occurs soon after formation of marine organic matter, or that transported material contains negligible amounts of alkenones. An isotopic mass balance calculation based on radiocarbon concentrations of co-occurring sediment components leads us to estimate that transported sediment components contain up to 12% of fossil organic matter that is free of or very poor in alkenones, but nevertheless appears to consist of a mixture of fresh and eroded fossil material. Considering all available constraints to characterize transported material, our results show that although focussing factors calculated from bulk sediment 230Thxs inventories may allow useful approximations of bulk redeposition, they do not provide a unique estimate of the amount of each laterally transported sediment component. Furthermore, our findings provide evidence that the occurrence of lateral sediment redistribution alone does not always hinder the use of multiple proxies but that individual sediment fractions are affected to variable extents by sediment focussing.

Investigations on fossil bivalves and wood from Seymour Island, Antarctica

Quasi-periodic variation in sea-surface temperature, precipitation, and sea-level pressure in the equatorial Pacific known as the El Niño - Southern Oscillation (ENSO) is an important mode of interannual variability in global climate. A collapse of the tropical Pacific onto a state resembling a so-called 'permanent El Niño', with a preferentially warmed eastern equatorial Pacific, flatter thermocline, and reduced interannual variability, in a warmer world is predicted by prevailing ENSO theory. If correct, future warming will be accompanied by a shift toward persistent conditions resembling El Niño years today, with major implications for global hydrological cycles and consequent impacts on socioeconomic and ecological systems. However, much uncertainty remains about how interannual variability will be affected. Here, we present multi-annual records of climate derived from growth increment widths in fossil bivalves and co-occurring driftwood from the Antarctic peninsula that demonstrate significant variability in the quasi-biennial and 3-6 year bands consistent with ENSO, despite early Eocene (~50 Mya) greenhouse conditions with global average temperature -10 degrees higher than today. A coupled climate model suggests an ENSO signal and teleconnections to this region during the Eocene, much like today. The presence of ENSO variation during this markedly warmer interval argues for the persistence of robust interannual variability in our future greenhouse world.

Mineralogy, lithology and stable isotope composition of diagenetic siliceous rocks from ODP Leg 129 sites

Ocean Drilling Program Leg 129 recovered chert, porcellanite, and radiolarite from Middle Jurassic to lower Miocene strata from the western Pacific that formed by different processes and within distinct host rocks. These cherts and porcellanites formed by (1) replacement of chalk or limestone, (2) silicification and in-situ silica phase-transformation of bedded clay-bearing biosiliceous deposits, (3) high-temperature silicification adjacent to volcanic flows or sills, and (4) silica phase-transformation of mixed biosiliceous-volcaniclastic sediments. Petrologic and O-isotopic studies highlight the key importance of permeability and time in controlling the formation of dense cherts and porcellanites. The formation of dense, vitreous cherts apparently requires the local addition and concentration of silica. The influence of permeability is shown by two examples, in which: (1) fragments of originally identical radiolarite that were differentially isolated from pore-water circulation by cement-filled fractures were silicified to different degrees, and (2) by the development of secondary porosity during the opal-CT to quartz inversion within conditions of negligible permeability. The importance of time is shown by the presence of quartz chert below, but not above, a Paleogene hiatus at Site 802, indicating that between 30 and 52 m.y. was required for the formation of quartz chert within calcareous-siliceous sediments. The oxygen-isotopic composition for all Leg 129 carbonate- and Fe/Mn-oxide-free whole-rock samples of chert and porcellanite range widely from d18O = 27.8 per mil to 39.8 per mil vs. V-SMOW. Opal-CT samples are consistently richer in 18O (34.1 per mil to 39.3 per mil) than quartz subsamples (27.8 per mil to 35.7 per mil). Using the O-isotopic fractionation expression for quartz-water of Knauth and Epstein (1976) and assuming d18Opore water = -1.0 per mil, model temperatures of formation are 7°-26°C for carbonate-replacement quartz cherts, 22°-25°C for bedded quartz cherts, and 32°-34°C for thermal quartz cherts. Large variations in O-isotopic composition exist at the same burial depth between co-existing silica phases in the same sample and within the same phase in adjacent lithologies. For example, quartz has a wide range of isotopic compositions within a single breccia sample; d18O = 33.4 per mil and 28.0 per mil for early and late stages of fracture-filling cementation, and 31.6 per mil and 30.2 per mil for microcrystalline quartz precipitation within enclosed chert and radiolarite fragments. Similarly, opal-CT d101 spacing varies across lithologic or diagenetic boundaries within single samples. Co-occurring opal-CT and chalcedonic quartz in shallowly buried chert and porcellanite from Sites 800 and 801 have an 8.7 per mil difference in d18O, suggesting that pore waters in the Pigafetta Basin underwent a Tertiary shift to strongly 18O-depleted values due to alteration of underlying Aptian to Albian-Cenomanian volcaniclastic deposits after opal-CT precipitation, but prior to precipitation of microfossil-filling chalcedony.

Station characteristics and copepod fatty alcohol, fatty acid and lipid composition during MSM02/4 in Fram Strait

The biodiversity of pelagic deep-sea ecosystems has received growing scientific interest in the last decade, especially in the framework of international marine biodiversity initiatives, such as Census of Marine Life (CoML). While a growing number of deep-sea zooplankton species has been identified and genetically characterized, little information is available on the mechanisms minimizing inter-specific competition and thus allowing closely related species to co-occur in the deep-sea pelagic realm. Focussing on the two dominant calanoid copepod families Euchaetidae and Aetideidae in Fram Strait, Arctic Ocean, the present study strives to characterize ecological niches of co-occurring species, with regard to vertical distribution, dietary composition as derived from lipid biomarkers, and trophic level on the basis of stable isotope signatures. Closely related species were usually restricted to different depth layers, resulting in a multi-layered vertical distribution pattern. Thus, vertical partitioning was an important mechanism to avoid inter-specific competition. Species occurring in the same depth strata usually belonged to different genera. They differed in fatty acid composition and trophic level, indicating different food preferences. Herbivorous Calanus represent major prey items for many omnivorous and carnivorous species throughout the water column. The seasonal and ontogenetic vertical migration of Calanus acts as a short-cut in food supply for pelagic deep-sea ecosystems in the Arctic.

(Table 2) Radiocarbon data (fMC; fraction modern carbon and conventional 14C ages) of foraminifera, TOC and alkenones of cores ME0005A-24JC, Y69-71P and MC16

Paired radiocarbon measurements on haptophyte biomarkers (alkenones) and on co-occurring tests of planktic foraminifera (Neogloboquadrina dutertrei and Globogerinoides sacculifer) from late glacial to Holocene sediments at core locations ME0005-24JC, Y69-71P, and MC16 from the south-western and central Panama Basin indicate no significant addition of pre-aged alkenones by lateral advection. The strong temporal correspondence between alkenones, foraminifera and total organic carbon (TOC) also implies negligible contributions of aged terrigenous material. Considering controversial evidence for sediment redistribution in previous studies of these sites, our data imply that the laterally supplied material cannot stem from remobilization of substantially aged sediments. Transport, if any, requires syn-depositional nepheloid layer transport and redistribution of low-density or fine-grained components within decades of particle formation. Such rapid and local transport minimizes the potential for temporal decoupling of proxies residing in different grain-size fractions and thus facilitates comparison of various proxies for paleoceanographic reconstructions in this study area. Anomalously old foraminiferal tests from a glacial depth interval of core Y69-71P may result from episodic spillover of fast bottom currents across the Carnegie Ridge transporting foraminiferal sands towards the north.

Impacts of groundwater discharge at Myora Springs (North Stradbroke Island, Australia) on the phenolic metabolism of eelgrass, Zostera muelleri, and grazing by the juvenile rabbitfish, Siganus fuscescens

Myora Springs is one of many groundwater discharge sites on North Stradbroke Island (Queensland, Australia). Here spring waters emerge from wetland forests to join Moreton Bay, mixing with seawater over seagrass meadows dominated by eelgrass, Zostera muelleri. We sought to determine how low pH / high CO2 conditions near the spring affect these plants and their interactions with the black rabbitfish (Siganus fuscescens), a co-occurring grazer. In paired-choice feeding trials S. fuscescens preferentially consumed Z. muelleri shoots collected nearest to Myora Springs. Proximity to the spring did not significantly alter the carbon and nitrogen contents of seagrass tissues but did result in the extraordinary loss of soluble phenolics, including Folin-reactive phenolics, condensed tannins, and phenolic acids by ?87%. Conversely, seagrass lignin contents were, in this and related experiments, unaffected or increased, suggesting a shift in secondary metabolism away from the production of soluble, but not insoluble, (poly)phenolics. We suggest that groundwater discharge sites such as Myora Springs, and other sites characterized by low pH, are likely to be popular feeding grounds for seagrass grazers seeking to reduce their exposure to soluble phenolics.

Juvenile–adult tree associations in a continental Mediterranean ecosystem: no evidence for sustained and general facilitation at increased aridity

Question: How do tree species identity, microhabitat and water availability affect inter- and intra-specific interactions between juvenile and adult woody plants? Location: Continental Mediterranean forests in Alto Tajo Natural Park, Guadalajara, Spain. Methods: A total of 2066 juveniles and adults of four co-occurring tree species were mapped in 17 plots. The frequency of juveniles at different microhabitats and water availability levels was analysed using log-linear models. We used nearest-neighbour contingency table analysis of spatial segregation and J-functions to describe the spatial patterns. Results: We found a complex spatial pattern that varied according to species identity and microhabitat. Recruitment was more frequent in gaps for Quercus ilex, while the other three species recruited preferentially under shrubs or trees depending on the water availability level. Juveniles were not spatially associated to conspecific adults, experiencing segregation from them inmany cases. Spatial associations, both positive and negative, were more common at higher water availability levels. Conclusions: Our results do not agree with expectations from the stressgradient hypothesis, suggesting that positive interactions do not increase in importance with increasing aridity in the study ecosystem. Regeneration patterns are species-specific and depend on microhabitat characteristics and dispersal strategies. In general, juveniles do not look for conspecific adult protection. This work contributes to the understanding of species co-existence, proving the importance of considering a multispecies approach at several plots to overcome limitations of simple pair-wise comparisons in a limited number of sites.

Growth rate and TRI5 gene expression profiles of Fusarium equiseti strains isolated from Spanish cereals cultivated on wheat and barley media at different environmental conditions

Fusarium equiseti is a toxigenic species that often contaminates ce real crops from diverse climatic regions such as Northern and Southern Europe. Previous results suggested the existence of two distinct populations within this species with differences in toxin pro file which largely corresponded to North and South Europe (Spain). In this work, growth rate profiles of 4 F. equiseti strains isolated from different cereals and distinct Spanish regions were determined on wheat and barley based media at a range of temperatures (15, 20, 25, 30, 35 and 40 °C) and water potentialregimens(−0.7,−2.8,−7.0,and −9.8MPa,correspondingto 0.99,0.98,0.95 and 0.93aw values).Growth was observed at all temperatures except at 40 °C, and at all the solute potential values except at−9.8 MPa when combined with 15 °C. Optimal growth was observed at 20– 30 °C and −0.7/−2.8 MPa. The effect of these factors on trichothecene biosynthesis was examined on a F. equiseti strain using a newly developed real time RT-PCR protocol to quantify TRI5 gene expression at 15, 25 and 35 °C and −0.7, −2.8, − 7.0 and −9.8 MPa on wheat and barley based media. Induction of TRI5 expression was detected between 25 and 35 °C and −0.7 and − 2.8 MPa, with maximum values at 35 °C and −2.8 MPa being higher in barley than in wheat medium. These results appeared to be consistent with a population well adapted to the present climatic conditions and predicted scenarios for Southern Europe and suggested some differences depending on the cereal considered. These are also discussed in relation to other Fusarium species co-occurring in cereals grown in this region and to their significance for prediction and control strategies of toxigenic risk in future scenarios of climate change for this region.

Natural regeneration in mediterranean mixed forests: effect of environmental conditions on the physiological optimum niche of the species = Regeneración natural en masas mixtas mediterráneas: efecto de las condiciones ambientales en el nicho óptimo fisiológico de las especies

Los montes Mediterráneos han experimentado múltiples cambios en las últimas décadas (tanto en clima como en usos), lo que ha conducido a variaciones en la distribución de especies. El aumento previsto de las temperaturas medias junto con la mayor variabilidad intra e inter anual en cuanto a la ocurrencia de eventos extremos o disturbios naturales (como periodos prolongados de sequía, olas de frío o calor, incendios forestales o vendavales) pueden dañar significativamente al regenerado, llevándolo hasta la muerte, y jugando un papel decisivo en la composición de especies y en la dinámica del monte. La amplitud ecológica de muchas especies forestales puede verse afectada, de forma que se esperan cambios en sus nichos actuales de regeneración. Sin embargo, la migración latitudinal de las especies en busca de mejores condiciones, podría ser una explicación demasiado simplista de un proceso mucho más complejo de interacción entre la temperatura y la precipitación, que afectaría a cada especie de un modo distinto. En este sentido tanto la capacidad de adaptación al estrés ambiental de una determinada especie, así como su habilidad para competir por los recursos limitados, podría significar variaciones dentro de una comunidad. Las características fisiológicas y morfológicas propias de cada especie se encuentran fuertemente relacionadas con el lugar donde cada una puede surgir, qué especies pueden convivir y como éstas responden a las condiciones ambientales. En este sentido, el conocimiento sobre las distintas respuestas ecofisiológicas observadas ante cambios ambientales puede ser fundamentales para la predicción de variaciones en la distribución de especies, composición de la comunidad y productividad del monte ante el cambio global. En esta tesis investigamos el grado de tolerancia y sensibilidad que cada una de las tres especies de estudio, coexistentes en el interior peninsular ibérico (Pinus pinea, Quercus ilex y Juniperus oxycedrus), muestra ante los factores abióticos de estrés típicos de la región Mediterránea. Nuestro trabajo se ha basado en la definición del nicho óptimo fisiológico para el regenerado de cada especie a través de la investigación en profundidad del efecto de la sequía, la temperatura y el ambiente lumínico. Para ello, hemos desarrollado un modelo de predicción de la tasa de asimilación de carbono que nos ha permitido identificar las condiciones óptimas ambientales donde el regenerado de cada especie podría establecerse con mayor facilidad. En apoyo a este trabajo y con la idea de estudiar el efecto de la sequía a nivel de toda la planta hemos desarrollado un experimento paralelo en invernadero. Aquí se han aplicado dos regímenes hídricos para estudiar las características fisiológicas y morfológicas de cada especie, sobre todo a nivel de raíz y crecimiento del tallo, y relacionarlas con las diferentes estrategias en el uso del agua de las especies. Por último, hemos estudiado los patrones de aclimatación y desaclimatación al frio de cada especie, identificando los periodos de sensibilidad a heladas, así como cuellos de botella donde la competencia entre especies podría surgir. A pesar de que el pino piñonero ha sido la especie objeto de la gestión de estas masas durante siglos, actualmente se encuentra en la posición más desfavorable para combatir el cambio global, presentado el nicho fisiológico más estrecho de las tres especies. La encina sin embargo, ha resultado ser la especie mejor cualificada para afrontar este cambio, seguida muy de cerca por el enebro. Nuestros resultados sugieren una posible expansión en el rango de distribución de la encina, un aumento en la presencia del enebro y una disminución progresiva del pino piñonero a medio plazo en estas masas. ABSTRACT Mediterranean forests have undergone multiple changes over the last decades (in both climate and land use), which have lead to variations in the distribution of species. The expected increase in mean annual temperature together with the greater inter and intra-annual variability in extreme events and disturbances occurrence (such as prolonged drought periods, cold or heat waves, wildfires or strong winds) can significantly damage natural regeneration, up to causing death, playing a decisive role on species composition and forest dynamics. The ecological amplitude for adaptation of many species can be affected in such a way that changes in the current regeneration niches of many species are expected. However, the forecasted poleward migration of species seeking better conditions could be an oversimplification of what is a more complex phenomenon of interactions among temperature and precipitation, that would affect different species in different ways. In this regard, either the ability to adapt to environmental stresses or to compete for limited resources of a single species in a mixed forest could lead to variations within a community. The ecophysiological and morphological traits specific to each species are strongly related to the place where each species can emerge, which species can coexist, and how they respond to environmental conditions. In this regard, the understanding of the ecophysiological responses observed against changes in environmental conditions can be essential for predicting variations in species distribution, community composition, and forest productivity in the context of global change. In this thesis we investigated the degree of tolerance and sensitivity that each of the three studied species, co-occurring in central of the Iberian Peninsula (Pinus pinea, Quercus ilex and Juniperus oxycedrus), show against the typical abiotic stress factors in the Mediterranean region. Our work is based on the optimal physiological niche for regeneration of each species through in-depth research on the effect of drought, temperature and light environment. For this purpose, we developed a model to predict the carbon assimilation rate which allows us to identify the optimal environmental conditions where regeneration from each species could establish itself more easily. To obtain a better understanding about the effect of low temperature on regeneration, we studied the acclimation and deacclimation patterns to cold of each species, identifying period of frost sensitivity, as well as bottlenecks where competition between species can arise. Finally, to support our results about the effect of water availabilty, we conducted a greenhouse experiment with a view of studying the drought effect at the whole plant level. Here, two watering regimes were applied in order to study the physiological and morphological traits of each species, mainly at the level of the root system and stem growth, and so relate them to the different water use strategies of the species. Despite the fact that stone pine has been the target species for centuries, nowadays this species is in the most unfavorable position to cope with climate change. Holm oak, however, resulted the species that is best adapted to tolerate the predicted changes, followed closely by prickly juniper. Our results suggest a feasible expansion of the distribution range in holm oak, an increase in the prickly juniper presence and a progressive decreasing of stone pine presence in the medium term in these stone pine-holm oak-prickly juniper mixed forests.

Physical stress and diversity-productivity relationships: The role of positive interactions

If environmental stress provides conditions under which positive relationships between plant species richness and productivity become apparent, then species that seem functionally redundant under constant conditions may add to community functioning under variable conditions. Using naturally co-occurring mosses and liverworts, we constructed bryophyte communities to test relationships between species diversity (1, 2, 4, 8, 16, 24, or 32 species) and productivity under constant conditions and when exposed to experimental drought. We found no relationship between species richness and biomass under constant conditions. However, when communities were exposed to experimental drought, biomass increased with species richness. Responses of individual species demonstrated that facilitative interactions rather than sampling effects or niche complementarity best explained results—survivorship increased for almost all species, and those species least resistant to drought in monoculture had the greatest increase in biomass. Positive interactions may be an important but previously underemphasized mechanism linking high diversity to high productivity under stressful environmental conditions.

The comparative effectiveness of rodents and dung beetles as local seed dispersers in Mediterranean oak forests

The process of seed dispersal of many animal-dispersed plants is frequently mediated by a small set of biotic agents. However, the contribution that each of these dispersers makes to the overall recruitment may differ largely, with important ecological and management implications for the population viability and dynamics of the species implied in these interactions. In this paper, we compared the relative contribution of two local guilds of scatter-hoarding animals with contrasting metabolic requirements and foraging behaviours (rodents and dung beetles) to the overall recruitment of two Quercus species co-occurring in the forests of southern Spain. For this purpose, we considered not only the quantity of dispersed seeds but also the quality of the seed dispersal process. The suitability for recruitment of the microhabitats where the seeds were deposited was evaluated in a multi-stage demographic approach. The highest rates of seed handling and predation occurred in those microhabitats located under shrubs, mostly due to the foraging activity of rodents. However, the probability of a seed being successfully cached was higher in microhabitats located beneath a tree canopy as a result of the feeding behaviour of beetles. Rodents and beetles showed remarkable differences in their effectiveness as local acorn dispersers. Quantitatively, rodents were much more important than beetles because they dispersed the vast majority of acorns. However, they were qualitatively less effective because they consumed a high proportion of them (over 95%), and seeds were mostly dispersed under shrubs, a less suitable microhabitat for short-term recruitment of the two oak species. Our findings demonstrate that certain species of dung beetles (such as Thorectes lusitanicus), despite being quantitatively less important than rodents, can act as effective local seed dispersers of Mediterranean oak species. Changes in the abundance of beetle populations could thus have profound implications for oak recruitment and community dynamics.

Community attributes determine facilitation potential in a semi-arid steppe

Studies on positive plant–plant relations have traditionally focused on pair-wise interactions. Conversely, the interaction with other co-occurring species has scarcely been addressed, despite the fact that the entire community may affect plant performance. We used woody vegetation patches as models to evaluate community facilitation in semi-arid steppes. We characterized biotic and physical attributes of 53 woody patches (patch size, litter accumulation, canopy density, vegetation cover, species number and identity, and phylogenetic distance), and soil fertility (organic C and total N), and evaluated their relative importance for the performance of seedlings of Pistacia lentiscus, a keystone woody species in western Mediterranean steppes. Seedlings were planted underneath the patches, and on their northern and southern edges. Woody patches positively affected seedling survival but not seedling growth. Soil fertility was higher underneath the patches than elsewhere. Physical and biotic attributes of woody patches affected seedling survival, but these effects depended on microsite conditions. The composition of the community of small shrubs and perennial grasses growing underneath the patches controlled seedling performance. An increase in Stipa tenacissima and a decrease in Brachypodium retusum increased the probability of survival. The cover of these species and other small shrubs, litter depth and community phylogenetic distance, were also related to seedling survival. Seedlings planted on the northern edge of the patches were mostly affected by attributes of the biotic community. These traits were of lesser importance in seedlings planted underneath and in the southern edge of patches, suggesting that constraints to seedling establishment differed within the patches. Our study highlights the importance of taking into consideration community attributes over pair-wise interactions when evaluating the outcome of ecological interactions in multi-specific communities, as they have profound implications in the composition, function and management of semi-arid steppes.

Drought and grazing combined: Contrasting shifts in plant interactions at species pair and community level

The combined effects of drought stress and grazing pressure on shaping plant–plant interactions are still poorly understood, while this combination is common in arid ecosystems. In this study we assessed the relative effect of grazing pressure and slope aspect (drought stress) on vegetation cover and soil functioning in semi-arid Mediterranean grassland–shrublands in southeastern Spain. Moreover, we linked these two stress factors to plant co-occurrence patterns at species-pair and community levels, by performing C-score analyses. Vegetation cover and soil functioning decreased with higher grazing pressure and more south-facing (drier) slopes. At the community level, plants at south-facing slopes were negatively associated at no grazing but positively associated at low grazing pressure and randomly associated at high grazing pressure. At north-facing slopes, grazing did not result in a shift in the direction of the association. In contrast, analysis of pairwise species co-occurrence patterns showed that the dominant species Stipa tenacissima and Anthyllis cytisoides shifted from excluding each other to co-occurring with increasing grazing pressure at north-facing slopes. Our findings highlight that for improved understanding of plant interactions along stress gradients, interactions between species pairs and interactions at the community level should be assessed, as these may reveal contrasting results.