The relationship between remotely-sensed spectral heterogeneity and bird diversity is modulated by landscape type.

To identify areas of high biodiversity and prioritize conservation efforts, it is crucial to understand the drivers of species richness patterns and their scale dependence. While classified land cover products are commonly used to explain bird species richness, recent studies suggest that unclassified remote-sensed images can provide equally good or better results. In our study, we aimed to investigate whether unclassified multispectral data from Landsat 8 can replace image classification for bird diversity modeling. Moreover, we also tested the Spectral Variability Hypothesis. Using the Atlas of Breeding Birds in the Czech Republic 2014-2017, we modeled species richness at two spatial resolutions of approx. 131 km2 (large squares) and 8 km2 (small squares). As predictors of the richness, we assessed 1) classified land cover data (Corine Land Cover 2018 database), 2) spectral heterogeneity (computed in three ways) and landscape composition derived from unclassified remote-sensed reflectance and vegetation indices. Furthermore, we integrated information about the landscape types (expressed by the most prevalent land cover class) into models based on unclassified remote-sensed data to investigate whether the landscape type plays a role in explaining bird species richness. We found that unclassified remote-sensed data, particularly spectral heterogeneity metrics, were better predictors of bird species richness than classified land cover data. The best results were achieved by models that included interactions between the unclassified data and landscape types, indicating that relationships between bird diversity and spectral heterogeneity vary across landscape types. Our findings demonstrate that spectral heterogeneity derived from unclassified multispectral data is effective for assessing bird diversity across the Czech Republic. When explaining bird species richness, it is important to account for the type of landscape and carefully consider the significance of the chosen spatial scale.

Sampling and modelling rare species: Conceptual guidelines for the neglected majority.

Biodiversity conservation faces a methodological conundrum: Biodiversity measurement often relies on species, most of which are rare at various scales, especially prone to extinction under global change, but also the most challenging to sample and model. Predicting the distribution change of rare species using conventional species distribution models is challenging because rare species are hardly captured by most survey systems. When enough data are available, predictions are usually spatially biased towards locations where the species is most likely to occur, violating the assumptions of many modelling frameworks. Workflows to predict and eventually map rare species distributions imply important trade-offs between data quantity, quality, representativeness and model complexity that need to be considered prior to survey and analysis. Our opinion is that study designs need to carefully integrate the different steps, from species sampling to modelling, in accordance with the different types of rarity and available data in order to improve our capacity for sound assessment and prediction of rare species distribution. In this article, we summarize and comment on how different categories of species rarity lead to different types of occurrence and distribution data depending on choices made during the survey process, namely the spatial distribution of samples (where to sample) and the sampling protocol in each selected location (how to sample). We then clarify which species distribution models are suitable depending on the different types of distribution data (how to model). Among others, for most rarity forms, we highlight the insights from systematic species-targeted sampling coupled with hierarchical models that allow correcting for overdispersion and spatial and sampling sources of bias. Our article provides scientists and practitioners with a much-needed guide through the ever-increasing diversity of methodological developments to improve the prediction of rare species distribution depending on rarity type and available data.

Integrating field and satellite data for spatially explicit inference on the density of threatened arboreal primates.

Spatially explicit models of animal abundance are a critical tool to inform conservation planning and management. However, they require the availability of spatially diffuse environmental predictors of abundance, which may be challenging, especially in complex and heterogeneous habitats. This is particularly the case for tropical mammals, such as nonhuman primates, that depend on multi-layered and species-rich tree canopy coverage, which is usually measured through a limited sample of ground plots. We developed an approach that calibrates remote-sensing imagery to ground measurements of tree density to derive basal area, in turn used as a predictor of primate density based on published models. We applied generalized linear models (GLM) to relate 9.8-ha ground samples of tree basal area to various metrics extracted from Landsat 8 imagery. We tested the potential of this approach for spatial inference of animal density by comparing the density predictions for an endangered colobus monkey, to previous estimates from field transect counts, measured basal area, and other predictors of abundance. The best GLM had high accuracy and showed no significant difference between predicted and observed values of basal area. Our species distribution model yielded predicted primate densities that matched those based on field measurements. Results show the potential of using open-access and global remote-sensing data to derive an important predictor of animal abundance in tropical forests and in turn to make spatially explicit inference on animal density. This approach has important, inherent applications as it greatly magnifies the relevance of abundance modeling for informing conservation. This is especially true for threatened species living in heterogeneous habitats where spatial patterns of abundance, in relation to habitat and/or human disturbance factors, are often complex and, management decisions, such as improving forest protection, may need to be focused on priority areas.

