气候及环境的变化与物种的地理分布和空间分布关系密切[1]，对物种的分布范围和繁殖发育产生深刻影响，因此，生物的气候需求及其与物种地理分布间的关系成为制定物种保护策略的重要依据。目前，模型分析被广泛应用于物种居群的高适生性分布区预测，其中MaxEnt模型可用来揭示物种居群或未知居群的分布点[2]，已被广泛用于濒危物种的适生区预测、气候环境适宜性、以及用于物种保护优先性评价等方面。如金丝猴[3]、丹顶鹤[4]等动物，伯乐树(Bretschneidera sinensis)、蒙古扁桃(Amygdalus mongolica)、太白红杉(Larix chinensis)、秦岭冷杉(Abieschensiensis)、独叶草(Kingdoniauniflora)[5-9]，及印度东北地区特有濒危冬青树(Ilex khasiana)、高山杨梅(Canacomyrica monticola)[10-11]等植物。

毛红椿(Toona ciliata var. pubescens)属楝科香椿属植物，被誉为“中国桃花心木”，1991年被列为国家二级保护濒危物种，主要分布于我国华南地区，如江西、湖南、湖北、安徽、广东、广西、四川、贵州、海南、浙江和云南等省，最北端位于安徽南部泾县的汀溪、苏红。目前毛红椿大面积天然林分已不复存在，多以斑块状零星分布为主[12]。相关研究表明天然毛红椿自然更新不良[13-16]，以根系萌蘖繁殖为主[12,17]；可能是由于生境破碎化导致其自然生殖的隔绝，遗传多样性普遍较低[13,18-19]。但目前，关于毛红椿地理分布格局，以及气候环境因子与其地理适宜性潜在分布的关系尚未深入研究。因此，本文利用MaxEnt模型预测其潜在自然分布，分析国内外分布格局，揭示毛红椿的生境需求及适宜的环境因子，为建立核心保护区域及其保护措施提供依据。

1 数据收集与研究方法

1.1 数据的收集与处理

1.1.1 毛红椿分布数据获取 主要取自全国科考调查、期刊文献、植物学类专著，以及本研究于江西省境内的调查成果，确定当前22个分布点的数据资料(表1)。

1.1.2 气候数据的获取 环境数据的19个气候因子(表2)来自网址http://www.worldclim.org/。

当时的我，对于图书分类编目一无所知。无奈之下，在图书馆借了一本《俄文编目手册》和《图书分类法》，在家里整整啃了两天。元旦后一上班，我就按照书上阐述的分类原则和图书著录格式，着手俄文图书的分编工作。同时向其他同志学习打蜡纸、油印卡片的方法和技巧。经过一年的边学边做，我慢慢忘却了走下三尺讲台的“失落感”，渐渐适应并爱上了图书编目工作。第二年我就被评上了校先进工作者，并在《图书情报工作》期刊上先后发表了《全苏科技情报中心现状》和《苏联科学院图书馆现状》两篇文章。撰写文章使我认识到了图书情报工作在科研工作中的重要地位与作用。同时也认识到，只要刻苦钻研、勤奋努力，在图书馆工作也是可以有所作为的。

1.2 数据处理和气候环境变量选取

1.2.1 数据处理 采用MaxEnt3.3.3k软件预测毛红椿当前的潜在适宜区。将地理分布数据以物种名称、经纬度的方式输入，存储成后缀名为csv文档格式。分布区图像采自中国地图、世界地图，并对气象数据进行影像的配准、裁剪和叠加等运行产生。

Guo X H,Zhou S D,Chen R J,et al.Survey of Toona ciliata var. pubescens community at Ting county in Anhui province[J].Journal of Anhui Normal University(Natural Science Edition),1983,27(2):63-71.

 

表1 毛红椿潜地理分布点基本情况

 

Tab.1 Basic information of the potential distribution of T. ciliata var. pubescens population

  

