{'fraction_non_alphanumeric': 0.030792065212269706, 'fraction_numerical': 0.016279508080855674, 'mean_word_length': 4.938304420817012, 'pattern_counts': {'":': 0, '<': 8, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Nesting success of red-winged blackbirds (Agelaius phoeniceus) in marshes in an anthropogenic landscape

Recent analyses show significant population declines in many abundant avian species, especially marsh-nesting species including the red-winged blackbird (RWBL). Hypothesized causes include reduced nesting success resulting from changing land-use patterns and exposure to contaminants. Our goal was to test the hypothesis that landscape and nest characteristics as well as exposure to polychlorinated biphenyls (PCBs) correlate with nesting success. From 2008 to 2014, we measured clutch size, egg and nestling mass, hatching and fledging success and daily survival of 1293 RWBL nests from 32 marshes in the Hudson River valley of New York. Using generalized linear effect and survival models, we found that: (i) Julian date was negatively related to hatching success and clutch size but positively related to egg mass; (ii) nest height was negatively related to hatching success; (iii) nestling mass decreased with increased nest density and distance to edges; (iv) fledging success was significantly lower in nests closer to the ground that were far from water; and (v) clutch size and daily survival were higher in nests farther from water. Results showed that nesting success was correlated with variables associated with flooding, population density and predation and provided no support for the predicted negative effects of PCB exposure.

This study also was intended to test potential effects of PCBs on nesting success relative to other factors that have been shown to affect RWBL nesting success. These factors include within-season timing of nesting [63,64], proximity to wooded edges and buildings [21,65], but see [29,47,66], nest height [42] and nest density and distribution within marshes [42][43][44][45]67,68]. We predicted that nest predation rates would decrease with proximity to open water (a proxy for water depth) [65] and with nest density (increased group defence) [45] and would increase with proximity to human habitations or forest edges [29,47,69]. Increasing nest height should also provide protection against flooding [42]. We further predicted that measures of nesting success that control for effects of nest predation and flooding (e.g. clutch size, egg mass, hatching success, nestling mass and fledging success) would decrease with Julian date [70], nest density (density-dependent reductions in food supply) [42], and with increasing exposure to PCBs.

Study area
This study was conducted at 32 marsh (primarily cattail (Typha)-dominated) study sites during six field seasons, 2008-2014 (figure 1; electronic supplementary material, figures S2-S24 showing maps of all marshes). Each marsh covers a distinct contiguous area, of which many appear to be fragments of historical marsh areas. The investigated marshes ranged in size from 0.12 to 8.1 ha (table 1). PCB exposure from a source of contamination on the upper Hudson River has been reported and may have the potential to affect the nesting productivity of RWBL. In 2002, 40 red-winged blackbird eggs from nests found within the upper Hudson River floodplain were collected and analysed for PCBs [54]. Sampling results show that PCB concentrations in eggs ranged from 0.08 to 34.8 ppm (figure 2), with the lowest average concentration found south of the Mohawk River [54]. For this reason, marshes situated within four distinct regions were chosen to represent different levels of exposure to waters and sediments contaminated with PCBs (table 1). (i) Hudson River Region: six of the marshes were directly connected to the Hudson River where they were exposed throughout the year to sediments contaminated with PCBs from the Superfund site on the upper Hudson River. (ii) Mohawk River Region: eight additional marshes were selected from the nearby Mohawk River where nests were also exposed directly to river sediments, but not to PCBs from the Superfund site, which was on the upper Hudson River. A series of dams on the two rivers prevent any mixing of water from the two rivers in the study area. The nearest known PCB source in the Mohawk River is over 100 miles upstream from where it flows into the Hudson River [71]. PCBs were not detected in most water column and sediment samples collected from this portion of the Mohawk River suggesting that nest data collected from nearby marshes could represent reference conditions [71][72][73] . The Mohawk River Region was used as the reference area for the nests located directly on the border of the Hudson River because there were no significant marshes located on the Hudson River upstream of the source of PCBs. (iii) Hudson Floodplain Region: 13 marshes in the floodplain of the upper Hudson River were chosen because they were exposed to Hudson River sediment deposition during spring flooding events. These sites ranged from 6 to 41 km along the river downstream from the Superfund site. (iv) Hudson Valley Region: five additional marshes were selected in the Hudson valley that were above the 100-year floodplain elevation such that they were not exposed to PCBs from the Superfund site. These latter marshes were considered as controls for the Hudson Floodplain marshes. Within each of these regions, marshes were selected at random from those identified by satellite images and topographic maps showing floodplain limits, with the additional constraints that they had to be accessible and have landowners' permission to conduct research on their property.
Marshes were all located in landscapes dominated by intensive human land uses such as row-crops and pastures or hayfields, roads and yards surrounding buildings (table 1, electronic supplementary material, figures S1-S23). Some were bordered by floodplain forest or rows of trees and shrubs growing at the edges of the marshes, which formed distinct edges. Other marshes, especially those along the rivers also had distinct edges bordering open water that was too deep to sustain emergent aquatic vegetation such as cattails. Marshes bordering agricultural fields, pastures and hayfields did not form distinct, permanent edges because water levels varied greatly throughout the study causing marshes to expand and contract. For this reason, we restricted our analyses of distances to edges to woody borders and to distances to buildings as an indicator of residential effects on a marsh [29,47]. Nests were mostly in cattails (Typha spp.; 74%), although a few were in other plants such as common royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 buttonbush (Cephalanthus occidentalis; 7.7%), sedge (Carex spp.; 4.4%), invasive purple loosestrife (Lythrum salicaria; 3.9%) and other vegetation types.

Field methods
We used the methods described in Martin  primarily on accessibility, safety concerns and field staff. Marshes were searched along transects that were roughly 5 m apart. Two observers walked slowly along these transects searching for all active nests. This method was designed to find all nests that were active on the day of searching. Each marsh was searched twice-once early in the season in May and once late in the season in June. Nest searching was only conducted during days with no precipitation and temperatures greater than 10°C to minimize losses caused by interference from investigators. Because we were studying many marshes spread over a large area (figure 1), we were only able to visit nests every 3-4 days, which introduced some possible error in our estimates of hatching and fledging success [74]. When nests were found, we recorded the nest height above the ground or the water surface and the plant species in which it was located. Each nest was mapped using GPS (stabilized readings accuracy +/− 3 m). Nests were flagged 2-3 m from the nest location to avoid attracting the attention of predators [75]. Each time a nest was visited, we recorded its contents. We tried to weigh the eggs as early in the incubation period as possible given the 3-4 day interval between visits. If the nest appeared to have a complete clutch (3-4 eggs; see below), we weighed the eggs on the first visit and only weighed them again if an additional egg appeared. If the nest was found empty or with a smaller clutch (1-2 eggs), we waited until the next visit to weigh the eggs. Hatching success was  royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 estimated from the number of newly hatched nestlings and the presence of unhatched eggs in nests at least 1-2 days after hatching based on the feather development of the nestlings. We did not use nests found during the nestling stage to estimate hatching success. If eggs disappeared during incubation prior to hatching, we did not count these as unhatched eggs as they could have represented partial predation or brood reduction by the parent. Nestlings were weighed and measured (tarsus length and wing cord length) and their feather development described. We tried to measure the nestlings as late in the nestling period as possible but before they were sufficiently developed that they would fledge when we approached (6-9 days after hatching). Based on nests with known hatching dates, we knew that flight feather development (sheathed, partially sheathed, unsheathed) was a reliable indicator that nestlings were 6-9 days posthatching (unpublished data).
Fledging success estimates were based on the proportion of nests in which young survived long enough to fledge and there was no evidence of predation in the late nestling period. If nests were empty on the last check and they were far enough along in development to possibly have fledged, we searched the area around the nest for fledglings and for females carrying food [71]. If we did not find any evidence of fledglings, we then checked the nest for signs of likely fledging, such as a build-up of faeces in the nest [71]. If nests were damaged, we considered this as evidence of predation.
We used the standards proposed by the Ornithological Council [76,77] regarding the humane treatment of birds. The only animal contact in this study occurred when we weighed and measured nestlings, a procedure that resulted in no harm to the individuals and lasted less than 5 min per nestling. All work was done under permits from NYSDEC (LCP #1457) and US Bird Banding Laboratory (#23373).

Predictor variables
Predictor variables included specific nest characteristics, marsh habitat characteristics and estimated PCB exposure, all of which we had hypothesized to affect nesting success. Nest characteristics included nest height, date the nest was found and dates both egg and hatchling mass were collected. Marsh habitat characteristics included distance to woody edge, nest density and distance to open water and buildings. Distance measures were derived using spatial analyses. Collected GPS coordinates for each nest were compared with aerial photos and existing geographical information system (GIS) databases that included land use, hydrology, roads and wetlands. Distance to edge was calculated as the distance from each nest to the nearest woodland or human structure. Distance to open water was calculated as the distance from each nest to the closest unvegetated open water of ponds or rivers.
We used nesting density as a proxy for population density of RWBLs in marshes. We derived nest density for each marsh using a time-to-event nest abundance estimator that accounts for variability in detection and survival probability based on site (marsh) and nesting stage (incubation or nestling) [78,79]. Based on published studies, the average duration of the incubation period was assumed to be 11 days [80]. Models used to derive annual abundance estimates for each marsh were calculated for each survey year using the nestAbund package in program R [78]. Final models were chosen based on Akaike's information criteria (AIC) and varied by year. Annual abundance estimates for each marsh were used in conjunction with the estimated marsh areas to provide annual nest density estimates for each marsh. When a marsh was not able to be searched in its entirety in a given year, annual density estimates were adjusted to account for the reduction in area searched (electronic supplementary material, table S1).
PCB exposure categories were defined by the maximum sediment/soil total PCB concentration obtained from existing samples collected within the contiguous marsh areas located in the upper Hudson River floodplains. When no existing samples were present, the maximum concentration of the nearest existing soil sample was used. The upper Hudson River floodplains have been heavily sampled for PCBs. When a sampling location was not present within a given marsh, nearby samples were available within an average distance of 79 m (range: 0.6 to 279 m, standard deviation: 78 m). Marshes from the Mohawk River and Hudson Valley Regions were located more than 2300 m from PCB sampling locations and were considered reference locations. Given the uncertainty in associated PCB concentrations resulting from several marshes without collocated PCB samples, PCB data is treated as an index of exposure using an ordinal variable based on nest PCB categories of: (i) reference; (ii) 0-1 ppm PCB; (iii) 1-3 ppm PCB; and (iv) greater than 3 ppm PCB, corresponding to a range of potential PCB effects from no effect to high potential effect (figure 1). Measured PCB concentrations from eggs collected within the Hudson River floodplain [81] correspond somewhat royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 with the PCB categories derived for nearby marshes included in this study but show no obvious trend with distance from the PCB source (figure 2 and electronic supplementary material, figure S1).

Statistical methods
We used generalized linear mixed-effects models (function glmer in the R package lme4) [82] and generalized linear models (function glm in the R package stats) [83] to assess relationships between environmental and habitat variables and five reproductive metrics: (i) clutch size, (ii) egg mass, (iii) hatching success ( proportion of eggs that hatched, weighted by the total clutch size), (iv) nestling mass, and (v) fledging success ( proportion of nestlings that fledged, weighted by the total brood size). Each of these five non-independent response variables was analysed separately, so we adjusted the alpha value to 0.01. Clutch size was not highly variable (1190 of 1300 clutches had either three or four eggs), so we treated it as a bivariate factor with clutches of up to three eggs scored as a 0 and clutches of four or more eggs scored as a 1. We used a binomial distribution to analyse clutch size as well as hatching and fledging success for each individual egg and nestling, and a Gaussian distribution for egg and nestling mass. Egg mass was squared prior to analysis to better fit a Gaussian distribution. We assessed fit of models with Gaussian distributions by visual inspection of plots of model residuals versus fitted values, sample versus theoretical quantiles (q-q plot), and frequency distributions of model residuals, as well as plots of model residuals against each fixed effect. Similar evaluations are not effective for assessment of binomial models, but we balanced all predictors through transformations and scaling to improve model fit.
We started with similar global models for all response variables, which included variables that accounted for effects of seasonality or timing (Julian day of year for onset of incubation or number of days since the onset of incubation), water depth (distance to open water), flooding (distance to open water and nest height), predation (nest density, distance to open water and nest height), and PCB exposure (nest PCB category). Initial analyses of hatching success, fledging success and clutch size included year, treated as a categorical factor, as a random intercept to control for unmeasured amongyear variation in environmental conditions that might influence reproductive success. Year was retained in the hatching and fledging success models, where intra-class correlation (ICC) suggested that observations within Year could be correlated. Analyses of egg mass and nestling mass included nest ID as a random intercept, to account for non-independence of measurements of eggs and offspring from the same nest. All initial analyses included marsh region (Hudson Floodplain, Hudson River, Hudson Valley and Mohawk) to account for unmeasured regional variation that might influence reproductive success. Region was retained in the hatching and fledging success models based on ICC results suggesting correlation among results by region.
A series of 24 a priori models was developed for each response variable to account for potential interactions (table 2). We treated PCB level as ordinal data for analyses, which facilitated consideration of the effects of PCBs that might depend on other environmental factors (i.e. inclusion in interaction effects). Distance to edge and distance to open water were right-skewed and so were log-transformed prior to analyses. Nest density had a slightly right-skewed distribution that was better improved through square-root transformation. We scaled all continuous fixed effects (after any required transformations) to a mean of zero and standard deviation of one, which allows for more direct comparison of effect sizes among fixed effects. Results were plotted with unscaled data, for ease of interpretation. Multi-collinearity among continuous fixed effects was not an issue (all pairwise |r| < 0.31).
We included scaled calendar day of the onset of incubation as a covariate in all global models. RWBLs may produce multiple broods per year [80], but the expectation is that reproductive investment might decline across the breeding season (table 2). We included nest height as a covariate in analyses of hatching and fledging success because nest height might alter exposure to some sources of nest failure (e.g. flooding or predation). There was little variability in nest height for nests within 5 m of the surface and the majority of nests known to be flooded were built less than or equal to 0.5 m of the surface. Therefore, to best represent risk from flooding, we classified nest height as a categorical variable scored low for nests less than or equal to 0.5 m and high for nests greater than 0.5 m from the surface. Eggs were weighed at differing days across the approximately 14-day incubation period, and so we included the scaled day of incubation as a covariate to control for water loss, which reduces egg mass over the course of incubation. Lastly, because nestlings were weighed at different ages (range 5-12 days old), we included a scaled estimate of nestling age as a covariate in analyses of royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 nestling mass. All covariates were included only as main effects and not considered in any interaction terms. Initial analyses of hatching success, fledging success and clutch size included year, treated as a categorical factor, as a random intercept to control for unmeasured among-year variation in environmental conditions that might influence reproductive success. Inclusion of year improved model fit for all global models, as indicated by a lower Akaike's information criterion corrected for Table 2. List of fixed variables included in 24 a priori models used to evaluate nesting success response data collected from marshes in 2008-2014 in various regions along the upper Hudson River and Mohawk River. Random variables were also included in each model to offset unmeasured variation among-years and by region and to account for non-independence of eggs and offspring from the same nests. , distance to open water (dWater), nest height categorized as less than or equal to 0.5 m or greater than 0.5 m (NestHt), nest density (Density), distance to edge (dEdge), and PCB. For the egg mass and nestling mass models, additional fixed variables were included in all models, to account for age of the egg or nestling (incubation day and nestling age, respectively).
royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 small sample size (AICc) for global models with versus without (using a generalized linear model with a binomial distribution) this random effect. Analyses of egg mass and nestling mass included nest ID as a random intercept, to account for non-independence of measurements of eggs and offspring from the same nest.
To calculate model-averaged parameter estimates (function model.ave in the R package MuMln) [84], we included model results for each nest success response variable that were ranked within 2 AICc of the topranked model (i.e. lowest AICc). Model-averaged results did not differ qualitatively from analyses that used the top-ranked model, although effects with weak support in the top-ranked model (i.e. 0.05 < p < 0.10) received less support in model-averaged results. Given the relatively large sample size for all analyses, these marginal effects are unlikely to be biologically important. We used the top-ranked model to generate plots to illustrate any effects identified as important following model averaging.
We used three different datasets for analyses of hatching and fledging success: one that included all nests, one that excluded nests that failed due to predation during incubation or the nestling stage (depending on the response variable being analysed), and a final one that excluded depredated nests as well as nests that failed due to flooding during the relevant breeding stage. Nests that failed during incubation were not included in analyses of fledging success. Predictor variables were adjusted to exclude those that corresponded directly with the excluded nests (table 2). For example, the fixed variables that corresponded with predation were not included when predated nests were removed from the dataset. Comparisons of results across these three complementary analyses provide insight into the influence of flooding and predation on associations between reproductive success and habitat characteristics (i.e. the fixed effects). We analysed the remaining response variables (clutch size, egg mass, nestling mass) using only one dataset, which included all data regardless of eventual nest outcome.
We evaluated daily nest survival as a function of the predictor variables using the RMark package [85]. RMark acts as an interface for the MARK program [86]. The same set of 24 a priori models evaluated for the generalized linear and mixed effect models were included in the analysis (table 2). Random variables were not included in the daily nest survival models; however, Year was included as a categorical variable for all models to account for annual variability due to weather conditions and other unmeasured among-year variability in environmental conditions. Nest age, the number of days since the onset of incubation, was also included as a continuous variable in all models to account for known variability in daily survival over the nesting period [87,88]. As with the previous models, distance to edge and distance to open water were log-transformed and nest density was square-root transformed prior to analysis.
Competing daily survival models were ranked based on AICc. Modelled coefficients of predictor variables included in models ranked within 2 AICc of the topped ranked model were considered to have strong support when their confidence intervals did not include zero [89]. Potential effects of highly supported predictor variables were evaluated based on the highest ranked model in which they were included using mean values of continuous variables and categorical variables that best represent mean daily survival from egg incubation through fledging.

Results
Reproductive metrics and nest characteristics used as fixed effects are summarized by region in tables 3 and 4, respectively. All relationships between reproductive metrics and fixed effects are reported with an effect estimate (β) ± standard error derived from full model averaging across all models ranked within 2 AICc of the top-ranked model (table 5). Effect estimates can be compared among continuous fixed effects within a model because they were scaled (mean = 0, s.d. = 1) prior to analysis.

Hatching success
In analysis of the full dataset (N = 1293 nests), which retained data from all nests, including those that failed due to predation or flooding during the incubation stage, hatching success declined as the breeding season progressed (β = −0.16 ± 0.03, z = 4.73, p < 0.001, table 6 and figure 5a). Hatching success was also lower in nests that were below 0.5 m from the ground (β = −0.24 ± 0.07, z = 3.50, p < 0.001, table 6 and figure 5b). There is a trend for the number of hatchlings to increase with distance to water (β = 0.10 ± 0.05, z = 1.99, p = 0.047, table 6). No other fixed effects or interactions retained in models ranked within 2 AICc of the top-ranked model predicted hatching success (table 5).
When analyses were repeated with nests that failed due to predation during the incubation stage excluded from analyses (N = 1019 nests retained, 274 depredated nests excluded), results were similar to those from the full dataset. Hatching success was lower in low nests relative to higher nests (β = −0.45 ± 0.09, z = 4.71, p < 0.001, table 6 and figure 6).
When the analysis was repeated with both nests that were depredated as well as nests that failed due to flooding during the incubation stage excluded (N = 971, 48 flooded and 275 depredated nests excluded), results differed somewhat from those of the analysis that only excluded depredated nests. Hatching success no longer differed between low and high nests, as is expected if lower nests are more likely to fail due to flooding (table 6).

Nestling mass
As expected, nestling mass increased with nestling age (N = 1446 nestlings from 512 nests, β = 1.90 ± 0.17, z = 10.96, p < 0.001, table 6 and figure 7a). Nestling mass also increased with PCB exposure (β = 0.50 ± 0.19, z = 2.70, p = 0.007, table 6 and figure 7b), but decreased with distance to edge (wooded or human   royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266   Table 6. Fixed effect estimates (with standard error) resulting from full model averaging across all best-performing models (ΔAICc < 2) to explain variation in nest success metrics in red-winged blackbirds. Random effect intercepts were used to account for non-independence of groups. Random effects are included in models where the intra-class correlation coefficient (ICC) suggests that observations within the random variables may be correlated.  Clutch size was treated as a binomial variable for analyses, with clutches less than or equal to 3 coded as 0 and greater than or equal to 4 coded as 1. c Day of year is the calendar day of the onset of incubation (with 1 May = 121). d Distance to water is the shortest distance to open water in wetlands or the river. e Incubation day is the estimated day of the approximately 14-day incubation period when eggs were weighed, included as a covariate to control for water loss. f Nest height is categorical, with nests described as low (less than or equal to 0.6 m) or high (greater than 0.6 m). g Distance to edge is the shortest distance to woods or human structures. h Nestling age is the estimated age of the nestlings when they were weighed, included as a covariate to control for growth over nestling development.
royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 15     royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 effects or interactions retained in models ranked within 2 AICc of the top-ranked model predicted fledging success (table 5). When analyses were restricted to exclude nests that failed during the nestling phase due to predation (N = 778 nests retained, 144 depredated nests excluded) or failed due to predation and flooding (N = 767   Figure 6. Model-predicted relationships between hatching success in red-winged blackbird nests ( proportion of eggs that hatched) and nest characteristics when known predated nests were removed from the analysis. When compared with analyses that included all nests, hatching success was still lower in nests that were built within 0.5 m of the ground surface. Nests higher than 0.5 m displayed greater hatching success. Model predictions are derived by specifying varying levels of the fixed effects illustrated in the figure, while holding all other fixed effects in the top-ranked model constant. Continuous variables were set at their mean value, and the categorical factor (nest height) was specified as the proportion of nests that were low (i.e. nests less than or equal to 0.5 m off the ground). See main text for statistical support.

Daily nest survival
Daily nest survival was best predicted by the model that included the variables year, nest age, day of year and distance to open water (table 7). The average rate of daily nest survival for all years is 0.96. This survival rate was best represented by daily nest survival in 2014, which was used as the base year to present model results for all other significant parameters (figure 9). The confidence intervals of both nest age and distance to water did not include zero, indicating strong support for their potential effects on daily nest survival (table 8). As expected, daily nest survival decreased with nest age (β = −0.05 ± 0.01, table 8). Conversely, daily nest survival increased with increased distance to water (β = 0.17 ± 0.08, table 8). Three additional models were within 2 AICc of the top-ranked model predicting daily nest survival (table 7). The second and third ranked models were variants of the top-ranked model. The fourth ranked model also included PCB and the interaction of PCB and distance to water; however, confidence intervals of both parameters included zero suggesting a lack of support for their effects. (c) ( d) Figure 7. Mass of red-winged blackbird nestlings increased with nestling age (a) and local PCB level (b) and decreased with distance to edges (c) and nest density (d ). The lines illustrate model-predicted relationships, and shaded areas illustrate 95% confidence interval. Model predictions are derived by specifying levels of the predictor variables in each figure (x-axis spans range of the data), while holding all other continuous fixed effects in the top-ranked model constant at their mean value. The categorical factor (nest height) was specified as the proportion of nests that were low (i.e. nests less than or equal to 0.5 m off the ground). See main text for statistical support.

Discussion
We found little evidence that characteristics associated with the human-dominated landscapes in the region were negatively correlated with nesting productivity of RWBLs. There were no edge-related reductions of nesting success, a reliable index of fitness [90], even for the few nests located immediately adjacent to wooded edges where predator communities may be different [69]. By contrast, we found that nestling mass was higher when nests were closer to wooded edges or buildings. As expected, nestling mass was also greater in areas with lower nest density, an effect that has been observed in other species nesting in human-modified landscapes [91]. Additionally, we found no evidence that exposure to PCBs was negatively correlated with nesting success of RWBLs in the region. In fact, PCB exposure was positively related to nestling mass; however, it is likely that the increase in nestling mass was not a direct result of PCB exposure, but instead the result of an  unmeasured variable correlated with areas of higher PCB exposure. Nest height was related to hatching success because lower nests were more likely to be flooded during the incubation phase. This was evident with the continued negative effect of nest height in models that analysed hatching success in nests not affected by predation, which might also be affected by nest height. Nesting success declined with date within the season, a reflection of decreasing clutch size and hatching success, but egg mass increased.

Edge effects
The relative lack of strong effects of distances to open water and wooded edges (or human structures) on nesting productivity in this study is generally consistent with the other studies that have examined edge effects in RWBLs [21,65]. Only nestling mass was found to have a significant relationship with distance to wooded edges. These results are consistent with Whittingham & Robertson [92] who found that parental provisioning was higher when parents foraged off-territory in woodland habitat. The increase in clutch size and daily survival rate with increased distance to water was not consistent with predictions based on protection provided by deeper, open water. Picman et al. [65] found that water depth, not distance to royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 edge, had a significant effect on predation rates in marshes, which they attributed to deeper water providing protection against mammalian nest predators. A similar mechanism may be responsible for the observed decrease in fledging success in nests built close to the ground when associated with increased distance to water. We did not, however, measure water depth under nests during our study and therefore do not have strong tests of this hypothesis. Our results also differ from those of Vierling [29], who found that nesting sites close to human habitation had significantly lower nesting success than those in other landscape contexts. A review of the many studies of RWBL nesting success [42] found no consistent effects of open water and water depth on nesting success of RWBLs. Higher predation rates in nests adjacent to forests might reflect vulnerability to forest-based predators such as some snakes and raccoons (Procyon lotor), which rarely visit deep water [93], but we found no consistent negative effects of distance to woody vegetation or buildings. Our results add to the growing number of studies that show few obvious edge effects on nesting success of birds in wetlands and grassy habitats [11].

Exposure to PCBs
Nesting success of RWBLs was not reduced in marshes with relatively high PCB concentrations compared with marshes with low PCB concentrations and with presumably limited or no PCB exposure. Nestling mass increased significantly with likely exposure to PCBs; given the lack of known mechanisms by which PCBs would increase nestling mass, it is possible that this reflects differences in the quality of food available in marshes. Although results were based on an index of potential PCB exposure rather than measured PCB concentrations in eggs or nestlings, results are consistent with other studies that have shown no effects of PCB exposure on nesting success [6,33,51,52,56,94]. RWBL foraging patterns are highly variable and increase the uncertainty of PCB exposure to marsh-nesting RWBLs that frequently forage in upland locations [95]. Adult RWBL are known to feed their young predatory aquatic odonates [35,[37][38][39] that are assumed to bioaccumulate toxins; feeding patterns, however, are highly variable. Although RWBLs are considered to be a species of intermediate sensitivity to PCBs [60], at least some eggs contain high concentrations of PCBs in the region (electronic supplementary material, figure S1; [81]). Additional studies that examine nesting success in conjunction with egg analysis would provide additional support for the lack of potential effects of PCBs on RWBL nest success [35,[37][38][39]. The lack of effects of likely PCB exposure on nesting success does not preclude possible adverse effects associated with other environmental contaminants, but it does show the continued need for field studies to understand the relative effects of toxins compared with other factors [96] including anthropogenic factors such as vehicle or structure collisions, wetland alteration, predation, or other myriad human effects [97]. Table 8. Effect estimates (with standard error) resulting from the best-performing model to explain variation in daily nest survival rates in red-winged blackbirds. Years were included as a dummy variable to assess effects between years. All analyses were conducted using a logistic exposure model in RMark. Parameters estimates with confidence intervals that do not include zero are indicated in bold.

. Nest height
Our data suggest that the primary risk associated with nesting low over the water is associated with flooding, which was primarily observed where water-level management occasionally led to spikes in water levels during the nesting season (e.g. Lions Park on the Mohawk River). Overall, only 4% of the nests were flooded. In fact, in nests closer to the ground, fledging success was higher in areas closer to water. Any increases in flooding frequency during the nesting season brought about by changing climates (e.g. higher frequency of severe summer rainstorms) [98,99] or by altered hydrological cycles caused by changes in drainage or water-level management in impoundments could increase losses to flooding during the breeding season.

Within-season variation
A decline in clutch size, hatching success and daily survival rate as the season progresses has been widely documented in birds and other vertebrates [70,100] and has been documented specifically in RWBLs [35,64]. Interestingly, the smaller late-season clutches had eggs that were, on average, larger than early season eggs. This pattern has been widely documented in vertebrates [63] and might represent a way in which RWBLs compensate for smaller clutches by laying eggs with more supplies to combat declining food supplies as marshes drain or to give their fledglings a greater chance of surviving in the face of competition with fledglings that left the nest earlier in the breeding season. Alternatively, it might be a response to increasing day length during the breeding season, providing more time for follicular growth during egg development [101].

Nesting density
Changing population densities of RWBLs should have contrasting effects on nest predation rates and on the growth and development of young. Nesting success should increase if mobbing effectively deters nest predators [45], but productivity from nests that escape predation should be lower if there are negative density-dependent reductions in food within marshes or in the nearby upland habitats [42]. If RWBLs are forced to nest at increasingly high population densities as a result of decreasing habitat availability, then increasing human land use could lead to lower nest predation rates, but reduced productivity from nests that survive. Alternatively, decreased regional populations may lead to decreased nesting density within marshes potentially leading to increased body condition of nestlings within surviving nests, but also increased nest predation rates. Our results support these predicted effects as we documented declining nestling mass in areas of higher density. Our results suggest that, although not significant at α = 0.01, increased nest density may have increased productivity through a reduction in potential predation. Avian nest predators can be deterred by mobbing behaviours [42]. In RWBLs, this is especially evident when defending against marsh wrens (Cistothorus palustris), which are prominent RWBL nest predators [43] and were present in most marshes (unpublished data). Although probably not as effective, mobbing behaviour has also been observed for mammalian nest predators such as mink [102] and raccoons [103], both of which are common in and around the Hudson River and its tributaries.

Possible additional factors
Other factors that have been listed as possible causes of reproductive failure that we have not studied here include other toxins such as methyl mercury and pesticides; reduced insect populations; proliferation of invasive species; decreasing marsh size; loss of upland nesting habitat; climate change; and direct mortality to pest control efforts [1]. The other factors listed can be related to human disturbance. In addition, suburban and rural landscapes have been shown to act as population sinks for RWBL [29] suggesting that measured nesting characteristics may have variable effects with differing levels of human disturbance. It is notable that because all marshes in the study area were located in rural locations with moderate levels of human influence, there is no measured range of human influence or reference areas to directly assess potential effects of increased or decreased human disturbance. Our marshes are relatively free of invasive species and there were no obvious differences in nesting success in the few nests we found in invasive plants (unpublished data). We had no data royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 9: 220266 from upland habitats because our study was designed originally to focus on nests that would be exposed to PCBs in the floodplain. For these variables, our analysis offers no insights.

Conclusion
Several variables including nest height and distance to open water or edges, both alone and in combination, had statistically significant effects on RWBL nesting success. However, few of the effects we documented were associated with large decreases in nesting success. Also, PCBs were not found to be negatively correlated with nesting success in the upper Hudson River. These results suggest that the long-term population declines of RWBLs in this region [2] may be caused by other factors such as outright loss of marsh habitat [22], climate change including an increase in summer flooding [99,[104][105][106], and intensification of agriculture, which can include loss of landscape heterogeneity [12] and earlier harvesting of hayfields where they also nest [13,28,36,107]. Alternatively, these declines could be caused by factors outside the breeding season, including killing of birds at roosting or feeding aggregation sites.
Ethics. We used the standards proposed by the Ornithological Council [76,77] regarding the humane treatment of birds.
All work was done under permits from NYSDEC (LCP #1457) and US Bird Banding Laboratory (#23373).
The data are provided in electronic supplementary material [110].