A wide variety of taxonomic groups depend on the conservation of the wetlands that harbour them (Gibbs 2000, Quesnelle et al. 2013). The widespread deterioration of wetlands around the world has had a negative effect on the associated bird communities (Ma et al. 2010, Maclean et al. 2013) resulting in a decrease in species richness and abundance and is especially detrimental to scarcer species (Finlayson et al. 1992, Parsons et al. 2002, Paracuellos 2008). Among these birds are the marsh passerines, a polyphyletic group that, due to its high degree of specialization and dependence on the conservation of wetlands, has a large number of severely threatened species (Tanneberger et al. 2009, Trnka et al. 2014).
Marsh passerines have species-specific preferences for different wetland habitats, including the plant community and its spatial structure as well as the area of flooded surface (Báldi and Kisbenedek 1999, Poulin et al. 2002, Martínez-Vilalta et al. 2002). These interspecific differences mean that management techniques applied in wetlands may affect the abundance of each species in different and even opposite ways, which can potentially create conflicts over the conservation of several threatened species (Mitchell et al. 2006, Ma et al. 2010). For example, winter cutting of reed beds causes a decrease in the abundance of some species that do not find suitable nesting sites in the cut areas (Poulin and Lefebvre 2002), but in contrast, this technique does not seem to negatively affect, and may even benefit, the populations of other species (Vadász et al. 2008, Tanneberger et al. 2009). Actually, wetland management is rarely focused entirely on biodiversity conservation, which is often in conflict with other socio-economic factors (Birol et al., 2009, Downard et al. 2014) given the wide diversity of resources and ecosystem services provided by wetlands (Zedler and Kercher 2005, De Groot et al. 2012). Taking this issue into account, studies on habitat requirements of marsh birds are needed to achieve optimal wetland management and take appropriate conservation measures (Milsom et al. 2000, Poulin et al. 2002, Smart et al. 2006).
The National Park and Biosphere Reserve of Tablas de Daimiel (Spain) is one of the most important wetlands on the Iberian Peninsula for the conservation of many marsh passerines. Despite being protected, this wetland has suffered great disturbances in its hydrologic regime during the past decades as a result of human activities (Álvarez-Cobelas et al. 2001, Sánchez-Carrillo et al. 2016). This has led to artificial water management in the park to ensure there is a minimum flooded area. Water transfers, the creation of three dams, and groundwater pumping from emergency wells are some of the measures taken for this purpose, although they appear to have limited effectiveness (Sánchez-Carrillo and Álvarez-Cobelas 2010, Aguilera et al. 2013). The overexploitation of the aquifer for agriculture continues to cause an alarming reduction in the wetland hydroperiod and the area of flooded surface (Llamas 1988, Álvarez-Cobelas et al. 2001), reaching critical situations such as the one that occurred in 2009 when the aquifer dried out completely and ignited by self-combustion (Moreno et al. 2010, Aguilera et al. 2013). These large fluctuations in water level can have negative impacts on the populations of threatened marsh passerines that breed in this area (Desgranges et al. 2006, Sánchez-Carrillo and Angeler 2010). In addition, variations in the water regime have caused significant changes in its plant community (Álvarez-Cobelas et al. 2001, Sánchez-Carrillo et al. 2016) resulting in a shift in the availability of suitable habitats for passerines (Poulin et al. 2002, Báldi 2005, Jiménez et al. 2015) that affects the abundance and potential recovery of populations of endangered species.
We analyzed the effect of water regime variation in the protected wetland of Tablas de Daimiel over the past decade on threatened marsh passerines to understand how water management can influence their conservation. Our hypotheses, based on previous studies, suggest that (1) problems in the water regime of the wetland may affect the population trends of the marsh passerines breeding in the national park. Although we expect (2) differences in the population trends of the study species due to (3) the differences of species-specific habitat requirements, and (4) declining water levels result in the disappearance of suitable habitats for some species and their populations will consequently be affected in a negative way.
The National Park of Tablas de Daimiel (Ciudad Real, Spain, 39º 08' 36'' N, 03º 42' 18'' W, 606 m a.s.l.), is one of the largest continental wetlands in the Iberian Peninsula (up to 1800 ha of flooded surface) and part of the Mancha Húmeda Biosphere Reserve (Fig. 1). The water regime of the wetland depends on the confluence of the Cigüela River with the discharge zone of the La Mancha Occidental aquifer (5500 km²), located in the Guadiana River Basin. Historically, the flooded area was maintained mainly by the discharge of the aquifer, but due to its over-exploitation, the wetland may at present even act as a recharge zone (Aguilera et al. 2013). This may lead to hydrological crises such as the one that occurred in 2006 when there was a severe drought eventually leading to the total desiccation of the wetland in 2009. In 2010, however, the park recovered its entire flooded area due to the discharge of the aquifer caused by the significant rise of the phreatic level that occurred after a large amount of rain during that year.
During the period of study (2012-2019) the flooded area of the wetland fluctuated strongly depending on the amount of groundwater reaching it from the Guadiana river basin. The general trend in the size of flooded area over the years was negative: y = -110.48 ± 223996.01, t = -3.15, R² = 0.56, P = 0.02 (Linear Regression model was conducted from the total flooded area data, in hectares, of the national park on the first weekend of May each year). Although in the first years the flooded surface was near its maximum, the flooded surface drastically declined throughout the last five years of the study and in the last year (2019) only 675 hectares were flooded, which is equivalent to a third of the total capacity of the wetland.
The wetland plant community is characterized by large areas of helophytic flora. The common reed (Phragmites australis) is the dominant species and it usually covers large monospecific areas. There are also patches of other species such as bulrush (Typha domingensis) or saw sedge (Cladium mariscus). In the driest areas, where the surface is not permanently flooded, grasslands predominate and the tamarisks (Tamarix canariensis and Tamarix gallica) constitute the arboreal strata.
The target species of this study were the Great Reed Warbler (Acrocephalus arundinaceus), the Moustached Warbler (Acrocephalus melanopogon), the Reed Bunting (Emberiza schoeniclus witherbyi), the Savi's Warbler (Locustella luscinioides), and the Bearded Tit (Panurus biarmicus). The criteria used for the selection of these species were that (1) the species were included in the Spanish Catalogue of Threatened Species and other national inventories of threatened fauna, (2) their population size in this wetland was important for the viability of this species in a regional and/or national context (e.g. Madroño et al. 2004, Monrós et al. 2018), and (3) that there was evidence of differences in the habitat requirements of the selected species (e.g. Martínez-Vilalta et al. 2002).
With the aim of estimating the population size and habitat selection preferences of the focal species, listening point censuses (Bibby et al. 1992) were carried out during the first weeks of May (based on the reproductive phenology of these birds) for 8 consecutive years, starting in 2012. Each year more than a hundred listening points were established over two consecutive days. We restricted sampling to the first and last three hours of the day when the peak of activity occurs and therefore it is easier to detect these species.
Prior to the censuses, the wetland was divided into several routes, which would be repeated every year, in order to cover the maximum surface area and the widest representation of habitats possible (Fig. 1). Five-minute listening points were established every 300 m of the route, and once each point was completed, the number of individuals of each target species detected was recorded. The number of listening points varied slightly between years due to issues related to the wetland conditions (such as changes in accessibility to certain areas due to changes in the flooded surface or the cover of marsh vegetation). We checked that such small variation produced no major effects on the number of detected individuals in other years.
Special attention was given to the census of the Reed Bunting. During the 5 minute sampling period, a recording of the territorial song of a male was played twice for one minute, followed by one minute of silence, resulting in the sequence: silence-playback-silence-playback-silence. This methodology was used according to the methodology of the Spanish national census for the species (Atienza 2006), because this species requires such a specific protocol (which does not affect the detection of the other species) in order to obtain a precise population estimate, which is of vital importance due to the critical situation that the taxon is going through.
Imperfect detection was taken into account when creating the census design. Target species detection was maximized by carrying out the censuses in the first and last 3 hours of the day in early May, in accordance with suggestions in the detection probability study carried out by Jimenez et al (2015) in the Tablas de Daimiel National Park with 3 of the 5 species in the study. These conditions were the reason for concentrating the greater part of the sampling effort (~130 listening points in an area of approximately 1800 hectares) in only two consecutive days. Census methodology for the Reed Bunting was also different from the others as it seems that this species is mostly detected visually, unlike other species that are easily detected by their song (Moskát and Báldi 1999, Jimenez et al 2015).
In addition, a characterization of the habitat was made at each listening point. To do this, we recorded the percent cover of each of six categories within a radius of 25 meters around the listening point. Four categories referred to vegetation type, including three dominated by the species Cladium mariscus, Phragmites australis, and Typha domingensis and a non-specific category called "DrierVeg", which included vegetation characteristic of non-floodable zones within the wetland, such as Tamarix sp. or grasses. The remaining two categories corresponded to either “OpenWater" or "DrySurface" within the floodable zone (i.e. bare mud).
Population trends of the different species were determined with Generalised Additive Models (GAMs) using the abundance detected of each species at each listening point in the different study years (Fewster et al. 2000). Models were constructed with the “gam” function of the mgcv package (Wood 2016). GAM models have previously been used to assess the long-term population trends of a wide range of bird groups, including waterbirds (Fewster et al 2000, Clarke et al 2003, Atkinson et al 2006).
To identify habitat requirements and to check which variables may play a stronger effect on the abundance of each species, we used Generalized Linear Mixed Models (GLMM), with a negative binomial distribution of the response variable and a log-link function. The cover of Cladium mariscus (saw sedge), Phragmites australis (common reed), Typha domingensis (bulrush), “DrierVeg”, “OpenWater” and “DrySurface” were used as fixed effect variables in the models, as they were the main habitat elements found in the wetland and have been used in previous similar studies (Martínez-Vilalta et al. 2002, Poulin et al. 2002, Vera et al. 2011). The total flooded area (“TFloodedArea”) of the wetland of the corresponding sampling year, was incorporated into the models as a random variable. Collinearity was tested prior to analysis, but no variable was discarded since no high degree of correlation was detected between the studied variables (R<0.7) (Dormann et al. 2013). These models were constructed in R 3.5.2 (R Core Team 2018) with the function “glmer.nb” of the lme4 package (Bates et al. 2015).
The best models resulting from the combination of variables were automatically selected according to the Akaike Information Criterion (AIC) using the “dredge” function of the MuMIn package (Barton 2020). For each species, the best model and those that did not exceed the lowest AIC value by two units were considered, and the weight (w) value was taken into account to calculate the relative importance of each of the variables (RIV) in the set of selected models (Burnham and Anderson, 2002). Finally, the function “tab_models” of the sjPlot package (Lüdecke et al. 2020) was used to obtain the descriptors for each of the selected models, including the RM² and RC² values to show the percentage of the model variance that is explained by the fixed effect variables and the random effect variables (Bates et al. 2015, Mestre et al. 2019).
To obtain an overall view of the relationships between the studied marsh passerines and between the available habitat types through time, we used the constrained ordination technique of Redundancy Analysis (RDA), which allows relating species assemblages to environmental gradients (Legendre and Legendre 2012). RDA was performed with the software Canoco5 (ter Braak and Šmilauer 2012).
The surveyed areas of the National Park of Tablas de Daimiel were mostly dominated by open water (Fig. 2). Open water was greatest in 2014, when it represented around 50% of the sampled wetland area, but declined as the years passed, along with the decrease in the total flooded area. Decline in water levels caused the appearance of dry surfaces within the usually flooded area (≈18% of the surface) during the last year.
As for vegetation, common reed was the dominant species in the wetland. It showed maximum extension values of around 40% of the sampled area at the beginning of the study period and decreased slightly in recent years. The area covered by bulrush was small at the beginning of the study period (≈10%) but increased over the years, although with large fluctuations, at times occupying up to 25% of the sampled area. The amount of saw sedge remained relatively constant at cover values close to 5% during the study period. Finally, the areas occupied by vegetation of non-flooded zones increased considerably in recent years after the initial drop at the beginning of the study.
The studied species experienced distinct and even opposite population trends during the period from 2012 to 2019 (Fig. 3).
The Great Reed Warbler was the only species with a substantial population size at the beginning of the study period. Its population peaked in 2013 with an average of 1.56 individuals per point count (ind./point) detected. The species experienced a severe population decline of 90% in 2019 when the mean number of Great Reed Warblers per point was 0.14. In contrast, the Reed Bunting started from a very low number of individuals (0.10 ind./point) and its populations progressively increased until 2019, quadrupling its population in 8 years (0.39 ind./point).
The Moustached Warbler and the Savi’s Warbler had similar population trends throughout the study period. Starting from a low number of individuals, especially for the Moustached Warbler (0.03 ind./point), both species experienced a considerable population increase until 2014, when the population growth stopped. Between 2014 and 2016 there was a two-year period of stability in the Savi’s Warbler populations and of population decline in the Moustached Warbler, preceding a new period of increase until 2018. In that year there is again a change in the population trend and the populations of both species fall considerably (from 0.40 to 0.10 ind./point in the Moustached Warbler and from 0.94 to 0.79 ind./point in Savi’s Warbler).
Finally, the Bearded Tit had strong population growth until 2017 (from 0.23 to 2.06 ind./point) when the trend changed drastically and it suffered a sharp decline that reduced the number of individuals by more than 50% (0.86 ind./point in 2019).
The results obtained indicate that the Moustached Warbler, Savi’s Warbler and Bearded Tit avoided vegetation from non-floodable areas within the wetland such as grasslands and tamarisk formations (RIV=1). This habitat element, however, does not seem to have had a major influence on the presence of Reed Bunting (RIV=0.5) and Great Reed Warbler (RIV=0.38) populations (Table A1.1, Fig. 4).
The Reed Bunting was underrepresented in habitat dominated by large bulrush coverage but benefited from the presence of saw sedge formations (Table A1.1, Fig. 4). High reed cover, on the other hand, seemed to negatively affect the Moustached Warbler (Table A1.1, Fig. 4, Fig. 5). For this species, the proportion of saw sedge was also generally higher in used than in available habitat, but this relationship was weaker than in the previous case, with a RIV= 0.55 in the GLMMs. Saw sedge, in addition, also seems to have a positive effect on Savi’s Warbler numbers (Fig. 4, Table A1.1).
The Bearded Tit, besides avoiding areas with non-floodable vegetation, also avoided the dry areas of the floodable zone of the wetland resulting from the decrease in water surface during the last study year. These dry surfaces, on the other hand, appeared to be used by the Reed Bunting (Table A1.1, Fig. 4).
The amount of open water was an important habitat element for every species' abundance, generally affecting positively except for Moustached Warbler, which apparently avoided this kind of habitat.
For the Great Reed Warbler the RM² value represented between 4 and 11% of the RC² value obtained in the best GLMMs (Table A1.1), so the fixed effect variables considered explained a very small proportion of the variance of this species compared to the random variable (total flooded area). In contrast, for the Moustached Warbler, Savi’s Warbler and Bearded Tits, RM² accounted for 70-80% of the total variance explained by the models, suggesting a greater strength of the fixed effect variables on the abundance of these species. For the Reed Bunting, with RM² representing around 55% of RC², the percentage of variance explained by the fixed and random effects variables was more balanced.
The RDA results (Fig. 5) indicate a strong relationship between the Great Reed Warbler population variability and total flooded area of the wetland, so as the negative effects of time progression (Year) on the species populations. The other species were all ordered on the opposite part of axis 1. In addition, the Reed Bunting was ordered on the positive part of axis 2 since its populations increased with time and benefited from the presence of saw sedge and dry surfaces in the floodable zone. The Bearded Tit, on the other hand, was ordered on the extreme negative part of axis 2, due to the strong negative effect that dry areas and non-floodable zones had on its populations, and the positive effect of open water surfaces on this species, in agreement with the GLMM results. Finally, Savi’s and Moustached Warbler remained in the intermediate part of axis 2, apparently negatively affected by the presence of non-floodable vegetation type and large reed extensions (although this was not entirely corroborated by the other analysis).
Our results, in agreement with similar studies carried out in other locations, show that the different species that form the marsh passerine community of Tablas de Daimiel have species-specific habitat requirements (Poulin et al. 2002, Martínez-Vilalta et al. 2002, Trnka and Prokop 2006, Jiménez et al. 2015).
The total flooded surface area of the wetland was the most important variable explaining the abundance of the Great Reed Warbler. Widely flooded areas, especially with certain fluctuation of water levels, promote insect population explosions and the development of reeds with optimal stem thickness that allow the Great Reed Warbler to perch and build its nests (Jedraszko-Dabrowska 1992; Graveland 1998, Poulin et al. 2002). More flooded area in the wetland leads to a greater amount of available open water which was positively related to Great Reed Warbler abundances, according to the GLMM results. This is in agreement with previous studies showing that this species is present at higher densities near the vegetation edges, especially in reed bed formations close to free water, where it also builds its nests (Rolando and Palestrini 1989, Graveland 1998, Prokešová and Kocian 2004, Trnka et al. 2009). Despite the large number of studies that conclude that there is a close relationship between this species and reeds (Rolando and Palestrini 1989, Martínez-Villalta et al. 2002, Méro et al. 2015), our results do not show a special preference by this species for reed over open water flooded areas. This may be indicating that the condition of reed plants, such as thickness, height or stem density, which depend upon flooding levels (Graveland 1998), might be more important for the Great Reed Warbler than the amount of reed cover itself. Finally, despite the strong dependence on the flooded surface found, our results indicate that Great Reed Warbler did not avoid dry vegetation from non-floodable areas within the wetland, as other species did. This may be because a large part of the wetland population of this species can be found in the south-western part of the national park, where the flooded surface and aquatic vegetation are very close to the boundary roads where grasslands and tamarisks are common, and where the Great Reed Warblers can occasionally be detected.
Our results show a clear positive relationship between the Reed Bunting and saw sedge formations, in agreement with previous surveys (Martínez-Villalta et al. 2002, Bigas and Copete 2004, Vera et al. 2009, Jiménez et al. 2015). The presence of dry, unvegetated surfaces in the wetland floodplain also favored its populations. This seems to be related to the way the Reed Bunting feeds and moves; in addition to climbing through the marsh vegetation, it frequently moves on the ground (Martínez-Villalta et al. 2002). Vera et al. (2011, 2014) found the presence of bulrush patches, rather than saw sedge, was a determining factor in the distribution of this Reed Bunting subspecies in Spanish wetlands (Vera et al. 2011, 2014). However, we found they seemed to avoid that vegetation type. These contradictory results could be explained by the difference in the scale of analysis of both studies. Vera et al.'s (2011) study was at a wetland scale within the distribution area of the species while our research was conducted at a microhabitat scale within the wetland. Perhaps the presence of bulrush patches is indicative of some feature of the wetland that favors the presence of the species in that particular locality, but on a smaller scale, Reed Buntings preferentially use saw sedge instead of bulrush formations.
Previous studies on the habitat requirements of Moustached Warbler’s Spanish populations, including Tablas de Daimiel, indicate that the species avoids large monospecific reed extensions (Martínez-Villalta et al. 2002, Castany and López 2005) and favors other helophytes that are more sensitive to fluctuations in the water level such as saw sedge (Jiménez et al., 2015). Our results support these conclusions, indicating a negative effect of high reed cover on the presence of the species and a positive effect, though weaker than in the previous case, of the saw sedge formations. Due to their nesting and feeding habits, it has been found that this marsh passerine inhabits areas close to flooded surfaces (Bibby 1982, Fessl and Hoi 1996, Trnka and Prokop 2006), in agreement with the negative effects of dry vegetation on its abundance in our study. The GLMM results, however, indicated that open water surfaces negatively affected the abundance of the species. The Moustached Warbler builds its nests at the base of the helophytic vegetation (0-0.5m), much lower than other species such as the Great Reed Warbler (Leisler 1991a, b). The large fluctuations in the water levels of the wetland during the study period, which could facilitate predation or cause flooding of the nests, could be reasons why this species avoided the low areas near wide open water surfaces and perhaps preferred small, more sheltered, waterlogged areas.
The Savi’s Warbler, like the Reed Bunting and Moustached Warbler, showed a preference for habitats with high cover values of saw sedge. The same results were found by Martínez-Villalta et al. (2002) in the Ebro Delta (Spain). This preference, however, varies across its geographical distribution, suggesting that the structure formed by vegetation, specifically the presence of an extensive basal stratum of straw litter, is more important than the plant species itself (Neto 2006, Bergner and Gezelius 2013). This preference has been suggested to be related to its mode of foraging near water surfaces (Bibby and Lunn 1982, Bergner and Gezelius 2013), which in turn agrees with our results on the positive effects of open water on the abundance of this species. The preference of the species for areas with saw sedge, a vegetation type of flooded areas with reduced fluctuations, and for open water areas, is supported by the negative effects of dry terrestrial vegetation on the abundance of the species found in our study. In contrast, other authors conclude that the species has a preference for dry areas where it also feeds (Martínez-Villalta et al 2002, Trnka and Prokop 2006).
Bearded Tits were absent in the driest areas of the wetland, where bare mud and vegetation of non-floodable surfaces appeared, while their abundances increased in open flooded areas. Dependence on open water by the Bearded Tit has been previously recorded and explained by its feeding and nesting mode (Bibby and Lunn 1982, Poulin et al. 2002, Brichetti and Grattini 2008). Despite this apparent dependence of the species on flooded areas, no relation was found between the total flooded area of the park and the abundance of the Bearded Tit. This may be the result of a temporary mismatch derived from the small initial population size, caused by the water crisis in the previous years and the sudden recovery of the total flooded area of the wetland in only one year. Adverse climatic conditions during particularly severe winters, such as heavy flooding, can have a significant effect on breeding populations of this species (Campbell et al 1996, Wilson and Peach 2006). Thus, the large floods recorded in the wetland in late 2010, together with other factors such as the absence of suitable habitat in previous years, could be responsible for the Bearded Tit’s small initial population size. Furthermore, no vegetation type was selected by the species in the study period, in agreement with the study by Jiménez et al (2015) carried out in Tablas de Daimiel, where they suggested that this species may be more generalist than other sympatric species such as the Reed Bunting and the Moustached Warbler.
Due to the conservation problems associated with the hydrology of this wetland, between 2006 and 2010 the Tablas de Daimiel National Park was completely dry. Under these conditions, the aquatic vegetation dependent on regular water levels was substantially decimated compared to the reed beds. Open water areas were non-existent and terrestrial vegetation began to replace the helophytic vegetation due to the absence of floods (Sánchez-Carrillo et al 2016). These transformations in the available habitat may have been responsible for the small population size of most species at the beginning of the study period according to our results, and as has been observed in other marsh passerine species in the same wetland (Jiménez et al. 2018).
After a particularly rainy autumn in 2010 and with help from artificial water management, the wetland suddenly recovered its entire flooded surface. The great fluctuations in the water levels caused the reed to benefit over other species of helophytes such as the saw sedge (Sánchez-Carrillo et al 2016) and also to have a considerable stem thickness (Poulin et al 2002). The sudden inundation of the wetland may have had negative effects on the abundance of some marsh passerines as mentioned above, but on the other hand, it also led to the appearance of wide-open water surfaces and to abundance peaks in insects. According to our results and the existing literature (Graveland 1998, Poulin et al. 2002, Prokešová and Kocian 2004), such conditions might have been optimal for the Great Reed Warbler’s populations, which at the beginning of the study period presented maximum values of abundance.
The presence of areas with constant water levels led to the appearance of extensions of other helophyte vegetation such as saw sedge and bulrush, which replaced reeds under certain conditions (Álvarez-Cobelas et al. 2001, Jiménez et al. 2015). These changes in the wetland vegetation structure may have benefited the populations of marsh passerines associated with these formations, including Reed Bunting, Moustached Warbler and Savi’s Warbler (Martínez-Villalta et al. 2002, Bigas and Copete 2004, Jiménez et al. 2015), which started to become more abundant. Due to the re-flooding, the dry areas and terrestrial vegetation that occupied the inner surface of the wetland disappeared and were occupied by open water, which would in turn, benefit the Bearded Tit populations. All these endangered birds appear to have benefited from the re-establishment of the wetland’s water regime and the changes in habitat structure resulting from it. However, as the years went by, the flooded area began to decrease again.
The strong reduction of the flooded surface seems to be the main cause of the sharp population decline suffered by the Great Reed Warbler. Even though some areas remained flooded, the decrease in water levels, especially pronounced between 2017 and 2019, caused the disappearance of wide open water surfaces, the appearance of dry areas, and an increase in terrestrial plant cover within the wetland. Those factors affected the abundance of the Bearded Tit, Moustached Warbler, and Savi’s Warbler and can explain population declines observed in the last study years. The population of the Reed Bunting, however, was most likely sustained by saw sedge plasticity, which may continue to grow even when the surface begins to dry out (Sánchez-Carrillo et al 2016) and by the appearance of the dry areas that seem to be used by the species for feeding (Martínez-Villalta et al. 2002).
We took precautions regarding bird detectability through our monitoring design, and after evaluating the data we consider these precautions sufficient to analyze population trends using the GAM models. Although problems with imperfect detectability can be addressed with statistical analyses that correct for consequential biases (Kéry and Schmidt 2008) to better estimate population abundances of the studied species, to focus conservation efforts on those having reduced densities, further research should take detectability issues into consideration.
Despite the numerous management strategies applied, the problems associated with the wetland water regime in the Tablas de Daimiel National Park remain unsolved (Sánchez-Carrillo et al 2021). During the study period, the flooded surface of the wetland went from maximum capacity to only a third in 5 years. The consequences on the marsh passerine community were recorded in the later years of the study. If the water level problems are not fixed, the park will dry out entirely again, which will lead to the collapse of the populations of the studied endangered species and many more.
Tablas de Daimiel is one of the most important wetlands at a regional and national level for the conservation of marsh passerines. For instance, in the first national census of the endangered Reed Bunting in 2005, its population in the park was the second largest in Spain, and it has been considered important for the viability of the species (Kvist et al. 2011). However, 10 years after the survey the number of breeding pairs had been reduced by almost 90%. The population increase recorded in the last years of this study is an important change of trend for the conservation of this species (Atienza 2006, Monrós et al. 2018) but if the current situation persists, its populations will fall again due to the disappearance of its suitable habitats (Sánchez-Carrillo et al 2016). Las Tablas de Daimiel is also the most important wetland for the metapopulation of the Bearded Tit in Castilla la Mancha, hosting the main population of this marsh passerine in the Iberian Peninsula due to the small population sizes in the northern wetlands of the region and the dramatic decline suffered in the Mediterranean ones (López-Iborra and Monrós 2004, Belenguer Barrionuevo et al. 2016). Reliable management strategies that would guarantee the viability of the populations in the park are required to ensure the conservation of these species.
The results of the present research highlight the need to ensure the long-term viability of the populations of threatened marsh passerines. Water management should focus on ensuring a minimum flooded area that would allow the development of helophyte vegetation, the availability of open water, and also prevent terrestrial vegetation from spreading over the wetland. The growth of terrestrial vegetation in periods of drought seems to be detrimental to several species in this study, and management by grazing, mowing, and burning seems to be effective solutions to this kind of issue (Vulink et al. 2000, Middleton et al. 2006). In addition, it is recommended that avoiding an excess of stability in the water level to recreate natural fluctuations that allow for explosions of insect populations, the appearance of dry zones, and changes in the composition and structure of the marsh vegetation, will benefit the populations of the passerine species (Poulin et al. 2002, Timmermans et al. 2008, Beemster et al 2010).
When making decisions about which management strategies to adopt, it should be kept in mind that potential conflicts between the conservation of coexisting endangered species may appear as a result of conservation actions focused on individual species. Many factors, such as interactions between species or differences in their habitat requirements, may cause certain management techniques to affect species in the same system in different or even opposite ways (Morrison et al. 1996, Thirgood et al. 2000, Raimondi et al. 2015). For example, previous works have found potential conflicts between the conservation of some marsh passerines, such as Bearded Tits and other waterbirds, resulting from water management focused on individual species (Baldi 2005, Wilson and Peach 2005).
Potential conflicts derived from differences in habitat requirements can be identified between our study species. For instance, dry areas on the floodable surface had a positive effect on the abundance of Reed Bunting but a negative effect on the Bearded Tit. As previously said, the Great Reed Warbler is strongly related to reed formations (Martínez-Villalta et al. 2002, Méro et al. 2015), but some species were dependent on the development of other helophyte plants such as saw sedge and even avoided monospecific reed extensions. Management focused on creating heterogeneous wetlands could be the solution to these possible conservation conflicts derived from species-specific habitat requirements. In fact, the heterogeneity of these types of environments appears to be a relevant cause of the richness and abundance of marshland birds (Beemster et al 2010). This demonstrates the importance of conducting studies of this type to carry out integrated management actions and ultimately ensure that the conservation measures applied are as effective as possible (Poulin et al. 2002).
Marsh passerines have species-specific habitat requirements, which means that changes in the wetland system may affect them in different and even opposite ways. Because of the complexity of biological communities, integrative studies of communities are necessary to carry out appropriate management actions and avoid indirect effects of conservation measures focused on a single species.
Unsolved problems related to the wetland water regime in the Tablas de Daimiel National Park threaten its community of endangered birds. Reduction of aquifer overexploitation and restoration of the natural water regime should be considered a priority to guarantee the conservation of these species that, in spite of having specific habitat requirements, ultimately depend on the wetland hydrology. On the other hand, management focused on promoting habitat heterogeneity within the same wetland may be the best option to address species-specific differences in habitat selection processes and thus minimize potential conflicts.
ACKNOWLEDGMENTS
We thank the Tablas de Daimiel National Park for providing the necessary authorization and logistics to carry out the censuses. This work was supported by the Fundación Biodiversidad of the Ministerio para la Transición Ecológica y el Reto Demográfico of Spain and by a FPU (Formación de Profesorado Universitario) grant (FPU/02579) from the Ministerio de Ciencia, Innovación y Universidades of Spain to I.A. Finally, we would like to thank the national park staff and the external volunteers who participated in the censuses during the years of the study: Edu, Carlos, Rubén, Marcial, Pablo, Pau, Lydia, Joaquín, Claudio, Felipe, Escudero, Dani, Manolo, Luís, Crescencio, Borondo, Álex, Manuel, Jordi, Alejandro, Chulen, Santi, Jesús, Pedro, Cruces, Consuegra.
Aguilera, H., S. Castaño, L. Moreno, M. E. Jiménez-Hernández, and A. De la Losa. 2013. Model of hydrological behaviour of the anthropized semiarid wetland of Las Tablas de Daimiel National Park (Spain) based on surface water-groundwater interactions. Hydrogeology Journal 21:623-641. https://doi.org/10.1007/s10040-012-0950-3
Álvarez-Cobelas, M., S. Cirujano, and S. Sánchez-Carrillo. 2001. Hydrological and botanical man-made changes in the Spanish wetland of Las Tablas de Daimiel. Biological Conservation 97:89-98. https://doi.org/10.1016/S0006-3207(00)00102-6
Atienza, J. C. 2006. El escribano palustre en España. I Censo Nacional. 2005. SEO/BirdLife, Madrid, Spain.
Atkinson, P. W., G. E. Austin, M. M. Rehfisch, H. Baker, P. Cranswick, M. Kershaw, J. Robinson, R. H.W. Langston, D. A. Stroud, C. Van Turnhout, and I. M.D. Maclean. 2006. Identifying declines in waterbirds: The effects of missing data, population variability and count period on the interpretation of long-term survey data. Biological Conservation 130(4): 549-559. https://doi.org/10.1016/j.biocon.2006.01.018
Báldi, Á. 2005. The importance of temporal dynamics of edge effect in reedbed design: a 12-year study on five bird species. Wetlands Ecology and Management 13:183-189. https://doi.org/10.1007/s11273-004-0313-5
Báldi, Á, and T. Kisbenedek. 1999. Species-specific distribution of reed-nesting passerine birds across reed-bed edges: effects of spatial scale and edge type. Acta Zoologica Academiae Scientiarum Hungaricae 45(2):97-114.
Barton, K. 2020. Mu-MIn: Multi-model inference. R Package Version 1.43.17. https://cran.r-project.org/web/packages/MuMIn
Bates, D., M. Maechler, B. Bolker, and S. Walker. 2015. Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67(1):1-48. https://doi.org/10.18637/jss.v067.i01
Beemster, N., E. Troost, and M. Platteeuw. 2010. Early successional stages of Reed Phragmites australis vegetations and its importance for the Bearded Reedling Panurus biarmicus in Oostvaardersplassen, The Netherlands. Ardea 98(3):339-354. https://doi.org/10.5253/078.098.0308
Belenguer Barrionuevo, R., G. M. López-Iborra, J.I. Dies, and J. Castany i Alvaro. 2016. Dramatic decline of the bearded reedling, Panurus biarmicus, in Spanish Mediterranean wetlands. Animal Biodiversity and Conservation 39(1):17-27. https://doi.org/10.32800/abc.2016.39.0017
Bergner, A. and L. Gezelius. 2013. Territory vegetation structure and habitat preferences of the Savi’s Warbler Locustella luscinioides in lake Tåkern, Sweden. Ornis Svecica 23:72-80. https://doi.org/10.34080/os.v23.22575
Bibby, C. J. 1982. Studies of West Paleartic birds: 184 Moustached Warbler. Britanic Birds 75(8):346-359.
Bibby, C. J. and J. Lunn. 1982. Conservation of reed beds and their avifauna in England and Wales. Biological Conservation 23:167-186. https://doi.org/10.1016/0006-3207(82)90074-X
Bibby, C. J., N. D. Burgess, D. A. Hill, and S. H. Mustoe. 1992. Bird census techniques. Academic Press, London, England.
Bigas, D. and J. L. Copete, J. 2004. Repicatalons Emberiza schoeniclus. Pages 546-547 in J. Estrada, V. Pedrocchi, L. Brotons, and S. Herrando, editors. Atles dels ocells nidificants de Catalunya 1999-2002. Lynx Edicions, Barcelona, Spain.
Birol, E., N. Hanley, P. Koundouri, and Y. Kountouris. 2009. Optimal management of wetlands: Quantifying trade-offs between flood risks, recreation, and biodiversity conservation. Water Resources Research 45:w11426. https://doi.org/10.1029/2008WR006955
Brichetti, P. and N. Grattini. 2008. Distribuzione, consistenza ed evoluzione delle popolazioni di basettino Panurus biarmicus nidificanti in Italia nel periodo 1980-2006. Avocetta 32:47-53.
Burnham, K. and D. Anderson. 2002. Models selection and multi-model inference. Springer, Berlin, Germany.
Campbell, L., J. Cayford, and D. Pearson. 1996. Bearded Tits in Britain and Ireland. British Birds 89:335-346.
Castany, J. and G. López. 2006. El carricerín real en España. I Censo Nacional (2005). SEO/BirdLife, Madrid.
Clarke, E. D., L. B. Spear, M. L. McCracken, F. F. C. Marques, D. L. Borchers, S. T. Buckland, and D. G. Ainley. 2003. Validating the use of generalized additive models and at-sea surveys to estimate size and temporal trends of seabird populations. Journal of Applied Ecology 40:278-292. https://doi.org/10.1046/j.1365-2664.2003.00802.x
De Groot, R., L. Brander, S. van der Ploeg, R. Costanza, F. Bernard, L. Braat, M. Christie, N. Crossman, A. Ghermandi, L. Hein, S. Hussain, P. Kumar, A. McVittie, R. Portela, L. C. Rodriguez, P. ten Brink, and P. van Beukering. 2012. Global estimates of the value of ecosystems and their services in monetary units. Ecosystem Services 1:50-61. https://doi.org/10.1016/j.ecoser.2012.07.005
Desgranges, J-L., J. Ingram, B. Drolet, J. Morin, C. Savage, and D. Borcard. 2006. Modelling wetland bird response to water level changes in the Lake Ontario - St. Lawrence River hydrosystem. Environmental Monitoring and Assessment 113:329-365. https://doi.org/10.1007/s10661-005-9087-3
Dormann, C. F., J. Elith, S. Bacher, C. Buchmann, G. Carl, G. Carré, J. R. García Marquéz, B. Gruber, B. Lafourcade, P. J. Leitão, T. Münkemüller, C. McClean, P. E. Osborne, B. Reineking, B. Schröder, A. K. Skidmore, D. Zurell, and S. Lautenbach. 2013. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 36:27-46. http://doi.org/10.1111/j.1600-0587.2012.07348.x
Downard, R., J. Endter-Wada, and K. M. Kettenring. 2014. Adaptive wetland management in an uncertain, and changing arid environment. Ecology, and Society 19(2):23. http://doi.org/10.5751/ES-06412-190223
Fessl, B., and H. Hoi. 1996. The significance of a two part song in the moustached warbler (Acrocephalus melanopogon). Ethology Ecology & Evolution 8:265-278
Fewster, R., S. Buckland, G. Siriwardena, S. Baillie and J. Wilson. 2000. Analysis of population trends for farmland birds using generalized additive models. Ecology 81(7):1970-1984. https://doi.org/10.2307/177286
Finlayson, C. M., G. E. Hollis, and T. J. Davis, editors. 1992. Managing Mediterranean wetlands and their birds. IWRB, Slimbridge. England.
Gibbs, J. P. 2000. Wetland loss and biodiversity conservation. Conservation Biology 14:314-317. https://doi.org/10.1046/j.1523-1739.2000.98608.x
Graveland, J. 1998. Reed die-back, water level management and the decline of the great reed warbler Acrocephalus arundinaceus in The Netherlands. Ardea 86(2):187-201.
Jedraszko-Dabrowska, D. 1992. Reeds as construction supporting Great Reed Warbler (Acrocephalus arundinaceus L.) and Reed Warbler (Acrocephalus scirpaceus Herm.) nests. Ekologia Polska 39:229-242.
Jiménez, J., J. M. Hernández, J. Feliú, M. Carrasco, and R. Moreno-Opo. 2018. Breeding in a dry wetland: demographic response to drought in the Common Reed Warbler Acrocephalus scirpaceus. Ardeola 65(2):247-259. https://doi.org/10.13157/arla.65.2.2018.ra4
Jiménez, J., R. Moreno-Opo, M. Carrasco, and J. Feliu. 2015. Estimating the abundance and habitat selection of conservation priority marsh-dwelling passerines with a double-observer approach. Ardeola 62(2):269-281. https://doi.org/10.13157/arla.62.2.2015.269
Kéry, M. and B. R. Schmidt. 2008. Imperfect detection and its consequences for monitoring for conservation. Community Ecology 9(2):207-2016. https://doi.org/10.1111/geb.12138
Kvist, L., S. Ponnikas, E. J. Belda, I. Encabo, E. Martínez, A. Onrubia, J. M. Hernández, P. Vera, J. M. Neto, and J. S. Monrós. 2011. Endangered subspecies of the Reed Bunting (Emberiza schoeniclus witherbyi and E. s. lusitanica) in Iberian Peninsula have different genetic structures. Journal of Ornithology 152:681-693. https://doi.org/10.1007/s10336-011-0646-9
Legendre, P. and L. Legendre. 2012. Numerical ecology. 3rd English edition. Elsevier Science BV, Amsterdam, Netherlands.
Llamas, M. R. 1988. Conflicts between wetland conservation and groundwater exploitation: two case histories in spain. Environmental Geology and Water Sciences 11(3):241-252. https://doi.org/10.1007/BF02574814
Leisler, B. 1991a. Acrocephalus melanopogon. Pages 217-251 in U.N. Glutz von Blotzheim, and K.M. Bauer, editors. Handbuch der Vögel Mitteleuropas. Aula, Wiesbaden, Germany.
Leisler, B. 1991b. Acrocephalus arundinaceus. Pages 488-539 in U.N. Glutz von Blotzheim, and K.M. Bauer, editors. Handbuch der Vögel Mitteleuropas. Aula, Wiesbaden, Germany.
López-Iborra, G. and J. Monrós, J. 2004. Bigotudo Panurus biarmicus. Pages 341-344 in A. Madroño, C. González, and J. C. Atienza, editors. Libro Rojo de las Aves de España. Ministerio de Medio Ambiente (MIMAM), Madrid, Spain.
Lüdecke, D., A. Bartel, C. Schwemmer, C. Powell, and A. Djalovski. 2020. sjPlot: Data visualization for statistics in social science. R Package Version 2.8.4. https://CRAN.R-project.org/package=sjPlot
Ma, Z., Y. Cai, B. Li, and J. Chen. 2010. Managing wetland habitats for waterbirds: an international perspective. Wetlands 30:15-27. https://doi.org/10.1007/s13157-009-0001-6
Maclean, I. M. D., J. P. Bird, and M. Hassall. 2013. Papyrus swamp drainage and the conservation status of their avifauna. Wetlands Ecology and Management 22:115-127. https://doi.org/10.1007/s11273-013-9292-8
Madroño, A., C. González, and J. A. Atienza. 2004. Libro Rojo de las Aves de España. Dirección General para la Biodiversidad - SEO/BirdLife, Madrid.
Martínez-Vilalta, J., A. Bertolero, D. Bigas, J-Y. Paquet, and A. Martínez-Vilalta. 2002. Habitat selection of passerine birds nesting in the Ebro delta reedbeds (NE Spain): management implications. Wetlands 22(2):318-325. https://doi.org/10.1672/0277-5212(2002)022[0318:HSOPBN]2.0.CO;2
Méro, T. O., A. Zuljevic, K. Varga, and S. Lengyel. 2015. Habitat use and nesting success of the Great Reed Warbler (Acrocephalus arundinaceus) in different reed habitats in Serbia. The Wilson Journal of Ornithology 127(3):477-485. https://doi.org/10.1676/13-200.1
Mestre, A., R. Poulin, R. D. Holt, M. Barfield, J. C. Clamp, G. Fernandez‐Leborans, and F. Mesquita‐Joanes. 2019. The interplay of nested biotic interactions, and the abiotic environment regulates populations of a hypersymbiont. Journal of Animal Ecology 88:1998-2010. https://doi.org/10.1111/1365-2656.13091
Middleton, B. A., B. Holsten, and R. van Diggelen. 2006. Biodiversity management of fens and fen meadows by grazing, cutting and burning. Applied Vegetation Science 9:307-316. https://doi.org/10.1111/j.1654-109X.2006.tb00680.x
Milsom, T. P., S. D. Langton, W. K. Parkin, S. Peel, J. D. Bishop, J. D. Hart, and N. P. Moore. 2000. Habitat models of bird species’ distribution: an aid to the management of coastal grazing marshes. Journal of Applied Ecology 37:706-727. https://doi.org/10.1046/j.1365-2664.2000.00529.x
Mitchell, L. R., S. Gabrey, P. P. Marra, and R. M. Erwin. 2006. Impacts of marsh management on coastal marsh birds’ habitats. Studies in Avian Biology 32:155-175.
Monrós, J. S., J. M. Neto, F. Arcos, P. Vera, and E. J. Belda. 2018. El escribano palustre en España y Portugal. Población reproductora en 2015 y método de censo. Universitat de València-Universitat Politécnica de València-SEO/BirdLife. Madrid. https://doi.org/10.31170/0059
Moreno, L., M. E. Jiménez-Hernández, H. Aguilera, P. Jiménez, and A. de la Losa. 2010. The 2009 smouldering peat fire in Las Tablas de Daimiel National Park (Spain). Fire Technology 47:519-538. https://doi.org/10.1007/s10694-010-0172-y
Morrison, D. A., R. T. Buckney, and B. J. Bewick. 1996. Conservation conflicts over burning bush in South-Eastern Australia. Biological Conservation 76:167-175. https://doi.org/10.1016/0006-3207(95)00098-4
Moskát, C. and A. Báldi. 1999. The importance of edge effect in line transect census applied in marshland habitats. Ornis Fennica 76:33-40.
Neto, J. M. 2006. Nest site selection and predation in Savi’s Warblers Locustella luscinioides. Bird study 53:171-176. https://doi.org/10.1080/00063650609461430
Paracuellos, M. 2008. Effects of long-term habitat fragmentation on a wetland bird community. Revue d'Écologie (La Terre et La Vie) 63:227-238.
Parsons, K. C., S. C. Brown, R. M. Erwin, H. A. Czech, and J. C. Coulson. 2002. Introduction: managing wetlands for waterbirds: integrated approaches. Waterbirds 25(suppl. 2):1-4.
Poulin, B. and G. Lefebvre. 2002. Effect of winter cutting on the passerine breeding assemblage in French Mediterranean reedbeds. Biodiversity and Conservation 11(9):1567-1581. https://doi.org/10.1023/A:1016866116220
Poulin, B., G. Lefebvre, and A. Mauchamp. 2002. Habitat requirements of passerines and reedbed management in southern France. Biological Conservation 107:315-325. https://doi.org/10.1016/S0006-3207(02)00070-8
Prokešová, J. and L. Kocian. 2004. Habitat selection of two Acrocephalus warblers breeding in reed beds near Malacky (Western Slovakia). Biologia, Bratislava 59(5):637-644.
Quesnelle, P. E., L. Fahrig, and K. E. Lindsay. 2013. Effects of habitat loss, habitat configuration and matrix composition on declining wetland species. Biological Conservation 160:200-208. https://doi.org/10.1016/j.biocon.2013.01.020
R Core Team. 2018. R: a language and environment for statistical computing. version 3.5.2. R Foundation for Statistical Computing, Vienna, Austria. www.r-project.org.
Raimondi, P., L. J. Jurgens, and M. T. Tinker. 2015. Evaluating potential conservation conflicts between two listed species: sea otters and black abalone. Ecology 96(11):3102-3108. https://doi.org/10.1890/15-0158.1
Rolando, A. and C. Palestrini. 1989. Habitat selection and interspecific territoriality in sympatric warblers at two Italian marshland areas. Ethology Ecology & Evolution 1:169-183. https://doi.org/10.1080/08927014.1989.9525521
Sánchez-Carrillo, S. and M. Álvarez-Cobelas. 2010. Climate and hydrologic trends: climate change versus hydrologic overexploitation as determinants of the fluctuating wetland hydrology. Pages 45-83 in S. Sánchez-Carrillo and D. G. Angeler, editors. Ecology of threatened semi-arid wetlands: long-term research in Las Tablas de Daimiel. Springer, Heidelberg, Germany.
Sánchez-Carrillo, S. and D. G. Angeler. 2010. Fish and Avian Communities: A Testimony of Wetland Degradation. Pages 197-212 in S. Sánchez-Carrillo and D. G. Angeler editors. Ecology of threatened semi-arid wetlands: long-term research in Las Tablas de Daimiel. Springer, Heidelberg, Germany.
Sánchez-Carrillo, S., M. Álvarez-Cobelas, S. Cirujano, M. Carrasco-Redondo, and A. Díaz-Cambronero. 2016. Long-term observations as a key tool for the assessment of environmental changes in Las Tablas de Daimiel: LTER-Daimiel. Ecosistemas 25(1):4-8. http://doi.org/10.7818/ECOS.2016.25-1.02
Sánchez-Carrillo, S., M. Álvarez-Cobelas, M. Merino-Ibarra, J. Ramírez-Zierold, and J. A. Morguí. 2021. Long-term nutrient dynamics in Las Tablas de Daimiel reveal the wetland has undergone enormous functional changes during the last 38 years (1980-2018). Limnetica 40(1). https://doi.org/10.23818/limn.40.11
Smart, J., J. A. Gill, W. J. Sutherland, and A. R. Watkinson. 2006. Grassland-breeding waders: identifying key habitat requirements for management. Journal of Applied Ecology 43:454-463. https://doi.org/10.1111/j.1365-2664.2006.01166.x
Tanneberger, F., C. Tegetmeyer, M. Dylawerski, M. Flade, and H. Joosten. 2009. Commercially cut reed as a new sustainable habitat for the globally threatened Aquatic Warbler. Biodiversity Conservation 18:1475-1489. https://doi.org/10.1007/s10531-008-9495-0ter
ter Braak, C. J. F. and P. Šmilauer. 2012. Canoco reference manual and user's guide: software for ordination, version 5.0. Microcomputer Power, Ithaca, United States of America.
Thirgood, S., S. Redpath, I. Newton, and P. Hudson. 2000. Raptors and Red Grouse: Conservation conflicts and management solutions. Conservation Biology 14:95-104. https://doi.org/10.1046/j.1523-1739.2000.99013.x
Timmermans, S. T. A., S. S. Badzinski, and J. W. Ingram. 2008. Associations between breeding marsh bird abundances and Great Lakes hydrology. Journal of Great Lakes Research 34:351-364. https://doi.org/10.3394/0380-1330(2008)34[351:ABBMBA]2.0.CO;2
Trnka, A. and P. Prokop. 2006. Reedbed structure and habitat preference of reed passerines during the post-breeding period. Biologia, Bratislava 61(2):225-230.
Trnka A., P. Batáry, and P. Prokop. 2009. Interacting effects of vegetation structure and breeding patterns on the survival of Great Reed Warbler Acrocephalus arundinaceus nests. Ardea 97(1):109-116. https://doi.org/10.5253/078.097.0113
Trnka, A., V. Peterková, P. Prokop, and P. Batáry. 2014. Management of reedbeds: mosaic reed cutting does not affect prey abundance and nest predation rate of reed passerine birds. Wetlands Ecology and Management 22:227-234. https://doi.org/10.1007/s11273-013-9325-3
Vadász, C., Á. Német, C. Biró, and T. Csörgő. 2008. The effect of reed cutting on the abundance and diversity of breeding passerines. Acta Zoologica Academiae Scientiarum Hungaricae 54(1):177-188.
Vera, P., E. J. Belda, L. Kvist, S. I. Encabo, M. Marín Villora, and J. S. Monrós. 2014. Habitat preferences for territory and nest settlement of the endangered eastern Iberian reed bunting Emberiza schoeniclus witherbyi. Ardeola 61(1): 97-110. http://doi:10.1357/arla.61.1.2014.97
Vera, P., E. J. Belda, L. Kvist, S. I. Encabo, and J. S. Monrós. 2011. Habitat preference of endangered Eastern Iberian Reed Buntings Emberiza schoeniclus witherbyi. Bird Study 58(3):238-247. http://doi.org/10.1080/00063657.2011.585393
Vera, P., M. Marín., E. J. Belda, and J. S. Monrós. 2009. Estructura y composición del nido del Escribano Palustre Iberoriental Emberiza schoeniclus witherbyi. Revista Catalana d’Ornitologia 25:43-48.
Vulink, J. T., H. J. Drost, and L. Jans. 2000. The influence of different grazing regimes on Phragmites and shrub vegetation in the well-drained zone of a eutrophic wetland. Applied Vegetation Science 3:73-80.
Wilson, J. and W. Peach. 2006. Impact of an exceptional winter flood on the population dynamics of bearded tits (Panurus biarmicus). Animal Conservation 9(4):463-473. https://doi.org/10.1111/j.1469-1795.2006.00063.x
Wood, S. 2016. mgcv: Mixed GAM computation vehicle with automatic smoothness estimation. R Package version 1.8-31. https://CRAN.R-project.org/package=mgcv
Zedler, J. B. and S. Kercher. 2005. Wetland resources: status, trends, ecosystem services, and restorability. Annual Review of Environment and Resources 30:39-74. https://doi.org/10.1146/annurev.energy.30.050504.144248