地区Region居群编号与地点Numberandlocationofpopulations纬度Latitude经度Longitude参考文献References江西省龙南县LongnancountyinJiangxiprovince1．九连山自然保护区冷水坑24°34′114°27′2．九连山自然保护区石塘24°36′114°30′3．九连山自然保护区虾公堂主沟24°31′114°27′4．九连山自然保护区廖屋坑24°33′114°26′江西省宜丰县YifengcountyinJiangxiprovince5．官山自然保护区28°30′114°29′刘军[13]6．官山自然保护区28°37′114°35′温强等[19]江西省抚州市FuzhoucityinJiangxiprovince7．马头山自然保护区27°40′117°08′温强等[19]江西省井冈山市JinggangshancityinJiangxiprovince8．井冈山自然保护区26°35′114°14′江西省黎川县LichuancountyinJiangxiprovince9．岩泉森林公园27°04′117°04′刘军[13]云南省宾川县BinchuancountyinYunnanprovince10．宾川力角乡25°02′100°16′刘军[13]云南省师宗县ShizongcountyinYunnanprovince11．五龙乡拔菌子山24°21′103°42′刘军[13]云南元谋县YunmoucountyinYunnanprovince12．元谋县25°17′101°49′刘军[13]安徽泾县JingxiancountyinAnhuiprovince13．苏红乡桃岭小泥坑30°26′115°36′郭新狐等[12]14．汀溪乡大坑30°33′118°37′郭新狐等[12]四川攀西地区PanxiregionsinSichuanprovince15．德昌县惠民乡27°21＇102°16′武华卫等[16]16．德昌县小高乡27°19＇102°19′武华卫等[16]17．延边县惠民乡26°45＇101°49′武华卫等[16]18．延边县渔门镇26°52＇101°38′武华卫等[16]19．米易县垭口镇26°52＇102°6′武华卫等[16]20．米易县丙谷镇26°48＇102°6′武华卫等[16]浙江仙居县XianjucountyinZhejiangprovince21．俞坑自然保护区28°48＇120°32′刘军等[18]浙江省遂昌县SuichangcountyinZhejiangprovince22．九龙山自然保护区28°19＇118°49′刘军等[18]

 

表2 气候环境变量描述

 

Tab.2 Description of climate environmental variable

  

