Predicting Habitat Utilization and Extent of Ecosystem Disturbance by an Increasing Herbivore Population

Herbivory can lead to shifts in ecosystem state or changes in ecosystem functioning, and recovery from herbivory is particularly slow in disturbance-sensitive ecosystems such as arctic tundra. Herbivore impacts on ecosystems are variable in space and time due to population fluctuations and selective utilization of habitats; thus there is a need to accurately predict herbivore impacts at the landscape scale. The habitat utilization and extent of disturbance caused by increasing populations of pink-footed geese (Anser brachyrhynchus) foraging in the high arctic tundra of Svalbard were assessed using a predictive model of the population’s habitat use. Pink-footed geese arrive in Svalbard in early spring when they forage for belowground plant parts; this foraging (called grubbing) can cause vegetation loss and soil disturbance. Surveys of the extent and intensity of grubbing were carried out to develop predictive models that were subsequently tested against data collected during the following year from different areas. Both habitat type at a particular point and the amount of preferred fen habitat in the surrounding area were powerful predictors of grubbing likelihood and the developed model correctly classified over 69% of validation observations with an AUC of 0.75. Pink-footed geese showed a strong preference for wetter habitats within low-lying landscapes. Extrapolation of the predictive model across the archipelago showed that a maximum potential area of 2300 km² (3.8% of the archipelago) could be disturbed by grubbing. Thus, increasing populations of geese may cause large-scale vegetation loss and soil disturbance in arctic ecosystems.     

Extrapolating herbivore-induced carbon loss across an arctic landscape

A high proportion of the global soil carbon stock is stored in tundra soils. However, populations of arctic-breeding migratory geese including pink-footed geese, Anser brachyrhynchus, are increasing due to agricultural changes and conservation measures in their wintering grounds. Foraging by these geese, which is widespread in extent, reduces the quantity of carbon stored in arctic tundra ecosystems. Here, the potential carbon loss caused by foraging pink-footed geese is modelled across the high-arctic archipelago of Svalbard, combining field experiments, habitat maps and published spatial models of foraging. The carbon loss caused by foraging geese was estimated three growing seasons following perturbation allowing for some recovery to take place. The carbon loss caused by 1-year worth of grubbing was estimated to be 1,700 tonnes, or 37 kg per goose. A total of over 340,000 tonnes of carbon could be affected given an unlimited increase in goose population. Estimated losses were mostly from wetter habitats, which are both carbon rich and highly selected for by foraging geese. The across-landscape carbon loss caused by geese is not great in magnitude in comparison to expected climate-driven carbon losses; however, it is locally severe, and demonstrates how migratory connectivity links processes, such as agricultural change and conservation measures in temperate Europe with carbon dynamics in the high arctic.     

Modelling pink-footed goose (Anser brachyrhynchus) wintering distributions for the year 2050: potential effects of land-use change in Europe

Feeding on farmland by overwintering populations of pink-footed geese ( Anser brachyrhynchus ) conflicts with agricultural interests in Northern Europe. In order to forecast the potential future of this conflict, we used generalized linear models to relate the presence and absence of pink-footed geese to variables describing the contemporary landscape, and predicted their future distributions in relation to two land-use scenarios for the year 2050. One future scenario represented a global, economically orientated world (A1) and the other represented a regional, environmentally concerned world (B2). The probability of goose occurrence increased within cropland and grassland, and could be explained by their proximity to coast, elevation, and the degree of habitat closure. Predictions to future scenarios revealed noticeable shifts in the suitability of goose habitat evident at the local and regional scale in response to future shifts in land use. In particular, as grasslands and croplands give way to unsuitable land-use types (e.g. woody biofuel crops, increased urbanization, and forest) under both future scenarios, our models predicted a decrease in habitat suitability for geese. If coupled with continued goose population expansion, we expect that the agricultural conflict will intensify under the A1 and particularly the B2 scenarios.     

Taxonomy of the bean goose-pink-footed goose

