Skip to main content
  • Research article
  • Open access
  • Published:

The luxury effect beyond cities: bats respond to socioeconomic variation across landscapes

Abstract

Background

The luxury effect describes the positive relationship between affluence and organism diversity or activity in urban ecosystems. Driven by human activities, the luxury effect can potentially be found at a broader scale across different landscapes. Previously, the luxury effect relationship has been established within a city for two bat species, the red bat (Lasiurus borealis) and the evening bat (Nycticeius humeralis). We examined landscape-scale patterns of bat activity distribution—using empirical data for seven bat species for the luxury effect. We also identified bat-land cover associations for each species. Across North Carolina, USA, we used the mobile transect survey protocol of the North American Monitoring Program to record bat activity at 43 sites from 2015 to 2018. We collected land cover and income data at our transect locations to construct generalized linear mixed models to identify bat-land cover and bat-income relationships.

Results

We found that across landscapes, activity of the red bat and the evening bat was positively correlated to income independent of land cover, consistent with previous single-city results. We found a negative relationship between hoary bat (Lasiurus cinereus) activity and income. All seven species had specific land cover associations. Additionally, we found a positive interaction term between income and evergreen forest for the red bat and a positive interaction term between income and woody wetland for hoary bat.

Conclusions

Our results demonstrated that the luxury effect is an ecological pattern that can be found at a broad spatial scale across different landscapes. We highlight the need for multi-scale ecology studies to identify the mechanism(s) underlying the luxury effect and that the luxury effect could cause inequity in how people receive the ecosystem services provided by bats.

Background

The term “luxury effect” describes a functional relationship that links affluence and biodiversity in urban ecosystems [1]. Hope et al. [1] suggested that in urban ecosystems the traditional resource-diversity relationship has been modified by human activities. Instead of natural resources being limiting factors, financial resources that interact with land use, legacy effects, and other sociocultural factors shape the biodiversity patterns in urban areas [1]. The luxury effect has been demonstrated worldwide across different plant and animal species in many urban areas (e.g. [2,3,4,5,6]).

There are several factors that could potentially drive the luxury effect in an urban area (see [7] for a comprehensive review). Affluence, generally measured by median household income, predicts where people live and how people manage their property. These human activities directly affect the physical structure, microclimate, and vegetation of different parts of an urban area [8,9,10,11], which can consequently impact animal distributions [12,13,14]. Human activities are also influenced by the existing environmental heterogeneity and inequity in an urban area and may amplify the heterogeneity and inequity over time [15,16,17,18]. Therefore, the core driver of the luxury effect is the difference in human activities caused by the socioeconomic differences.

The luxury effect is found in cities likely because urban ecosystems have the highest human density and are highly impacted by human activities [19]. As the world is urbanizing, it is important to examine whether the luxury effect of urban ecosystems could be found in other landscapes and whether it is evident at a broader scale. At the global scale, studies have shown that terrestrial biodiversity hotspots are in low-income developing countries [20,21,22]. However, many biodiversity hotspots such as the Amazon rainforest or the south-central region of China have very limited human presence and still contain primary vegetation [20]. In contrast, there has been no empirical study to examine the luxury effect in landscapes that all have been heavily modified by humans and compare the magnitude of this effect across landscapes.

Differences in human activities caused by socioeconomic differences can be found in landscapes other than the urban landscape. In the US, the socioeconomic background of farmers influences whether they would participate in conservation-compatible practices and programs on farmlands [23]. In Australian agroecosystems, whether farmers participate in conservation schemes for farming practices is influenced by complex social factors [24]. Piemonti et al. [25] demonstrated that watershed conservation and management actions could be influenced by landowners’ socioeconomic background. Even in protected areas, various land use activities are driven by socioeconomic factors [26], such as exurban housing developments in areas close to nature for the wealthy [27, 28] or cattle grazing on mesic habitats driven by rent economics [29]. Therefore, it is reasonable to expect that the luxury effect in urban ecosystems can scale up spatially and be found across landscapes.

Previously, a single-city study of bat distribution patterns indicated bats responded to urban affluence variations. Li and Wilkins [4] found that in a medium-sized city (Waco, Texas, USA), two tree-roosting bat species, the eastern red bat (Lasiurus borealis) and the evening bat (Nycticeius humeralis) had a higher probability of presence in neighborhoods with higher median household income. Further, in this study these two species responded to vegetation coverage, which was consistent with their habitat preference [30]. However, contradictory to studies that examine urban vegetation and affluence (e.g. [8, 31]), vegetation coverage was not correlated to median household income in Waco. The independence of vegetation coverage and median household income, and that both factors affect bat distribution [4], suggests that the luxury effect is more complicated than a simple altering of habitat availability [1].

Bats are able to explore large areas and are associated with different land cover types [32, 33]. In France, a landscape scale study showed that four bat species avoid intensive agriculture whereas the effect of urban land cover varies by species [34]. In the USA, two bat species positively respond to urban land cover [35]. In Italy, the lesser horseshoe bat (Rhinolophus hipposideros) distribution at the landscape scale is associated broad-leaved forests [36]. In addition to land cover types, socioeconomic differences among landscapes elements could potentially be perceived by bats and affect distribution.

The objective of our study was to examine if the luxury effect could be found in bats at the landscape scale using empirical data collected across multiple urban centers. We also aimed to identify species-specific bat-land cover associations and test whether the luxury effect and bat-landcover association were two independent processes. Because there are only a few bat species in most parts of the US and the occurrence of these species is common, alpha diversity is less informative than the activity of bats [7]. Hence, we selected seven common North American bat species as the target to examine the luxury effect in species-specific bat activity. The target species were the big brown (Eptesicus fuscus), eastern red (L. borealis), hoary (L. cinereus), silver haired (Lasionycteris noctivagans), evening (Nycticeius humeralis), tricolored (Perimyotis subflavus), and Mexican free-tailed (Tadarida brasiliensis) bat. Specifically, we hypothesized that across different landscapes, the luxury effect would be found in L. borealis and N. humeralis but not in other species based on previous single-city results [4]. We predicted that both species would positively respond to affluence as measured by median household income. Furthermore, we expected no correlation between affluence and any specific type of land cover and species-specific bat-land cover associations.

Results

Between 2015 and 2018, we collected 10,899 bat passes that met identification criteria. From these bat passes we identified 805 passes of E. fuscus, 2609 passes of L. borealis, 158 passes of L. cinereus, 585 passes of L. noctivagans, 1857 passes of N. humeralis, 1016 passes of P. subflavus, and 350 passes of T. brasiliensis. Among the 43 sites, deciduous forest was the most common land cover type followed by evergreen forest and pasture/hay (Table 1). Of all 15 land cover types, 13 were present in all sites (Table 1). High intensity developed land was absent in one site and emergent herbaceous wetland absent in seven sites. The average income measured at the 43 sites was $45,882.42 (US dollar). The range for income was $32,075.15 to $65,285.81 (all US dollar, Table 1).

