Source: UNIV OF WISCONSIN submitted to NRP
LANDSCAPE AND LOCAL DRIVERS OF ROOSTING HABITAT FOR ENDANGERED FOREST-DWELLING BATS USING REMOTE SENSING
Sponsoring Institution
National Institute of Food and Agriculture
Project Status
COMPLETE
Funding Source
Reporting Frequency
Annual
Accession No.
1014010
Grant No.
(N/A)
Cumulative Award Amt.
(N/A)
Proposal No.
(N/A)
Multistate No.
(N/A)
Project Start Date
Oct 4, 2017
Project End Date
Sep 30, 2020
Grant Year
(N/A)
Program Code
[(N/A)]- (N/A)
Recipient Organization
UNIV OF WISCONSIN
21 N PARK ST STE 6401
MADISON,WI 53715-1218
Performing Department
Forest and Wildlife Ecology
Non Technical Summary
This project is designed to provide land managers with guidance for optimizing forest management to recover forest-roosting bats following the white-nosed syndrome epizootic. White-nose syndrome is a fungal disease that causes mass mortality in bats that hibernate in caves or mines. In Wisconsin, two species of cave bats roost in buildings and bat houses during summer, while another two species primarily use forests for day roosts and foraging habitat. These forest-dwelling species, the northern long-eared bat (Myotis septentrionalis; NLEB) and the eastern pipistrelle (Perimyotis subflavus; EAPI), depend on forest resources for critical life-history processes that drive population growth. Conservationists expect a wave of local extirpations among cave-hibernating bats because of the disease, and consequently most hibernating bats are on State or Federal watch lists with varying levels of regulatory classification (e.g. species of concern, threatened, endangered). Recovery of these populations and restoration of the ecological benefits provided by bats will require land managers and conservationists to manage forests for range of factors related to habitat, addressing both mortality sources and optimizing recruitment to achieve positive population growth following the inevitable die-off.Unfortunately, basic life-history information on Wisconsin's threatened bats is relatively unknown. In the context of recovery, managers have little information on the spatial relationships between large essentially permanent winter hibernacula (caves, mines) and dispersed summer roosting sites that are temporally ephemeral (dead trees in an appropriate state of decay). Bats give birth and provision their young while using summer roosting habitat hence quality summer habitat is critically important for recruitment and population growth. Forest management to help recover post-epizootic bat populations could include the retention or creation of dead trees, but how close to hibernacula do they need to be? How soon after tree death is a dead tree useful for roosting and how long does a dead tree provide a suitable roosting resource? How does ecological context affect the potential use of a roost tree?This project takes an interdisciplinary approach to answering these and related questions by using remote sensing of forests (first discipline) and wildlife ecology (second discipline) to mount a rigorous multiscale analysis of the location of roosting habitat relative to landscape position of known bat hibernacula and suitable forest conditions (likely a juxtaposition of bat foraging habitat and presence of suitable dead trees in Wisconsin).
Animal Health Component
100%
Research Effort Categories
Basic
(N/A)
Applied
100%
Developmental
(N/A)
Classification