The bean goose Anser fabalis and the pink-footed goose A. brachyrhynchus breed in the tundra and taiga zones of Eurasia and eastern Greenland, and the taxonomy of the group based on morphology has been controversial. We investigated the phylogenetic relationships within the bean goose–the pink-footed goose complex using mitochondrial control region sequences of 199 individuals collected from the breeding areas in the Palaearctic and Eastern Nearctic. We found three mitochondrial clades geographically distributed to (1) Greenland, Iceland and Svalbard (A. brachyrhynchus), (2) the eastern taiga zone (former subspecies A. fabalis middendorffii), and (3) the western taiga and the tundra zone (subspecies A. fabalis rossicus, serrirostris and fabalis). MtDNA phylogeny suggests that morphological affinities between the taxa, e.g. in the bill structure, result from convergent evolution due to adaptation to similar habitats. Although a latitudinal cline in morphology was observed, clear phylogenetic discontinuities exist in the taiga and tundra zones supporting a species status for brachyrhynchus and middendorffii.     

First de novo whole genome sequencing and assembly of the pink-footed goose

Annotated genomes can provide new perspectives on the biology of species. We present the first de novo whole genome sequencing for the pink-footed goose. In order to obtain a high-quality de novo assembly the strategy used was to combine one short insert paired-end library with two mate-pair libraries. The pink-footed goose genome was assembled de novo using three different assemblers and an assembly evaluation was subsequently performed in order to choose the best assembler. For our data, ALLPATHS-LG performed the best, since the assembly produced covers most of the genome, while introducing the fewest errors. A total of 26,134 genes were annotated, with bird species accounting for virtually all BLAST hits. We also estimated the substitution rate in the pink-footed goose, which can be of use in future demographic studies, by using a comparative approach with the genome of the chicken, the mallard and the swan goose. A substitution rate of 1.38 × 10^? 7 per nucleotide per generation was obtained when comparing the genomes of the two closely-related goose species (the pink-footed and the swan goose). Altogether, we provide a valuable tool for future genomic studies aiming at particular genes and regions of the pink-footed goose genome as well as other bird species.     

Colonization history of the high-arctic pink-footed goose Anser brachyrhynchus

Population structure and phylogeography of the pink-footed goose, Anser brachyrhynchus Baillon 1833, was studied using mtDNA control region sequences (221 bp) from 142 individuals. Present breeding areas of the species in Greenland, Iceland, and Svalbard were largely covered by ice during the late Pleistocene. In pairwise comparisons phiST estimates showed significant differentiation among eastern and western populations, whereas sampling localities within both areas were not differentiated. The mtDNA data indicate that the populations have separated recently (less than 10 000 years ago) and present breeding areas were colonized from one refugial population. The levels of haplotype and nucleotide diversity were approximately five times higher for the eastern population compared to the western population and suggest that the latter was colonized by a subset of eastern birds. Time to the most recent common ancestor of the species is 32 000-46 000 years, i.e. the present mtDNA variation of the pink-footed goose has accumulated during the last 0.1 My. Estimates of the long-term female effective population size (5400-7700 for the eastern population) imply that the refugial population of the pink-footed goose has been large. Tundra habitats were more extensive in cold periods of the late Pleistocene than today and may have sustained population sizes that allowed the accumulation of extant genetic polymorphism. It is not probable that the postulated small refugial areas in the high latitudes had a significant role in maintaining this diversity.     

Predator behaviour and predation risk in the heterogeneous Arctic environment

1. Habitat heterogeneity and predator behaviour can strongly affect predator-prey interactions but these factors are rarely considered simultaneously, especially when systems encompass multiple predators and prey. 2. In the Arctic, greater snow geese Anser caerulescens atlanticus L. nest in two structurally different habitats: wetlands that form intricate networks of water channels, and mesic tundra where such obstacles are absent. In this heterogeneous environment, goose eggs are exposed to two types of predators: the arctic fox Vulpes lagopus L. and a diversity of avian predators. We hypothesized that, contrary to birds, the hunting ability of foxes would be impaired by the structurally complex wetland habitat, resulting in a lower predation risk for goose eggs. 3. In addition, lemmings, the main prey of foxes, show strong population cycles. We thus further examined how their fluctuations influenced the interaction between habitat heterogeneity and fox predation on goose eggs. 4. An experimental approach with artificial nests suggested that foxes were faster than avian predators to find unattended goose nests in mesic tundra whereas the reverse was true in wetlands. Foxes spent 3.5 times more time between consecutive attacks on real goose nests in wetlands than in mesic tundra. Their attacks on goose nests were also half as successful in wetlands than in mesic tundra whereas no difference was found for avian predators. 5. Nesting success in wetlands (65%) was higher than in mesic tundra (56%) but the difference between habitats increased during lemming crashes (15%) compared to other phases of the cycle (5%). Nests located at the edge of wetland patches were also less successful than central ones, suggesting a gradient in accessibility of goose nests in wetlands for foxes. 6. Our study shows that the structural complexity of wetlands decreases predation risk from foxes but not avian predators in arctic-nesting birds. Our results also demonstrate that cyclic lemming populations indirectly alter the spatial distribution of productive nests due to a complex interaction between habitat structure, prey-switching and foraging success of foxes.     