Table 1 Mean, standard deviation, and range for land cover (percentage in 5 km buffer) and income (US dollar) used to examine bat-land cover/bat-income relationships

For species-specific bat activity–land cover relationships, we found that E. fuscus activity was positively correlated to all four types of urban development as well as deciduous forest (Table 2, Fig. 1). Negative correlations were found between E. fuscus activity and evergreen forest, shrub, grassland/herbaceous, cultivated crops, and both types of wetlands (Table 2, Fig. 1). Lasionycteris noctivagans showed similar association patterns to E. fuscus except that L. noctivagans did not respond to grassland/herbaceous land cover. Tadarida brasiliensis had positive correlations with urban development, and the same negative bat activity–land cover relationships as L. noctivagans. Additionally, T. brasiliensis activity decreased as barren land cover increased (Table 2, Fig. 1).

Table 2 Generalized linear mixed regression model coefficient 95% confidence interval for bat-land cover, and bat-income relationships; significant relationships in italic
Fig. 1
figure 1

Generalized linear mixed model regression coefficient 95% confidence intervals (based on 1000 rounds of bootstrapping) in relation to 0 (0 as red vertical dash line; non-overlapping means significant in yellow solid line, overlapping means non-significant in blue dash line) for bat-land cover relationships; EPFU, Eptesicus fuscus; LABO, Lasiurus borealis; LACI, L. cinereus; LANO, Lasionycteris noctivagans; NYHU, Nycticeius humeralis; PESU, Perimyotis subflavus; TABR, Tadarida brasiliensis

In contrast to the three species above, L. borealis and P. subflavus, had negative relationships with all 4 types of urban development land cover (Table 2, Fig. 1). Activity of both species also decreased with increased barren land cover. As evergreen forest, mixed forest, and/or shrub land cover increased, so did L. borealis activity. However, there was a negative relationship between emergent herbaceous wetland land cover and L. borealis activity. The only land cover type that increased P. subflavus activity was woody wetland (Table 2, Fig. 1). Lasiurus cinereus and N. humeralis activity was not associated with urban development land cover. Lasiurus cinereus activity increased with deciduous forest cover and decreased with evergreen forest, shrub, grassland/herbaceous, pasture/hay, and both types of wetland cover (Table 2, Fig. 1). Nycticeius humeralis activity decreased with pasture/hay cover and increased with evergreen forest, shrub, or woody wetland cover (Table 2, Fig. 1).

As hypothesized, there was a positive relationship between bat activity and income for L. borealis and N. humeralis (Fig. 2). We also found a negative correlation between income and L. cinereus activity (Fig. 2). No bat-income relationship was found for the other species. In the post hoc modeling for L. borealis and N. humeralis, we examined interactions between income and land cover type that had a positive relationship on bat activity. We found a positive interaction term between income and evergreen forest for L. borealis activity (Table 3), suggesting the positive effect of income on L. borealis activity was stronger in areas with more evergreen forest. No interaction term was found for other land cover–income pairs for L. borealis and N. humeralis. We examined interactions between income and land cover variables that had a negative relationship with L. cinereus activity. We found a positive interaction term between income and woody wetland for L. cinereus activity (Table 3), suggesting L. cinereus activity would decrease more as income increased in areas with more woody wetlands.

Fig. 2
figure 2

Generalized linear mixed model regression coefficient 95% confidence intervals (based on 1000 rounds of bootstrapping) in relation to 0 (0 as red vertical dash line; non-overlapping means significant in yellow solid line, overlapping means non-significant in blue dash line) and scatter plot with estimated significant regression trend line/range for bat-income relationships; EPFU, Eptesicus fuscus; LABO, Lasiurus borealis; LACI, L. cinereus; LANO, Lasionycteris noctivagans; NYHU, Nycticeius humeralis; PESU, Perimyotis subflavus; TABR, Tadarida brasiliensis

Table 3 Post hoc generalized linear mixed regression model coefficient 95% confidence interval for income and land cover interaction term for bat activity (significant terms in italic); other variables in the models were reported in Additional file 2

Discussion

Our study showed that bat-land cover associations at the landscape scale were consistent with previous studies at finer scales. Both L. borealis and N. humeralis are tree roosting species [37, 38] associated with tree cover at the local scale (e.g. [4, 39, 40]). Herein, we found the same relationships between bat activity and certain types of forest cover at a much broader scale. Moreover, N. humeralis can additionally roost in buildings whereas L. borealis roosts exclusively in trees [37, 38, 41]. Thus, it was not surprising that we found L. borealis negatively associated with urban development whereas the response of N. humeralis to urban development was neutral.

Interestingly, our study also showed that the luxury effect could be found across landscapes. After accounting for spatial and temporal relationships, two previously studied species, L. borealis and N. humeralis, showed the same positive response to affluence as found in a single-city study [4]. However, certain ecological traits might not be sufficient to explain why the luxury effect could be found across different landscapes beyond urban ecosystems. This is because income was not correlated with any land cover type across our study area. Furthermore, there was a positive interaction term between income and evergreen forest cover for L. borealis, suggesting that the luxury effect was more profound in areas with more evergreen forest. This interaction, consistent with Li and Wilkins [4], suggests that income and land cover habitat might function through different mechanisms to shape the bat activity patterns.

Certain mechanisms have been suggested for the luxury effect. For sedentary organisms, such as plants, human actions determined by affluence can directly alter the physical environment, such as water availability, which may favor certain ecological traits [1, 6, 7]. Similarly, water sources could potentially affect more mobile species. Our land cover data has a 30-m resolution [42], which did not capture small water bodies such as small creeks, ponds, or swimming pools. The importance of both natural and artificial water source to bats has been established (e.g. [43,44,45,46]) and wealthier people prefer to purchase properties that include water bodies for various reasons [27, 28, 47]. Thus, we speculate that water availability, especially small water bodies might contribute to the luxury effect found at the broad scale.

In addition to water sources, food or prey sources, which may be influenced by multiple environmental variables have been shown to affect predators’ distribution [48]. Studies of highly mobile meso-carnivores demonstrated that wealthier neighborhoods had higher vegetation coverage to support more small mammals, which were the prey for urban coyotes [14, 49]. Similarly, the importance of insect diversity and abundance to bats have been well documented (e.g. [50, 51]). Even though very few studies have compared insect diversity and abundance across landscapes with varying affluence, it is evident that high affluence organic farming promotes high insect diversity and abundance [23, 24, 52, 53]. Thus, we also speculate that insect food sources might also contribute to the luxury effect. Future studies should examine the luxury effect on insects at broad spatial scales.

Both water and food sources could potentially contribute to the luxury effect we found for L. borealis and N. humeralis. However, it is puzzling that we also found L. cinereus activity negatively associated with income, independent of land cover. In Li and Wilkins [4], L. cinereus activity was only recorded on the outskirts of the city, thus was not tested for the luxury effect. Lasiurus cinereus is a solitary tree foliage roosting species [54], which is consistent with its association with deciduous forests in our study. Lasiurus cinereus also prefers open space to forage [41, 54]. However, none of these ecological traits explains the association with low income areas. One possible explanation could be that L. cinereus responded to nighttime illumination. Regional economic activities have been found to be correlated to nighttime illumination [55]. How bats responded to nighttime illumination have been found for certain species (e.g. [56, 57]), a study on nighttime light effects on L. cinereus may help to understand our result.

