Study site
We monitored nest boxes installed by the Louisiana Department of Wildlife and Fisheries (LDWF) primarily for WODU use, although Hooded Mergansers (Lophodytes cucullatus) and BBWD frequently nest in them. These nest boxes were located in Iberville Parish, Louisiana, on the Atchafalaya National Wildlife Refuge at two portions of Sherburne Wildlife Management Area known as “North Farm” and “South Farm” (30.4166, -91.5186), which are managed as moist soil impoundments for waterfowl (Fig. 1). South Farm is open to the public, while North Farm is closed except for youth lottery hunts. Nest boxes were located on the side of levees and accessed by all-terrain vehicles or trucks. Sherburne is one constituent site of a larger research project on cavity-nesting ducks; nests have been monitored weekly and ducks have been captured/marked in boxes since 2020.
Historically, all nest boxes at our site were single units mounted on a pole. To increase sample size and commensurate with other LDWF wood duck monitoring sites, in February 2021 we converted all units to duplexes (two boxes on either side of the pole). Thus, for our study in 2022, we monitored 40 duplex-style nest boxes. We categorized nest boxes erected prior to 2021 as “old” and nest boxes erected in February 2021 as “new;” new nest boxes were identical to the existing (old) boxes (Fig. 2). 20 nest boxes were located at South Farm and 20 at North Farm; all duplexes were outfitted with a conical baffle predator guard (Fig. 2).
Nest monitoring
We visited nest boxes at approximately seven-day intervals. We numbered each egg with a permanent marker in the order it appeared in the nest [17], determined the species of each egg [4], and determined whether nests were active or inactive during each visit. We classified a nest as active if we observed a bird incubating, new eggs were laid since the previous visit, or incubation progressed [20]. We classified inactive nests as abandoned, depredated, or successful. We considered a nest abandoned if egg laying or incubation were discontinued without sign of depredation. We classified a nest as depredated if egg laying or incubation ceased and eggs were destroyed or missing. Successful nests survived to hatch ≥ 1 egg; we counted the number of unhatched eggs and egg membranes to determine the number of eggs that hatched. We also categorized each nest based on the apparent presence or absence of conspecific brood parasitism (CBP). CBP parasitized nests received > 1 egg per day during the laying stage [21] and/or additional eggs following day 4 of incubation, as our species initiate incubation about 3 days prior to laying the last egg [13]. We also identified interspecific brood parasitism (IBP); IBP nests received ≥ 1 egg of a different species.
During nest visits, we attempted to capture the incubating individual by covering the nest box entrance and removing the duck from the box through the side door. When caught, we checked the individual for any markers, including passive integrated transponder (PIT) tags and leg bands. We detected PIT tags in the hand using a handheld reader running the device along the dorsal side of the bird several times. We banded individuals using the appropriate-sized federal band for their species and embedded 2×12 mm, 125kHz PIT tags subcutaneously between the scapulae with a 12-gauge needle and implanter syringe; we closed the injection site with surgical adhesive.
As the hatch date approached, we checked each nest every 1–3 days to ensure we were there to process ducklings because they typically leave the nest within 24 hours of hatching [1]. For each duckling, we determined their species by plumage and bill shape (Fig. 3) and marked them with PIT tags in the same way as adults.
We collected data under U.S. Fish and Wildlife Service banding permit #06669 and Special Use Permit 43612-20-04; LDWF state collecting permits WDP-20-037 and WDP-21-060, and Wildlife Management Area Permit WL-Research-2020-03; protocols were approved by the Louisiana State University Institutional Animal Care and Use Committee Protocol A2019-27. This study complies with the ARRIVE guidelines for observational research.
RFID readers
We created stationary PIT tag readers using a custom radio frequency identification (RFID) circuit board [22]. Each unit was equipped with two loop-style antennas, so that one circuit board could be used on a duplex-style nest box with antennas encircling the nest box entrances. We configured the units as a simple data logger that recorded the date, time, and alphanumeric ID of each individual’s PIT tag as they entered or left the nest box.
We created the antennas using 26-gauge copper magnet wire and a circular jig to ensure the wire was wrapped in a uniform circle with an 11 cm diameter. PIT tags were particularly sensitive to antenna inductance, and deviation from optimal inductance resulted in failure of the antennas to detect PIT tags. We tuned each antenna to 1.2 mH using a digital inductance meter; each antenna required 67–68 turns of the copper wire to reach that inductance. We found that different brands of wire varied in the number of turns needed to achieve 1.2 mH, so each brand was tuned for the number of turns needed before creating a batch of antennas. Once the antenna wire was coiled, we wrapped it snugly with electrical tape (this step was important to keep consistent, as the tightness of the wire coil can change the antenna’s inductance) and coated it with Plasti-Dip, leaving the free ends (hereafter, leads) of the antenna exposed. Then, we used a fine-grit sandpaper to remove the insulation from the last ~ 1 cm of the antenna leads so that they could be connected to the circuit board. We attached the antenna leads to the circuit board’s screw clamps with a small flathead screwdriver and covered any exposed areas of the wire that remained with electrical tape to prevent them from touching one another. We marked “antenna 1” with a small piece of colored tape so the two could be easily distinguished when installing the readers in the field (Fig. 4).
After attaching the antennas to the circuit board, we installed a 3V CR1025 battery to the back of the circuit board as a backup clock battery and loaded firmware to each board using the Arduino Integrated Development Environment (Arduino LLC, Scarmagno, Italy). While installing the firmware, we ran a PIT tag through each antenna to ensure they both worked and assigned each circuit board a numeric ID that is displayed as the filename of data recorded to the SD card.
To prepare each unit for the field, we first ran the solar panel cable through a waterproof container so that the circuit board was protected from water while the solar panel was mounted on top of the nest box. We used dry boxes from Outdoor Products initially but switched to Ziploc Twist ‘n Loc containers due to cost and availability. For both container types, we drilled a hole to run the solar panel cable through and sealed it using waterproof epoxy putty. We then connected the cable to the battery pack (Voltaic V25 USB battery pack) using an adapter, connected the battery pack to the circuit board using a mini-USB cable, and inserted a FAT-formatted SD card into the circuit board. We ran the antenna lead wires through the lids without drilling a hole through the side of the container (Fig. 4).
We installed each unit on duplex-style nest boxes in the field, using large zip ties to attach the waterproof container to the pole underneath the nest boxes. To attach the antennas, we drilled 4 holes around the next box entrance and zip tied the antenna around it (Fig. 2). We then attached the cable extending from the waterproof container to the solar panel (Voltaic 3.5-watt/6 volt) and attached it to the top of the nest box using wood screws. We installed each solar panel facing south for optimal sun exposure (Fig. 2). After installation, we turned on the circuit board and confirmed that it was working properly. During weekly nest visits, we changed the SD cards in each unit and confirmed that the unit was still working properly.
Post-processing and analysis
The original dataset generated by the RFID readers had numerous duplicate scans generated from an individual sitting at the nest box entrance for an extended period of time. To remove these unnecessary detections, we removed consecutive scans within 15 seconds of one another, keeping the first one. We assigned each BBWD individual as “paired” or “unpaired,” with pairs defined as a male and female BBWD captured in the same nest box while the same nest was active. Since BBWD likely form long-term pair bonds [23, 24] and we observed previously confirmed pairs visiting the same nest boxes on RFID data, this includes BBWD pairs that were confirmed as such when captured in 2020–2022. We assigned each BBWD pair an identification number for analysis. Since only WODU females incubate the nest, we identified individuals by their PIT tag identification number.
We categorized PIT tag detections into prospecting, laying, and incubation periods for each BBWD pair and WODU female. We defined each period based on empirical field data: the prospecting period is defined as the time period when a pair/female does not have an active nest of their own (they may be laying parasitically or prospecting boxes for nest sites), the laying period as the interval between nest initiation and onset of incubation, and the incubation period as the time between onset of incubation and the nest hatching or being abandoned/depredated.
We conducted all analyses in R version 4.3.2 [25] using a Bayesian framework and Poisson underlying distributions. We used a generalized linear model to predict BBWD and WODU preference for old vs. new nest boxes during the prospecting period:
c i ~ Poisson(λi),
log(λi) = α + β1 * (nest box agei * speciesi)
where (ci) is the estimate number of visits to nest box i. We used the log link function to evaluate the relationship between our fixed effects and the expected count (λi). Our fixed effects consisted of an interaction between nest box age and the species. We ran this model twice: once to predict the number of unique new and old nest boxes visited by BBWD and WODU, and once to predict the total number of visits to new and old boxes by both species during the prospecting period. Additionally, we used a second model:
c i ~ Poisson(λi),
log(λi) = α + β1 * speciesi
to evaluate the relationship between species and the total number of unique nest boxes visited throughout the entire breeding period, regardless of nest box age. We ran this model an additional two times: once to predict the duration of the prospecting period, and once to predict the number of nest boxes visited during the prospecting period. We fitted models using package brms in R [26] (Bürkner 2017) and used the default uninformative prior distributions for fixed effects. We ran 4 Markov Chain Monte Carlo (MCMC) chains of 5,000 iterations, discarding 2,500 iterations during the warm-up period. We ensured that MCMC chains converged by examining trace plots and using the Gelman-Rubin statistic [27] (\(\widehat{\text{R}}\) value < 1.05). For all models, 95% credible intervals (CI) that did not overlap zero indicated significance. We confirmed model convergence by examining \(\widehat{\text{R}}\) values and trace plots.