Finding food in a highly seasonal landscape: where and how pink footed geese Anser brachyrhynchus forage during the Arctic spring

Food accessibility and availability in the highly seasonal Arctic landscape can be restricted by snow cover and frozen substrate, particularly in early spring. Therefore, to determine how a long distance migratory herbivore forages in such a landscape, pink‐footed geese Anser brachyrhynchus at an early spring feeding area in Svalbard were studied. Birds arrived in mid May when extensive snow cover restricted habitat availability. Geese fed in all habitats, but the highest densities occurred in wet tundra. However, prolonged snow lie restricted access to wet areas compared to dry and mesic habitats. Above ground biomass was very low in all habitats; yet sizeable amounts occurred below ground. In line with this, the majority of birds (86%) grubbed for below ground plant storage organs such as stem bases and rhizomes. Wet habitat contained greater quantities of edible and lower amounts of inedible below ground material (ratio 1:0.3) than dry (ratio 1:9) or mesic (ratio 1:4) areas. Although foraging in wet habitat prevented geese from encountering high proportions of inedible plant parts, forage species characteristic of this habitat, such as Dupontia grasses and the rush Eriophorum scheuchzeri, were more difficult to extract than food plants typical of drier habitats such as the forb Bistorta vivipara. Hence, we suggest that wet areas are preferred by pink‐footed geese, but the prolonged snow lie there made it necessary to use less preferred but much more abundant drier habitats, which experienced earlier snowmelt and indeed accommodated more than half of all goose foraging recordings.     

Habitat correlates of the Eurasian otter Lutra lutra recolonizing Central Poland

The increase in Eurasian otter Lutra lutra populations in their natural range and recolonization processes are recently observed in several European countries. We address the process of otter recolonization and habitat utilization in Central Poland over 14 years. Field surveys in 1998 and 2007 documented increase in occurrence of the species. The frequency of positive sites denoted 15 % in 1993, 38 % in 1998, and 89 % in 2007. Otter occurrence at study sites was positively affected by river width while negatively affected by presence of buildings at the site and river regulation. During the most intensive colonization process in the 1990s, the habitat preferences of the otter did not change. However, the sites inhabited by otters after 1998 were characterized by lower river width and tree cover and were more often located on regulated river sections, suggesting change in habitat tolerance during expansion. The otter abundance in transformed habitats is a result of increasing population numbers and the necessity to inhabit suboptimal sections of watercourses. Thus, it seems that presence–absence data for otter populations cannot be considered a reliable indicator of habitat quality, being depended of the population density.     

Predictive modelling to aid the regional-scale management of a vertebrate pest

Extensive resources are allocated to managing vertebrate pests, yet spatial understanding of pest threats, and how they respond to management, is limited at the regional scale where much decision-making is undertaken. We provide regional-scale spatial models and management guidance for European rabbits (Oryctolagus cuniculus) in a 260,791 km² region in Australia by determining habitat suitability, habitat susceptibility and the effects of the primary rabbit management options (barrier fence, shooting and baiting and warren ripping) or changing predation or disease control levels. A participatory modelling approach was used to develop a Bayesian network which captured the main drivers of suitability and spread, which in turn was linked spatially to develop high resolution risk maps. Policy-makers, rabbit managers and technical experts were responsible for defining the questions the model needed to address, and for subsequently developing and parameterising the model. Habitat suitability was determined by conditions required for warren-building and by above-ground requirements, such as food and harbour, and habitat susceptibility by the distance from current distributions, habitat suitability, and the costs of traversing habitats of different quality. At least one-third of the region had a high probability of being highly suitable (support high rabbit densities), with the model supported by validation. Habitat susceptibility was largely restricted by the current known rabbit distribution. Warren ripping was the most effective control option as warrens were considered essential for rabbit persistence. The anticipated increase in disease resistance was predicted to increase the probability of moderately suitable habitat becoming highly suitable, but not increase the at-risk area. We demonstrate that it is possible to build spatial models to guide regional-level management of vertebrate pests which use the best available knowledge and capture fine spatial-scale processes.     

