Quantifying avian relative abundance and ecosystem service value to identify conservation opportunities in the Midwestern U . S .

Grassland birds are highly imperiled because of historical habitat loss and ongoing conversion of grasslands to agricultural and urban land uses. Therefore, prioritizing and further justifying conservation action in remaining grasslands is critical to protecting what remains. Grassland bird conservation has focused on identifying and protecting large grassland complexes referred to as Grassland Bird Conservation Areas (GBCAs). We identified and classified GBCAs in a region highly impacted by both agricultural and urban land conversion using previously developed methods. Then, we extended the analysis to include estimated relative abundance of five grassland focal species in each GBCA. Models of relative abundance were built using eight years of monitoring data collected by citizen scientists. Finally, we quantified the value of ecosystem services provided by each GBCA. There were nearly 55,000 ha of grassland habitats in the Chicago Metropolitan Region that met GBCA criteria, 33% (18,415 ha) of which were protected. Proportion of abundance in protected versus unprotected areas was similar for Bobolink (Dolichonyx oryzivorus; 46%), Grasshopper Sparrow (Ammodramus savannarum; 52%), and Sedge Wren (Cistothorus platensis; 48%), whereas, Henslow’s Sparrow (Ammodramus henslowii; 75%) had a higher proportion of relative abundance in protected GBCAs and Eastern Meadowlark (Sturnella magna) had lower proportions (37%). GBCAs provisioned just under $900 million annually in ecosystem services, 73% of which is because of flood control. Outputs of this comprehensive approach will inform grassland bird conservation by providing detailed information about the value for birds and people of grassland habitats. Utilisation de la mesure de l'abondance relative aviaire et de la valeur des écoservices pour déterminer les occasions de conservation dans le Midwest étasunien RÉSUMÉ. Les oiseaux de prairie sont grandement menacés en raison des pertes d'habitat subies dans le passé et de la conversion continuelle des prairies à des fins d'agriculture ou d'urbanisation. Ainsi, la priorisation et la justification des mesures de conservation des prairies restantes sont critiques si on veut les protéger. Les mesures actuelles de conservation des oiseaux de prairie reposent surtout sur l'identification et la protection de vastes complexes de prairies, aussi nommés Grassland Bird Conservation Areas (GBCA). Nous avons identifié et classé les GBCA d'une région fortement touchée par la conversion des terres à des fins d'agriculture et d'urbanisation au moyen de méthodes élaborées précédemment. Nous avons ensuite élargi l'analyse pour inclure l'abondance relative de cinq espèces d'oiseaux de prairie estimée dans chaque GBCA. Les modèles d'abondance relative ont été construits à partir de huit années de données récoltées par des citoyens participant à des projets scientifiques. Finalement, nous avons mesuré la valeur des écoservices rendus par chaque GBCA. Près de 55 000 ha de milieux de prairie dans la région métropolitaine de Chicago correspondaient aux critères d'un GBCA, dont 33 % (18 415 ha) étaient protégés. L'abondance dans les régions protégées se comparait à celle dans les régions non protégées pour le Goglu des prés (Dolichonyx oryzivorus; 46 %), le Bruant sauterelle (Ammodramus savannarum; 52 %) et le Troglodyte à bec court (Cistothorus platensis; 48 %), tandis qu'elle était plus élevée pour le Bruant de Henslow (Ammodramus henslowii; 75 %) et plus faible pour la Sturnelle des prés (Sturnella magna; 37 %) dans les GBCA protégés. Les écoservices fournis par les GBCA se chiffraient tout juste sous la barre des 900 millions de dollars, dont 73 % dans la régulation des inondations. Les résultats de cette approche compréhensive peuvent contribuer à la conservation des oiseaux de prairie en fournissant des renseignements détaillés sur la valeur qu'ont les milieux de prairie sur les oiseaux et les humains.


