Ecology and Evolution. 2018;8:10489–10496. | 10489 www.ecolevol.org Received:14April2017 | Revised:24February2018 | Accepted:24March2018 DOI: 10.1002/ece3.4514 ORIGINAL RESEARCH Multiple components of environmental change drive populations of breeding waders in seminatural grasslands Karsten Laursen1 | Javier Balbontín2 | Ole Thorup3 | Henrik Haaning Nielsen4 | Tommy Asferg1 | Anders Pape Møller5 ThisisanopenaccessarticleunderthetermsoftheCreativeCommonsAttributionLicense,whichpermitsuse,distributionandreproductioninanymedium, providedtheoriginalworkisproperlycited. ©2018TheAuthors.Ecology and EvolutionpublishedbyJohnWiley&SonsLtd. 1DepartmentofBioscience,Aarhus University,Rønde,Denmark 2DepartamentodeZoología,Facultadde Biología,Sevilla,Spain 3AmphiConsult,Ribe,Denmark 4AvifaunaConsult,Vesløs,Denmark 5EcologieSystématiqueEvolution, CNRS,UniversitéParis-Sud,AgroParisTech, UniversitéParis-Saclay,Orsay,France Correspondence KarstenLaursen,InstituteofBioscience, AarhusUniversity,Rønde,Denmark. Email:
[email protected] Funding information AageV.JensenNaturfond Abstract Environmentsarerapidlychangingduetoclimatechange,landuse,intensiveagriculture,andtheimpactofhuntingonpredatorpopulations.Here,weanalyzedlong- termdatarecordedduring1928–2014onthesizeofbreedingpopulationsofwaders attwolargenaturereservesinDenmark,VejlerneandTipperne,todeterminethe effectsofcomponentsofenvironmentalchangeonbreedingpopulationsofwaders. Environmentalvariablesandcountsofwadersweretemporallyautocorrelated,and weusedgeneralizedleastsquare(GLS)byincorporatingthefirst-orderautoregressivecorrelationstructureintheanalyses.Weattemptedtopredicttheabundanceof wadersforshort-termtrendsfortwonaturereserves(35years)andforlong-term trendsforonenaturereserve(86years),usingprecipitation,temperature,nutrients, abundanceoffoxesVulpes vulpes,areagrazed,andnumberofcattle.Therewasevidenceofimpactsofnutrients,climate(long-termchangesintemperatureandprecipitation),grazing,mowing,andpredationonbirdpopulations.Weusedstandard randomeffectsmeta-analysesweightedby(N–3)toquantifythesemeaneffects. Therewasnosignificantdifferenceineffectsizeamongspecies,whilemeaneffect sizedifferedconsistentlyamongenvironmentalfactors,andtheinteractionbetween effectsizeforspeciesandenvironmentalfactorswasalsosignificant.Thus,environmentalfactorsaffectedthedifferentspeciesdifferently.Meaneffectsizewasthe largestat+0.20forrain,+0.11fortemperature,−0.09forfoxabundance,and−0.03 fornumberofcattle,whiletherewasnosignificantmeaneffectforfertilizer,area grazed, and year. Effect sizes for two short-term time series from Tipperne and Vejlernewerepositivelycorrelatedaswereeffectsizesforshort-termandlong-term timeseriesatTipperne.Thisimpliesthatenvironmentalfactorshadconsistenteffectsacrosslargetemporalandspatialscales. KEYWORDS climatechange,environmentalchange,fertilizer,landuse,long-termstudies,nutrients, precipitation,studymethods,temperature
10490 | LAURSEN Et AL. 1 | INTRODUCTION Breeding populations of waders are declining across western Europe(Thorup,2006;Roodbergen,vanderWert,&Hötker,2012; Robinson,Morrison,&Baillie,2014).Thereasonsforthesedeclines includeintensifiedagriculturalpractice,reclamationofcoastalhabitats,increasingpredationpressure,humandisturbance,andclimate change (Beintema & Müskens, 1987; Wilson, Ausden, & Milsom, 2004;Smart,Gill,Sutherland,&Watkinson,2006;Holm&Laursen, 2009;Roodbergenetal.,2012;Stephensetal.,2016).Demographic studiesshowthatthedeclinesaremainlycausedbypoorchicksurvivalratherthanadultsurvival(Roodbergen,Klok,&Schekkerman, 2008; Roodbergen etal., 2012). Low breeding success is caused bylossesofneststopredationandflooding(Hötker&Segebade, 2000;VandePoletal.,2010;Bellebaum&Bock,2009;Thorup& Koffijberg,2016).Studiesofwaderpopulationsduringrecentyears havefocusedonfarmlandrevealingcomplexinteractionsbetween agricultural practice and climate, affecting the survival of young (Kleijn etal., 2010; Schroeder etal., 2012). Young waders have tofindtheirownfood,andduetohighenergyrequirementsand poorabilitiestosaveenergy,theyoperatewithinnarrowenergetic margins, for example, a balance between quantity and quality of food (Schekkerman &Visser,2001;Maier,2013).Tothisarrayof parametersinfluencing waderpopulations,weincludedchange in nutrient load in the environment as a proxy for primary productivityandamountofbenthosandthusindirectlycarryingcapacity forbreedingbirds(Philippartetal.,2007;Møller,Flensted-Jensen, Laursen,&Mardal,2015).Thus,nutrientloadcanbeanimportant parameter,althoughastraightforwardrelationshipbetweennutrientsandamountofinvertebratefoodforyoungcannotbeexpected (Schekkerman & Beintema, 2007). Outside the breeding season duringmigrationandatwinteringsites,severalwaderspeciesstage and forage in marine estuaries and along coasts that are influencedbynutrients(Vitousek,Mooney,Lubchenko,&Melillo,1997; Windolf,Blicher-Mathiesen,Carstensen,&Krovang,2012). Different components of environmental change such as climatechange,landuse,andfisheriesarecurrentlybeingdescribed asdeterminantsofpoorreproductiveperformance,reducedsurvivorship,anddecliningpopulationtrendsforwadersandotherbird populations (Frederiksen, Wanless, Harries, Rothery, & Wilson, 2004; Schroeder etal.,2012;Maier,2013).However, fewstudies havedemonstratedforclimatechangethatdemographicvariables determine population size (Dunn & Møller, 2014; Robinson etal., 2014;Stephensetal.,2016).Attemptstomakeintegratedanalyses oftherelativecontributionsofmultiplefactorsaccountingforsuch trendsarescarce.Anumberofstudieshaverecentlyinvestigated theeffectsofclimatechangeandlanduseonpopulationsize(Møller, Flensted-Jensen,&Mardal,2007;Pimm,2009;Eglington&Pearce- Higgins,2012;Mantyka-Pringle,Martin,&Rhodes,2012;Jørgensen etal.,2015;Martin,vanDyck,Dendoncker,&Titeux,2013;Møller& Laursen,2015).Thisleavesanumberofadditionalfactorsinneedof study,includingpopulationchangesandeffectsofindustrializedagriculturewithhighlevelsoffertilizeruse.Sincesuchanalysesquickly includemanypredictors,theyareresourcedemandingandonlyfeasiblewhenbasedonlong-termstudies.Unfortunately,mosttimeseriesinecologicalresearchareshortandonlyrarelyexceed50years, therebyoftenpreventinginclusionofallorevenmostcrucialpredictorsandcertainlynotinclusionofinteractionsamongvariables. Theaimsofthisstudywere(a)toquantifytheeffectsizeforthe impactofenvironmentalconditionsonthepopulationsizeofseven waderspecies;(b)toanalyzetheeffectsofclimate,nutrients,land use,andpredatorabundanceonsizeoflocalbreedingpopulations ina86-yeartimeseriesofwaderbirdcommunitiesatTipperneand twoshort-termtimeseriesof35yearsatTipperneandVejlerne;(c) totestforconsistencyineffectsizebetweenoneshort-termand onelong-termtimeseriesatTipperneandforoneshort-termtime seriesofeffectsizesatTipperneandoneshort-termtimeseriesat Vejlerne.Wedidsobyanalyzingthebreedingabundanceofseven wader species during1928–2014 atthenaturereserveTipperne, Denmark, and during 1978–2014 at the nature reserve Vejlerne, Denmark. 2 | METHODS 2.1 | Sites, wader species, and study periods We studied breeding waders at two seminatural grassland sites (TipperneandVejlerne)inwesternDenmarkextensivelyfarmed,for example,grazedbycattleatlowdensity,latemowing,andwithout directuseoffertilizer(seesitedescriptioninSupportinginformation AppendixS1).Wefocusedonthewaderbirdcommunitycomposed of seven species: Oystercatcher Haematopus ostralegus, lapwing Vanellus vanellus,black-tailedgodwitLimosa limosa,redshankTringa totanus, dunlin Calidris alpina, ruff Philomachus pugnax, and avocet Recurvirostra avocetta. Data on bird numbers are available at Tipperneduring1928–2014andatVejlerne1978–2014,exceptfor 2004andforoystercatcherduring2004–2014(seeFig.1).Theyare allmigratorybirdswiththeearlybreedingspeciesarrivinginMarch (oystercatcher,lapwing),whiletheotherspeciesarriveinApril.The twostudyareasareNATURA-2000sitesthatarestrongholdsfor breedingwadersinDenmark(Thorup,2004). 2.2 | Census methods CensusmethodsaredescribedbyMøller(1983),Thorup(1998),and Kjeldsen(2008).Theyhavebeenadjustedduringthestudyperiod duetochangeinvegetationheight.AtTipperneduring1928–1957, nestsweresearchedintensively.During1958–1964,nestsearches weresupplementedbymappingbirdsgivingalarmcallstoidentify territories. In 1965–1985, breeding birds giving alarm calls were mapped together with nests. From 1986, breeding waders were mappedfromadistancewithatelescope(60×magnification)supplementedwithmappingeventsofwarningbirds.Thesemethods weregroupedinthreecategoriesandenteredinthestatisticalanalysesasafactorwiththefollowinglevels:nestsearch,nestsearch
| 10491 LAURSEN Et AL. andterritorymapping,andterritorymappingandtelescopeuse.A detaileddescriptionofthecensusmethodsisgiveninSupporting informationAppendixS1. 2.3 | Climate Climate estimated as long-term change in precipitation and temperaturetogetherwithgroundwaterlevelisknowntoaffectbreeding wader populations (Thorup, 1998; Bellebaum & Bock, 2009; Schroederetal.,2012;Maier,2013).Weusedmeanofdailymaximumairtemperature(°C)inAprilandsumofmonthlyprecipitation forMarch, April, and May (mm, datafrom DanishMeteorological Institute). We used data on climate measured at a local coastal weatherstationatVestervig.GroundwateratTipperneandwater level at the outlet sluice at Vejlerne were measured, but not includedintheenvironmentalanalysisduetosignificantcorrelations (p<0.05)withprecipitationinmostspringmonths. 2.4 | Nutrients Data on fertilizer from farmland in Denmark were estimated for 1928–2014 as the annual amount of outlet of total-N (tons) to coastalwaters(Conleyetal.,2007)andupdatesbyHansen(2011) andThodsenetal.(2016).Eachupdatewascalibrated totheformerlevel.Thisdatasetwasusedforthelong-termstatisticalanalysesforTipperne.From1989,dataonnitrogen(μg/L)concentration werecollectedinRingkøbingFjord(surroundingTipperne)andfor nitrogen in Limfjorden (adjacent to Vejlerne) as a part of the nationalmonitoringprogramNOVANA(datafromEnvironmentCentre Ringkøbing;bothsites;Hansen,2015).Thesedatasetswereused fortheshort-termstatisticalanalysesatTipperneandVejlerne. 2.5 | Management and predation Habitatmanagementasgrazingandmowingisimportantformaintaining populations of breeding waders while predators reduce breedingsuccess(Thorup,1998;Bellebaum &Bock,2009;Maier, 2013). At Tipperne, we have data on the number of cattle, area grazed(ha),andtherelativeareamowedsince1931.Duetoqualitativeinformationonareamowed,itwasclassifiedas0(0%–5%ofthe meadowsmowed),1(6%–25%mowed),2(26%–50%mowed),or3 (51%–75%mowed).AtVejlerne,dataonthenumberofcattle,the areagrazed,andtheareamowedwereavailable. We analyzed the effects of foxes on nest predation (Thorup, 1998;Kjeldsen,2008).Informationaboutfoxabundancewastaken from the official Danish bag statistics, which goes back to 1941. Fromthisyearupto1972,weusedthenationalbagsizeoffoxesasa proxyforfoxabundance(annualnumberoffoxesshotinDenmark). From1973to2014,informationexistsaboutfoxbagatthecountry level(annualnumberoffoxesshotpercountysurroundingTipperne andVejlerne),andweusedthislocalinformationtoaccountfordifferences between regions. Environmental data for Tipperne and VejlerneareshowninSupportinginformationTablesS1andS2. 2.6 | Statistical methods Weusedgeneralizedleastsquare(GLS)thatallowserrorstobe correlatedandtohaveunequalvariances.WetestedfortemporalautocorrelationforeachwaderspeciesbuildingtwoGLS, FIGURE1 Numberofbreedingpairsofsevenwaderspeciesat (a)Tipperneduring1928–2014and(b)Vejlerneduring1978–2014 (a) (b) FIGURE2 Boxplotsofeffectsizeforimpactsofdifferent environmentalconditionsonpopulationsizeandpopulation trendsofwaders.Theboxplotsshowmedians,quartiles,5and95 percentiles,andextremevalues –0.8 –0.4 0.0 0.4 0.8 Effect size Area grazed No. of cattle Amount No. of foxes AmountTemperature Year fertilizer rain
10492 | LAURSEN Et AL. oneof themincorporatingthe first-orderautoregressivecorrelation structure for year and the other one without it. We comparedbothmodelsusingANOVAandlikelihood-ratiotest (LRT) to evaluate which model performed better (Pinheiro & Bates,2000).Allanalysesshowedevidenceoftemporalautocorrelation(allp <0.01).Hence,weusedGLSthatincorporate the first-order autoregressive correlation structure for year toevaluate foreachwaderspecies seven predictors explaining variation in population size across years. The predictors weretemperatureinApril,totalprecipitationfromFebruaryto May,amountofnutrients,numberofcattle,totalareagrazed, numberoffoxes,andyear.Countnumbersforeachwaderspecieswerelog-transformedpriortoanalyses(exceptforlapwing countsthatshowedadistributionclosetoaGaussiancurve).All predictorswerestandardizedtoameanofzeroandSD =1to ensurethatallvariablesbeingatthesamescaleandtheinterceptbeinginterpretable.Residualsofeachmodelwerevisually inspectedfordeviationfornormalityusingnormalQQplots, andheteroskedasticity wasalsovisuallyinspected withplots ofresidualsagainstthefittedvalues.Webuiltmodelsforlong- termandshort-termdataatTipperneandshort-termdataat Vejlerne. We employed the libraries “nlme” (Pinheiro, Bates, DebRoy,&Sarkar,2018)andusedRversion3.3.3(RCoreTeam 2017). WeestimatedeffectsizesasPearson’sproduct–momentcorrelationcoefficientsbyusingstandardconversions(Rosenthal,1994). Theseeffectsizeanalysesonlyincludedfewobservations,implying thatthestatisticalpowerofanyspecificanalysisislow.Weadopted Cohen’s(1988)recommendationsforthemagnitudeofeffectsbeing small(Pearsonr=0.10,explaining1%ofthevariance),intermediate (r=0.30,explaining9%ofthevariance),orlarge(r=0.50,explaining 25%ofthevariance). 3 | RESULTS 3.1 | Effect size, environmental change, and species at Tipperne during 1928–2014 Overall effect size weighted by (N–3) was on average +0.033 (SE =0.016), N =147, marginally differing from zero (t =2.02, df =146, p =0.045). There was no significant difference in effectsizebetweenthetwosites(F =0.254,df =1,145,r2 =0.00, p =0.615(Figure2)). Effectsizedidnotdiffersignificantlyamongspecies(F =1.08, df =6, 140, r2 =0.0003, p =0.377). The four largest effect sizes were for foxes reducing the short-term abundance of avocet at Vejlerne(r =−0.552,t =2.90,df =26,p =0.0074),rainincreasingthe long-term abundance of lapwing at Tipperne (r =+0.427, t =3.58, df =66, p =0.0077), temperature increasing the long-term abundanceofgodwitatTipperne(r =+0.412,t =3.44,df =66,p =0.001), andfoxesreducingtheshort-termabundanceoflapwingatVejlerne (r =−0.339,t =3.42,df =26,p =0.0021).Short-termeffectsizefor godwitatVejlerneincreasedwithnutrientswithalargeeffectof +0.468, t =2.476, df =26, p =0.020, and avocet in the long term at Tipperne likewise had a strong positive effect size of +0.313, t =2.655,df =66,p = 0.0099. Effect size differed significantly among environmental variables (Figure3; F =5.797, df =6, 140, r2 =0.16, p<0.0001). The interactionbetweenspeciesandenvironmentalvariableswasalso statistically significant (Figure3; F =3.415, df =36, 140, r2 =0.44, p<0.0001). Thus, not all six environmental variables had similar effectsonthedifferentspecies.Precipitationhadanintermediate significantpositiveeffectonabundance,whiletemperaturehada significantpositiveeffect,foxesasignificantnegativeeffect,numberofcattleasignificantnegativeeffect,andareagrazedasignificantpositiveeffect(Figure3;Table1). FIGURE3 (a)Relationshipbetweeneffectsizesfordifferent charactersbetweenshort-termstudiesatTipperneandlong-term studiesatTipperne.Thelinesshowthepositivelinearregressions. (b)Relationshipbetweeneffectsizesfordifferentcharacters betweenshort-termstudiesatTipperneandshort-termstudiesat Vejlerne –0.8 –0.4 0.0 0.4 0.8 –0.8 –0.4 0.0 0.4 0.8 Short-term effect sizes Vejlerne Short-term effect sizes Tipperne (b) – 0.8 – 0.4 0.00.4 0. 8 – 0.8 – 0.4 0.0 0.4 0.8 Long-term effect sizes Tipperne Short-term effect sizes Tipperne (a)
| 10493 LAURSEN Et AL. 3.2 | Correlations between effect sizes at different sites and in periods of different duration For the data from Tipperne, there was a positive correlation between short- and long-term effect sizes (Figure 3a; F =10.187, df =1,47,r2 =0.16,p =0.0025,estimate(SE) =0.334(0.105),effect size+0.42).Thus,populationchangeshavebeenparallelatshortand at long terms.Similarly, there was a positivecorrelation between short-termeffectsizeatTipperneandshort-termeffectsatVejlerne (Figure3b; F =12.260, df =1, 47, r2 =0.19, p =0.0010, estimate (SE) =0.391(0.112),effectsize+0.45).Therefore,specificenvironmentalfactorshadsimilareffectsacrossspatialscales. 4 | DISCUSSION Themainfindingsofthislong-termstudyofpopulationtrendsof breedingwadersattwonaturereservesinDenmarksince1928and 1978weresignificantimpactsofmultiplecomponentsofenvironmentalchange.Inthepresentstudy,wetestedforconsistencyin effectsizesforenvironmentalvariableshypothesizedtoaffectpopulationtrendsofwaders.Wefocusedonasuiteofenvironmental changeparametershypothesizedtoactonbreedingwaderpopulations,althoughthesespeciesaremigrantsthatspendthenonbreedingseasonalongtheEastAtlanticFlywayfromtheWaddenSeain thenorthtoWestAfricainthesouth(Bønløkkeetal.,2006).The long-termtimeseriesdatingbackto1928allowedustotestforheterogeneity in strength of the relationship between abundance of breedingwadersandmultiplecomponentsofenvironmentalchange. Inparticular,wewereabletorankenvironmentalcomponentsin terms of importance. We found no significant difference among speciesineffectsize,butwedocumentedsignificantheterogeneityamongenvironmentalfactors,andtheseeffectsvariedamong species.Importantly, effectsizesfor thetwo shorttime seriesat TipperneandVejlernewerepositivelycorrelated,andthatwasalso thecasewhenrelatingeffectsizesfortheshorttermandthelong termatTipperne.Thus,effectsizesforenvironmentalfactorswere consistentacrosstemporalandspatialscales. Climatechangeiscurrentlyconsideredamajordeterminantof phenology, demography, distribution, and population trends (reviewsinMøller,Fiedler,&Berthold,2010;Pearce-Higgins&Green, 2014;Stephensetal.,2016).However,thereisonlylittleempirical evidencesuggestingthatpopulationsizeofbirdsisimpactedbyclimatechange(Dunn&Møller,2014;Stephensetal.,2016).Arrival dateandearlybreedinginwadersareclearlyaffectedbywarmer temperature in spring, although altered by agricultural practice (Petersen, Meltofte,&Tøttrup, 2012;Schroederetal.,2012).We showedconsiderablefluctuationsinbreedingpopulationsofwaders at Tipperne since 1928, explained by spring temperature and in particular precipitation. We documented a mean effect size of +0.197(SE =0.040)forspringprecipitationandaweakereffectsize of+0.110(0.044)forspringtemperature.Sincebothtemperature andprecipitationhaveincreased,thisimpliesthatthechangesinclimatearethosethatimprovedtheabundanceofbreedingwadersthe most.Temperatureandprecipitationareknowntoaffectreproductioninwaders(Schroederetal.,2012),andpopulationsresponding positivelytoclimatechangearelikelytoincrease(Møller,Rubolini,& Lehikoinen,2008;Stephensetal.,2016). Landusehaschangedconsistentlyinthetwostudysitesdueto management.Interestingly,wehaddetailedinformationonthearea grazedandthereforeonitsrelativeimportanceforbreedingwaders. Grazingreducedtheheightofthevegetation,sloweddownoreven reversedsuccessionalprocesses,andincreasedarthropodbiomass, which in turn had beneficial effects on growth of young (Norris etal.,1998;Eglingtonetal.,2008;Maier,2013).Inaddition,vegetationheightaffectstheabilityofwaderstoobservepredatorsata distance,butalsopreventsuseofsuitablenestingsitesandaccess topreferredforaginghabitatsforyoungduetotallswards(Cramp, 1985;Kleijnetal.,2010).Wefoundevidenceofaweakmeaneffectofthenumberofcattleonpopulationsizeofdifferentspecies ofwadersatthetwostudysites.However,theweakmeaneffect ofnumberofcattleof−0.032(0.026)impliesthatthismanagement toolinfacthadadetrimentaleffectratherthanabeneficialeffect onthenumberofbreedingwaders.Incontrast,theareagrazedhad ameaneffectof +0.031(0.030).This impliesthat thenumberof breedingwadersimprovedmarginallywiththeareagrazedbycattle. Weanalyzedtheeffectsoffoxesonpopulationsizeofbreeding waders that mainly act through effects on nest predation andtoasmallerextentonpredationonadultbirds(Beintema& Müskens, 1987; Meisner etal., 2014). We expected that foxes wouldnegativelyimpactpopulationtrendsofwaders.Indeed,we foundanexpectedmeaneffectof−0.106(0.050).Atotalof16 Variable Mean SE Wilcoxon W N p Year −0.014 0.053 −9,330 21 0.13 Amountof fertilizer +0.024 0.033 4,075 21 0.51 Areagrazed +0.031 0.030 33,871 21 <0.0001 No.ofcattle −0.032 0.026 −47,365 21 <0.0001 No.offoxes −0.086 0.041 −76,427 21 <0.0001 Temperature +0.110 0.044 82,291 21 <0.0001 Precipitation +0.197 0.040 19,652 21 <0.0001 TABLE1 Mean(SE)effectsizesfor environmentalfactorsatTipperneand Vejlerneandtestsfordifferencefromthe nullhypothesisofzeroinWilcoxon matched-pairssigned-ranktests
10494 | LAURSEN Et AL. outof21effectsizeswerenegative.Differencesineffectsoffox predationovertimeandbetweensitescouldbetheconsequence ofoutbreaksofthesarcopticmangediseasecausedbytheskin- dwellingmiteSarcoptes scabiei.Scabiescanreducefoxpopulations locallyandnationallywitheffectsonthebreedingwaterbirdcommunity(Forchhammer&Asferg,2000;Clausen&Kahlert,2010). There was a negative impact of foxes on waders, as shown by Roodbergenetal.(2012). Fertilizer increases nutrient availability, primary productivity, andsubsequenteffectsontheabundanceofanimalsathighertrophiclevelsinthemarineenvironment(Phillippartetal.,2007;Møller etal.,2015).Meaneffectsizeoffertilizeronwaderswasonly+0.024 (0.033),whichwasnotsignificantlydifferentfromzero.Here,we documentedtworelativelylargeeffects.Wefoundthateffectsize forgodwitatVejlerneincreasedwithfertilizerusewithalargeeffectsizeof+0.468,t =2.476,df =26,p =0.020,andavocetinthe long-termatTippernelikewisehadastrongpositiveeffectsizeof +0.313,t =2.655,df =66,p =0.0099.TipperneandthesurroundingRingkøbingFjordarestronglyimpactedbydrainagefromSkjern River,whichtransportshugeamountsofnutrientsfromlargesurroundingagriculturalareas(Petersenetal.,2008). Thisstudyhassignificantfutureprospects.First,weonlyanalyzed the main factors that most likely affect populations and populationtrendsofbreedingwaders.Clearly,itwouldbeinterestingtoassesstowhichextentinteractionsbetweenfactorsmay impactpopulationsandpopulationtrends.Second,thisstudyhas importantimplicationsforconservationintermsofmanagement ofnaturereserves,managementofpopulationsofpredatorssuch asfoxes,andmanagementofvegetationasaffectedbygrazing. Clearly,mostnaturereserves,oftensmallinsize,willbeimpacted byanumberoffactorsrelatedtoenvironmentalchange,butalso influenced by surrounding areasthat aremuch larger.Wehave shownherethateffectsizesforshort-andlong-termeffectsare positivelycorrelatedandthatdifferentfactorsaffectthesizeof breedingwaderpopulations. In conclusion, we have analyzed unique long-term data on breeding populations of waders with the longest time series starting in 1928. When analyzing the relationship between the abundanceofbreedingwadersandmultiplecomponentsofenvironmental change, we found significant evidence of effects. Whiletherewerepositivecorrelationsbetweeneffectsizesfor short-termstudiesatdifferentsites,orlong-andshort-termeffect sizesatthesamesite,thesecorrelationswererelativelystrongaccountingfor10%–15%ofthevariance.Thesefindingsemphasize theimportanceofconsideringheterogeneityinfuturemonitoring schemes. ACKNOWLEDGMENTS WethankJettePoulsenEngholm,EnvironmentCenterRingkøbing, for providing data for nutrients. O. Gordo provided constructive criticism.ThestudyreceivedfinancialsupportfromAageV.Jensen Naturfond.Wedeclarenoconflictofinterests. AUTHORS CONTRIBUTION APMconceivedtheproject;OTandHHNcollectedthedata;APM andKLanalyzedthedata;KL,APM,OT,HHN,andTAwrotethe manuscript. DATA ACCESSIBILITY WesuggesttostorethedataatTreeBASE,ifthepaperisaccepted. ORCID Karsten Laursen http://orcid.org/0000-0002-0380-8991 Anders Pape Møller http://orcid.org/0000-0003-3739-4675 REFERENCES Beintema,A.J.,&Müskens,G.J.D.M.(1987).Nestingsuccessofbirds breedinginDutchagriculturalgrasslands.Journal of Applied Ecology, 24,743–758.https://doi.org/10.2307/2403978 Bellebaum, J., & Bock, C. (2009). Influence of ground predators and waterlevelsonLapwingVanellus vanellusbreedingsuccessintwo continentalsites.Journal of Ornithology,150,221–230.https://doi. org/10.1007/s10336-008-0341-7 Bønløkke,J.,Madsen,J.J.,Thorup,K.,Pedersen,K.T.,Bjerrum,M.,& Rahbek, C. (2006). The Danish Bird Migration Atlas. Copenhagen, Denmark: Rhodos and Zoological Museum, University of Copenhagen. Clausen,P.,&Kahlert,J.(2010).Ynglefugle i Tøndermarsken og Margrethe Kog 1975-2009.FagligrapportfraDMUnr.778.(InDanish). Cohen,J.(1988).Statistical power analysis for the behavior science,2nded. Hillsdale,NJ:LawrenceErlbaum. Conley,D.J.,Carstensen,J.,Ærtebjerg,G.,Christensen,P.B.,Dalsgaard, T.,Hansen,J.L.,&Josefson,A.B.(2007).Long-termchangesandimpactofhypoxiainDanishcoastalwaters.Ecological Applications,17(5 Supplement),165–184.https://doi.org/10.1890/05-0766.1 Cramp,S.(Ed.)(1985).The birds of the Western Palearctic,Vol.4.Oxford, UK:OxfordUniversityPress. Dunn,P.O.,&Møller,A.P.(2014).Changesinbreedingphenologyand population size of birds. Journal of Animal Ecology, 83, 729–739. https://doi.org/10.1111/1365-2656.12162 Eglington,S.M.,Gill,J.A.,Bolton,M.,Smart,M.,Sutherland,W.J.,& Watkin,A.R.(2008).Restorationofwetfeaturesforbreedingwadersonlowgrassland.Journal of Applied Ecology,45,305–314. Eglington,S.M.,&Pearce-Higgins,J.W.(2012).Disentanglingtherelativeimportanceofchangesinclimateandland-useintensityindrivingrecentbirdpopulationtrends.Public Library of Science ONE,7(3), e30407. Forchhammer, M. C., & Asferg, T. (2000). Invading parasites cause a structuralshiftinredfoxdynamics.Proceedings of the Royal Society of London B,267,779–786.https://doi.org/10.1098/rspb.2000.1071 Frederiksen,M.,Wanless,S.,Harries,M.P.,Rothery,P.,&Wilson,L.J. ( 2 0 0 4 ) . T h e r o l e o f i n d u s t r i a l f i s h e r i e s a n d o c e a n o g r a p h i c i n t h e d e c l i n e ofNorthSeablack-leggedkittiwakes.Journal of Applied Ecology,41, 1129–1139.https://doi.org/10.1111/j.0021-8901.2004.00966.x Gill,J.A.,Norris,K.,Potts,P.M.,Gunnarsson,T.G.,Atkinson,P.W.,& Sutherland,W.J.(2001).Thebuffereffectandlarge-scalepopulationregulationinmigratorybirds.Nature,412,436–438.https://doi. org/10.1038/35086568 Goss-Custard,J.D.,&Stillman,R.A.(2008).Individual-basedmodelsand themanagementofshorebirdspopulations.Natural Resource Modelling, 21,3–71.https://doi.org/10.1111/j.1939-7445.2008.00006.x
| 10495 LAURSEN Et AL. Gunnarsson,T.G.,Gill,J.A.,Petersen,A.,Appleton,G.F.,&Sutherland, W. J. (2005). A double buffer effect in a migratory shorebird population. Journal of Animal Ecology, 74, 965–971. https://doi. org/10.1111/j.1365-2656.2005.00994.x Holm, T. E., & Laursen, K. (2009). Experimental disturbance by walkers affects behaviour and territory density of nesting Black-tailed Godwit Limosa limosa. Ibis, 151, 77–87. https://doi. org/10.1111/j.1474-919X.2008.00889.x Hötker,H.,&Segebade,A.(2000).Effectsofpredationandweatheron thebreedingsuccessofAvocetsRecurvirostra avosetta. Bird Study,47, 91–101.https://doi.org/10.1080/00063650009461163 Jørgensen,P.S.,Böhning-Gaese,K.,Thorup,K.,Tøttrup,A.P.,Chylarecki, P.,Jiguet,F.,…Rahbek,C.(2015).Continent-scaleglobalchangeattribution in European birds – combining annual and decadal time scale.Global Change Biology,22,530–543. Kjeldsen, J. P. (2008). Breeding birds in Vejlerne. Dansk Ornitologisk Forenings Tidsskrift,102,1–238. Klaasen, M., Lindström, Å., Meltofte, H., & Piersma, T. (2001). Arctic waders are not capital breeders. Nature, 413, 794. https://doi. org/10.1038/35101654 Kleijn,D.,Schekkerman,H.,Dimmers,W.,vanKats,R.J.M.,Melman, D.,&Teunissen,W.A.(2010).AdverseeffectsofagriculturalintensificationandclimatechangesonbreedinghabitatqualityofBlack- tailedGodwitsLimosa l. limosaintheNetherlands.Ibis,152,475–486. https://doi.org/10.1111/j.1474-919X.2010.01025.x Laursen,K.,&Hald,A.B.(2012).IdentificationofBlack-tailedGodwit’s Limosa limosa breeding habitat by botanical and environmental indicators. Journal of Ornithology, 153, 1141–1152. https://doi. org/10.1007/s10336-012-0845-z Laursen,K.,&Møller,A.P.(2013).Long-termchangesinnutrientsand mussel stocks are related to number of breeding Eider Somateria mollissimaatalargeBalticcolony.Public Library of Science One,9(4), e95851. Maier, M. (2013). Managing Mainland Salt Marshes for Breeding Birds. Thesis.Oldenburg,Germany:UniversityofOldenburg. Mantyka-Pringle,C.S.,Martin,T.G.,&Rhodes,J.R.(2012).Interactions betweenclimateandhabitatlosseffectsonbiodiversity:Asystematicreviewandmeta-analysis.Global Change Biology,18,1239–1252. https://doi.org/10.1111/j.1365-2486.2011.02593.x Martin,Y.,vanDyck,H.,Dendoncker,N.,&Titeux,N.(2013).Testing insteadofassumingtheimportanceoflandusechangescenariosto modelspeciesdistributionsunderclimatechange.Global Ecology & Biogeography,22,1204–1216.https://doi.org/10.1111/geb.12087 McClave,J.T.,&Sincich,T.(2003).Statistics,9thed.NewJersey,NJ: Prentice-Hall. Meisner,K.,Sunde,P.,Clausen,K.K.,Clausen,P.,Fælled,C.C.,&Hoelgaard, M.(2014).Foragingecologyandspatialbehavioroftheredfox(Vulpes vulpes)inawetgrasslandecosystem.Acta Theriologica,59,377–389. https://doi.org/10.1007/s13364-014-0178-9 Møller,A.P.(1983).Methods for monitoring bird populations in the Nordic countries.Oslo,Norway:NordicCouncilofMinisters. Møller,A.P.,Fiedler,W.,&Berthold,P.(Eds.)(2010).Effects of climate change on birds.Oxford,UK:OxfordUniversityPress. Møller, A. P., Flensted-Jensen, E., Laursen, K., & Mardal, W. (2015). Fertilizerleakagetothemarineenvironment,ecosystemeffectsand populationtrendsofwaterbirdsinDenmark.Ecosystems,18,30–44. https://doi.org/10.1007/s10021-014-9810-4 Møller,A.P.,Flensted-Jensen,E.,&Mardal,W.(2007).Agriculture, fertilizersandlifehistoryofacoastalseabird.Journal of Animal Ecology,76,515–525.https://doi.org/10.1111/j.1365-2656.2007. 01235.x Møller,A.P.,&Laursen,K.(2015).Reversibleeffectsoffertilizeruseon populationtrendsofwaterbirdsinEurope.Biological Conservation, 184,389–395.https://doi.org/10.1016/j.biocon.2015.02.022 Møller, A. P., Rubolini,D.,&Lehikoinen, E. (2008). Populationsofmigratorybirdspeciesthatdidnotshowaphenologicalresponseto climatechangearedeclining.Proceedings of the National Academy of Sciences of the United States of America,105,16195–16200.https:// doi.org/10.1073/pnas.0803825105 Norris, K., Brindley, E., Cook, T., Babbs, S., Brown, C. F., & Taxley, R. (1998). Is the density of Redshank Tringa totanusnestingonsaltmarshesinGreatBritaindecliningduetochangesingrazingmanagement? Journal of Applied Ecology, 35, 621–634. https://doi. org/10.1046/j.1365-2664.1998.355339.x Pearce-Higgins, J.W.,& Green,R.E.(2014).Birds and climate change: Impact and conservation responses. Cambridge, UK: Cambridge UniversityPress.https://doi.org/10.1017/CBO9781139047791 Petersen,J.K.,Hansen,J.W.,Laursen,M.B.,Clausen,P.,Carstensen, N. J., & Conley, D. J. (2008). Regime shift in a coastal marine ecosystem. Ecological Applications, 18, 497–510. https://doi. org/10.1890/07-0752.1 Petersen,T.L.,Meltofte,H.,&Tøttrup,A.P.(2012).AdvancedspringarrivalofavianmigrantsonTipperne,westernDenmark,during1929– 2008. Dansk Ornitologisk Forenings Tidsskrift,106,65–72. Petersen, J. W., & Petersen, D. L. J. (2011). Marine Områder 2010, NOVANA.ScientificreportfromDCEno.6.Aarhus,Denmark:Aarhus University(InDanish). Philippart,C.J.M.,Beukema,J.J.,Cadée,G.C.,Dekker,R.,Goedhart,P. W.,vanIperen,J.M.,…Herman,P.M.J.(2007).Impactsofnutrients reductiononcoastalcommunities.Ecosystems,10,95–118. Pimm,S.L.(2009).Climatedisruptionandbiodiversity.Current Biology, 19,R595–R601.https://doi.org/10.1016/j.cub.2009.05.055 Pinheiro, J., & Bates, D. (2000). Mixed-effects models in S and S-PLUS. New York, NY: Springer-Verlag. https://doi. org/10.1007/978-1-4419-0318-1 Pinheiro,J.,Bates,D.,DebRoy,S.,&Sarkar,D.&RCoreTeam.(2018). nlme:Linear and Nonlinear Mixed Effects Models.Rpackageversion 3.1-131.1,https://CRAN.R-project.org/package=nlme RCoreTeam.2017).R: A language and environment for statistical computing.Vienna,Austria:RFoundationforStatisticalComputing.https:// www.R-project.org/ Robinson,R.A.,Morrison,C.A.,&Baillie,S.R.(2014).Integratingdemographicdata:Towardsaframeworkformonitoringwildlifepopulationsatlargespatialscales.Methods in Ecology and Evolution,5, 1361–1372.https://doi.org/10.1111/2041-210X.12204 Roodbergen,M.,Klok, C.K., &Schekkerman,H.(2008). The ongoing declineofthebreedingpopulationofblack-tailedgodwitsLimosa l. limosainTheNetherlandsisnotexplainedbychangesinadultsurvival.Ardea,96,207–218.https://doi.org/10.5253/078.096.0206 Roodbergen,M.,vanderWert,B.,&Hötker,H.(2012).RevealingthecontributionsofreproductionandsurvivaltotheEurope-widedecline inmeadowsbirds:Reviewandmeta-analysis.Journal of Ornithology, 153,53–74.https://doi.org/10.1007/s10336-011-0733-y Rosenthal,S.(1994).Parametricmeasuresofeffectsize.InH.Cooper,& L.V.Hedges(Eds.),The handbook of research synthesis(pp.231–244). NewYork,NY:RussellSageFoundation. SASInstituteInc.(2012).JMP version 10.0.Cary,NC,SASInstituteInc.. Schekkerman,H.,&Beintema,A.J.(2007).Abundanceofinvertebrates andforagingsuccessofBlack-tailedGodwitLimosa limosachicksin relation to agricultural grassland management. Ardea, 95, 39–45. https://doi.org/10.5253/078.095.0105 Schekkerman, H., & Visser, G. H. (2001). Prefledging energy requirements in shorebirds: Energetic implications of self-feeding precocial development. The Auk, 188, 944–957. https://doi. org/10.1642/0004-8038(2001)118[0944:PERISE]2.0.CO;2 Schroeder, J., Piersma, T., Groen, N. M., Hooijmeijer, J. C. E. W., Kentie, R., Lourenco, P. M., … Both, C. (2012). Reproductive timing and investment in relation to spring warming and advancing
10496 | LAURSEN Et AL. agriculturalschedules.Journal of Ornithology,153,327–336.https:// doi.org/10.1007/s10336-011-0747-5 Smart, J., Gill, J. A., Sutherland, W. J., & Watkinson, A. R. (2006). Grassland-breedingwaders:Identifyingkeyhabitatrequirementsfor management.Journal of Applied Ecology, 43, 454–463. https://doi. org/10.1111/j.1365-2664.2006.01166.x Stephens,P.A.,Mason,L.R.,Green,R.E.,Gregory,R.D.,Sauer,J.R., Alison,J.,…Willis,S.G.(2016).Consistentresponseofbirdpopulations to climate change on two continents. Science, 352, 84–87. https://doi.org/10.1126/science.aac4858 Thodsen, H., Windolf, J., Rasmussen, J., Bøgestrand, J., Larsen, S.E., Tornbjerg,H.,…Wiberg-Larsen,P.(2016).Vandløb 2015.Videnskabelig rapportfraDCEnr.206.AarhusUniversitet.(InDanish). Thorup,O.(1998).ThebreedingbirdsonTipperne1928–1992.Dansk Ornitologisk Forenings Tidsskrift,92,1–192. Thorup, O. (2004). Status of populations and management of Dunlin Calidris alpina, Ruff Philomachus pugnax and Black-tailed Godwit Limosa limosainDenmark.Dansk Ornitologisk Forenings Tidsskrift,98, 7–20. Thorup, O.(Ed.)(2006).Breeding Waders in Europe 2000.International WaderStudies14.Slimbridge,UK:InternationalWaderStudyGroup. Thorup,O.,&Koffijberg,K.(2016).Breeding success in the Wadden Sea 2009–2012. A Review. Ecosystem No. 36.Wilhelmshaven,Germany: CommonWaddenSeaSecretariat. VandePol,M.,Ens, B.J.,Heg, D., Brouwer,L.,Krol,J.,Maier,M.,… Koffijberg,K.(2010).Dochangesinthefrequency,magnitudeand timingofextremeclimaticeventsthreatenthepopulationviabilityof coastalbirds?Journal of Applied Ecology,47,720–730. Vitousek,P.M.,Mooney,H.A.,Lubchenko,J.,&Melillo,J.M.(1997). Human domination of Earth’s ecosystems. Science, 277, 494–499. https://doi.org/10.1126/science.277.5325.494 Wilson,A.M.,Ausden,M.,&Milsom,T.P.(2004).Changesinbreeding waderpopulationonlowlandwetgrasslandsinEnglandandWales: Causesandpotentialsolutions.Ibis,146(Suppl.2),32–40.https://doi. org/10.1111/j.1474-919X.2004.00371.x Windolf,J.,Blicher-Mathiesen,G.,Carstensen,J.,&Krovang,B.(2012). Changesinnitrogenloadstoestuariesfollowingimplementationof governmentalactionplansinDenmark:Apairedcatchmentandestuaryapproachforanalyzingregionalresponses.Environmental Science and Policy,24,24–33.https://doi.org/10.1016/j.envsci.2012.08.009 SUPPORTING INFORMATION Additional supporting information may be found online in the SupportingInformationsectionattheendofthearticle. How to cite this article:LaursenK,BalbontínJ,ThorupO, HaaningNielsenH,AsfergT,MøllerAP.Multiplecomponents ofenvironmentalchangedrivepopulationsofbreeding wadersinseminaturalgrasslands.Ecol Evol. 2018;8:10489– 10496. https://doi.org/10.1002/ece3.4514