Better up,worse down: bidirectional consequences of three decades of climate change on a relict population of Erebia cassioides

Orophilous species are often unable to escape the consequences of climate change because mountains are surrounded by unsuitable habitats. Among them, several endemic species belonging to the genus Erebia Dalman (Lepidoptera, Nymphalidae, Satyrinae) can be considered as key species to assess the risk of biodiversity loss of mountain habitats. The aim of this paper is to measure changes that have occurred in the altitudinal distribution of Erebia cassioides on the Pollino Massif (Southern Italy) during the last 37 years. Sixteen sites sampled in 1975 have been resampled after about three decades (2004, 2012). In 1975 56 % of the sampled population inhabited sites above and 44 % sites below the treeline, while in 2004 and 2012 99 % of the population were observed above the treeline. Furthermore, we observed an uphill shift of 180 m in the barycentre altitude of the species distribution and an unexpected increased density of the population above the treeline which led to a range reduction coupled to population increase of E. cassioides. This pattern contrasts with the usually observed one that couples habitat reduction to population decreasing. The reason for the observed pattern is unclear, but the implication for conservation strategies could be important if confirmed for other species. In fact, during coming decades local extinctions as a consequence of climate change might be fewer and more delayed than expected, and relict populations of cold adapted species could be preserved for a longer time span within optimal habitat refugia.     

Towards a Solution to the Goose-Agriculture Conflict in North Norway, 1988–2012: The Interplay between Policy,Stakeholder Influence and Goose Population Dynamics

This paper presents results from a multidisciplinary study of a negotiation process between farmers and wildlife authorities which led to an agricultural subsidy scheme to alleviate conflicts between agriculture and geese in Norway. The Svalbard-breeding population of pink-footed geese Anser brachyrhynchus has increased considerably over the last decades and conflicts with farmers have escalated, especially at stopover sites in spring when geese feed on newly sprouted pasture grass. In Vesterålen, an important stopover site for geese in North Norway, farmers deployed scaring of geese at varying intensity dependent on the level of conflict during 1988–2012. We assessed the efficiency of a subsidy scheme established in 2006, in terms of its conflict mitigation, reflected in a near discontinuation of scaring activities. The presence of pink-footed geese was analysed in relation to scaring intensity, the total goose population size and the increasing occurrence of another goose species, the barnacle goose Branta leucopsis. Scaring significantly affected the number of geese staging in Vesterålen, both in absolute and relative terms (controlling for total population size). The geese responded immediately to an increased, and reduced, level of scaring. Despite the establishment of the subsidy scheme, the number of pink-footed geese has recently declined which is probably caused by the increasing number of barnacle geese. For the farmers, the subsidy scheme provides funding that reduces the economic costs caused by the geese. Sustaining a low level of conflict will require close monitoring, dialogue and adaptation of the subsidy scheme to cater for changes in goose population dynamics.     

Effects of environmental and anthropogenic drivers on Amur tiger distribution in northeastern China

We examined environmental and anthropogenic factors drive range loss in large mammals, using presence data of Amur tigers opportunistically collected between 2000 and 2012, and anthropogenic and environmental variables to model the distribution of the Amur tiger in northeastern China. Our results suggested that population distribution models of different subregions showed different habitat factors determining tiger population distribution patterns. Where farmland cover was over 50 km² per pixel (196 km²), distance was within 15 km to the railway in Changbaishan and road density (length per pixel) increased in Wandashan, the relative probability of Amur tiger occurrence exhibited monotonic avoidance responses; however, where distance was within 150 km of the Sino-Russia border, the occurrence probability of Amur tiger was relatively high. We analyzed the avoidance or preference responses of Amur tiger distribution to elevation, snow depth and Viewshed. Furthermore, different subregional models detected a variety of spatial autocorrelation distances due to different population clustering patterns. We found that spatial models significantly improved model fits for non-spatial models and made more robust habitat suitability predications than that of non-spatial models. Consequently, these findings provide useful guidance for habitat conservation and management.     

