{'fraction_non_alphanumeric': 0.031707769330613755, 'fraction_numerical': 0.025254960133506397, 'mean_word_length': 4.749680170575693, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Forest edge creates small-scale variation in reproductive rate of sika deer

Habitat edges are considered to have an important role in determining the abundance of deer in forest landscapes, but to our knowledge there are few lines of evidence indicating that forest edge enhances the vital rate of deer. We examined pregnancy of female sika deer in Boso peninsula, central Japan, and explored how forest edges, food availability in forests, and local population density influence the pregnancy rate of sika deer. Local deer density was estimated by the number of fecal pellets, and food availability in forests was estimated by combining GIS data of vegetation distribution and the relationship between vegetation biomass and local deer density. Forest edge length was also determined by GIS data. Model selection was performed with multiple logistic regression analyses using the AIC to find the best model for accounting for the observed variation in pregnancy rates of the deer. Multiple logistic regression analysis showed that the length of forest edge had a positive effect on the pregnancy rate of females, whereas food availability in forests and local deer density had little effect. This forest edge effect was detected in a 100–200-m radius from deer captured locations, indicating that deer pregnancy is primarily determined by habitat quality within a 10-ha area. This result was confirmed by tracking females with GPS telemetry, which found that the core areas of the home range were less than 12 ha. The positive effect of edges and the lack of density dependence could be a result of high plant productivity in open environments that produces forages not depleted by high deer densities. Our results support the view that land management is the cause of the current problem of deer overabundance.


Introduction
Recent studies have shown that deer populations have spatial structures comprising units that differ in densities and vital rates (Coulson et al. 1997;Focardi et al. 2002;Pettorelli et al. 2003;Zannese et al. 2006), most of which are formed by heterogeneous landscape structure. One of the important elements causing forest heterogeneity is the forest edges that abut clearcuts, farmlands, and roads, which provides abundant and high-quality foliage (e.g., Lyon and Jensen 1980;Alverson et al. 1988;Decalesta 1997). Several lines of evidence suggest that edges are high-quality habitats for deer; local population densities were higher (e.g., Wahlstrom and Kjellander 1995;Hemami et al. 2005), home range sizes were smaller (Kie et al. 2002;Said and Servanty 2005), and utilization rates by deer were higher than would be expected by random (e.g., Moe and Wegge 1994;Tufto et al. 1996;Lamberti et al. 2006) in edge habitats. However, preference for a particular habitat (or vegetation type) does not necessarily lead to a higher reproductive output. For instance, Juncus-dominated marsh was the highest ranked resource for red deer based on relative use versus availability, but it had little effect on reproductive success (McLoughlin et al. 2006). Although there is an increasing number of studies showing the linkage between habitat quality and individual reproductive output in deer (e.g., Conradt et al. 1999;Nilsen et al. 2004;Pettorelli et al. 2005;McLoughlin et al. 2006McLoughlin et al. , 2007), very few studies have demonstrated that higher utilization of forest edges by deer translates into increased reproductive output (but see McLoughlin et al. 2007). Examining the extent to which forest edges affect deer reproduction is crucial for predicting population dynamics and implementing appropriate management plans in forested landscapes.
The sika deer (Cervus nippon) population in the Boso peninsula, central Japan, has been increasing over the past 30 years and is becoming a serious threat to biodiversity and in forest ecosystems (Miyashita et al. 2004;Suzuki et al. 2007) as well as to agricultural activities (Chiba Prefecture 2004). The landscape of Boso peninsula is a mosaic of fine-scale elements including broadleaved forests, Japanese cedar plantations, and agricultural fields. Our previous work has documented that fecal nitrogen levels, a measure of food quality, were positively correlated with forest edge length within a landscape, yet negatively correlated with local deer density ). Accordingly, aside from the landscape structure, local deer density is likely to affect reproduction of deer through food resource depression, as reported in other studies on deer (e.g., Gaillard et al. 2000;Focardi et al. 2002;Keyser et al. 2005;Stewart et al. 2005). To determine if such effects are manifested in the reproductive rate, we first analyzed factors determining the pregnancy rate in female sika deer, with particular attention to forest edges and food availability in the forest understory. Because we estimated food availability using the relationship between deer density and understory plant biomass, it mainly represents the effect of deer density through resource depression. Our analysis also revealed the spatial scale at which significant patterns in deer reproductive rate emerge. To demonstrate the validity of this scale, we examined the home range sizes of freeranging female deer using GIS telemetry. We expected that any positive influence of forest edge on pregnancy rate should be more pronounced at the female home range scale rather than at a larger landscape scale.