INTRODUCTION
The United States has lost 98% of its original tallgrass prairie as the result of urban development and vast agricultural expansion (Samson and Knoff 1994), and conversion of grasslands in general continues to outpace restoration (Lark et al. 2015).Grassland bird populations have been among the fastest and most consistently declining suite of species in North America over the past 40 years.Of the 46 North American grassland-breeding birds, 55% have shown significant declines and 48% are of conservation concern including four that are federally endangered in the United States (NABCI 2009).An indicator for grassland birds based on 24 grassland-obligate species dropped 40% from the baseline value set in 1968 (NABCI 2009).In recent years this indicator has stabilized, albeit at historically low levels (NABCI 2014), but recent increases in grassland conversion within the contiguous U.S. to row crops (1-5% annually; Wright and Wimberley 2013), primarily corn (Lark et al. 2015), may lead to further declines.
Grassland bird conservation has included the protection of large tracts and complexes of grassland habitat referred to as Grassland Bird Conservation Areas (GBCAs; e.g., Johnson et al. 2010).This framework was developed to prioritize conservation given the sensitivity of many grassland birds to habitat fragmentation (Herkert 1994, Herkert et al. 2003, but see Benson et al. 2013), but GBCA designation does not by itself equate with any level of protection.Designations of GBCAs have been based on expert opinion (CEC and TNC 2005) as well as quantitative geospatial criteria (Sample andMossman 1997, Johnson et al. 2010).GBCAs include a core area and surrounding matrix and are further classified into types (I-IV).The most valuable GBCAs (Type I) have the largest core area and proportion of grassland in the surrounding matrix.Core patches and their surrounding matrix can be thought of as conservation or management units, for example, when a core is a protected area.GBCAs have been identified in the Prairie Pothole Region of the U.S. (Johnson et al. 2010), and in the states of Missouri and Iowa (Jacobs et al. 2005(Jacobs et al. , 2012)).Surprisingly, bird observations are not used to identify or rank GBCAs because the identification process occurs at scales larger than available observational data.Thus, using empirical data to test whether GBCAs are effective in conserving grassland birds is necessary to ensure they effectively contribute to prioritization of conservation actions (Howard et al. 2014, Veloz et al. 2015).
In addition to wildlife conservation, grasslands provision valuable ecosystem services.Ecosystem services, or the benefits provided to people by natural resources, enhance the quality of life for people worldwide (Birch et al. 2010, Bateman et al. 2013, Arkema et al. 2015, Ferraro et al. 2015).For example, forest restoration has been found to increase carbon storage and tourism in Latin America and Southeast Asia (Birch et al. 2010, Ferraro et al. 2015); coastal protection has increased fishery and recreational revenues while decreasing damages caused by storms in Latin America and the U.S. (Arkema et al. 2015, Reddy et al. 2016); and targeted land use planning has increased urban green space and bird species diversity in the United Kingdom (Bateman et al. 2013).
Quantifying ecosystem services across a landscape is increasingly used as a method to demonstrate the economic value associated with the protection of natural areas (Polasky et al. 2011, Arkema et al. 2015, Hamel et al. 2015).Additional studies have illustrated the link between ecosystem services and biodiversity conservation (Nelson et al. 2009, Chaplin-Kramer et al. 2015, Reuchlin-Hugenholtz and McKenzie 2015).It can be valuable for land managers to incorporate the benefits provided by ecosystem services when making decisions regarding land use (Polasky et al. 2011).
We (1) developed regional models of relative abundance for grassland birds based on bird survey data and environmental measures, (2) used those models and landscape configuration attributes to designate GBCAs for regional populations, and (3) estimated the ecosystem services provided by this network of GBCAs.This analysis aimed to inform land management decision making in the Chicago Metropolitan Region by identifying currently unprotected areas that are valuable for both grassland birds and ecosystem services.The work used a multiyear systematic survey of avian abundance collected by citizen scientists to build models of avian relative abundance for five grassland bird species covering a range of habitat requirements: Bobolink (Dolichonyx oryzivorus), Sedge Wren (Cistothorus platensis), Henslow's Sparrow (Ammodramus henslowii), Eastern Meadowlark (Sturnella magna), and Grasshopper Sparrow (Ammodramus savannarum).We estimated relative abundance based on landcover composition and configuration, soils, and vegetation condition.We then classified existing grasslands into four GBCA types and summarized estimates of relative abundance and ecosystem service value from the Chicago Wilderness Green Infrastructure Vision (GIV, described below) for each GBCA.