Population densities and habitat use of the golden jackal (Canis aureus) in farmlands across the Balkan Peninsula

The main objective of this study was to analyze the habitat use and population densities of the golden jackal in four countries across lowland regions of the Balkan Peninsula, known as the core area of the species' distribution in Europe. Using indirect (acoustic) method for detecting territorial golden jackals, we surveyed jackal presence and densities on 331 monitoring sites in four countries, covering area an of 4,296 km² in total during April and May 2007–2012. We used GIS to assess landscape and environmental characteristics in a 2-km circular buffer (12.6 km²) around calling stations. Average population density of golden jackals in the study areas ranged between 0.6 and 1.1 territorial groups/10 km² (mean ± SE, 0.6?±?0.06 groups/10 km²), with several high-density areas with up to 4.8 territorial groups/10 km². Analysis of habitat use showed that for both jackal occurrence and number of jackal groups, the only significant parameter was the interaction between country and intensity of agriculture, indicating that jackals adapt their habitat selection patterns in relation to the habitat availability. We observed that selection of the more suitable habitats (shrub–herbaceous vegetation/heterogeneous agricultural vegetation) increased with lower proportion of these habitat types in the study area. Our study confirms high habitat plasticity of the golden jackal and offers explanation for its recent range expansion, which might be connected with the land use changes during the last decades in the Balkan Peninsula.     

Spring feeding by pink-footed geese reduces carbon stocks and sink strength in tundra ecosystems

Tundra ecosystems are widely recognized as precious areas and globally important carbon (C) sinks, yet our understanding of potential threats to these habitats and their large soil C store is limited. Land‐use changes and conservation measures in temperate regions have led to a dramatic expansion of arctic‐breeding geese, making them important herbivores of high‐latitude systems. In field experiments conducted in high‐Arctic Spitsbergen, Svalbard, we demonstrate that a brief period of early season belowground foraging by pink‐footed geese is sufficient to strongly reduce C sink strength and soil C stocks of arctic tundra. Mechanisms are suggested whereby vegetation disruption due to repeated use of grubbed areas opens the soil organic layer to erosion and will thus lead to progressive C loss. Our study shows, for the first time, that increases in goose abundance through land‐use change and conservation measures in temperate climes can dramatically affect the C balance of arctic tundra.     

Integrating aspects of ecology and predictive modelling: implications for the conservation of the leopard cat (Prionailurus bengalensis) in the Eastern Himalaya

An understanding of species ecology is vital for effective conservation, particularly if the species forms an important constituent of the lesser mammal guild and regulates small mammal and bird populations. As the ecological role of the leopard cat (Prionailurus bengalensis) in the intricate eastern Himalayan habitats is not known, we assessed the site occupancy, detection probability and activity pattern of leopard cats in Khangchendzonga Biosphere Reserve, India, based on sign surveys and camera trapping. The estimated site occupancy was 0.352?±?0.061 and detection probability was 0.143?±?0.0484. Occupancy modelling indicated low elevation, high rodent abundance and tree cover as best predictors for the occupancy of leopard cat. Diet based on analysed scats revealed murids as the most dominant prey (89.2 %). Information based on photographic captures indicated that the leopard cat exhibited a nocturnal activity pattern (peak activity between 0200–0300 hours), which coincided with its principal prey (revealed through diet analysis), but mainly contradicted with other sympatric competitors, hence indicating a temporal partitioning of resources among them. Ecological niche factor analysis indicated that the leopard cat exhibits high global marginality (1.32) and low global tolerance (0.275). The habitat suitability map for leopard cats showed majority of the habitat as unsuitable (1,959.44 km²) and predicted only 164.54 km² areas of lower temperate forests as moderate to highly suitable. As highly suitable habitats of the leopard cat are in close proximity to villages, conflict issues are a major threat and therefore need to be addressed in conservation program for this felid.     

Establishment of a new population by a single gravid colonist: implications for Hyalella biogeography and speciation

