Review: advances in in situ and satellite phenological observations in Japan

To accurately evaluate the responses of spatial and temporal variation of ecosystem functioning (evapotranspiration and photosynthesis) and services (regulating and cultural services) to the rapid changes caused by global warming, we depend on long-term, continuous, near-surface, and satellite remote sensing of phenology over wide areas. Here, we review such phenological studies in Japan and discuss our current knowledge, problems, and future developments. In contrast with North America and Europe, Japan has been able to evaluate plant phenology along vertical and horizontal gradients within a narrow area because of the country’s high topographic relief. Phenological observation networks that support scientific studies and outreach activities have used near-surface tools such as digital cameras and spectral radiometers. Differences in phenology among ecosystems and tree species have been detected by analyzing the seasonal variation of red, green, and blue digital numbers (RGB values) extracted from phenological images, as well as spectral reflectance and vegetation indices. The relationships between seasonal variations in RGB-derived indices or spectral characteristics and the ecological and CO2 flux measurement data have been well validated. In contrast, insufficient satellite remote-sensing observations have been conducted because of the coarse spatial resolution of previous datasets, which could not detect the heterogeneous plant phenology that results from Japan’s complex topography and vegetation. To improve Japanese phenological observations, multidisciplinary analysis and evaluation will be needed to link traditional phenological observations with “index trees,” near-surface and satellite remote-sensing observations, “citizen science” (observations by citizens), and results published on the Internet.


