Patch and landscape responses of bird abundance to fragmentation in agroecosystems of east-central Argentina

Forest fragmentation in agroecosystems is linked to declines of avian species worldwide. Agriculture has greatly reduced native forest cover in east-central Argentina. Assessing the influence of fragmentation on forest bird populations is vital to inform reliable conservation and management strategies for the Espinal region of Argentina. We determined the relationships of vegetation structure within native forest patches, size and shape of these patches (patch scale), composition and spatial configuration (at landscape scale) to relative abundance of 17 forest bird species during austral fall-winter and spring-summer seasons. Birds were sampled from 2007–2009 in 45 forest patches within three landscape mosaics (30 x 30 km) representing a gradient of native forest fragmentation in east-central Argentina. We used an information-theoretic approach and approximated model inference to examine the effect of predictor environmental variables at two spatial scales on patterns of forest bird abundance. Density of trees within forest patches was the main predictor of bird abundance at the patch scale. Amount of forest and spatial configuration were also important. The abundance of several bird species was greater in patches with high density of trees and landscapes characterized by more forest cover and larger patches in close proximity. We found two main avian response patterns to forest fragmentation and patch characteristics. This information contributes reliable knowledge for the design of conservation measures of agricultural landscapes in the Espinal region of Argentina. Abondance d'oiseaux dans les îlots et les paysages à la suite de la fragmentation dans des agrosystèmes du centre-est de l'Argentine RÉSUMÉ. La fragmentation forestière dans les agrosystèmes est liée à la diminution d'espèces aviaires partout dans le monde. L'agriculture a grandement réduit la superficie de forêts indigènes dans le centre-est de l'Argentine. L'évaluation de l'influence de la fragmentation sur les populations d'oiseaux forestiers est cruciale pour que les gestionnaires puissent choisit des stratégies de conservation et de gestion fiables dans la région de l'Espinal en Argentine. Nous avons déterminé les relations de la structure de la végétation dans les îlots forestiers, de la taille et de la forme de ces îlots (échelle des îlots), de la composition et de la configuration spatiale (à l'échelle du paysage) avec l'abondance relative de 17 espèces d'oiseaux forestiers durant les saisons automne-hiver et printempsété austraux. Les oiseaux ont été échantillonnés de 2007 à 2009 dans 45 îlots forestiers se trouvant dans 3 mosaïques paysagères (30 x 30 km) représentatives d'un gradient de fragmentation de la forêt indigène dans le centre-est de l'Argentine. Nous avons utilisé une approche fondée sur la théorie de l'information et un modèle d'inférence flou pour examiner l'effet des variables environnementales explicatives sur la tendance de l'abondance des oiseaux forestiers aux deux échelles spatiales. La densité d'arbres dans les îlots forestiers était la principale variable explicative de l'abondance à l'échelle des îlots. La superficie de forêt et la configuration spatiale étaient aussi importantes. L'abondance de plusieurs espèces d'oiseaux était supérieure dans les îlots où la densité d'arbres était élevée et dans les paysages caractérisés par une plus grande superficie de couvert forestier et par de plus grands îlots près les uns des autres. Nous avons trouvé deux types principaux de réaction des oiseaux à la fragmentation forestière et aux caractéristiques des îlots. Ces résultats contribuent à une connaissance plus éclairée en vue d'élaborer des mesures de conservation de paysages agricoles dans la région de l'Espinal en Argentine.


