Distribution , habitat , and conservation status of the near-threatened Japanese Paradise-Flycatcher ( Terpsiphone atrocaudata periophthalmica ) on Lanyu , Taiwan

Distribution, habitat, and conservation status of the near-threatened Japanese Paradise-Flycatcher (Terpsiphone atrocaudata periophthalmica) on Lanyu, Taiwan Thorsten Späth , Mei-Ling Bai , Lucia L. Severinghaus 3 and Bruno Andreas Walther 4 Georg-August-University-Göttingen, Conservation Biology/Workgroup on Endangered Species, Göttingen, Germany, Formosa Natural History Information Ltd., Taipei, Taiwan, Biodiversity Research Center, Academia Sinica, Taipei, Taiwan, Master Program in Global Health and Development, College of Public Health, Taipei Medical University, Taipei, Taiwan


INTRODUCTION
Island bird species and populations are especially threatened by extinction (Johnson and Stattersfield 1990), and habitat destruction and invasive species are the two main causes (Veitch andClout 2002, Reaser et al. 2007).Therefore, there is a growing focus on the value of conserving island populations as independent conservation units (Kier et al. 2009, Pruett et al. 2017).Because conservation and management of species is dependent on knowledge of a species' distribution, ecology, and population size (Sutherland et al. 2004), an understanding of the main factors determining a species' presence will assist conservation managers in prioritizing actions and increasing the efficacy of any applied strategy (Sutherland et al. 2004).The East Asian region has a large number of island bird species that remain little known and understudied (Collar et al. 2001, Brazil 2009, Ando et al. 2014).
The near-threatened conservation status of the Japanese Paradise-Flycatcher appears to be based mostly on older Japanese studies (Appendix 1) and an assumed habitat loss on the wintering grounds (BirdLife International 2015).However, BirdLife International (2015) proposed that more studies and careful monitoring are needed in order to (1) determine its current distribution and abundance on both the breeding and wintering range and (2) increase knowledge of its habitat requirements.To further improve our knowledge of this species, we (1) assessed the species' current conservation status using the available literature, (2) used the species' ecological preferences to determine the geographic distribution on Lanyu, (3) estimated Lanyu's population size, and (4) estimated the habitat change on Lanyu over the last half century.

Conservation status assessment
To assess the current conservation status of the Japanese Paradise-Flycatcher, we gathered information from all recent publications and reports about the species' status in Japan, Korea, Taiwan, and the Philippines.Publications were found by (1) searching Google, Google Scholar, and Web of Science using the appropriate search terms (e.g., English and Latin species name, country names, etc.); (2) checking standard references (e.g., Brazil 1991, Brazil 2009, Coates et al. 2006, Severinghaus et al. 2010, 2017) and websites (e.g., BirdLife International 2015, HBW Alive 2016) for information and further references; (3) emailing any researcher whose email was given in these references and asking for further information and publications (any researcher who responded is mentioned in the Acknowledgements).

Study site
Lanyu (21°58′-22°06′N, 121°29′-121°35′E), Taitung County, Taiwan, is located 64 km east off the southern tip of Taiwan.It is a 46.3 km² mountainous volcanic island of mainly wooded habitats, generally shallow soil, and mostly steep topography.The tallest mountain is Mt.Hongtoushan at 552 m a. s. l., and nearly 10 peaks are > 400 m (Fig. A1.1).The island is bounded by narrow strips of coastal flatlands (Fig. A1.4), and seven villages are located within these strips.The mean annual temperature is 22.7º C, and the mean annual rainfall is approximately 3000 mm with an average of about 220 rainy days (Taiwan Central Weather Bureau).Because of year-round heavy winds and frequent typhoons (usually several per year), the slopes are largely covered by dense shrubs and stunted trees while the occurrence of taller forests is confined to wind-sheltered plains and valleys (Chao et al. 2010).Although the most protected forests can reach a height of 20 meters, the forest canopy on slopes is usually around eight meters and, depending on wind exposure, even lower.Much of the taller forests of the coastal flatlands have been transformed into agricultural fields or frequently burned grasslands, while other accessible forests are structurally altered for the purpose of fruit production or construction wood (T.S., personal observation).

