Integrating habitat- and species-based perspectives for wetland conservation in lowland agricultural landscapes

Wetlands are among the most endangered ecosystems worldwide with multiple direct and indirect stressors, especially in human-altered areas like intensive agricultural landscapes. Conservation management and efforts often focus on species diversity and charismatic taxa, but scarcely consider habitats. By focusing on a complex formed by 107 permanent wetlands at 18 Natura 2000 sites in the Emilia-Romagna region (northern Italy), the patterns of habitats of conservation concern were investigated and the concordance with threatened species patterns was analysed. Wetlands were characterised in terms of morphology, connectivity, land use and management as drivers of assemblage and richness patterns of habitats. Our results showed a strong concordance between the distribution and richness patterns of both habitats and threatened taxa (birds, mammals, amphibians, reptiles, fish, invertebrates, and plants). Thus, habitats seem an effective proxy of species patterns. The variables related with perimeter, environmental heterogeneity and presence of water bodies were the most important ones associated with habitat richness patterns. The presence of aquatic systems (measured as the percentage of wetland area occupied by an aquatic surface) and their position in the hydrographic network were associated mostly with habitats distribution. Low richness wetlands (in habitat terms) were not complementary as no new habitat types were supported. The results stressed the relevance of wetlands with wide water body perimeters composed of diverse systems as being key for biodiversity conservation in a simplified agricultural matrix. Integrating habitat- and species-based perspectives seems a promising field and may provide a rapid assessment tool to acquire effective information for wetlands conservation and assessment.


Introduction
Wetlands are critical ecosystems for the biodiversity they support and the ecosystem services they provide (De Groot et al. 2012;Junk et al. 2013). At the same time however, these diversified ecosystems suffer numerous and multiple impacts, even in the climate change context (Cížková et al. 2013;Junk et al. 2013). This is especially true for lowland agricultural landscape, where wetlands and small standing-water ecosystems are residual hotspots of local biodiversity and ecosystem functioning, despite their progressive and global disappearance (e.g. Zedler 2003;Bolpagni et al. 2019). Hence, in heavily exploited agricultural landscapes, wetlands represent unique "natural islands" in a dense anthropogenic matrix (Lastrucci et al. 2010;Bolpagni and Piotti 2016).
Measuring and monitoring biodiversity are crucial tasks for wetlands conservation and management, also for designating effective protected areas and ecological restoration plans. However, evaluating the conservation status and monitoring the biodiversity of wetlands is a complex matter. Surveying biodiversity may present certain difficulties, such as limited financial resources and/or time, and requiring considerable expertise for different components. For these reasons, surrogates, proxies or indicators, herein used as synonymous, as species surrogates (e.g. Sánchez-Fernández et al. 2006;Caro 2010) or environmental surrogates (e.g. Faith 2003; Beier and de Albuquerque 2015) have been tested or suggested in numerous geographic areas and ecosystems (for complete definitions, see Hunter Jr et al. 2016). Given the strict relationship linking terrestrial, riparian and aquatic habitats in wetlands, different kinds of indicators have been tested including plants (Rogers et al. 2012;Chamberlain and Brooks 2016), invertebrates (Kashian and Burton 2000;Ormerod et al. 2010;Guareschi et al. 2015a) and vertebrates (e.g. Frederick et al. 2009;Guareschi et al. 2015b). However, many authors have stressed contrasting findings after using surrogates or context-specific responses (e.g. Araújo et al. 2001;Lindenmayer et al. 2015;Pierson et al. 2015). To date, conservation and restoration efforts often focus on specific plants or charismatic animals (e.g. species red lists), and a limited integration between both approaches has been recorded (McAlpine et al. 2016).
In global terms, the Ramsar convention on wetlands (Ramsar Bureau 2000) designates sites of international importance (more than 2300 worldwide) following specific requisites based on animals, plants or ecological communities in general, plus specific criteria on vertebrates (fish and waterbirds) in particular. In Europe, the Natura 2000 (N2000) network represents the key tool for biodiversity conservation, following Habitats Directive 92/43/ EC (HD) and Birds Directives 79/409/EC and 2009/147/EC (BD). Nationwide, each country also identifies and selects specific local taxa or habitats of conservation concern (e.g. in Italy by Bolpagni et al. 2010;Genovesi et al. 2014).
Vegetation-based indicators have been used to assess wetland restoration and conservation priorities (Matthews et al. 2009;Benavent-González et al. 2014;Angiolini et al. 2017), by considering that several animal communities can be affected by vegetation structures/ systems (Schwab et al. 2002;Zellweger et al. 2017). In fact, vegetation plays a pivotal role in shaping physical habitats, especially in aquatic contexts where aquatic plants act as engineer species deeply modifying the colonized environments (Pierobon et al. 2010;Ribaudo et al. 2018).
For this reason, vegetation has been used in environmental legislation as a mainstay to define habitats of conservation priority (e.g. sensu HD). In the biodiversity assessment and conservation planning context, habitats are increasingly being tested for species conservation (Bunce et al. 2013), and their relevance has already been specifically stressed for marine (Dalleau et al. 2010) and terrestrial ecosystems (Oliver et al. 2004). Nevertheless, little is known about their roles in inland wetlands, mainly in lowland agricultural landscapes (but see an example in a Mediterranean wetland network, Angiolini et al. 2017). Indeed, the relationship between livestock, farming, irrigation and nature conservation is a critical and crucial topic, with open debate in biological conservation (Bolpagni and Piotti 2016). The need for conservation efforts in densely populated areas has been stressed by Ricketts and Imhoff (2003) in North America as being just as important as preserving pristine regions. Moreover, the relevance of agri-environment schemes and agricultural matrices for constructive biodiversity conservation has been pointed out in Europe and elsewhere (Vandermeer and Perfecto 2007;Batáry et al. 2015;Bolpagni and Piotti 2016).
In order to contribute to the biological conservation of wetlands and to provide pivotal information for protected area management and future definitions, we herein integrated habitat-and species-based perspectives to better understand the contribution of habitats (sensu HD) to wetland biodiversity in a heavily exploited agricultural landscape in northern Italy.
This work specifically aimed to study and test: (i) if European and national habitats of conservation interest (HCI hereinafter) can act as a proxy of species of conservation concern (SCC hereinafter), in both richness and composition terms; (ii) the main environmental descriptors associated with the richness and distribution of HCI; and (iii) the role of protected areas with low HCI richness values as potential complementary sources of representativeness and conservation interest. Our hypotheses were that HCI would serve as indicators of wealth in rare and protected species (SCC), and that morphological and ecological drivers, related mainly to the morphological complexity of wetlands, could explain the observed patterns.
To address the first point, we assess whether HCI could be used as rapid tools for assessing wetland biotic integrity and biodiversity in a heavily exploited landscape, and also for their versatility as ecological surrogates, while exploring the other two topics provides applied information for wetlands management and conservation in protected area networks.