High habitat richness reduces the risk of tick-borne encephalitis in Europe: A multi-scale study.

Background: The natural transmission cycle of tick-borne encephalitis (TBE) virus is enhanced by complex interactions between ticks and key hosts strongly connected to habitat characteristics. The diversity of wildlife host species and their relative abundance is known to affect transmission of tick-borne diseases. Therefore, in the current context of global biodiversity loss, we explored the relationship between habitat richness and the pattern of human TBE cases in Europe to assess biodiversity's role in disease risk mitigation. Methods: We assessed human TBE case distribution across 879 European regions using official epidemiological data reported to The European Surveillance System (TESSy) between 2017 and 2021 from 15 countries. We explored the relationship between TBE presence and the habitat richness index (HRI1) by means of binomial regression. We validated our findings at local scale using data collected between 2017 and 2021 in 227 municipalities located in Trento and Belluno provinces, two known TBE foci in northern Italy. Findings: Our results showed a significant parabolic effect of HRI on the probability of presence of human TBE cases in the European regions included in our dataset, and a significant, negative effect of HRI on the local presence of TBE in northern Italy. At both spatial scales, TBE risk decreases in areas with higher values of HRI. Interpretation: To our knowledge, no efforts have yet been made to explore the relationship between biodiversity and TBE risk, probably due to the scarcity of high-resolution, large-scale data about the abundance or density of critical host species. Hence, in this study we considered habitat richness as proxy for vertebrate host diversity. The results suggest that in highly diverse habitats TBE risk decreases. Hence, biodiversity loss could enhance TBE risk for both humans and wildlife. This association is relevant to support the hypothesis that the maintenance of highly diverse ecosystems mitigates disease risk.

Grassland vertical height heterogeneity predicts flower and bee diversity: an UAV photogrammetric approach.

The ecosystem services offered by pollinators are vital for supporting agriculture and ecosystem functioning, with bees standing out as especially valuable contributors among these insects. Threats such as habitat fragmentation, intensive agriculture, and climate change are contributing to the decline of natural bee populations. Remote sensing could be a useful tool to identify sites of high diversity before investing into more expensive field survey. In this study, the ability of Unoccupied Aerial Vehicles (UAV) images to estimate biodiversity at a local scale has been assessed while testing the concept of the Height Variation Hypothesis (HVH). This hypothesis states that the higher the vegetation height heterogeneity (HH) measured by remote sensing information, the higher the vegetation vertical complexity and the associated species diversity. In this study, the concept has been further developed to understand if vegetation HH can also be considered a proxy for bee diversity and abundance. We tested this approach in 30 grasslands in the South of the Netherlands, where an intensive field data campaign (collection of flower and bee diversity and abundance) was carried out in 2021, along with a UAV campaign (collection of true color-RGB-images at high spatial resolution). Canopy Height Models (CHM) of the grasslands were derived using the photogrammetry technique "Structure from Motion" (SfM) with horizontal resolution (spatial) of 10 cm, 25 cm, and 50 cm. The accuracy of the CHM derived from UAV photogrammetry was assessed by comparing them through linear regression against local CHM LiDAR (Light Detection and Ranging) data derived from an Airborne Laser Scanner campaign completed in 2020/2021, yielding an [Formula: see text] of 0.71. Subsequently, the HH assessed on the CHMs at the three spatial resolutions, using four different heterogeneity indices (Rao's Q, Coefficient of Variation, Berger-Parker index, and Simpson's D index), was correlated with the ground-based flower and bee diversity and bee abundance data. The Rao's Q index was the most effective heterogeneity index, reaching high correlations with the ground-based data (0.44 for flower diversity, 0.47 for bee diversity, and 0.34 for bee abundance). Interestingly, the correlations were not significantly influenced by the spatial resolution of the CHM derived from UAV photogrammetry. Our results suggest that vegetation height heterogeneity can be used as a proxy for large-scale, standardized, and cost-effective inference of flower diversity and habitat quality for bees.