As Ackley et al. [2] suggested, affluence might be a holistic index of habitat structure that characterizes a combination of many different aspects of the environment. Thus, future studies should focus on understanding what exactly bats perceive as habitat. It is also important to point out that the income range in this study was not particularly large. The mobile transect survey we used in our study required a certain length of a transect and a constant sampling speed [58] and thus prevented sampling in some urban areas with the highest income. Additionally, our study region lacked extreme low-income areas. Future studies should consider other sampling methods to include high-income urban areas and expand the spatial scale to include areas with lower oncome. Lastly, based on the intermediate disturbance theory, the relationship between bat activity and income could be non-linear [59, 60], which is also a topic for future studies.

For species that did not show linear association with income, our study provides a broad scale examination of their land cover associations, which are consistent with each species ecology based on fine scale studies. For example, P. subflavus, is tree roosting species that uses wetlands and forages on aquatic insects [44, 61,62,63]. Our analysis shows the same association with woody wetland at the landscape scale. Our results also support the importance of functional guilds in land cover selection for bats [36, 64,65,66]. For example, open space foragers, such as E. fuscus or T. brasiliensis showed associations with urban cover whereas tree roosting species generally associated with forest cover. Furthermore, our results showed that no single land cover type could benefit all seven species, suggesting the importance of conserving different land cover types for the conservation of bats. The species-specific bat activity-land cover associations also call for more fine scale ecological studies to answer questions such as why L. noctivagans showed a positive association with deciduous forests but a negative one with evergreen forest. Lastly, our study analyzed the composition of land cover types and it is important to note that the spatial configuration of land cover types also affects the distribution of organisms [67]. Future studies should address how connectivity (or fragmentation) of certain types of land cover affects bat activity.

Conclusion

Our study demonstrated that the luxury effect in urban ecosystems can be found beyond cities. It is an ecological pattern that is evident across different landscapes at a broad spatial scale. The luxury effect may function through a mechanism or mechanisms different from simple bat-land cover associations. Fine scale studies are needed to better understand the specific habitat requirements of bats, especially due to the rapid urbanization and increase of anthropogenic environments. We also suggest that the luxury effect should be examined in different taxa, worldwide. The consistency of the luxury effect between spatial scales we found in the state of North Carolina, USA and in one city, Waco, Texas, USA, suggests that the luxury effect is broad and scalable. Importantly, the luxury effect has cascading impacts on both the environment and on humans [7]. There is an enormous value of bats to humans because of the ecosystem services they provide such as insect predation, seed dispersal, and pollination [68,69,70,71]. Potentially, the ecosystem services provided by bats may be unequally received by people with different levels of income, as found in urban vegetation, due to the luxury effect [3, 31, 72].

Methods

Study region and sampling site selection

We recorded bat activity across North Carolina (NC), USA (Fig. 3a). All seven target species have statewide distributions [35]. North Carolina is approximately 140,000 km2 with three geographical regions (west to east): mountains, piedmont, and coastal plain. According to US Census data in 2016, the median household income of NC is approximately $50,000 (the 38th in the US), lower than the national value (approx. $59,000). The two largest metropolitan areas (Charlotte, NC and Raleigh, NC) are in the piedmont region. Smaller urban centers (e.g. Asheville, NC and Wilmington, NC) are present in the mountains or coastal plain regions, functioning as regional hubs and tourism destinations [73].

Fig. 3
figure 3