Study area
This study focussed mostly on a complex formed by 107 permanent inland wetlands at 18 N2000 sites in the Emilia-Romagna region, northern Italy ( Fig. 1, details available in Supplementary Material ESM1). These freshwater wetlands are located in a ~ 5600 km 2 heavily exploited, lowland agricultural landscape (within the range 4-40 m a.s.l.) through the Po River plain. This area is one of the most productive in the northern Hemisphere, and it is characterised by a marked imbalance between fertilisation and the uptake capacity of crops to result in widespread water contamination by phosphorus and nitrogen (Viaroli et al. 2018). The lowland Po plain presents a temperate climate (Köppen climate classification Cfa, Peel et al. 2007), with mean annual precipitation of ~ 700 mm and mean annual temperature of 13-14 °C.
Both habitats of European (EU) and national (ITA) relevance were considered to obtain the final overall HCI value (HCI codes and descriptions available in the results section, Table 1). Regarding SCC, attention was paid to the birds belonging to Annex I of the BD 1 3 All 107 wetlands were characterised to identify the role played by morphology, connectivity, land use and management (20 environmental variables; Table 2) to determine the patterns and richness of HCI. Data about HCI were obtained from the Emilia-Romagna Regional Habitat Map (updated in 2015), and from specific field surveys carried out in the summer and autumn of 2016 to verify the spatial patterns and local representativeness of HCI. Data on SCC were obtained from the Emilia-Romagna Region N2000 database (available at https ://ambie nte.regio ne.emili a-romag na.it/it/parch i-natur a2000 /rete-natur a-2000/habit at-e-speci e-di-inter esse-europ eo).

Statistical analysis
Firstly, the relationships among the abundance of national, European and total HCI were tested by Pearson correlations. Then, the potential concordance patterns of HCI and SCC were analysed in two complementary ways: richness and composition patterns.
Linear regressions and correlations were used to investigate the relationship between the overall richness patterns of habitats and taxa of conservation interest. The latter was possible by considering the general set of the N2000 sites (n = 18), which comprise 120 wetlands as data were available on a protected area scale (N2000, list and names available in ESM1), but not for every single water body.    Weighted classical multidimensional scaling (the "wcmdscale" function in vegan, also known as a weighted principal coordinates analysis) was used to summarise the composition patterns of HCI and SCC on the N2000 protected site scale (n = 18). Analyses were performed with a matrix based on presence-absence data, using the Jaccard Index (Legendre and Anderson 1999). A Procrustean analysis, accompanied by a permutation test (n = 9999), was applied to evaluate the degree and significance of the concordance between the ordinations obtained by respectively considering HCI and SCC. A Procrustean rotation analysis is regarded as a robust method for concordance analyses (Peres-Neto and Jackson 2001) and is frequently used to study ecosystem and community patterns (e.g. Virtanen et al. 2009;Guareschi et al. 2015b;Slimani et al. 2019). The statistic obtained is a Procrustes correlation r that derives from the symmetric Procrustes residual m 2 (r = √1 − m 2 ).
Having investigated the concordance patterns, the study focussed on the environmental variables associated with HCI assemblages and richness. Despite the biodiversity relevance of temporary wetlands and watercourses (e.g. Zacharias and Zamparas 2010), for the subsequent statistical analysis, attention was paid to the main pool of the 107 permanent systems at the 18 N2000 to manage a homogenous dataset. This allowed tests to be done with the complete list of environmental variables, including those related with water surface.
To avoid any collinearity among the 20 environmental drivers, a selection based on variance inflation factor (VIF with "th value" fixed = 2) was used to reduce the number of variables in the modelling procedure (Naimi et al. 2014).
Generalised linear models (GLM; Crawley 1993) were applied to investigate which environmental descriptors were associated with the richness of the HCI in the studied area. GLM analyses were carried out by considering a Poisson error distribution and a log-link function. Outliers, overdispersion and independence of residuals (Shapiro-Wilk test) were assessed following Fox and Weisberg (2011) and Zuur et al. (2009).
Non-metric multidimensional scaling (NMDS with function 'metaMDS') was used to assess the main patterns in the HCI assemblage structure. The analysis was done with a habitat matrix based on presence-absence data, using the Jaccard Index. The threshold at which ordination was not considered reliable was set at 0.20. Linear fittings, using the vegan function 'envfit', were performed between the selected environmental variables and the ordination outputs to identify the environmental factors associated with HCI distribution. The significance of the fitted variables was assessed by a permutation procedure (9999 permutations).
Finally, low richness areas (HCI ≤ 2) were explored to assess their role in habitat representativeness and were directly compared with richer areas (HCI > 2). All the analyses were performed with the R 3.5.1 software of the R statistical environment (R Core Team 2018) with packages "vegan" (Oksanen et al. 2018), "psych" (Revelle 2018), "usdm" (Naimi et al. 2014) and "car" (Fox and Weisberg 2011).
The habitat of national concern "Pa" (reed-beds, riparian formations of Phragmition) and EU habitat code 3270 (Rivers with muddy banks with vegetation of Chenopodion rubri pp and Bidention pp) were the commonest habitats, and were present in 99 and 87 permanent wetlands, respectively, and in 17 of the 18 N2000 protected sites (Fig. 2).
The correlation analysis stressed a very strong relationship (Pearson correlation r = 0.96, and Spearman rank correlation r = 0.97, both p < 0.01, n = 107, Fig. 3a) between EU habitats and the overall HCI value (HCI = EU + ITA). Lower values were obtained between the EU and national habitats (Spearman rank r = 0.63, p < 0.01). The same correlation values were obtained when considering the entire pool of 120 wetlands.
A linear model between species and habitats of conservation concern (HCI) stressed a strong relationship between their richness patterns at N2000 scale (adjusted R-squared = 0.693, p < 0.0001, n = 18, Fig. 3b). The same findings were obtained with the correlation analysis (Pearson r = 0.84, p < 0.001). This was also confirmed by testing the mean number of HCI when considering the 107 wetlands (adjusted R-squared = 0.384, p = 0.0036, n = 107).
According to the protest analysis, the assemblage of SCC and HCI was significantly concordant across the N2000 protected sites (r value = 0.758, m 2 statistic = 0.425, p value = 0.0001, n = 18). The strong correlation detected between SCC and HCI in both  Table 1) 1 3 terms of richness and distribution patterns (aim 1) allowed the multivariate and GLM analyses for HCI to be applied and to fully pursue the other research aims (2-3).