INTRODUCTION
Habitat loss and fragmentation as consequences of land use changes are major drivers of global biodiversity loss (Baillie et al. 2004, Renfrew andRibic 2008).These processes in turn may modify the composition, structure, or function of landscapes (Fraterrigo et al. 2009, Magrach et al. 2011).Further, fragmentation may alter the original configuration of a landscape, generating isolated patches immersed in a frequently hostile matrix (Forman 1995, McComb 1999, Fahrig 2003).Both processes can negatively influence the abundance, movement, and Address of Correspondent: Noelia C. Calamari, Ruta 11, km 12.7, 3101, Oro Verde, Entre Ríos, Argentina, calamari.noelia@inta.gob.ardispersal of birds (Mazerolle and Villard 1999, Zurita and Bellocq 2007, Prugh et al. 2008).This may negatively influence demography and long-term persistence of bird populations, and affect key ecological roles of ecosystem function (Lundberg and Moberg 2003, Hanski 2011, Ye et al. 2013).
Forest ecosystems around the world are among the most affected by habitat loss and fragmentation, in many cases driven by the expansion and intensification of agricultural activities (Geist and Lambin 2002, Hoekstra et al. 2005, Gasparri and Grau 2009).In recent decades, the greatest loss of forests worldwide have occurred in the tropics, particularly in South America and Africa where annual rates of forest loss of 0.43% and 0.54%, respectively, have been reported (FAO 2016).With a loss of 297,000 hectares of forest/year (1.1%/year; FAO 2016), Argentina ranked as one of the top 10 countries with the greatest annual net loss of forest area during 2010 to 2015.In the province of Entre Ríos, located in east-central Argentina, the agricultural frontier has greatly expanded, despite soils poorly suited for production (Wilson 2008).Landscape composition and spatial configuration of native Espinal forest have changed during the past century.Forest cover in the province of Entre Ríos decreased from 56% in 1914 to approximately 13.8% at present (Muñoz et al. 2005, De la Fuente and Suárez 2008, Aizen et al. 2009).This has resulted in a subsequent loss of avian diversity in the region (Schrag et al. 2009, Gavier-Pizarro et al. 2012).
Bird species respond differently to environmental changes in fragmented landscapes, exhibiting both negative and positive effects (Antongiovanni andMetzger 2005, Uezu et al. 2005).Generally, negative effects are characterized by population declines from increased mortality, changes in resource availability or nesting sites (Lopez de Casenave et al. 1998, Zanette and Jenkins 2000, Bennett 2003).Fragmentation may also alter bird community structure through reductions of forest interior species, followed by an increase of generalist species associated with edge habitats (Merriam and Wegner 1992).Further, the increased isolation of remnant forest patches may compromise movement of species with reduced dispersal ability with potential negative consequences when remnant forest patches are immersed in a matrix of unsuitable habitat (Lens et al. 2002(Lens et al. , Şekercioglu et al. 2002)).Consequently, spatial scale plays a prominent role in the response of bird species to fragmentation (Holland and Bennett 2009).
The spatial scales at which birds respond to environmental characteristics may have conservation and management implications (Holland and Bennett 2009).For example, within patch, i.e., tree density and height, and patch-scale, i.e., patch size, features have been linked to population abundance and community structure of forest birds (Forman 2006, Skórka et al. 2016).Alternatively, factors characterizing the surrounding matrix can determine the main effects on bird community structure and composition (Gascon et al. 1999, Ribic and Sample 2001, Horn and Koford 2006).These hierarchal, i.e., patch to landscape, responses reflect the various environmental components required by forest birds to maximize fitness and complete their annual cycle (Manning et al. 2004, Thornton et al. 2011).Assessing the degree of sensitivity or tolerance at patch and landscape scales can provide insight into the effects of agriculture expansion on native forest birds.If we compare and integrate agricultural landscapes that represent a gradient of changes that occurred during the fragmentation process, trading space for time, we can determine bird response to agricultural expansion (dos Anjos 2006, Şekercioglu andSodhi 2007).
Agricultural expansion is expected to continue in the province of Entre Ríos with consequent alterations to native forest cover and configuration.Therefore, it is important to assess the influence of this multiscale process on one of the better indicators of ecological conditions, namely, native forest birds (Canterbury et al. 2000, O'Connell et al. 2000).In this study we analyzed the effect of forest fragmentation at patch and landscape scales on birds in east-central Argentina.Our objectives were to examine relationships between relative abundance of birds and environmental characteristics at patch and landscape scales and identify predictor variables influencing bird abundance.We predicted response by native forest birds to fragmentation in our study landscape would be a function of degree of habitat specialization and would differ by individual species along an autoecological gradient from sensitive to tolerant.Further, we predicted relative abundance of forest bird species would reflect variability of site-specific differences in patch characteristics.We used predictive habitat models to relate bird abundance to environmental characteristics of fragmented forest landscapes to assess their use for conservation planning (Mazerolle and Villard 1999, Henle et al. 2004, Thornton et al. 2011).

Study area
We conducted our study in the departments of Paraná, Diamante, Nogoyá, and La Paz of Entre Ríos province located in eastern Argentina.Entre Ríos is found within the Espinal ecoregion (Cabrera 1994).The Espinal ecoregion is characterized by semixerophytic forests dominated by hardwood tree species such as Prosopis affinis, Acacia caven, Geoffroea decorticans, and Celtis ehrenbergiana (Lewis et al. 2009, Matteucci 2012).Mean annual temperature is 18-20 °C and precipitation ranges from 800 to 1000 mm.Topography of the landscape is mostly flat plains interspersed with smooth undulations.The current state of Espinal forest remnants is the outcome of a process of secondary succession.Anthropogenic disturbances originated by indigenous peoples approximately 2000 years ago, but deforestation substantially intensified during the last 200 years.The most important changes in vegetation have occurred since the beginning of the 20th century, mainly from cattle ranching followed by industrial agriculture, mostly soybean (Muñoz et al. 2005, Engler and Vicente 2009, Matteucci 2012).
We placed a grid with a cell size of 30 x 30 km in each of the provincial departments encompassing our study area using a Geographic Information System.We classified all cells, using visual interpretation, according to the proportion of native forest cover into three categories (> 50%, 15-50%, < 5 %), then randomly selected one cell for each category.The selected landscape mosaics represented a gradient of forest fragmentation resulting from agricultural expansion processes (Fig. 1).Mosaic I (31°20'31.5"S, 59°26'43.2''W) included 54% native forest cover interspersed with agricultural patches.Mosaic II (31°34'52.2''S, 60°04'29.3''W) was dominated by an agricultural matrix with 23% forest cover connected by corridors of riparian forest.Finally, mosaic III (32° 02'2.6''S, 60°18'29.9''W) was dominated by intensive crop agriculture with interspersed forest fragments representing 4% of the total area.Each landscape mosaic was subdivided into nine nonoverlapping 10 x 10 km cells using ArcGIS version 9.2 (Esri, Redlands, CA, USA).We randomly selected six cells, hereafter termed sublandscapes.Within each sublandscape we randomly selected forest patches in three size categories (< 7 ha, 8-20 ha, > 21 ha), 15 in each mosaic for a total of 45 patches.A forest patch was defined as a contiguous area covered by native forest that differed from its surroundings.Average size of these 45 native forest patches (± standard error) was 37.16 ± 5.25 ha.However, mean patch area was 253.95 (± 40.03), 64.36 (± 12.85), and 15.12 (± 4.19) for mosaics I, II, and III, respectively.Following a first visit to each patch, we planned our bird and vegetation structure sampling scheme.

Bird community
We sampled forest birds during the austral seasons of fall-winter 2007 and 2008 (May to July) and spring-summer 2007-2009 (November to January).We selected 19 diurnal, terrestrial bird species to derive abundance estimates (Table A1.1).Selected taxa were common species representative of the Espinal forest region (Narosky and Yzurieta 2010).Further, selected species illustrated a diversity of foraging guilds, e.g., insectivores and granivores, and habitat use preferences of varying forest strata (understory, midstory, canopy).Common and taxonomic names of bird species followed De la Peña and Rumboll (1998).
We used 10-minute circular point counts (50-m radius) to estimate relative abundance, expressed as number of birds per point per patch (Savard andHooper 1995, Bibby et al. 2000).We randomized distance between transects and placed points along transects perpendicular to the edge of each forest patch.The number of points on a given transect was proportional to patch size (range = 4-16 points per patch) and points were located at a minimum distance of 100 m to decrease the possibility of double counting (Bibby et al. 2000;Fig. 1).Surveys were conducted between 0600-1000 h (spring-summer), 0700-1100 h (fall-winter) or 1500-1900 h (spring-summer), 1400-1800 h (fall-winter), with one patch sampled per morning-afternoon period.We visited each point once per season.Surveys were not conducted during rainy or windy days (wind ≥ 20 km/h).During the first set of surveys, the number of observations was greater in the morning block ( = 253, SE = 69.30)compared to the afternoon ( = 181, SE = 29.79).Therefore, we kept the same block observation in the first two surveys and reversed the block schedule in the last two surveys.This decision was made to avoid masking possible seasonal influences on bird detections (Bibby et al. 2000).We attempted to estimate detection functions for all bird species in native forest patches based on distance data and following Buckland et al. (2001).However, for most species this approach was not feasible primarily because of an insufficient number of independent detections to fit the detection function, thus violating a key assumption of distance sampling.Therefore, abundance estimates were reported as relative abundance ± standard error.

Patch-scale variables
We selected vegetation structure variables to measure characteristics within forest patches based on the literature (García Del Rey and Cresswell 2005, Johnston and Holberton 2009).We used the dot-quadrant method (Cottam and Curtis 1956) at each point to quantify tree density (D_Arb) and shrub density (D_arb) as number of trees/shrubs per hectare.We also documented additional measures collected from the tree closest to the center of the plot in each quadrant, including tree height (A_Arb) using a clinometer (Romahn de la Vega et al. 1994), and diameter at breast height (DBH).These measures were averaged over all sampling points in a patch to obtain a single value representing patch scale vegetation structural variables.
We considered additional patch level variables known to affect bird populations (i.e., breeding habitat), community structure, or spatial and temporal distribution (Estades 2001, Graham and Blake 2001, Suarez-Rubio and Thomlinson 2009).These included patch area (PA), total edge (TE), shape index (SI), perimeter-area ratio (PARA), fractal dimension (FD), linear distance of patch to nearest settlement (DistLoc), latitude (X), and longitude (Y) for the central point of each forest patch (Gauss-Kruger, zone 5).The last two variables were included in the analysis to account for potential landscape scale spatial patterns masking relationships between bird abundance and patch characteristics, and account for potential correlation of other variables responding to the same landscape scale gradients.
Variables were measured using Quickbird images (available in Google Earth™ http://earth.google.com)and the Patch Analyst extension for ArcGIS (Rempel et al. 2012).These described the amount and configuration of each forest patch beyond those that described environmental characteristics within patches.Possible correlations among patch scale variables were tested using a Pearson's matrix and discarded variables that showed high levels of correlation (r ≥ 0.70, P < 0.05; Table 1).Images were geometrically corrected using a first-degree polynomial model given study area topography.Images were orthorectified using a digital elevation model (DEM, 30-m spatial resolution) and nearest neighbor as a resampling method.Vector topology was reconstructed using a PAEK algorithm (polynomial approximation with exponential Kernel) with a 50 m tolerance (Menéndez and La Roca 2007), then vectors rasterized (5-m cell size) in ArcGIS.
Landsat™ satellite images were classified with support vector machines (SVM), a nonparametric supervised classification algorithm (van der Linden et al. 2009).The SVM approach distinguishes classes by fitting separating hyperplanes in feature space based on training samples (Foody and Mathur 2004).After training sites were digitized, statistical characterizations were derived and 10 land use-cover types identified: water bodies (W), corn (C), sunflower (S), soybean (S), sorghum (So), other crops (Oc) included annual and perennial pastures and fields with weeds, introduced forest (Pf), native forest (Fn), flooded vegetation (Vf) and urban-bare soil (URS).
We trained the SVM using 250 randomly selected pixels of each class.Although the SVM approach was originally developed for binary classifications, we worked with multiclass images overcoming that problem through a one-against-one approach that applied a set of individual classifiers to all possible class pairs and performed a majority vote to assign the winning class (Huang et al. 2002, Pal andMather 2005).Further, we determined magnitude of penalty given to misclassified training data (C) and width of the Gaussian kernel function selected (γ).Once optimal parameters were found, we used the resulting SVM to classify the images (Janz et al. 2007, Kuemmerle et al. 2008).We eliminated pixel groups smaller than 0.81 ha (3 x 3 pixels) representing classification artifacts or areas without ecological importance to the scale of our study.We validated classification results with 100 randomly selected points per land class type using Quickbird images and ground truthing (Congalton and Green 2009).Overall classification accuracy was 82% ( 2007) and 84% (2008; Table A2.1).
We estimated forest fragmentation in each sublandscape by calculating metrics that quantified amount and edge of forest cover, namely total forest area in each mosaic (FA, A3.1), number of forest patches (NP, A3.2), mean patch area (PA, A3.3), and edge density (ED, A3.4).Additionally these metrics described shape complexity of forest patches with measures of shape index (SI, A3.5), fractal dimension index (FD, A3.6), and perimeterarea ratio (PARA, A3.7).Finally, they addressed spatial configuration of the forest via mean distance (Euclidean) to nearest neighbor patch (ENN, A3.8), patch cohesion index (COH, A3.9), and aggregation index (AI, A3.10).These are considered robust metrics for fragmentation measurements and allow the definition of patterns (Wang et al. 2014).When constructing our models, we averaged metrics for 2007 and 2008 as there was no difference among years (P > 0.05).We used area of each cover type, contagion (CO), interspersion and juxtaposition (IJI), and richness (PR) to characterize the surrounding matrix on each sublandscape.All landscape-scale metrics were calculated using FRAGSTATS version 4 (McGarigal et al. 2012).

Analysis
We used Kruskal Wallis tests (P > 0.05) to evaluate differences at patch and landscape scale variables between mosaics.The relationships between the environmental metrics selected to characterize forest patches indicated variables PA and TE were highly correlated (Table 1).Therefore, we selected PA for model building given its wide use in other studies of species-habitat relationships (Munguía-Rosas and Montiel 2014).Variables describing forest patches PARA and FD were also strongly correlated.Therefore, we included only FD in our models because this index reflects shape complexity across a range of spatial scales (patch sizes; McGarigal et al. 2012; Table 1).We discarded X and Y given these were autocorrelated, although were uncorrelated with any of the other explanatory variables.Therefore, this discards any possible spatial patterns of patch scale variables along the landscape fragmentation gradient.Further, we used a correlation approach to examine relations among cover type metrics at the landscape scale.Because forest cover landscape metrics were correlated (r ≥ 0.70, P < 0.05; Table 2), we used a principal components analysis (PCA) to obtain independent measurements of these variables (McCune et al. 2002).The resulting two vectors accounted for 86% of the variance in the original data matrix (PC1 63%, PC2 23%).Axis PC1 represented a gradient of sublandscapes dominated by large patches of connected forest with less edge and simpler shapes to sublandscapes with reduced forest cover characterized by smaller patches of forest with greater edge and complex shapes.Axis PC2 represented a gradient of covariation among sublandscapes with low number of large patches to sublandscapes with greater number of smaller patches and finally, sublandscapes characterized by low patch number of small size.Relative abundance of each bird species per season was modeled as a function of patch and landscape variables.We built two sets of models independently.These included models where patch variables were considered, and models with landscape variables only.We selected this approach given differences in sample sizes, with 45 patches for the first set of models and 18 (six sublandscapes in each of three mosaics) at the landscape scale.
Based on this criterion, each model included only one to three explanatory variables.
The response variable (relative bird abundance) and explanatory variables (patch and landscape metrics) were examined graphically and analytically for multicollinearity and deviation from normality using Infostat (Di Rienzo et al. 2010).We then used generalized linear models (GLM; Guisan et al. 2002, Seavy et al. 2005) with a Gaussian error structure (log-transformed when condition was not met) to examine the influence of environmental variables on relative abundance of each bird species using package glm in program R ver.2.14 (R Development Core Team 2013).We first constructed models for each single variable, followed by models with sets of two and three variables within each model set (Table A4.1).We used this exploratory approach, rather than a more restricted set of a prior models (Fletcher and Koford 2002) given the absence of prior information on the explanatory power of our variables.We used Akaike's information criteria adjusted for small sample size (AICc) to compare model performance and considered Δ AICc scores ≤ 2 to be competitive models (Burnham and Anderson 2002).
For each spatial scale analyzed we used model averaging and sum of weights of competitive models where a variable was present (∑ωi; Burnham and Anderson 2002) to evaluate individual variable performance at each scale (patch and landscape).
Variables with a value of ∑ωi close to 1 were the most important predictor variables.Additionally, we used a hierarchical partitioning analysis to calculate the percent variance of the full model explained by each variable when all other variables were included in the model.For estimates of avian abundance, we fitted all possible models based on different combinations of the explanatory variables.For each fitted model the variable of interest was removed before refitting the model.Variable importance was calculated as the average change in R² when a variable was removed from all fitted models (MacNally 2002).Finally, we assessed spatial autocorrelation using randomized semivariograms of model residuals (Isaaks and Srivastava 1989).
There was no evidence of a spatial autocorrelation effect in our models.

Bird species abundance
We removed Greyish Saltator (Saltator coerulescens) and Shortbilled Canastero (Asthenes baeri), from subsequent analysis because of the low number of detections.Relative abundance estimates differed for most species by season (fall-winter, springsummer) and landscape mosaic (Fig. 2

Mosaic description
As expected, tree density was greater in mosaic I (= 355 ± 36 indiv/ ha) and mosaic II ( = 451 ± 25 indiv/ha) compared with mosaic III ( = 166 ± 29 indiv/ha, Table A5.1).However, tree height was greater and statistically significant only in mosaic II ( = 5.40 ± 0.18 m) compared with mosaic I ( = 4.44 ± 0.18 m).Shrub density did not differ between mosaics (mosaic I = 528 ± 84 indiv/ha, mosaic II = 462 ± 79 indiv/ha and mosaic III = 289 ± 44 indiv/ha).Forest amount and configuration differed (P < 0.05) between mosaics.Mosaic I was dominated by continuous forest cover while mosaic III had the lowest amount of forest (Table A5.1).Mosaic I was mainly characterized by 54% native forest, 38% annual crops, whereas mosaic II had 23% native forest, 66% annual crops, and finally mosaic III had 4% native forest and 77% annual crops.Additionally, mosaic II had the greater number of forest patches indicating an intermediate stage in the fragmentation process, whereas mosaic III had the lowest number of patches; its mean area was also the lowest and its mean shape of patches was the most complex (Table A5.1).
We observed that some variables presented associations among spatial scales, despite evaluating relationships between bird abundance and environmental metrics at both patch and landscape scales using independent sets of models (Table A6.1).
On the one hand, these associations indicated that patches with greater tree height were found in landscapes characterized by a greater number of patches.On the other hand, the shape complexity of the 45 sampled patches was positively associated with the shape pattern at landscape scale.Finally, the landscapes characterized by greater number of forest patches (included in axis PC2) were also characterized by the more disaggregated land use-cover types (Table A6.1).

Patch and landscape scales bird response
Although the response to forest fragmentation varied among avian species, we identified two main patterns: (1) species sensitive to vegetation structure conditions within patch and forest fragmentation (e.g., Narrow-billed Woodcreeper, Scimitar-billed Woodcreeper, and Brown Cacholote (Pseudoseisura lophotes), and (2) tolerant species (e.g., Picui Ground Dove, Saffron Finch, and Masked Gnatcatcher).
Vegetation structural variables were incorporated in most individual bird models (Table 3).Number of trees per hectare was important when determining abundance of bird species during both seasons.Variables related to tree size (DBH and A_Arb) were not as important as tree density (Table 4).Tree height had a strong and positive effect on several species including Sootyfronted Spinetail (Synallaxis frontalis), Great Antshrike, Creamybellied Thrush, Golden-billed Saltator (Saltator aurantiirostris) and a negative effect on White-crested Tyrannulet (Serpophaga subcristata), Saffron Finch, and Monk Parakeet.Shrub density had a clear negative effect on only three species; Rufous-bellied Thrush (Turdus rufiventris), Grayish Baywing, and Picui Ground Dove.The abundance of Narrow-billed Woodcreeper, Scimitarbilled Woodcreeper, Great Antshrike, and Creamy-bellied Thrush increased proportionally with density of trees within forest patches.Conversely, Grayish Baywing, Picui Ground Dove, Black-capped Warbling Finch, and Saffron Finch were more abundant in patches with lesser tree densities (Table 3).
Variables characterizing patch shape complexity, such as shape index (SI) and fractal dimension (FD) were incorporated in most individual bird models during both seasons (Table 3).Furthermore, relative importance based on model weights of SI (range = 0.28-0.98)and FD (range = 0.24-1.00)was greater for most bird species (Table 4).For example, Scimitar-billed Woodcreeper, Narrow-billed Woodcreeper, Grayish Baywing, Red-crested Cardinal, and Brown Cacholote were more abundant in patches of simpler shape and lesser ratio of edge to forest interior.Conversely, Picui Ground Dove, Great Antshrike, Saffron Finch, Sooty-fronted Spinetail, Masked Gnatcatcher, and Golden-billed Saltator were more abundant in patches of greater shape complexity and greater ratio of edge to forest interior.Proximity to population settlements did not have a clear effect on any bird species, except for Monk Parakeet that was clearly favored by closeness to settlements (Table 3).
Birds also responded to amount and spatial configuration of forest at landscape scale as forest fragmentation metrics were present in the best models (Table 3).Variables summarized by PC1 exhibited greater mean relative importance (0.68) than variables summarized by PC2 (0.41; Table 4).Furthermore, PC1 was found in most species landscape scale models.The relationship of individual bird species to PC1 varied depending on their affinity to forest habitat, ranging from species that feed, reproduce, and nest within the forest, i.e., forest specialists, to open habitat generalists that use the forest, but also the surrounding matrix.Conversely, PC2 was selected in a smaller number of models, some of these included bird species that exhibited a positive response to landscapes with greater number of patches of native forest, e.g., Grayish Baywing, Great Antshrike, and Rufous-bellied Thrush, while other species such as Black-capped Warbling Finch, Creamy-bellied Thrush, Whitecrested Tyranulet, and Monk Parakeet responded negatively (Table 3).
Finally, the surrounding matrix influenced many bird species such as the Golden-billed Saltator, who exhibited a positive effect of interspersion and juxtaposition (IJI), and whose abundance was greater in landscapes where patch types were close together.Conversely, other species were negatively influenced by IJI including the Scimitar-billed Woodcreeper and White-crested Tyranulet (Table 3).On the other hand, abundance of Rufousbellied Thrush and Monk Parakeet was positively associated to contagion, and was greater in landscapes where all patch types were aggregated, while Narrow-billed Woodcreeper, Sootyfronted Spinetail, and Creamy-bellied Thrush had a negative response to CO (Table 3).Finally, the abundance of Grayish Baywing, Red-crested Cardinal, and Brown Cacholote was greater in landscapes with fewer types of different land use-cover (PR; Table 3).Nevertheless, relative importance (∑ωi) of landscape predictor variables describing the matrix was low (< 0.5; Table 4).

DISCUSSION
Forest structure, amount, and spatial configuration successfully explained patterns of bird abundance in fragmented landscapes of the Espinal region in east-central Argentina.Although some species exhibited no apparent response or a positive response to fragmentation, a subset of bird species responded negatively to these patterns at both patch and landscape scales.Thus, potential changes in bird communities in east-central Argentina could be driven by reduced abundance of forest sensitive species.This in turn may enhance the potential for population declines, including local extinctions, as landscapes become increasingly fragmented and simplified (Boulinier et al. 2001).These results were  consistent with other published reports (Gehring and Swihart 2003, Cunningham and Johnson 2006, Thornton et al. 2011) on avian responses to forest fragmentation.Additionally, our results highlighted the importance of evaluating species-habitat relationships at multiple, i.e., patch and landscape, scales because factors varied in importance depending on the spatial scale.For instance, at the patch scale, vegetation structure within patches was comparatively more important than patch shape and size, while at the landscape scale the amount and spatial configuration of native forest cover were comparatively more important than composition and configuration of all cover types present in the landscape.
Vegetation structure within patches was an important predictor for many species in our study, similarly to patterns reported for other regions of the world (Martin 1998, Forman 2006).Tree density within patches was strongly associated (positively or negatively) with 15 out of 17 bird species in our study during both fall-winter and spring-summer.Increased density of trees would benefit cavity nesting species like the Narrow-billed Woodcreeper and Scimitar-billed Woodcreeper by increasing the availability of potential nesting and foraging sites, suggesting these species may be associated to particular forest patch conditions (Nadkarni and Longino 1990, Berg et al. 1994, Cockle et al. 2012).Interestingly, these results differed from studies conducted in North America where an inverse relationship to tree density has been reported for cavity nesting birds (Flack 1976, Brawn 1988).Most North American cavity nesters are more abundant in mature forests dominated by larger trees (and lower tree density), probably because older forests tend to have more dead trees and dead branches on live trees where cavities may be excavated or utilized.Many other species are secondary cavity nesters and rely on abandoned woodpecker nests or natural cavities.In the Espinal forest, bird species that nest in tree cavities are secondary users, Table 4. Relative importance of spatial metrics at patch and landscape scales based on sum of model weights (∑ωi, for each predictor variable shows the sum of Akaike weights for all possible models in which the predictor variable was incorporated at each level) and hierarchical variance partitioning (HP, percent of variance of the full model explained by each variable) for bird species sampled during the austral fallwinter and spring-summer, in the province of Entre Ríos, Argentina.taking advantage of holes made by medium-sized woodpeckers.The existing Espinal forest is mostly secondary, where stands commonly include trees of different sizes and ages.In summary, the forest patch characteristics and historical land use, e.g., proportion of shrub cover, degradation, livestock ranching, and extraction of mature trees, may explain the lack of an observed inverse relationship between tree density and availability of nest cavities.Regarding other species like Monk Parakeet and Picui Ground Dove, preference for patches characterized by reduced tree density could facilitate the detection of potential nest predators (Peris andAramburú 1995, Eberhard 1998).
Most species in our study did not respond to forest patch size contrary to other studies that reported a response (positive or negative) between bird abundance and patch size (Robinson et al. 1995, Estades 2001, Lindenmayer et al. 2002, Ferraz et al. 2007).This suggests ecological processes manifested at certain spatial scales may be of greater importance.For example, patch size response may be influenced by species-specific characteristics such as life history strategies or attributes of the landscape.Similarly, abundant species may be present in small patches (MacNally et al. 2000), other species may increase in abundance in proportion with patch size (Vickery et al. 1994, Johnson andIgl 2001), while others may exhibit no definite pattern (Andrén 1994, Estades andTemple 1999).In our study, the relatively weak influence of patch size could be related to the absence of clear statistical differences between landscapes.However, this relationship was identified, i.e., PC1, at the landscape scale, suggesting large tracts of continuous forest in mosaic 1 could still support vulnerable bird species or those sensitive to fragmentation.Conversely, these species were not abundant in landscapes characterized by relatively small forest patches surrounded by widespread agricultural development.These patterns suggest ecological thresholds, i.e., response, may exist for species requiring particular forest habitat conditions (Luck 2005, Betts et al. 2010), such as the Narrow-billed Woodcreeper, Scimitar-billed Woodcreeper, or Brown Cacholote.
Patch shape complexity was consistently important, a result similar to other published reports regarding its relationship to bird abundance (Hawrot andNiemi 1996, Magrach et al. 2011).Simpler patch shapes in our study area favored cavity nesting birds like Woodcreepers and species like Red-crested Cardinal and Suiriri Flycatcher.Patch shape is determined by the interaction between area and perimeter, which in turn determines the amount of core habitat in a fragmented landscape (Laurance andYensen 1991, Collinge 1996).We did not quantify core area, given the extensive information in the published literature regarding core area differences among species (Watson et al. 2004, Zurita et al. 2012).Further, edge effect relationships to bird species were not part of our study objectives.Additionally, edge was not clearly defined in the Espinal forest compared to other forested biomes of southern South America, such as the Chaco forests of Argentina (Lopez de Casenave et al. 1998, Sosa 2008, Banks-Leite et al. 2010, Zurita et al. 2012).The absence of a clear edge in Espinal forest may be related to the semixerophytic nature of these forests, dominated by plant species adapted to harsh environmental conditions typical of edges, and the influence of cattle grazing in Entre Ríos province (Cano et al. 1980, Sabattini et al. 2002).
Previous information on bird use of forest edge versus interior differed among published studies.Sosa (2008) identified Masked Gnatcatcher as a forest interior species, yet in our study this was more abundant in landscapes characterized by small forest patches with little to no core habitat.Conversely, Golden-billed Saltator was reported as exclusively associated with forest edge (Sosa 2008).However, it was less abundant in study landscapes characterized by small patches and therefore lesser forest interior area.Furthermore, Dardanelli et al. (2006) differentiated use of forest interior and edge by some of the same species in our study and reported contrasting patterns to results by Sosa (2008).These differences may be due to the fact that while these other studies were also conducted in the Espinal region, the composition of tree species differed from our study area.Weather conditions in our study region also differed from Sosa (2008).Further, forest management practices in relation to fire and extractions during the first half of the 20th century resulted in changes to vegetation structure (Dussart et al. 2011).Most importantly, spatial scale may explain the differences among our study and Sosa (2008), given the area of this previous study was 9344 ha, equivalent to just one of 18 landscapes (10,000 ha) in our study.
Landscape scale factors related to amount and spatial configuration of forest were more important than other cover types, indicating bird species in our sample responded to patterns associated to forest habitat loss and fragmentation.Although many published studies highlight the quality of the matrix when determining composition and abundance of species within patches (Laurance 1991, Gascon et al. 1999, Kotze and Samways 1999, Cook et al. 2002), our results suggested the surrounding matrix exhibited a limited influence on bird abundance.A possible explanation could be that many of the selected bird species use forest patches both for feeding and nesting, which could limit their use of the matrix.On the other hand, the movements and resource selection patterns of species using both forest patches and surrounding matrix are unknown for the Espinal, so changes in land use, i.e. soybeans one year, followed by corn, may not influence avian distribution patterns.Also, bird species that actively forage in the matrix, e.g., Picui Ground Dove, may have access to food equally provided by different land use types, resulting in weak responses to quantity and spatial arrangement of crops that characterized the matrix in our study.Finally, there may be additional variables that were not incorporated or the manner in which some of the incorporated variables were measured may have failed to capture any potential influences on the bird community.Some species in our study, e.g., Scimitar-billed Woodcreeper, Narrow-billed Woodcreeper, Brown Cacholote, Creamy-bellied Thrush, and Red-crested Cardinal, were more abundant in landscapes dominated by closely arranged, interconnected forest patches.This may reflect the benefit of increased landscape connectivity and access to available patch resources (Andrade andMarini 2001, Lees andPeres 2008).On the other hand, Picui Ground Dove, Grayish Baywing, Great Antshrike, and Saffron Finch were more abundant in landscapes characterized by greater forest fragmentation.These species are known generalists commonly found in modified agricultural landscapes.Their preference for landscapes characterized by isolated forest patches likely reflects their greater ecological plasticity in resource use (Fraser andStutchbury 2004, Thornton et al. 2011).Our twoscale approach to assess bird responses to environmental characteristics highlighted the differential use of landscape elements by a representative sample of birds in the Espinal region.Hence, negative impacts of loss and fragmentation of native forest may not necessarily be offset by local measures alone, such as retaining large trees.Landscape scale measures describing structure, amount, and spatial configuration of forest should be considered to ensure the regional persistence of bird species.This information contributes reliable knowledge to integrate in the design of conservation measures for agricultural landscapes in the Espinal region of Argentina.

CONCLUSIONS
In the absence of land use regulations agricultural expansion will continue in the province of Entre Ríos, with increasing alteration of native forest cover, spatial configuration, and structure within patches.Our results provide useful and reliable information for the development of management and land conservation plans.As such, these may serve to provide guidance on landscape conservation design for the fragmented agroecosystems of the region.At the patch scale, retaining larger trees in forest patches will help maintain the abundance of forest species while allowing agricultural activities inside regular lots.At the landscape scale, limiting the conversion of forests to agricultural lands, avoiding further fragmentation of forest patches, and designating areas where large forest fragments are protected and maintained will benefit forest specialist species.To prevent major negative effects on biodiversity, additional conservation efforts are necessary.Some of these are currently in progress, including the National Law to protect Argentine forests.However, its effectiveness will depend on decision makers and land owners given the lack of control on land use policies.If soybean production and cattle ranching continue to expand across the Espinal forest ecosystem without land use planning, the negative effects of these activities on avian biodiversity will likely continue to increase.
Responses to this article can be read online at: http://www.ace-eco.org/issues/responses.php/1222Table A4.1 Set of plausible models describing the relationship of within patch, patch and landscape variables with bird abundance for each season.Models are shown with their AICc values and Akaike weight (ωi).Best set of models were those having ΔAICc ≤ 2 (Burnham and Anderson 2002)

Fig. 1 .
Fig. 1.Map of Argentina indicating location of Entre Ríos province and study mosaics.Detail of sublandscape with highlighted forest patch (dark green areas) and outline of bird sampling points within forest patch.

Table 1 .
Correlation matrix for forest patch scale variables.D_Arb = Density of trees, D_arb = Density of shrubs, A_Arb = Tree height, DBH = Diameter at breast height, PA = Forest patch area, TE = Total edge, SI = Shape index, PARA = Perimeter-area ratio, FD = Fractal dimension, DistLoc = Lineal distance from forest patch to nearest settlement.

Table 2 .
Correlation matrix of spatial metrics for amount and spatial configuration of forest cover of sublandscapes (10 x 10 km).FA = Total forest area (ha), NP = Number of patches, ED = Edge density, PA = Mean patch area, SI = Shape index, FD = Fractal dimension, PARA = Area-perimeter ratio, ENN = Mean euclidean nearest neighbor patch distance, COH = Patch cohesion index, AI = Aggregation index.

Table 3 .
Best generalized linear models describing bird abundance in forest patches during the austral fall-winter and spring summer, in the province of Entre Ríos, Argentina.Parentheses indicate a negative relationship to bird abundance.
† Codes of bird species listed in TableA1.1.
. ΔAICc= AICi -AICmin values.Parentheses indicate a negative relationship to bird abundance.Codes of bird species listed in Table A1.1.
Table A5.1 Mosaic description based on patch and landscape variables during 2007 and 2008 in the province of Entre Ríos, Argentina.Values are mean ± SE in (). and b indicate statistical differences of variables between mosaics.Means with a common letter are not significantly different (P > 0.05) by Kruskal Wallis test. a