Introduction
To accurately evaluate changes in the spatial and temporal variation of ecosystem functioning (e.g., evapotranspiration, photosynthesis) and services (e.g., regulating, cultural) under the rapid climate changes that are occurring due to global warming, long-term, continuous, and broad-scale phenological observations are required (Richardson et al. 2013). The spatial and temporal variations of the timing of leaf flush and leaf fall affect biogeochemical processes such as CO 2 uptake during photosynthesis and the emission of biogenic volatile organic compounds (Peñuelas et al. 2009;Polgar and Primack 2011;Keenan et al. 2014) and carbon, water, and heat budgets through the effects on evapotranspiration and physical properties such as latent and sensible heat, albedo, and the aerodynamic roughness of the land surface (Peñuelas et al. 2009;Polgar and Primack 2011;Schwartz et al. 2013). A phenological mismatch between plants and their animal pollinators and consumers may increase the risk of losses of biodiversity and supporting services (SCBD 2010;Polgar and Primack 2011;Kudo 2014). Under the future warming that is predicted to occur, the time lag between the blossoming date of species such as cherry (Prunus spp.) and the dates of associated festivals may disrupt tourism activities. To prevent the loss of cultural services such as flower tourism, decision-makers must develop adaptation plans (Sparks 2014).
Many previous studies examined the spatial and temporal variation of phenology based on fixed-point visual observations of Bindex trees^by experts (Linderholm 2006;Ogawa-Onishi and Berry 2013). Such observations have the merit of permitting long-term detailed phenological records. However, this method is labor-intensive (Richardson et al. 2007), the detection of phenological stages is subjective, the numbers of species and individuals that can be surveyed are limited, the observations fail to cover the full range of the species (they are mainly distributed in cities and low-altitude areas), and the observations may be affected by urbanization (e.g., the heat island effect; Ohashi et al. 2012). To solve these problems, we can take advantage of in situ (near-surface) and satellite remote sensing, which permit more objective, long-term, continuous phenological observations of many individuals over a wide area (Linderholm 2006;Morisette et al. 2009).
Long-term continuous phenological observations with digital cameras and spectral radiometers have been conducted since around 2003 for a range of ecosystems, mainly by American and Japanese networks. In the USA, this has been done by the PhenoCam network (http://phenocam.sr.unh.edu/ webcam/. Accessed 12 August 2015; Richardson et al. 2007); in Japan, it has been done by the Phenological Eyes Network (PEN; http://www.pheno-eye.org. Accessed 12 August 2015; Nasahara and Nagai 2015). However, these observations cover a limited area. Satellite remote sensing of vegetation indices, such as the data from the Advanced Very High Resolution Radiometer (AVHRR) sensor onboard National Oceanic and Atmospheric Administration (NOAA) satellites since 1981, can detect the spatial and temporal variation of phenology at large scales (Delbart et al. 2006;Buitenwerf et al. 2015). Furthermore, the Moderate-resolution Imaging Spectroradiometer (MODIS) sensors onboard the Terra and Aqua satellites have provided phenological data since 2000, and the VEGETATION sensor onboard the Satellite Pour l'Observation de la Terre (SPOT) have provided data since 1999 at point, regional, and continental scales, with a resolution of 500 to 1000 m (Delbart et al. 2006;Nagai et al. 2015b).
Japan includes subarctic to subtropical zones, and 70 % of the land is mountainous. The mountainous region in central Japan ranges widely in altitude (by 2000 to 3000 m) over relatively short distances (50 km). Japan's maritime climate is strongly influenced by global atmospheric patterns such as the El Niño-Southern Oscillation and the Arctic Oscillation (JMA 2014). These geographical characteristics create large regional differences and diversity and high spatial variation of vegetation (Natural Environmental Information GIS 1999). For this reason, we require not only phenological observations based on index trees but also quantitative satellite data that cover many ecosystems and plant species over a wide area, with low labor costs, to obtain and process the data.
In this review, we describe plant phenological observations in Japan based on remote sensing and discuss our current knowledge, problems, and future developments based on these studies. The geography of Japan allows for evaluation of the spatial and temporal variation of plant phenology along vertical (altitude) and horizontal (latitude) gradients over short distances (Doi and Takahashi 2008;Matsumoto 2010;Nagai et al. 2015b). Nagai et al. (2015b) found that the spatial and temporal variations of the timing of the start and end of the deciduous forest's growing season were strongly affected by both vertical and horizontal gradients. In addition, the Japanese Ministry of the Environment, Forestry Agency and other organizations have surveyed the vegetation cover and forest structure throughout the country for many years (Natural Environmental Information GIS 1999). This geographical and social background suggests the possibility of taking advantage of a large body of evidence to accurately evaluate the spatial and temporal variation of plant phenology. In this review, we will not discuss previous agricultural and horticultural studies because the effects of human manipulations may conceal the effects of climate change.

Background for Japanese phenological observations
The dramatic lifestyle changes since World War II (rapid urbanization) drastically reduced opportunities to experience nature. However, citizens continue to enjoy viewing of spring cherry blossoms (hanami; Aono and Kazui 2008), which has been recorded since the ninth century AD, and maple leaves (Acer spp.) in autumn (momiji-gari; Aono and Tani 2014), which has been recorded since the tenth century AD, because these activities are part of Japan's traditional culture. Recently, meteorological and tourist services, city offices, and tourist associations have published phenological information such as flowering and leaf coloring dates at scenic spots, resorts, temples, shrines, and parks on their Web sites. The Internet has dramatically increased the availability of such phenological information (Nagai et al. 2015a). Social networks such as Twitter and Facebook also let citizens share phenological information.
To find and experience these seasonal changes, citizens travel to scenic spots to view flowers and leaf coloring. Although technological crop cultivation (e.g., in greenhouses) and imported foods have undermined the sense of seasonality, citizens nonetheless receive nature's bounty (supporting ecosystem services) by purchasing seasonal vegetables and fruits. Seasonal events such as the timing of flowering and appearance of insects and birds have been important in Japanese literature, including classical poetry (waka) and haiku, since ancient times. Japanese citizen scientists have actively conducted phenological observations. BIkimono Log,^managed by the Biodiversity Center of Japan, is a major system for collecting and sharing biological information (http://ikilog.biodic.go.jp. Accessed 12 August 2015). These cultural, historical, and social backgrounds have created favorable conditions for Japanese phenological observations.

Phenological observations using index trees
In this section, we will briefly summarize previous phenological studies based on observations of index trees. Ogawa-Onishi and Berry (2013) provide detailed descriptions of Japanese studies. Typical data include the timing of flowering, leaf flush, leaf coloring, and leaf fall that have been published by meteorological agencies (Matsumoto et al. 2003;Doi and Takahashi 2008;Ibáñez et al. 2010;Matsumoto 2010). This dataset has been compiled at weather stations throughout Japan since 1953(JMA 1985. In addition, long-term continuous observations of these parameters have been conducted at research forests and agricultural experiment stations (Miller-Rushing et al. 2007;Fujisawa and Kobayashi 2010).
Japan stretches more than 3000 km, from 45°N to 20°N and from 123°E to 154°E, and comprises nine climatic divisions (Suzuki 1962). This allows mapping of the spatial characteristics of the timing of flowering, leaf flush, leaf coloring, and leaf fall in each season and examination of the changes in timing in response to climate change (JMA 1985). By way of example, Fig. 1 shows the spatial distribution of the first flowering date of cherry (Prunus × yedoensis; somei-yoshino) in 2014, and Table 1 provides details of these data. Global warming has advanced blooming of this cherry by 5 days from 1953 to 2013 (JMA 2014) and leaf coloring of Ginkgo biloba by 4 days from 1953 to 2000 (Matsumoto et al. 2003), but has delayed leaf coloring of Acer palmatum by 18 days from 1953 to 2013 (JMA 2014) and leaf fall of G. biloba by 8 days from 1953 to 2000 (Matsumoto et al. 2003). Unlike the case in North America and Europe, Japan's delay in autumn phenology has been clearer than the advance in spring phenology (Matsumoto et al. 2003;Ibáñez et al. 2010;Matsumoto 2010).

Japanese phenological observation networks
Recently, automatic fixed-point phenological observations using near-surface remote-sensing techniques such as digital cameras and spectral radiometers have been conducted to support research and outreach activities. Typical networks include the Internet Nature Information System (INIS) in national parks across the country, managed by the Biodiversity Center of Japan (http://www.sizenken.biodic.go.jp. Accessed 12 August 2015; Ide and Oguma 2010), PEN at ecosystem study sites (http://www.pheno-eye.org. Accessed 12 August 2015; Nasahara and Nagai 2015), and Cyberforest, which mainly supports university forest research (http://landscape. nenv.k.u-tokyo.ac.jp/cyberforest/Welcome.html. Accessed 12 August 2015; Saito et al. 2004). These networks have published daily phenological images on their Web sites. In addition, PEN has conducted automatic fixed-point spectral measurements of the canopy surface in forests and grasslands to provide vegetation index data (Nasahara and Nagai 2015).
The main aims of these systems were to introduce national parks to citizens through outreach activities (Ide and Oguma 2010), permit ground truthing of satellite remote-sensing data (Nasahara and Nagai 2015), and improve environmental education (Saito et al. 2004). INIS targets the landscape scale (200 to 1000 m), whereas PEN and Cyberforest aim for coverage from shoot to landscape scales (20 cm to 200 m; Saito et al. 2004;Ide and Oguma 2010;Nasahara and Nagai 2015). Figure 2 summarizes the locations of the sites in each phenological observation network.

Phenological observations by digital cameras and spectral radiometers Digital cameras
We can detect differences in plant phenological characteristics among ecosystems and tree species by analyzing the temporal variations in the red, green, and blue (RGB) digital numbers extracted from daily digital images (Richardson et al. 2007;Ide and Oguma 2010). These analyses have been conducted by INIS and PEN in various ecosystems (Ide and Oguma 2010;Nasahara and Nagai 2015). At some study sites, continuous daily phenological images exist for more than 10 years. Target ecosystems have been typical Japanese vegetated landscapes. Accordingly, the studies provide ground truthing data for phenological observations obtained by satellite remote sensing. However, some typical vegetation types such as oak (Quercus serrata) and mangrove (Rhizophoraceae spp.) forests have not yet been analyzed. In addition, Cyberforest has not conducted this analysis. Table 2 summarizes the previous studies that have been conducted in each Japanese ecosystem.
The year-to-year variation in the timing of flowering, leaf flush, and leaf fall for the whole canopy and for each tree species in deciduous forests can be detected by analyzing RGB time series (Ide and Oguma 2010;Nagai et al. 2011a;Inoue et al. 2014). The seasonal variation of these values correlates well with that of gross primary production (GPP) in this forest , deciduous coniferous forest (larch (Larix kaempferi) forest; Ide and Oguma 2011), and evergreen coniferous forest (Japanese cedar forest; Saitoh et al. 2012;Nagai et al. 2013a). Nagai et al. (2013a) examined differences between the seasonal variation in RGB-derived indices between a Japanese closed canopy evergreen coniferous forest (Japanese cedar, Cryptomeria japonica) and an Alaskan open canopy evergreen coniferous forest (black spruce, Picea mariana). However, insufficient comparative studies have been conducted in Japan. Sonnentag et al. (2012) compared the RGB-derived indices extracted from images taken from several types of digital cameras. They found that the seasonal patterns of RGB-derived indices were likely to be more relevant to detecting changes in plant phenology than were particular RGB-derived indices values. This means that the seasonal patterns of the changes in RGB-derived indices could be used to detect phenological transitions without having to account for the characteristics of different types of digital cameras. However, additional examinations of the unique and common characteristics of the seasonal variation of RGB-derived indices are required for Japan and elsewhere.
Very steep terrain, with large altitude changes (e.g., 900 m) over short distances (e.g., 14 km), exists in the mountainous region of central Japan (Muraoka and Koizumi 2009). Nagai et al. (2015b) reported that along a vertical gradient in deciduous forests between 35°N and 46°N, leaf flush began 1 day later per 38.6 to 68.5 m increase in altitude and leaf fall began 1 day earlier per 78.1 to 142.9 m increase. This suggests the possibility that studying a small river basin in steep terrain would let researchers evaluate the sensitivity of plant phenology to climate change using a Bspace for time substitution,^a lthough more such studies are needed to confirm this. Studies conducted over a wide range of latitudes (horizontal distance) are also required to account for the spatial distribution of tree species and populations, genetic diversity, and phenological plasticity (Doi et al. 2010). However, studies over a wide range of latitudes include more uncertainties because of the greater variation in topographic and environmental characteristics.
In forest ecosystems, digital cameras have generally been installed at the top of CO 2 flux or ecosystem observation towers (Richardson et al. 2007;Nasahara and Nagai 2015). Nagai et al. (2013b) found that the seasonal patterns of RGBderived indices extracted from upward images taken from the forest floor detected canopy phenology almost as well as downward images taken from the top of a tower in a closed canopy deciduous broad-leaved forest. Eliminating the requirement to install towers would allow such studies to be conducted at more sites for the same research budget.
Plant phenology is mainly affected by physiological events such as the onset of flowering or changes of leaf biomass and the amount of leaf pigments (Sims and Gamon 2002); Nagai et al. (2011a) found that RGB-derived indices can sensitively detect these changes based on changes of leaf color during leaf flush (from yellowish to dark green) and leaf coloring (from yellowish green to yellow or red). This seasonal variation can also be detected by spectral reflectance data and vegetation indices measured by spectral radiometers. In particular, the  green-red vegetation index (GRVI), which is based on the visible red and green reflectance bands, was better able to detect the canopy surface color changes than the normalizeddifference vegetation index (NDVI) and the enhanced vegetation index (EVI), which are based on the visible red and blue and near-infrared reflectance bands (Motohka et al. 2010;Ide and Oguma 2011;Nagai et al. 2014a).
However, automatic fixed-point observations that target the whole canopy using spectral radiometers cannot detect differences in the seasonal changes among tree species. Nakaji et al. (2011) found that the seasonal variation in phenology of each deciduous broad-leaved tree species could be detected by measurements with a three-band multispectral camera. Saitoh et al. (2012) found that the seasonal variation in whole canopy However, it would be prohibitively difficult to perform such observations for multiple tree species and individuals at multiple points because of the high cost of the instruments and the labor-intensive analysis. Further examinations of the relationship between RGB-derived indices and spectral characteristics for a range of tree species will be required to confirm that RGB-derived indices can be used as proxies for spectral reflectance and vegetation indices. In addition, examination of the relationships between seasonal variation of RGB-derived indices, leaf traits, and ecophysiological parameters for various tree species (Noda et al. 2015b) would provide useful direct evidence for analysis of GPP values (Ide and Oguma 2011). Digital cameras can capture phenological observations at shoot, canopy, and landscape scales. However, the long distance between the camera and the target vegetation that is required to cover larger areas may prevent detailed phenological observations. Nagai et al. (2014b) found that the year-to-year variation in the timing of full blooming of Prunus sargentii could be detected at a canopy scale, but detailed continuous flowering phenology could not. Inoue et al. (2014) reported that the threshold RGB-derived index values for detection of the timing of leaf flush and leaf fall could not be defined uniquely by analyzing phenological images taken with a fisheye lens. The problem resulted from the different conditions of target individuals; the center of the images captured the canopy surface of a tree, whereas the perimeter captured the side of the canopy. These results suggest that it may be necessary to install multiple digital cameras at a study site to capture enough usable images at different spatial scales. In addition, the problem resulted from a requirement to identify particular RGB-derived index values (i.e., thresholds) for determining phenological transitions. Curve fitting approaches that can detect the inflection point for RGB-derived indices during a phenological change may help to solve this problem (Ide and Oguma 2013;Henneken et al. 2014;Klosterman et al. 2014). To continuously observe detailed flowering, fruiting, leaf flush, leaf coloring, and leaf fall phenology, it will be necessary to improve phenological observations at the shoot scale, as has been done by Cyberforest (Saito et al. 2004), and comparative analyses of phenological images acquired at different spatial scales will be required.
To conserve biodiversity, the Japanese government developed and revised the national biodiversity strategy in 2012 (MEJ 2012). This strategy has actively encouraged learning by students through hands-on experience with environmental studies. Such learning is now included in the course of study developed by Japan's Ministry of Education, Culture, Sports, Science and Technology (http://www.mext.go.jp). For instance, phenological observations of plants such as morning glory (Ipomoea spp.) are a valued part of science education at Japanese primary schools. To encourage learning through environmental studies, the online phenological data published on the Web site of Cyberforest may be helpful. Phenological observations by digital cameras are also helpful because they can replace or complement labor-intensive phenological observations using index trees conducted in university research forests and exhibition gardens (Fujimoto 2007;Sato 2014; the tree phenological observation network conducted by the Research Forests Council of Japanese universities, http://www.forest.kyushu-u.ac.jp/phenology/. Accessed 12 August 2015). Ide and Oguma (2013) demonstrated the usefulness of phenological observations by digital cameras in alpine ecosystems, where accessibility is difficult but remarkable environmental changes are occurring (Kudo 2014). These phenological observations by digital cameras may encourage students to pursue their interests through both university education and lifelong environmental learning. Strengthening cooperation with citizen scientists in a project such as the BMonitoring Sites 1000 Project,^conducted by the Ministry of the Environment's Biodiversity Center of Japan (Kudo 2014; http://www.biodic.go.jp/moni1000/. Accessed 12 August 2015), may also develop environmental studies by encouraging activities related to phenological observation networks implemented using digital cameras.

Spectral radiometers
Spectral reflectance and vegetation indices measured by fixedpoint spectral radiometers can detect the characteristics of plant phenology in secondary deciduous broad-leaved forest (Motohka et al. 2010;Nagai et al. 2010a), deciduous coniferous larch forest (Motohka et al. 2010), evergreen coniferous cypress (Chamaecyparis obtusa) forest (Nakaji et al. 2008) and Japanese cedar forest , and grassland (Motohka et al. 2010). However, fewer ecosystems and study sites have been used for measurements with fixed-point spectral radiometers than for fixed-point digital cameras. Nasahara and Nagai (2015) used spectral radiometers to detect the mean tree phenology within 2 to 10 m of the device. To support ground truthing for satellite remote-sensing data is problematic because the footprint of the in situ spectral measurement is much smaller than that of the satellite data (typically ≥30 m and as great as 1000 m). NDVI, EVI, and GRVI values in deciduous forests and grasslands increase rapidly during leaf flush and then become saturated before the vegetation achieves its maximum leaf area; they then decrease gradually during the leaf coloring and leaf fall periods (Nagai et al. 2010a;Motohka et al. 2010). Nagai et al. (2010a) examined the relationship between NDVI and EVI values observed by a spectral radiometer and the corresponding phenological images in a deciduous broad-leaved forest. NDVI values of 0.6 to 0.7 and EVI values of 0.4 to 0.5 potentially detected the onset of leaf flush, but neither index could detect the timing of leaf fall. This problem resulted from differences in the characteristics of the pattern and timing of leaf coloring and leaf fall among tree species (Nagai et al. 2014a). Nagai et al. (2014a) examined the relationship between seasonal variation of spectral reflectance and vegetation indices and spatial and temporal variations of leaf litter in a deciduous broad-leaved forest. GRVI, which observes only the seasonal color variation of the canopy surface, detected the timing of leaf fall more accurately than NDVI and EVI. Motohka et al. (2010) examined the relationship between seasonal variation of spectral reflectance and vegetation indices and daily phenological images in several deciduous forests and a grassland. The first date on which GRVI was greater than zero in the spring and the first date on which GRVI was less than zero in the autumn corresponded to the start and end of the growing season, respectively. This ground truthing information indicated a general ability to detect the start and end of the growing season in deciduous forests and a grassland by thresholdbased analysis using GRVI (Motohka et al. 2010).
The seasonality of NDVI and EVI in evergreen coniferous forests was much smaller than that in deciduous forests (Nakaji et al. 2008;Nagai et al. 2012). In contrast, GRVI could detect the seasonal variation of leaf coloring at the canopy surface of a Japanese cedar forest in winter (winter reddening; Nagai et al. 2012). The seasonal variation of leaf color at the canopy surface in Japanese cedar forests was caused by changes in leaf pigments (Han et al. 2003). This seasonality correlated well with the seasonal variation in leaf photosynthetic capacity . Detecting such changes is important because evergreen coniferous forests account for 30 % of Japan's total forested area (Japan FAO Association 1997).
Many previous studies indicated that seasonal variation of satellite EVI values was correlated with that of GPP in various ecosystems (Rahman et al. 2005;Huete et al. 2006). These reports suggest that GPP could be directly evaluated by satellite EVI. To test this possibility, Nagai et al. (2010b) and Muraoka et al. (2013) examined the relationships between EVI values observed with a spectral radiometer and GPP obtained from CO 2 flux tower measurements and between various vegetation indices (including EVI) observed with a spectral radiometer and photosynthetic capacity estimated by ecophysiological measurements and ecosystem modeling in the same deciduous broad-leaved forest. GPP increased exponentially with increasing EVI in spring; in contrast, GPP showed a linear decrease with decreasing EVI in autumn.
Previous studies used 8-, 10-, and 16-day composite satellite EVI data that comprised the best-quality data, with no effects of cloud contamination or atmospheric noise within these periods (Rahman et al. 2005;Huete et al. 2006). Estimated GPP values based on the composite satellite data were closer to the potential GPP values (i.e., the GPP under ideal conditions) observed under clear skies than to estimates of the actual GPP values under clear skies or rainy and cloudy conditions (Nagai et al. 2010b). During the leaf flush and leaf fall periods, high-quality Terra/MODIS NDVI data with no effects of cloud contamination or atmospheric noise were available in Japan for only 3 to 7 days per month (Nagai et al. 2011b). For this reason, GPP values estimated from satellite Japanese EVI data may overestimate the real values.

Phenological observations by satellite remote sensing
In this section, we will briefly summarize previous phenological research based on satellite remote sensing. Our previous research found fewer phenological studies in Japan than in North America, Europe, and mainland Asia (Nagai et al. 2015b). This may be because the 8-km spatial resolution of the Global Inventory Modeling and Mapping Studies data (http://staff.glcf.umd.edu/sns/branch/htdocs.sns/data/gimms/. Accessed 12 August 2015) is too coarse to account for the heterogeneous Japanese plant phenology that results from the complex topography and vegetation cover. In contrast, the Terra and Aqua/MODIS GRVI, with 500-m spatial resolution, could detect the spatial variation in the timing of the start and end of the growing season along both vertical and horizontal gradients in Japan (Nagai et al. 2015b). Figure 3 shows the resulting spatial distribution of the timing of the start and end of the growing season in 2014.
To accurately evaluate the sensitivity of plant phenology to climate change, it is necessary to evaluate the year-to-year spatial variation of the timing of leaf flush and leaf fall over a wide area and a long time period. Although Nagai et al. (2015b) performed a decadal-scale analysis for Japan, they only showed the mean spatial and temporal timing of the start and end of the growing season during the study period for horizontal zones (2°latitude intervals) and vertical zones (200-m intervals). Previous research suggests that in a temperate climate, the sensitivity of plant phenology to climate change is clearer in southern parts of the region than in northern parts (Shen et al. 2014;Wang et al. 2015). In Japan, many river basins are divided by mountain ranges, creating small climatic zones close to each other but with different geographical characteristics. These divisions often occur at the boundaries between local authorities that are responsible for tracking phenological data. As a result, obtaining a detailed evaluation of the sensitivity of plant phenology to climate change would require the participation of Japan's 47 prefectural and city governments. Hadano et al. (2013) calibrated a phenological model for prediction of the timing of leaf flush by using both spatial variation of this timing detected from MODIS GRVI data with a 500-m spatial resolution and meteorological data at weather stations throughout Japan. They predicted the spatial and temporal variation of the timing of leaf flush under future climate change. Previous phenological models were based on in situ  Fig. 4 Summary of the relationship between the amount and spatial resolution of the available datasets for phenological observations. Data obtained by citizen scientists tends to be openly published on the Internet phenological datasets at multiple points (Xu and Chen 2013), and Hadano et al. (2013) suggested that it may be possible to accurately detect the sensitivity of plant phenology to climate change by integrating an analysis of satellite data with such a phenological model. In a temperate zone, the timing of flowering and leaf flush was estimated based on the number of chilling hours in winter and on the number of warming hours before bud burst in spring (Kim et al. 2012;Luedeling et al. 2013). Examining the relationship between the spatial variation of the timing of flowering and leaf flush predicted by phenological models and the timing of leaf flush detected by satellite observations may improve our understanding of the sensitivity of plant phenology to climate change at regional to continental scales. Nagai et al. (2015a) examined the relationship between the spatial variation of Landsat-8 GRVI values with a 30-m spatial resolution and leaf coloring information published on Web sites on 1 day during the leaf coloring period in 2014. Ground truthing information for plant phenology obtained from weather stations, research forests, agricultural experiment stations, and phenological observation sites with digital cameras is not sufficient because of the scale mismatch between the satellite and in situ data. Further examination of the relationship between satellite-detected plant phenology and phenological information published on the Internet is required to account for the footprint of the satellite data. Sekizawa et al. (2015) detected the evacuation-induced land cover changes following the Fukushima Daiichi nuclear power plant disaster in 2011 by analyzing seasonal patterns of vegetation indices derived from MODIS data with a 250-m spatial resolution. Plant phenology may change significantly in response to changes in land cover or land use, ecological succession, and invasion by nonnative species. In Japan, progressive abandonment of cultivation and succession to bamboo forests has occurred in mid-to low-elevation mountains and suburban regions (MEJ 2012). Further studies will be required to accurately evaluate the relationship between the spatial and temporal variability of satellite-observed vegetation indices and those of vegetation cover at scales ranging from individual sites to whole landscapes.

Conclusions and future developments
In the many previous studies described in this review, phenological observations were obtained using index trees for a few tree species, near-surface remote-sensing observations for a few ecosystems or tree species, and satellite remote sensing at a regional scale. The differences in these approaches resulted from the different goals of the studies and differences in the available data. Suzuki (2015) noted that many previous studies of remote-sensing data were mainly from a technical perspective and that the ecological meaning By the Phenological Eyes Network (Nasahara and Nagai 2015) c By the MODIS spectral radiometer sensor aboard the Terra satellite of the signals was not sufficiently considered. In addition, many previous studies only used the limited datasets observed by scientists and other experts rather than the larger amount of data provided by citizen scientists. This may have prevented researchers from obtaining phenological data at a synoptic scale that covered multiple species over wide areas.
To solve these problems, researchers should attempt a multidisciplinary analysis that links index trees, near-surface remote sensing, satellite remote sensing, citizen science, and open-access information published on the Internet. Figure 4 illustrates the relationship between the amounts and spatial resolutions of the available phenological datasets. Japanese citizens are very interested in plant phenology, and their interest provides data at a small personal scale that can be linked with the abundant data at wider regional and national scales. Establishing this link will require the development of databases of phenological observations that will reveal ecosystem changes from the leaf scale to national scales (Muraoka and Koizumi 2009). Each source of data has advantages and problems (Table 3), but because these differ among the data sources, the advantages of one dataset may compensate for another dataset's limitations.
Two studies from outside Japan provide examples of successful multidisciplinary studies that linked remote-sensing observations with index tree observations or citizen science. Delbart et al. (2015) examined the relationship between the start of the growing season in Canada, detected by SPOT/VEGETATION vegetation index data, and the timing of flowering and leaf flush reported by citizen scientists. The phenological information observed by the citizen scientists provided useful ground truthing for the satellite data. Kim et al. (2012) examined the relationships among the timing of cherry blooming at multiple locations in Korea, the start of the growing period detected from a MODIS vegetation index, and a model for prediction of the timing of leaf flush. Phenological information from index tree observations provided useful data for a multidisciplinary phenological study in both cases.
Phenological studies tend to be local because of the specific and diverse goals of the researchers. Unfortunately, this also means that there have been no direct comparisons between Japanese plant phenology and that in other countries. On the other hand, many Japanese studies were conducted both before and after World War II. However, those studies were published in Japanese. As a result, insufficient Japanese phenological information has been made available to the international research community (Ogawa-Onishi and Berry 2013). Recently, huge amounts of phenological information have become available on the Internet in Japanese, which means that it remains unavailable to researchers who cannot read Japanese. By way of example, Table 1 summarizes the first flowering date of cherry near 58 weather stations in Japan published by the Japan Meteorological Agency. To solve these issues, we propose the following challenging but important tasks: (1) We should conduct more multidisciplinary phenological studies in Japan through collaboration among researchers with different goals and between researchers and citizen scientists, with the aim of improving our understanding of phenology and permitting comparisons among observation networks.
(2) We should evaluate both the unique characteristics of Japanese phenology and the characteristics shared with vegetation in other regions through international collaborations and through comparisons with results from foreign observation networks. (3) We should publicize Japanese language phenological data in the international research literature so that the data will be more accessible to the international research community. Of course, to accomplish this it will be necessary to obtain the necessary legal permissions to share this data, and it will also be necessary to find volunteers willing to translate the information into English and other languages. By implementing these recommendations, Japanese phenological studies will become increasingly powerful.