LiDAR GEDI derived tree canopy height heterogeneity reveals patterns of biodiversity in forest ecosystems.

The "Height Variation Hypothesis" is an indirect approach used to estimate forest biodiversity through remote sensing data, stating that greater tree height heterogeneity (HH) measured by CHM LiDAR data indicates higher forest structure complexity and tree species diversity. This approach has traditionally been analyzed using only airborne LiDAR data, which limits its application to the availability of the dedicated flight campaigns. In this study we analyzed the relationship between tree species diversity and HH, calculated with four different heterogeneity indices using two freely available CHMs derived from the new space-borne GEDI LiDAR data. The first, with a spatial resolution of 30 m, was produced through a regression tree machine learning algorithm integrating GEDI LiDAR data and Landsat optical information. The second, with a spatial resolution of 10 m, was created using Sentinel-2 images and a deep learning convolutional neural network. We tested this approach separately in 30 forest plots situated in the northern Italian Alps, in 100 plots in the forested area of Traunstein (Germany) and successively in all the 130 plots through a cross-validation analysis. Forest density information was also included as influencing factor in a multiple regression analysis. Our results show that the GEDI CHMs can be used to assess biodiversity patterns in forest ecosystems through the estimation of the HH that is correlated to the tree species diversity. However, the results also indicate that this method is influenced by different factors including the GEDI CHMs dataset of choice and their related spatial resolution, the heterogeneity indices used to calculate the HH and the forest density. Our finding suggest that GEDI LIDAR data can be a valuable tool in the estimation of forest tree heterogeneity and related tree species diversity in forest ecosystems, which can aid in global biodiversity estimation.

The Spectral Species Concept in Living Color.

Biodiversity monitoring is an almost inconceivable challenge at the scale of the entire Earth. The current (and soon to be flown) generation of spaceborne and airborne optical sensors (i.e., imaging spectrometers) can collect detailed information at unprecedented spatial, temporal, and spectral resolutions. These new data streams are preceded by a revolution in modeling and analytics that can utilize the richness of these datasets to measure a wide range of plant traits, community composition, and ecosystem functions. At the heart of this framework for monitoring plant biodiversity is the idea of remotely identifying species by making use of the 'spectral species' concept. In theory, the spectral species concept can be defined as a species characterized by a unique spectral signature and thus remotely detectable within pixel units of a spectral image. In reality, depending on spatial resolution, pixels may contain several species which renders species-specific assignment of spectral information more challenging. The aim of this paper is to review the spectral species concept and relate it to underlying ecological principles, while also discussing the complexities, challenges and opportunities to apply this concept given current and future scientific advances in remote sensing.

Terra and Aqua satellites track tiger mosquito invasion: modelling the potential distribution of Aedes albopictus in north-eastern Italy.

