A high-resolution prediction system for birch pollen in Sapporo

https://doi.org/10.1016/j.agrformet.2020.108229Get rights and content

Highlights

  • A high-resolution birch-tree map was created by a field survey and street viewing.

  • Birch pollen prediction system represented the emission, advection and deposition.

  • The prediction was succeed to reproduce the observation at a site.

  • Birch pollen concentration is highly correlated with birch tree density.

Abstract

This study aims to develop a high-resolution forecast model for airborne birch pollen concentration in Sapporo. Using results that shows that the bulk of pollen grains was deposited near the source plant, we created a high-resolution birch tree density map from a field survey along public roads and footpaths and Street View on Google Maps. The map was linearly extrapolated to cover an area 10 km from the city centre, including forests in the south and southwest. A new atmospheric dispersion deposition model for birch pollen was developed to represent wind advection, turbulent mixing, and convection in the atmospheric boundary layer, and gravitational settling and scavenging with precipitation. Pollen emission was also modelled as a function of the birch tree density and surface air temperatures. The onset of the pollen season was based on the pollen observations by a sampler at Hokkaido Institute of Public Health in Sapporo, and the end of pollen season was assumed to be 40 days after its onset. Comparison with observations suggested that the hindcast simulation with the model reproduced the daily variations in the dry deposition well in most years from 2001 to 2011. The geographical distribution of birch pollen deposition was highly correlated with the birch tree density map, even on windy days.

Introduction

Hokkaido is the northernmost main island of Japan (Fig. 1a) with a population of about 5.3 million, and is an area with a high prevalence of birch pollinosis, which causes springtime allergy-related symptoms, such as stuffy nose, watery eyes and sneezing, and asthmatic reactions (Gotoda et al. 2001; Baba and Nakae 2008). About one-third of the population is concentrated in Sapporo, the capital of Hokkaido Prefecture. Residential and commercial areas extend north and east over the plain about 10 km from the centre, and there are mountainous areas to the south and west, with Mt. Teine in the west, Mt. Moiwa in the southeast, and Mt. Yakeyama in the south (Fig. 1b). Because birch is a shade-intolerant species, there are clusters of Betula platyphylla var. japonica, B. ermanii, and B. maximowicziana around the summit and in the ridgelines of the mountains. Moreover, there are birch trees, mostly B. platyphylla var. japonica, planted along major roads and in city parks administered by the Greenery Promotion Department of Sapporo City. Given that we expected that the bulk of pollen would be deposited near the source plant (Bricchi et al. 2000; Tampieri et al. 1977; Wright 1953), the instantaneous amount of airborne birch pollen was expected to be sensitive to the density of birch trees (Skjøth et al. 2008). Hence, the prediction of daily or hourly amount of pollen in Sapporo requires primarily a high-resolution map of birch tree distribution, which would help pollinosis patients take preventive action (Peden and Reed 2010). However, no map has been compiled because the trees are distributed irregularly in roads, parks, estates, and forests owned by various organisations and individuals.