Richness and composition patterns of HCI and environmental variables
The GLM analysis explained 60.2% of HCI richness deviance (n = 107 wetlands, AIC = 421.7) with the considered set of environmental variables (Table 3). The water Fig. 3 a Correlation patterns (Spearman rank correlation and p < 0.01) among European (EUhab), national (ITAhab) and total habitats (HCI), upper panel; b Linear regression plot between HCI and species of conservation interest (SCC richness = 4.395HCI -2.086; lower panel) body mean perimeter (wb_per_mean), the percentage of wetland area occupied by aquatic areas (%_wb, with a negative sign), and the number of water bodies in each wetland (num_wb) were the significant environmental predictors associated with the HCI richness patterns (complete details in Table 3).
For HCI assemblage, the best two-dimensional solution ordination presented a final stress value of 0.126 (Fig. 4). Assemblage was principally affected again by the percentage of wetland area occupied by aquatic areas (%_wb) and by the hydrographic networks in terms of the distance and connectivity from other aquatic systems (dist_10 W, dist_riv). All the most significant environmental variables presented similar r 2 values (0.10-0.13; p < 0.01) (see Table 4).  Table 2 for details of variable codes and Table 4 for the complete results)

Comparing high and low HCI richness wetlands
Ninety of the 107 wetlands presented more than two HCIs (mean = 5.8, max = 9.0), while 17 wetlands were classified as poor in habitat richness terms (HCI ≤ 2). In both rich and poor wetlands, the national habitat "Pa" (Phragmition) and the EU habitat code 3270 (Chenopodion rubri pp and Bidention pp vegetation) were the most widespread in abundance and percentage terms, and "Pa" was particularly common in both cases (see Table 5 for details). EU habitats code 3140 (Hard oligo-mesotrophic waters with benthic vegetation of Chara spp.) and 91E0 (Alluvial forests with Alnus glutinosa and Fraxinus excelsior) were the rarest HCI as they appeared in only one (BB) and two (BO, and PT) protected sites, respectively. Accordingly, they were detected only in the richest group. The low richness group did not provide any exclusive or rare HCI (Table 5). Finally, even the few non-perennial wetlands excluded from the other statistical analysis did not provide any complementary or rare habitats (data not shown).