BACKGROUND: The continuing spread of the Asian tiger mosquito Aedes albopictus in Europe is of increasing public health concern due to the potential risk of new outbreaks of exotic vector-borne diseases that this species can transmit as competent vector. We predicted the most favorable areas for a short term invasion of Ae. albopictus in north-eastern Italy using reconstructed daily satellite data time series (MODIS Land Surface Temperature maps, LST). We reconstructed more than 11,000 daily MODIS LST maps for the period 2001-09 (i.e. performed spatial and temporal gap-filling) in an Open Source GIS framework. We aggregated these LST maps over time and identified the potential distribution areas of Ae. albopictus by adapting published temperature threshold values using three variables as predictors (0°C for mean January temperatures, 11°C for annual mean temperatures and 1350 growing degree days filtered for areas with autumnal mean temperatures > 11°C). The resulting maps were integrated into the final potential distribution map and this was compared with the known current distribution of Ae. albopictus in north-eastern Italy. RESULTS: LST maps show the microclimatic characteristics peculiar to complex terrains, which would not be visible in maps commonly derived from interpolated meteorological station data. The patterns of the three indicator variables partially differ from each other, while winter temperature is the determining limiting factor for the distribution of Ae. albopictus. All three variables show a similar spatial pattern with some local differences, in particular in the northern part of the study area (upper Adige valley). CONCLUSIONS: Reconstructed daily land surface temperature data from satellites can be used to predict areas of short term invasion of the tiger mosquito with sufficient accuracy (200 m pixel resolution size). Furthermore, they may be applied to other species of arthropod of medical interest for which temperature is a relevant limiting factor. The results indicate that, during the next few years, the tiger mosquito will probably spread toward northern latitudes and higher altitudes in north-eastern Italy, which will considerably expand the range of the current distribution of this species.

Diversity of European habitat types is correlated with geography more than climate and human pressure.

Habitat richness, that is, the diversity of ecosystem types, is a complex, spatially explicit aspect of biodiversity, which is affected by bioclimatic, geographic, and anthropogenic variables. The distribution of habitat types is a key component for understanding broad-scale biodiversity and for developing conservation strategies. We used data on the distribution of European Union (EU) habitats to answer the following questions: (i) how do bioclimatic, geographic, and anthropogenic variables affect habitat richness? (ii) Which of those factors is the most important? (iii) How do interactions among these variables influence habitat richness and which combinations produce the strongest interactions? The distribution maps of 222 terrestrial habitat types as defined by the Natura 2000 network were used to calculate habitat richness for the 10 km × 10 km EU grid map. We then investigated how environmental variables affect habitat richness, using generalized linear models, generalized additive models, and boosted regression trees. The main factors associated with habitat richness were geographic variables, with negative relationships observed for both latitude and longitude, and a positive relationship for terrain ruggedness. Bioclimatic variables played a secondary role, with habitat richness increasing slightly with annual mean temperature and overall annual precipitation. We also found an interaction between anthropogenic variables, with the combination of increased landscape fragmentation and increased population density strongly decreasing habitat richness. This is the first attempt to disentangle spatial patterns of habitat richness at the continental scale, as a key tool for protecting biodiversity. The number of European habitats is related to geography more than climate and human pressure, reflecting a major component of biogeographical patterns similar to the drivers observed at the species level. The interaction between anthropogenic variables highlights the need for coordinated, continental-scale management plans for biodiversity conservation.

rasterdiv-An Information Theory tailored R package for measuring ecosystem heterogeneity from space: To the origin and back.

Ecosystem heterogeneity has been widely recognized as a key ecological indicator of several ecological functions, diversity patterns and change, metapopulation dynamics, population connectivity or gene flow.In this paper, we present a new R package-rasterdiv-to calculate heterogeneity indices based on remotely sensed data. We also provide an ecological application at the landscape scale and demonstrate its power in revealing potentially hidden heterogeneity patterns.The rasterdiv package allows calculating multiple indices, robustly rooted in Information Theory, and based on reproducible open-source algorithms.

Priority list of biodiversity metrics to observe from space.

