基于遥感数据量化黄河三角洲刺槐林的边缘效应

刘思思, 王红, 车纯广, 王鑫

南京林业大学学报(自然科学版) ›› 2025, Vol. 49 ›› Issue (6) : 37-46.

PDF(4043 KB)
PDF(4043 KB)
南京林业大学学报(自然科学版) ›› 2025, Vol. 49 ›› Issue (6) : 37-46. DOI: 10.12302/j.issn.1000-2006.202411008
专题报道Ⅰ: 第十四届海峡两岸森林经理研讨会专题(执行主编 李凤日 曹林)

基于遥感数据量化黄河三角洲刺槐林的边缘效应

作者信息 +

Quantification of the edge effect of Robinia pseudoacacia plantation in the Yellow River delta based on remote sensing data

Author information +
文章历史 +

摘要

【目的】边缘效应是指在不同生态系统交界处,由于生态因子的差异及相互作用,引发系统内某些组分显著变化的现象。这一效应广泛存在于森林、湿地、草原等生态系统中。研究森林边缘效应对生态保护、营林管理、气候变化应对和碳储量计算具有重要意义。【方法】以黄河三角洲刺槐人工林为研究对象,利用UAV-LiDAR和Sentinel-2遥感数据,分别提取7个结构指标和3个生理指标,量化分析人工林三维结构和生态功能的边缘效应,并探讨不同边缘类型的边缘效应特点。【结果】①边缘效应对刺槐人工林分结构产生了负面影响。边缘区域的林分树高降低1~5 m、覆盖度减少7%~53%、叶面积指数(LAI)下降10%~77%,林分密度下降6%~52%,同时其冠幅面积和体积也减小。②边缘效应对森林生态功能也产生了负面影响,如削弱了植被的光合作用,其吸收的光合有效辐射(FAPAR)降低8%~50%,叶绿素含量下降9%~72%,冠层含水量也下降8%~50%。③不同边缘类型的影响范围不一样,河流边缘的影响范围可达150 m,而道路边缘影响强度较小,范围约为30 m。【结论】边缘效应对黄河三角洲滨海防护林的植被结构和生态系统功能产生了多维度影响。本研究有助于滨海防护林管理和生态保护,帮助应对气候变化与区域环境退化。

Abstract

【Objective】The edge effect refers to the phenomenon that some components within the ecosystems change significantly due to the differences and interactions of ecological factors at the junction of different ecosystems. This phenomenon is ubiquitously observed across diverse ecosystems including forests, wetlands and grasslands. The investigation of forest edge effects holds substantial importance for advancing ecological conservation efforts, enhancing forestry management practices, addressing climate change impacts, and refining carbon storage estimations.【Method】This study quantified edge effects in Robinia pseudoacacia plantations in the Yellow River delta using UAV-LiDAR and Sentinel-2 data. Seven structural and three physiological indices were analyzed to assess 3D structural and ecological functional variations, with distinct edge effect patterns identified across different edge types.【Result】(1) The edge effect significantly altered vegetation structural parameters, with observed reductions of 1-5 m in tree height, 7%-53% in vegetation coverage, 10%-77% in leaf area index (LAI), and 6%-52% in tree density. Canopy structural features, including projected area and volume, also exhibited marked decreases. (2) Ecological functions were substantially compromised at forest edges, manifested by impaired photosynthetic capacity with 8%-50% reduction in fraction of absorbed photosynthetically active radiation (FAPAR), 9%-72% decline in chlorophyll content, and 8%-50% decrease in canopy water content. (3) Edge influence zones showed distinct spatial patterns, with river edges demonstrating an extensive impact range of 150 meters, contrasting with the relatively limited 30-meter influence zone of road edges.【Conclusion】The edge effect exerts multidimensional influences on both structural characteristics and ecosystem functionalities of coastal shelterbelts in the Yellow River Delta. This research provides critical insights for enhancing coastal shelterbelt management and conservation strategies, facilitating the optimization of sustainable forest plantation development, and contributing to climate change mitigation and regional environmental restoration efforts.

关键词

滨海防护林 / 刺槐人工林 / 边缘效应 / 植被结构 / 生态系统功能 / 黄河三角洲

Key words

coastal shelter forest / Robinia pseudoacacia plantation / edge effect / vegetation structure / ecosystem function / Yellow River delta

