We tested whether prior-summer and next-summer ozone and warmth product slopes differ for recorded bird richness. We used a published US bird survey panel and four region-period samples in 1995-2023. Models controlled route-observer pairs, region-year changes, climate and survey conditions. Only the early western contrast had a four-test corrected interval below zero: -0.946 species per 10 ppb per 2 C, with interval -1.564 to -0.329. The other three intervals included zero. Removing each state kept this contrast negative but could greatly widen uncertainty. A common time-specific pattern did not repeat. Future exposure cannot cause an earlier count, but shared causes and serial exposure correlation can affect both slopes. Monitor error, future-record selection, detection bias and few western clusters limit inference. This is not a causal or equivalence test. The bundle includes code, aggregate estimates, a figure, exact input hashes and a verified retrieval script. External record-level data remain in cited signed publications.
Filters: arxiv: q-bio.PE (Populations and Evolution) ×
We checked direct support for air pollution comparisons between beneficial invertebrates and significant pests in the public data from Ryalls et al. (2024). We selected 581 records from 89 publications. Only 32 records (5.5 percent) had both roles within the same pollutant, publication and broad outcome type. These formed four groups from two publications: a 1997 US study and a 2022 UK study. Ozone had two matched publications, nitrogen oxides had one, and sulfur dioxide and particulate matter had none. The UK paper supplied both pollutant comparisons. Descriptive ozone contrasts were +0.050 for US abundance, -0.136 for US survival and -0.657 for UK abundance, in beneficial minus pest log response ratios. An exact fixed effects calculation shows that only these two papers inform a role coefficient after separate intercepts for each publication and outcome pair are removed. No region and publication-era stratum had more than one matched paper. Publication years do not give field observation periods. Thus this release cannot support an independent matched test across regions and time. This audit does not test the validity of the original random effects model or show equal pollution effects. Broad outcome matches still compare different taxa and measurements. Climate interactions were not tested. Code, source hashes, aggregate evidence, a figure and independent numerical checks are included.
We tested whether changes in scored plant species change the annual ozone and soil moisture relation with recorded leaf injury. We used 4,271 site-years in two eastern US regions and two periods in 1994-2010. We compared an equal-species injury rate for all scored species with the same rate for species scored in every retained year at each site. The four paired product-slope changes were +0.97, +0.36, -2.78 and -0.13 injury percentage points per 10 ppm-h ozone per moisture SD. All four 95 percent intervals included zero. Fixed sets changed mean injury rates, but gave no clear repeated slope change. The outcomes differ from a prior binary injury study, so this test does not directly explain its slope. Selected persistent species, few clusters, variable plant counts, unknown missed injury and estimated exposures limit inference. One smaller sensitivity fit gave a negative variance estimate. Code, source hashes, aggregate results, independent numerical checks and a figure support reproduction from a verified source publication. Source records remain outside this package.
Question: did the apparent temperature sensitivity (days per °C) of spring phenology at Japan Meteorological Agency (JMA) stations change between 1953–1986 and 1987–2025, do heat-sum models with the same temperatures predict the change, and which model part must change to reproduce it? This is an independent replication of [citation 1], [citation 3] and [citation 2], with a new decomposition. Data: JMA records (regular species) for 8 spring events at 33–82 stations; daily station temperature from GHCN-Daily, with 6,415 gap station-months filled from JMA daily pages and each month anchored to the official JMA monthly mean. Method: per-station OLS slope of event day on the 60-day pre-event mean temperature, for observed days and for days from forcing-only (thermal-time) and two chill-forcing models fitted by grid search with one requirement per station. Uncertainty: year resampling shared by all stations and a permutation test over random year partitions. Results: for cherry first flowering, cherry full bloom, ginkgo budburst and wisteria flowering, the observed sensitivity became 0.73–1.04 days/°C stronger, while the forcing-only model predicted −0.08 to +0.14. Observed minus model: −0.86 (95% year CI −1.50 to −0.27), −0.72 (−1.32 to 0.02), −1.06 (−1.73 to −0.11) and −1.06 (−1.61 to −0.39); permutation p 0.002, 0.017, 0.013, 0.005 (Holm over 8 events: 0.016, 0.085, 0.078, 0.035). Mulberry, azalea, dandelion and ume showed no significant difference. The result held for 30–90-day windows, detrending, control for earlier temperature, official monthly temperature and a neighbour-station instrument. Its size depends on the split year, and moving-window curves are not monotonic. New: a period-specific heat requirement reproduced none of the change (−0.06 to +0.13 days/°C), so a requirement change (for example, from less chilling) cannot explain it in these models. A period-specific base temperature reproduced 36–120% of it; the shift was identified only for wisteria (−5 to +1 °C; year CI of shift 1 to 15 °C). Also: the forcing-only model explains weaker cherry sensitivity at warmer stations (observed gradient +0.41, model +0.36 days/°C per °C), but not the absent gradient for ginkgo, wisteria and mulberry; chill terms cut cherry cross-validated error by 12–17%. Limits: one urban reference plant per station, simple grid models, most events end in 2020, cause not identified.
Question: did the apparent temperature sensitivity (days per °C) of spring phenology at Japan Meteorological Agency (JMA) stations change between 1953–1986 and 1987–2025, do heat-sum models with the same temperatures predict the change, and which model part must change to reproduce it? This is an independent replication of [citation 1], [citation 3] and [citation 2], with a new decomposition. Data: JMA records (regular species) for 8 spring events at 33–82 stations; daily station temperature from GHCN-Daily, with 6,415 gap station-months filled from JMA daily pages and each month anchored to the official JMA monthly mean. Method: per-station OLS slope of event day on the 60-day pre-event mean temperature, for observed days and for days from forcing-only (thermal-time) and two chill-forcing models fitted by grid search with one requirement per station. Uncertainty: year resampling shared by all stations and a permutation test over random year partitions. Results: for cherry first flowering, cherry full bloom, ginkgo budburst and wisteria flowering, the observed sensitivity became 0.73–1.04 days/°C stronger, while the forcing-only model predicted −0.08 to +0.14. Observed minus model: −0.86 (95% year CI −1.50 to −0.27), −0.72 (−1.32 to 0.02), −1.06 (−1.73 to −0.11) and −1.06 (−1.61 to −0.39); permutation p 0.002, 0.017, 0.013, 0.005 (Holm over 8 events: 0.016, 0.085, 0.078, 0.035). Mulberry, azalea, dandelion and ume showed no significant difference. The result held for 30–90-day windows, detrending, control for earlier temperature, official monthly temperature and a neighbour-station instrument. Its size depends on the split year, and moving-window curves are not monotonic. New: a period-specific heat requirement reproduced none of the change (−0.06 to +0.13 days/°C), so a requirement change (for example, from less chilling) cannot explain it in these models. A period-specific base temperature reproduced 36–120% of it; the shift was identified only for wisteria (−5 to +1 °C; year CI of shift 1 to 15 °C). Also: the forcing-only model explains weaker cherry sensitivity at warmer stations (observed gradient +0.41, model +0.36 days/°C per °C), but not the absent gradient for ginkgo, wisteria and mulberry; chill terms cut cherry cross-validated error by 12–17%. Limits: one urban reference plant per station, simple grid models, most events end in 2020, cause not identified.
We tested whether curved previous-season ozone and temperature responses changed their interaction with recorded bird richness. We reused 30,986 surveys in four eastern and western US samples in 1995-2024. The four-term spline estimates ranged from -0.491 to +0.208 recorded species per 10 ppb ozone and 2 C. All four corrected intervals included zero. No paired coefficient change was clear. The early western corrected interval was below zero with six spline terms, but included zero with four. This threshold result depends on curve choice; it does not establish that curves explain the association. Selected routes, few western state clusters, coarse exposure and detection bias limit inference. This is a model check on reused observations, not a causal climate test. Code, the verified input panel, source hashes, aggregate results and numerical checks are included.
We tested whether low May-June rainfall strengthened the link between May-June maximum temperature and April to mid-June ozone at US bird survey routes. Four region and period models used 32,956 accepted surveys in 1995-2024. Models controlled route effects, Bird Conservation Region by year effects, and separate temperature and rainfall curves. The change in the ozone-temperature slope per 50 mm rainfall deficit ranged from -0.008 to +0.088 ppb per C. All four corrected intervals included zero. State by year controls gave a negative early eastern estimate, but this did not repeat and its two-way cluster interval included zero. Complete exposure records covered 20-34 percent of accepted western surveys and 64-66 percent of eastern surveys. The result tests exposure structure in a biodiversity sample. It does not test bird responses, daily extremes, or a causal climate effect. Different seasonal windows, changing monitors, coarse weather, sample selection and few western state clusters limit inference. Code, source hashes, aggregate results and a verified-parent retrieval script are included.
Artificial light at night can decrease flying insects, so it can decrease the food of birds that catch insects in the air. We tested this prediction with North American Breeding Bird Survey counts and a harmonized DMSP/VIIRS night-light series. We used four independent samples: the DMSP era (light 1992-94 to 2011-13, counts 1993-97 to 2010-14) and the VIIRS era (light 2014-16 to 2022-24, counts 2015-18 to 2022-25), each east and west of 100 W (827-1,560 routes per sample). A route-level long-difference model related the change in light to the change in counts of aerial insectivores (swallows, swifts, nightjars, tyrant flycatchers) minus non-insectivores, with BCR fixed effects and observer and survey-date controls. The contrast per doubling of (1 + DN) was +5.4% (95% interval -4.9 to +16.8), +0.3% (-9.7 to +11.3), +2.8% (-2.7 to +8.5) and -4.9% (-10.7 to +1.3). The signs did not agree and no test passed Holm correction. In the eastern samples, the intervals exclude relative decreases larger than about 5%. As a secondary result, other insectivores (mostly forest songbirds) decreased relative to non-insectivores where DMSP-era light increased in the east (-11.6%, -16.2 to -6.6; permutation p = 0.001). This result persisted in five checks but did not repeat in the VIIRS-era east (+0.4%, -2.6 to +3.5). In annual route-year panels, future light predicted counts as well as past light. Light change is a marker of development. Thus, these data cannot separate light from habitat loss, traffic noise or lower detection of songbirds near roads. The analysis does not measure deaths during migration at lit buildings. Code, derived route tables, source hashes and figures are included.
We tested two explanations for the shift from acid-bark lichens (acidophytes) to nitrogen-loving lichens (nitrophytes) in Great Britain: high ammonia-type nitrogen (NHx) deposition, or the decrease of sulphur (SOx). We used 77,729 British Lichen Society recording visits from GBIF (1960-2025) and the van Herk indicator lists. Poisson models counted indicator species per visit, with list length as the effort control, and compared changes from 1986-1995 across 2,582 10-km squares with CBED deposition. Bark nitrophytes increased to about 1.5 times the 1987-1989 level in England and in Scotland and Wales. SO2-tolerant acidophytes decreased to 0.30 and 0.45. The nitrophyte gain did not follow NHx: +0.035 log units per SD (95% CI -0.029 to 0.099) in 2005-2012, against a national gain of about 0.29. Acidophytes decreased less, not more, at high NHx. SO2-tolerant acidophytes decreased most where SOx decreased most (-0.158 per SD; repeated in both regions and in 2013-2025). SO2-sensitive acidophytes and nitrophytes had no clear SOx term; common non-indicator bark lichens had the largest SOx gain. The planned model gave a negative NHx term for nitrophytes, but this came from fewer churchyard surveys after 1995; models that compare visits of the same survey type remove it. Limitations: opportunistic records without substrate, 10-km exposure, correlated NHx, NOx and SOx, CBED ends in 2012, and indicator-only download after 1991 because of GBIF rate limits.
We tested whether an ozone by warming contrast in harebell flower width changed in the same direction from early to late summer at two ozone targets. We used four monthly visits in one North Wales experiment in 2021. The changes were -0.1224 log units for 80 versus 30 ppb and +0.0659 for 110 versus 30 ppb. The latter became negative when September was omitted. Each treatment had one dome. Treatment columns added no rank to chamber indicators. Fixed chamber offsets cancel in the seasonal contrast, but errors that change across months do not. Exact bounds include zero when allowed cell errors reach 1.54% and 0.83%, respectively. These are assumed error limits, not measured error rates or confidence limits. Missing flowers, selected width measurements and unequal achieved exposures limit inference. The results do not identify a causal warming effect or replicate across regions or years. Code, source hashes, aggregate results, checks and a figure are included. External source records remain outside the publication.
We tested whether prior seasonal warmth increased the prior ozone slope for unusually low recorded bird richness. We used a published US Breeding Bird Survey panel. Low richness was at or below the route-specific 20th percentile for routes with at least 10 eligible surveys. Four separate East and West, early and late samples used 28,404 surveys. Models controlled route-observer and region-year effects, survey conditions and weather. Interaction estimates ranged from 0.11 to 1.48 percentage points per 5 ppb ozone and 1 C. All four primary adjusted intervals included zero. Under this model, all upper bounds were below a 5-point reference size. A 10th-percentile check gave a positive early western signal that did not repeat. This result does not establish safety or a causal effect. Threshold choice, route selection, coarse exposure data, few western state clusters and detection bias limit inference. The package gives the exact published input, code, estimates, a figure and numerical checks.
We tested whether deep snow had a larger positive link to ground than tree bird reports on the same feeder checklists. We used 125,580 selected North American checklists from 2001-2005 and 2021-2024. Fixed region-specific species pools came from EltonTraits 1.0 and early reporting frequency. Models removed site-winter and state or province by week effects. They controlled reported effort and two-day minimum temperature. Ground-minus-tree contrasts were +0.63 and +1.03 percentage points in the East. Their four-test intervals were +0.02 to +1.24 and +0.54 to +1.51. West contrasts were +0.35 and +0.59; both four-test intervals included zero. Eastern secondary checks kept positive intervals. Deleting any one species kept all four point estimates positive. In later East, the difference included both more ground-pool reports and fewer tree-pool reports. This agrees with a regional food-access prediction but does not identify its cause. Strict trait pools, missing snow reports, selected sites, coarse weather, and group-specific detection limit inference. The study does not measure true occupancy, population gain, pollution effects, or a causal climate-change effect. Source records and transformed record-level data remain external. The package contains code, aggregate evidence, figures, exact source hashes, and independent checks.
We tested whether a current-year ozone and soil moisture product had a different leaf-injury slope from the next-year product. We used 3,717 consecutive-year records from a verified public-data panel in two eastern US regions and two periods from 1994 to 2010. Joint models removed site and year effects and controlled survey timing, plant and species counts, and temperature. Current-minus-future slopes were +0.096, +0.044, +0.102 and -0.051 injury probability per 10 ppm-h ozone per one moisture SD. All four two-way cluster 95 percent intervals included zero. Curved exposure terms and before-visit exposures also gave no clear contrast. A closer-monitor early Northeast result had a positive marginal interval but not a four-test interval. Thus this analysis did not establish a distinct current-year product slope. It does not establish equal slopes or disprove ozone damage. Selected repeat visits, few time clusters, coarse exposure, unknown missed injury and reuse of rounded source records limit inference. Code, aggregate evidence and verified-source retrieval instructions are included.
We tested whether separate responses to pond gains and losses improve predictions of regional duck shares. We used ten taxa and 24 fixed prairie regions from the public waterfowl breeding survey. Each model used previous bird shares, previous May pond levels, and current pond change. We fitted the models on 23 years before 2000 and tested them on 18 complete years in 2000–2019. Separate gain and loss slopes decreased mean squared Hellinger error by 0.013%. The nominal 95% test-year block range was −0.372% to +0.365%. A check with 300 training refits and test-year draws gave a range of −1.861% to +0.701%. Thus separate slopes gave no clear mean gain. In contrast, the symmetric pond-history model decreased error by 20.54% against bird history alone. Taxon, country, and time-split results had mixed signs. A smooth curved response also gave no clear mean gain. This is a retrospective model test on previously studied survey data. It does not establish equal ecological responses or a causal effect of ponds. The package preserves source data, code, provenance, predictions, uncertainty draws, figures, and reproduction instructions.
We tested whether partial sharing of pond coefficients improves predictions of duck distribution. We used ten duck taxa and 24 fixed prairie regions from the public waterfowl breeding survey. We chose the sharing strength with years before 2000. We then fitted models to 23 annual transitions before 2000 and tested 18 years in 2000–2019. The selected ridge penalty gave weak sharing. Mean squared Hellinger error was 0.025952, compared with 0.025936 for separate coefficients. The error decrease was −0.063%. A nominal 95% test-year block range was −0.102% to −0.024%. With training and test resamples, the range was −0.113% to +0.004%. That range held the selected penalty fixed and omitted survey error. All eight training-block deletion checks repeated selection and gave no mean test gain. One common coefficient gave 2.020% greater error than separate coefficients. The tested partial-sharing rule did not improve mean prediction. Its observed cost was small. The result applies to current survey shares for taxa with long training histories. It does not test scarce taxa or identify a causal habitat effect. The package preserves the source archive, code, provenance, all predictions, uncertainty draws, figures, and an independent arithmetic audit.
We tested whether short ice-season intervals occurred together more often at fixed Finnish lake sites. We used open NSIDC ice dates and selected 30 sites with at least 25 valid paired dates in each of 1961–1990 and 1991–2020. A short interval was below the site's observed early 10th percentile. Across 435 fixed pairs, mean joint frequency rose from 4.41% to 29.56%. The change was +25.15 percentage points, with a nominal 95% common five-year block interval of +7.88 to +43.77. An exact split gave +14.78 points for the product of local event frequencies and +10.37 for the excess in shared years. This excess depends on local event rates. When each site instead supplied its shortest 10% in each period, the joint measure changed by +0.46 points, with an interval of −2.02 to +3.78. Thus, more shared short seasons do not establish stronger alignment of locally rare events. Exact bounds for missing records retain an increase in full-calendar joint frequency at the fixed historical limits. The package preserves source files, code, all results, figures and independent checks. The selected observations do not establish a cause or a species response.
We tested the sensitivity of prior-season ozone slopes for recorded bird richness to an omitted factor. A published US bird survey panel gave negative slopes in four East and West, early and late samples after weather and survey controls. All four adjusted intervals included zero. A hypothetical factor with equal partial R2 values of 0.77-3.44% for ozone and richness can remove each slope. The pooled slope was -0.472 species per 10 ppb (95% interval -0.745 to -0.200). Its equal-strength threshold was 2.17%. Measured survey and weather blocks gave smaller calibration bias bounds than the pooled slope. These thresholds describe linear regression sensitivity. They do not show that an omitted factor exists or that ozone has no effect. Monitor coverage, measurement error, detection and selected surveys limit inference. The periods share some routes. Code, aggregate results, source hashes, a figure and a second numerical fit are included.
We tested whether prior seasonal ozone and warmth jointly predicted a rise in bird count concentration. We rebuilt counts from the USGS 2026 release and reused a verified exposure panel. Models used route-observer and state-year effects, survey controls and state-cluster errors. Four East and West, early and late samples in 1995-2024 gave mixed interaction signs. All four corrected intervals included zero. The pooled sample used 25,608 surveys on 1,993 routes. Its estimate was -0.118 concentration points per +10 ppb ozone and +5 C, with a 95 percent interval from -0.314 to +0.077. A closer-monitor check was negative. The tested data did not establish a repeated rise. Count concentration describes recorded birds, not true abundance. Coarse exposures, detection bias, selected routes, source-version changes and few western clusters limit inference. Code, exact derived evidence, source hashes, checks and a figure are included.
We tested whether precipitation before May adds information about May duck distribution. We used ten duck taxa and 24 prairie strata from the Waterfowl Breeding Population and Habitat Survey. The reference model used previous duck shares, previous pond levels and fixed spatial preferences. The added term used October–April precipitation from NASA POWER at one point per stratum. We fitted separate models on 16 target years in 1982–1999 and kept them fixed. On 18 test years in 2000–2019, added precipitation decreased mean squared Hellinger error by 1.57%. The nominal 95% test-year block interval was −5.67% to +8.46%. Five taxon point gains were positive and five were negative. A shorter training period gave a larger gain, so the result depends on the fit history. Adding current May pond change decreased error by 20.69% under the same primary comparison. This test did not establish a mean gain from the precipitation term. The weather data are a retrospective reanalysis. We did not test real-time data availability, other weather models or causal habitat effects. The package preserves data, source records, code, predictions, figures and repeat checks.
We tested whether breeding duck distributions were more concentrated in years with fewer May ponds. We used ten taxa, 24 fixed prairie survey regions, and 44 complete annual surveys from 1974 to 2019. The response was mean concentration from squared bird shares divided by region area shares. With a linear time trend, half as many ponds predicted a rise of 0.167 in concentration. The nominal 95% calendar-HAC interval was 0.065 to 0.268. A block interval was 0.083 to 0.251. An exact partition gave fitted changes of 0.096 in the shared distribution and 0.071 in differences between taxa. The largest positive taxon estimates were for green-winged teal, redhead, canvasback, and scaup. The US-only estimate had the opposite sign, and an entropy measure had an interval that included zero. An approximate survey-error check reduced the main estimate to 0.119. The intervals omit joint survey error. The result applies to estimated regional distributions under the tested model. It does not prove an effect of drought or measure use of individual wetlands. The package preserves source archives, provenance, code, results, figures, and repeat instructions. Separate calculations passed, and a full repeat gave identical hashes for 16 numerical files.