Monitoring global biodiversity from space through remotely sensing geospatial patterns has high potential to add to our knowledge acquired by field observation. Although a framework of essential biodiversity variables (EBVs) is emerging for monitoring biodiversity, its poor alignment with remote sensing products hinders interpolation between field observations. This study compiles a comprehensive, prioritized list of remote sensing biodiversity products that can further improve the monitoring of geospatial biodiversity patterns, enhancing the EBV framework and its applicability. The ecosystem structure and ecosystem function EBV classes, which capture the biological effects of disturbance as well as habitat structure, are shown by an expert review process to be the most relevant, feasible, accurate and mature for direct monitoring of biodiversity from satellites. Biodiversity products that require satellite remote sensing of a finer resolution that is still under development are given lower priority (for example, for the EBV class species traits). Some EBVs are not directly measurable by remote sensing from space, specifically the EBV class genetic composition. Linking remote sensing products to EBVs will accelerate product generation, improving reporting on the state of biodiversity from local to global scales.

Spectral variation versus species beta-diversity at different spatial scales: a test in African highland savannas.

Few studies exist that explicitly analyse the effect of grain, i.e. the sampling unit dimension, on vascular plant species turnover (beta-diversity) among sites. While high beta-diversity is often a result of high environmental heterogeneity, remotely sensed spectral distances among sampling units may be used as a proxy of environmental gradients which spatially shape the patterns of species turnover. In this communication, we aimed to (i) test the potential relation between spectral variation and species beta-diversity in a savanna environment and to (ii) investigate the effect of grain on the achieved patterns. Field data gathered by the BIOTA Southern Africa biodiversity monitoring programme were used to model the relation between spectral variation and species turnover at different spatial grains (10 m x 10 m and 20 m x 50 m). Our results indicate that the overall fit was greater at the larger grain size, confirming the theoretical assumption that using a lower grain size would generally lead to a higher noise in the calculation of species turnover. This communication represents one of the first attempts at relating beta-diversity to spectral variation, while incorporating the effects of grain size in the study. The results of this study could have significant implications for biodiversity research and conservation planning at a regional or even larger spatial scale.

A grid-based map for the Biogeographical Regions of Europe.

BACKGROUND: Biogeographical units are widely adopted in ecological research and nature conservation management, even though biogeographical regionalisation is still under scientific debate. The European Environment Agency provided an official map of the European Biogeographical Regions (EBRs), which contains the official boundaries used in the Habitats and Birds Directives. However, these boundaries bisect cells in the official EU 10 km × 10 km grid used for many purposes, including reporting species and habitat data, meaning that 6881 cells overlap two or more regions. Therefore, superimposing the EBRs vector map over the grid creates ambiguities in associating some cells with European Biogeographical Regions. NEW INFORMATION: To provide an operational tool to unambiguously define the boundaries of the eleven European Biogeographical Regions, we provide a specifically developed raster map of Grid-Based European Biogeographical Regions (GB-EBRs). In this new map, the borders of the EBRs are reshaped to coherently match the standard European 10 km × 10 km grid imposed for reporting tasks by Article 17 of the Habitats Directive and used for many other datasets. We assign each cell to the EBR with the largest area within the cell.

Geographical Distribution of Ljungan Virus in Small Mammals in Europe.

Ljungan virus (LV), which belongs to the Parechovirus genus in the Picornaviridae family, was first isolated from bank voles (Myodes glareolus) in Sweden in 1998 and proposed as a zoonotic agent. To improve knowledge of the host association and geographical distribution of LV, tissues from 1685 animals belonging to multiple rodent and insectivore species from 12 European countries were screened for LV-RNA using reverse transcriptase (RT)-PCR. In addition, we investigated how the prevalence of LV-RNA in bank voles is associated with various intrinsic and extrinsic factors. We show that LV is widespread geographically, having been detected in at least one host species in nine European countries. Twelve out of 21 species screened were LV-RNA PCR positive, including, for the first time, the red vole (Myodes rutilus) and the root or tundra vole (Alexandromys formerly Microtus oeconomus), as well as in insectivores, including the bicolored white-toothed shrew (Crocidura leucodon) and the Valais shrew (Sorex antinorii). Results indicated that bank voles are the main rodent host for this virus (overall RT-PCR prevalence: 15.2%). Linear modeling of intrinsic and extrinsic factors that could impact LV prevalence showed a concave-down relationship between body mass and LV occurrence, so that subadults had the highest LV positivity, but LV in older animals was less prevalent. Also, LV prevalence was higher in autumn and lower in spring, and the amount of precipitation recorded during the 6 months preceding the trapping date was negatively correlated with the presence of the virus. Phylogenetic analysis on the 185 base pair species-specific sequence of the 5' untranslated region identified high genetic diversity (46.5%) between 80 haplotypes, although no geographical or host-specific patterns of diversity were detected.

