Study sites and study species
The study islands (locally known as chars) were located in the Ganges-Padma River floodplain of Rajshahi (Char Majardia) and in the Brahmaputra-Jamuna River floodplain of Kurigram (Char Sakahati), Bangladesh referred to hereafter as Rajshahi and Kurigram (Fig. 1). The total area of the Rajshahi site covered 19.4 sq km, and was approximately 1.2 km from the mainland, while the Kurigram site covered 35.2 sq km and was about 3.4 km from the mainland, based on 2018 Google Earth data. The two sample islands are about 170 km apart. The study areas were chosen to represent conditions along the two rivers and particularly for the presence of Bristled Grassbird (Jahan et al. 2022). These sites are roughly near the eastern edge of the Bristled Grassbird range and roughly in the center of the mainland Asian range for Striated Grassbird (BirdLife International 2022). Both areas are subject to seasonal flooding during the monsoon (July–October). During the bird breeding season, when surveys were conducted (May to July), wet, inundated areas were scarce. Land-cover of both areas was characterized by floodplain grasslands containing a mosaic of patches of both short grasses (Cynodon sp.) and tall grasses (Saccharum sp.), with a mix of forbs and bushes. Forbs and bushes included Solanum sp., Tamarisk sp., Ipomoea carnea, Lippia nodiflora, Persicaria sp. and Xanthium strumarium, and were usually mixed with grasses. For the two woody species:
Ipomoea carnea, considered invasive to these tropical grasslands (Akter and Zuberi, 2009), was more abundant (based on individual plants/sample point) in Kurigram (0.58 stems per sample point) than at Rajshahi (0.07stems/point); while Tamarisk sp. was more abundant in Rajshahi (3.88 stems/point) than at Kurigram (0.75 stems/point). Ipomoea carnea can reach heights of up to 5 m and occurs in patches or clumps amid other vegetation. Tamarisk sp. also can reach heights of up to 5 m and occurs in extensive patches or small patches amid other vegetation. Saccharum sp. mostly grows up to 3 m tall and occur in dense and more extensive patches. Croplands and small villages were also characteristic of the landscape. Differences in crop varieties were noted between islands. In Kurigram site a wide range of crops were planted including rice, lentils, jute, maize, sesame, dhaincha (Sesbania bispinosa), berries (Ziziphus mauritiana), pointed gourd, bitter gourd, chili, and other vegetables; while in Rajshahi site mostly rice was planted with small areas of sesame and jute. Kurigram site supported year-round human occupancy, in villages located on the island, while Rajshahi site had no human occupants within the study area, although a village named Nobinagar was located at the south edge of the island, adjacent to the study area, close to the border with India (part of the island is located in India, Fig. 1). Livestock grazing (cows Bos indicus, goats Capra hircus, sheep Ovis aries, and water buffaloes Bubalus bubalis) and grass (Saccharum sp.) harvesting by local people takes place at both sites before the onset of the monsoon. Grazing intensity was observed to be higher in Rajshahi compared to Kurigram in 2019 (Fig. 1).
The Striated Grassbird Megalurus palustris toklao is typically present in grassy floodplains throughout the country and is considered to be a common and widespread resident (IUCN Bangladesh 2015). Bristled Grassbird Schoenicola striatus, formerly described as common in parts of Bangladesh in the nineteenth century (Ali and Ripley 1987), went un-recorded in the country for over a century. It was rediscovered in northeast Bangladesh in 2014, and considered to be rare (Haque and Tareq 2014, Khan et al. 2015). However, a recent study found that Bristled Grassbird was more widespread along the major rivers than previously recognized (Chowdhury et al. 2022). Both species breed during the summer and early monsoon season (March-July) (Jahan et al. 2022, Khan et al. 2015). While Striated Grassbird is a highly visible, year-round resident (IUCN Bangladesh, 2015), outside the breeding season Bristled Grassbird remains largely unrecorded probably due to a relative lack of surveys.
Survey methods
Abundance data for both grassbird species were collected during their breeding season, when birds were singing and most easily detected. Point-count surveys were conducted from 21 May to 2 July in 2018 and from 22 May to 16 July in 2019. A total of 382 point-counts were conducted (185 points in 2018 and 197 in 2019) by overlaying a grid of 250 x 250 m on 2018 Google Earth imagery of the study islands. The grid size was decided based on field surveys to minimize double counting the same birds in adjacent points. Based on image interpretation and field checking, grid points that fell on water, bare sand, woodland, forest plantations and human settlements were excluded, and the remaining points were in grassland patches and other open-country vegetation (bushes, reed, swamp, scrub, cultivation). Sample points were chosen randomly from among these, one per grid-square so that all sample points were ≥ 250 m apart. One observer (IJ) recorded individual Bristled Grassbirds and Striated Grassbirds seen or heard within 100-m fixed-radius point-counts (Hutto et al. 1986) of 10-min duration without a waiting period (Savard and Hooper 1995). Both species were noted as either heard, or seen in song-flight, or perching. Differences in song flight behaviour of both species were also recorded. Distances between the observer and birds were estimated directly using a rangefinder when birds were visible, otherwise in 5 m bands to 30 m, then 10 m bands to 100 m (Rosenstock et al. 2002). With points spaced 250 m apart, and birds only counted to a distance of 100 m at each point, the 50 m buffer minimized the likelihood of double-counting the same individual. Moreover, the observer also carefully noted the locations of individual birds seen at adjacent points to further avoid double-counting the same individuals. Surveys took place between 06h00 and 10h00 under good weather conditions (i.e., with no rain or strong winds). Weather conditions (wind and cloud cover) were also recorded. Wind velocity was estimated following the Beaufort scale, cloud cover (visual percentage approximation) and survey time and duration were recorded. Rajshahi was surveyed at 94 points which were sampled in the two survey years and an additional 12 points surveyed in cropland in 2019 (Fig. 1). Kurigram was surveyed at 91 points which were also sampled twice (once in each survey year). In total, the detection area of the survey points covered about 17% of the Rajshahi site and 8% of the Kurigram site.
Habitat measurements
In 2018 and 2019, plant species and height were recorded at 5 m intervals along a 100 m radius from each bird survey point in each of the four cardinal directions (see Jahan et al. 2022 for additional details). This gave a total of 80 vegetation samples per bird survey point. The plants were classified into four categories: Saccharum sp., Cynodon sp., forbs and bushes, or cropland. Bristled Grassbird and Striated Grassbird were observed using Saccharum sp., Ipomoea carnea, Tamarisk sp. forbs (> 3m height) and trees (> 7m heigt) as their perching substrate (Fig. 2). Of which Ipomoea carnea and Tamarisk sp., were tall and woody classed within the forbs and bushes category, were additionally classified as song posts. Crop and non-crop vegetation height (cm) was recorded as an average for each sample point within a 50 cm radius at 5 m intervals. Human use variables (grazing and harvesting) were measured only in 2019. Grazing intensity was based on a count of livestock dung (Sheidai-Karkaj et al. 2022) and grass harvesting intensity was based on a count of cut grass stems within a 1.5 m radius, every 5 m interval along the four cardinal direction sampling transects. Some changes in habitat were observed between 2018 and 2019: four points previously covered with Saccharum sp. or forbs and bushes in 2018 were bare land in 2019, and an additional five points previously covered with Saccharum sp. or forbs and bushes were inundated in 2019. Crop cultivation increased in 2019 (29 points) compared to 2018 (12 points). This comparison does not include an additional 12 points which were surveyed in cropland in 2019 in order to increase the sample of cultivated points. Vegetation heterogeneity was quantified based on the extent of different vegetation heights recorded at each bird survey point. Vegetation heterogeneity was summarized by calculating means and standard deviations of vegetation height measures within each bird survey plot (n = 80) and then averaged these across all plots (n = 382). Vegetation heterogeneity was considered as the extent to which short vegetation and tussocks of taller vegetation occur in proximity. Areas with a mix of both short and tall vegetation were considered to be more heterogeneous than areas with only tall or only short vegetation for example.
Analysis
Vegetation associations
Both male and female Bristled Grassbirds and Striated Grassbirds were recorded at each sample point. However, females of both grassbird species were relatively cryptic, whereas males were more readily detected, particularly in song flight, and when perched and singing. All the aural detections from song were also likely from males. Therefore, models assessing vegetation associations and density estimates were performed for only males of both species. To evaluate the occurrence of singing males in 100-m point counts of each species relative to habitat variables, we used binomial regression models in a generalized linear mixed models (GLMMs) framework, using R package glmmTMB (Brooks et al. 2017). Models were fitted as a binary response variable: detected = 1 and not detected = 0. Model sets were developed to conduct vegetation structural preference analysis and the effects of grazing and harvesting for Bristled Grassbird and Striated Grassbird combining both study islands and both years. We used the same set of randomly sampled points to replicate the survey in both years (apart from the 12 additional points in cropland in 2019). However, the impact of the differing extent of flooding and changes in vegetation on bird usage during the two years meant that these sets could reasonably be considered to be independent of one another. Two sets of models were developed: the first set of models included vegetation type with vegetation heterogeneity, height and song post, using data recorded in 2018 and 2019. The second set of models added grazing and harvesting intensity, these data were recorded in 2019 only. Vegetation categories were Saccharum sp., Cynodon sp., forbs and bushes, and cultivated crops. Use of tall, woody song posts (Ipomoea carnea and Tamarisk sp.) was analysed for both of the bird species. Both model sets contained each variable and combination of variables, in addition to null and global models. Collinearity among predictor variables was checked before analysis. Variables with correlations above 0.70 were not included in the same model per Dormann et al. (2013). AIC model comparison methods (Burnham and Anderson 2002) were used to evaluate the relative ability of each model to explain variation in occurrence for Bristled Grassbird and Striated Grassbird. All models were checked for the homogeneity of their residuals. Models were compared and ranked using Akaike’s information criterion corrected for small sample size (AICc) following Burnham and Anderson (2002). Models and their relative importance were also assessed using Arnold (2010) to avoid the inclusion of uninformative parameters. AICc weights (w) provided the relative weight of any particular model (Burnham and Anderson 2002). Model coefficients (β) with 95% confidence intervals (CI) from AICc selected models (Burnham and Anderson 2002, Arnold 2010) were reported, and predictor variables were considered to be strongly influential on response variables if the CI for β did not overlap zero (Bolker 2008, Arnold 2010). We also used model-averaged coefficients to make multi-model inferences when model-selection uncertainty was high (Burnham and Anderson 2002). All statistical analyses were carried out in the Program R version 4.1.0 (R Core Team, 2020) using glmmTMB, tidyverse, ggpubr, AICcmodavg, MuMIn, visreg packages.
Other habitat variables tested included the presence of standing water within 100 m of the survey point as an indicator of ‘wet’ habitat as a rough test of a hypothesis regarding the relative preferences for the two grassbird species for wetter versus drier habitats.
Grassbird density
Distance sampling was applied to estimate bird densities using Program DISTANCE (Buckland et al. 2001). Four estimators: uniform, half normal, hazard-rate, and negative exponential were employed to determine the detection model of best fit. Hazard-rate models were not considered when the output gave implausible shapes (Thomas et al. 2010). All detection models were run with the cosine adjustment. Models were evaluated using AICc (Buckland et al. 2001). The detection model with the lowest AIC score was considered the best fit and was used to estimate density of individuals within that specific habitat category (McCollum et al. 2018). Kolmogorov-Smirnov goodness-of-fit (GOF) test, were used to evaluated models with models having P > 0.05 considered well-fitted to the data (Buckland et al. 2001).