METHODS
The project area was defined as the Chicago Metropolitan Region, a 14,615 km² region including seven metropolitan counties (McHenry, Lake, Kane, DuPage, Cook, Kendall, and Will) surrounding Chicago, Illinois, USA.The regional population is 3.09 million people with roughly 1.04 million living in urban Chicago (U.S. Census Bureau 2010).The study area is a mix of high-density urban, suburban, and rural agricultural land uses.

Observation data
We built models of relative abundance for five obligate grassland bird species: Bobolink, Eastern Meadowlark, Henslow's Sparrow, Grasshopper Sparrow, and Sedge Wren.Annual five-minute, fixed-radius (75 m) point counts conducted between 2007 and 2014 were used to estimate avian abundance.All birds seen and heard were counted.Surveys were conducted in the morning (before 09:00) in the month of June.The majority of observations were collected by Bird Conservation Network (http://www.bcnbirds.org)volunteers.Point count data collected by the Lake County Forest Preserves, the Illinois Department of Natural Resources (Will County, IL), and the Illinois Natural History Survey's Critical Trends Assessment Program were also included, but constituted a small proportion (< 1%) of available data.These were 10-minute counts collected in time increments, and only the first five minutes were used.All surveys were repeated annually at the same point locations.Points were separated by 150 m and revisited in subsequent years.To reduce sampling bias in the dataset that could lead to a biased abundance model (Boria et al. 2014), we selected only one observation from each sampling location by choosing the year with the maximum abundance (selecting randomly between ties) to minimize the impact of imperfect detection by the volunteer surveyors.Maximum abundance has been shown to more closely reflect true abundance (Toms et al. 2006).We then sampled annual environmental data for the year in which each filtered observation was made.