Explicitly accounting for pixel dimension in calculating classical and fractal landscape shape metrics.

Different summarized shape indices, like mean shape index (MSI) and area weighted mean shape index (AWMSI) can change over multiple size scales. This variation is important to describe scale heterogeneity of landscapes, but the exact mathematical form of the dependence is rarely known. In this paper, the use of fractal geometry (by the perimeter and area Hausdorff dimensions) made us able to describe the scale dependence of these indices. Moreover, we showed how fractal dimensions can be deducted from existing MSI and AWMSI data. In this way, the equality of a multiscale tabulated MSI and AWMSI dataset and two scale-invariant fractal dimensions has been demonstrated.

Algorithmic foundation of spectral rarefaction for measuring satellite imagery heterogeneity at multiple spatial scales.

Measuring heterogeneity in satellite imagery is an important task to deal with. Most measures of spectral diversity have been based on Shannon Information theory. However, this approach does not inherently address different scales, ranging from local (hereafter referred to alpha diversity) to global scales (gamma diversity). The aim of this paper is to propose a method for measuring spectral heterogeneity at multiple scales based on rarefaction curves. An algorithmic solution of rarefaction applied to image pixel values (Digital Numbers, DNs) is provided and discussed.

Advances in Microclimate Ecology Arising from Remote Sensing.

Microclimates at the land-air interface affect the physiological functioning of organisms which, in turn, influences the structure, composition, and functioning of ecosystems. We review how remote sensing technologies that deliver detailed data about the structure and thermal composition of environments are improving the assessment of microclimate over space and time. Mapping landscape-level heterogeneity of microclimate advances our ability to study how organisms respond to climate variation, which has important implications for understanding climate-change impacts on biodiversity and ecosystems. Interpolating in situ microclimate measurements and downscaling macroclimate provides an organism-centered perspective for studying climate-species interactions and species distribution dynamics. We envisage that mapping of microclimate will soon become commonplace, enabling more reliable predictions of species and ecosystem responses to global change.

A suite of essential biodiversity variables for detecting critical biodiversity change.

Key global indicators of biodiversity decline, such as the IUCN Red List Index and the Living Planet Index, have relatively long assessment intervals. This means they, due to their inherent structure, function as late-warning indicators that are retrospective, rather than prospective. These indicators are unquestionably important in providing information for biodiversity conservation, but the detection of early-warning signs of critical biodiversity change is also needed so that proactive management responses can be enacted promptly where required. Generally, biodiversity conservation has dealt poorly with the scattered distribution of necessary detailed information, and needs to find a solution to assemble, harmonize and standardize the data. The prospect of monitoring essential biodiversity variables (EBVs) has been suggested in response to this challenge. The concept has generated much attention, but the EBVs themselves are still in development due to the complexity of the task, the limited resources available, and a lack of long-term commitment to maintain EBV data sets. As a first step, the scientific community and the policy sphere should agree on a set of priority candidate EBVs to be developed within the coming years to advance both large-scale ecological research as well as global and regional biodiversity conservation. Critical ecological transitions are of high importance from both a scientific as well as from a conservation policy point of view, as they can lead to long-lasting biodiversity change with a high potential for deleterious effects on whole ecosystems and therefore also on human well-being. We evaluated candidate EBVs using six criteria: relevance, sensitivity to change, generalizability, scalability, feasibility, and data availability and provide a literature-based review for eight EBVs with high sensitivity to change. The proposed suite of EBVs comprises abundance, allelic diversity, body mass index, ecosystem heterogeneity, phenology, range dynamics, size at first reproduction, and survival rates. The eight candidate EBVs provide for the early detection of critical and potentially long-lasting biodiversity change and should be operationalized as a priority. Only with such an approach can science predict the future status of global biodiversity with high certainty and set up the appropriate conservation measures early and efficiently. Importantly, the selected EBVs would address a large range of conservation issues and contribute to a total of 15 of the 20 Aichi targets and are, hence, of high biological relevance.