Conservation in a heavily exploited environmental matrix
According to previous studies, a strong concordance among taxonomic groups should be indicated by high r values (e.g. > 0.7, Heino 2010, and references within). Here, following our findings, the HCI patterns seemed to act as a good multipurpose proxy of threatened taxonomic species diversity (SCC) in richness and distribution terms, as we initially hypothesised. However, it should be taken into account that most of the SCC considered herein were waterbird species, whereas the other taxonomic groups were poorly represented. Based on this, further research is recommended to better understand and confirm the concordance patterns of HCI with other animal and/or plant species (e.g. multitaxon datasets), and ideally not only from a taxonomic point of view (multi perspective approach).  Table 5 Comparing low (HCI ≤ 2) and high (HCI > 2) habitat richness wetlands in terms of number, type and percentage of HCI (n. wetlands = 107) European and Italian habitats are displayed. The commonest habitats are shown in bold. See Table 1  The overall value of the 10 EU habitats of conservation interest (of a total value of 13 HCI) is not particularly high as it represents around 8% of all EU habitats that have been recognised for the Italian Continental Region (n = 83, Genovesi et al. 2014), but represents the bulk of the recorded HCI values (77%). Conservation in highly human-exploited ecosystems is a truly hard challenge for the different stakeholders involved, like conservationists, natural resource managers, environmental agencies, farmers, and society as a whole. The protected sites in similar matrices suffer multiple pressures, which imply the risk of biodiversity loss and trivial communities and, due to their geographical location, they can be greatly demanded by the general public (e.g. recreational services or improper uses). Current and future anthropogenic pressures have been demonstrated to have direct implications for protected areas and biodiversity conservation. In fact, the distance between protected areas and cities is predicted to dramatically reduce in some regions (McDonald et al. 2008). Moreover, the importance of the human footprint in shaping the distribution of numerous terrestrial and aquatic invaders have also been stressed (Gallardo et al. 2015), and heavily exploited environmental matrices like lowland agricultural and urban landscapes usually attain high levels of invasion (e.g. Chytrý et al. 2009). These issues represent challenges and threats to which conservation planning should respond (e.g. management and monitoring of reserves, Margules and Pressey 2000).