Territory size estimation
To estimate the territory size of male Japanese Paradise-Flycatchers, we selected five study plots that represented the various forest types found on the island and that were sufficiently accessible for frequent and intensive observations (some steep slopes are virtually inaccessible).From early March to late September 2009 and 2010, we censused each plot for the flycatcher's presence using conspecific playback at control points if we had not made spontaneous observations within the first few minutes after reaching the plot.The majority of these territories were occupied from late March, and no individuals left the island before September.Our observations also indicated that territorial borders became less distinct and that vocal activity or response to playback decreased as the breeding season progressed.We therefore chose April and May in 2010 as the period of highest territorial and vocal activity during which to carry out our islandwide census.
In 2009, we also captured flycatchers using mist nets at the same five study plots and color-banded each individual.We could furthermore distinguish uncaptured males in these study plots by differences in body characteristics, such as the length of tail feathers and aspects of plumage coloration on the throat, breast belly, mantle, greater coverts, and tail feathers.We then mapped out individual territories within each study plot between late March and late May.Once a week, each study plot was searched in the morning; points of encounters were recorded on a map based on high-resolution satellite images, and the bird was followed until we lost sight of it.These points were then used to plot the trajectory of the bird.The outermost points were later connected as a minimum convex polygon using Hawth's Tools for ArcGIS.The mean and standard deviation of 32 mapped territories was 1.16 ± 0.43 hectares (some mentioned in Severinghaus and Bai 2009, unpublished data cited in Appendix 1).

