The European Union (EU) Horizon 2020 Coordination and Support Action ESMERALDA aimed at developing guidance and a flexible methodology for Mapping and Assessment of Ecosystems and their Services (MAES) to support the EU member states in the implementation of the EU Biodiversity Strategy’s Target 2 Action 5. ESMERALDA’s key tasks included network creation, stakeholder engagement, enhancing ecosystem services mapping and assessment methods across various spatial scales and value domains, work in case studies and support of EU member states in MAES implementation. Thus ESMERALDA aimed at integrating various project outcomes around four major strands: i) Networking, ii) Policy, iii) Research and iv) Application. The objective was to provide guidance for integrated ecosystem service mapping and assessment that can be used for sustainable decision-making in policy, business, society, practice and science at EU, national and regional levels. This article presents the overall ESMERALDA approach of integrating the above-mentioned project components and outcomes and provides an overview of how the enhanced methods were applied and how they can be used to support MAES implementation in the EU member states. Experiences with implementing such a large pan-European Coordination and Support Action in the context of EU policy are discussed and recommendations for future actions are given.

Premise Different cytotypes of a species may differ in their morphology, phenology, physiology, and their tolerance of extreme environments. We studied the ecological niches of two subspecies of Saxifraga rosacea with different ploidy levels: the hexaploid Central European endemic subspecies sponhemica and the more widely distributed octoploid subspecies rosacea. Methods For both cytotypes, we recorded local environmental conditions and mean plant trait values in populations across their areas of distribution, analyzed their distributions by niche modeling, studied their performance at two transplant sites with contrasting conditions, and experimentally tested their cold resistance. Results Mean annual temperature was higher in hexaploid than in octoploid populations and experiments indicated that frost tolerance of the hexaploid is lower than that of the octoploid. Reproduction of octoploids from Central Europe was higher than that of hexaploids at a transplant site in subarctic Iceland, whereas the opposite was true in temperate Luxembourg, indicating adaptation of the octoploids to colder conditions. Temperature variables were also most important in niche models predicting the distribution of the two cytotypes. Genetic differences in survival among populations were larger for the octoploids than for the hexaploids in both field gardens, suggesting that greater genetic variability may contribute to the octoploid's larger distributional range. Conclusions Our results support the hypotheses that different cytotypes may have different niches leading to spatial segregation, and that higher ploidy levels can result in a broader ecological niche and greater tolerance of more extreme conditions.

We examined recruitment, survival, life cycle and fecundity of two metallicolous (M, on metalliferous calamine soils) and two non‐metallicolous (NM, on normal soils) populations of Thlaspi caerulescens in Belgium and Luxemburg. In each population, permanent plots were monitored over two reproductive seasons. In M populations, plots were located in two contrasting environments (grass versus grove) in order to test the influence of vegetation cover on life strategy. Our results show that the monocarpic life cycle is dominant in all populations of T. caerulescens. However the length of the pre‐reproductive period varies from several months (winter annuals) to 1 year or more (perennials), and is partly related to plant origin (M versus NM). Most plants growing in metalliferous environments were annuals, whereas NM plants were mostly perennials. These differences in life cycle were related to differences in survival during summer, which was better in NM than in M populations. Within each M population, different survival conditions and life cycles were observed according to vegetation cover. Plants growing in grass areas were mostly annuals and had a low survival rate in summer whereas grove plants were mostly perennials and survived better in summer. Our results suggest the selection of stress avoiders (shortening of life cycle) in M populations of T. caerulescens but only for individuals growing in grass areas. Summer survival seems to play a key role in selection of life strategy in T. caerulescens.

Summary1. The Bayesian inference of demographic parameters under an Isolation-Migration (IM) model of population evolution offers a major improvement over previously available approaches. This method is implemented in a popular program, IMa, widely used in population genetic studies. 2. While the robustness of the method to deviations of the IM model has previously been evaluated, we assess the performance of the program with two populations when the model used to generate the analysed data meets the assumptions of the IM model completely; the goal is to identify the conditions under which the method works best. Overall, we test eighteen sets of conditions and analyse AE 500 simulated data sets, for a total of over 200,000 hours of analyses using a large computer cluster. 3. Although we find clear differences in quality estimates among models, the best ranges of demographic parameter values to infer accurate estimates differ among parameters. Divergence time is best estimated in the absence of gene flow and when population sizes are large compared to divergence time. In contrast, the classic population parameter ϑ (= 4Nl) is best estimated, for the two current populations, when divergence time is large compared to population size, with or without migration. The parameter is always poorly estimated in the case of the ancestral population. While it is possible to distinguish between scenarios with or without gene flow, estimating the extent of gene flow, when different from 0, is associated with relatively high error rates. In general, increasing the number of loci or the sample size reduces the variance and credible interval of the estimates, and only for the migration rate, it slightly improves the accuracy of the estimate as well. Increasing the prior distribution range of a parameter can dramatically increase that of its posterior distribution. Surprisingly, differences are highlighted among the estimates inferred from sequences generated by different simulation programs, especially for the simulation program SIMDIV. 4. Overall, the performances of the method shown here probably reflect the limitation of the method in general and/or of the historical information contained in DNA sequence data.