Contrasting changes in reproductive success of birds in cities with reduced and increased human disturbance during the COVID-19 lockdown

Gábor Seress (  seressg@almos.uni-pannon.hu ) MTA-PE Evolutionary Ecology Research Group, University of Pannonia, Veszprém Krisztina Sándor Behavioral Ecology Research Group, Center for Natural Sciences, University of Pannonia, Veszprém Ernő Vincze MTA-PE Evolutionary Ecology Research Group, University of Pannonia, Veszprém Ivett Pipoly MTA-PE Evolutionary Ecology Research Group, University of Pannonia, Veszprém Boglárka Bukor Behavioral Ecology Research Group, Center for Natural Sciences, University of Pannonia, Veszprém Nóra Ágh MTA-PE Evolutionary Ecology Research Group, University of Pannonia, Veszprém András Liker MTA-PE Evolutionary Ecology Research Group, University of Pannonia, Veszprém


Introduction
Urbanization often has negative effects on wildlife and the dramatic expansion of urban areas is a major factor responsible for the severe decline of animal populations worldwide (Maxwell et al. 2016, Ibáñez-Álamo et al. 2017, Chaudhary et al. 2018. One of the key features of urban environments is the ubiquitous presence and activity of humans that generate various stressors for urban wildlife, including chemical, acoustic, and light pollutions associated with vehicle tra c and direct disturbances caused by pedestrians and recreational activities in urban green spaces (Forman 2014).
To impede the spreading of the coronavirus disease (COVID-19), many countries went into lockdown during the spring of 2020, resulting in greatly reduced human mobility and tra c in cities and towns across the world for several months, a phenomenon coined as 'anthropause' (Rutz et al. 2020). Although observations suggest that urban wildlife quickly responded to the anthropause, for example by roaming in cities more freely than before (Rutz et al. 2020) or by adjusting their behavior to the new circumstances (Derryberry et al. 2020, Jiang et al. 2020, Soh et al. 2021, quantitative assessments of the ecological impacts are still largely lacking. In the northern temperate zone, the duration of the lockdown coincided with the breeding season of many animal populations, hence providing excellent opportunities to investigate the tness consequences of human activities in urban wildlife populations. Under usual circumstances, these impacts are typically inferred from comparisons between populations living in habitats with different levels of urbanization or along urban-to-rural gradients (McDonnell and Hahs 2008), because conducting experiments on relevant spatial and temporal scales in cities is challenging (but see Marzluff et al. 2016). Although some components of the urban environments can be effectively manipulated by experiments, e.g. food availability (Seress et al. 2020), noise levels (Mulholland et al. 2018), or light pollution (Welbers et al. 2017), other factors, such as road tra c and the ubiquitous presence of humans and the associated direct and indirect disturbances, are much more di cult to control in cities. For this reason, we have limited and mostly correlative evidence on the impacts of these anthropogenic effects on the tness of urban animals.
Thus, in this study, we exploit the opportunity that the unprecedented conditions during the lockdown offered, and use a quasi-experimental approach to study the impacts of human activity on the breeding success of urban animals. To do so, we compare the reproductive performance of a widespread urban adapter bird, the great tit (Parus major) between areas where the lockdown induced markedly different effects on human activities. In a long-term study of this model system, we already showed that urban great tits have strongly reduced breeding success and nestling development relative to forest populations (Seress et al. 2018(Seress et al. , 2020. One main factor responsible for this disparity in breeding success is the limited availability of arthropod-rich nestling food (especially caterpillars) in urban areas (Seress et al. 2018(Seress et al. , 2020, although several other mechanisms may contribute to this phenomenon. For example, increased levels of tra c-related air pollution may affect urban birds detrimentally either via their health status (Salmón et al. 2018) or by reducing the abundance of their arthropod food supply (Jones and Leather 2012). Also, persistent exposure to anthropogenic noise may induce chronic stress (Kleist et al. 2018) and promote nest abandonment (Strasser and Heath 2013), and can acoustically mask parent-offspring communication and other relevant auditory cues in birds' environment resulting in decreased tness (Schroeder et al. 2012). Finally, high levels of human disturbance associated with, for example, human recreational activity in urban green spaces may also disrupt parents' nestling feeding patterns or constrain their access to foraging sites (Fernández-Juricic andTellería 2000, Merrall and. Therefore, here we test the prediction that the large scale (city-wide) and persistent reduction in human activities due to the lockdown during the birds' breeding period should facilitate their reproductive success in urban areas (Corsini et al. 2017). To do so, we combined spatial and temporal comparisons.
First, we compared the birds' reproductive performance between a city (Veszprém) and a nearby mature forest area (Szentgál) that were affected differently by the movement restrictions during the spring of 2020, and compared this to 2019 as a reference (a year with no lockdown). While in 2020 the lockdown resulted in a signi cantly lower number of humans counted around the urban nests in Veszprém (compared to a long-term reference period), forest nests remained relatively undisturbed, similarly to the preceding years (Fig. 1a). Therefore, we predicted that the anthropause would yield positive effects on the birds' reproductive success in the urban but not in the forest area, resulting in a reduced difference between the forest and the city in 2020 compared to their difference in 2019.
Second, we compared the birds' reproductive success between two cities, Veszprém and Budapest (the capital city of Hungary), which responded in strikingly different ways to the lockdown. Interestingly, while the lockdown reduced human density in the public green areas of Veszprém by 44% (see above), we recorded a 40% increase in human presence around nest sites in the urban parks of Budapest, where residents used green spaces for recreational activities more frequently than in the preceding year (the levels of human presence were similar in the two urban areas in 2019; Fig. 1b). Building on this observation, we hypothesized that if human disturbance is indeed a major factor affecting reproductive success (Hut uss and Dingemanse 2019), then birds should have a better breeding performance in Veszprém relative to Budapest in 2020 than in the pre-pandemic year (2019).

Data collection
Data on the birds' reproductive success and the number of humans present at nest sites were collected as part of a long-term, ongoing monitoring project in Hungary, in which we investigate the impacts of urbanization on populations of great tits. The great tit is an insectivorous passerine bird that is widespread across the Western Palearctic, occupies both urban and forest habitats, readily accepts nestboxes, and shares many important ecological traits with other tit or chickadee species also occurring in urban habitats (del Hoyo et al. 2008). These traits make this species an ideal model organism for studying the effects of the anthropause on wildlife in different environments.
Study sites. We monitored breeding great tit populations and also collected human presence data in two urban areas and at one forest study site. In one of the urban sites, Veszprém (47°05′17.29″N, 17°54′29.66″E; human population: c. 56,000; the monitoring scheme started in 2013), the nestboxes were placed in public green spaces (public parks, university campuses, a bus station, and a cemetery) that are surrounded by built-up areas and roads, and experience frequent anthropogenic disturbance. At the other urban size, Budapest (47°30'27.4"N,19°01'03.4"E; the capital city of Hungary, human population: c. 1.75 million; the monitoring scheme started in 2019), the nestboxes were placed in two public urban parks, located c. 400 m from each other in the city core area and separated by high-tra c roads. The parks are freely accessible to residents and are heavily embedded within the urban matrix. At both urban sites, most of the nestboxes are distributed along paths or walking trails. Even though the two cities greatly differ in their size and human population, our urban study plots in both cities have similar general characteristics: these are surrounded by built-up areas, are at a similar distance (c. 3-4 km) from the nearest forested areas (for Veszprém, this is the forest at Vilma-puszta: 47°05′06.7″N, 17°51′51.4″E; for Budapest, this is the forest at Normafa: 47°30'27.7"N 18°57'51.1"E), and nests also experienced a similar level of human disturbance in the pre-COVID reference period (Fig. 1b). The forest site, Szentgál (47°06'39.75"N, 17°41'17.94"E; the monitoring scheme started in 2013), is a mature woodland, dominated by beech (Fagus sylvatica) and hornbeam (Carpinus betulus), located 3 km away from the nearest human habitation (Szentgál, human population: c. 2.800), c. 20 km and 110 km away from Veszprém and Budapest, respectively. There are no paved roads in the forest, and the area is relatively free from human disturbance although it experiences occasional hunting and logging activity.
Human presence around nests. To quantify human presence at our study sites for 2020 and the reference years we counted the number of humans (motorized vehicles excluded) during each nest check, for 30 sec, in the proximity of the nestboxes (for similar approach see Corsini et al. 2019). The number of humans was recorded within a 50-m radius of the nestboxes between 2013-2018 (Veszprém, Szentgál), and within a 15-m radius distance in 2019-2020 (all sites). We changed the counting distance in 2019 due to methodological reasons following Corsini et al. 2019. However, to be able to compare the human presence data of 2020 in Veszprém and Szentgál to that recorded in earlier years, in 2020 we performed the counts with both the 15-m and the 50-m radius distances at these two sites. Thus, for 2020 in Veszprém, we have human presence data both for the 50-m and the 15-m radius areas that were used in the forest-city and the between-cities comparisons, respectively (see below). For each year and study site, we used human presence data only from seasonally rst broods (de ned below), and only from nests where there were either already eggs or nestlings in the nest, resulting in 9.4 ± 3.6 (mean ± SD) observations per brood which is a reliable indicator of human presence (Corsini et al. 2019).
Birds' reproductive success. We monitored nestboxes each year at least twice a week from mid-March to early June to record laying and hatching dates, clutch size, hatching success, and the number of nestlings in active great tit nests. We ringed nestlings at day 14-16 post-hatch (i.e. a few days before edging; hatching day of the rst chick = day 1) with a numbered metal ring and also recorded their body mass (to the nearest 0.1 g), tarsus length (to the nearest 0.1 mm and following Svensson's 'alternative' method; Svensson 1992) and wing length (from the bend of the wing to the longest primary; to nearest 1 mm). Shortly after edging we carefully examined the nest material to identify and count the number of chicks that died after ringing that we included in the calculation of nestling survival.
From these data we determined clutch size (the maximum number of eggs observed in a brood), the number of edglings, and nestling survival (the proportion of edged young/hatched chicks). The number of edglings (i.e. the number of young edged successfully) was calculated as the number of chicks ringed minus the number of chicks found dead in the nest after the ringing. We involved only seasonally rst breeding attempts, de ned as follows. In our study system breeding great tits are captured on their nests and receive a unique combination of colour rings. Active nests are also routinely equipped with a small, concealed video camera enabling us to reliably identify over 80% of breeding individuals each year (Seress et al. 2017). Thus, relying on this setup, we considered a clutch as a rst breeding attempt of a pair if it was initiated before the date of the rst egg laid in the earliest second clutch at that site by an individually identi able (i.e. colour-ringed) female that successfully raised her rst clutch (i.e. edged at least one young) in that year.
Air pollution and meteorological conditions. To describe the levels of tra c-related air pollution (nitrogen dioxide [NO 2 ], nitrogen oxides [NO X ] and ozone [O 3 ]) and the meteorological conditions (temperature and precipitation) at the two urban study sites (Veszprém and Budapest), we used data provided by the Hungarian Air Quality Monitoring Network and the Hungarian Meteorological Service, respectively. To better understand which aspect of the anthropause might have affected great tits' breeding success we thus assessed if the lockdown affected air pollution levels differently at the two urban study sites (compared to 2019), or if weather conditions showed different uctuations between 2019 and 2020 at the two cities. For more details on the statistical analyses and results, see ESM: section 1.

Statistical analyses
The duration of the o cial restrictions on human mobility (lockdown) spanned between 28 March-4 May in Veszprém (Julian date: 88-125) and 28  in Budapest. During this period people were allowed to leave their homes e.g. to run essential errands including individual sport and recreational activities in public green spaces, although with keeping at least 1.5 meters from each other (social distancing). Very importantly, from the point of view of our study, the period of movement restrictions had almost completely overlapped with the seasonally rst breeding attempts (from egglaying to edging) of great tits at both urban sites. The date of laying the 1 st egg (Julian date, mean ± SD) in Veszprém was 94.2 ± 6.4, while in Budapest 97 ± 7.8; the date of chick ringing and measuring in Veszprém was 128 ± 5.3, while in Budapest: 133 ± 9.1. Thus, we decided not to exclude any rst broods based on the date in order to maximize our sample size. Similarly, the period from which we involved human presence data was also strongly overlapped with the duration of the movement restrictions in both cities. Therefore, the Julian dates of the rst and the last human count at each nest, respectively, were as follows: in Veszprém 87-108 (median: 100) and 121-142 (median: 132), while in Budapest 87-128 (median: 98) and 118-155 (median: 128).
Human presence around nests. In accordance with our rst objective (forest-city comparisons), we explored if the lockdown in 2020 caused any changes in human disturbance around the great tit nests. To do so, we compared the number of humans (50-m radius of the nests) between 2020 and the 2013-2018 reference period, separately for the forest (Szentgál) and urban (Veszprém) study sites. Note that in 2019, we did not collect data on human presence within a 50-m radius at Veszprém and Szentgál (see above: Data collection), therefore 2019 was not included in the reference period of this analysis. We, however, also compared human presence in Veszprém between 2019 and 2020 using the 15-m radius data which indicates a change that is consistent with the differences found using the 50-m radius data (detailed below).
First, we built generalized linear mixed-effects (GLM, lme4 R package) models with Poisson error distribution with the number of humans as the response variable, including year as a xed factor and nestbox ID as random factor to control for non-independence of the data. Next, we extracted the mean values (least-squares means; package emmeans (Lenth et al. 2020)) and associated standard errors for each year as estimated by the model. We computed the mean of these yearly mean estimates for the 2013-2018 reference period (i.e. calculated a single overall mean describing the whole reference period) and compared this long-term mean to the mean estimate of 2020 by calculating the linear contrast between them (with the 'contrast' function of the emmeans package), and expressed linear contrasts as 2020 minus the reference period.
For our second objective (between-cities comparisons), we compared the changes in human disturbance around the nestboxes at the two urban study sites, Veszprém and Budapest, using the number of humans recorded within the 15-m radius of the active nests in 2019 and 2020. We analysed the data from Budapest and Veszprém separately and built generalized linear mixed-effects models with Poisson error distribution with the number of humans (15-m radius of the nests) as the response variable, including year as a xed factor and nestbox ID as random factor to control for non-independence.
Birds' reproductive success. We used data from 2019 (reference; for justi cation see below in this section) and 2020 (lockdown). First, we constructed separate linear models to analyse each component of reproductive success (response variables), and for the forest-city and the between-cities comparisons. We used linear models (LM) for clutch size and the number of edglings, generalized linear models (GLM, with quasi-binomial error distribution) for nestling survival, and linear-mixed effects models (LME) for nestling body size traits (body mass, tarsus length, and wing length). Models on nestling body size traits contained nestlings' age at ringing as a confounding variable (three-level factor: 14, 15, or 16 d of age) and brood ID as a random factor to control for the non-independence of chicks raised in the same brood. Finally, these models always contained a habitat (Veszprém or Szentgál) × year (2019 or 2020) interaction term for forest-city comparisons and a city (Budapest or Veszprém) × year (2019 or 2020) interaction term for between-cities comparisons. We checked assumptions of residuals' normality and homogeneity of variance by inspecting the residuals plots which were respected for all models.
Next, to test the prediction for our rst objective (forest-city comparisons), we extracted the mean values (least-squares means) and associated standard errors of each response variable for each habitat × year combination as estimated by the linear model's interaction. Then, from these estimates, we calculated habitat contrasts, i.e. the mean forest-city difference (forest minus urban) for each year (i.e. for 2019 and 2020), and compared the mean habitat contrast for the 2019 reference year to the mean habitat contrast of 2020; for similar approach see Vincze et al. 2019, Seress et al. 2020 For our second objective (between-cities comparisons), we followed the same procedure as for the forestcity comparisons (detailed above) except that here we compared the differences between cities (Budapest minus Veszprém) in 2020 and 2019. These full models (i.e. for the forest-city and betweencities comparisons) are presented in Table S1-S2 (ESM: section 2).
In our study, we chose 2019 as a reference year for multiple reasons. First, because this was temporally the closest year without a lockdown. Second, because for Budapest we have monitoring data only from 2019-2020, using 2019 and 2020 in all analyses makes the results more comparable. Finally, although we have monitoring data from a total of eight years (2013-2020) for Szentgál (forest site) and Veszprém (urban site), for the forest-city comparisons we did not include years before 2019 in the reference because we noticed a negative trend in birds' reproductive success throughout the study years (Fig. S3). This trend was especially apparent in the forest population, and may have reduced the forest-city difference by the end of the study period. Indeed, 2019 and 2020 were amongst the poorest years and resulted in a very similar reproductive success between both years within both habitats (Fig. S3). Because such temporal trend may have confounded the comparisons of 2020 with earlier years, to take account for its effect, and to further justify our approach of using 2019 as the reference year, we conducted additional analyses on birds' reproductive success by comparing both 2019 and 2020 to the 2013-2018 long-term reference period. We predict that if 2019 and 2020 are similarly affected by the decreasing trend in reproductive success than then the differences between the long-term reference period and 2019 and 2020, respectively, should be similar. For the details of these long-term forest-city comparisons see ESM: section 3 and Table S3).
Clutches that failed before reaching the incubation stage (due to predation or desertion; i.e. nal clutch size was uncertain), suffered complete mortality due to weather (nestbox fall from the tree), and cases when complete or partial clutch or brood loss may have occurred due to the monitoring process (e.g. when a nest box fell or when chick mortality occurred soon after capturing a parent on the nest) were excluded from all analyses. In the analyses investigating the number of edglings, edging success, and nestling body size traits we involved nests only where at least one nestling hatched, and excluded broods that were involved in a food-supplementation experiment (as treatment group) during the nestling rearing period in 2017 (Seress et al. 2020). All analyses were implemented in the R 3.5.2 software environment (R Core Team 2016).

Results
Human presence around nests. Regarding our rst objective (forest-city comparisons), our results indicated that, for Veszprém, the number of humans around the nests was signi cantly lower in 2020 than in the 2013-2018 reference period (linear contrast, mean ± SE= -0.98 ± 0.06, z-ratio= -17.132, P< 0.001, n = 2807 observations; Fig. 1a). In contrast, in Szentgál, we detected no such differences between 2020 and the previous years (linear contrast, mean ± SE= 0.52 ± 0.44, z-ratio= 1,167, P= 0.243, n = 2368 observations; Fig. 1a), as the number of humans registered within the 50-m radius of the nestboxes remained very low throughout the rst-brood period similarly to the reference years (at Szentgál 98.4 % (in 2020(at Szentgál 98.4 % (in ) and 99.0 % (2013(at Szentgál 98.4 % (in -2018 of the censuses indicated no humans around the nestboxes).
For our second objective (between-cities comparisons), we found that the number of humans around the nests in 2020 was signi cantly lower in Veszprém (b ± SE= -0.48 ± 0.13, z-value= -3.761, P < 0.001, number of observations: 695) and signi cantly higher in Budapest (b ± SE= 0.42 ± 0.10, z-value= 4.103, P< 0.001, number of observations: 411) than in 2019 (Fig. 1b). Parameter estimates are expressed as 2020 compared to 2019, i.e. a negative value indicates a lower number of humans in 2020.
Birds' reproductive success. Regarding our rst objective and contrary to our expectation, we found no signi cant changes from 2019 to 2020 in the forest-city differences in any components of great tits' breeding success and nestling body size traits ( Fig. 2; Table 1a), and the site × year interaction was also non-signi cant for any of the response variables (Table S1). Combined together, these results mean that the lockdown period in 2020 was not associated with an increased reproductive success in the city relative to the forest population. Our additional results for the forest-city comparisons using the long-term data (2019 and 2020 vs. 2013-2018) revealed that the forest-city differences in the birds' reproductive success between 2019 and the long-term reference period was similar to those between 2020 and the long-term reference period (Table S3). This latter result indicates that (1) the long-term temporal trend in reproductive success is unlikely to be related to the lockdown-effect, and (2) justi es further the similarity between 2019 and 2020, despite that lockdown occurred only in 2020.
For our second objective, the between-cities comparisons, we found signi cant changes in the betweencities differences in nestling body size traits (as measured by body mass, tarsus length, and wing length) between 2019 and 2020 ( Fig. 2a-c; Table 1b), in line with our prediction. These changes were driven by birds' differing response to the lockdown in the two cities: while in Budapest nestlings reached signi cantly smaller body mass and size in 2020 than in 2019, in Veszprém we found no differences in nestlings' body size traits between the two years (signi cant site × year interaction in all nestling traits, see Table S2). In contrast, we found no similar temporal changes in the between-cities differences for clutch size, the number of edglings, or nestling survival (Table 1b; non-signi cant site × year interactions, see Table S2).

Discussion
Our results provide two key implications. First, as suggested by comparisons between a forest and a smaller city (Veszprém), the lockdown did not eliminate urban birds' disadvantage in reproductive success relative to a non-urban population, despite the considerably reduced human activities in the city during the breeding season of 2020. In theory, this lack of effect might be related to the relatively low levels of pre-lockdown human disturbance (tra c intensity, pedestrian density) in smaller cities like Veszprém. However, this explanation is not likely, as in 2019 the number of humans recorded around the nests in Veszprém was very similar to that recorded in the urban parks of Budapest (Fig. 1b), although the behavior and activity of people (hence the intensity of disturbance) may differ between the two cities. An alternative and in our opinion more likely explanation for the lack of lockdown-effect in Veszprém is that the difference in reproductive success between the city and the nearby forest site (Szentgál) is mainly driven by differences in food availability with relatively low contribution from other anthropogenic effects. This interpretation is strongly supported by the considerably lower abundance of caterpillars (the principal nestling food for great tits) at our urban compared to forest sites (Seress et al. 2018), and especially by a food supplementation experiment that eliminated much of the forest-city differences in both nestling body size traits and nestling survival (Seress et al. 2020).
Second, nestlings' body size traits were strongly reduced during the lockdown in the parks of Budapest where recreational activities signi cantly increased human presence (hence probably the intensity of disturbance to birds) in the proximity of nests compared to 2019. Importantly, this hindered nestling development might have posed tness consequences for birds in Budapest, as reduced pre-edging body mass predicts lower rst-year survival in many songbird species including great tits (Magrath 1991). These results provide one of the rare large-scale quasi-experimental supports for the hypothesis that human disturbance can pose considerable tness costs for animals even at habitats where high level of human presence is common, and even for such a successful urban-dweller species like the great tit. The negative impacts on nestlings' body size traits may be caused by the adverse effects of increased disturbance in several, non-exclusive ways. Human disturbance can act as a stressor resulting in elevated stress hormone levels (e.g. glucocorticoid concentrations) in animals, helping them to cope with stressful conditions at the cost of suppressing other activities like reproduction or parental care (Wing eld and Sapolsky 2003) as it was demonstrated in other urban-breeding bird species (Strasser and Heath 2013).
Similarly, increased human activity and its by-products (e.g. noise) can also have subtle impacts on birds' behaviour. For example, human disturbance can disrupt birds' foraging patterns, hindering parentoffspring communication (Schroeder et al. 2012), constraining birds' access to preferred foraging sites (Fernández-Juricic and Tellería 2000), or reducing parents' foraging success, which has also been shown in other species (Suraci et al. 2019, Merrall and. These impacts can be especially strong in food-limited environments. In line with this idea, a recent meta-analysis found a massive decline in the abundance and diversity of arthropods, especially for lepidopterans, in temperate zone cities (Fenoglio et al. 2020). Under such conditions, a sudden and lasting increase in another stressor (e.g. human disturbance) may further limit the birds' feeding success and hinder the already impaired nestling development. Note that the opposite may not be necessarily true -as illustrated by that great tits at Veszprém could not realize higher breeding success in 2020 (compared to earlier years) under signi cantly reduced human disturbance around their nests. Finally, in contrast with nestlings' body size, we found no temporal changes in the between-cities differences for clutch size, the number of edglings, or nestling survival. These latter results might be explained, at least to some extent, by the relatively short period elapsed from the beginning of the lockdown to the initiation of the egg-laying (on average c. 6 days (Veszprém) and 9 days (Budapest); see Methods) that might have reduced the potential impacts of the changes in human disturbance on clutch size (and perhaps also on edgling number).
Besides inducing changes in human activity, the lockdown had other consequences that could have affected urban birds' breeding by different means, although we argue that these may be less likely in our case. For example, the reduction in vehicular road tra c during the lockdown was paralleled by reduced concentrations of air pollutants, resulting in improved air quality worldwide (Salma et al. 2020) that could have bene cial effects on urban birds directly e.g. by affecting their oxidative stress physiology (Salmón et al. 2018) or indirectly (via their arthropod food supply). However, while in Budapest the levels of NO 2 and NO X (two reliable indicators of tra c-related urban air pollution, measured during birds' rst-brood period) were substantially reduced during the lockdown compared to 2019 (Fig. S1) or longer reference periods (Salma et al. 2020), in Veszprém the levels of both pollutants were similar in 2019 and 2020 (Fig.  S1). Note that these different changes in air pollution in the two cities would predict an opposite change in birds' reproductive performance of what we detected in our urban populations. Decreased vehicle tra c can, however, result in actually elevated levels of another air pollutant, the ground-level O 3 . Indeed, O 3 concentrations were found to be strongly increased in several cities during the lockdown due to the unprecedented reductions of NO X emissions (Sicard et al. 2020). High levels of O 3 have several harmful impacts on birds' reproductive success, and a recent long-term, continental-scale study demonstrated negative associations between bird abundance and O 3 concentrations in the United States . Even though in our study system the levels of O 3 were higher in Veszprém than in Budapest (both in 2019 and 2020), its concentrations did not change signi cantly in either city during the lockdown when compared to the 2019 reference period (Fig. S1). Therefore, building on these results, we conclude that it is not likely that the lockdown-associated changes in the levels of the studied air pollutants can explain the differences in birds' reproductive success between the cities. Finally, we can also exclude that differences in meteorological conditions at our two urban study sites biased our results because both daily temperature and the amount of daily precipitation (measured during birds' breeding period) changed similarly at the two urban study sites between 2019 and 2020 (ESM: section 1.2 and Fig. S2).
In conclusion, our study demonstrates that reductions in human presence during the lockdown had little effect on great tits' reproductive performance in urban habitats, where other factors, e.g. food availability, are known to strongly limit the birds' reproduction. On the other hand, our study also showed that the lockdown could result in strongly increased human presence in some urban areas to a level that has detrimental tness consequences for urban-breeding birds, even if the lockdown otherwise considerably reduced other negative anthropogenic effects like air and noise pollution. Taken together, these results suggest that, in ecological systems like ours, it would need more than just a few months of 'anthropause' to offset the detrimental impacts of urbanization, as the positive impacts resulting from the reduced human mobility very likely need more time to build-up and take effect on urban wildlife.
Declarations ACKNOWLEDGEMENTS: We are thankful to all who participated in eld data collection in 2013-2020. GS is thankful to Prof. Ádám Miklósi because without his provocative attitude this manuscript may have never been written. Table   Table 1. Changes in population differences in the reproductive success of great tits between 2019 (reference period) and 2020 (breeding during the lockdown). Population differences were calculated either (a) between a forest and an urban habitat (Szentgál minus Veszprém) or (b) between two cities (Budapest minus Veszprém). Negative contrasts indicate smaller between-population differences during the 2020 lockdown than in 2019. Model estimated contrast ± SE values are on the original scales except for nestling survival (the proportion of fledged young/hatched chicks), which is given on the log odds ratio scale. Statistically significant results (P < 0.05) are highlighted in bold. See Table S1 and Table S2 for statistical details of the models used for the forest-city and between-cities contrast calculations, respectively.
(a) Comparisons between forest and urban habitats diff. in 2020 -diff.    <p>The impacts of the lockdown on nestlings&rsquo; body mass, tarsus length, and wing length in urban and forest populations of great tits. Black colour indicates comparisons between a forest (Szentg&aacute;l, not affected by the lockdown) and a nearby urban site (Veszpr&eacute;m, where the lockdown decreased human activity around the nests), whereas orange colour shows comparisons between Veszpr&eacute;m and the capital city, Budapest (where the lockdown increased human activity at the study sites). Panels show the differences between these populations for 2020 (when breeding