Environmental data layers
We assembled 22 environmental layers with a resolution of 100 m (or 0.1 x 0.1 km or 1 hectare grid cell) for Lanyu (Table A1.1).The two normalized differenced vegetation index (NDVI) layers were generated from a satellite image taken by FORMOSAT-2 in April 2007 (Center for Space and Remote Sensing Research; http://www.csrsr.ncu.edu.tw/) with an 8 m multispectral resolution.We calculated the NDVI for each 8 x 8 m pixel and summarized the mean and standard variation of the NDVI values in each 1 hectare grid cell (Fig. A1.2).
A digital terrain model (DTM) with 40 m resolution generated by the Aerial Survey Office (http://www.afasi.gov.tw/) was used to calculate four topographical data layers, namely elevation, aspect, slope, and solar irradiation for each 1 hectare grid cell.
These four data layers were then used to generate the data layers 3-13 (Table A1.1).
Based on the manual interpretation of aerial photographs taken in 1948 by the U.S. Airforce and in 2006 taken by the Aerial Survey Office, we created two layers of land cover types using ArcGIS that included the following categories: (1) roads; (2) bare coastal land (beach, rock); (3) other built-over land (buildings, parking spaces, etc.); (4) farmland or grassland; (5) hilltop shrub; (6) stunted forest; (7) tall forest (either mixed mature forests or forests dominated by mature Pometia pinnata; Fig. A1.3), (8); other types, e.g., landslides, inland water, etc. Stunted forest contained both low secondary forests and primary shrubby woodland on wind exposed slopes because these two types could not be differentiated in the aerial photographs.These two data layers were used to generate the data layers 14-22 (Table A1.1).To characterize the land cover types of each grid cell, we calculated the percentage value of each type for each cell (data layers 14-21 in Table A1.1), and finally used the same data to generate a categorical data layer, which corresponded to the dominant land cover type of each grid cell (data layer 22 in Table A1.1).Only this data layer 22 of the dominant land cover types was generated for both 1948 and 2006 in order to assess habitat change.

Modeling species distribution
Because we had established that male territories were about one hectare, we used ArcGIS to establish 1 hectare (100 m x 100 m) grid cells as sampling units across the entire island.Any grid cells that were not entirely covered by land surface were a priori excluded.We then used Hawth's Tools for ArcGIS to randomly select 224 grid cells among the remaining grid cells, which represent about 5% of Lanyu's total land surface.Because of our limited resources and the island's rugged terrain that renders many parts inaccessible, we a priori chose this 5% threshold, which was the maximum number of grid cells that we could reasonably survey during the period of the flycatcher's highest territorial activity, i.e., April and May.
Each grid cell was visited once between sunrise and noon during April and May 2010.Once the center of the grid cell had been reached, we waited for spontaneous calls or songs for up to three minutes; if none were detected, we used playback consisting of recordings of calls and songs of the local population and played them up to three times at three minute intervals to elicit a response.The species was defined to be present if it was detected visually or acoustically during this time, and was defined to be absent from the respective grid cell if it could not be heard or seen during this one visit.This protocol for detecting presence was based on our preliminary field work.During our regular controls of the five study plots but also when checking additional transect lines, we used playback at fixed control points to check for the presence or absence of the territorial male if it could not be detected by visual or acoustic detection.We found that there was no significant effect of time of day on the playback response, and that we never failed to detect the territorial male even if the playback was conducted only in the grid cell's center.Therefore, we are confident that our method of detection was reliable.
However, we are aware that absence is much harder to establish than presence, and that some of the absence records in Fig. 2 may be false absences.For this reason, we chose a modeling software that only needs presence records, namely Maxent.Therefore, our absence records were not used further in the modeling analyses described below.We avoided multicollinearity and overfitting in our Maxent distribution models by reducing our original set of 21 continuous variables, i.e. environmental data layers, to eight variables (Table A1.1).We assessed collinearity by constructing a correlation matrix for the 21 variables based on Spearman's correlation coefficient (r s ), and we removed all but one variable if r s was > 0.80 between two or more variables.In each case, we kept only the variable with the highest percent contributions across the 48 global Maxent models which we initially ran.
With this reduced set of eight continuous variables and one categorical variable (dominant land cover type), we generated suitability distributions for the Japanese Paradise-Flycatcher across Lanyu using the maximum entropy algorithm implemented in the software Maxent (Phillips et al. 2006, Phillips and Dudík 2008, Elith et al. 2011).We used Maxent with the logistic output of probabilities.Maxent finds the probability distribution of maximum entropy subject to constraints imposed by the information available from the observed distribution of the species and environmental conditions across the study area.Thereby, Maxent transforms environmental variables into feature vectors and then uses entropy as the means to generalize specific observations of the species' presence; therefore, it does not require absence points within its theoretical framework.
Recommended default values for the Maxent modeling procedure were used for the convergence threshold (10-5), maximum number of iterations (500), data being randomly divided into 70% training and 30% testing data, and cross-validation using jackknife resampling.The selection of "features" (environmental variables or functions thereof) was also carried out automatically, following default rules dependent on the number of presence records.
In order to obtain the most suitable model for our data, we generated 48 models by varying the following modeling features: 1.The number of background points was varied at four settings: 500, 1000, 2500, and 5000.
In this way, we generated 4 x 4 x 3 = 48 models for comparison.First, we generated 48 models using all 22 environmental variables, but only to score the heuristic estimates of each variable's relative model contribution and average each variable's mean contribution over the entire 48 models.This average percent contribution was then used to eliminate highly correlating variables.We then generated another 48 models using the same combination of modeling features but with only the remaining nine variables (eight continuous variables plus one categorical variable).
We then evaluated the support for each of these 48 competing models using the likelihood-based methods based on the information theoretic approach proposed in Burnham and Anderson (2002).We ranked models based on the Akaike's Information Criterion adjusted for small sample size (AICc), and we accepted only those models with an AICc difference (ΔAICc) < 2 as having "substantial" support (sensu page 70 in Burnham and Anderson 2002; see also Majić et al. 2011, Grabowska-Zhang et al. 2012, Hong et al. 2016).These are also the models with the largest Akaike weights (Wi) and the smallest evidence ratios (Burnham and Anderson 2002).
To choose a final model, the best model with the lowest ΔAICc = zero can be used (e.g., Botero-Delgadillo et al. 2015a,b).However, we chose to adopt the strategy of Burnham and Anderson ( 2002) whereby all models with substantial support should be considered as having value.Therefore, we generated an ensemble model (sensu Araújo and New 2007) as our final model by combining all the models with substantial support (ΔAICc < 2).Although there are also several ways of combining models (M.Araújo, 2012, personal communication), we simply calculated the mean of all the models that were included into the final ensemble model.
Maxent also generates response curves, which show how each environmental variable affects the prediction.These response curves reflect the dependence of predicted suitability both on the selected variable and on dependencies induced by correlations between the selected variable and other variables (Phillips 2017).We examined these response curves to check model fit and how the environmental variable affects the prediction.
As recommended by Botero-Delgadillo et al. (2015a,b), model significance was also evaluated with threshold-dependent binomial probability tests applied in each replicate, and model performance was evaluated with the mean values and standard

Population size estimate
To estimate population size, we divided the size of the entire area deemed suitable in the binary distribution map generated using the minimum training presence threshold by the mean estimated size of a male's territory, namely 1.16 ± 0.43 hectares.It should be noted that this estimation assumes that every territory is occupied.

Habitat change assessment
Our results below demonstrate that the Japanese Paradise-Flycatcher occurred almost exclusively in hilltop shrub, stunted forest, and tall forest.For our historical comparison of habitat change, we therefore lumped these three land cover types into one type called "forest."We then compared the forest area of 1948 to the forest area of 2006 using data layer 22 (Table A1.1) to assess the amount of change in the area of habitat that is potentially suitable for the Japanese Paradise-Flycatcher.

Conservation status assessment
A brief assessment of the current conservation status of the Japanese Paradise-Flycatcher in Japan, Korea, Taiwan, and the Philippines is given in Appendix 1. From the available but relatively sparse information, it appears that the Japanese population is now stable after a decline between the 1970s and 1990s, that the Korean population is probably stable or slightly increasing (although some experts contended that the species has been decreasing), the Taiwanese population is stable, while the status of the Philippine population is unknown although it was certainly present in 2006 and 2007.

Modeling species distribution
The field work in 2010 established 120 presence and 104 absence grid cells (Fig. 2).Almost all presence grid cells were at some distance from the coastline, i.e., almost no singing males were found in the coastal lowlands.Away from coastal areas, the presence grid cells were distributed relatively evenly across the entire island, from north to south and east to west.The means, variations, and ranges of the environmental values within the presence and absence grid cells are given in Table A1.2.
The presence grid cells covered an elevational range from 16-399 m with a mean of 172 m (Table A1.2, Fig. 3A).The area of tall forests within the presence grid cells ranged from 0-100%, with a mean percentage of 26% (Fig. 3B).Furthermore, 115 out of the 120 presence grid cells (or 96%) ranged from 0.68-0.77for mean NDVI, with a mean of 0.74 (Fig. 3C).The dominant land cover types for the presence grid cells were farm/grassland (5 grid cells), hilltop shrub (6 grid cells), stunted forest (77 grid cells), and tall forest (32 grid cells).Therefore, grid cells containing either mostly stunted forest (Fig. A1.5) or mostly tall forest (Fig. 3B) represented 91% of all the presence grid cells.Furthermore, grid cells containing mostly one of the three "forest" types (namely, hilltop scrub, stunted forest, or tall forest) represented 96% of all the presence grid cells.To summarize, the Japanese Paradise-Flycatcher was detected mostly at midelevations of 50-300 m, preferred the areas of Lanyu with the highest NDVI values, and was found almost exclusively in forest, and predominantly in stunted or tall forest.
We ran Maxent with the nine selected environmental variables and used them to construct 48 models.Of these 48 models, five models had a ΔAICc < 2 (Table A1.3).
For these five "best" models, the mean regularized training gains were 0.49 ± 0.04 (n = 5 * 20 replicates = 100 replicates, same below), the mean values of the threshold-based training omission  error rates were 0.07 ± 0.02 (n = 100), and the threshold-based test omission error rates were 0.17 ± 0.05 (n = 100), thereby confirming that these models performed reasonably well.All cross-validation replicates of these models were statistically significant according to the threshold-dependent binomial probability test (all P < 0.05).
Among these 48 models, the three variables that were almost consistently chosen as the variables with the highest heuristic estimates of their relative model contribution were mean elevation, percentage of tall forest, and mean NDVI, having a combined contribution of 54.4% (Table 1).Although these three variables did not yield the highest training gain achieved with all but the regarded variable, these three variables achieved the highest training gain achieved with only the regarded variable, with one exception, namely dominant land cover type (Table 1).
We combined these five models into an ensemble model (Fig. 4).Using the minimum training presence threshold, we then converted this ensemble model into a binary map of the estimated ecological niche of the Japanese Paradise-Flycatcher on Lanyu Island (Fig. 5).This binary map covers 1203 out of a total of 4632 grid cells, or 26.0% of the island's surface (equaling 12.03 km²).
After modeling environmental suitability as a function of these nine variables, the individual response curves generated by Maxent for the three most important environmental variables indicated that environmental suitability for the Japanese Paradise-Flycatcher was highest at midelevations of about 30-350 m (Fig. 6A) and increased with the percentage of tall forest within the grid cell (Fig. 6B) and with NDVI (Fig. 6C).These results generated by Maxent therefore correspond reasonably well with the results from the presence grid cells (Fig. 3).

Population size estimate
Because the average size of a male's territory was 1.16 ± 0.43 hectares, we extrapolated that the area deemed suitable in the binary distribution map (namely 12.03 km²) could hold 1037.1 territories.If we use the lower and upper estimates given by the standard deviation (0.73 and 1.59, respectively), we can calculate an upper and lower limit of the number of male territories of 1647.9 territories and 756.6 territories, respectively.

Habitat change assessment
Our historical comparison of the availability of "forest" revealed that very small areas of forest were lost between 1948 and 2006 (specifically, 0.66 km² or 1.43% of the island's area).Almost all of these areas are located near the coast (Fig. A1.6).During the same time period, relatively large areas of forest were gained (specifically, 8.23 km² or 17.77% of the island's area).These areas are distributed across the entire island, but especially in the areas just south of the island's center (Fig. A1.7).The net gain of forest areas was therefore 7.57 km² (or 16.35% of the island's area).Some of these gained forest areas are located within our binary presence model (Fig. A1.8), and this gained area covers 3.76 km² (or 8.13% of the island's area).Assuming that the Japanese Paradise-Flycatcher's habitat preference for the three "forest" categories did not change between 1948 and 2006, and that other important variables, e.g., territory size, also did not change, this area of potentially suitable habitat gained from 1948 to 2006 could support an additional 324.3 male territories, or 31.27% of the entire 1037.1 territories.

DISCUSSION
Our field work indicates that males of the periophthalmica subspecies of the Japanese Paradise-Flycatcher found on Lanyu Island prefer to establish territories in midelevational forest habitats that include both stunted and tall forests and that have relatively high NDVI values.Using the 120 presence locations from our field work, we produced an ensemble model of the flycatcher's ecological niche that estimates that such suitable habitat extends over about 12.0 km² or about 26% of Lanyu's area.Because our field work also established that the average size of a male's territory is 1.16 ± 0.43 hectares, we could extrapolate that approximately 1037 territories may exist on Lanyu Island, with upper and lower limits of approximately 1648 territories and 757 territories, respectively.However, we note that any such population size estimate is associated with various uncertainties, such as that not all potential territories may be occupied, that there may be annual variations in population size, and that the use of different thresholds (Nenzén and Araújo 2011) than the one we used (minimum training presence threshold) would have resulted in different estimates of the area of suitable habitat.Nevertheless, even with these uncertainties, our estimate is a substantial increase over previous estimates of < 100 breeding pairs (Brazil 2009) and < 500 individuals (Fang 2005).
Our results can also be used to design conservation measures for the Lanyu population because our results suggest that the Japanese Paradise-Flycatcher prefers relatively wet midelevational forest habitats.Currently, the Lanyu population is probably not threatened in its preferred habitat because limited land use conversion has been occurring in Lanyu, with an overall net gain of 7.57 km² of forest habitats between 1948 and 2006, which, all other things being equal, should have resulted in a population increase of around 30%.Despite the fact that the local population rejected a plan to turn parts of the island into a national park in the 1980s and 1990s (Huang 1997), there are currently no socioeconomic factors that would drive forest destruction or conversion.Therefore, continuous monitoring and maintenance of this suitable forest habitat should ensure the long-time survival of this species, assuming no stochastic catastrophic events, such as typhoons, or long-term changes, such as climate change.
Consequently, the question arises whether the Japanese Paradise-Flycatcher should remain in the "near-threatened" category.The main reason given by BirdLife International ( 2015) is that declines were noted in parts of Japan's breeding range that were presumably caused mostly by habitat loss and degradation within the wintering grounds.However, as outlined in Appendix 1, this assessment is based mostly on relatively old studies and only from parts of Japan.Since then, its status in Japan appears unchanged since the 1990s, including the Ryukyu population (H.Higuchi, 2015, personal communication).Our study suggests that the Lanyu population is also relatively safe, and despite some yearly fluctuations (L.L. S., unpublished data), it may even have increased over the last half century because of habitat expansion, although there are no historical data on population numbers to confirm or reject this supposition.
Finally, the Korean population may have been increasing, which could be due to a northward expansion due to climate change (e.g., Kwon et al. 2014, Wu andShi 2016) and successful reforestation, although other experts considered this to be an artificial effect derived from increased sampling effort or even contended that the species has been decreasing (Appendix 1).If the Philippine population is also found to be stable (pending further study), a case could be made for down-listing the species to "least concern."However, there remains an urgent need to find out more about this species' migration routes and wintering grounds because habitat loss and degradation in the wintering grounds were suspected to have caused the original decline of the Japanese populations (BirdLife International 2015).With rainforest habitats continuing to decline rapidly across Southeast Asia (Wilcove et al. 2013, Walther et al. 2016), it is likely that the Japanese Paradise-Flycatcher is being impacted by these environmental changes.Given rapid environmental changes, the continuous reassessment of the conservation status of Taiwan's birds (Walther et al. 2011, Wu et al. 2014, Lin et al. 2016) and East Asia's birds (Collar et al. 2001, Kirby et al. 2008) must remain a research priority.
Consequently, we have four main recommendations concerning the T. a. periophthalmica subspecies: 1.Because there is no reliable population estimate for the Philippine population (Appendix 1), a study of the size of the Philippine breeding population should be a conservation priority.

2.
Further monitoring and protection of Taiwan's breeding population should be a conservation priority for Taiwanese researchers and authorities.

3.
A study of the species' migration and wintering grounds is overdue, including changes and threats to suitable habitats in these grounds.

4.
A study to determine if this subspecies actually deserves full species status (e.g., McGregor 1907, Alcasid 1965) should be a priority for taxonomists.
Indeed, further taxonomic and genetic studies are warranted for all the subspecies (or even all of its independent conservation units) of the Japanese Paradise-Flycatcher because some or all of its subspecies may be sufficiently distinct to merit full species status.If, for example, T. a. periophthalmica was to be elevated to species status, its conservation status might be "nearthreatened" or even "vulnerable" (IUCN 2012) because of its existence in only five known locations (Lanyu and the four Batanes islands), small geographic range (~12 km² on Lanyu and unknown in Batanes), and its relatively small population size (~1000 singing males on Lanyu and unknown in Batanes).Under the IUCN criteria, species with either small populations (< 10,000 individuals) or small areas of occupancy (< 2000 km²) may be classified as vulnerable, regardless of the trajectory of their populations.T. a. periophthalmica would fulfil both of these criteria.
Even if T. a. periophthalmica was not elevated to species status, these two island populations could justifiably be considered independent conservation units, especially given that they are the only populations within either Taiwan or the Philippines.Island populations as independent conservation units have been proposed for other island bird populations (e.g., Dudaniec et al. 2011, Garcia-del-Rey et al. 2013, Ando et al. 2014, Forcina et al. 2014, Pruett et al. 2017).Conservation units have also been called evolutionarily significant units or management units (cf.Moritz 1994, Rayner et al. 2010).However, these studies have routinely included genetic analyses, which were not part of our study, again emphasizing the need for more genetic and taxonomic studies of this species.
The Lanyu and Batanes islands are regularly subjected to devastating super typhoons (e.g., Fritz 2016), with no information about the impact of these typhoons on local bird populations (although see Hong et al. 2016 for an example from Taiwan's mainland).This lack of information about the recent fate of these island populations further underlines the need for more continuous field work, especially in Batanes, but also in Lanyu.
We therefore recommend that the governments of Taiwan and the Philippines should support future research on bird populations on islands that could be designated as independent conservation units.We further recommend that more continuous monitoring of such populations is financed, because the only continuous field work on the Lanyu population is the one reported in this study which lasted for only two breeding seasons.Although the recently established Taiwan Breeding Bird Survey (Ko et al. 2015) is a step in the right direction, more long-term monitoring targeted at specific species and populations (e.g., Lin et al. 2007) is required. Responses

LITERATURE CITED
In Korea, the Ministry of the Environment designated the Japanese Paradise-Flycatcher as an Endangered Bird Species Level II and reported records up to 37°N between 1997and 2005(Kim et al. 2010a).The distribution has been moving northwards with several recently confirmed records over 37°N (H. Kim, in litt. 2016).A recent status report published by Birds Korea (2014) categorized the species as an amber species which means it has (1) a global status of near-threatened, and (2) its national status assessed by NIBR (2011,2012) is "vulnerable."This report also stated that the species "is increasingly recorded in the ROK [= South Korea], with recent breeding confirmed north to at least 35°N," and both its historical and recent population trends are increasing.The breeding population was estimated to be between 100 and 999 individuals.However, Kim et al. (2010c)  Within the administrative boundaries of Taiwan, T. a. periophtalmica is known to breed exclusively on Lanyu, and its conservation status is given as category 2 "rare and valuable" (Severinghaus et al. 2010, 2017, WCAT 2013, Lin et al. 2016).Fang (2005) estimated the population as fewer than 500 individuals and warned that clearance for agriculture, roads and construction is reducing suitable habitat   A1.2) are depicted in two colours from lowest (blue) to highest (red) in 10 equal intervals, and the 100m x 100m grid cells are bounded by grey lines.A1.2) are depicted in two colours from lowest (blue) to highest (red) in 10 equal intervals, and the 100m x 100m grid cells are bounded by grey lines. .The category 'other types' is not depicted because it covered a very small area (Table A1.2).Grey areas are non-forest areas (i.e., not categorized as either hilltop shrub, stunted forest, or tall forest).

Fig. 3 .
Fig. 3. Frequency distributions depicting the variation in (A) mean elevation, (B) percentage of tall forest, and (C) mean normalized differenced vegetation index (NDVI) associated with the 120 presence grid cells for the Japanese Paradise-Flycatcher (Terpsiphone atrocaudata).Note that the bimodal distribution in Fig. 3B results from the species' preference for stunted forest (left peak; see also Fig. A1.2) and tall forest (right peak).

Fig. 4 .
Fig. 4. The ensemble distribution model of the Japanese Paradise-Flycatcher (Terpsiphone atrocaudata) across Lanyu, Taiwan, shown as a probability surface (see text for details).Probability values are depicted in two colors from lowest (blue) to highest (red) in 10 equal intervals.

Fig. 5 .
Fig. 5.The binary distribution model of Japanese Paradise-Flycatcher (Terpsiphone atrocaudata) across Lanyu, Taiwan, using the minimum training presence threshold to turn the probability surface depicted in Fig. 4 into a binary map with modeled presence (black, n = 1203) and absence (grey, n = 3429) grid cells.

Fig. 6 .
Fig. 6.The probability of presence of the Japanese Paradise-Flycatcher (Terpsiphone atrocaudata) in relation to spatial variation in (A) mean elevation (in m); (B) percentage of tall forest; and (C) mean normalized differenced vegetation index (NDVI; unitless) in Lanyu Island, Taiwan.Each individual response curve was generated by Maxent, with the curve representing the Maxent model created by using only the respective variable.Here, we used the curves generated for the best model (ΔAICc = 0).The vertical axis represents the probability of presence from the logistic output of Maxent model.The red line shows the average of the 100 replicate runs, while the blue lines shows ± 1 standard deviation.

Fig. A1. 2 .
Fig. A1.2.Gridded maps of the important environmental data layers of Lanyu used in the distribution modelling of the Japanese Paradise-Flycatcher.This map depicts NDVI mean; the values (cf.Table

Fig. A1. 3 .
Fig. A1.3.Gridded maps of the important environmental data layers of Lanyu used in the distribution modelling of the Japanese Paradise-Flycatcher.This map depicts % tall forest; the values (cf.Table

Fig. A1. 4 .
Fig. A1.4.Gridded maps of the important environmental data layers of Lanyu used in the distribution modelling of the Japanese Paradise-Flycatcher.This map depicts dominant land cover type: the colours refer to roads (black), bare coastal lands (greenish brown), other built-over land (yellow), farmland or grassland (grey), hilltop shrub (greenish yellow), stunted forest (orange beige), and tall forest (light beige).The category 'other types' is not depicted because it covered a very small area (TableA1.2).

Fig. A1. 5 .
Fig. A1.5.Frequency distribution depicting the variation in percentage of stunted forest associated with the 120 presence grid cells for the Japanese Paradise-Flycatcher.Note that the right peak in this Figure largely corresponds to the left peak in Fig. 3B.

Fig. A1. 8 .
Fig. A1.8.Binary distribution model of Japanese Paradise-Flycatcher as depicted in Fig. 5 of main text, but here, the predicted presence grid cells are overlaid with the forest areas gained between 1948 to 2006 (yellow areas in Fig. A1.7).

Table 1 .
Heuristic estimates of the relative contribution and relative performance of nine environmental variables to modeling the distribution of the Japanese Paradise-Flycatcher (Terpsiphone atrocaudata) on Lanyu Island, Taiwan.The variables are ordered by their mean percent contribution.NDVI = normalized differenced vegetation index.Values represent the model's training gain achieved with all but the regarded variable ± 1 standard deviation.§ Values represent the model's training gain achieved with only the regarded variable and excluding the remaining eight variables ± 1 standard deviation. ‡ to this article can be read online at: http://www.ace-eco.org/issues/responses.php/1167overall logistics, material, and additional support by field assistants was only possible because of substantial funding from Forestry Bureau of Council of Agriculture of Executive Yuan and Academia Sinica.
(Kim et al. 2010becies "has rapidly declined due to habitat loss (i.e.deforestation and industrialization)."Thisapparentcontradiction in the assessment of the species' trend is likely due to several factors: (1) increased sampling effort in both intensity and geographic coverage because of (i) increased survey efforts (C.-Y.Choi in litt.2016) and (ii) the "increased number of birders and bird photographers during [the last] 20 years [which] may be A 2010 survey of Jeju Island, which is a 1,848 km 2 island of the southern tip of Korea, found a total of 124 individuals in evergreen forest between 83 and 1106 meters a. s. l.(Kim et al. 2010b).Jeju's broad-leaved forests have increased in area and quality, and many parts of them are protected within a national park and thus serve as a stronghold of the species inKorea (C.-Y.Choi, in litt.2016).
over 100 times more than before" (H.Kim, in litt.2016);(2) successful reforestation after the Korean War has steadily increased suitable forest habitat in South Korea while North Korea continues to suffer from deforestation (C.-Y.Choi, H. Kim, in litt.2016); (3) however, locally, detrimental pressures may negatively affect the species.In conclusion, the general impression is that the species trend is stable or slightly increasing while spreading northwards; however, this impression is somewhat clouded by the increased observer density and the lack of reliable long-term data on a national basis.C.-Y. Choi (in litt.2016) even asserted that "any nationwide assessment on the trend of the species is unreliable and probably biased at this stage."

Table A1 .
3. Model set for modeling the ecological niche of the Japanese Paradise-Flycatcher on Lanyu Island, Taiwan, containing different combinations of the number of background points (BPs), coefficient