a North American Bat Monitoring Program (NABat) grid cells selected for the mobile transect survey across North Carolina with the 2016 American Community Survey 5-year estimates of median household income shown at census block group scale (darker shade in blue indicates higher level of income, white indicates no data, income data available at https://www.census.gov/programs-surveys/acs/data/data-via-ftp.html). b A demonstration of a survey site, including one NABat grid cell, one mobile transect, one 5 km buffer, and 100 random points for extracting income data

We followed the North American Bat Monitoring Program (NABat) grid cell framework to select sites for field sampling. Detailed NABat grid cell framework and sampling site selection criteria were described in Loeb et al. [58] and Li and Kalcounis-Rueppell [35]. The NABat grid cell framework divides North America into 10 km by 10 km grid cells. All grid cells are ranked with a generalized random-tessellation stratified survey design algorithm, which allows subsampling of NABat grid cells to form a set of sites that are spatially balanced yet randomized [58]. We treated one NABat grid cell as a sampling site. In 2015, we selected 32 sites for field sampling. In each year from 2016 to 2018 we re-evaluated site availability based on logistic constraints and added sites to replace sites that were no longer available for sampling. We had 39, 33, 34 sites for 2016 to 2018 field seasons respectively. Overall, we had 43 unique sampling sites across the state of NC that represented different landscapes (Fig. 3a). Within the 43 sties we collected 276 individual samples.

Bat mobile transect survey and acoustic identification

To sample bat activity at a site, we used the NABat mobile transect survey protocol. At a sampling site, we mapped a 30–35 km non-overlapping transect before the field season. Details of transect road selection and an example of the transect spatial layout are presented in Li and Kalcounis-Rueppell [35]. Briefly, we used Anabat SD2 bat detectors and related car mount accessories (Titley Scientific, Australia) for the NABat mobile transect survey protocol to record bat acoustic activity. All sampling was in June and July of each year. We only sampled on nights with no rain or heavy fog, and with wind speeds less than 10 km/h. The transect survey started 45 min after sunset and was driven at a constant speed of 32 km/h and followed all local traffic rules. Each transect was driven twice a season within 7 nights. We also coordinated transect sampling dates at sites across years for consistency. For example, if a site was sampled early in June the first year, we endeavored to sample it early in June in subsequent years. For each transect sample, we recorded date, temperature, relative humidity, wind speed, and percent cloud cover at the beginning and the end of the survey. We used the mean of each variable collected at the beginning and the end of the survey (except for date) as survey covariates for the transect sample.

All acoustic recordings were processed with Analook (Version 4.2g, Titley Scientific, Australia) for manual species identification, as in Li and Kalcounis-Rueppell [35]. We only included recordings with at least three complete and clear bat echolocation pulses for species identification and defined each qualified recording as a bat pass. The first author compared all recordings based on bat echolocation pulse characteristics (measured by AnaLook) suggested by O’Farrell et al. [74] and Kunz and Parsons [75] with a known bat echolocation library described in Li and Kalcounis-Rueppell [35] to identify species. The species identification process was conservative. Certain species pairs, such as E. fuscus and L. noctivagans, or L. borealis and N. humeralis, can produce similar echolocation pulses. For these pairs, we did not identify a recording to species unless the unique characteristics of a species were found in multiple pulses. Further, we did not differentiate L. borealis from L. seminolus (Seminole bats) as these two species cannot be identified acoustically [76]. We only identified the seven target species. For each species, the dependent variable for statistical analyses was bat activity measured as number of bat passes per transect. We did not standardize bat activity by survey time because the amount of time to drive a transect was uniform (approximately 60 min). Further, all dependent variables and statistical analyses were species specific and we did not make any comparisons across species because there are species biases in acoustic methods of bat surveys [77].

Land cover and affluence data collection

We used land cover as an indicator of general habitat structure experienced by bats [78, 79]. Our land cover data source was the National Land Cover Database 2016 (NLCD 2016, [42]). Using ArcMap (10.4.1, ESRI, California), we generated a 5 km radius buffer around each transect. We selected 5 km as the buffer radius because it is the common nightly active range of most target bat species [30, 35, 80, 81] and we were interested in examining the luxury effect at the landscape scale, not the local scale [40, 67, 82]. We used ArcMap to extract land cover raster images from NLCD 2016 and used FRAGSTATS [83] to calculate the percentage of each land cover type within the buffer as the indicator of the amount of each land cover type available. All land cover types described in Yang et al. [42] were considered. A preliminary multi-collinearity test of all land cover types showed that all four developed land cover types (open space, low, medium, high intensities) were highly correlated (all correlation coefficients > 0.6, all variance inflation factors > 5, Additional file 1; [84]).

To avoid multi-collinearity, we modeled each land cover type separately. Further, we opted to examine all land cover types without combining highly correlated land cover types because different types of urban land cover may function as specific habitat structure for bats. For example, high intensity developed land which includes urban center office buildings might function as roosting habitat for bats [85, 86], and open space developed land which includes parks might function as foraging habitat [87, 88].

We used median household income as the indicator of affluence [1, 7]. For our median household income variable, we used the 2016 American Community Survey (ACS) 5-year estimates published by the US Census (https://www.census.gov/programs-surveys/acs/data/data-via-ftp.html), because we needed income data from all areas regardless of the community size (Fig. 3a; [89]). Median household income was available at the census block group scale from the ACS 5-year estimates. Because most census blocks did not align with the 5 km buffer and some blocks only overlay with the buffer by a small portion, we decided to use random points to calculate the weighted average median household income for each buffer. At each bat sampling site, within the 5 km radius buffer of the transect, we used ArcGIS to generate 100 spatially random points (Fig. 3b) to extract 100 median household income measurements for a site. We explicitly excluded census block groups with no data during this process. We calculated the mean of 100 measurements and used the mean as the “income” variable to represent the affluence of each bat sampling transect. Lastly, we calculated the variance inflation factor for income with land cover and found no collinearity (all correlation coefficients < 0.3, all variance inflation factors < 2, Additional file 1).

Statistical analyses

All analyses were done in R (version 3.4.2, [90]). We modeled each independent variable separately then explored possible interaction terms, post hoc, for selected land cover type and income pairs. We found that for all seven species, bat activity variance was different from the mean (all variance to mean ratios > 4). Therefore, we modeled bat activity with a negative binomial distribution [91, 92]. We also tested all possible spatial and temporal autocorrelations in bat activity in the preliminary analysis, and adjusted models accordingly [93].

Temporal autocorrelation of bat activity could occur between the two sample nights at a site within a season. Since all sites had two nights within a season, we first used Wilcoxon signed-rank test to compare all night one with all night two, independent from survey covariates. We found a difference between night one and night two (p < 0.05 for all species). Next, we used generalized linear models as the preliminary analysis to if the survey covariates would affect bat activity. We found that all species except for L. noctivagans were significantly affected by temperature (Additional file 2), consistent with previous work in the area [94]. Other survey covariates mentioned above had no effect on bat activity. Therefore, we included temperature as a covariate for all bat activity-land cover/income models except for L. noctivagans.

All target bat species may migrate regionally during the winter and cause local population size difference between summers [95]. Hence, we also conducted preliminary analyses with generalized linear models to evaluate if year had an effect on bat activity. We found that year was a significant variable for all bat species (Fig. 4, Additional file 2). Thus, we used generalized linear mixed-effects models to model bat activity-land cover/income relationship with year as a random effect using R package lme4 [84]. We further tested if the effect of temperature was nested within years by testing the interaction term between year and temperature. We found that for E. fuscus, N. humeralis, and P. subflavus, the effect of temperature varied between years (Additional file 2). Thus, for these three species, temperature was modeled as a covariate nested within the random effect of year. In this way, temperature was treated as a random slope and the relationship between bat activity and temperature was allowed to vary in slope with each year.

Fig. 4
figure 4

Box plot of species-specific bat activity (number of bat passes per transect) in relation to sampling years; all species responded to year significantly; EPFU, Eptesicus fuscus; LABO, Lasiurus borealis; LACI, L. cinereus; LANO, Lasionycteris noctivagans; NYHU, Nycticeius humeralis; PESU, Perimyotis subflavus; TABR, Tadarida brasiliensis

To test for potential spatial autocorrelation, we used the null generalized linear mixed-effects model including only temperature and year for all species but L. noctivagans, and only year for L. noctivagans, to calculate model residuals for Moran’s I test using R package ape v5.2 [96, 97]. Using p < 0.05 as the significance criterion, we found spatial autocorrelation in E. fuscus, L. borealis, L. cinereus, L. noctivagans, and N. humeralis. For these five species, we calculated a spatial autocovariate using the distance-band based neighbor scheme and neighbor weights of inverse distance with R package spdep v0.8-1 [14, 96, 98]. We included the spatial autocovariate to account for spatial autocorrelation in the models for these species. Although other modeling techniques can address spatial autocorrelation [96, 99] we used this approach to support the negative binomial distribution for our bat activity variable (see [100]).

For a more robust inference, we applied bootstrapping to the generalized linear mixed-effects models and generated confidence intervals for independent variables as suggested by Bates et al. [101]. We resampled 1000 times for each generalized linear mixed-effects model. As we modeled each independent variable separately, our models took approximately 6–7 degrees of freedom (one independent variable of interest, covariates temperature, spatial autocovariate, random effect year of 4 categories). Thus, our sample size of 276 was adequate for generalized linear mixed-effects models and bootstrapping [102]. From the bootstrapping process, we generated the 95% confidence interval for land cover and income. If the 95% confidence interval did not overlap with 0, the independent variable had a significant relationship with bat activity. For post hoc modeling, we applied the same modeling method to the interaction term between income and each significant land cover type variable for each species where income was significant. We only conducted post hoc modeling for land cover type variables that had the same direction of effect, positive or negative, as income. For example, if income was found positively correlated with a species’ activity, we conducted post hoc modeling for the interaction term between income and each land cover type that was also positively correlated to bat activity. In a post hoc model, income and a selected land cover type and their interaction term were included (Additional file 2).

Availability of data and materials

The bat datasets generated during and/or analyzed during the current study are available at the North American Bat Monitoring Program repository (https://www.nabatmonitoring.org/). Complied bat/land cover/income data are available from the corresponding author on reasonable request.

Abbreviations

ACS:

American Community Survey

EPFU:

big brown bat Eptesicus fuscus

LABO:

Eastern red bat Lasiurus borealis

LACI:

hoary bat Lasiurus cinereus

LANO:

silver haired bat Lasionycteris noctivagans

NABat:

North American Bat Monitoring Program

NC:

North Carolina

NLCD:

National Land Cover Database

NYHU:

evening bat Nycticeius humeralis

PESU:

tri-colored bat Perimyotis subflavus

TABR:

Mexican free-tailed bat Tadarida brasiliensis

References

  1. Hope D, Gries C, Zhu W, et al. Socioeconomics drive urban plant diversity. Proc Natl Acad Sci. 2003;100:8788–92. https://doi.org/10.1073/pnas.1537557100.

    Article  CAS  PubMed  Google Scholar 

  2. Ackley JW, Wu J, Angilletta MJ, et al. Rich lizards: how affluence and land cover influence the diversity and abundance of desert reptiles persisting in an urban landscape. Biol Conserv. 2015;182:87–92. https://doi.org/10.1016/j.biocon.2014.11.009.

    Article  Google Scholar 

  3. Clarke LW, Jenerette GD, Davila A. The luxury of vegetation and the legacy of tree biodiversity in Los Angeles, CA. Landsc Urban Plan. 2013;116:48–59. https://doi.org/10.1016/j.landurbplan.2013.04.006.

    Article  Google Scholar 

  4. Li H, Wilkins KT. Patch or mosaic: bat activity responds to fine-scale urban heterogeneity in a medium-sized city in the United States. Urban Ecosyst. 2014. https://doi.org/10.1007/s11252-014-0369-9.

    Article  Google Scholar 

  5. van Heezik Y, Freeman C, Porter S, Dickinson KJM. Garden size, householder knowledge, and socio-economic status influence plant and bird diversity at the scale of individual gardens. Ecosystems. 2013;16:1442–54. https://doi.org/10.1007/s10021-013-9694-8.

    Article  Google Scholar 

  6. Wang H-F, Qureshi S, Knapp S, et al. A basic assessment of residential plant diversity and its ecosystem services and disservices in Beijing, China. Appl Geogr. 2015;64:121–31. https://doi.org/10.1016/j.apgeog.2015.08.006.

    Article  Google Scholar 

  7. Leong M, Dunn RR, Trautwein MD. Biodiversity and socioeconomics in the city: a review of the luxury effect. Biol Lett. 2018;14:20180082. https://doi.org/10.1098/rsbl.2018.0082.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Grove JM, Locke DH, O’Neil-Dunne JPM. An ecology of prestige in New York City: examining the relationships among population density, socio-economic status, group identity, and residential canopy cover. Environ Manag. 2014;54:402–19. https://doi.org/10.1007/s00267-014-0310-2.

    Article  Google Scholar 

  9. Larsen L, Harlan SL. Desert dreamscapes: residential landscape preference and behavior. Landsc Urban Plan. 2006;78:85–100. https://doi.org/10.1016/j.landurbplan.2005.06.002.

    Article  Google Scholar 

  10. Meyer WB. Bringing hypsography back in: altitude and residence in American cities. Urban Geogr. 1994;15:505–13. https://doi.org/10.2747/0272-3638.15.6.505.

    Article  Google Scholar 

  11. Niemelä J, Breuste JH, Guntenspergen G, et al. Urban ecology: patterns, processes, and applications. Oxford: Oxford University Press; 2011.

    Book  Google Scholar 

  12. Belaire JA, Westphal LM, Minor ES. Different social drivers, including perceptions of urban wildlife, explain the ecological resources in residential landscapes. Landsc Ecol. 2016;31:401–13. https://doi.org/10.1007/s10980-015-0256-7.

    Article  Google Scholar 

  13. Fuller RA, Warren PH, Armsworth PR, et al. Garden bird feeding predicts the structure of urban avian assemblages. Divers Distrib. 2008;14:131–7. https://doi.org/10.1111/j.1472-4642.2007.00439.x.

    Article  Google Scholar 

  14. Wine S, Gagné SA, Meentemeyer RK. Understanding human–coyote encounters in urban ecosystems using citizen science data: what do socioeconomics tell us? Environ Manag. 2015;55:159–70. https://doi.org/10.1007/s00267-014-0373-0.

    Article  Google Scholar 

  15. Boone CG. Environmental justice as process and new avenues for research. Environ Justice. 2008;1:149–54. https://doi.org/10.1089/env.2008.0530.

    Article  Google Scholar 

  16. Cohen M, Baudoin R, Palibrk M, et al. Urban biodiversity and social inequalities in built-up cities: new evidences, next questions. The example of Paris, France. Landsc Urban Plan. 2012;106:277–87. https://doi.org/10.1016/j.landurbplan.2012.03.007.

    Article  Google Scholar 

  17. O’Neill MS, Jerrett M, Kawachi I, et al. Health, wealth, and air pollution: advancing theory and methods. Environ Health Perspect. 2003;111:1861–70.

    Article  Google Scholar 

  18. Ripplinger J, Franklin J, Collins SL. When the economic engine stalls—a multi-scale comparison of vegetation dynamics in pre- and post-recession Phoenix, Arizona, USA. Landsc Urban Plan. 2016;153:140–8. https://doi.org/10.1016/j.landurbplan.2016.05.009.

    Article  Google Scholar 

  19. Donihue CM, Lambert MR. Adaptive evolution in urban ecosystems. Ambio. 2015;44:194–203. https://doi.org/10.1007/s13280-014-0547-2.

    Article  PubMed  Google Scholar 

  20. Myers N, Mittermeier RA, Mittermeier CG, et al. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853.

    Article  CAS  Google Scholar 

  21. Newbold T, Hudson LN, Hill SLL, et al. Global effects of land use on local terrestrial biodiversity. Nature. 2015;520:45–50. https://doi.org/10.1038/nature14324.

    Article  CAS  Google Scholar 

  22. Turner WR, Brandon K, Brooks TM, et al. Global biodiversity conservation and the alleviation of poverty. Bioscience. 2012;62:85–92. https://doi.org/10.1525/bio.2012.62.1.13.

    Article  Google Scholar 

  23. Lambert DM, Sullivan P, Claassen R, Foreman LF. Conservation-compatible practices and programs: who participates?. Washington, DC: United States Department of Agriculture, Economic Research Service; 2006.

    Google Scholar 

  24. Page G, Bellotti B. Farmers value on-farm ecosystem services as important, but what are the impediments to participation in PES schemes? Sci Total Environ. 2015;515–516:12–9. https://doi.org/10.1016/j.scitotenv.2015.02.029.

    Article  CAS  PubMed  Google Scholar 

  25. Piemonti AD, Babbar-Sebens M, Luzar EJ. Optimizing conservation practices in watersheds: do community preferences matter? Water Resour Res. 2013;49:6425–49. https://doi.org/10.1002/wrcr.20491.

    Article  Google Scholar 

  26. Martinuzzi S, Radeloff VC, Joppa LN, et al. Scenarios of future land use change around United States’ protected areas. Biol Conserv. 2015;184:446–55. https://doi.org/10.1016/j.biocon.2015.02.015.

    Article  Google Scholar 

  27. Kirk RW, Bolstad PV, Manson SM. Spatio-temporal trend analysis of long-term development patterns (1900–2030) in a Southern Appalachian county. Landsc Urban Plan. 2012;104:47–58. https://doi.org/10.1016/j.landurbplan.2011.09.008.

    Article  Google Scholar 

  28. Esparza A. The exurbanization process and rural housing markets. In: Lapping M, Marcouiller D, Furuseth O, editors. Rural housing, exurbanization, and amenity-driven development. London: Routledge; 2011.

    Google Scholar 

  29. Runge CA, Plantinga AJ, Larsen AE, et al. Unintended habitat loss on private land from grazing restrictions on public rangelands. J Appl Ecol. 2018. https://doi.org/10.1111/1365-2664.13271.

    Article  Google Scholar 

  30. Kunz TH, Fenton MB. Bat ecology. Chicago: University of Chicago Press; 2006.

    Google Scholar 

  31. Schwarz K, Fragkias M, Boone CG, et al. Trees grow on money: urban tree canopy cover and environmental justice. PLoS ONE. 2015;10:e0122051. https://doi.org/10.1371/journal.pone.0122051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Fenton MB. Science and the conservation of bats. J Mammal. 1997;78:1–14.

    Article  Google Scholar 

  33. Voigt CC, Kingston T. Bats in the anthropocene: conservation of bats in a changing world. New York: Springer International Publishing; 2016.

    Book  Google Scholar 

  34. Azam C, Viol IL, Julien J-F, et al. Disentangling the relative effect of light pollution, impervious surfaces and intensive agriculture on bat activity with a national-scale monitoring program. Landsc Ecol. 2016;31:2471–83. https://doi.org/10.1007/s10980-016-0417-3.

    Article  Google Scholar 

  35. Li H, Kalcounis-Rueppell M. Separating the effects of water quality and urbanization on temperate insectivorous bats at the landscape scale. Ecol Evol. 2018;8:667–78. https://doi.org/10.1002/ece3.3693.

    Article  PubMed  Google Scholar 

  36. Ducci L, Agnelli P, Febbraro MD, et al. Different bat guilds perceive their habitat in different ways: a multiscale landscape approach for variable selection in species distribution modelling. Landsc Ecol. 2015;30:2147–59. https://doi.org/10.1007/s10980-015-0237-x.

    Article  Google Scholar 

  37. Shump KA, Shump AU. Lasiurus borealis. Mamm Species Arch. 1982;183:1–6. https://doi.org/10.2307/3503843.

    Article  Google Scholar 

  38. Watkins LC. Nycticeius humeralis. Mamm Species. 1972. https://doi.org/10.2307/3503945.

    Article  Google Scholar 

  39. Dixon MD. Relationship between land cover and insectivorous bat activity in an urban landscape. Urban Ecosyst. 2011;15:683–95. https://doi.org/10.1007/s11252-011-0219-y.

    Article  Google Scholar 

  40. Gallo T, Lehrer EW, Fidino M, et al. Need for multiscale planning for conservation of urban bats. Conserv Biol. 2018;32:638–47. https://doi.org/10.1111/cobi.13047.

    Article  PubMed  Google Scholar 

  41. Ammerman LK, Hice CL, Schmidly DJ. Bats of Texas. College Station: Texas A&M University Press; 2012.

    Google Scholar 

  42. Yang L, Jin S, Danielson P, Homer C, Gass L, Bender SM, et al. A new generation of the United States National Land Cover Database: Requirements, research priorities, design, and implementation strategies. ISPRS J Photogramm Remote Sens. 2018;146:108–23.

    Article  Google Scholar 

  43. Korine C, Adams R, Russo D, et al. Bats and water: anthropogenic alterations threaten global bat populations. In: Voigt CC, Kingston T, editors. Bats in the anthropocene: conservation of bats in a changing world. Cham: Springer International Publishing; 2016. p. 215–24.

    Chapter  Google Scholar 

  44. Parker KA, Springall BT, Garshong RA, et al. Rapid increases in bat activity and diversity after wetland construction in an urban ecosystem. Wetlands. 2018. https://doi.org/10.1007/s13157-018-1115-5.

    Article  PubMed  Google Scholar 

  45. Salvarina I. Bats and aquatic habitats: a review of habitat use and anthropogenic impacts. Mammal Rev. 2016;46:131–43. https://doi.org/10.1111/mam.12059.

    Article  Google Scholar 

  46. Stahlschmidt P, Pätzold A, Ressl L, et al. Constructed wetlands support bats in agricultural landscapes. Basic Appl Ecol. 2012;13:196–203. https://doi.org/10.1016/j.baae.2012.02.001.

    Article  Google Scholar 

  47. Sirmans GS, MacDonald L, Macpherson DA, Zietz EN. The value of housing characteristics: a meta analysis. J Real Estate Finance Econ. 2006;33:215–40. https://doi.org/10.1007/s11146-006-9983-5.

    Article  Google Scholar 

  48. Tsuji Y, Ito TY, Kaneko Y. Variation in the diets of Japanese martens Martes melampus. Mammal Rev. 2019;49:121–8. https://doi.org/10.1111/mam.12147.

    Article  Google Scholar 

  49. Magle SB, Lehrer EW, Fidino M. Urban mesopredator distribution: examining the relative effects of landscape and socioeconomic factors. Anim Conserv. 2016;19:163–75. https://doi.org/10.1111/acv.12231.

    Article  Google Scholar 

  50. Blakey RV, Kingsford RT, Law BS, Stoklosa J. Floodplain habitat is disproportionately important for bats in a large river basin. Biol Conserv. 2017;215:1–10. https://doi.org/10.1016/j.biocon.2017.08.030.

    Article  Google Scholar 

  51. Threlfall CG, Law B, Banks PB. Sensitivity of insectivorous bats to urbanization: implications for suburban conservation planning. Biol Conserv. 2012;146:41–52. https://doi.org/10.1016/j.biocon.2011.11.026.

    Article  Google Scholar 

  52. Bengtsson J, Ahnström J, Weibull A-C. The effects of organic agriculture on biodiversity and abundance: a meta-analysis. J Appl Ecol. 2005;42:261–9. https://doi.org/10.1111/j.1365-2664.2005.01005.x.

    Article  Google Scholar 

  53. Wickramasinghe LP, Harris S, Jones G, Vaughan Jennings N. Abundance and species richness of nocturnal insects on organic and conventional farms: effects of agricultural intensification on bat foraging. Conserv Biol. 2004;18:1283–92.

    Article  Google Scholar 

  54. Shump KA, Shump AU. Lasiurus cinereus. Mamm Species. 1982. https://doi.org/10.2307/3503878.

    Article  Google Scholar 

  55. Doll CNH, Muller J-P, Morley JG. Mapping regional economic activity from night-time light satellite imagery. Ecol Econ. 2006;57:75–92. https://doi.org/10.1016/j.ecolecon.2005.03.00.

    Article  Google Scholar 

  56. Russo D, Cistrone L, Libralato N, et al. Adverse effects of artificial illumination on bat drinking activity. Anim Conserv. 2017. https://doi.org/10.1111/acv.12340.

    Article  Google Scholar 

  57. Schoeman MC. Light pollution at stadiums favors urban exploiter bats. Anim Conserv. 2015. https://doi.org/10.1111/acv.1222o.

    Article  Google Scholar 

  58. Loeb SC, Rodhouse TJ, et al. A plan for the North American Bat Monitoring Program (NABat). 2015.

  59. McKinney ML. Effects of urbanization on species richness: a review of plants and animals. Urban Ecosyst. 2008;11:161–76. https://doi.org/10.1007/s11252-007-0045-4.

    Article  Google Scholar 

  60. Yuan ZY, Jiao F, Li YH, Kallenbach RL. Anthropogenic disturbances are key to maintaining the biodiversity of grasslands. Sci Rep. 2016;6:22132. https://doi.org/10.1038/srep22132.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Fujita MS, Kunz TH. Pipistrellus subflavus. Mamm Species. 1984;228:1–6. https://doi.org/10.2307/3504021.

    Article  Google Scholar 

  62. Kalcounis-Rueppell MC, Payne VH, Huff SR, Boyko AL. Effects of wastewater treatment plant effluent on bat foraging ecology in an urban stream system. Biol Conserv. 2007;138:120–30. https://doi.org/10.1016/j.biocon.2007.04.009.

    Article  Google Scholar 

  63. Neece BD, Loeb SC, Jachowski DS. Variation in regional and landscape effects on occupancy of temperate bats in the southeastern US. PLOS ONE. 2018;13:e0206857. https://doi.org/10.1371/journal.pone.0206857.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Jung K, Threlfall CG. Trait-dependent tolerance of bats to urbanization: a global meta-analysis. Proc R Soc B Biol Sci. 2018. https://doi.org/10.1098/rspb.2018.1222.

    Article  Google Scholar 

  65. Luck GW, Smallbone L, Threlfall C, Law B. Patterns in bat functional guilds across multiple urban centres in south-eastern Australia. Landsc Ecol. 2013;28:455–69. https://doi.org/10.1007/s10980-012-9842-0.

    Article  Google Scholar 

  66. Threlfall CG, Law B, Penman T, Banks PB. Ecological processes in urban landscapes: mechanisms influencing the distribution and activity of insectivorous bats. Ecography. 2011;34:814–26. https://doi.org/10.1111/j.1600-0587.2010.06939.x.

    Article  Google Scholar 

  67. Turner M, Gardner R, O’Neill R. Landscape ecology in theory and practice: pattern and process. New York: Springer; 2001.

    Google Scholar 

  68. Boyles JG, Cryan PM, McCracken GF, Kunz TH. Economic importance of bats in agriculture. Science. 2011;332:41–2.

    Article  Google Scholar 

  69. Kunz TH, de Torrez EB, Bauer D, et al. Ecosystem services provided by bats. Ann N Y Acad Sci. 2011;1223:1–38. https://doi.org/10.1111/j.1749-6632.2011.06004.x.

    Article  PubMed  Google Scholar 

  70. Russo D, Bosso L, Ancillotto L. Novel perspectives on bat insectivory highlight the value of this ecosystem service in farmland: research frontiers and management implications. Agric Ecosyst Environ. 2018;266:31–8. https://doi.org/10.1016/j.agee.2018.07.024.

    Article  Google Scholar 

  71. Saldaña-Vázquez RA, Castaño JH, Baldwin J, Pérez-Torres J. Does seed ingestion by bats enhance germination? A new meta-analysis 15 years later. Mammal Rev. 2019;49:201–9. https://doi.org/10.1111/mam.12153.

    Article  Google Scholar 

  72. Brooks KR, Kelley W, Amiri S. Social equity of street trees in the pedestrian realm. Pap Appl Geogr. 2016;2:216–35. https://doi.org/10.1080/23754931.2015.1121163.

    Article  Google Scholar 

  73. North Carolina Department of Commerce. Travel economic impact model. 2013. https://www.nccommerce.com/tourism/research/economic-impact/teim. Accessed 19 Oct 2018.

  74. O’Farrell MJ, Miller BW, Gannon WL. Qualitative identification of free-flying bats using the Anabat detector. J Mammal. 1999;80:11–23.

    Article  Google Scholar 

  75. Kunz TH, Parsons S. Ecological and behavioral methods for the study of bats. Baltimore: Johns Hopkins University Press; 2009.

    Google Scholar 

  76. Loeb SC, Post CJ, Hall ST. Relationship between urbanization and bat community structure in national parks of the southeastern US. Urban Ecosyst. 2009;12:197–214. https://doi.org/10.1007/s11252-008-0075-6.

    Article  Google Scholar 

  77. Russo D, Voigt CC. The use of automated identification of bat echolocation calls in acoustic monitoring: a cautionary note for a sound analysis. Ecol Indic. 2016;66:598–602.

    Article  Google Scholar 

  78. Adamo M, Tarantino C, Tomaselli V, et al. Expert knowledge for translating land cover/use maps to general habitat categories (GHC). Landsc Ecol. 2014;29:1045–67. https://doi.org/10.1007/s10980-014-0028-9.

    Article  Google Scholar 

  79. Tomaselli V, Dimopoulos P, Marangi C, et al. Translating land cover/land use classifications to habitat taxonomies for landscape monitoring: a Mediterranean assessment. Landsc Ecol. 2013;28:905–30. https://doi.org/10.1007/s10980-013-9863-3.

    Article  Google Scholar 

  80. Barclay RMR. Long- versus short-range foraging strategies of hoary (Lasiurus cinereus) and silver-haired (Lasionycteris noctivagans) bats and the consequences for prey selection. Can J Zool. 1985;63:2507–15. https://doi.org/10.1139/z85-371.

    Article  Google Scholar 

  81. Norberg UM. Vertebrate flight: mechanics, physiology, morphology, ecology and evolution. New York: Springer; 1990.

    Book  Google Scholar 

  82. Miguet P, Jackson HB, Jackson ND, et al. What determines the spatial extent of landscape effects on species? Landsc Ecol. 2015;31:1177–94. https://doi.org/10.1007/s10980-015-0314-1.

    Article  Google Scholar 

  83. McGarigal K, Cushman S, Ene E. FRAGSTATS v4: spatial pattern analysis program for categorical and continuous maps. Computer software program produced by the authors at the University of Massachusetts, Amherst. 2012.

  84. Zuur A, Ieno EN, Walker N, et al. Mixed effects models and extensions in ecology with R. New York: Springer; 2009.

    Book  Google Scholar 

  85. Li H, Wilkins KT. Selection of building roosts by Mexican free-tailed bats (Tadarida brasiliensis) in an urban area. Acta Chiropterologica. 2015;17:321–30. https://doi.org/10.3161/15081109ACC2015.17.2.007.

    Article  Google Scholar 

  86. Neubaum DJ, Wilson KR, O’shea TJ. Urban maternity-roost selection by big brown bats in Colorado. J Wildl Manag. 2007;71:728–36. https://doi.org/10.2193/2005-684.

    Article  Google Scholar 

  87. Avila-Flores R, Fenton MB. Use of spatial features by foraging insectivorous bats in a large urban landscape. J Mammal. 2005;86:1193–204. https://doi.org/10.1644/04-MAMM-A-085R1.1.

    Article  Google Scholar 

  88. Silva de Araújo MLV, Bernard E. Green remnants are hotspots for bat activity in a large Brazilian urban area. Urban Ecosyst. 2015. https://doi.org/10.1007/s11252-015-0487-z.

    Article  Google Scholar 

  89. U.S. Census Bureau. A compass for understanding and using American Community Survey data: what researchers need to know. 2009.

  90. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2017.

    Google Scholar 

  91. Burnham KP, Anderson DR. Model selection and multi-model inference: a practical information-theoretic approach. New York: Springer; 2002.

    Google Scholar 

  92. Frühwirth-Schnatter S, Frühwirth R, Held L, Rue H. Improved auxiliary mixture sampling for hierarchical models of non-Gaussian data. Stat Comput. 2009;19:479. https://doi.org/10.1007/s11222-008-9109-4.

    Article  Google Scholar 

  93. Gutzwiller KJ, Riffell SK. Using statistical models to study temporal dynamics of animal—landscape relations. In: Bissonette JA, Storch I, editors. Temporal dimensions of landscape ecology. New York: Springer; 2007. p. 93–118.

    Chapter  Google Scholar 

  94. Schimpp SA, Li H, Kalcounis-Rueppell MC. Determining species specific nightly bat activity in sites with varying urban intensity. Urban Ecosyst. 2018;21:541–50. https://doi.org/10.1007/s11252-018-0737-y.

    Article  Google Scholar 

  95. Krauel JJ, McCracken GF. Recent advances in bat migration research. In: Adams RA, Pedersen SC, editors. Bat evolution, ecology, and conservation. New York: Springer; 2013. p. 293–313.

    Chapter  Google Scholar 

  96. Bivand RS, Pebesma EJ, Rubio VG. Applied spatial data: analysis with R. New York: Springer; 2008.

    Google Scholar 

  97. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90. https://doi.org/10.1093/bioinformatics/btg412.

    Article  CAS  PubMed  Google Scholar 

  98. Augustin NH, Mugglestone MA, Buckland ST. An autologistic model for the spatial distribution of wildlife. J Appl Ecol. 1996;33:339–47. https://doi.org/10.2307/2404755.

    Article  Google Scholar 

  99. Beale CM, Lennon JJ, Yearsley JM, et al. Regression analysis of spatial data. Ecol Lett. 2010;13:246–64. https://doi.org/10.1111/j.1461-0248.2009.01422.x.

    Article  PubMed  Google Scholar 

  100. Lee D, Rushworth A, Napier G. Spatio-temporal areal unit modeling in R with conditional autoregressive priors using the CARBayesST package. J Stat Softw. 2018. https://doi.org/10.18637/jss.v084.i09.

    Article  Google Scholar 

  101. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. ArXiv14065823 Stat. 2014.

  102. Quinn GP, Keough MJ. Experimental design and data analysis for biologists. New York: Cambridge University Press; 2002.

    Book  Google Scholar 

  103. Sikes RS (2016) 2016 Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. J Mammal 97:663–688. https://doi.org/10.1093/jmammal/gyw078.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We sincerely thank B. Cherry, R. Purkey, E. Britzke, J. Patten, M. Whitby, J. Nawrocki, M. Frazer, B. Saul, W. Ricks, L. Raleigh, S. Bryan, H. Cobos, S. Cameron, M. Morse, D. Bradley, C. Hickman, M. LaVoie, J. Paris, D. Blatny, E. Corey, D. Allen, A. Boynton, K. Caldwell, J. Carpenter, A. Melford, B. Skiba, K. Weeks, A. Bledsoe, K. Confortin, T. Porter, A. Matteson, S. Scott, A. Wagner, L. Zarecky, N. Allen, C. Kendall, M. Vindigni, K. Livengood, C. Corben, B. Morgan, E. Travis, and B. Sherrill for assistance in collecting bat recordings. We deeply appreciate the lodging accommodations provided by B. Redd, C. Redd, M. Vann, T. Vann, G. Parker, and K. Parker Sr. during the field seasons. We thank J. Nelson for his support with GIS software and B. Springall, R. Petric, and B. McLean for input on the manuscript. We appreciate the thoughtful comments of two anonymous reviewers that improved the presentation of this work.

Funding

This work was supported by the U.S. Fish and Wildlife Service Competitive State Wildlife Grant Program [Grant Number SC-U2-F14AP00958] and North Carolina Wildlife Resources Commission [Grant Number WM-0321]. Both grants specified that the North American Bat Monitoring Program (NABat, [58]) should be used as the overarching framework for field site selection and bat acoustic monitoring. The funding bodies were not involved in the design of the research question, field data collection, analysis and interpretation of data, or writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

HL and MKR were responsible for funding acquisition and conceiving the project. HL and KP conducted field data collection, data analysis, and visual presentation. HL wrote the first draft of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Han Li.

Ethics declarations

Ethics approval and consent to participate

All fieldwork followed the American Society of Mammalogist’s guidelines for use of wild animals in research [103] and was approved by the North Carolina Wildlife Resources Commission (permit number 15-SC00162, 16-SC00162, 17-SC00162, 18-SC00162).

Consent to publish

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Independent variable (land cover and income) multi-collinearity analysis result. Correlation scatter plots on the lower left of the graph and correlation coefficient between paired independent variables on the upper right of the graph; land cover codes are—cls_21_5km: Developed open space; cls_22_5km: Developed low intensity; cls_23_5km: Developed medium intensity; cls_24_5km: Developed high intensity; cls_31_5km: Cultivated crops; cls_41_5km: Pasture/hay; cls_42_5km: Barren land; cls_43_5km: Grassland/herbaceous; cls_52_5km: Shrub; cls_71_5km: Deciduous forest; cls_81_5km: Evergreen forest; cls_82_5km: Mixed forest; cls_90_5km: Emergent herbaceous wetlands; cls_95_5km: Woody wetlands.

Additional file 2.

Regression model results. Preliminary regression model results on the effects of temperature and year on worksheet 1; interaction model results on income and selected land cover on worksheet 2.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, H., Parker, K.A. & Kalcounis-Rueppell, M.C. The luxury effect beyond cities: bats respond to socioeconomic variation across landscapes. BMC Ecol 19, 46 (2019). https://doi.org/10.1186/s12898-019-0262-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12898-019-0262-8

Keywords