Richness and distribution HCI patterns
The richness model explained a large amount of deviance with the considered environmental predictors. It would appear that at least some variables were significantly associated with the HCI patterns in the studied sites.
On the one hand, high values of the mean perimeter of water bodies within a given wetland (wb_per_mean) and the presence of diverse waterbodies (in number terms) were positively associated with the highest HCI richness values. Both variables can indicate the extent of the riparian, and therefore ecotonal, zones in each wetland and, consequently, overall environmental (i.e. morphological) heterogeneity. The relevance of habitat heterogeneity has already been stressed by Shi et al. (2010) as being important for vascular plant species richness in Chinese wetlands, and by Báldi (2008) for invertebrate richness in Hungarian protected areas, but it is also a matter of some other contrasting findings (e.g. Palmer et al. 2010). Aquatic-terrestrial ecotones have been recognised as crucial zones because they act as natural filters, a source of organic carbon, and as connection lines for energy and materials fluxes (Décamps and Naiman 1990;Naiman et al. 2002).
On the other hand, the percentage of wetland area occupied by just aquatic areas was negatively associated with the HCI richness patterns. This should not clash with what has been discussed above. In fact, a massive extension of water surface can occupy a large wetland area without providing specific habitat and wide environmental heterogeneity. Accordingly, it would seem that numerous and diverse, but not extended, aquatic systems (i.e. water bodies) would benefit HCI richness by promoting heterogeneous riparian and ecotonal zones. This is indeed true for heavily exploited watersheds where the quality level of surface waters is expected to be generally poor. Similar evidence has been collected in a comparable system of wetlands located along the Oglio River (Lombardy region), a left tributary of the Po River (Bolpagni et al. 2013;Piotti 2015, 2016). These authors verified that aquatic and amphibian plant diversity was driven more by site features (natural vs. artificial), and by their morphological complexity (i.e. presence of wide ecotones) rather than by hydrology.

3
HCI distribution was associated with only a few significant predictors, which were not particularly strong in r 2 terms, but differed from those associated with HCI richness, except for the variable "percentage of wetland area occupied by aquatic areas". In this case, this predictor was mostly coupled with the descriptors related to the hydrographic network, such as distance to other aquatic systems (i.e. rivers and wetlands).
Overall, the traditional variables, widely considered in conservation planning like wetland area or wetland perimeter, were not the first stressed predictors in both cases (richness and distribution) despite some indirect relations with some significant variables being observed. This seems to underline the relevance of other variables (not just the commonest ones) associated with HCI patterns, at least in the considered highly exploited environmental matrix.
On the considered regional scale, the wetlands with low habitat richness seem quite trivial as they never bring new or complementary HCI. However, they still play a relevant role in habitat redundancy and resilience by ensuring the presence of habitats of interest in the study area.

Conservation implications and final remarks
In wetland biodiversity and monitoring context, habitats (in richness and distribution terms) can be considered a proxy for species of conservation interest (as a strictly taxonomic measure). Furthermore, HCI help to improve our capability to understand the real and/or potential conservation value of small natural features, like small wetlands in agricultural landscapes (Hunter Jr 2017;Bolpagni et al. 2019). Accordingly, they could be used as rapid tools for assessing the biotic integrity of wetlands and should also be tested under other geographical and environmental conditions. Significantly, only three invertebrate taxa are cited as SCC in the 18 studied protected sites, namely L. dispar (Lepidoptera), G. bilineatus (Coleoptera), and A. pallipes (Decapoda), along with one aquatic fern (M. quadrofolia, Marsileaceae). This is probably because there are only a few invertebrates (especially aquatic insects) and macrophytes, comprising non-vascular species, included on official Conservationist Lists (e.g. HD, Annex II). In fact, the European Union's conservation efforts are taxonomically biased towards vertebrates (Mammides 2019) which testifies that more research and efforts on non-target (e.g. invertebrates) and non-charismatic taxa is needed (Filz et al. 2013;Guareschi et al. 2015c;Habel et al. 2019). Specifically, periodic monitoring of these rare taxa is recommended (e.g. population persistence). Similarly, a recent ecological systematic review of the ecosystem value of small standing-water ecosystems, Bolpagni et al. (2019), verified the existence of clear narrative trends and wide knowledge gaps across geographical areas, biological components and target issues (e.g. competition, environmental drivers, human drivers).
Our research is an ambitious study into a lowland matrix exploited for agriculture and impacted by urbanization that provides useful information for its biodiversity conservation. Overall, the use of HCI as a rapid assessment tool may assist environmental managers and conservationists with limited resources. This is especially true considering that data on the presence of HCI are already available and freely accessible (at least considering the EU N2000 network, see Bunce et al. 2013). If not, HCI data seem easily achievable as their detection and effective monitoring during fieldwork is time-and cost-saving compared to in-depth taxonomic efforts and specific surveys for numerous taxa (e.g. Gigante et al. 2016). Integrating habitat-and species-based perspectives seems a promising field that 1 3 may provide a rapid assessment tool to acquire effective information for wetlands management and conservation.