Environmental data
Environmental datasets used to predict the relative abundance of the five focal species included landcover, soils, and vegetative condition.We extracted annual landcover from 2007-2014 for the study area from the National Agricultural Statistics Service (USDA) CropScape web service (http://www.cropscape.com;Han et al. 2012).
The 30-meter resolution dataset classifies cropland on an annual basis and utilizes National Landcover Database (Homer et al. 2015) classes for nonagricultural areas.We consolidated cropland and landcover classes into 12 classes (Table 1).Landcover data were summarized in multiple size and configuration metrics at each survey point, including the landcover class and patch size, proportion of each consolidated landcover class within one kilometer, and distance to other landcover classes.Also characterized at each point were annual vegetation condition, i.e., greenness measured as normalized difference vegetation index (NDVI), in mid-June, extracted from the VegScape web service (Han et al. 2012), and soil hydrologic group, based on the gridded SSURGO soils dataset (NRCS-USDA [date unknown]).The resulting environmental dataset included 33 predictor variables.Survey time and date were also included as predictors, resulting in a total of 35 variables.All environmental data was summarized at each point location using the raster package (Hijmans 2015) in program R (R Core Team 2014).

Abundance modeling
Estimating absolute densities or population sizes from point counts is a challenge (Efford and Dawson 2009).The best available methods require distance estimates or time-of-detection data for each observation, or multiple-observers per count (Farnsworth et al. 2005, Johnson 2008).This information is then used to estimate species-specific detection probabilities that are used to adjust estimated abundances for imperfect detection.Low detection is less of an issue for highly vocal species, such as Western Meadowlark, but can lead to underestimation of density in species that sing from the ground, such as the Grasshopper Sparrow (Leston et al. 2015).
The point-count data available for this analysis lacked distance or time-of-detection information.Therefore, we opted to model relative abundance across sites instead of true population size.
We generated models of relative abundance for each species using Boosted Regression Trees (BRTs).BRTs are a machine learning algorithm with high predictive performance able to capture nonlinear relationships between environmental and response variables and interactions among predictor variables (Elith et al. 2008).We built models based on techniques outlined in Elith et al. (2008) using the gbm.step function in the package dismo (Hijmans et al. 2015).We used a Poisson model for count data and set tree complexity to five.Learning rate was adjusted (range: 0.005-0.01)such that at least 1000 trees were used to create each final model (Elith et al. 2008).We used a stepwise variable removal procedure (Elith et al. 2008) to reduce the 33 environmental variables to a smaller set for each species.The stepwise procedure was run on the full dataset using ten-fold cross-validation to assess model performance.We eliminated all predictor variables that showed no sign of improving predictive performance by plotting performance against the number of variables included in the model and identifying the most parsimonious variable set that also maximized performance.The reduced variable set was used in subsequent analyses.
We modeled avian mean relative abundance as the performanceweighted average of 20 BRT models generated for each species (as in Barker et al. 2014).We assessed the performance of each model by cross-validation using spatially stratified subsampling.Spatially stratified subsampling generates robust performance measures because model training and test data are spatially independent (Bahn and McGill 2013).We used semivariograms to explore spatial autocorrelation in our environmental dataset and determined that on average spatial autocorrelation leveled off at distances greater than five kilometers.We then created a novel procedure for building model training and test datasets in which a 5-km resolution grid was overlaid on the study area and all sampling points within randomly selected 5-km grid cells were assigned to either the training or test dataset at a ratio of 5:

Relative abundance prediction
We predicted avian relative abundance as a continuous surface throughout the study area using environmental data from 2014.
Relative abundance was estimated as a performance-weighted average of the 20 prediction grids using the Pearson's correlation as the performance measure.We scaled relative abundance to the area of each 30-m grid cell by dividing by the area of each 75-m point count radius and multiplying by the area of each grid cell.

Landscape analysis
We identified patches of suitable grassland habitat appropriate for landscape-scale conservation from modeled abundances (described above) and landcover.We generated patches of suitable habitat from areas classified as herbaceous grasslands, wheat, other forest and shrub, herbaceous wetlands, and woody wetlands in the Cropland Data Layer (Han et al. 2014) by clumping adjacent grid cells into patches.Most grassland birds avoid woody habitat.However, the other forest and shrub and woody wetlands classes were included because comparisons with locally delineated vegetation maps suggested their classification accuracies were low and because they represented a relatively low proportion (< 5%) of the landscape.Furthermore, these areas are often converted to grasslands through land management and therefore represent a conservation opportunity.
We used criteria for minimum patch size and modeled bird abundance to further refine the patch list.Many grassland birds show evidence of being area sensitive (Herkert 1994, Johnson and Igl 2001, Herkert et al. 2003, Ribic et al. 2009).Yet, recent studies synthesizing information across larger geographic scales suggest that area sensitivity may differ across regions (Ribic et al. 2009).Benson et al. (2013) found little relationship between patch size and nest success in the Midwestern U.S. Walk et al. (2010) demonstrated the importance of small grassland patches for some nesting grassland birds, such as the Eastern Meadowlark.Given the highly fragmented nature of grassland habitats in the study area, we chose to include minimum patch size of 12 ha to avoid overlooking smaller sites that may still be of use to grassland birds.We then removed those patches with modeled relative abundances failing to meet a minimum threshold specific to each species, based on estimated densities from distance sampling conducted previously at a subset of sites.We determined these thresholds by taking the estimated average density found by Buxton (2014) and converting it to the number of birds estimated to be within an area of 12 ha to be consistent with our minimum patch size criterion.We then applied this minimum threshold to all patches, regardless of size.For example, Buxton (2014) found an average density of about 0.95 Bobolink per ha, which translates to a threshold of 11.4 Bobolink across a 12-ha area.We then only included patches with at least 11.4 Bobolink, or another species' minimum threshold.Our rationale assumed that small patches of below-average density were of lower conservation priority and that larger patches supporting this minimum population size were of interest.Thus, the criteria for including small patches was high, while the criteria for including large patches was low.We assumed that large patches with small estimated populations still represented a conservation interest because they could be enhanced through management.

Compiling ecosystem services
Green infrastructure refers to the network of natural areas that not only support native species but also provide ecosystem services such as water filtration, flood control, and carbon storage.
We This study focused on four ecosystem services: flood control, groundwater recharge, water purification, and carbon storage (Table 2).We summarized these particular services because these are the only ecosystem services having valuation grids for the region prepared by The Conservation Fund and thus can be quantified using the GIV Summary Tool.
The valuation of these ecosystem services considered several factors, some of which included avoided cost (e.g., from property damages or public health problems that would have occurred in the absence of those services), replacement cost (i.e., the cost of man-made systems that would emulate those services), and factor income (e.g., improved water quality increases commercial fisheries, which increase fishermen incomes; Farber et al. 2002).

Protected areas analysis
We determined how much of each patch is currently protected and unprotected by overlaying each set of patches with a layer representing protected lands (sourced from The Field Museum in Chicago) in the region.We calculated the proportion that is protected and unprotected for each patch and multiplied these proportions by the patch area, bird relative count, and ecosystem service values.This enabled us to identify the largest patches of unprotected habitat and their corresponding values for relative bird abundance and ecosystem services.

Grassland Bird Conservation Area analysis
We classified patches of suitable habitat into GBCAs using criteria modified from Johnson et al. (2010).Classifications were based on patch size and proportion of grassland cover in the matrix surrounding the patch (Table 3).We added a fourth GBCA type to account for the smaller patches found in the Chicago metropolitan region's highly fragmented landscape.Each patch was classified as a Type I, II, III, or IV core based on the area and width of the patch.Next, we calculated the proportion of compatible habitat encompassed by each patch and its matrix, the area enclosed by a 1.6-km buffer surrounding the patch.Type level was assigned based on the lowest criteria met for either patch size or proportion of grassland cover.For example, a 300+ ha patch with only 32% grassland within its matrix would be classified as a Type II GBCA.Modeled relative abundance and ecosystem services valuations were then aggregated by GBCA types and protection status.These included an affinity for high proportions of grass/hay and avoidance of croplands at the landscape scale across species, and preferences for hydric soil types and proximity to water for some species (e.g., Sedge Wren).
Nearly 55,000 ha of grassland habitats in the Chicago Metropolitan Region met GBCA criteria (Fig. A1.3); 33% (18,415 ha) of which are protected (Fig. 1).The majority of habitat was classified as Type III and IV GBCAs.Type I GBCAs, the largest areas surrounded by high proportions of grasslands, represent the smallest proportion of GBCA-assigned habitat in the region (3.5%).A Type I GBCA was more likely to be protected than not, whereas Type III and IV GBCAs were more likely to be unprotected and Type II was approximately half protected (Fig. 1).Protected GBCAs host relatively fewer Eastern Meadowlark (37%) than unprotected GBCAs, whereas Henslow's Sparrow (75%) had higher relative abundance in protected GBCAs.Proportions of abundance in protected versus unprotected areas were similar for Bobolink (46%), Grasshopper Sparrow (52%), and Sedge Wren (48%).Across species, a bird located in a Type I or II GBCA was more likely to be on protected land than one found on a Type III or IV GBCA.Area-adjusted indices of relative abundance (i.e., a density index) showed higher per-ha densities in protected than unprotected lands.Density indices were similar across GBCAs for the most abundant species: Bobolink and Eastern Meadowlark.Density indices for less common species were higher in Type I and II GBCAs.
Grassland GBCAs provisioned just under $900 million annually in ecosystem services based on GIV valuations (Table 4).The majority (73%) of that value came from flood control followed by groundwater recharge (18%), water purification (9%), and carbon sequestration (< 1%).As with avian abundance, ecosystem services provisioned in Type I and II GBCAs are more likely to be from protected lands than those in Types III and IV.Fifty-five percent of the total valuation for flood control was protected, 42% for groundwater recharge, 47% for water purification, and 46% for carbon sequestration.   3 for definitions of GBCA types.White areas represent incompatible habitat.

DISCUSSION
The Chicago Metropolitan Region hosts four Type I GBCAs, 91 Type II GBCAs, and hundreds of Type III and IV GBCAs (Fig. 2).Nearly 33% of these areas are protected.These GBCAs provided important habitat for all five grassland birds modeled.The value of these sites extends beyond their support of wildlife populations to include nearly $900 million in ecosystem services.
The most valuable of these is flood control, comprising 73% of the total valuation.Estimates of avian relative abundance and ecosystem services demonstrate the value of these natural habitats and provide justification for protecting currently unprotected areas in the region as well as habitat restoration.
Quantifying value is an initial step in building a robust regional grassland conservation plan.The outputs summarized above characterize the importance of grassland habitats to birds and people in the Chicago Metropolitan Region.On average across species and ecosystem services, less than half of relative abundances and ecosystem services value are currently protected.
Remaining unprotected lands are therefore a high-value conservation opportunity.Prioritizing among habitat patches could be done through a spatial optimization routine, treating birds and ecosystem services as conservation objectives (e.g., Ban et al. 2009).Alternatively, priorities could be set through coordination across local land management agencies.Currently, the National Audubon Society is working with county land managers to develop a regional grassland conservation vision based on these data.3 for definitions of GBCA types.
One unexpected result was the high proportion of relative abundance associated with Type III and IV GBCAs.For the most part, this reflects differences in the total area; Type III and IV GBCAs are more than twice the area than Type I and II (Fig. 2).Furthermore, density indices are lower in Type III and IV GBCAs.
Thus, the return on investment, which is often used to prioritize conservation decisions (Withey et al. 2012), for Type III and IV GBCAs may be lower, at least in terms of relative bird abundance for less common species of high conservation need.
One conceptual model for GBCAs in the northern tallgrass prairie targeted a core patch size of 800 ha surrounded by an additional 1000 ha of grasslands within a 1.6-km buffer (Fitzgerald et al. 1998).Multiple, adjacent Type I GBCAs located within the Midewin National Tallgrass Prairie when combined, would reach this goal.For example, the largest single patch in our analysis (682 ha), also located in Midewin, is surrounded by 2300 ha of grasslands within its 1.6-km matrix buffer.Thus, the Midewin grassland effectively meets the criteria of a premier GBCA.Most (54%) Type I and II GBCAs are already protected.Therefore, increasing connectivity among these areas, and protecting unprotected Type III and IV GBCAs (73% of which are unprotected) could create additional 1000-ha grasslands with the required matrix of surrounding grasslands.Current land-use change trends are in the opposite direction (Lark et al. 2015), but this could change locally or regionally through federal farm policy.The goal of protected large GBCAs in the Chicago Metropolitan Region is attainable and is desirable for grassland bird conservation and the added benefits of provisioned ecosystem services.
We utilized citizen science data to characterize relative abundance across grassland patches.The use of citizen science data to build models of bird-habitat relationships is relatively new.Comparisons of models built using citizen-collected data versus structured surveys conducted by trained professionals suggest that citizen science is a useful proxy (Munson et al. 2010, Callaghan andGawlik 2015).We used standard-effort point counts to model relative abundance.In a separate analysis, we converted abundance outputs to densities and compared them to site-level density estimates based on distance sampling for a handful of grasslands in the study area (Buxton 2014).We found our estimates biased high for abundant species (Bobolink and Eastern Meadowlark) and biased low for rare species (Grasshopper Sparrow, Henslow's Sparrow, and Sedge Wren).A bias high for abundant species may suggest a tendency of volunteers to overestimate counts of highly detectable and abundant species.For example, loud vocalizations (e.g., Eastern Meadowlark) may have led observers to include individuals that were > 75 m away.A bias low for less common species suggests that low detection probabilities are an issue.This is further supported by the presence of detection variables, such as time of day and survey date, in final models (Appendix Figs.A1.4a-e).Therefore, modeling detection probability would likely improve estimates of relative abundance for less common species.In spite of these biases, models of relative abundance did reflect known species-habitat relationships and is therefore a useful measure of relative value across management units (i.e., GBCAs).
Modeled abundance-habitat relationships were consistent with the literature.As expected, proportion of grass and hay was a top variable for all five species.Previous study has found evidence of sensitivity to decreasing grassland patch size for Bobolink, Eastern Meadowlark, Grasshopper Sparrow, and Henslow's Sparrow (Herkert 1994).Patch size was retained in models for these four species.Drivers of abundance among Bobolink, Eastern Meadowlark, Grasshopper Sparrow, and Henslow's Sparrow were similar.All included some measure of grass/hay, development, and cropland.Meanwhile, Sedge Wren stands out from the group with key predictors relating to water, soils, and variables associated with detection.Previous study suggests that vegetation characteristics drive Sedge Wren occurrence (Herkert 1994).Low performance of the Sedge Wren model, in particular, may reflect the absence of information on local vegetation.Vegetation condition, as measured by NDVI, was not in the top five predictors for any species.NDVI has been used elsewhere to model avian-habitat relationships (Gottschalk et al. 2005), although additional information found in raw reflectance data may improve performance over NDVI (Shirley et al. 2013).Here, landscape composition provided greater information than NDVI for predicting avian abundance.We explored modeling with raw reflectance data in preliminary analyses, but the models underperformed relative to models based on classified landcover and the availability of cloud-free images during the breeding season reduced the available dataset.
The use of classified landcover based on remote sensing in predictive modeling as compared to local vegetation measures introduces additional errors in abundance models.The Cropland Data Layer (Han et al. 2014) has high classification accuracies for agriculture.For example, classification accuracy for corn is 98%.However, natural vegetation classifications are based on the National Landcover Database (Homer et al. 2015), for which accuracy of grasslands and haylands is likely low.Unfortunately, no formal accuracy assessment exists for the National Landcover Database.Refining local landcover data and augmenting it with information on grassland vegetation composition and management, e.g., enrollment in the Conservation Reserve Program, would likely improve model performance.
Ecosystem service valuations from the Chicago Wilderness GIV estimate that grassland habitats provide hundreds of millions of dollars of services annually, largely for flood control.Ecosystem service valuations are helpful for communicating the benefits to society from protecting natural areas (Guerry et al. 2015); however, valuations are still an estimate and have associated errors.In particular, an approach called "benefit transfer" takes a valuation from one location and transfers it to another location with similar landcover (Plummer 2009), but the validity of the transfer is never verified.This characterizes the ecosystem service valuation used here and may lead to overestimation of monetary values if estimates are made in areas of high known value, e.g., along waterways, and transferred to upland areas elsewhere.Furthermore, the dependence of this method on landcover means that ecosystem service valuations are subject to the same underlying errors of the landcover datasets.
Our valuation of ecosystem services associated with grassland bird habitat is limited to four particular ecosystem services, but birds themselves are known to provide ecosystem services (Wenny et al. 2011), such as seed dispersal (Howe and Smallwood 1982, Nathan and Muller-Landau 2000, Cordeiro and Howe 2003), pollination (Feinsinger et al. 1982, Kelly et al. 2004), or pest control for adjacent agriculture (Mols and Visser 2002, Van Bael et al. 2003, Perfecto et al. 2004).Future studies quantifying the monetary value associated with such services will generate further support for GBCAs as highly valuable habitats worthy of protection for the benefit of birds and people.
Responses to this article can be read online at: http://www.ace-eco.org/issues/responses.php/902

Fig. 1 .
Fig. 1.Type I-IV Grassland Bird Conservation Areas (GBCAs) in the Chicago Metropolitan Region.Type I GBCAs cover a total of 1933 ha; Type II GBCAs cover 11,938 ha; Type III GBCAs cover 25,343 ha; and Type IV GBCAs cover 15,778 ha.See Table3for definitions of GBCA types.White areas represent incompatible habitat.

Fig. 2 .
Fig. 2. Hectares of protected and unprotected grassland habitat classified as Type I-IV Grassland Bird Conservation Areas (GBCAs).See Table3for definitions of GBCA types.

Table 1 .
Environmental data used to predict relative avian abundance including consolidated landcover classes and derived variables, vegetative condition, and soil data.
(Weber 2014)he Chicago Metropolitan Agency for Planning's (CMAP) GIV Summary Tool (available at http://datahub.cmap.illinois.gov)tocalculateecosystemservice values for each patch that met the habitat, minimum area, and minimum abundance criteria.This tool uses dollar values (based on existing literature regarding ecosystem services) estimated by The Conservation Fund, which used studies based in the Chicago Metropolitan region when possible(Weber 2014).For example,Pimentel et al.
(1997)estimated a water quality value of $118.30associated with each hectare of grassland per year (2006; $138.91/ha/year in 2014).Also, Costanza et al. (1997) estimated a gas and climate regulation value of $9.50 associated with each hectare of grassland per year (2006; $11.16/ha/year in 2014).

Table 2 .
List of ecosystem services included in valuation by Green Infrastructure Vision (GIV) 2.3(Allen et al. 2014).
decreasing performance.Model performance was a function of the proportion of nonzero counts (i.e., prevalence) for each modeled species (F = 20.11,p < 0.0001), suggesting that the paucity of observations available for some species limited our ability to model relative abundance.Overall model performance was low to moderate across species; however, variable response curves suggest that underlying processes are being captured.

Table 4 .
Total ecosystem service value ($2014/ac/year) associated with both protected and unprotected portions of each Grassland Bird Conservation Area type.