Detecting and Tracking Communal Bird Roosts in Weather Radar Data

The US weather radar archive holds detailed information about biological phenomena in the atmosphere over the last 20 years. Communally roosting birds congregate in large numbers at nighttime roosting locations, and their morning exodus from the roost is often visible as a distinctive pattern in radar images. This paper describes a machine learning system to detect and track roost signatures in weather radar data. A significant challenge is that labels were collected opportunistically from previous research studies and there are systematic differences in labeling style. We contribute a latent variable model and EM algorithm to learn a detection model together with models of labeling styles for individual annotators. By properly accounting for these variations we learn a significantly more accurate detector. The resulting system detects previously unknown roosting locations and provides comprehensive spatio-temporal data about roosts across the US. This data will provide biologists important information about the poorly understood phenomena of broad-scale habitat use and movements of communally roosting birds during the non-breeding season.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 2

page 4

page 7

11/30/2021

Synthetic weather radar using hybrid quantum-classical machine learning

The availability of high-resolution weather radar images underpins effec...
04/25/2014

Improving weather radar by fusion and classification

In air traffic management (ATM) all necessary operations (tactical plani...
07/01/2019

Short-term prediction of Electricity Outages Caused by Convective Storms

Prediction of power outages caused by convective storms which are highly...
10/11/2020

Weather Radar in Nepal: Opportunities and Challenges in Mountainous Region

Extreme rainfall is one of the major causes of natural hazards (for exam...
07/15/2021

CMU-GPR Dataset: Ground Penetrating Radar Dataset for Robot Localization and Mapping

There has been exciting recent progress in using radar as a sensor for r...
11/14/2019

Convolutional Neural Network for Convective Storm Nowcasting Using 3D Doppler Weather Radar Data

Convective storms are one of the severe weather hazards found during the...
10/01/2019

MASS-UMAP: Fast and accurate analog ensemble search in weather radar archive

The use of analogs - similar weather patterns - for weather forecasting ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

The US weather radar network offers an unprecedented opportunity to study dynamic, continent-scale movements of animals over a long time period. The National Weather Service operates a network of 143 radars in the contiguous US. These radars were designed to study weather phenomena such as precipitation and severe storms. However, they are very sensitive and also detect flying animals, including birds, bats, and insects [17]. The data are archived and available from the early 1990s to present, and provide comprehensive views of a number of significant biological phenomena. These include broad-scale patterns of bird migration, such as the density and velocity of all nocturnally migrating birds flying over the US in the spring and fall [11, 26], and the number of migrants leaving different stopover habitats [5, 4]. Radars can also detect phenomena that can be matched to individual species, such as insect hatches and departures of large flocks of birds or bats from roosting locations [32, 16, 6].

Figure 1: (a) Illustration of roost exodus. (b) A radar traces out cone-shaped slices of the atmosphere (left), which are rendered as top-down images (center). This image from the Dover, DE radar station at 6:52 am on Oct 2, 2010 shows at least 8 roosts. Several are shown in more detail to the right, together with crops of one roost from five consecutive reflectivity and radial velocity images over a period of 39 minutes. These show the distinctive expanding ring and “red-white-green” diverging velocity patterns.

This paper is about the development of an AI system to track the roosting and migration behavior of communally roosting songbirds. Such birds congregate in large roosts at night during portions of the year. Roosts may contain thousands to millions of birds and are packed so densely that when birds depart in the morning they are visible on radar. Roosts occur all over the US, but few are documented, and many or most are never witnessed by humans, since birds enter and leave the roosts during twilight. Those who do witness the swarming behavior of huge flocks at roosts report an awe-inspiring spectacle111See: https://people.cs.umass.edu/~zezhoucheng/roosts.

We are particularly interested in swallow species, which, due to their unique behavior, create a distinctive expanding ring pattern in radar (see Fig. 1). During fall and winter, the Tree Swallow (Tachycineta bicolor) is the only swallow species present in North America, so roost signatures provide nearly unambiguous information about this single species, which is highly significant from an ecological perspective. Swallows are aerial insectivores, which are rapidly declining in North America [12, 22]. Information about the full life cycles of these birds is urgently needed to understand and reverse the causes of their declines [21]. Radar data represents an unparalleled opportunity to gather this information, but it is too difficult to access manually, which has greatly limited the scope (in terms of spatial/temporal extent or detail) of past research studies [18, 2].

We seek an AI system to fill this gap. An AI system has the potential to collect fine-grained measurements of swallows—across the continent, at a daily time scale—from the entire 20-year radar archive. This would allow scientists to study trends in population size and migration behavior in relation to changes in climate and habitat availability, and may provide some of the first insights into where and when birds die during their annual life cycle, which is critical information to guide conservation efforts.

The dramatic successes of deep learning for recognition tasks make it an excellent candidate for detecting and tracking roosts in radar. We develop a processing pipeline to extract useful biological information by solving several challenging sub-tasks. First, we develop a single-frame roost detector system based on Faster R-CNN, an established object detection framework for natural images 

[24]. Then, we develop a tracking system based on the “detect-then-track” paradigm [25] to assemble roosts into sequences to compute meaningful biological measures and improve detector performance. Finally, we use auxiliary information about precipitation and wind farms to reduce certain sources of false positives that are poorly represented in our training data. In the final system, 93.5% of the top 521 tracks on fresh test data are correctly identified and tracked bird roosts.

A significant challenge in our application is the presence of systematic differences in labeling style. Our training data was annotated by different researchers and naturalists using a public tool developed for prior research studies [18]. They had different backgrounds and goals, and roost appearance in radar is poorly understood, leading to considerable variability, much of which is specific to individual users. This variation makes evaluation using held out data very difficult and inhibits learning due to inconsistent supervision.

We contribute a novel approach to jointly learn a detector together with user-specific models of labeling style. We model the true label as a latent variable, and introduce a probabilistic user model for the observed label conditioned on the image, true label, and user. We present a variational EM learning algorithm that permits learning with only black-box access to an existing object detection model, such as Faster R-CNN. We show that accounting for user-specific labeling bias significantly improves evaluation and learning.

Finally, we conduct a case study by applying our models to detect roosts across the entire eastern US on a daily basis for fall and winter of 2013-2014 and cross-referencing roost locations with habitat maps. The detector has excellent precision, detects previously unknown roosts, and demonstrates the ability to generate urgently needed continent-scale information about the non-breeding season movements and habitat use of an aerial insectivore. Our data points to the importance of the eastern seaboard and Mississippi valley as migration corridors, and suggests that birds rely heavily on croplands (e.g., corn and sugar cane) earlier in the fall prior to harvest.

2 A System to Detect and Track Roosts

Radar Data

We use radar data from the US NEXRAD network of over 140 radars operated by the National Weather Service [8]. They have ranges of several hundred kilometers and cover nearly the entire US. Data is available from the 1990s to present in the form of raster data products summarizing the results of radar volume scans, during which a radar scans the surrounding airspace by rotating the antenna 360 at different elevation angles (e.g., 0.5, 1.5) to sample a cone-shaped “slice” of airspace (Fig. 1b). Radar scans are available every 4–10 minutes at each station. Conventional radar images are top-down views of these sweeps; we will also render data this way for processing.

Standard radar scans collect 3 data products at 5 elevation angles, for 15 total channels. We focus on data products that are most relevant for detecting roosts. Reflectivity is the base measurement of the density of objects in the atmosphere. Radial velocity uses the Doppler shift of the returned signal to measure the speed at which objects are approaching or departing the radar. Copolar cross-correlation coefficient is a newer data product, available since 2013, that is useful for discriminating rain from biology [27]. We use it for post-processing, but not training, since most of our labels are from before 2013.

Roosts

A roost exodus (Fig. 1a) is the mass departure of a large flock of birds from a nighttime roosting location. They occur 15–30 minutes before sunrise and are very rarely witnessed by humans. However, roost signatures are visible on radar as birds fly upward and outward into the radar domain. Fig. 1b, center, shows a radar reflectivity image with at least 8 roost signatures in a area. Swallow roosts, in particular, have a characteristic signature shown in Fig. 1b, right. The center row shows reflectivity images of one roost expanding over time. The bottom row shows the characteristic radial velocity pattern of birds dispersing away from the center of the roost. Birds moving toward the radar station (bottom left) have negative radial velocity (green) and birds moving away from the radar station (top right) have positive radial velocity (red).

Annotations

We obtained a data set of manually annotated roosts collected for prior ecological research [18]. They are believed to be nearly 100% Tree Swallow roosts. Each label records the position and radius of a circle within a radar image that best approximates the roost. We restricted to seven stations in the eastern US and to month-long periods that were exhaustively labeled, so we could infer the absence of roosts in scans with no labels. We restricted to scans from 30 minutes before to 90 minutes after sunrise, leading to a data set of 63691 labeled roosts in 88972 radar scans. A significant issue with this data set is systematic differences in labeling style by different researchers. This poses serious challenges to building and evaluating a detection model. We discuss this further and present a methodological solution in Sec. 3.

Related Work on Roost Detection and Tracking

There is a long history to the study of roosting behavior with the radar data, almost entirely based on human interpretation of images [32, 18, 2]

. That work is therefore restricted to analyze only limited regions, short-time periods or coarse-grained information about the roosts. chilson2019 [chilson2019] developed a deep-learning image classifier to identify radar images that contain roosts. While useful, this provides only limited biological information. Our method locates roosts within images and tracks them across frames, which is a substantially harder problem, and important biologically. For example, to create population estimates or locate roost locations on the ground, one needs to know where roosts occur within the radar image; to study bird behavior, one needs to measure how roosts move and expand from frame to frame. Our work is the first machine-learning method able to extract the type of higher-level biological information presented in our case study (Sec. 

5). Our case study illustrates, to the best of our knowledge, the first continent-scale remotely sensed observations of the migration of a single bird species.

Methodology

Our overall approach consists of four steps (see Fig. 2): we render radar scans as multi-channel images, run a single-frame detector, assemble and rescore tracks, and then postprocess detections using other geospatial data to filter specific sources of false positives.

Figure 2: Detection and tracking pipeline. A final step (not shown) uses auxiliary data to filter rain and wind-farms.

Detection Architecture

Our single-frame detector is based on Faster R-CNNs [24]. Region-based CNN detectors such as Faster R-CNNs are state-of-the-art on several object detection benchmarks.

A significant advantage of these architectures comes from pretraining parts of the network on large labeled image datasets such as ImageNet 

[9]. To make radar data compatible with these networks, we must select only 3 of the 15 available channels to feed into the RGB-based models. We select the radar products that are most discriminative for humans: reflectivity at 0.5, radial velocity at 0.5 degrees, and reflectivity at 1.5. Roosts appear predominantly in the lowest elevations and are distinguished by the ring pattern in reflectivity images and distinctive velocity pattern. These three data products are then rendered as a image in the “top-down” Cartesian-coordinate view (out to 150km from the radar station) resulting in a 3-channel image. The three channel images are fed into Faster R-CNN initialized with a pretrained VGG-M network [7]. All detectors are trained for the single “roost” object class, using bounding boxes derived from the labeled dataset described above.

Although radar data is visually different from natural images, we found ImageNet pretraining to be quite useful; without pretraining the networks took significantly longer to converge and resulted in a 15% lower performance. We also experimented with models that map 15 radar channels down to 3 using a learned transformation. These networks were not consistently better than ones using hand-selected channels. Models trained with shallower networks that mimic handcrafted features, such as those based on gradient histograms, performed 15-20% worse depending on the architecture. See the supplementary material on the project page for details on these baseline detection models.

We defer training details to Sec. 3, where we discuss our approach of jointly learning the Faster R-CNN detector together with user-models for labeling style.

Roost Tracking and Rescoring

Associating and tracking detections across frames is important for several reasons. It helps rule out false detections due to rain and other phenomenon that have different temporal properties than roosts (see Sec. 5

). Detection tracks are also what associate directly to the biological entity—a single flock of birds—so they are needed to estimate biological parameters such as roost size, rate of expansion, location and habitat of first appearance, etc. We employ a greedy heuristic to assemble detections from individual frames into tracks 

[25], starting with high scoring detections and incrementally adding unmatched detections with high overlap in nearby frames. Detections that match multiple tracks are assigned to the longest one. After associating detections we apply a Kalman smoother to each track using a linear dynamical system model for the bounding box center and radius. This model captures the dynamics of roost formation and growth with parameters estimated from ground-truth annotations. We then conduct a final rescoring step where track-level features (e.g., number of frames, average detection score of all bounding boxes in track) are associated to individual detections, which are then rescored using a linear SVM. This step suppresses false positives that appear roost-like in single frames but do not behave like roosts. Further details of association, tracking and rescoring can be found in the supplementary material.

Postprocessing with Auxiliary Information

In preliminary experiments, the majority of high-scoring tracks were roosts, but there were also a significant number of high-scoring false positives caused by specific phenomena, especially wind farms and precipitation (see Sec. 5). We found it was possible to reliably reject these false positives using auxiliary information. To eliminate rain in modern data, we use the radar measurement of copolar cross-correlation coefficient, , which is available since 2013 [27]. Biological targets have much lower

values than precipitation due to their high variance in orientation, position and shape over time. A common rule is to classify pixels as rain if

 [10]. We classify a roost detection as precipitation if a majority of pixels inside its bounding box have . For historic data one may use automatic methods for segmenting precipitation in radar images such as [19]. For wind farms, we can use recorded turbine locations from the U.S. Wind Turbine Database [15]. A detection is identified as a wind farm if any turbine from the database is located inside its bounding box.

3 Modeling Labeling Styles

Figure 3: Labeling style variation leads to inaccurate evaluation and suboptimal detectors. All of these detections (pink boxes) are misidentified as false positives because of insufficient overlap with annotations of one user (green boxes) with a tight labeling style. Label variation also hurts training and leads to suboptimal models.

Preliminary experiments revealed that systematic variations in labeling style were a significant barrier to training and evaluating a detector. Fig. 3 shows example detections that correctly locate and circumscribe roosts, but are classified as false positives because the annotator used labels (originally circles) to “trace” roosts instead of circumscribing them. Although it is clear upon inspection that these detections are “correct”, with 63691 labels and a range of labeling styles, there is no simple adjustment to accurately judge performance of a system. Furthermore, labeling variation also inhibits learning and leads to suboptimal models. This motivates our approach to jointly learn a detector along with user-specific models of labeling style. Our goal is a generic and principled approach that can leverage existing detection frameworks with little or no modification.

Latent Variable Model and Variational EM Algorithm

To model variability due to annotation styles we use the following graphical model:

where is the image, represents the unobserved “true” or gold-standard label, is the user (or features thereof), and is the observed label in user ’s labeling style. In this model

  • [leftmargin=*, itemsep=0pt]

  • is the detection model, with parameters

    . We generally assume the negative log-likelihood of the detection model is equal to the loss function of the base detector. For example, in our application,

    , the loss function of Faster R-CNN.222Faster R-CNN includes a region proposal network to detect and localize candidate objects and a classification network to assign class labels. The networks share parameters and are trained jointly to minimize a sum of several loss functions; we take the set of all parameters as and the sum of loss functions as .

  • is the forward user model for the labeling style of user , with parameters . In our application, much of the variability can be captured by user-specific scaling of the bounding boxes, so we adopt the following user model: for each bounding-box we model the observed radius as where is the unobserved true radius and is the user-specific scaling factor. In this model, the bounding-box centers are unmodified and the user model does not depend on the image , even though our more general framework allows both.

  • is the reverse user model. It is determined by the previous two models, and is needed to reason about the true labels given the noisy ones during training. Since this distribution is generally intractable, we use instead a variational reverse user model , with parameters . In our application, , which is another user-specific rescaling of the radius.

We train the user models jointly with Faster R-CNN using variational EM. We initialize the Faster R-CNN parameters by training for 50K iterations starting from the ImageNet pretrained VGG-M model using the original uncorrected labels. We then initialize the forward user model parameters using the Faster R-CNN predictions: if a predicted roost with radius overlaps sufficiently with a labeled roost (intersection-over-union > 0.2) and has high enough detection score (> 0.9), we generate a training pair where

is the labeled radius. We then estimate the forward regression model parameters as a standard linear regression with these pairs.

After initialization, we repeat the following steps (in which is an index for annotations):

  1. [leftmargin=*,itemsep=0pt]

  2. Update parameters of the reverse user model by minimizing the combined loss . The optimization is performed separately to determine the reverse scaling factor for each user using Brent’s method with search boundary and black-box access to .

  3. Resample annotations on the training set by sampling for all , then update by training Faster R-CNN for 50K iterations using the resampled annotations.

  4. Update by training the forward user models using pairs , where

    is the radius of the imputed label.

Formally, each step can be justified as maximizing the evidence lower bound (ELBO) [1] of the log marginal likelihood with respect to the variational distribution . Steps 1, 2, and 3 maximize the ELBO with respect to , , and , respectively. Steps 1 and 2 require samples from the reverse user model; we found that using the maximum a posteriori instead of sampling is simple and performs well in practice, so we used this in our application.

The derivation is presented in the supplementary material. It assumes that is a structured label that includes all bounding boxes for an image. This justifies equating with the loss function of an existing detection framework that predicts bounding boxes simultaneously for an entire image (e.g., using heuristics like non-maximum suppression). This is important because it is modular. We can use any detection framework that provides a loss function, with no other changes. A typical user model will then act on (a set of bounding boxes) by acting independently on each of its components, as in our application.

We anticipate this framework can be applied to a range of applications. More sophisticated user models may also depend on the image to capture different labeling biases, such as different thresholds for labeling objects, or tendencies to mislabel objects of a certain class or appearance. However, it is an open question how to design more complex user models and we caution about the possibility of very complex user models “explaining away” true patterns in the data.

Related Work on Labeling Style and Noise

jiang2017face [jiang2017face] discuss how systematic differences in labeling style across face-detection benchmarks significantly complicate evaluation, and propose fine-tuning techniques for style adaptation. Our EM approach is a good candidate to unify training and evaluation across these different benchmarks. Prior research on label noise [13] has observed that noisy labels degrade classifier performance [23] and proposed various methods to deal with noisy labels [30, 14, 3, 33, 28], including EM [20]. While some considerations are similar (degradation of training performance, latent variable models), labeling style is qualitatively different in that an explanatory variable (the user) is available to help model systematic label variation, as opposed to pure “noise”. Also, the prior work in this literature is for classification. Our approach is the first noise-correction method for bounding boxes in object detection.

4 Experiments

Dataset

We divided the 88972 radar scans from the manually labeled dataset (Sec. 2) into training, validation, and test sets. Tab. 1 gives details of training and test data by station. The validation set (not shown) is roughly half the size of the test set and was used to set the hyper-parameters of the detector and the tracker.

Evaluation Metric

To evaluate the detector we use established evaluation metrics for object detection employed in common computer vision benchmarks. A detection is a true positive if its overlap with an annotated bounding-box, measured using the intersection-over-union (IoU) metric, is greater than 0.5. The mean average precision (MAP) is computed as the area under the precision-recall curve. For the purposes of evaluating the detector we mark roosts smaller than

in a radar image as difficult and ignore them during evaluation. Humans typically detect such roosts by looking at adjacent frames. As discussed previously (Fig. 3), evaluation is unreliable when user labels have different labeling styles. To address this we propose an evaluation metric (“+User”) that rescales predictions on a per-user basis prior to computing MAP. Scaling factors are estimated following the same procedure used to initialize variational EM. This assumes that the user information is known for the test set, where it is only used for rescaling predictions, and not by the detector.

Results: Roost Detector and User Model

Tab. 1 shows the performance of various detectors across radar stations. We trained two detector variants, one a standard Faster R-CNN, and another trained with the variational EM algorithm. We evaluated the detectors based on whether annotation bias was accounted for during testing (Tab. 1, “+User”).

The noisy annotations cause inaccurate evaluation. A large number of the detections on KDOX are misidentified as negatives because of the low overlap with the annotations, which are illustrated in Fig. 3, leading to a low MAP score of 9.1%. This improves to 44.8% when the annotation biases are accounted for during testing. As a sanity check we trained and evaluated a detector on annotations of a single user on KDOX and found its performance to be in the mid fifties. However, the score was low when this model was evaluated on annotations from other users or stations.

The detector trained jointly with user-models using variational EM further improves performance across all stations (Tab. 1, “+EM+User”), with larger improvements for stations with less training data. Overall MAP improves from 44.2% to 45.5%. To verify the statistical significance of this result, we drew 20 sets of bootstrap resamples from the entire test set (contains

images), computed the MAP of the model trained with EM and without EM on each set. The mean and standard error of MAPs for the model trained with EM are 45.5% and 0.12% respectively, while they are 44.4% and 0.11% for the model trained without EM.

Station Test Train R-CNN +User +EM+User
KMLB 9133 19998 47.5 47.8 49.2
KTBW 7195 16382 47.3 50.0 50.8
KLIX 4077 10192 32.4 35.1 35.7
KOKX 1404 2994 23.2 27.3 29.9
KAMX 860 1898 29.9 30.8 31.6
KDOX 639 902 9.1 44.8 50.2
KLCH 112 441 32.1 39.8 43.1
entire 23.7k 53.6k 41.0 44.2 45.5
Table 1: Roost detection MAP for detector variants.

Results: Tracking and Rescoring

After obtaining the roost detections from our single-frame detector, we can apply our roost tracking model to establish roost tracks over time. Fig. 4 shows an example radar sequence where roost detections have been successfully tracked over time and some false positives removed. We also systematically evaluated the tracking and rescoring model on scans from the KOKX station. For this study we performed a manual evaluation of the top 800 detections before and after the contextual rescoring. Manual evaluation was necessary due to human labeling biases, especially the omission of labels at the beginning or end of a roost sequence when roost signatures are not as obvious. Fig. 4, middle panel, shows that the tracking and rescoring improves the precision across the entire range of . Our tracking model also enables us to study the roost dynamics over time (see Sec. 5 and Fig. 4 right panel).

Figure 4: Left: tracking example, with raw detections (top) and track (bottom). Transient false positives in several frames lead to poor tracks and are removed by the rescoring step. Middle: precision@k before and after rescoring. Right: Roost radius relative to time after sunrise.

5 Case study

Figure 5: Tree Swallow fall migration in 2013. The color circles show detected roost locations with each half-month period. The location of each roost is determined by the center of the first bounding box in the track, when the airborne birds are closest to their location on the ground. Faint gray triangles show radar station locations.

We conducted a case study to use our detector and tracker to synthesize knowledge about continent-scale movement patterns of swallows. We applied our pipeline to radar scans collected from 86 radar stations in the Eastern US (see Figure 5) from October 2013 through March 2014. During these months, Tree Swallows are the only (fall/winter) or predominant (early spring) swallow species in the US and responsible for the vast majority of radar roost signatures. This case study is therefore the first system to obtain comprehensive measurements of a single species of bird across its range on a daily basis. We ran our detector and tracking pipeline on all radar scans from 30 minutes before sunrise to 90 minutes after sunrise. We kept tracks having at least two detections with a detector score of 0.7 or more, and then ranked tracks by the sum of the detector score for each detection in the track.

Error Analysis

There were several specific phenomena that were frequently detected as false positives prior to post-processing. We reviewed and classified all tracks with a total detection score of 5 or more prior to postprocessing (678 tracks total) to evaluate detector performance “in the wild” and the effectiveness of post-processing. This also served to vet the final data used in the biological analysis. Tab. 2 shows the number of detections by category before and after post-processing. Roughly two-thirds of initial high scoring detections were swallow roosts, with another 5.6% being communal roosts of some bird species.

Pre Post Pre Post
Swallow roost 454 449 Other roost 38 38
Precipitation 109 5 Clutter 22 21
Wind farm 47 0 Unknown 8 8
Table 2: Detections by type pre- and post-filtering with auxiliary data. Post-processing effectively removes false positives due to precipitation and wind farms.

The most false positives were due to precipitation, which appears as highly complex and variable patterns in radar images, so it is common to find small image patches that share the general shape and velocity pattern of roosts (Fig. 6, fourth column). Humans recognize precipitation from larger-scale patterns and movement. Filtering using nearly eliminates rain false positives. The second leading source of false positives was wind farms. Surprisingly, these share several features of roosts: they appear as small high-reflectivity “blobs” and have a diverse velocity field due to spinning turbine blades (Fig. 6 last column). Humans can easily distinguish wind farms from roosts using temporal properties. All wind farms are filtered successfully using the wind turbine database. Since our case study focuses on Tree Swallows, we marked as “other roost” detections that were believed to be from other communally roosting species (e.g. American Robins, blackbirds, crows). These appear in radar less frequently and with different appearance (usually “blobs” instead of “rings”; Fig. 6, fifth column) due to behavior differences. Humans use appearance cues as well as habitat, region, and time of year to judge the likely species of a roost. We marked uncertain cases as “other roost”.

Figure 6: Some detections visualized on the reflectivity (top) and radial velocity (bottom) channels of different scans. The first three columns show swallow roost detections while the next three columns show detections due to rain, roosts of other species, and windmills.

Migration and Habitat Use

Fig. 5 shows swallow roost locations and habitat types for five half-month periods starting in October to illustrate the migration patterns and seasonal habitat use of Tree Swallows. Habitat assignments are based on the majority habitat class from the National Land Cover Database (NLCD) [29] within a area surrounding the roost center, following the approach of [2] for Purple Martins. Unlike Purple Martins, the dominant habitat type for Tree Swallows is wetlands (38% of all roosts), followed by croplands (29%). These reflect the known habits of Tree Swallows to roost in reedy vegetation—either natural wetlands (e.g. cattails and phragmites) or agricultural fields (e.g., corn, sugar cane) [31].

In early October, Tree Swallows have left their breeding territories and formed migratory and pre-migratory roosts throughout their breeding range across the northern US [31]. Agricultural roosts are widespread in the upper midwest. Some birds have begun their southbound migration, which is evident by the presence of roosts along the Gulf Coast, which is outside the breeding range. In late October, roosts concentrate along the eastern seaboard (mostly wetland habitat) and in the central US (mostly cropland). Most of the central US roosts occur near major rivers (e.g., the Mississippi) or other water bodies. The line of wetland roosts along the eastern seaboard likely delineates a migration route followed by a large number of individuals who make daily “hops” from roost to roost along this route [32]. By early November, only a few roosts linger near major water bodies in the central US. Some birds have left the US entirely to points farther south, while some remain in staging areas along the Gulf Coast [18]. By December, Gulf Coast activity has diminished, and roosts concentrate more in Florida, where a population of Tree Swallows will spend the entire winter.

Widespread statistics of roost locations and habitat usage throughout a migratory season have not previously been documented, but are enabled by our AI system to automatically detect and track roosts. Our results are a starting point to better understand and conserve these populations. They highlight the importance of the eastern seaboard and Mississippi valley as migration corridors, with different patterns of habitat use (wetland vs. agricultural) in each. The strong association with agricultural habitats during the harvest season suggests interesting potential interactions between humans and the migration strategy of swallows.

Roost Emergence Dynamics

Our AI system also enables us to collect more detailed information about roosts than previously possible, such as their dynamics over time, to answer questions about their behavior. Fig. 4 shows the roost radius relative to time after sunrise for roosts detected by our system. Roosts appear around 1000 seconds before sunrise and expand at a fairly consistent rate. The best fit line corresponds to swallows dispersing from the center of the roost with an average airspeed velocity of (unladen).

6 Conclusion

We presented a pipeline for detecting communal bird roosts using weather radar. We showed that user-specific label noise is a significant hurdle to doing machine learning with the available data set, and presented a method to overcome this. Our approach reveals new insights into the continental-scale roosting behavior of migratory Tree Swallows, and can be built upon to conduct historical analysis using 20+ years of archived radar data to study their long-term population patterns in comparison with climate and land use change.

Acknowledgements

This research was supported in part by NSF #1749833, #1749854, #1661259 and the MassTech Collaborative for funding the UMass GPU cluster.

References

  • [1] D. M. Blei, A. Kucukelbir, and J. D. McAuliffe (2017) Variational inference: a review for statisticians. Journal of the American Statistical Association 112 (518). Cited by: §3.
  • [2] E. S. Bridge, S. M. Pletschet, T. Fagin, P. B. Chilson, K. G. Horton, K. R. Broadfoot, and J. F. Kelly (2016) Persistence and habitat associations of purple martin roosts quantified via weather surveillance radar. Landscape ecology 31 (1). Cited by: §1, §2, §5.
  • [3] C. E. Brodley and M. A. Friedl (1999) Identifying mislabeled training data.

    Journal of artificial intelligence research

    .
    Cited by: §3.
  • [4] J. J. Buler and D. K. Dawson (2014) Radar analysis of fall bird migration stopover sites in the northeastern us. The Condor 116 (3). Cited by: §1.
  • [5] J. J. Buler and R. H. Diehl (2009) Quantifying bird density during migratory stopover using weather surveillance radar. IEEE Transactions on Geoscience and Remote Sensing 47 (8). Cited by: §1.
  • [6] J. J. Buler, L. A. Randall, J. P. Fleskes, W. C. Barrow Jr, T. Bogart, and D. Kluver (2012) Mapping wintering waterfowl distributions using weather surveillance radar. PloS one 7 (7). Cited by: §1.
  • [7] K. Chatfield, K. Simonyan, A. Vedaldi, and A. Zisserman (2014) Return of the devil in the details: delving deep into convolutional nets. arXiv preprint arXiv:1405.3531. Cited by: §2.
  • [8] T. D. Crum and R. L. Alberty (1993) The WSR-88D and the WSR-88D operational support facility. Bulletin of the American Meteorological Society 74 (9). Cited by: §2.
  • [9] J. Deng, W. Dong, R. Socher, L. Li, K. Li, and L. Fei-Fei (2009) ImageNet: A large-scale hierarchical image database. In CVPR, Cited by: §2.
  • [10] A. M. Dokter, P. Desmet, J. H. Spaaks, S. van Hoey, L. Veen, L. Verlinden, C. Nilsson, G. Haase, H. Leijnse, A. Farnsworth, W. Bouten, and J. Shamoun-Baranes (2018) bioRad: biological analysis and visualization of weather radar data. Ecography. Cited by: §2.
  • [11] A. Farnsworth, B. M. Van Doren, W. M. Hochachka, D. Sheldon, K. Winner, J. Irvine, J. Geevarghese, and S. Kelling (2016) A characterization of autumn nocturnal migration detected by weather surveillance radars in the northeastern USA. Ecological Applications 26 (3). External Links: ISSN 1939-5582 Cited by: §1.
  • [12] K. C. Fraser, B. J. M. Stutchbury, C. Silverio, P. M. Kramer, J. Barrow, D. Newstead, N. Mickle, B. F. Cousens, J. C. Lee, D. M. Morrison, T. Shaheen, P. Mammenga, K. Applegate, and J. Tautin (2012) Continent-wide tracking to determine migratory connectivity and tropical habitat associations of a declining aerial insectivore. Proceedings of the Royal Society of London B: Biological Sciences. External Links: Document, ISSN 0962-8452, Link, http://rspb.royalsocietypublishing.org/content/early/2012/10/22/rspb.2012.2207.full.pdf Cited by: §1.
  • [13] B. Frénay and M. Verleysen (2014) Classification in the presence of label noise: a survey.

    IEEE transactions on neural networks and learning systems

    25 (5), pp. 845–869.
    Cited by: §3.
  • [14] A. Ghosh, H. Kumar, and P. Sastry (2017) Robust loss functions under label noise for deep neural networks.. In AAAI, pp. 1919–1925. Cited by: §3.
  • [15] B. D. Hoen, J. E. Diffendorfer, J. T. Rand, L. A. Kramer, C. P. Garrity, and H.E. Hunt (2019) United states wind turbine database. U.S. Geological Survey, American Wind Energy Association, and Lawrence Berkeley National Laboratory data release: USWTDB V1.3.. Note: https://eerscmap.usgs.gov/uswtdb Cited by: §2.
  • [16] J. W. Horn and T. H. Kunz (2008) Analyzing nexrad doppler radar images to assess nightly dispersal patterns and population trends in brazilian free-tailed bats (tadarida brasiliensis). Integrative and Comparative Biology 48 (1). Cited by: §1.
  • [17] T. H. Kunz, S. A. Gauthreaux Jr, N. I. Hristov, J. W. Horn, G. Jones, E. K. Kalko, R. P. Larkin, G. F. McCracken, S. M. Swartz, R. B. Srygley, et al. (2008) Aeroecology: probing and modeling the aerosphere. Integrative and comparative biology 48 (1), pp. 1–11. Cited by: §1.
  • [18] A. J. Laughlin, D. R. Sheldon, D. W. Winkler, and C. M. Taylor (2016-10) Quantifying non-breeding season occupancy patterns and the timing and drivers of autumn migration for a migratory songbird using doppler radar. Ecography 39 (10), pp. 1017–1024. External Links: ISSN 1600-0587, Document, Link Cited by: §1, §1, §2, §2, §5.
  • [19] T. Lin, K. Winner, G. Bernstein, A. Mittal, A. M. Dokter, K. G. Horton, C. Nilsson, B. M. Van Doren, A. Farnsworth, F. A. La Sorte, S. Maji, and D. Sheldon (2019)

    MistNet: Measuring historical bird migration in the US using archived weather radar data and convolutional neural networks

    .
    Methods in Ecology and Evolution 10 (11), pp. 1908–1922. External Links: Document Cited by: §2.
  • [20] V. Mnih and G. E. Hinton (2012) Learning to label aerial images from noisy data. In ICML, Cited by: §3.
  • [21] NABCI (2012) The State of Canada’s Birds, 2012. Environment Canada, Ottawa, Canada. Cited by: §1.
  • [22] S. Nebel, A. Mills, J. D. McCracken, and P. D. Taylor (2010) Declines of aerial insectivores in North America follow a geographic gradient. Avian Conservation and Ecology. Cited by: §1.
  • [23] D. F. Nettleton, A. Orriols-Puig, and A. Fornells (2010)

    A study of the effect of different types of noise on the precision of supervised learning techniques

    .
    Artificial intelligence review. Cited by: §3.
  • [24] S. Ren, K. He, R. Girshick, and J. Sun (2015) Faster R-CNN: Towards real-time object detection with region proposal networks. In NIPS, Cited by: §1, §2.
  • [25] X. Ren (2008) Finding people in archive films through tracking. In CVPR, Cited by: §1, §2.
  • [26] J. Shamoun-Baranes, A. Farnsworth, B. Aelterman, J. A. Alves, K. Azijn, G. Bernstein, S. Branco, P. Desmet, A. M. Dokter, K. Horton, S. Kelling, J. F. Kelly, H. Leijnse, J. Rong, D. Sheldon, W. Van den Broeck, J. K. Van Den Meersche, B. M. Van Doren, and H. van Gasteren (2016) Innovative visualizations shed light on avian nocturnal migration. PLoS ONE 11 (8). Cited by: §1.
  • [27] P. M. Stepanian, K. G. Horton, V. M. Melnikov, D. S. Zrnić, and S. A. Gauthreaux (2016) Dual-polarization radar products for biological applications. Ecosphere 7 (11). Cited by: §2, §2.
  • [28] D. Tanaka, D. Ikami, T. Yamasaki, and K. Aizawa (2018) Joint optimization framework for learning with noisy labels. In CVPR, Cited by: §3.
  • [29] N. USGS (2011) Land cover (2011 edition, amended 2014), national geospatial data asset (ngda) land use land cover, 2011, editor. 2011. US Geological Survey. Cited by: §5.
  • [30] B. Van Rooyen, A. Menon, and R. C. Williamson (2015) Learning with symmetric label noise: the importance of being unhinged. In NIPS, Cited by: §3.
  • [31] D. W. Winkler, K. K. Hallinger, D. R. Ardia, R. J. Robertson, B. J. Stutchbury, and R. R. Cohen (2011) Tree Swallow (Tachycineta bicolor), version 2.0. In The Birds of North America, P. G. Rodewald (Ed.), Cited by: §5, §5.
  • [32] D. W. Winkler (2006) Roosts and migrations of swallows. Hornero 21 (2), pp. 85–97. Cited by: §1, §2, §5.
  • [33] T. Xiao, T. Xia, Y. Yang, C. Huang, and X. Wang (2015) Learning from massive noisy labeled data for image classification. In CVPR, Cited by: §3.