引用本文

导出引用
刘思思, 王红, 车纯广, . 基于遥感数据量化黄河三角洲刺槐林的边缘效应[J]. 南京林业大学学报(自然科学版). 2025, 49(6): 37-46 https://doi.org/10.12302/j.issn.1000-2006.202411008
LIU Sisi, WANG Hong, CHE Chunguang, et al. Quantification of the edge effect of Robinia pseudoacacia plantation in the Yellow River delta based on remote sensing data[J]. Journal of Nanjing Forestry University (Natural Sciences Edition). 2025, 49(6): 37-46 https://doi.org/10.12302/j.issn.1000-2006.202411008
中图分类号: S727   

参考文献

[1]
王伯荪, 彭少麟. 鼎湖山森林群落分析:Ⅹ.边缘效应[J]. 中山大学学报(自然科学版), 1986, 25(4):31-38.
WANG B S, PENG S L. Analysis on the forest communities of Dinghushan Guangdong:Ⅹ.Communities edge effect[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 1986, 25(4):31-38.
[2]
PFEIFER M, LEFEBVRE V, PERES C A, et al. Creation of forest edges has a global impact on forest vertebrates[J]. Nature, 2017, 551(7679):187-191.DOI:10.1038/nature24457.
[3]
FRANKLIN C M A, HARPER K A, CLARKE M J. Trends in studies of edge influence on vegetation at human-created and natural forest edges across time and space[J]. Canadian Journal of Forest Research, 2021, 51(2):274-282.DOI:10.1139/cjfr-2020-0308.
[4]
MEEUSSEN C, GOVAERT S, VANNESTE T, et al. Structural variation of forest edges across Europe[J]. Forest Ecology and Management, 2020,462:117929.DOI:10.1016/j.foreco.2020.117929.
[5]
MURCIA C. Edge effects in fragmented forests:implications for conservation[J]. Trends in Ecology & Evolution, 1995, 10(2):58-62.DOI:10.1016/S0169-5347(00)88977-6.
[6]
ANDERSON J A, MCCLEAN C J, SIM S, et al. Weak edge effects on trees in Bornean rainforest remnants bordering oil palm[J]. Biotropica, 2022, 54(4):921-932.DOI:10.1111/btp.13115.
[7]
DANTAS DE PAULA M, GROENEVELD J, HUTH A. The extent of edge effects in fragmented landscapes:insights from satellite measurements of tree cover[J]. Ecological Indicators, 2016, 69:196-204.DOI:10.1016/j.ecolind.2016.04.018.
[8]
LI J P, ZHAO C Y, PENG Y Y, et al. Edge effects on tree growth and species diversity in forests of different types and ages[J]. Polish Journal of Ecology, 2018, 66(3):239.DOI:10.3161/15052249pje2018.66.3.004.
[9]
BLANCHARD G, BARBIER N, VIEILLEDENT G, et al. UAV-Lidar reveals that canopy structure mediates the influence of edge effects on forest diversity,function and microclimate[J]. Journal of Ecology, 2023, 111(7):1411-1427.DOI:10.1111/1365-2745.14105.
\n \n \n Widespread forest loss and fragmentation dramatically increases the proportion of forest areas located close to edges. Although detrimental, the precise extent and mechanisms by which edge proximity impacts remnant forests remain to be ascertained.
[10]
MAEDA E E, NUNES M H, CALDERS K, et al. Shifts in structural diversity of Amazonian forest edges detected using terrestrial laser scanning[J]. Remote Sensing of Environment, 2022,271:112895.DOI:10.1016/j.rse.2022.112895.
[11]
MUTANGA O, DUBE T, OMER G. Assessing edge effect on the spatial distribution of selected forest biochemical properties derived using the Worldview data in Dukuduku forests,south Africa[J]. African Journal of Ecology, 2019, 57(3):314-326.DOI:10.1111/aje.12604.
[12]
BAUER L, HUTH A, BOGDANOWSKI A, et al. Edge effects in Amazon forests:integrating remote sensing and modelling to assess changes in biomass and productivity[J]. Remote Sensing, 2024, 16(3):501.DOI:10.3390/rs16030501.
[13]
HARPER K A, MACDONALD S E, BURTON P J, et al. Edge influence on forest structure and composition in fragmented landscapes[J]. Conservation Biology, 2005, 19(3):768-782.DOI:10.1111/j.1523-1739.2005.00045.x.
[14]
BUSS J, DABROS A, HIGGINS K L, et al. Comparison of edge effects from well pads and industrial roads on mixed upland boreal forest vegetation in Alberta[J]. Plant Ecology, 2024, 225(4):331-343.DOI:10.1007/s11258-023-01393-3.
[15]
MARCHAND P, HOULE G. Spatial patterns of plant species richness along a forest edge:what are their determinants?[J]. Forest Ecology and Management, 2006, 223(1/2/3):113-124.DOI:10.1016/j.foreco.2005.10.064.
[16]
ARAGÓN G, ABUJA L, BELINCHÓN R, et al. Edge type determines the intensity of forest edge effect on epiphytic communities[J]. European Journal of Forest Research, 2015, 134(3):443-451.DOI:10.1007/s10342-015-0863-5.
[17]
SMITH I A, HUTYRA L R, REINMANN A B, et al. Piecing together the fragments:elucidating edge effects on forest carbon dynamics[J]. Frontiers in Ecology and the Environment, 2018, 16(4):213-221.DOI:10.1002/fee.1793.
Forest fragmentation is pervasive throughout the world's forests, impacting growing conditions and carbon (C) dynamics through edge effects that produce gradients in microclimate, biogeochemistry, and stand structure. Despite the majority of global forests being <1 km from an edge, our understanding of forest C dynamics is largely derived from intact forest systems. Edge effects on the C cycle vary by biome in their direction and magnitude, but current forest C accounting methods and ecosystem models generally fail to include edge effects. In the mesic northeastern US, large increases in C stocks and productivity are found near the temperate forest edge, with over 23% of the forest area within 30 m of an edge. Changes in the wind, fire, and moisture regimes near tropical forest edges result in decreases in C stocks and productivity. This review explores differences in C dynamics observed across biomes through a trade‐offs framework that considers edge microenvironmental changes and limiting factors to productivity.
[18]
HAMBERG L, LEHVÄVIRTA S, KOTZE D J. Forest edge structure as a shaping factor of understorey vegetation in urban forests in Finland[J]. Forest Ecology and Management, 2009, 257(2):712-722.DOI:10.1016/j.foreco.2008.10.003.
[19]
ORDWAY E M, ASNER G P. Carbon declines along tropical forest edges correspond to heterogeneous effects on canopy structure and function[J]. Proceedings of the National Academy of Sciences of the United States of America, 2020, 117(14):7863-7870.DOI:10.1073/pnas.1914420117.
Nearly 20% of tropical forests are within 100 m of a nonforest edge, a consequence of rapid deforestation for agriculture. Despite widespread conversion, roughly 1.2 billion ha of tropical forest remain, constituting the largest terrestrial component of the global carbon budget. Effects of deforestation on carbon dynamics in remnant forests, and spatial variation in underlying changes in structure and function at the plant scale, remain highly uncertain. Using airborne imaging spectroscopy and light detection and ranging (LiDAR) data, we mapped and quantified changes in forest structure and foliar characteristics along forest/oil palm boundaries in Malaysian Borneo to understand spatial and temporal variation in the influence of edges on aboveground carbon and associated changes in ecosystem structure and function. We uncovered declines in aboveground carbon averaging 22% along edges that extended over 100 m into the forest. Aboveground carbon losses were correlated with significant reductions in canopy height and leaf mass per area and increased foliar phosphorus, three plant traits related to light capture and growth. Carbon declines amplified with edge age. Our results indicate that carbon losses along forest edges can arise from multiple, distinct effects on canopy structure and function that vary with edge age and environmental conditions, pointing to a need for consideration of differences in ecosystem sensitivity when developing land-use and conservation strategies. Our findings reveal that, although edge effects on ecosystem structure and function vary, forests neighboring agricultural plantations are consistently vulnerable to long-lasting negative effects on fundamental ecosystem characteristics controlling primary productivity and carbon storage.
[20]
LAURANCE W F, NASCIMENTO H E M, LAURANCE S G, et al. Habitat fragmentation,variable edge effects,and the landscape-divergence hypothesis[J]. PLoS One, 2007, 2(10):e1017.DOI:10.1371/journal.pone.0001017.
[21]
LAURANCE W F, CAMARGO J L C, FEARNSIDE P M, et al. An Amazonian rainforest and its fragments as a laboratory of global change[J]. Biological Reviews, 2018, 93(1):223-247.DOI:10.1111/brv.12343.
\n We synthesize findings from one of the world's largest and longest‐running experimental investigations, the\n B\n iological\n D\n ynamics of\n F\n orest\n F\n ragments\n P\n roject (\n BDFFP\n ). Spanning an area of\n ∼\n 1000 km\n 2\n in central\n A\n mazonia, the\n BDFFP\n was initially designed to evaluate the effects of fragment area on rainforest biodiversity and ecological processes. However, over its 38‐year history to date the project has far transcended its original mission, and now focuses more broadly on landscape dynamics, forest regeneration, regional‐ and global‐change phenomena, and their potential interactions and implications for\n A\n mazonian forest conservation. The project has yielded a wealth of insights into the ecological and environmental changes in fragmented forests. For instance, many rainforest species are naturally rare and hence are either missing entirely from many fragments or so sparsely represented as to have little chance of long‐term survival. Additionally, edge effects are a prominent driver of fragment dynamics, strongly affecting forest microclimate, tree mortality, carbon storage and a diversity of fauna.\n
[22]
PÖPPERL F, SEIDL R. Effects of stand edges on the structure,functioning,and diversity of a temperate mountain forest landscape[J]. Ecosphere, 2021, 12(8):e03692.DOI:10.1002/ecs2.3692.
[23]
WANG Z J, YANG J Y. Urbanization strengthens the edge effects on species diversity and composition of woody plants in remnant forests[J]. Forest Ecosystems, 2022,9:100063.DOI:10.1016/j.fecs.2022.100063.
[24]
YUAN C M, GENG Y F, CHAI Y, et al. Response of lianas to edge effects in mid-montane moist evergreen broadleaved forests in the Ailao Mountains,SW China[J]. Biodiversity Science, 2016, 24(1):40-47.DOI:10.17520/biods.2015248.
[25]
NGUYEN T A, EHBRECHT M, CAMARRETTA N. Application of point cloud data to assess edge effects on rainforest structural characteristics in tropical Sumatra,Indonesia[J]. Landscape Ecology, 2023, 38(5):1191-1208.DOI:10.1007/s10980-023-01609-x.
In forestry, edge zones created by forest degradation and fragmentation are more susceptible to disturbances and extreme weather events. The increase in light regime near the edge can greatly alter forest microclimate and forest structure in the long term. In this context, understanding edge effects and their impact on forest structure could help to identify risks, facilitate forest management decisions or prioritise areas for conservation.\n
[26]
SILVA JUNIOR C H L, ARAGÃO L E O C, ANDERSON L O, et al. Persistent collapse of biomass in Amazonian forest edges following deforestation leads to unaccounted carbon losses[J]. Science Advances, 2020, 6(40):eaaz8360.DOI:10.1126/sciadv.aaz8360.
[27]
MELIN M, HINSLEY S A, BROUGHTON R K, et al. Living on the edge:utilising lidar data to assess the importance of vegetation structure for avian diversity in fragmented woodlands and their edges[J]. Landscape Ecology, 2018, 33(6):895-910.DOI:10.1007/s10980-018-0639-7.
[28]
VEPAKOMMA U, KNEESHAW D D, DE GRANDPRÉ L.Influence of natural and anthropogenic linear canopy openings on forest structural patterns investigated using LiDAR[J]. Forests, 2018, 9(9):540-559.DOI:10.3390/f9090540.
[29]
GULCI S, AKAY A E, OGUZ H, et al. Assessment of the road impacts on coniferous species within the road-effect zone using NDVI analysis approach[J]. Fresenius Environmental Bulletin, 2017, 26(2a): 1654-1662.
[30]
HARPER K A, MACDONALD S E, MAYERHOFER M S, et al. Edge influence on vegetation at natural and anthropogenic edges of boreal forests in Canada and Fennoscandia[J]. Journal of Ecology, 2015, 103(3):550-562.DOI:10.1111/1365-2745.12398.
\n\n\nAlthough anthropogenic edges are an important consequence of timber harvesting, edges due to natural disturbances or landscape heterogeneity are also common. Forest edges have been well studied in temperate and tropical forests, but less so in less productive, disturbance‐adapted boreal forests.
[31]
HONG M M, WANG J L, HAN B M. Spatial and temporal pattern changes and driving forces:analysis of salinization in the Yellow River delta from 2015 to 2020[J]. Journal of Resources and Ecology, 2022, 13(5):786-796.DOI:10.5814/j.issn.1674-764x.2022.05.004.
[32]
CHEN H Y, ZHAO G X, LI Y H, et al. Monitoring the seasonal dynamics of soil salinization in the Yellow River delta of China using Landsat data[J]. Natural Hazards and Earth System Sciences, 2019, 19(7):1499-1508.DOI:10.5194/nhess-19-1499-2019.
[33]
WU C S, LIU G H, HUANG C. Prediction of soil salinity in the Yellow River delta using geographically weighted regression[J]. Archives of Agronomy and Soil Science, 2017, 63(7):928-941.DOI:10.1080/03650340.2016.1249475.
[34]
ZHAO X Q, GUO Q H, SU Y J, et al. Improved progressive TIN densification filtering algorithm for airborne LiDAR data in forested areas[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2016, 117:79-91.DOI:10.1016/j.isprsjprs.2016.03.016.
[35]
BOHN F J, HUTH A. The importance of forest structure to biodiversity-productivity relationships[J]. Royal Society Open Science, 2017, 4(1):160521.DOI:10.1098/rsos.160521.
\n While various relationships between productivity and biodiversity are found in forests, the processes underlying these relationships remain unclear and theory struggles to coherently explain them. In this work, we analyse diversity–productivity relationships through an examination of forest structure (described by basal area and tree height heterogeneity). We use a new modelling approach, called ‘forest factory’, which generates various forest stands and calculates their annual productivity (above-ground wood increment). Analysing approximately 300 000 forest stands, we find that mean forest productivity does not increase with species diversity. Instead forest structure emerges as the key variable. Similar patterns can be observed by analysing 5054 forest plots of the German National Forest Inventory. Furthermore, we group the forest stands into nine forest structure classes, in which we find increasing, decreasing, invariant and even bell-shaped relationships between productivity and diversity. In addition, we introduce a new index, called optimal species distribution, which describes the ratio of realized to the maximal possible productivity (by shuffling species identities). The optimal species distribution and forest structure indices explain the obtained productivity values quite well (\n R\n 2\n between 0.7 and 0.95), whereby the influence of these attributes varies within the nine forest structure classes.\n
[36]
BROCKERHOFF E G, BARBARO L, CASTAGNEYROL B, et al. Forest biodiversity, ecosystem functioning and the provision of ecosystem services[J]. Biodiversity and Conservation, 2017, 26(13):3005-3035.DOI:10.1007/s10531-017-1453-2.
[37]
HADDAD N M, BRUDVIG L A, CLOBERT J, et al. Habitat fragmentation and its lasting impact on earth's ecosystems[J]. Science Advances, 2015, 1(2):e1500052.DOI:10.1126/sciadv.1500052.
Urgent need for conservation and restoration measures to improve landscape connectivity.
[38]
ZHAO Z, LI W, CIAIS P, et al. Fire enhances forest degradation within forest edge zones in Africa[J]. Nature Geoscience, 2021,14:479-483.DOI:10.1038/s41561-021-00763-8.
[39]
SUNG W G, LEE D K, JIN Y H. Analyzing difference of urban forest edge vegetation condition by land cover types using spatio-temporal data fusion method[J]. Journal of Environmental Impact Assessment, 2018, 27(3): 279-290. DOI:10.14249/eia.2018.27.3.279.
[40]
宋音, 王红, 路开宇, 等. 基于CCA方法的黄河三角洲不同健康刺槐林的土壤属性研究[J]. 江西农业学报, 2017, 29(10):48-53.
SONG Y, WANG H, LU K Y, et al. Study on soil properties of Robinia pseudoacacia forest under different health conditions in Yellow River delta based on canonical correspondence analysis[J]. Acta Agriculturae Jiangxi, 2017, 29(10):48-53.DOI:10.19386/j.cnki.jxnyxb.2017.10.11.
[41]
WANG H, LU K Y, PU R L. Mapping Robinia pseudoacacia forest health in the Yellow River delta by using high-resolution IKONOS imagery and object-based image analysis[J]. Journal of Applied Remote Sensing, 2016, 10(4):045022.DOI:10.1117/1.jrs.10.045022.
[42]
TAYLOR P J, NUBERG I K, HATTON T J. Enhanced transpiration in response to wind effects at the edge of a blue gum (Eucalyptus globulus) plantation[J]. Tree Physiology, 2001, 21(6):403-408.DOI:10.1093/treephys/21.6.403.
In Australia, tree planting has been widely promoted to alleviate dryland salinity and one proposed planting configuration is that of strategically placed interception belts. We conducted an experiment to determine the effect of tree position in a belt on transpiration rate. We also assessed how much the effect of tree position can be explained by advection and environmental conditions. Daily transpiration rates were determined by the heat pulse velocity technique for four edge and 12 inner trees in a 7-year-old Tasmanian blue gum (Eucalyptus globulus) plantation in South Australia. Various climatic variables were logged automatically at one edge of the plantation. The relationship between daily sap flow and sapwood area was strongly linear for the edge trees (r2 = 0.97), but only moderately correlated for the inner trees (r2 = 0.46), suggesting an edge effect. For all trees, sap flow normalized to sapwood area (Qs) increased with potential evaporation (PE) initially and then became independent as PE increased further. There was a fairly close correlation between transpiration of the edge and inner trees, implying that water availability was partially responsible for the difference between inner and edge trees. However, the ratio of edge tree to inner tree transpiration differed from unity, indicating differences in canopy conductance, which were estimated by an inverse form of the Penman-Monteith equation. When canopy conductances were less than a critical value, there was a strong linear relationship between Qs of the edge and inner trees. When canopy conductances of the edge trees were greater than the critical value, the slope of the linear relationship was steeper, indicating greater transpiration of the edge trees compared with the inner trees. This was interpreted as evidence for enhancement of transpiration of the edge trees by advection of wind energy.
[43]
ALMEIDA D R A, STARK S C, SCHIETTI J, et al. Persistent effects of fragmentation on tropical rainforest canopy structure after 20 yr of isolation[J]. Ecological Applications, 2019, 29(6):e01952.DOI:10.1002/eap.1952.
Assessing the persistent impacts of fragmentation on aboveground structure of tropical forests is essential to understanding the consequences of land use change for carbon storage and other ecosystem functions. We investigated the influence of edge distance and fragment size on canopy structure, aboveground woody biomass (AGB), and AGB turnover in the Biological Dynamics of Forest Fragments Project (BDFFP) in central Amazon, Brazil, after 22+ yr of fragment isolation, by combining canopy variables collected with portable canopy profiling lidar and airborne laser scanning surveys with long‐term forest inventories. Forest height decreased by 30% at edges of large fragments (>10 ha) and interiors of small fragments (<3 ha). In larger fragments, canopy height was reduced up to 40 m from edges. Leaf area density profiles differed near edges: the density of understory vegetation was higher and midstory vegetation lower, consistent with canopy reorganization via increased regeneration of pioneers following post‐fragmentation mortality of large trees. However, canopy openness and leaf area index remained similar to control plots throughout fragments, while canopy spatial heterogeneity was generally lower at edges. AGB stocks and fluxes were positively related to canopy height and negatively related to spatial heterogeneity. Other forest structure variables typically used to assess the ecological impacts of fragmentation (basal area, density of individuals, and density of pioneer trees) were also related to lidar‐derived canopy surface variables. Canopy reorganization through the replacement of edge‐sensitive species by disturbance‐tolerant ones may have mitigated the biomass loss effects due to fragmentation observed in the earlier years of BDFFP. Lidar technology offered novel insights and observational scales for analysis of the ecological impacts of fragmentation on forest structure and function, specifically aboveground biomass storage.
[44]
DELJOUEI A, SADEGHI S M M, ABDI E, et al. The impact of road disturbance on vegetation and soil properties in a beech stand,Hyrcanian forest[J]. European Journal of Forest Research, 2018, 137(6):759-770.DOI:10.1007/s10342-018-1138-8.

基金

国家自然科学基金项目(31971579)

编辑: 郑琰燚
PDF(4043 KB)

Accesses

Citation

Detail

段落导航
相关文章

/