价值核算体系融合物质流分析的沱江流域生态系统生产总值分配研究
Study on Distribution of Tuojiang River Basin Gross Ecosystem Product by Value Accounting System Combined with Material Flow Analysis
-
Keywords:
- gross ecosystem product /
- accounting /
- material flow analysis /
- allocation /
- ecological compensation
-
繁育系统在影响植物表型和进化路线方面有着重要作用,主要包括花部形态性状、开放式样、交配系统、传粉者类别和传粉频率等,其中,交配系统可以控制植物种群数量以及种群遗传变异情况,在繁育系统中起着核心作用[1-2]。相关研究认为:植物繁育特性(如开花传粉、交配系统、种子散布方式和进化历史)对其遗传结构、遗传变异及分化和基因流等具有重要影响[3]。植物授粉系统随着环境和授粉条件的变化而变化,其繁育系统决定着植物的分布和扩散能力,也影响着植物的进化潜力。有效传粉是植物果实和种子结实的首要条件,决定着植物的生存潜力,其中有效传粉者起着关键作用[4]。了解繁育系统的生殖特性,对于分析极小种群野生植物致危原因以及生存于恶劣环境导致传粉受限的物种有着重要的作用[5]。
亚洲地区发现的杜鹃属(Rhododendron)植物多达900余种,且一半以上分布在中国;中国74%以上的杜鹃花为特有种,以云南、贵州和四川最为丰富,构成了杜鹃属植物的现代分布和演化中心[6-7]。有学者对极小种群野生植物杜鹃属的物种遗传多样性进行了研究,如宽杯杜鹃(R. sinofalconeri)[8]和红萼杜鹃(R. meddianum)[9]。近年来,杜鹃属植物繁育系统与环境关系的研究备受关注,大部分杜鹃花属植物的繁育系统倾向于异交的复合交配方式[10-14],而宽杯杜鹃[15]和长梗杜鹃(R. longipedicellatum)[16]的繁育系统为自交和异交的混合交配系统,但有关波叶杜鹃繁育系统的研究尚未见报道。
波叶杜鹃(R. hemsleyanum)是杜鹃属常绿灌木,花冠白色宽钟状,边缘波状,其分布海拔低,生长适应性强,有较高的观赏价值,是一种极具开发应用潜力的野生花卉资源。波叶杜鹃仅分布于四川峨眉山,该地处于旅游开发区,其野生资源受人为干扰严重,分布区域狭窄并多伴有人为建筑,生境破碎化,种群数量稀少。在《中国高等植物受威胁物种名录》[17]《中国生物多样性红色名录》[18]《IUCN 红色名录》[19]和《The Red List of Rhododendrons》[20]中,波叶杜鹃被评估为极度濒危物种。经过多次实地考察,课题组发现波叶杜鹃在峨眉山仅有6个分布单元,现存个体数163株,种群更新较弱,是典型的极小种群野生植物。本研究对四川峨眉山的野生波叶杜鹃资源现状进行深入调查,观察其花部形态和昆虫访花行为,结合杂交指数、花粉—胚珠比和人工控制授粉等试验揭示其繁育系统特征,并提出其有效合理的保护措施,为其作为乡土园林树种推广和极小种群杜鹃属植物资源的可持续利用提供理论基础。
1. 材料与方法
1.1 研究区概况
研究区位于四川省峨眉山市峨眉山长老坪(N29º57′16.6′′,E103º36′63.9′′,海拔1 449 m)和石笋沟(N29º55′88.8′′ ,E103º37′11.1′′ ,海拔1 191 m)。该区为中亚热带季风气候,年均温17 ℃,年平均降水量1 555 mm,空气相对湿度达80%以上[21];坡度较陡;土壤多为黄壤和棕壤,少见黄棕壤[22]。波叶杜鹃群落结构复杂,伴生植物乔木层以樟科(Lauraceae)为主,伴有日本柳杉(Cryptomeria japonica)、鞘柄木(Torricellia tiliifolia)、杨叶木姜子(Litsea populifolia)和锐齿槭(Acer sikkimense),灌木层以柃属(Eurya)为主,伴有勾儿茶(Berchemia sinica)和箭竹(Fargesia spathacea)等植物[23]。
1.2 研究材料
供试波叶杜鹃为研究区生长的野生波叶杜鹃。
1.3 研究方法
1.3.1 花期物候与花部特征观察
2021年5—6月对峨眉山长老坪和石笋沟的波叶杜鹃野外居群进行群体花期物候观察。随机挑选波叶杜鹃10株记录其始花、终花以及花期持续时间和开花同步性;在波叶杜鹃花期随机挑选10株,每株均随机挂牌花序3个,每天观察和统计花序新增的开花数量和开放持续时间;在盛花期,观察并统计波叶杜鹃花部形态特征并随机挂牌花朵30朵,记录单花花期时间和形态变化,记录始花、终花时间及花期同步指数[24]。
$ {花}{期} {同}{步}{指}{数}=\dfrac{1}{N}\left(\dfrac{1}{{T}_{{\rm{K}}}}\right)\displaystyle\sum _{{\rm{J=K}}}^{N}{R}_{{\rm{J}}\ne {\rm{K}}}。 $
式中:花期同步指数范围为0~1;N为种群个体总数;TK为单个花期时间;RJ为个体K和个体J的花期重复时间。
1.3.2 花粉活力和柱头可授性检测
波叶杜鹃开花当日花粉活力测定:随机选取当日开花的单花15朵,从9: 00—21: 00每2 h取花粉1次,将其涂抹于载玻片,用5.0%蔗糖与0.5%氯化三苯基四氮唑(TTC)的混合液进行染色,放入培养皿(内置湿润滤纸),25 ℃、无光条件静置2 h,于显微镜下观察染色情况,重复3次。染色花粉数与观察花粉总数比值的百分数为波叶杜鹃的花粉活力值。
波叶杜鹃花期的花粉活力测定:标记花蕾期(主要指开花前1~2 d)和开花后1~8 d的花朵30朵,每天9: 00随机选取单花3朵测其花粉活力,重复3次。测定方法同开花当日花粉活力的测定。
波叶杜鹃开花过程中柱头可授性测定:每天11: 00随机选取花蕾期和开花后1~8 d的单花3朵,共选取30朵,参照DAFNI[25]的方法,将柱头置于凹面载玻片上,滴加反应液(V3%过氧化氢: V水: V 1%联苯胺 = 11∶22∶4)后观察柱头是否染色。根据颜色和气泡判定柱头的可授性,具体标准为:柱头接触反应液后无气泡和颜色出现,为不具有可授性,记为0;柱头接触反应液后产生1~6个气泡且未出现蓝色反应,为部分具有可授性,记为1;柱头1/2面积产生气泡且有微弱的蓝色反应出现,为具有可授性,记为2;柱头2/3面积产生大量气泡且蓝色加深,为可授性较强,记为3;整个柱头几乎都有气泡产生且有明显的蓝色反应,为可授性强,记为4。
1.3.3 估算杂交指数和花粉—胚珠比
于波叶杜鹃盛花期,随机挑选其花序30个(每株单株3个花序,共10株单株),观察其花部特征,参照DAFNI等[26]的标准(花冠直径、花部特征和柱头可授性)计算其杂交指数,并评估波叶杜鹃的繁育系统。
随机挂牌标记已开放但是花药还未开裂的花蕾30个(每株单株采集3个花序,共10株单株),将其采集到室内,待花药开裂后将其制作成花药悬浮液,于体式显微镜下观察并计数其花粉量。摘除波叶杜鹃花瓣、花萼、花冠和雄蕊,将子房在解剖镜下用解剖刀切开,找出胚珠,记录胚珠数。参照CRUDEN[27]的标准评估波叶杜鹃的繁育系统,即花粉—胚珠比为244.7~2 588.0时,繁育系统为兼性异交,花粉—胚珠比越小则近交程度越强,花粉—胚珠比越大则远交程度越强。
1.3.4 授粉方式研究
选取峨眉山长老坪的波叶杜鹃开花植株,根据TAMURA等[28]的方法对其进行6种不同授粉方式的处理:(1)未开花前去雄蕊,套硫酸纸袋处理(孤雌生殖);(2)花蕾期去雄蕊,不套任何袋(传粉媒介介导的异花授粉);(3)整个单花花期套网袋(自花授粉);(4)花蕾期去雄蕊并套硫酸纸袋,开花到第3 天时用同株但不同花的花粉授粉并套硫酸纸袋(同株异花授粉);(5)花蕾期去雄套袋,第3 天使用间隔一定距离(10 m)且单花开放的植株花粉授粉,授粉完成后套袋(异株异花授粉);(6)自然对照。所有处理挂牌标记,种子成熟时统计结实率。结实率=成熟果实数/处理花朵数×100%。
1.3.5 访花昆虫及其行为观察
2021年5月中旬(波叶杜鹃盛花期)进行连续3 d访花昆虫观察,记录访花昆虫种类、数量和行为等。用捕虫网收集访花昆虫并置于75%酒精中保存,用于后续鉴定。根据STOUT[29]判定有效访花昆虫的标准(即具有携带花粉的能力又能触碰柱头的昆虫),明确波叶杜鹃的有效传粉者。
2. 结果与分析
2.1 花期物候与花部特征
波叶杜鹃为顶生总状伞形花序,每个花序花朵数7~16朵,花朵展示醒目。波叶杜鹃5月初始花出现,5月中旬为盛花期, 6月初花期结束,单花开花6~10 d,单花开放进程可划分为花蕾期、盛花期和凋谢期(图1)。单株开花同步性较高,花期同步指数为0.72±0.31,呈多峰曲线式开花进程,开花数量于始花出现后第6天达到峰值(图2)。
波叶杜鹃花冠呈白色宽钟状,直径30.17~76.47 mm,7裂,裂片扁圆形,顶端有缺刻;雄蕊不等长,平均14~16 mm,花药白色长椭圆形。波叶杜鹃花部主要形态指标见表1。花蕾期时,波叶杜鹃柱头被花瓣包裹,花朵顶部淡红,剖开花朵后发现柱头白色未成熟,柱头未分泌黏液;盛花期时,花瓣变为白色,完全打开,柱头弯曲呈拱形,有黏液,雄蕊和雌蕊分离;凋谢期时,花瓣边缘卷曲呈褐色,柱头颜色逐渐变为褐色,有黏液。
表 1 波叶杜鹃花部特征Table 1. Floral syndrome of Rhododendron hemsleyanum指标
indexes最小值
minimum最大值
maximum平均值
mean标准差
standard deviation变异系数/%
coefficient of variation花瓣长/mm petal length 16.11 28.57 23.54 2.59 11.00 花瓣宽 /mm petal width 15.59 24.01 18.11 3.14 17.33 花冠直径/mm flower diameter 30.17 76.47 43.05 4.11 9.55 柱头宽 /mm stigma diameter 3.01 4.79 4.34 0.92 21.20 花柱长/mm style length 39.34 52.07 45.29 3.79 8.37 花梗长/mm pedicel length 26.55 46.59 34.31 4.17 12.15 花丝长/mm filament length 19.17 28.03 24.59 4.77 19.40 花萼长/mm calyx length 1.12 1.76 1.47 0.74 50.34 花药长/mm anther length 2.87 4.52 3.62 1.01 27.90 子房长/mm ovary length 6.21 8.27 7.14 3.14 44.00 子房宽/mm ovary width 4.14 7.69 5.29 1.02 19.28 柱头与花药最短距离/mm
shortest distance of stigma to anthers16.74 41.20 28.43 3.57 12.56 2.2 花粉活力及柱头可授性
由图3可知:开花当天,波叶杜鹃的花粉活力在11: 00最高(90.25%),之后呈逐渐降低趋势,17: 00时最低(48.67%),之后缓慢上升;花蕾期的花粉活力为87.11%,随着花朵开放和花药开裂,花粉活力上升,花粉活力在开花第1天最大(91.05%),之后花粉活力逐日下降,花药和柱头颜色逐渐变黑,随着花粉散出其活力逐渐下降,至凋谢期时(花期第9~10天)花粉活力仍保持在约30%,可见,自然状态下波叶杜鹃能长期保持较高的花粉活力。此外,在花期第1和2天时,波叶杜鹃柱头不产生黏液,未出现呈蓝色的颜色反应,不具有可授性;在开花后3~7 d,柱头与反应液反应逐渐剧烈,蓝色逐渐明显,气泡增多;在第6和7 天柱头与反应液反应剧烈,柱头产生气泡最多,可授性最强;在第8和9天柱头与反应液反应逐渐减弱,蓝色变浅的同时气泡变少;开花第10天,野外观察到波叶杜鹃花瓣和雄蕊开始脱落,柱头可授性结果显示无颜色反应和未产生气泡,此时波叶杜鹃失去柱头可授性。
2.3 杂交指数与花粉—胚珠比
按照DAFNI等 [26]的标准估算波叶杜鹃杂交指数。波叶杜鹃花冠直径为(43.05±4.11) mm>6 mm,记为3;其同时具有雌蕊和雄蕊,其中雄蕊先熟,记为1;开花时柱头弯曲且位置高于花药,产生空间分离现象,记为1。因此,波叶杜鹃杂交指数为3+1+1=5,可推断波叶杜鹃繁育系统部分自交亲和。
波叶杜鹃单花花粉数量为416 289±5 125.71,单花胚珠数量为649.51±67.45,因此,其花粉—胚珠比为640.93±65.52,根据CRUDEN[27]对繁育系统的评判标准,判断波叶杜鹃繁育系统为兼性异交。
2.4 不同授粉方式对波叶杜鹃结实率影响
由表2可知:自然对照组果实17个,结实率56.67%,说明波叶杜鹃能在峨眉山野生居群下自然结实,但结实率较差;去雄套袋果实数和结实率均为0,说明波叶杜鹃不存在孤雌生殖现象;去雄不套袋果实数和结实率均略低于自然对照组,说明波叶杜鹃可以在虫媒帮助下授粉结实,且基本能达到自然条件下的结实率;自花授粉能产生2个果实,结实率为6.67%,判断波叶杜鹃能通过自交结实,具有一定程度的自交亲和性;异株异花授粉和同株异花授粉的果实数和结实率均明显高于自然对照组,进一步说明波叶杜鹃异交传粉机制有效,主要以异交为主,其结实受传粉者限制。
表 2 不同授粉方式处理后波叶杜鹃的结实率Table 2. Fruit set of R. hemsleyanum under different pollination treatments授粉方式
pollination treatments花序数量
number of inflorescences授粉花数量
number of flowers果实数量
number of fruits结实率/%
fruit set自然对照 natural pollination 10 30 17 56.67 去雄套袋 emasculated and bagged 10 30 0 0 去雄不套袋 emasculation is not bagged 10 30 15 50.00 自花授粉 autonomous self-fertilization 10 30 2 6.67 同株异花授粉 strains cross-pollination 10 30 23 76.67 异株异花授粉 different strains cross-pollination 10 30 26 86.67 2.5 访花昆虫及其行为
波叶杜鹃访花昆虫主要包括蜜蜂科(Apidae)、蚁科(Formicidae)、食蚜蝇科(Syrphidae)和蚜总科(Aphidoidea) (表3)。
表 3 波叶杜鹃访花昆虫Table 3. Flower visiting insects of R. hemsleyanum目
order科
family种
species是否为有效传粉者
effective pollinator or not膜翅目 Hymenoptera 蜜蜂科 Apidae 中华蜜蜂 A. cerana 是 yes 排蜂 A. dorsata 是 yes 蚁科 Formicidae 蚁科某种 Formicidae sp. 否 no 双翅目 Diptera 食蚜蝇科 Syrphidae 灰带管蚜蝇 Eristalis cerealis 否 no 黑带食蚜蝇 Episyphus balteatus 否 no 同翅目 Homoptera 蚜总科 Aphidoidea 蚜总科某种 Aphidoidea sp. 否 no 中华蜜蜂访花前先在空中盘旋,之后停留在某一植株上访花,落到波叶杜鹃花冠上时头部伸入花腔内,通过其细长的口器吸食花蜜或花粉,转动身体,后足更换着力点,行动较为缓慢,在同一朵花上停留数分钟,在此过程中身体沾到大量花粉,能有效地为波叶杜鹃传粉(图4a~c)。排蜂靠近波叶杜鹃花朵时,悬空震动翅膀,前足先接触花冠,然后采粉,随之晃动使得花粉粘在其头部和身体上,用其足部刮刷,将花粉聚集到其花粉筐中(图4d~f);当排蜂进入另一朵花时,就将其携带的花粉授到柱头上,能有效地为波叶杜鹃传粉。而黑带食蚜蝇(Episyphus balteatus) 、蓟马(Thrips)和蚂蚁等小型的昆虫主要以盗取花蜜为主,身体各部均不携带花粉且不接触柱头。因此,判定波叶杜鹃白天的有效传粉者是中华蜜蜂和排蜂。此外,观察发现传粉者的传粉活动可能受自然因素的影响较大,波叶杜鹃开花时期,温度较高且天气晴朗的时候访花者频率较高,主要是由于波叶杜鹃花朵展示醒目并散发芳香味道,吸引传粉者;而盛花期时阴雨天较多,空气湿度低,温度较低,影响访花者活动频率。
3. 讨论
3.1 波叶杜鹃花粉活力及花部特征
本研究中,波叶杜鹃花冠为大型花冠,开花呈伞形花序,开花时花朵白色,花香浓郁,花瓣不折叠,开花量较大且同步较高,具有蜜腺,增加花展示面积的同时可吸引中华蜜蜂和排蜂传粉。开花式样通过吸引传粉者影响花粉散布方式,进而影响植物的交配系统[30]。波叶杜鹃雄蕊和雌蕊成熟时期不同,根据观察发现其雄蕊先熟;其花粉活力在开花当天11: 00达到峰值(90.27%),低于宽杯杜鹃 (花粉活力最高达92.80%)[15]和长梗杜鹃 (花粉活力最高达92.18%)[31],高于锈叶杜鹃 (R. siderophyllum,花粉活力最高达83.95%)[13]和大喇叭杜鹃 (R. excellens,花粉活力最高达84.43%)[14]。波叶杜鹃花粉活力保持在50%以上的时间有7 d,花朵开放时期的花粉活力均在40%以上,推断波叶杜鹃具有大量花粉及高效保持花粉活力的机制以保证其传粉成功。凋谢期时,波叶杜鹃的花粉活力仍可保持在29%以上,推断这可能是由于盛花期未成熟花粉逐渐成熟并累积到凋谢期所致。
本研究表明:波叶杜鹃柱头可授性最强时,花粉活力仍保持在40%以上;野外观察发现:雄蕊和雌蕊有一段短暂的同时期表达时段,但该时期波叶杜鹃花柱弯曲呈拱形,花柱悬垂,导致雄蕊和雌蕊产生空间异位现象,根据自花授粉结果,推断雌雄异位现象有效避免了自花传粉的可能性。两性植物保证异交成功的花部特征是花朵雌雄异位现象,可以避免自交和有性功能干扰的机制[32-33];同时波叶杜鹃花柱的悬垂形状还有效避免了雨水对柱头花粉和表面黏液的冲刷,这与前人研究结果[34-35]一致。开展遗传育种时,作为父本的波叶杜鹃花粉活力在开花当天11: 00达到最高值,此时采集花粉进行人工授粉试验可能最为有利。
3.2 波叶杜鹃繁育系统和授粉方式
根据DAFNI等[26]和CRUDEN[27]的划分标准,波叶杜鹃繁育系统类型为兼性异交,且自然条件下其结实率达56.67%。然而,在野生居群中极少发现波叶杜鹃实生苗,说明该物种存在自然更新障碍,这可能与其生境土层浅薄和郁闭度高等立地条件导致的种子萌发困难有关。FREE[36]和READER[37]研究表明大多数杜鹃花科植物都具有自交亲和性,波叶杜鹃自花授粉结实率6.67%,判断其能通过自交结实。野外调查发现:波叶杜鹃在开花第3~7天时花药逐渐散落且花粉呈线状;开花第8~10天,花冠开始枯萎下垂,花托开始与基部分离,导致花瓣和雄蕊向花柱偏移,花粉黏附在花柱上实现部分自交,故判断波叶杜鹃具有一定程度的自交亲和性,与马雄杜鹃 (R. maxiongense)[12]、锈叶杜鹃[13]和大喇叭杜鹃[14]的繁育系统一致。本研究推断波叶杜鹃繁育系统由自交向异交演化的可能性较大,因为波叶杜鹃同长梗杜鹃[31]、牛皮杜鹃 (R. aureum)[38]等大部分杜鹃花一样进化出柱头呈拱形悬垂状的特征,通过雄蕊和雌蕊空间异位避免自花授粉现象。波叶杜鹃分布在陡壁悬崖边和沟谷地,种间距离较远,传粉媒介较为缺乏,异交受到限制,部分自交亲和可能是为适应群落特征而产生的补偿作用。频繁的远交和近交都会导致其繁育系统适合度下降,因此,最适合的交配方式是实现适度的杂交[39],波叶杜鹃的花药和柱头在空间上异位,这将大大降低自花授粉,防止自交衰退,增强子代的适应能力。本研究认为波叶杜鹃部分自交亲和,这与LEE等[40]的观点一致,即:起源于较小奠基者种群的物种,在分布区域狭窄、生境破碎化和种群数量较小时,为了种群繁育经常会发生近交现象。自交物种的环境适应能力及形成新种的潜力较弱,自交不适宜作为长期生殖策略[41]。将波叶杜鹃选作母本时,其柱头可授性在花朵开花后第6和7天最强,此时通过人工异花授粉的方式可提高波叶杜鹃自然条件下异交的概率,促进不同种群间的基因交流。
3.3 访花昆虫
国内已报道大喇叭杜鹃[14]、长梗杜鹃[16]和大叶金顶杜鹃(R. prattii Franch.)[42]的有效传粉昆虫分别为中华蜜蜂、熊蜂和食蚜蝇类。本研究中,波叶杜鹃白天的主要访花者为中华蜜蜂和排蜂,且中华蜜蜂访问对象偏向初开期花朵,此时花粉活力强且散出量大,能保证传粉的时效性。蜜蜂访问频次可能与其耐热性、耐寒性、抗逆性和花粉传递效率等相关[43]。日间蝴蝶并没有访问过波叶杜鹃,这种差异可能是由于生境的环境气候条件所致[44]。据观察,蓟马和蚂蚁等小型的昆虫主要盗取花蜜,根据BASNETT等[39]提出“盗蜜者和传粉者授粉能力可能相同”的观点,推测蓟马和蚂蚁盗取花蜜的行为可能有利于波叶杜鹃花部特征的进化。HUANG等[45]研究表明:白色芳香杜鹃花的日间传粉者主要是蜜蜂和蝴蝶,夜间则会吸引蛾类授粉者。然而,由于夜间环境条件恶劣加之分布地紧邻悬崖峭壁,未开展夜间传粉观察,未来有条件需开展波叶杜鹃夜间有效传粉者的观察研究。在通过人工繁育等方法扩大波叶杜鹃群体数量的同时,应在花期结合释放有效访花昆虫的方式,全面提高波叶杜鹃在自然条件下的异交结实率。
4. 结论
波叶杜鹃繁育系统是倾向于异交的混合交配系统,自然条件下有效传粉者为中华蜜蜂和排蜂。波叶杜鹃花部特征表现出对虫媒传粉的适应性,在单株和群体水平上波叶杜鹃具有较高的开花同步性,可吸引较多访花者,产生较高的结实率,其雄蕊和雌蕊的空间异位可在一定程度上避免自交。波叶杜鹃野生分布地处于旅游开发区,受人为干扰程度较强,可通过人工繁育等方式扩大其种群,也可在花期释放有效昆虫等提高其自然条件下异交的概率。
-
图 1 研究区位置示意图
注:四川省行政地图审图号:川S(2021)00050号,四川省标准地图水系版审图号:川S(2021)00068号;底图无修改。
Figure 1. Location diagram of study area
Note: Approval number of Sichuan administrative map is Sichuan S (2021) 00050, approval number of water system regarding Sichuan standard map is Sichuan S (2021) 00068; the base map has not been modified.
表 1 流域生态系统生产总值核算方法
Table 1 Accounting method of gross ecosystem product of river basin ecosystem
一级指标
the 1st level indicator二级指标
the 2nd level indicator核算公式
accounting formula公式说明
formula description物质产品
material products农业、林业、畜牧业和渔业产品;水资源
products of agriculture, forestry, animal husbandry and fishery; water resources$V_{\rm{M}}=V_{{\rm{G D P}}} \times \dfrac{S_{\text {流域 }}}{S} $ VM为物质产品价值量,亿元;VGDP为该区域GDP值,亿元;S流域为该区域内流域面积,km2;S为该区域行政面积,km2。 调节服务
regulation services水源涵养
water conservation$\begin{aligned} & Q_{{\rm{w r}}}=\sum_{i=1}^n\left(P_i-R_i-E T_i\right) \times {A_i} \\& V_{\mathrm{wr} }=Q_{\mathrm{wr} } \times C_{\mathrm{we} }\end{aligned}$ Qwr为总水源涵养量,m3/a;Pi为产流降雨量, mm/a;Ri为地表径流量,mm/a;ETi为蒸散发量,mm/a;Ai为i类生态系统面积,m2;i为研究区第i类生态系统类型;n为研究区生态系统类型;Vwr为水源涵养总价值量,m3/a;Cwe为水资源交易市场价格,元/m3。 土壤保持
soil conservationQsr=R×K×L×S×(1−C1×P)
Vsr=Vsd+Vdpd
Vsd=λ×(Qsr/ρ) ×C2$V_{{\rm{d p d}}}=\displaystyle\sum_{i=1}^n Q_{{\rm{s r}}} \times C_i \times P_i $ Qsr为土壤保持量;R为降雨侵蚀力因子;K为土壤可蚀性因子;L为坡长因子;S为坡度因子;C1为植被覆盖和管理因子;P为水土保持措施因子;Vsr为生态系统土壤保持价值,元/a;Vsd为减少泥沙淤积价值,元/a;Vdpd为减少面源污染价值,元/a;C2为单位水库清淤工程费用,元/m3;ρ为土壤容重,t/m3;λ为泥沙淤积系数;i为土壤中氮磷等营养物质数量,i=1,2,...,n;Ci为土壤中氮磷等营养物质的纯含量,%;Pi为处理成本,元/m3。 洪水调蓄
flood regulation
and storageCfm=Crc+Clc
Vfm=Cfm×CweCfm为洪水调蓄量,m3/a;Crc为库塘洪水调蓄量, m3/a;Clc为湖区洪水调蓄量,m3/a;Vfm为生态系统洪水调蓄价值,元/a;Cwe为水库单位库容的建设成本,元/m3。 碳固定
carbon fixation$Q_{\mathrm{CSO}_2}=\displaystyle\sum_{i=1}^n A_i \times S_ { i}$
Vcf= QCO2×CcQCSO2为土壤固碳量,t/a;Ai为不同系统中的土壤面积,hm2;Si为不同系统实测土壤固碳量,t/(hm2∙a); Vcf为生态系统固碳价值,元/a;QCO2为生态系统固碳总量,t/a;Cc为碳价格,元/t。 空气净化
air purificationAPQ= QSO2+QNOx+Cdusts
APV=QSO2×CSO2+QNOx×CNOx+Qdusts×CdustsAPQ为大气污染物排放总量,t;QSO2为二氧化硫排放量,t;QNOx为氮氧化物排放量,t;Qdusts为颗粒物排放量,t;APV为空气净化价值,元;CSO2为二氧化硫治理成本,元/t;CNOx为氮氧化物治理成本,元/t;Cdusts为颗粒物治理成本,元/t。 水质净化
water purification$V_{\rm{w}}=\displaystyle\sum_{i=1}^3 c_i \times Q_i$ Vw为生态系统水质净化的价值,元/a;ci为治理第i种污染物的成本,元/t;Qi为第i种污染物的排放量,t。 气候调节
climate regulationEwe=Ew×q×103/3600+Ew×y Ewe为生态系统蒸发消耗的能量,kW/h;Ew为年蒸散量,mm;q为挥发潜热,即蒸发1 g水所需热量,J; y为加湿器将1 m3水转化为蒸汽的耗电量,kW/h。 文化服务
cultural services自然景观
natural landscape$ N_t=\displaystyle\sum_{i=1}^n N_{t i}$ $ V_r=\displaystyle\sum_{j=1}^J N_j \times T C_j$
TCj=Tj×Wj+Cj
Cj=Ctc,j+Clf,j+Cef,jNt为旅客总人数;Nti为第i个旅游区的人数;n为旅游区个数,i=1,2,3...n;Vr为被核算地点的休闲旅游价值,元/a;Nj为每年j地到核算地区旅游的总人数;j为被核算地点旅游的游客所在区域;TCj为来自j地的游客的平均旅行成本,元/人;Tj为来自j地的游客用于旅游地点的平均时间,d/人;Wj为来自j地的游客的当地平均日工资,元/人;Cj为来自j地的游客花费的平均直接旅行费用,元/人;Ctc,j为游客从j地到核算区域的交通费用,元/人;Clf,j为食宿花费,元/人;Cef,j为门票费用,元/人。 -
[1] 靳乐山, 刘晋宏, 孔德帅. 将 GEP 纳入生态补偿绩效考核评估分析[J]. 生态学报, 2019, 39(1): 24. DOI: 10.5846/stxb201809252081. [2] HUANG L, HE C L, WANG B. Study on the spatial changes concerning ecosystem services value in Lhasa River Basin, China[J]. Environmental Science and Pollution Research, 2022, 29(5): 7827. DOI: 10.1007/s11356-021-16245-8.
[3] ZHANG Z M, GAO J F, FAN X Y, et al. Response of ecosystem services to socioeconomic development in the Yangtze River Basin, China[J]. Ecological Indicators, 2017, 72: 481. DOI: 10.1016/j.ecolind.2016.08.035.
[4] YU H J, XIE W, YANG L, et al. From payments for ecosystem services to eco-compensation: conceptual change or paradigm shift?[J]. Science of the Total Environment, 2020, 700: 134627. DOI: 10.1016/j.scitotenv.2019.134627.
[5] 陈默, 林育青, 张建云, 等. 水生态系统生产总值核算体系及应用[J]. 水资源保护, 2023, 39(1): 234. DOI: 10.3880/j.issn.1004-6933.2023.01.030. [6] 欧阳志云, 朱春全, 杨广斌, 等. 生态系统生产总值核算: 概念、核算方法与案例研究[J]. 生态学报, 2013, 33(21): 6747. DOI: 10.5846/stxb201310092428. [7] 石敏俊, 陈岭楠. GEP核算: 理论内涵与现实挑战[J]. 中国环境管理, 2022, 14(2): 5. DOI: 10.16868/j.cnki.1674-6252.2022.02.005. [8] COSTANZA R, D'ARGE R, GROOT R, et al. The value of the world’s ecosystem services and natural capital[J]. Nature, 1997, 387(6630): 253. DOI: 10.1038/387253a0.
[9] 谢高地, 鲁春霞, 成升魁. 全球生态系统服务价值评估研究进展[J]. 资源科学, 2001, 23(6): 5. DOI: 10.3321/j.issn:1007-7588.2001.06.002. [10] 邓娇娇, 常璐, 张月, 等. 福州市生态系统生产总值核算[J]. 应用生态学报, 2021, 32(11): 3835. DOI: 10.13287/j.1001-9332.202111.018. [11] 罗芳, 潘安, 陈忠升, 等. 四川省土地利用变化对生态系统服务价值的影响研究[J]. 云南农业大学学报(自然科学), 2021, 36(4): 734. DOI: 10.12101/j.issn.1004-390X(n).202010031. [12] 宋昌素, 欧阳志云. 面向生态效益评估的生态系统生产总值 GEP 核算研究: 以青海省为例[J]. 生态学报, 2020, 40(10): 3207. DOI: 10.5846/stxb202004260999. [13] 牛乐, 王娇月, 郗凤明, 等. 基于生态系统服务付费模型的福州市区域生态补偿额度核算[J]. 应用生态学报, 2021, 32(11): 3805. DOI: 10.13287/j.1001-9332.202111.009. [14] 白玛, 卓嘎, 肖燚, 等. 基于生态系统生产总值核算的习水县生态保护成效评估[J]. 生态学报, 2020, 40(2): 499. DOI: 10.5846/stxb201812202767. [15] 牟雪洁, 王夏晖, 张箫, 等. 北京市延庆区生态系统生产总值核算及空间化[J]. 水土保持研究, 2020, 27(1): 265. DOI: 10.13869/j.cnki.rswc.2020.01.037. [16] 陈梅, 纪荣婷, 刘溪, 等. “两山”基地生态系统生产总值核算与“两山”转化分析: 以浙江省宁海县为例[J]. 生态学报, 2021, 41(14): 5899. DOI: 10.5846/stxb202008032025. [17] 高玉娟, 张莹. 国有林区生态生产总值GEP核算及生态贡献度研究[J]. 林业经济问题, 2020, 40(2): 173. DOI: 10.16832/j.cnki.1005-9709.2020.02.008. [18] 陈宗铸, 雷金睿, 吴庭天, 等. 国家公园生态系统生产总值核算: 以海南热带雨林国家公园为例[J]. 应用生态学报, 2021, 32(11): 3883. DOI: 10.13287/j.1001-9332.202111.010. [19] OUYANG Z Y, SONG C S, ZHENG H, et al. Using gross ecosystem product (GEP) to value nature in decision making[J]. Proceedings of the National Academy of Sciences, 2020, 117(25): 14593. DOI: 10.1073/pnas.1911439117.
[20] 赵筱青, 苗培培, 普军伟, 等. 抚仙湖流域土地利用变化及其生态系统生产总值影响[J]. 水土保持研究, 2020, 27(2): 291. DOI: 10.13869/j.cnki.rswc.2020.02.041. [21] PASTOR A V, TZORAKI O, BRUNO D, et al. Rethinking ecosystem service indicators for their application to intermittent rivers[J]. Ecological Indicators, 2022, 137: 108693. DOI: 10.1016/j.ecolind.2022.108693.
[22] 刘家倩, 陆建忠, 田晴, 等. 干旱对鄱阳湖流域植被生态系统 GEP 的影响[J]. 华中师范大学学报(自然科学版), 2022, 56(3): 513. DOI: 10.19603/j.cnki.1000-1190.2022.03.019. [23] 马元希, 于德永, 王欣烨, 等. 黄河青海流域生态系统服务价值评估方法研究及时空变化分析[J]. 青海民族大学学报(社会科学版), 2022, 48(2): 70. DOI: 10.3969/j.issn.1674-9227.2022.02.009. [24] 耿翔燕, 葛颜祥, 张化楠. 基于重置成本的流域生态补偿标准研究: 以小清河流域为例[J]. 中国人口资源与环境, 2018, 28(1): 140. DOI: 10.12062/cpre.20170720. [25] 姚婧, 肖宇婷, 樊敏, 等. 沱江流域主要污染负荷预测及时空分布特征[J]. 中国环境监测, 2022, 38(3): 105. DOI: 10.19316/j.issn.1002-6002.2022.03.12. [26] 吴之见, 杜思敏, 黄云, 等. 基于生态系统生产总值核算的生态保护成效评估: 以赣南地区为例[J]. 生态学报, 2022, 42(16): 1. DOI: 10.5846/stxb202104181007. [27] 樊敏, 肖宇婷, 姚婧, 等. 沱江流域污染负荷与生态系统服务价值空间关系解析研究[J]. 长江流域资源与环境, 2022, 31(6): 1354. DOI: 10.11870/cjlyzyyhj202206016. [28] 彭立, 邓伟, 黄佩, 等. 四川盆地多重生态系统服务景观指数评价与服务簇识别[J]. 生态学报, 2021, 41(23): 9328. DOI: 10.5846/stxb202009062319. [29] 陈效逑, 乔立佳. 中国经济—环境系统的物质流分析[J]. 自然资源学报, 2000, 15(1): 17. DOI: 10.11849/zrzyxb.2000.01.003. [30] 邓伟. 环境保护税税额地方确定权的实施现状与制度完善[J]. 税务研究, 2022(6): 33. DOI: 10.19376/j.cnki.cn11-1011/f.2022.06.021. [31] 李倩娜, 唐洪松, 胡艳. 沱江流域生态补偿模式演变、实践困境及优化对策[J]. 环境保护, 2022, 50(19): 24. DOI: 10.14026/j.cnki.0253-9705.2022.19.003. [32] 杨亮洁, 杨晓蓉, 杨永春. 城市内生竞争力与外生竞争力耦合协调研究: 以成渝城市群为例[J]. 人文地理, 2021, 36(6): 76. DOI: 10.13959/j.issn.1003-2398.2021.06.009. [33] 王问苈, 时培豪, 黄庆华. 成渝地区双城经济圈绿色创新效率、演变趋势与影响因素研究[J]. 西南大学学报(自然科学版), 2022, 44(3): 172. DOI: 10.13718/j.cnki.xdzk.2022.03.020. [34] 郑小强, 蒲泱竹. 成渝双城经济圈产业结构升级与环境效率: 基于非动态面板门槛模型[J]. 软科学, 2021, 35(11): 58. DOI: 10.13956/j.ss.1001-8409.2021.11.10. [35] 朱靖, 余玉冰, 王淑. 岷沱江流域水环境治理绩效综合评价方法研究[J]. 长江流域资源与环境, 2020, 29(9): 1995. DOI: 10.11870/cjlyzyyhj202009010. [36] 张玲玲, 车力, 王宗志, 等. 长江经济带水污染与经济脱钩关系演变及驱动因素研究[J]. 水利水电技术, 2021, 52(12): 47. DOI: 10.13928/j.cnki.wrahe.2021.12.005. [37] 温泉, 马迎群, 秦延文, 等. 成渝地区中小河流水生态环境保护存在的问题与对策[J]. 环境工程技术学报, 2022, 12(2): 493. DOI: 10.12153/j.issn.1674-991X.20210639. [38] 蔡绍洪, 谷城, 张再杰. 长江经济带绿色发展水平测度及时空演化特征[J]. 华东经济管理, 2021, 35(11): 25. DOI: 10.19629/j.cnki.34-1014/f.210519005. [39] ALHARBI O M L, BASHEER A A, KHATTAB R A. Health and environmental effects of persistent organic pollutants[J]. Journal of Molecular Liquids, 2018, 263: 442. DOI: 10.1016/j.molliq.2018.05.029.
[40] 邓元杰, 侯孟阳, 贾磊, 等. 基于生态系统服务价值评估的长征沿线革命老区生态补偿策略[J]. 应用生态学报, 2022, 33(1): 159. DOI: 10.13287/j.1001-9332.202201.019. [41] 陈影, 文传浩, 沈体雁. 成渝地区双城经济圈绿色发展效率评价及时空演变研究[J]. 长江流域资源与环境, 2022, 31(5): 1137. DOI: 10.11870/cjlyzyyhj202205017. -
期刊类型引用(6)
1. 侯秋梅,周艳,杨朔,张佳佳,周洪英,吴洪娥. 不同香水月季染色体倍性鉴定及其繁育系统分析. 西北植物学报. 2024(04): 654-661 . 百度学术
2. 熊文艳,吴雅文,白天,解玮佳,董燕,张敬丽. 露珠杜鹃繁育系统及杂交亲和性. 植物研究. 2024(03): 370-379 . 百度学术
3. 李薇,邱利文,王毅敏,彭刚志,高晗,李争艳,刘磊. 疏花水柏枝残存居群开花特性及繁育系统. 植物科学学报. 2024(03): 296-303 . 百度学术
4. 柳文,杨琳祥,木正春,马永鹏,谢孟,万有名,李正红,马宏. 极小种群野生植物波叶杜鹃组培快繁体系的构建. 云南农业大学学报(自然科学). 2024(04): 99-109 . 百度学术
5. 俞群,阮淑明,施清华,黄晓春,曾华浩,高伟. 珍稀植物银粉蔷薇开花生物学特性与繁育系统. 西北植物学报. 2024(12): 1937-1945 . 百度学术
6. 周宏香,邵凤侠,张旻桓,廖菊阳,文亚峰,刘艳,张峥,李雨彤,廉欣欣,刘鲲鹏. 三种杜鹃开花生物学特性研究. 植物科学学报. 2024(06): 703-716 . 百度学术
其他类型引用(1)