Study area
The Boso Peninsula is located in Chiba Prefecture,central Japan (35°N,140°E). This region has a warm temperate climate, with mean temperatures of 25°C in mid-summer and 5°C in mid-winter. Annual precipitation is 2,000-2,400 mm, with almost no snowfall in winter. Broadleaved evergreen forests (Castanopsis sieboldii and Quercus spp.) and coniferous plantations (Cryptomeria japonica and Chamaecyparis obtusa) dominate the vegetation. The forest understory vegetation consists of evergreen shrubs (e.g., Aucuba japonica, Eurya japonica), forbs (e.g., Rubus buergeri), sedges, and various fern species.
Sika deer were once restricted to a small area of 40 km 2 . The population began to increase in the early 1970s, and its distribution has continued to expand, covering a current area of 440 km 2 (Chiba Prefecture 2004). The deer density varies locally, with southern areas generally having higher densities ([20 deer/km 2 ).

Deer pregnancy
We examined the pregnancy status of female sika deer culled from mid-February to mid-March in 2005. Age and fetus presence were determined for 61 females that were obtained from a 340 km 2 region that accounted for approximately three-fourths of the entire distribution area. Additionally, kidney fat index as determined by the ratio of fat mass and fat-free kidney mass multiplied by 100 (Riney 1955) was measured for 56 females to estimate deer condition. The culling location of each individual was recorded on a map (1:50,000). Deer age was determined by examining the annual layers in the tooth cementum ([2 years old) (Scheffer 1950), and by the degree of replacement of milk teeth by permanent teeth (0-2 years old) (Ohtaishi 1980). For the data analysis, only two age classes (yearling and adult) were used because the sika deer is known to show stable age-specific pregnancy rates when becoming adult (Takatsuki 2006).

Deer density index
Distributions of deer fecal pellets were investigated during the winters (December-January) from 1996 to 2005, when fecal decomposition rates were low. The density of fecal pellets [X (per 100 quadrats)] shows a high correlation with local deer density [D (km -2 )] as estimated using a block counting method (Maruyama and Nakama 1983), creating the following relationship: D ¼ 0:76 þ 0:079 X (r 2 = 0.731, n = 14, P \ 0.001; Chiba Prefecture 1998). Thus, the number of fecal pellets found on each census line can be regarded as an index of deer density. In 1996, a total of 92 transects (each 1 km long and 1.6 km apart) were placed to cover the entire distribution area of sika deer. Investigators counted all fecal pellets found within 200 quadrats along each transect; the quadrats were 1 9 1m and were placed at 5-m intervals. Approximately half of all transects were investigated each year alternately, except in 2005 when all transects were surveyed. We defined the deer density index (DDI, m -2 /year) of a line transect as the total number of fecal pellets found from 1996 to 2005 divided by the search area per transect (200 m 2 ) and by census time. Because each transect was 1 km in length, we considered that the local density estimated by the DDI comprised the home ranges of several females (see '' Results'' section).
The distribution of the DDI across the whole area of the deer distribution was estimated at 1-km resolutions, assuming that the DDI was at the geometric center of the 1 km-mesh and represents the mean value of that mesh. The DDI at each mesh was calculated as the inverse distance weighted interpolation from the DDI data within 4 km of the mesh center.
Whether pellet counts are a reliable proxy for deer density has been heavily debated (e.g., Fuller 1991; Morellet et al. 2007). However, the spatial gradient of pellet densities of the Boso deer population was fairly stable in 8 years (see Electronic Supplementary Material S1). If the pellet count did not reflect local densities, such consistent patterns could not have been obtained. Thus, DDI appears to capture the spatial gradient of deer density successfully in our case study.

Food availability
We estimated food availability for deer based on the empirical relationships between the biomass of palatable plants and the DDI, as well as the vegetation classification map on GIS. First, we conducted two series of surveys of the forest-floor vegetation: ''interior'' and ''edge'' surveys. For the interior survey, we placed five quadrats (2 9 2m ) in each of 30 cedar plantations and 35 broad-leaved forests that had varying deer densities and measured the coverage of palatable plants with heights of less than 2 m within each quadrat. The interior survey was conducted in February and September 2005. More detailed methods of this survey were described in Suzuki et al. (2007). The relationship between plant cover (F %) and DDI was described using a logistic equation : We converted plant coverage (F %) into dry biomass (W, g/m 2 ) by the following formula, This relationship was obtained by examining the projective cover and foliage dry biomass of palatable plants in 75 quadrats (2 9 2 m). For the edge survey, we established four sets of four quadrats (1 9 1 m) inside the forest, with each set of quadrats located at 0, 2, 5, and 10 m from the forest edge and measured the projective cover of palatable plants with a height of less than 2 m. This survey was conducted in 14 cedar plantations and 16 broad-leaved forests in August 2006. We initially estimated the biomass of palatable plants (W, g/m 2 ) at each quadrat using Eq. 2, and then fitted them to the following logistic equation, where DFE (m) represents the distance from the forest boundary, implying that plant biomass decreases from the boundary to the forest interior. Parameters of both Eqs. 1 and 3 were estimated separately for conifer plantations and broad-leaved forests in summer and winter. We regarded the total cover of evergreen plants as the cover in winter. Regression analyses were performed using R for Windows 2.2.0. Next, we classified landscape elements of the Boso Peninsula into four categories (open habitats, coniferous plantations, broad-leaved forests and other) based on the most recent vegetation/land use map available on J-IBIS (Ministry of the Environment 2000). Major land-use changes after the last update of J-IBIS was corrected based on an ASTER image taken in March 2004 (Earth Remote Sensing Data Analysis Center 2004). The grid resolutions of the GIS analyses and satellite images were 5 and 15 m, respectively. The area of conifer plantations and broadleaved forests were calculated for every 50-m mesh, which was used for the plant biomass calculation.
The plant biomass in a 50-m mesh (W 50 ) was calculated as: where W i and W e are the biomass in the interior parts ([10 m from the boundaries) and edges (\10 m from the boundaries), respectively. W i was computed from Eqs. 1 and 2, whereas W e was calculated by integrating Eq. 3 for 0-10 m of DEF, given a value of DDI. These computations were made on ArcGIS 8.3 (ESRI Co. Ltd) with the help of Hawth's Tools for Spatial Analyst (Beyer 2004).

Forest edge length
We calculated the length of the forest edge that was likely to influence the reproductive rates of sika deer in our system, because this proved to be the only influential landscape variable determining food quality estimated by fecal nitrogen levels ). The forest edge length was defined as the total length of forest perimeter abutting open habitats and roads with a width [2 m; this was calculated using the GIS data mentioned in the ''Food availability'' section.

Home range
Six females were caught in forested areas where deer densities were relatively high. Three females were equipped with telemetry collars (Lotek GPS-3300S, Ontario, Canada) in early summer (late June to early July), and another three females were equipped in late autumn (mid-November) of 2005. Animals tracked in the summer were located every 30 min, and their collars were retrieved after 4 weeks using a self-drop-off system, whereas individuals tracked in the winter were located every hour and their collars were retrieved after 17 weeks.

Factors determining pregnancy rate
We performed multiple logistic regression analyses to detect factors influencing deer pregnancy rates. The dependent variable was the presence/absence of a fetus for each female, and the independent variables were forest edge length, summer or winter food availability, local deer density, and age (yearling vs. adult). Because deer density has a high correlation with food availability, they were not included in the same regression model. The spatial extent at which independent variables should be extracted is not clear, so we generated a buffer circle with a given radius (100, 200, 300, 400, and 500 m) around each location where females had been caught and then calculated forest edge length (m/ha), food availability (vegetation biomass/ m 2 ) in forests, and the DDI within the buffer area using ArcGIS. We performed model selection for the multiple logistic regressions using the Akaike Information Criterion (AIC). Based on the AIC we also computed Akaike weights (w i ), which represent the probabilities that model i is the best model in the set of models considered. Additionally, we performed supplementary analyses separately for yearling and adult in the same way above, except that we used AICc for model selection criterion because of small sample sizes.
To assess the ability of the best model to discriminate between pregnant and non-pregnant females, we generated an ROC curve and calculated the area under the ROC curve (AUC). AUC values that were [0.7 indicate an acceptable discrimination capacity (Hosmer and Lemeshow 2000).
To examine the difference in kidney fat index between females with and without a fetus, a two-way ANOVA was conducted using the presence/absence of a fetus and age (yearling vs. adult) as two factors.

Home range
The home range of deer was computed using 95 and 50% kernel isopleths (Worton 1989), using reference bandwidth (Worton 1995). These estimates were computed using the R/adehabitat package (Calenge 2006).
To evaluate habitat preference by the deer, a habitat map of the home range area (MCP estimates) was produced using aerial photos taken in September 2003 and GIS software (ESRI ArcMap 8.3). The classification of habitats was the same as in food availability, except that forested areas within 5 m of the forest boundaries were also classified as grassland because of the high productivity of deer forage .
A deer preference index was calculated for each habitat type based on Aebischer et al. (1993). The preference index (p i ) is the proportional use (u i ) of the habitat type inversely weighted by the availability of the habitat type (a i ) (Aebischer et al. 1993): Here habitat availability was defined as the relative area of each habitat type within the home range, and proportional use was defined as the relative number of location points in each habitat type. A few location points for deer no. 6 fell on open water and were excluded from the analysis.

Relationship between plant biomass and deer density
Deer density showed a strong negative effect on the coverage of palatable plants in the forest-floor vegetation, both during summer and winter, and in both the conifer and broad-leaved forests (Table 1, parameter a 1 ). The effect of deer was similar in broad-leaved forests and in conifer plantations, as revealed by the similar slope values (parameter a 1 ), but palatable plants were generally more abundant in conifer forests, as shown in the intercepts (parameter a 0 ). After controlling for the effect of deer density, the biomass of palatable plants in summer showed a marked decrease with increasing distance from the forest boundaries. A similar, but weaker tendency was found in the winter (Table 2, parameter b 1 ).

Factors determining pregnancy rate
To explore the spatial extent at which independent variables could explain the pregnancy rates of sika deer, we show the model with the lowest AIC for each radius of buffer circle (Fig. 1) together with the null model (no variables). For all females combined, buffers with a radius of 100 and 200 m showed similar levels of the lowest AIC; these values were considerably lower than the AIC of the null model (Fig. 1a). However, the lowest AIC increased abruptly at 300 m and came close to the AIC of the null model. In buffers of 400 and 500 m radius, the null model had the lowest AIC. Accordingly, a spatial extent of 100 or 200 m appears to determine deer pregnancy. Supplementary analyses showed that yearlings and adults exhibited the lowest AIC in the buffer radius of 200 and 100 m, respectively (Fig. 1b, c). Table 3a shows the parameters of the five best models and the null model for the 200-m radius buffer (the best spatial extent). In all cases, ''forest edge length'' was always included and had an estimate that was approximately twice as large as its standard error. All other variables that are supposed to affect deer pregnancy, including food availability and deer density, exhibited estimated values that were equal to or lower than their standard errors, indicating that they had a very weak effect on the deer pregnancy rates. When analyzing separately for the age (Table 3b, c), forest edge always included in all models as well. Summer food was additionally included in the best model, but the sign was positive for yearlings and negative for adults. As the estimates of these variables were less than twice as large as their standard errors, they were not strong determinants for deer pregnancy. Figure 2 shows the relationship between the length of forest edge and the pregnancy of females in the model with a 200-m radius buffer; this model had the lowest AIC among all possible models. Almost all deer inhabiting areas with more than 100 m/ha of forest edge became pregnant. The AUC value of the regression was 0.732, indicating satisfactory discrimination capacity.

Home range
The home range of each female deer was estimated using more than 200 location points, indicating adequate sample sizes for kernel estimates (Seaman et al. 1999). The deer exhibited small home-range sizes and core areas both in summer and winter (Table 4), with 95 and 50% kernel estimates being less than 50 and 12 ha, respectively. The deer generally exhibited a higher preference for grassland compared to other landscape elements. Neither body mass nor age appeared to affect home range size and habitat use.

Discussion
Our most important finding is that increases in forest edge in a landscape led to an increased probability of pregnancy in female sika deer. This was demonstrated by the facts that (1) the best regression model included only the variable of forest edge length, and the difference in AICs between the best and null models was approximately nine, indicating substantial support for the best model (Burnham and Anderson 2002), (2) the estimate of forest edge length was more than twice as large as its standard error, whereas the estimates of other variables were equal to or smaller than their standard errors, and (3) the best model was able to satisfactorily discriminate between pregnant and nonpregnant females (AUC [ 0.7). As revealed in previous studies, forest edges generally have higher productivity and larger biomass than forest-floor habitats (e.g., Lyon and Jensen 1980;Alverson et al.1988;Kremsater and Bunnell 1992;Rea 2003). In our study area, agricultural fields adjacent to forests harbor large amounts of plant material (Takada et al. 2002), and clearcuts and roadsides are also rich in plant biomass (unpubl. data). Furthermore, our previous work revealed higher levels of fecal nitrogen in a landscape with more forest edge ). Thus, the high food availability in these habitats appeared to have enhanced nutritional condition, which then facilitated pregnancy rates in sika deer. This inference was supported by the positive correlation between the kidney fat index of females and the forest edge length in a landscape. Although edge length did not clearly affect the upper limit of nutritional condition, it did enhance the lower limit, which probably has increased the pregnancy rate of sika deer. Although numerous studies have demonstrated the role of edge habitats on home range size (Kie et al. 2002;Said and Servanty 2005), habitat use (e.g., Moe and Wegge 1994;Tufto et al. 1996;Lamberti et al. 2006), and local densities (e.g., Wahlstrom and Kjellander 1995;Hemami et al. 2005) of deer, to our knowledge, there is only one other study that provided evidence of the importance of edge habitat on deer reproduction. McLoughlin et al. (2007) recently showed a positive association between lifetime reproductive success of female roe deer and home ranges containing forest edges (meadows and road). The striking similarity of the two results is probably due to the sedentary nature common to sika deer and roe deer. The home range size of sika deer in our study area was only 20-50 ha, which is similar to the home range size of roe deer in European countries (e.g., Said and Servanty 2005;Focardi et al. 2006;Lamberti et al. 2006). Thus, fine-scale variation in habitat quality derived from the sedentary nature of sika deer and roe deer is likely to have caused a strong association between the length of forest edge and reproductive output. Social factors such as the high site fidelity of females may be responsible for the sedentarity (Porter et al. 1991(Porter et al. , 2004Pettorelli et al. 2001).
It is noteworthy that the spatial extent that appears to determine the deer pregnancy rate is approximately 10 ha or less because the regression model with a buffer radius of 200 m around the capture point of deer displayed the  lowest AIC, followed by the model with a 100-m radius buffer. It seems unlikely that spatial extents larger than 300 m were effective because models with the lowest AICs were null models in those extents. The effective spatial extent of 200 m was supported by the home range size of females as estimated using GPS telemetry. Although home range estimates by 95% kernel methods were larger than this extent (20-50 ha vs. 10 ha), the 50% kernel estimate was equivalent to or smaller than this (3-12 ha vs. 10 ha). These results suggest that the pregnancy of sika deer is largely determined by the habitat quality of the core area of its home range and that fine-scale variability in reproductive rates is ascribed to the small range size coupled with fine-scale landscape heterogeneity in this system. Contrary to our expectation, we found no evidence for a density-dependent reduction in deer pregnancy as estimated by food availability in forests (a function of local density) and local density per se, despite the fact that food availability in forests clearly decreased with increasing deer density. This is probably because food availability in open environments is much higher than that in the forest interior and is not depleted by deer foraging due to its high productivity. This suggests that vegetation in open environments provides ample food for deer, and the heavy pressure on forest vegetation by deer did not result in a negative feedback to deer reproduction. The large amount of annual rainfall (2,000-2,400 mm) and warm winter temperature (5°C) appear to be responsible for the high productivity in our study region. An alternative explanation for the absence of density dependence could be that the range of deer densities in our sample was not large enough to detect the negative effects of density. This is unlikely, however, because the range of local deer densities was 2-20 individuals/km 2 . Other studies also revealed that sika deer do not exhibit clear density-dependence unless they reach extremely high population densities leading to population eruption (Putman and Clifton-Bligh 1997;Kaji et al. 2004).
Our results have important implications for the management of sika deer. If deer depend on forest-floor vegetation for their food, the population density should be regulated in a density-dependent manner, and chronic high densities should not be possible. Actually, however, sika deer in this system appear to be sustained by food resources from open environments outside the forest interior, with no obvious decreases by consumption (''donor-controlled'' consumer-resource interaction, e.g., Polis and Strong 1996). This may create continuous negative effects on forest ecosystems even after forest-floor vegetation is severely damaged. This is likely to cause soil erosion and the deterioration of ecosystem function through changes in the physical and chemical properties of the soil (Binkley et al. 2003;Furusawa et al. 2003;Danell et al. 2006). To prevent such ecosystem change, we should increase the awareness of local governments and stakeholders as to how change in land use could lead to deer overabundance and associated negative effects on forests. Although this scenario has long been advocated by many researchers (white-tailed deer: Alverson et al. 1988, Decalesta 1997; black-tailed deer: Kremsater and Bunnell 1992; mule deer: Kie et al. 2002;muntjac: Hemami et al. 2005), it has now become more persuasive in light of our results. We would like to re-emphasize the statement by Schmitz and Sinclair (1997) that the problem of deer overabundance may be attributed to land management issues, rather than the result of a lack of population management.