Establishment of a new population is a significant event in a lineage, one that impacts geographic distribution, community dynamics, and speciation. We examined population establishment in experimental populations of Hyalella, an amphipod crustacean, by introducing a single gravid female into experimental habitats that differed in habitat size and the presence of an interspecific competitor, Physella snails. We found that over the course of 142 days, the length of a typical reproductive season, a single female gave rise to a viable population of descendants in 62 % of experimental habitats, with more than half of these viable populations containing more than 100 descendants. We found no effect of habitat size or interspecific competition on extinction probability or final abundance of descendants. Based on empirical results we estimated the probability that a single female will establish a lasting population, and examined this establishment probability under a broad range of overwinter survival scenarios. We found that establishment probability declines gradually with reductions in overwinter survival, such that establishment probability exceeds 33 % across a broad range of overwinter mortality rates. Our finding that a single female may have a reasonably high probability of establishing a new population helps to resolve the enigmatic biogeographic distribution of Hyalella amphipods, which appear to have very limited overland dispersal, yet are nearly ubiquitous in permanent freshwater habitats of North America.     

Small rodents in the shrub tundra of Yamal (Russia): Density dependence in habitat use?

Northern small rodents are well known for their population cycles which represent a key process for the functioning of arctic and boreal ecosystems. Habitat use often changes in the course of the cycle. Higher densities can either lead to spill-over into secondary habitats or to increased habitat specificity because of interspecific competition. Here we investigate whether voles in the shrub tundra of southern Yamal exhibit density dependent habitat use. Voles were trapped at the Erkuta Tundra Monitoring Site (N 68.2°, E 69.2°) in three characteristic habitats over five years covering all phases of the population cycle. Our analyses focused on the two most numerous species Microtus gregalis (52% of individuals caught) and M. middendorffii (36%). A small-scale spill-over effect was observed for M. gregalis, which increasingly used the open habitat adjacent to their preferred willow thickets at high abundance. At a larger scale no such effect was observed for the two Microtus species – a result which is explained by the overall moderate densities of voles and the large spatial extent of the primary habitat of M. middendorffii: moist moss dwarf shrub tundra.     

Lack of competition between barnacle geese Branta leucopsis and pink-footed geese Anser brachyrhynchus during the pre-breeding period in Svalbard

The feeding ecology of barnacle geese and pink-footed geese was studied in Sassendalen, Svalbard during the pre-nesting period (late May) to assess the potential for inter-specific competition. Barnacle geese fed almost exclusively (97%) by grazing above-ground plant material, mostly (79%) along snow edges in moss-mat habitats. Pink-footed geese fed mostly (93%) by excavating below-ground parts of plants, mostly (56%) away from snow and were more evenly distributed between habitat types. Barnacle goose faeces contained mostly (62%) moss, that of pink-footed geese mostly (48%) below-ground plant storage organs (especially Bistorta viviparum L.). Principal components analysis of dropping contents showed no overlap in species diet in allopatry or sympatry. There was little overlap in diet and feeding ecology of the two species at this pre-nesting feeding site. Hence, unless increased goose feeding densities affect future vegetation density and composition, under present circumstances, increasing numbers of either species is unlikely to affect foraging conditions for the other at this important stage in the annual cycle. However, such changes could have local density-dependent intra-specific effects.     

黑龙江省老爷岭南部穆棱林区马鹿种群分布数量及生境适宜性评价

开展野生动物种群分布数量和生境适宜性评价研究是制定物种科学保护决策的基础和关键。以东北马鹿(Cervus canadensis xanthopygus)为研究对象,于2015—2020年在黑龙江省老爷岭南部穆棱林区,采用大样方调查方法收集马鹿在雪地留下的足迹信息分析马鹿种群数量;通过相机监测、足迹链跟踪和大样方调查3种方法综合收集马鹿活动点信息,利用最大熵(MaxEnt)建模分析马鹿生境适宜性。马鹿种群分布数量研究结果表明,穆棱林区马鹿种群密度为(0.0645±0.009)只/km~2,种群数量为47—61只,主要分布在研究地区南部的和平、龙爪沟和共和林场。生境分析结果表明：人为干扰因子中,居民区对MaxEnt模型的贡献率为44%,马鹿主要在距居民区距离约5 km和10—15 km的区域活动,在偏僻的林间小道和乡村道路马鹿生境适宜性较高,其对模型贡献率分别为16.8%和10.2%;植被因子中,在距常绿针叶林和针阔混交林距离4 km范围内,随距离增加马鹿生境适宜性逐渐降低。生境适宜性分析结果表明,研究地区马鹿适宜生境和次适宜生境面积为87.09 km~2,仅占研究区域的10.39%,主要...