The prediction of pollen concentration has long been a subject of study and requires a dynamical model system that accounts for the total amount, flowering, emission, dispersion, and deposition. The total amount of pollen during a single season has often been evaluated by seasonal pollen index, a proxy for allergenic exposure of the population (D'Amato et al. 2007; Gioulekas et al. 2004). The magnitude of flowering varies between years and is affected by the environmental conditions and the flowering intensity of the previous year (Dahl and Strandhede 1996; Ranta et al. 2008; Stach et al. 2008a). The total amount of pollen exhibits a pronounced biennial rhythm in Sapporo (Kobayashi et al. 2013). The total amount of pollen is also highly correlated with the number of catkins (Ranta et al. 2008), and can be predicted the year before (Yasaka et al. 2009). Shirasaki et al. (2014) found a significant positive correlation with the sunshine hours in June of the preceding year. Birch flowering timing is driven by heating degree days relative to a base temperature of 3.5°C from 1 March (Larsson 1993; Linkosalo et al. 2010). The onset date ranges from 70 to 120 degree days in central Europe (Sofiev et al. 2013; Verstaeten et al. 2019), and a pollen calendar model has been developed for predicting daily concentrations of airborne pollen in Serbia (Šikoparija et al. 2018). High temperatures trigger pollen emission into the air (Clot 2001; Ribeiro 2003) and the diurnal cycle of temperatures and sunshine also affects pollen emission (Fehér and Járai-Komlódi 1998). Based on these results, a dynamical model system for pollen prediction was developed using a conventional transport-deposition model for aerosol particles (Sofiev et al. 2006a, 2013). The dynamical model predicted the transboundary pollen transport effectively (Siljamo et al. 2012), and Sofiev et al. (2015) performed the first ensemble modelling experiment on birch pollen in Europe. However, because up to 1% of released pollen grains can travel far from the source despite the large dry deposition velocity of around 1.2 cm s−1 (Sofiev et al. 2006b), the accuracy of pollen prediction strongly depends on information about plant distribution.

Inatsu et al. (2014) reported a simulation of daily birch pollen in Sapporo. They constructed a multiple regression model with two independent predictors of lagged pollen amount and diurnal temperature range (Kawashima and Takahashi 1995), based on observed data from 2001 to 2010 in Sapporo (Kobayashi et al. 1998). A hindcast simulation for 2011 reproduced the daily variation in pollen amount with a correlation coefficient with the observed data of 0.84. However, using the pollen amount from the previous day (Inatsu et al. 2014; Stach et al. 2008b) may not be practical because the current day's pollen data are not available in time to produce tomorrow's pollen forecast. Ritenberga et al. (2016) excluded this nudging predictor by preprocessing, including normalisation by the seasonal pollen index, using heat sum as a time axis, projecting meteorological variables to pollen concentrations, and removing the mean pollen season curve. Nevertheless, no matter how complex the system, statistical prediction only provides the pollen amount as a relative value that optimally fits training data at a particular site. Because most of the released pollen grains settle near the source, the pollen amount slightly away from the site must differ substantially from the statistically predicted amount at the site, although these two amounts may be strongly correlated with each other. Moreover, other aerosol studies (Neubauer et al. 2014) suggest that the spatial distribution, although dependent on the tree density distribution, may have a much larger effect than the temporal variation. Hence, the statistical models for a single site may not be suitable for evaluating the higher exposure risk around birch trees. Even though the absolute amount of airborne pollen grains is helpful for pollinosis patients, the spatial distribution has not been examined, probably due to the lack of a high-resolution birch tree map.

The purpose of this study is to develop a high-resolution prediction system for birch pollen in Sapporo, including mapping the tree density and modelling the airborne birch pollen transport and deposition. A crucial stumbling block for such a prediction is the lack of a high-resolution birch tree distribution map. About one-third of birch pollen is deposited within 800 m of the source (Bricchi et al. 2000); thus, we first create a map of Sapporo city centre with 400 m spatial resolution based on a field survey along public roads and footpaths and Street View on Google Maps (Section 3). Second, we originally develop a transport-deposition model driven by uniform meteorological data for the city (Section 4). The model includes all the components necessary for birch pollen prediction, but it should be as simple as possible to facilitate a practical forecast. Third, we perform a hindcast simulation from 2001 to 2011 in Sapporo and compare the result with the pollen observations at a site (Section 5). We also compare the spatial variation with the temporal variation in the simulation result, which has not been reported in previous studies. Finally, we discuss the improvement of our prediction system (Section 6).

Section snippets

Data

We used the data for airborne pollen amount between 2001 and 2011 observed at the Hokkaido Institute of Public Health (HIPH) in Sapporo, which is located at 43°05′N, 141°20′E (Fig. 1b). The pollen was collected on petrolatum-coated glass slides in a Durham sampler installed on the rooftop 16 m above ground level (AGL). Durham samplers have a comparable performance to Burkard samplers (Kishikawa et al. 2009). After the slides were stained with gentian violet, the number of pollen grains was

Birch tree map

The horizontal distribution of birch trees was mapped onto a 400-m horizontal grid within 10 km of Sapporo city centre (Fig. 1b), based on a field survey, a search by Google Street View and Google Maps and an extrapolation. The procedure is illustrated in Fig. 3. First, we generated about 36,000 1-ha blocks in the circle with MANDARA (Geographic Information System, http://ktgis.net/mandara/). An author (RY) manually counted the number of birch trees along the city roads in a block in the

Prediction system

The prediction system we developed contains three components: flowering time, emission, and transport and deposition.

Hindcast simulation

We performed the hindcast simulation with the model described in Section 4. The target years were 2001 to 2011, and the daily deposition amount of the simulation at the grid cell containing HIPH was compared with the daily amount of pollen grains (unit of cm−2 day−1) collected in the Durham sampler there (Fig. 8)3

Conclusion and Discussion

We developed a high-resolution prediction system for birch pollen in Sapporo. We created a birch tree density map within 10 km of Sapporo city centre in 16 ha blocks by a field survey along public roads and footpaths and using Google Maps Street View (Section 3). We developed a system consisting of components such as flowering time, emission, and transport and deposition (Section 4). To facilitate practical forecasting, the model was designed to be as simple as possible, based on phenological

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

We thank the two anonymous reviewers for providing us with insightful comments that helped improve this paper. We are deeply indebted to Hokkaido Agricultural Research Center of National Agriculture and Food Research Organization, Hokkaido Research Center of Forestry and Forest Products Research Institute, Hokkaido Agricultural Technical College of Hakkou Gakuen, and Sapporo Development Corporation for assistance with the field survey of birch tree distribution in their domain. The Greenery

References (42)

  • G D'Amato et al.

    Allergenic pollen and pollen allergy in Europe

    Allergy

    (2007)
  • J W Deardorff et al.

    Further results from a laboratory model of the convective planetary boundary layer

    Bound-Layer Meteor

    (1985)
  • Z Fehér et al.

    A new weather factor predicting airborne pollen concentration: Péczely's macrosynoptic weather types

    Aerobiologia

    (1998)
  • D Gioulekas et al.

    Allergenic pollen records (15 years) and sensitization in patients with respiratory allergy in Thessaloniki, Greece

    Allergy

    (2004)
  • H Gotoda et al.

    Spring time pollinosis and oral allergy syndrome in Sapporo

    Auris Nasus Larynx

    (2001)
  • M Inatsu et al.

    Statistical analysis on daily variations of birch pollen amount with climatic variables in Sapporo

    SOLA

    (2014)
  • M Inatsu et al.

    Toward predicting expressway closures due to blowing snow events

    Cold Region Science and Technology

    (2020)
  • S Kawashima et al.

    Modelling and simulation o mesoscale dispersion processes for airborne cedar pollen

    Grana

    (1995)
  • R Kishikawa et al.

    Japanese Cedar Airborne Pollen Monitoring by Durham's and Burkard Samplers in Japan ― Estimation of the Usefulness of Durham's Sampler on Japanese Cedar Pollinosis ―

    Global Environ Res.

    (2009)
  • S Kobayashi et al.

    Airborne pollen survey in North Ward, Sapporo in 1997

    Report of Hokkaido Institute of Public Health

    (1998)
  • S Kobayashi et al.

    Trends in annual counts of Betula pollen from six cities in Hokkaido

    Japanese J Palynology

    (2013)
  • Cited by (0)

    View full text