Knowledge Area (KA)Subject of Investigation (SOI)Field of Science (FOS)Percent
2160860107050%
1230860107050%
Goals / Objectives
The goal of this research is a validated predictive model of use of forests by bats in Wisconsin to help aid management for the recovery of bat populations following impacts White-Nose Syndrome disease. We hypothesize that: 1) large-scale selection of forest habitat will be driven by landscape context that reflects distances from established permanent hibernacula, forest cover types, and forest fragmentation; 2) fine-scale use of forest habitat (selection of roost trees) will be dependent on probability of large-scale selection and fine-scale measures of density of dead trees, as well local features such as soil type, and proximity to wetlands (because of their influence on arthropod prey).
Project Methods
We propose three discrete steps towards deriving a rigorous predictive model of summer roosting habitat for NLEBs and EAPIs: 1) Estimation of a landscape model of summer habitat selection in Wisconsin forests 2) Estimation of a roost-site selection model conditioned on landscape level selection, and 3) Evaluation and validation of landscape roost site selection models using independently generated data.Step one will be an evaluation of geo-referenced bat detections for NLEBs and EAPIs generated by the WI DNR's bureau of Natural Heritage (please see attached letter of support). The WI DNR database contains roughly 650 and 740 summer detections of NLEB and EAPIs (identified roost trees, capture records, acoustic detections) from throughout the state in addition to geo-referenced records of visits to forested sites where these bats were not found. These detections/non-detections will be treated as experimental observations and reflect mostly a pre-WNS condition. We will use Generalized Linear Models (GLMs) or occupancy modeling (MacKenzie et al. 2005) to model detection/non-detection probabilities as functions of landscape context and remotely sensed habitat features that vary at the scale of discrete forest stands (e.g. species composition, age, stocking density). Occupancy modeling has the advantage of enabling the analyst to estimate separately site-level probabilities of occupancy and detection but it requires repeated observations. Our decision to use occupancy modeling will depend on whether the database contains natural aggregations of observations (detections/non-detection) in space or time that can be treated as repeated observations and whether aggregated data produces a reasonable sample size and gives adequate coverage of Wisconsin's forests.Under a GLM analysis, observations that are within a 2 foraging radii for bats will be considered spatially dependent and we will use random effects to account for unmeasured spatial dependence. We will estimate 2 models, one for each species (NLEB and EAPI). Maps of habitat type, land cover, and fragmentation will be derived from the forthcoming 30-m Wiscland 2 product (dnr.wi.gov/maps/gis/datalandcover.html). Other data sources, such as the annual USDA National Agricultural Statistics Service (NASS) 30-m Cropland Data Layer (CDL) will also be used for cover type mapping (www.nass.usda.gov/Research_and_Science/Cropland/SARS1a.php).Step two will involve subsampling of archived geo-referenced bat observations used in step one to reflect a representative range of predicted landscape-level occupancy or occurence probabilities. The subsampled points will then be used to study the fine-scale influence of dead trees of various ages. We will use high resolution remotely sensed imagery from NAIP (1m, available statewide 2005, 2006 [2m], 2008, 2010, 2013, 2015, imagery before 2010 is true color) to measure the density (and potentially age) of dead trees within 2 foraging radii of each subsample point. We will segment imagery to identify likely trees (Fig 2a) and use machine learning algorithms such as support vector machines for classification. (Cortes and Vapnik 1995). Training data for live vs. dead segments will be developed from manual interpretation of NAIP images and validation will be conducted using field surveys (see below) and independent sources such as Google Earth (Achanta et al. 2012). Dead trees will be categorized into age classes on the basis of when they first appear in the imagery as dead trees. Again, we will use GLMs to examine the influence of dead trees on the probability of bat detection by species. We will accomplish this by modeling the influence of dead tree measurements on residuals for step one models (above) fit to the subsample data and by modeling probability of occupancy in the subsample data with predicted landscape detection rates as a covariate.Step three will involve using the high resolution remotely sensed imagery and the step two models to select an independent set of data points for model verification. We will then visit these sites and survey for bat presence with acoustic detectors using the DNR field protocols used to generate the modeling dataset (described below). DNR monitoring of WNS in known hibernacula and step 1 modeling will be used to guide study design so that field sampling will avoid areas where absence of bats is likely to be a result of disease. We will then fit derived (steps 1 and 2) and null models to the validation dataset and evaluate the validity of our derived models on the basis of comparisons to null models using Akaike's Information Criterion (Burnham and Anderson 2003). We will then refine models as needed and produce a "heat map" of Wisconsin forests showing how relative potential as summer roosting habitat varies with orientation to hibernacula and large- and small-scale ecological drivers.Data generated for the validation data set will incorporate information on the influence of WNS-affected hibernacula. We view this as a natural experiment for validating landscape level effects and we will evaluate models using additional covariates to account for year-specific detections of WNS in discrete hibernacula.Field sampling for bat occupancy and model validation: Once a sample site is identified, we will follow USFWS's range-wide summer survey guidelines for Indiana bat and NLEB (USFWS 2016). We will start occupancy assessments by using ultrasound detectors designed to record bat echolocation calls. Recorded calls of high enough quality can be identified to species. For non-linear projects, USFWS requires four detector nights per 123 acres. If target species NLEB or EAPI are recorded, we will continue assessments by netting for 42 net nights. If target species are captured, Wisconsin Department of Natural Resources will attach radio transmitters and track bats for the duration of the life of the transmitter. We will collect habitat characteristics from roost trees to which we track bats, further adding to the model.