气候环境变量Climateenvironmentalvariable气候环境变量Climateenvironmentalvariable年平均气温度Meanannualtemperature(Bio1)最冷季度平均温度Meantemperatureofcoldestquarter(Bio11)昼夜温差月均值Meandiurnalrange[Meanofmonthly(maxtemperature-mintemperature)(Bio2)年均降雨量Annualprecipitation(Bio12)昼夜温差与年温差比值Isothermality(Bio2/Bio7×100)(Bio3)最湿月降雨量Precipitationofwettestmonth(Bio13)温度变化方差SDoftemperatureseasonality(Bio4)最干月降雨量Precipitationofdriestmonth(Bio14)最热月份最高温Maxtemperatureofwarmestmonth(Bio5)湿度变化方差SDofhumidityseasonality(Bio15)最冷月份最低温Mintemperatureofcoldestmonth(Bio6)最湿季降雨量Precipitationofwettestquarter(Bio16)年温度变化范围Temperatureannualrange(Bio5-Bio6)(Bio7)最干季降雨量Precipitationofdriestquarter(Bio17)最湿季度平均温度Meantemperatureofwettestquarter(Bio8)最热季度平均降雨量Precipitationofwarmestquarter(Bio18)最干季度平均温度Meantemperatureofdriestquarter(Bio9)最冷季度平均降雨量Precipitationofcoldestquarter(Bio19)最热季度平均温度Meantemperatureofwarmestquarter(Bio10)

1.2.3 气候变量筛选和提取 受试者工作特性曲线(receiver operating characteristic curve，ROC曲线) 面积，以AUC值(area under roc curve)的大小评价模型的准确性，选取AUC>0.9评价本模型的预测效果。采用“Do jackknife to measure variable importance”功能检验和评价气候各环境因子的相对贡献率，筛选影响毛红椿地理分布的主导气候因子。

急性阑尾炎是在多种因素下所形成的炎性改变，大多因感染所引发。腹痛是其典型的临床症状，而单纯性阑尾炎则表现为持续性或阵发性的钝痛或胀痛，甚至波及双侧下腹部或中下腹，给患者的生活造成严重影响[7-9]。祖国医学认为阑尾炎多因饮食不节、寒温不适、忧思抑郁等因素所导致的肠道功能紊乱，同时伴随气血不和、传化不利等症状。临床主张行手术切除已病变的阑尾，辅以穴位按压及口服中药汤药等可大大缓解患者的临床症状，更好达到活血化瘀、清热解毒的功效，并可增加其机体的免疫能力，利于炎症的消除[10-12]。

2 结果与分析

2.1 ROC 曲线与预测精度

深灰色线代表训练数据的ROC曲线，浅灰色线为测试数据的ROC曲线，曲线下面积AUC 值越接近于1，说明预测物种“存在”的能力越高。图1中AUC值为 0.959±0.036，说明该模型模拟毛红椿的地理分布结果不是随机存在的，而是有规律的，表明气候因素对其分布有重要影响，且预测结果具有很高的可信度。

2.2 毛红椿潜在地理分布格局及其适宜性等级

理论上，存在内控缺陷的企业更有可能在其财务报表中出现重大错报，注册会计师将面临更高的审计风险，为了降低风险只有通过实施更多的审计程序降低检查风险，成本将增加，审计定价也将更高。由此，提出假设一：H1：审计定价与企业内控缺陷之间存在正相关关系。

  

图1 MaxEnt 模型的ROC 检验Fig.1 ROC test of MaxEnt model

  

图2 毛红椿在中国的潜在适生区地理分布Fig.2 Potential suitable distribution in China

  

图3 毛红椿在全球的潜在适生区地理分布Fig.3 Potential suitable distribution in the world

图2显示，毛红椿在我国的潜在分布格局呈连续性，适宜及其以上分布区集中，且高分布值(0.77～1.0)区域呈片状(除台湾最适分布区外)，其中，2个紧邻的最适分布区域分别是安徽和浙江的东南部，以及江西南部、福建、广东。而在国外的潜在生态分布区，则包括喜马拉雅山西北坡、印度东部、孟加拉国、缅甸、泰国和越南的连续狭长腹地林区，间断连接至朝鲜、韩国与日本的阔叶林带，且国外的最高分布值超过我国(图3)。全球高分布区之间存在地理空间分布的断裂，具有相对独立性，这种斑块状分布有利于保持适应各自环境的遗传稳定性，但也影响了高适区优良种群间的基因交流。模拟结果显示，毛红椿全球潜在分布既包括了已知的分布[12]，如新近报道的福建南平市武夷山和延平区等地[20]，也包括新增的我国台湾分布区，且大于实际已知的生长区，其中以我国境内分布的连续性较好(图2～3)。

Wu W W,Gu S S,Wu J,et al.Impact of climate change on distribution of breeding sites of Red-Crowned Crane in China[J].Journal of Ecology and Rural Environment,2012,28(3):243-248.

 

表3 气候变量重要性的Jackknife检验

 

Tab.3 Results of Jackknife test of climate environmental variable importance

  

变量Variable贡献率/%Percentcontribution随机分布重要值Permutationimportance变量Variable贡献率/%Percentcontribution随机分布重要值Permutationimportancebio1865．022．6bio-90．10．0bio421．43．9bio-160．10．1bio67．467．1bio-70．10．2bio191．51．1bio-130．10．0bio-151．50．0bio-100．00．0bio-141．21．2bio-170．00．0bio-80．60．3bio-30．00．0bio-10．62．1bio-120．00．0bio-20．31．3bio-50．00．0

  

图4 3个主要气候因子的响应曲线Fig.4 Response curves of 3 main climate factor

2.3 毛红椿潜在分布的关键气候因子

经Jackknife 检验、AUC检验及环境变量的贡献率分析，表明有3 个生物气候主导因子影响毛红椿的分布，分别为bio-18、bio-4和bio-6，其贡献率分别为 65.0%、21.4%和7.4%，对应的随机分布重要值分别为22.6、3.9、67.1(表3，图4～6)。即最热季度平均降雨量(Bio18)494.79～1 504 mm，最适值在717.16 mm，温度变化方差(Bio4)2 362.88～9 986.69，最适4 904.15，最冷月份最低温(Bio6)-10.35～10.44 ℃，最适-0.49 ℃(图4)。这些主导因子与我国目前已知分布区亚热带气候特征一致，可作为毛红椿引种及其种群保护的重要依据。

  

图5 生物气候环境变量对MaxEnt 模型的贡献值Fig.5 The contribution value of Biology climate environmental variable to MaxEnt model

  

图6 生物气候环境变量的AUC检验Fig.6 AUC test of biology climate environmental variable

3 结论与讨论

基于毛红椿已知的生态位MaxEnt模型，补充完善了已有的分布记录，较好地预测和识别毛红椿的潜在生态位，以及较高精度反映了该植物的分布和环境条件之间的匹配性，其地理分布的主要气候限制因子依次是最热季降水量、温度变化和极端最低温。根据目前调查与报道的毛红椿分布资料，我国主要分布区有江西、云南、四川、安徽、湖南、广西、广东、福建等省境的南方或西南林区，集中在中亚热带、南亚热带气候带，因而具有分布的连续性特性，也与毛红椿阳性喜光，喜温暖、湿润环境需求一致。据气象资料统计，江西毛红椿主要分布区为官山、马头山和九连山国家自然保护区等地，年平均气温10 ℃以上，1月份平均气温8.3 ℃，毛红椿保存和生长发育良好。安徽的泾县、宣城和宁国等县分布点，森林资源丰富的阔叶林小气候条件优良[12]，夏季平均降水量495～587 mm，最冷月份最低温-13.35～14.7 ℃，1月份平均气温2.6～3.0 ℃。另报道低温影响毛红椿自然更新繁衍[17]，以及移栽南京受到低温的限制[22]。可见，毛红椿经长期适应亚热带气候带环境，导致其对降雨、低温、温度变化的敏感性，与本文模拟得到的最热季度平均降雨量、温度变化和最冷月份最低温的3个主导因子结论相吻合，可以推断上述气象因子是决定毛红椿适宜分布及其生存环境的主要影响因子。

从全球区域来看，毛红椿生态空间分布格局整体连续，但高分布区的斑块化特征明显。这有利于生境内物种具有较高的生物多样性和生态系统稳定性，并维持较高的遗传多样性，这种斑块内生物多样性的维持机制是就地保护的重要依据和基础资料[21]。毛红椿现有残存分布及破碎化生境产生的生殖隔绝导致较高的遗传分化程度和偏低的遗传多样性[13,18]，并不利于该物种的保护和种质资源库构建。在本研究模拟的高分布区，与目前已知的遗传多样性较高的研究区域(江西、浙江、云南)[13,18,19]相吻合。从经典岛屿理论来看，毛红椿气候生物学高适宜的斑块化生境，对其物种多样性分布格局的形成可能发挥着重要作用。因此，需深入研究毛红椿各高分布区域，优先保护区域内适宜生境及优势群落，分析斑块生境的“面积效应”及可能存在的特殊变化规律。

参考文献：

[1] Walther G R,Post E,Convey P,et al.Ecological response to recent climate change[J].Nature,2002,416:389-395.

[16] 武华卫,杨柳璐,钟金,等.攀西地区毛红椿资源分布现状初步调查[J].四川林业科技,2015,36(5):37-40.

坑道工程口部伪装是工程防护重要研究方向之一，一些工程在施工过程中破坏了原有地貌，伪装时需要轻质仿石器材来模拟口部附近的真山石，造成口部是一些普通石块的假象，从而形成有效欺骗。随着园林行业的发展、仿石技术的进步，在可见光波段模拟山石已非难事。然而军事侦察主要波段除了可见光外还有近红外、热红外、雷达等波段[1-3]，而真山石的光学、热红外和雷达波段特征具有独特性，园林仿石很容易被揭露。所以，能起到伪装作用的仿石器材需要一些新技术、新材料才能达到对抗多频谱侦察的需要。

[3] 徐卫华,罗翀.MAXENT 模型在秦岭川金丝猴生境评价中的应用[J].森林工程,2010,26(2):1-4.

Xu W H,Luo C.Applicationof MAXENT model in Rhinopithecus roxellanae habitat assessment in Qinling Mountain[J].Forest Engineering,2010,26(2):1-4.

[4] 吴伟伟,顾莎莎,吴军,等.气候变化对我国丹顶鹤繁殖地分布的影响[J].生态与农村环境学报,2012,28(3):243-248.

党的十八大以来，云南省法院紧紧围绕“努力让人民群众在每一个司法案件中感受到公平正义”的目标，坚持司法为民、公正司法，充分发挥人民法院审判职能作用，为乡村经济发展、社会和谐稳定作出了积极贡献。2013年到2017年，131个基层法院审（执）结各类案件1460286件，占全省法院总结案数的75.15%。

马克·吐温是著名的批判现实主义作家，他的笔下描绘了很多形形色色的社会人物形象，写实与幽默的同时又常暗含嘲讽之意。目前对马克·吐温短篇小说的研究，基本可以概括为以下两个方面：

[5] 龚维,夏青,陈红锋,等.珍稀濒危植物伯乐树的潜在适生区预测[J].华南农业大学学报,2015,36(4):98-104.

Gong W,Xia Q,Chen H F,et al.Prediction of potential distributions of Bretschneidera sinensis，an rare and endangered plant species in China[J].Journal of South China Agricultural University,2015,36(4):98-104.

明朝首辅张居正爱摆谱，虽摆出了一时得意，却摆掉了荣华富贵。为了省亲时在家乡父老面前炫耀自己，他不仅仪仗队伍浩浩荡荡，还有100名火枪手护卫，并斥巨资做了32人抬大轿，精美绝伦，有客厅，有卧室，有厨房，还有金童玉女伺候，极尽奢侈之能事，且一路招摇，收礼无数。一餐百菜，尚嫌“无下箸处”。谱摆得太大，不仅大臣纷纷上书弹劾，皇帝也极不满意，张居正一死，就立即下令抄家，削尽其宫秩，迫夺生前所赐玺书、诰命，张居正家人或饿死或流放，惨不忍睹。

[6] 马松梅,聂迎彬,段霞,等．蒙古扁桃植物的潜在地理分布及居群保护优先性[J].生态学报,2015，35(9):2960-2966.

Ma S M,Nie Y B,Duan X,et al.The potential distribution and population protection priority of Amygdalus mongolica[J].Acta Ecologica Sinica,2015,35(9):2960-2966.

[7] 高蓓,卫海燕，郭彦龙,等.应用 GIS 和最大熵模型分析秦岭冷杉潜在地理分布[J].生态学报,2015,34(3):843-852.

Gao B,Wei H Y,Guo Y L et al.Using GIS and MaxEnt to analyze the potential distribution of Abies chensiensis[J].Acta Ecologica Sinica,2015,34(3):843-852.

[8] 徐军,曹博,白成科.基于MaxEnt 濒危植物独叶草的中国潜在适生分布区预测[J].生态学报,2015,34(12):3354 -3359.

Xu J,Cao B,Bai C K.Prediction of potential suitable distribution of endangered plant Kingdonia uniflora in China with MaxEnt[J].Chinese Journal of Ecology,2015,34(12):3354 -3359.

[9] 吴建国.气候变化对 7种保护植物分布的潜在影响[J].武汉植物学研究,2010,28(4):437-452.

Wu J G.Potential effects of climate change on the distribution of seven protected plants in China[J].Journal of Wuhan Botanical Research,2010,28(4):437-452.

潜在分布阈值是由训练集数据中10%的随机分布点产生，分布值处于0.3以下的地区视为小概率事件。因此，可将毛红椿适宜度划分3 个分布等级，分别由浅灰色、灰色和黑色表示，其中0.39～0.62为适宜分布区(浅灰色)，0.63～0.85为较适分布区(灰色)，0.86～以上最适分布区(黑色)(图2～3)。

[10] Kumar S,Stohlgren T J.Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticoal in New Caledonia[J].Journal of Ecology and Natural Environment,2009,1:94-98.

[11] Adhikari D,Barik S K,Upadhaya K．Habitat distribution modeling for reintroduction of Ilex khasiana Purk,a critically endangered tree species of northeastern India[J].Ecological Engineering,2012,40:37-43.

[12] 郭新弧,周树德,陈仁钧,等.安徽泾县毛红椿林的群落调查[J].安徽师大学报:自然科学版,1983,27(2):63-71.

1.2.2 模型运行 基于一系列环境变量产生一个连续栅格的预测分布图，其概率值从 0到 1，代表研究区域内对当地气候环境的适宜性。模型设置重复运行 15 次，每次运行随机选择 25%的分布信息点为测试集，剩余 75%的样本为训练集用于验证。

[13] 刘军.毛红椿天然居群遗传结构研究[D].北京:中国林业科学研究院,2007.

乔瞧把鞋踢了，就要朝水里跳，秀容月明叫了声“别忙，我先下去”。秀容月明自幼练武，身子健壮，不怕冷。他下了水，把手伸给乔瞧：“我拉着你，你慢点下来。”

使用Stata 11.0软件对所提取的资料信息进行Meta‐分析。本研究的结局指标为二分类变量，以比值比（OR）及其95%置信区间（95%CI）作为合并效应量。采用Q检验及I2检验对所纳入研究的异质性进行检验，如果P≤0.10或I2＞50%时，表明研究间异质性较大，将采用随机效应模型进行结局指标效应值的合并分析；反之则采用固定效应模型进行合并分析。采用漏斗图定性判断发表偏倚并采用Egger's检验定量判断发表偏倚，P＜0.05时，表明存在发表偏倚。

Liu J.Study on genetic structure of natural populations of Toona ciliata var. pubescens[D].Beijing:Forest Tree Genetics and Breeding,2007.

[14] 付方林.江西九连山毛红椿群落学特征研究[D].南昌:江西农业大学,2008.

[17] 萧运峰.安徽省毛红楝子的调查报告[J].植物生态学报,1983,7(2):152-157.

[15] 黄红兰.九连山自然保护区毛红椿种群生态学特征[D].南昌:江西农业大学,2012.

Huang H L.The population ecological characteristics of Toona ciliata var. pubescens in Jiulianshan national reserve[D].Nanchang:Jiangxi Agricultural University,2012.

[2] Elith J,Phillips S J,Hastie T,et al.A statistical explanation of MaxEnt for ecologists[J].Diversity and Distributions,2011,17:43-57.

“突突突……”11月的一天早晨，在监利县何王庙长江江豚省级自然保护区，一条铁驳船缓缓驶入河汊。协巡员周家喜把太阳帽扔进船舱，穿上救生衣和下水裤，跳下船去，从水里扯起一根小竹竿，带起一张10多米长的渔网。“江豚很有可能误闯进这种渔网中受到伤害。”周家喜将渔网从水中慢慢拖上岸，密封的下水裤和“江豚保护”救生衣让他感到闷热，额头很快冒出了汗珠。

Wu H W,Yang L L,Zhong J,et al.Preliminary investigation on the Distribution of Toona ciliata var. pubescens resources in Panxi areas[J].Journal of Sichuan Forestry Science and Technology,2015,36(5):37-40.

Fu F L.Study on the community characteristics of Toona ciliata var. pubescens in Jiulianshan national reserve[D].Nanchang:Jiangxi Agricultural University,2008.

Xiao Y F.Investigation on Toona ciliata var. pubescens from Anhui Province[J].Chinese Journal of Plant Ecology,1983,7(2):152-157.

[18] 刘军,陈益泰,罗阳富,等.毛红椿天然林群落结构特征研究[J].林业科学研究,2010,23(1):93-97.

Liu J,Chen Y T,Luo Y F,et al.Studyon community characters of Toona ciliata var. pubescens natural forest[J].Forest Research,2010,23(1):93-97.

[19] 温强,叶金山,周诚.江西毛红椿天然群体的遗传多样性分析[J].南方林业科技,2016,44(3):1-6.

Wen Q,Ye J S,Zhou C.Analysis of genetic diversity of natural populations of Toona ciliata var. pubescens in Jiangxi Province[J].South China Forestry Science,2016,44(3):1-6.

（2）进行心理健康教育：繁重的工作会导致施工人员不可避免的产生抑郁、焦躁等心理疾病，严重的影响其在施工过程中的安全性。因此本工程通过定期的进行视频教育的方式，警示施工人员不执行安全操作可能出现的严重后果，同时通过开展各种娱乐活动，加强施工人员良好工作情绪的培养，减少不必要的施工误操作，有效的提高了施工安全管理的水平。

[20] 陈文荣.延平区：首次发现国家二级保护植物“毛红椿”群落[J].福建林业,2016,32(6):14.

Chen W R.Yanping areas:first discovered community of Toona ciliata var. pubescens of the national second-grade protection of China[J].Fujian Forestry,2016,32(6):14.

2.1 本工程上人梯共2个，分别位于1#、2#煤仓之间南侧和2#、3#煤仓之间南侧，每相邻两个仓共用1个上人梯。

[21] 陈国科,马克平.生态系统受危险等级的评估标准与方法[J].生物多样性,2012,20(1):66-70.

Chen G K,Ma K P.Criteria and methods for assessing the threat status of ecosystem[J].Biodiversity Science,2012,20(1):66-70.

我踢着街上在小石子，那是什么时候的事？哦，扒锅街消失的第三年，我照旧在新家的街道上被我妈拿着铲子追着，只是再也没有人给我通风报信，通常我都会被揍得很惨。

[22] 宗世贤,陶金川,杨志斌,等.毛红椿的生态地理分布及其南京引种的初步观察[J].植物生态学与地植物学学报,1988,12(3):222-231.

Zong S X,Tao J C,Yang Z B,et al.Ecogeographical distribution of Toona ciliata var. pubescens and preliminary observation on its introduction in Nanjing[J].Acta Phytoecologica et Geobotanica Sinica,1988,12(3):222-231.