Estimates of ozone AOT40 from passive sampling in forest sites in South-Western Europe.

Weekly-fortnightly ozone (O3) concentrations measured by passive sampling at 81 forest monitoring plots in France, Italy, Spain and Switzerland over the period 2000-2002 were used to estimate the cumulative exposure index AOT40. The estimation method is based on a deterministic model which describes the O3 daily profile as a function of relative altitude (the difference between the altitude of the site and the lowest altitude within a 5 km radius) and the time of the day. Estimated AOT40 values (AOT40(e)) were evaluated against co-located automatic measurement stations and with 14 independent automatic stations located throughout Italy whose weekly mean O3 values were used to simulate passive samplers. AOT40 can be predicted by modelling passive sampling data (R2: 0.90; P<0.0001, SE of estimates: 3271 ppb h), although considerable deviations can occur for individual sites. Estimated AOT40 shows a distinct, significant latitudinal and altitudinal gradient. Taking the 3-year average as a whole, exceedance of critical level of 5000 ppb h occurs at 77-100% of the monitored sites, respectively.

Invasion by the Alien Tree Prunus serotina Alters Ecosystem Functions in a Temperate Deciduous Forest.

Alien invasive species can affect large areas, often with wide-ranging impacts on ecosystem structure, function, and services. Prunus serotina is a widespread invader of European temperate forests, where it tends to form homogeneous stands and limits recruitment of indigenous trees. We hypotesized that invasion by P. serotina would be reflected in the nutrient contents of the native species' leaves and in the respiration of invaded plots as efficient resource uptake and changes in nutrient cycling by P. serotina probably underly its aggressive invasiveness. We combined data from 48 field plots in the forest of Compiègne, France, and data from an experiment using 96 microcosms derived from those field plots. We used general linear models to separate effects of invasion by P. serotina on heterotrophic soil and litter respiration rates and on canopy foliar nutrient content from effects of soil chemical properties, litter quantity, litter species composition, and tree species composition. In invaded stands, average respiration rates were 5.6% higher for soil (without litter) and 32% higher for soil and litter combined. Compared to indigenous tree species, P. serotina exhibited higher foliar N (+24.0%), foliar P (+50.7%), and lower foliar C:N (-22.4%) and N:P (-10.1%) ratios. P. serotina affected foliar nutrient contents of co-occuring indigenous tree species leading to decreased foliar N (-8.7 %) and increased C:N ratio (+9.5%) in Fagus sylvatica, decreased foliar N:P ratio in Carpinus betulus (-13.5%) and F. sylvatica (-11.8%), and increased foliar P in Pinus sylvestris (+12.3%) in invaded vs. uninvaded stands. Our results suggest that P. serotina is changing nitrogen, phosphorus, and carbon cycles to its own advantage, hereby increasing carbon turnover via labile litter, affecting the relative nutrient contents in the overstory leaves, and potentially altering the photosynthetic capacity of the long-lived indigenous broadleaved species. Uncontrolled invasion of European temperate forests by P. serotina may affect the climate change mitigation potential of these forests in the long term, through additive effects on local nutrient cycles.