Progress 10/04/17 to 09/30/20

Outputs
Target Audience:The target audience include scientists who study the life history of bats and the effects of white-nosed syndrome. Other audiences include academic and agency scientists and managers interested in the conservation of North American wildlife, particularly bat. We also target lay audiences who are interested broadly in the conservation of bats. Changes/Problems: Nothing Reported What opportunities for training and professional development has the project provided?This project supported the training of a graduate student who earned a Master's of Science degree in Wildlife Ecology at the University of Wisconsin - Madison. How have the results been disseminated to communities of interest?The graduate student works for the Natural Heritage division of the Wisconsin Department of Natural Resources and has communicated her findings to managers directly. She is also an active member of The Wildlife Society and of a professional group focused on bat conservation in North America and has given several presentations at scientific meetings. Her thesis has been deposited in the Library at the University of Wisconsin -Madison and is available using library indexing services. The graduate student has two manuscripts drafted for submission to peer reviewed science journals, functionally the two chapters of the thesis. Both manuscripts are under review with co-authors and collaborators and will be submitted during 2021. What do you plan to do during the next reporting period to accomplish the goals? Nothing Reported

Impacts
What was accomplished under these goals? This analysis was based on a database of bat detections maintained by the Wisconsin Department of Natural Resources and remotely sensed imagery obtained from the P. Townsend lab (University of Wisconsin - Madison. Weused spatially thinned occurrence data from June and July 2010 to 2017 to build species distribution models. Using acoustic and capture records for tricolored bats and northern long-eared bats we developed distribution models for these species in Wisconsin and identified land cover and climatic variables contributing to habitat suitability. We then tested the performance of the habitat suitability models using an independent dataset, however due to low occurrence of tricolored bat and northern long-eared bat during post-White-Nose Syndrome (WNS) acoustic sampling (occurrences at only one site out of 16), we were unable to assess model performance post-WNS using these data. The tricolored bat distribution model performed better than the northern long-eared bat model, however models for the tricolored bat were driven mainly by average summer temperature which provides limited opportunity for managementfor this species. With significant declines in populations because ofWNS, the forest-dwelling northern long-eared bat (Myotis septentrionalis) is of particular interest to wildlife managers seeking methods to aid recovering populations. Using number of bats captured and time-to-detection as proxies for bat abundance, we investigated relationships between abundance of northern long-eared bats and snag density metrics from 30 mist-net survey sites in four forests in Wisconsin in 2014 prior to arrival of WNS in the state. We used high resolution imagery to classify land cover surrounding netting sites and extracted number of snag patches, clumpiness index, mean size of snag patches, and closest snag patch within 1 km from netting sites. I used Cox Proportional Hazards time-to-event survival modeling and Poisson regression to assess relationships between bat abundance and measures of snag proximity and abundance. I captured northern long-eared bats at 13 of the 30 netting sites, and sites ranged from 0-3 bats captured. Using high resolution imagery to classify land cover maps wequantified snag metrics at a landscape scale and investigate the role that snag metrics play on local abundance of northern long-eared bats in Wisconsin. Number of snag patches within 1 km and area delineation were the best predictors of time-to-event (first capture) and the counts of northern long-eared bats captured. However, the effect of number of snag patches was small and confidence intervals overlapped zero suggesting an uninformative parameter or insufficient sample sizes in this modeling. The top model using Coxph regression of number of snag patches and area delineation had a 29%chance of being the best approximating model given the data, however with this low weight, there exists model selection uncertainty. Top model weight was higher using Poisson regression (wi = 0.36), although model selection uncertainty was still present.

Publications

  • Type: Theses/Dissertations Status: Accepted Year Published: 2019 Citation: Kaarakka, H. M., 2019. Summer distribution and habitat suitability for two threatened bat species in Wisconsin. Thesis. University of Wisconsin - Madison, USA.


Progress 10/01/18 to 09/30/19

Outputs
Target Audience:The target audience for this research includes scientists, natural resources managers, and conservationists who are concerned with the ecology and conservation of native bat species. Changes/Problems: Nothing Reported What opportunities for training and professional development has the project provided?This project has provided research training for a master's degree student in the department of Forest and Wildlife Ecology at the University of Wisconsin - Madison. How have the results been disseminated to communities of interest?The recently completed Master's Thesis has been distributed to managers in the Wisconsin DNR and the graduate student has presented her work at several professional conferences. Manuscripts for peer review publication are being drafted. What do you plan to do during the next reporting period to accomplish the goals?We intend to prepare manuscripts for publication and address issues of model validation and model sensitivity.

Impacts
What was accomplished under these goals? We used maximum entropy modeling (MaxEnt) to estimate distributions and determine land cover and climatic variables influencing detections of two cave-dwelling bat species in Wisconsin, both of whom depend on forest roosting sites during summer. We used spatially thinned occurrence data from June and July 2010 to 2017 to build species distribution models. Tricolored bats were restricted in Wisconsin in summer to regions south of Wisconsin's floristic tension zone and were influenced primarily by a positive relationship with average daily summer temperature. Northern long-eared bats were found statewide, but areas of high habitat suitability in central and southwest Wisconsin aligned with major river corridors. Important predictor variables affecting northern long-eared bats included increasing amount of lowland forest within 1 km, higher summer precipitation, and reduced distance from major hibernation sites. Results from acoustic sampling at 16 locations in southern and western Wisconsin in 2019 resulted in tricolored bat recorded at one site, and no northern long-eared bats were recorded.

Publications

  • Type: Theses/Dissertations Status: Accepted Year Published: 2019 Citation: Kaarakka, H. M. 2019. Summer distribution and habitat suitability for two threatened bat species in Wisconsin. Thesis. University of Wisconsin - Madison. 114pp.


Progress 10/04/17 to 09/30/18

Outputs
Target Audience:Target audience is state, federal, and private sector biologists who are concerned about the effects of white-nosed syndrome on native bats. Changes/Problems: Nothing Reported What opportunities for training and professional development has the project provided?Professional development included support for a masters student in Forestry and Wildlife. Graduate students attended2018 Midwest Fish and Wildlife Conference and the 2018 North American Society for Bat Research annual meeting. How have the results been disseminated to communities of interest? Nothing Reported What do you plan to do during the next reporting period to accomplish the goals?We plan to finalize habitat suitability models and develop validation methods for these models. Snag classification maps will be finalized and survival models will be developed to investigate snag density in relation to occupancy by one species of bat. Models will be developed investigating finer scale habitat use by bats while foraging on the summerlandscape. Student plans to present resultsat least one meeting in 2019.

Impacts
What was accomplished under these goals? In the first year of the grant, a thesis proposal was drafted and the graduate student worked with the Department of Natural Resources to secure existing data for analysis. We began developing habitat suitability models and maps for two species of bats, and snag classification maps for one species of bat. The two species of bat are the Northern Long-Eared Bat (Myotis septentrionalis) and the Eastern Pipistrelle (Perimyotis subflavus). Both are resident in the sense that they don't migrate but rather hibernate in mass in caves where the WNS fungus can grow and transmit. Work accomplished under the project included field sampling for bat occurrence, compiling of the cover type and landscape data needed to support the summer habitat models and initial model fitting with evaluation of fit statistics and model diagnostics. Habitat suitability models are essentially multivariate comparisons of survey data where bats were detected using trapping or acoustic sampling to areas where bats have not been detected The model is then applied to the landscape generally to identify or classify bat habitat in parts of the landscape where surveys have not occurred. Heather has been using the DNR's detection data to fit alternate models when landscape and cover type features are used as predictors. She is then using model selection techniques and model diagnostics to identify those models that are optimal given the data available and the sampling constraints. In the wake of white-nose syndrome's devastation of bat populations in Wisconsin and North America, identification of surviving bat populations and critical habitat is important for conservation efforts. Habitat suitability models offer a method of determining important climatic and landscape factors that influence where bats are likely to befound. Results from ourmodeling will guide future surveys for these two bats in Wisconsinand focus conservation efforts ofsummer habitat that iscritical for reproduction and continuation of these species.

Publications