Spatiotemporal heterogeneity and driving factors of ecosystem service value in the Yellow River Basin
-
摘要:目的
在黄河流域生态恢复背景下，研究生态系统服务价值（ESV）的时空动态演变特征及其驱动因素，为可持续发展提供理论参考。
方法基于2001—2020年20期的黄河流域土地利用数据，计算黄河流域ESV，利用空间自相关和探索性时空数据分析，揭示黄河流域ESV的时空分异动态特征，并引入时空地理加权回归模型探究黄河流域ESV时空分异的驱动因素。
结果（1）2001—2020年黄河流域ESV整体呈上升趋势，生态系统服务中调节服务价值所占比重最大。不同土地利用类型中，草地和林地所提供的生态服务价值最高。（2）黄河流域ESV空间上存在显著正相关。大多数地区，局部空间自相关指标时间路径的相对长度和弯曲度较小，结构比较稳定，但空间依赖效应较弱。ESV移动方向中负向协同增长占比最大，呈低速增长特征。2001—2020年ESV时空凝聚概率超过80%，说明黄河流域ESV局部空间格局呈现相对锁定状态。（3）流域ESV驱动因素中，年降雨量主要起负向作用，年均气温和植被覆盖则主要起正向作用，而人口密度与夜间灯光指数的影响作用相对较弱。
结论黄河流域ESV增长现状相对稳定，短期内较难改变，在影响ESV变化的因素中自然因素起主导作用。
Abstract:ObjectiveThis paper aims to study the spatiotemporal dynamic evolution characteristics and driving factors of ecosystem service value (ESV) in the context of ecological restoration in the Yellow River Basin, and provide theoretical reference for sustainable development.
MethodBased on the land use data of the Yellow River Basin from 2001 to 2020, the ESV of the Yellow River Basin was calculated. Spatial autocorrelation and exploratory spatiotemporal data analysis (ESTDA) were used to reveal the dynamic characteristics of spatiotemporal differentiation of ESV in the Yellow River Basin. The spatiotemporal geographic weighted regression (GTWR) model was introduced to explore the driving factors of spatiotemporal differentiation of ESV in the Yellow River Basin.
Result(1) From 2001 to 2020, the ESV of the Yellow River Basin showed an overall upward trend, and the value of regulation services accounted for the largest proportion of ecosystem services. Among different land use types, grassland and forest land provided the highest ecological service value. (2) There was a significant positive correlation in the spatial distribution of ESV in the Yellow River Basin. In most regions, the relative length and curvature of the local spatial autocorrelation index time path were relatively small, and the structure was relatively stable, but the spatial dependence effect was weak. The proportion of negative synergistic growth in the direction of ESV movement was the highest, showing a low-speed growth characteristics. The probability of spatiotemporal aggregation of ESV from 2001 to 2020 exceeded 80%, indicating that the local spatial pattern of ESV in the Yellow River Basin was in a relatively locked state. (3) Among the driving factors of watershed ESV, annual rainfall mainly played a negative role, while annual average temperature and vegetation cover mainly played a positive role, while the impact of population density and nighttime light index was relatively weak.
ConclusionThe current situation of ESV growth in the Yellow River Basin is relatively stable and difficult to change in the short term. Natural factors play a dominant role in influencing ESV changes.
-
随着全球气候变暖，极端气候频现，干旱、盐害等非生物胁迫已成为威胁全球农林业生产和发展的一个重大挑战[1]。植物为了适应逆境环境并得以生存，进化出相关基因来维持细胞内结构和功能机制的稳定[2]。迄今为止，已发现许多与干旱和盐胁迫防御相关的功能基因，其中包括膜联蛋白基因。膜联蛋白是一个超级基因家族的多功能脂结合蛋白，到目前已从超过65个物种中分离和克隆了编码基因，包括真菌、原生生物、植物、高等脊椎动物和原核生物 [3]。Creutz等人[4]1978年首次从牛肾上腺髓质嗜铬细胞中分离纯化出一种蛋白质Annexin A7（synexin），这种蛋白促使嗜铬颗粒通过依赖Ca2+来结合一些膜或者内膜成分。而国内学者对植物膜联蛋白的研究较晚，第1个植物膜联蛋白是从悬浮培养的番茄 (Solanum lycopersicum) 细胞中发现和分离的[5]，且与动物膜联蛋白在序列上有较高的同源性[6]。随着对植物膜联蛋白的进一步了解，发现膜联蛋白广泛分布于各种植物中，如棉花（Gossypium spp.）[7]、烟草（Nicotiana tabacum）[8]、玉米（Zea mays）[9]、拟南芥（Arabidopsis thaliana）[10]、水稻（Oryza sativa）[11]等。不同植物中的膜联蛋白以基因家族的形式存在，其家族成员各自具有不同的功能和时空表达模式，表达具有组织特异性。据报道膜联蛋白参与重要的生物过程，如膜运输、细胞骨架组织、细胞内稳态和离子运输，具有Ca2+结合、离子通道和过氧化物酶活性等不同的生物学功能，在植物生长发育、信号转导、适应环境和响应非生物逆境胁迫中发挥着重要作用[12]。在渗透胁迫、脱落酸（ABA）或干旱等非生物胁迫时紫花苜蓿（Medicago sativa）中的膜联蛋白被激活上调表达[13]。拟南芥中膜联蛋白的表达能响应各种非生物胁迫（如盐胁迫、渗透胁迫和低温胁迫）而发生变化[14]。Cantero等人[15]评估了拟南芥中的8种膜联蛋白基因表达模式，并建立了它们在盐、干旱和其他非生物胁迫下的差异调控机制。拟南芥t-dna突变体表明，AtAnn1和AtAnn4在萌发过程中能应答渗透胁迫[16]。拟南芥AtAnn1基因是植物中研究较多的一个膜联蛋白，过表达AtAnn1导致耐旱性增强，而AtAnn1缺失突变株抗旱能力降低，对干旱更敏感[17]。印度芥菜（Brassica juncea）中，过表达AnnBj1分别增强了转基因烟草和棉花的耐旱性和耐盐性[18]。小麦 经低温胁迫可诱导膜联蛋白p39和p22.5的表达，嵌插在质膜上感受或转导钙信号的同时还调节胞质中的钙离子浓度适应胁迫环境[19]。虽然膜联蛋白在植物胁迫反应中发挥重要作用，但不同植物膜联蛋白的分子性能、生理和生物学功能是否有所不同，还属未知。
胡杨（Populus euphratica）适应非生物胁迫环境能力极强，其抗逆机制受到国内外广泛关注。胡杨耐旱性的生理及分子机制虽有较多的研究，但目前未有胡杨膜联蛋白与水分胁迫相关性的报道。本研究通过对比胡杨Annexin1转基因拟南芥、野生型、突变体（atann1）在渗透胁迫、干旱胁迫及复水后的表型及生理生化反应，明确PeAnn1在植物耐受水分胁迫中的作用，研究结果有助于揭示树木的抗旱性机制。
1. 材料与方法
1.1 实验材料
将来自新疆的胡杨1年生实生苗栽植在北京林业大学苗圃的温室，培养3个月后，利用250 mmol/L甘露醇处理12 h，在渗透胁迫处理0、4、8、12 h后采集叶片，进行RT-qPCR分析。
将野生型拟南芥（WT）、突变体（atann1）和2个过表达胡杨PeAnn1拟南芥（PeAnn1-OE1和PeAnn1-OE2）种子在超净台中用质量分数为1%的次氯酸钠消毒10 min后，无菌水冲洗5次，播种于1/2 MS固体培养基（含有10 g/L蔗糖和质量分数为0.3%的植物凝胶，pH 5.8）上，4 ℃低温暗培养春化2 d，取出置于22 ℃、16 h光照/8 h黑暗的人工气候培养箱中生长。生长7 d的幼苗移栽于草炭土和蛭石的混合比例为 2∶1(w∶w)的土壤中，转入22 ℃、16 h光照/8 h黑暗，55% ~ 65%相对湿度的温室中生长。
胡杨苗木培养，PeAnn1基因克隆，RT-qPCR分析，以及转基因植物的培育及分子检测参见文献[20]。RNA提取试剂盒购自康为世纪生物科技有限公司；cDNA反转录试剂盒购自普洛麦格生物技术有限公司；其余试剂由北京拜尔迪生物科技有限公司提供。引物合成由北京睿博兴科生物技术有限公司完成。RNA提取、纯化、RT-qPCR参见文献[21]。
1.2 PeAnn1序列与进化树分析
通过NCBI（http://blast.ncbi.nlm.nih.gov/Blast.cgi）在线获取不同植物Annexin同源基因序列。利用Mega 6软件（http://www.megasoftware.net/index.php）对胡杨PeAnn1、毛白杨（Populus tomentosa)PtAnn1、拟南芥AtAnn1、野生大豆（Glycine soja）GsANN和水稻OsAnn8等基因序列进行多重序列比对分析，并对不同物种同源Annexins构建系统进化树。
1.3 渗透胁迫下拟南芥耐旱性检测
1.3.1 种子萌发率
在1/2 MS固体培养基中添加不同浓度的甘露醇模拟渗透胁迫，甘露醇浓度梯度设为150、200、250、300 mmol/L 4个梯度，以不加甘露醇的培养基作为对照。拟南芥种子在超净台中用质量分数为1%的次氯酸钠消毒10 min后，无菌水冲洗5次，依次点播在各个培养基上，每个株系点播50粒种子，每种处理做3次重复。4 ℃低温暗培养2 d后置于22 ℃、16 h光照/8 h黑暗的恒温培养箱中培养，10 d后记录萌发率。
1.3.2 根 长
选取在1/2MS固体培养基中萌发5 d，长势相似的拟南芥幼苗，移栽到正常对照培养基和甘露醇浓度梯度150、200、250、300 mmol/L的固体培养基上，观察拍照并测量记录拟南芥竖直生长7 d后的根长。
1.4 渗透胁迫下转基因拟南芥H2O2含量的测定
采用荧光探针H2DCF-DA（molecular probe）染色测定H2O2含量[22]。将在1/2MS固体培养基上萌发、生长7 d的各株系拟南芥幼苗（WT、atann1、PeAnn1-OE1
和PeAnn1-OE2）移到含有250 mmol/L甘露醇的LMS培养基中，12 h后取出幼苗，将幼苗浸入10 μmol/L H2DCF-DA的探针溶液，室温避光孵育15 min。取出幼苗，用蒸馏水冲洗3 ~ 4次，去除表面残留探针染料。利用Leica SP5激光共聚焦显微镜（Wetzlar, 德国）测定H2DCF-DA绿色荧光，激发光波长488 nm，发射光波长510 ~ 530 nm，相对荧光强度由Image-Pro Plus version 4.5软件计算。 1.5 渗透胁迫下转基因拟南芥抗氧化酶活性的测定
采用紫外吸收法[23]测定过氧化氢酶（CAT）活性，以1 min内ΔA240减少0.01的酶量为1个过氧化氢酶活性单位（U）。采用NBT（氮蓝四唑）光化还原法[24]测定超氧化物歧化酶（SOD）活性，以抑制NBT光还原50%的酶量为1个酶活力单位（U）。采用愈创木酚法[25]测定过氧化物酶（POD）活性，以每分钟OD值变化（升高）0.01为1个酶活性单位（U）。按照考马斯亮蓝G250法进行蛋白含量的测定。
1.6 转基因拟南芥中抗逆相关基因表达水平的检测
将野生型拟南芥和转基因拟南芥播种于1/2MS培养基上进行培养，萌发5 d后长势相似的拟南芥移至250 mmol/L甘露醇的1/2MS培养基上进行处理。处理10 d后取样，利用TRIzon试剂（康为世纪）提取总RNA并反转录合成第1条cDNA，通过荧光定量PCR检测抗逆相关基因SOD（超氧化物歧化酶基因）、CAT（过氧化氢酶基因）、POD（过氧化物酶基因）的转录水平表达情况，以拟南芥AtACTIN2（At3g18780）为内参，各抗氧化酶基因荧光定量引物见表1。
表 1 本文实验中所用到的引物序列Table 1. Gene-specific primer sequences used in this study引物名称 Primer name 上游引物 Forward primer (5′−3′) 下游引物 Reverse primer (5′−3′) AtACTIN2 GGTAACATTGTGCTCAGTGGTGG AACGACCTTAATCTTCATGCTGC AtSOD AGGAAACATCACTGTTGGAGAT GAGTTTGGTCCAGTAAGAGGAA AtCAT AGGATCAAACTTTGAGGGGTAG CTTGTGGTTCCTGGAATCTACT AtPOD CGTGCCCTTCATATTGTTGG GACGCCATCAACAACGAGTC 1.7 自然干旱条件下拟南芥叶绿素荧光参数的测定
选取生长20 d，状况良好且一致的野生型、突变体和转基因株系通过控制浇水模拟土壤干旱，浇足水后当天作为干旱的第1天，设置5个重复。分别在正常浇水、干旱处理8 d和复水3 d后采样进行相关指标的测定。
1.7.1 叶绿素荧光参数
每次早上10：00开始测定，选取暗适应30 min后的拟南芥叶片。利用Junior PAM便携式荧光仪（WALZ，德国）测定叶片正常浇水，干旱处理8 d和复水3 d后的荧光参数，包括PSⅡ最大光量子效率（Fv/Fm）、实际光合量子产量（ΦPSⅡ）和相对电子传递速率（ETR）。各个参数直接从仪器导出，每个株系重复3 ~ 5次，取平均值。
1.7.2 叶绿素相对含量
利用日本Minolta公司生产的便携式叶绿素计（SPAD-502 Plus）原位测定叶绿素相对含量，以SPAD值表示[26]。测定时选取长势一致的叶片、叶缘和叶脉中间的部位。要注意避开叶脉集中的部位，同时适当遮挡直射阳光，以保证测量的准确性。每个株系重复3 ~ 5次，取平均值。
1.8 数据分析
实验数据采用Excel和SPSS 19.0软件进行数据统计、分析并绘图。用SPSS 19.0软件中的最小显著性差异法（LSD）进行单因素方差分析，用Duncan法进行多重比较分析，显著性水平设定为α = 0.05。
2. 结果与分析
2.1 渗透胁迫诱导胡杨PeAnn1基因表达
胡杨苗木经250 mmol/L甘露醇处理后，PeAnn1基因表达发生明显变化，PeAnn1基因相对表达量在渗透胁迫处理4 h后显著提高，之后下降，在渗透胁迫处理12 h后降至对照水平（图1）。
![]()
图 1 甘露醇处理对胡杨PeAnn1基因表达的影响胡杨1年生苗木经250 mmol/L甘露醇处理12 h，渗透胁迫处理0、4、8、12 h后采集叶片进行RT-qPCR分析，每个数值均为3次生物样本重复，*表示差异显著，P < 0.05。One-year-old seedlings of P. euphratica were subjected to 250 mmol/L mannitol and stressed for 12 hours. Then leaves were sampled after being treated for 0, 4, 8, 12 hours under osmotic stress and used for RT-qPCR analysis. Each column is the mean of three biological repeats. Asterisk represents significant differences at P < 0.05 level.Figure 1. Effects of mannitol stress on expression of PeAnn1 in Populus euphratica leaves2.2 胡杨PeAnn1序列与进化树分析
胡杨PeAnn1与毛白杨PtAnn1、拟南芥AtAnn1、大豆GsANN和水稻OsAnn8同源基因的相似度分别为85%、72.56%、56%、29.51%，胡杨PeAnn1与毛白杨PtAnn1的相似性最高（图2A）。系统进化树分析结果表明，PeAnn1与毛白杨PtAnn1的亲缘关系最为接近（图2B）。
![]()
图 2 胡杨PeAnn1基因序列比对与系统进化树分析A. 胡杨 PeAnn1基因与其他植物Annexin同源基因多重序列比对。B. 不同植物Annexin蛋白的系统进化树分析。胡杨PeAnn1（XM_011018027.1）、毛白杨PtAnn1（JX986594.1）、拟南芥AtAnn1（NM_103274.4）、野生大豆GsANN（GU474544.1）、水稻OsAnn8（LOC_Os09g20330.1）。A, multiple alignment of the deduced gene sequences of PeAnn1 with other Annexins from different plant species. B, phylogenetic tree analysis of Annexin proteins. P. euphratica Anneixn1（XM_011018027.1）, P. tomentosa PtAnn1（JX986594.1）, Arabidopsis thaliana AtAnn1（NM_103274.4）, Glycine soja GsANN（GU474544.1）, and Oryza sative OsAnn8（LOC_Os09g20330.1）.Figure 2. Multiple alignment of PeAnn1 gene sequences with other Annexins from different plant species and phylogenetic tree analysis of Annexin proteins2.3 渗透胁迫下拟南芥的萌发率
将各株系拟南芥种子播种于含有不同浓度甘露醇的1/2 MS培养基中培养，观察萌发率和生长状态。从图3A中可以看出，未经甘露醇处理时，各个基因型的拟南芥均正常生长，长势一致，无显著差异。但在含有甘露醇的培养基上，各株系的拟南芥的生长均受到抑制，且随着甘露醇浓度的升高，生长抑制作用逐渐加强。当甘露醇浓度达到250 mmol/L时，野生型、突变体种子萌发率急剧降低，降至75.8%、78.7%，转基因拟南芥OE1和OE2的种子受到的抑制作用更严重，萌发率降至69.0%、71.6%（图3B）。当甘露醇浓度达到300 mmol/L时，转基因株系OE1和OE2的种子萌发率分别为46.9%、48.9%，仅是对照组的47.8%、49.1%，与野生型和突变体间差异显著（P < 0.05, 图3B）。结果表明，过表达PeAnn1拟南芥株系OE1和OE2在甘露醇培养基上的长势弱于野生型和突变体，对渗透胁迫更敏感。
![]()
图 3 不同浓度甘露醇对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2 ）拟南芥萌发率的影响 A. 甘露醇对种子萌发生长的影响（拟南芥播种在含有不同浓度甘露醇的1/2MS培养基上，生长10 d后的照片）；B. 种子萌发率分析。不同字母表示在P < 0.05水平上差异显著，下同。A, effects of mannitol stress on seed germination rate (seeds of Arabidopsis thaliana are allowed to germinate on 1/2 MS medium supplemented with different concentrations of mannitol. Representative photographs are taken after 10 days of treatment); B, analysis on seed germination rate. Different letters denote significant differences at P < 0.05 level, the same below.Figure 3. Effects of different concentrations mannitol on seed germination rate of WT, atann1 and PeAnn1-transgenic Arabidopsis thaliana (PeAnn1-OE1 , PeAnn1-OE2) 2.4 渗透胁迫下拟南根长的变化
根长结果显示，在无甘露醇的对照培养基上，转基因拟南芥OE2根长生长低于其他株系（图4A、4B）。在含有不同浓度梯度（150、200、250、300 mmol/L）的甘露醇培养基上，各拟南芥的根长生长均受到抑制，且随着处理浓度的增加抑制作用逐渐加强，其长度逐渐缩短（图4A）。250 mmol/L浓度的甘露醇处理后，野生型和突变体株系根长分别降低了42.5%和46.3%，转基因株系OE1和OE2根的长度下降更显著，达到58.8%和57.7%（图4B）。当甘露醇的浓度达到300 mmol/L时，转基因OE1和OE2的根长分别降低了72.4%和70.7%（图4B），说明高浓度的甘露醇对拟南芥根生长产生极大的抑制作用。而渗透胁迫下野生型与突变体根长显著高于转基因拟南芥OE1和OE2，说明过表达PeAnn1影响了转基因株系的根长生长。根据甘露醇梯度实验结果，确定250 mmol/L甘露醇为拟南芥渗透胁迫处理的最适浓度，并进行H2O2含量、抗氧化酶活性和基因表达等指标的测定。
![]()
图 4 不同浓度甘露醇对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥根长生长的影响A. 甘露醇对根长生长的影响（拟南芥播种在1/2MS培养基萌发5 d后，移植到含有不同浓度甘露醇的1/2MS培养基垂直生长7 d后的照片）, 标尺为1 cm；B. 根系生长分析。A, effects of mannitol stress on root growth (seeds of Arabidopsis thaliana are allowed to germinate on 1/2 MS medium for 5 days, then transplanted to 1/2 MS medium supplemented with different concentrations of mannitol. Representative photographs are taken after 7 days of treatment). Scale bar = 1 cm; B, root growth analysis.Figure 4. Effects of different concentrations mannitol on root growth of WT, atann1 and PeAnn1-transgenicArabidopsis thaliana (PeAnn1-OE1, PeAnn1-OE2)2.5 渗透胁迫下拟南芥H2O2含量的变化
H2DCFDA是一种脂溶性荧光探针，具有细胞渗透性且自身无荧光的氧化反应指示剂。进入细胞后被细胞酯酶水解成非脂溶性的DCFH且不被质膜排出，之后在过氧化物酶的催化下可被H2O2氧化转变成为强荧光性的DCF，细胞内DCF的荧光强度与胞内H2O2水平呈显著正相关，因此被广泛用来检测细胞内活性氧（ROS）的含量和监测细胞氧化还原反应过程[27]。H2DCFDA的测定结果显示，在正常条件下，转基因株系OE1和OE2根细胞H2O2水平较高，与野生型、突变体拟南芥有显著差异。经250 mmol/L甘露醇渗透胁迫后，除转基因株系OE2外，各拟南芥株系的H2O2水平均有显著增加（图5），说明渗透胁迫下，植物体内积累了大量的活性氧，而且250 mmol/L甘露醇对转基因株系的作用更加明显：OE1和OE2 H2O2的水平分别是野生型拟南芥、突变体的1.19 ~ 1.36倍和1.37 ~ 1.55倍（图5），说明转基因拟南芥在渗透逆境下难以清除体内过多的H2O2，抗氧化防御系统能力弱于野生型拟南芥和突变体。
![]()
图 5 甘露醇胁迫下野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥根细胞H2O2荧光强度的变化 甘露醇（250 mmol/L）胁迫下拟南芥根细胞H2O2荧光强度的变化。H2O2 fluorescence intensity variation of Arabidopsis thaliana under mannitol stress (250 mmol/L).Figure 5. Effects of mannitol stress on H2O2 fluorescence intensity in root cells of wild-type (WT), atann1 and transgenic lines(PeAnn1-OE1, PeAnn1-OE2)2.6 渗透胁迫下拟南芥抗氧化物酶活性的变化
当植物处于不利于自身生长的条件时，活性氧增多，SOD、CAT、POD等抗氧化酶清除超氧阴离子和H2O2，从而保持活性氧的动态平衡，减轻ROS对植物的危害[28]。通过对250 mmol/L甘露醇胁迫处理后各株系中SOD活性的检测，发现甘露醇处理后各拟南芥体内的SOD活性均显著增加，野生型和突变体的SOD值变化幅度较大，分别较对照组上升了33.1%和52.8%，说明拟南芥幼苗的SOD对其在干旱胁迫下自由基的清除具有重要作用。转基因株系OE1和OE2的SOD活性增幅较低，且在渗透胁迫后，SOD活性显著低于野生型和突变体（图6A）。POD能清除对植物体有毒害作用的氧自由基，防止细胞衰老，提高植物的抗逆能力。由图6B可知，渗透胁迫后，各拟南芥的POD活性均呈不同程度的下降，其中野生型、转基因株系OE1和OE2与对照组间差异显著，分别下降了15.2%、33.4%、35.7%；突变体拟南芥无显著差异。CAT作为植物体内分解、清除H2O2的重要保护酶类之一。在250 mmol/L甘露醇胁迫条件下，水分的缺失诱导植物体内CAT活性提高，从而加速H2O2的分解，避免羟基、自由基的产生，对植物细胞膜的过氧化有一定减缓作用。图6C中结果显示，正常条件下各个拟南芥株系的CAT活性差异显著。胁迫处理后，野生型和突变体的CAT活性显著增加，而转基因株系OE1的CAT活性升高幅度较小，OE2的CAT活性显著降低，较处理前下降了29.3%。
![]()
图 6 甘露醇胁迫下野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥中SOD酶活性、POD酶活性和CAT酶活性的变化A. 甘露醇胁迫下SOD活性；B. 甘露醇胁迫下POD活性；C. 甘露醇胁迫下CAT活性。A, SOD activity under mannitol stress; B, POD activity under mannitol stress; C, CAT activity under mannitol stress.Figure 6. Effects of mannitol stress on the activities of SOD, POD and CAT in wild-type (WT), atann1 and transgenic lines (PeAnn1-OE1, PeAnn1-OE2)2.7 渗透胁迫下抗氧化酶基因表达量的变化
为了探究PeAnn1转基因株系在渗透胁迫下抗氧化酶活性下降的原因，本文检测了渗透胁迫后各基因型抗氧化酶基因的表达。图7A结果表明，渗透胁迫处理后，AtSOD基因的表达水平在野生型和突变体拟南芥中显著升高，与对照组相比分别上升了50%、38.5%；而在转基因株系OE1和OE2中显著下降。图7B中显示，渗透胁迫处理的野生型和突变体拟南芥AtPOD基因表达水平几乎没有变化，而转基因株系OE1和OE2与对照组间差异显著，AtPOD表达量分别下降了33.4%、35.7%。同样，在渗透胁迫后，各个拟南芥的AtCAT表达水平也均显著降低，转基因株系OE1和OE2的AtCAT基因表达水平下降幅度最大，分别较对照组下降了58.2%、66.1%（图7C）。综上可知，250 mmol/L甘露醇处理条件下转基因株系OE1和OE2体内AtSOD、AtPOD和AtCAT基因的表达水平显著均低于野生型和突变体（图7A~7C），说明过表达胡杨PeAnn1降低了拟南芥ROS清除酶基因的表达，导致不能清除体内过量的H2O2，抗旱性较弱。
![]()
图 7 甘露醇胁迫对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2 ）拟南芥抗氧化酶基因AtSOD、AtPOD和AtCAT表达的影响 A. 甘露醇胁迫下AtSOD基因表达；B. 甘露醇胁迫下AtPOD基因表达；C. 甘露醇胁迫下AtCAT基因表达；内参基因为AtACTIN2，n = 3。A, AtSOD expression under mannitol stress; B, AtPOD expression under mannitol stress; C, AtCAT expression under mannitol stress; reference gene: AtACTIN2, n = 3.Figure 7. Effects of mannitol stress on relative expression of antioxidant-enzyme genes (AtSOD, AtPOD and AtCAT) in wild-type (WT), atann1 and transgenic lines (PeAnn1-OE1, PeAnn1-OE2)![]()
图 8 干旱和复水后野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥的光合参数变化A. 干旱及复水后拟南芥的生长状况；B~D. 正常浇水、干旱8 d、复水3 d拟南芥的PSⅡ最大光量子效率（Fv/Fm）、相对电子传递效率（ETR）、实际光合量子产量（ΦPSⅡ）的变化。A, plant performance of Arabidopsis thaliana after drought stress and rewatering; B−D, maximum photon efficiency of PSⅡ（Fv/Fm）, relative electron transfer efficiency（ETR）, and actual photosynthetic quantum yield（ΦPSⅡ） of Arabidopsis thaliana under normal watering, 8 days of drought stress and 3 days of rewatering.Figure 8. Changes of photosynthetic parameters in wild-type (WT), atann1, and transgenic lines (PeAnn1-OE1, PeAnn1-OE2) under drought stress and rewatering2.8 干旱胁迫下拟南芥叶绿素荧光参数的变化
对各株系拟南芥进行土壤干旱处理，并测定叶绿素含量及荧光参数。干旱处理前，拟南芥生长状况没有明显差异，干旱胁迫8 d后植株失水枯黄，野生型和突变体拟南芥长势较好。复水3 d后，转基因株系死亡率较高（图8A）。
叶绿素荧光参数反映植物光反应中心活性和光合生理状况。胁迫条件下PSⅡ的最大光量子效率（Fv/Fm）的变化可反映PSⅡ的光抑制程度[29]。如图8B所示，在干旱胁迫前，野生型、突变体和转基因株系OE1和OE2的Fv/Fm为0.706 ~ 0.747，其中转基因株系OE2的Fv/Fm最低，与其他拟南芥株系有显著差异（P < 0.05）。干旱胁迫8 d后，各株系的Fv/Fm有不同程度的降低，转基因株系的Fv/Fm下降更为显著；复水3 d后，PSⅡ最大光化学效率均迅速恢复，但转基因株系的恢复幅度不及野生型和突变体（图8B）。说明PeAnn1转基因株系的光反应中心在干旱胁迫条件下受到逆境的伤害较大。
干旱处理前，转基因株系和突变体的相对电子传递效率（ETR）、实际光合量子产量（ΦPSⅡ）均低于野生型（图8C、8D）。干旱8 d后，转基因株系的ETR和ΦPSⅡ值均大幅度降低，明显高于野生型、突变体株系的下降幅度（图8C、8D）。复水3 d后，各个拟南芥株系的相对电子传递效率（ETR）和实际光合量子产量（ΦPSⅡ）均快速恢复，但转基因株系的恢复程度远低于野生型、突变体拟南芥（图8C、8D）。这说明，与光反应中心相似，PeAnn1转基因株系的捕光蛋白复合体在干旱胁迫下被严重破坏，电子传递受阻，光能转化效率降低。
2.9 干旱胁迫下拟南芥叶绿素含量的变化
叶绿素相对含量采用叶绿素仪（SPAD 仪）原位测定[22]。由图9可以看出，在干旱前，突变体atann1的SPAD值显著高于野生型和转基因株系OE1和OE2（P < 0.05）。干旱8 d后，SPAD值均显著降低，其中转基因株系下降幅度最大。复水后第3天，各株系拟南芥的SPAD值均有所增加，但仍远低于干旱前的水平，特别是转基因株系。这表明干旱胁迫使转基因拟南芥叶绿素含量显著下降。与转基因株系相比，突变体atann1能在干旱和复水期间维持较高的叶绿素含量。
![]()
图 9 干旱和复水后野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥的叶绿素相对含量拟南芥正常浇水、干旱8 d、复水3 d叶绿素SPAD值的变化。Chlorophyll SPAD values of Arabidopsis thaliana under normal watering, 8 days of drought stress and 3 days of rewatering.Figure 9. Changes of relative chlorophyll content in wild-type (WT), atann1 and transgenic plant (PeAnn1-OE1, PeAnn1-OE2) under drought stress and rewatering3. 讨 论
渗透胁迫处理诱导了胡杨叶片PeAnn1基因瞬时的上调表达（图1），表明PeAnn1参与了胡杨对水分胁迫的响应。为了明确胡杨PeAnn1在植物耐受水分胁迫中的作用，本文对野生型（WT）、Annexin1突变体（atann1）和PeAnn1转基因株系OE1和OE2进行渗透胁迫和干旱胁迫处理，研究各株系的表型差异。结果表明，甘露醇处理后，WT、atann1和转基因株系OE1和OE2的种子萌发率和根长均受到抑制，且抑制作用随甘露醇浓度提高而增强（图3、4）。在高渗透胁迫（250 ~ 300 mmol/L甘露醇）条件下，转基因植株的萌发率、根长显著低于野生型和突变体。干旱胁迫下，转基因株系的光合参数下降幅度高于野生型和突变体（图8、9），这些表型结果说明，过表达胡杨PeAnn1使转基因植株对水分胁迫的敏感性提高。
PeAnn1转基因株系对水分胁迫耐受性的下降与光合参数的下降有关。叶绿素含量、PSⅡ最大光量子效率（Fv/Fm）、相对电子传递效率（ETR）和实际光合量子产量（ΦPSⅡ）的高低反映植物的光合能力和生长状况[30]。SPAD 值是一项表征叶片光合活性的指标，它与单位叶面积的叶绿素含量及叶绿素的密度显著正相关[31]，因此，常用SPAD 值的大小来衡量叶片中叶绿素含量。干旱胁迫下SPAD显著降低（图9），说明在干旱条件下叶绿素蛋白质的合成受到抑制，叶绿素逐渐分解，光合色素含量下降。干旱还会破坏捕光蛋白复合体，使电子传递受阻，光能转化效率降低，导致Fv/Fm、ETR、ΦPSⅡ显著下降（图8）。复水后，光合参数有所回升，说明光系统并未完全受到损伤。但与野生型和突变体相比，转基因株系修复损伤能力最弱，这与其干旱期间光系统受损的程度有关（图8、9）。
在逆境胁迫下,植株体内会生成大量的超氧自由基和过氧化氢，SOD作为植物抗氧化系统的第一道防线，通过歧化超氧阴离子生成H2O2，生成的H2O2进一步被CAT、POD等分解消除，保护细胞免受氧化伤害[32]。渗透胁迫下各株系H2O2水平提高（图5），野生型和突变体维持较高抗氧化酶活性，有助于活性氧的清除（图6）。但转基因株系中SOD、CAT、POD酶活较低（图6），不能有效清除活性氧（图5），减轻氧化胁迫对细胞膜的伤害，难以保持光合系统的相对稳定。这与在印度芥菜[18]和小麦[19]中的研究结果不完全一致，表明不同物种膜联蛋白在抗氧化防御中的作用有所差异。众多研究表明，抗氧化酶活性与编码基因表达的表达量相关[33-34]。在渗透胁迫后PeAnn1转基因株系OE1和OE2抗氧化酶SOD、POD和CAT的活性和基因表达量均显著低于野生型和突变体（图6、7），表明过表达PeAnn1株系在干旱胁迫下不能及时清除体内大量的活性氧，破坏了细胞内蛋白质的结构和功能。
关于膜联蛋白在植物耐受水分逆境中的作用有不同的研究报道。Lee等人[16]用甘露醇模拟渗透胁迫，发现过表达AtAnn1表现强耐旱性。这与Konopka-Postupolska等人[17]报道的结果类似，但不同于Huh等人[35]研究的结果。PeAnn1基因序列与拟南芥AtAnn1有较大差别（图2），因此，PeAnn1在调控植物抗旱性中作用可能有所不同。王希等人[36]通过超表达GsANN发现转基因拟南芥提高了对盐胁迫和干旱胁迫的敏感性，却志群等人[37]研究水稻OsAnn8基因可能通过负调控方式参与水稻对干旱胁迫的响应。这些研究与本文结论类似。另外，Liao等[38]发现拟南芥MYB30转录因子调控植物胁迫响应的信号转导机制。拟南芥MYB30通过结合ANN1和ANN4基因的启动子，抑制ANN1和ANN4基因表达，从而影响ANNs介导的特异胞内钙信号的产生，以此调控植物的氧化胁迫和热胁迫响应过程。胡杨PeAnn1是否也与相关转录因子互作，通过调控胞内钙信号影响抗氧化酶基因的表达，有待深入研究。
综上所述，本文从生理和分子转录水平上揭示了胡杨PeAnn1基因对植物抗旱性的影响，研究结果有助于阐明胡杨的抗旱机制。
-
![]()
图 2 各地类提供生态系统服务价值（ESV）及各地类ESV变化生态贡献率
Figure 2. Ecosystem service value (ESV) provided by different regions and ecological contribution rates of ESV changes in different regions
![]()
图 3 2001—2020年黄河流域土地利用及ESV空间分布
Figure 3. Spatial distribution of land use and ESV in the Yellow River Basin from 2001 to 2020
![]()
图 4 2001—2020年黄河流域生态系统服务价值LISA图
Figure 4. LISA diagram of ESV in the Yellow River Basin from 2001 to 2020
![]()
图 5 2001—2020年LISA时间路径相对长度、弯曲度与移动方向
Figure 5. Relative length, curvature and moving direction of LISA time path from 2001 to 2020
![]()
图 6 黄河流域ESV驱动因素空间异质性
括号内为显著性比例，显著性水平比例为 GTWR 估计的参数值在显著性水平为 0.05 时通过检验的比例。The proportion of significance is shown in parentheses, and the proportion of significance level is the proportion of parameter values estimated by GTWR that pass the test at a significance level of 0.05.
Figure 6. Spatial heterogeneity of ESV driving factors in the Yellow River Basin
表 1 黄河流域生态系统服务价值系数 元/hm2
Table 1 Ecosystem service value (ESV) coefficients of the Yellow River Basin
CNY/ha 服务类别
Service type单项服务
Individual service林地
Forest land灌木林地
Shrub land草地
Grassland湿地
Wetland耕地
Cropland水域
Waters供给服务
Provision service食物生产
Food production420.39 295.83 155.70 794.06 1 728.26 622.80 原料生产
Material production980.90 669.51 217.98 778.50 389.25 186.84 水资源供给
Water resource supply513.81 342.54 124.56 4 032.60 −2 039.66 8 143.06 调节服务
Regulating service气体调节
Gas regulation3 222.97 2 195.36 794.06 2 958.28 1 385.72 747.36 气候调节
Climate regulation9 653.34 6 586.07 2 086.37 5 605.16 731.79 2 210.93 净化环境
Environmental purification2 802.58 1 992.95 685.08 5 605.16 217.98 4 452.99 水文调节
Hydrological regulation6 009.98 5 215.92 1 525.85 37 725.87 2 335.49 85 151.78 支持服务
Supporting service土壤保持
Soil conservation3 923.61 2 678.02 965.33 3 596.65 809.63 731.79 维持养分循环
Maintaining nutrient cycle295.83 202.41 77.85 280.26 249.12 62.28 生物多样性
Biodiversity3 581.08 2 444.47 871.91 12 253.51 264.69 1 992.95 文化服务
Cultural service美学景观
Aesthetic scenery1 572.56 1 074.32 389.25 7 364.56 124.56 1 541.42 总计 Total 32 977.05 23 697.40 7 893.94 80 994.61 6 196.83 105 844.20 
下载: 导出CSV
表 2 2001—2020年黄河流域ESV变化 108元
Table 2 Changes in ESV of the Yellow River Basin from 2001 to 2020
108 CNY 服务类别
Service
type单项服务
Individual
service2001 2005 2010 2015 2020 ESV 比例
Proportion/%ESV 比例
Proportion/%ESV 比例
Proportion/%ESV 比例
Proportion/%ESV 比例
Proportion/%供给服务
Provision
service食物生产
Food
production350.13 5.16 365.16 5.31 373.72 5.35 379.38 5.40 383.28 5.37 原料生产
Material
production214.50 3.16 218.57 3.18 222.55 3.19 224.26 3.19 227.33 3.19 水资源供给
Water resource
supply−181.01 −2.67 −197.89 −2.88 −204.93 −2.93 −210.78 −3.00 −212.11 −2.97 调节服务
Regulating
service气体调节
Gas regulation765.74 11.28 779.27 11.34 792.53 11.34 798.11 11.35 808.05 11.33 气候调节
Climate
regulation1 638.88 24.14 1 651.97 24.04 1 677.24 24.01 1 683.77 23.95 1 708.39 23.95 净化环境
Environmental
purification532.07 7.84 535.67 7.80 543.45 7.78 545.73 7.76 553.41 7.76 水文调节
Hydrological
regulation1 600.55 23.57 1 633.02 23.76 1 667.12 23.86 1 684.68 23.96 1 713.89 24.02 支持服务
Supporting
service土壤保持
Soil
conservation807.13 11.89 815.81 11.87 827.98 11.85 831.87 11.83 842.40 11.81 维持养分循环
Maintaining
nutrient cycle90.40 1.33 92.67 1.35 94.42 1.35 95.31 1.36 96.41 1.35 生物多样性
Biodiversity670.01 9.87 674.48 9.81 684.66 9.80 688.28 9.79 698.57 9.79 文化服务
Cultural
service美学景观
Aesthetic
scenery301.13 4.44 303.28 4.41 308.04 4.41 309.99 4.41 314.84 4.41 总计 Total 6 789.52 100.00 6 872.00 100.00 6 986.78 100.00 7 030.61 100.00 7 134.48 100.00 变化率 Rate of change/% 1.21 1.67 0.63 1.48
下载: 导出CSV
表 3 黄河流域ESV全局Moran’s I
Table 3 Global Moran’s I of ESV in the Yellow River Basin
项目 Item 2001 2005 2010 2015 2020 Moran’s I 0.640 0.644 0.651 0.645 0.650 Z 81.934 82.474 83.368 82.506 83.121 P < 0.001 < 0.001 < 0.001 < 0.001 < 0.001
下载: 导出CSV
表 4 黄河流域ESV Moran’s I转移概率矩阵与时空跃迁
Table 4 Moran’s I transition probability matrix and spatiotemporal transition of ESV in the Yellow River Basin
转移矩阵
Transfer matrix2001—2020 类型
Type比例
Proportion变迁
Transition凝聚
CondensationHH LH LL HL HH 0.103 0.009 0.016 0.002 Type0 0.779 0.197 0.801 LH 0.013 0.046 0.035 0.000 Type1 0.142 LL 0.005 0.013 0.618 0.006 Type2 0.055 HL 0.006 0.002 0.114 0.012 Type3 0.024
下载: 导出CSV
-
[1] 石涛. 黄河流域生态保护与经济高质量发展耦合协调度及空间网络效应[J]. 区域经济评论, 2020(3): 25−34. Shi T. Spatial correlation network and regional connected effect of coupling coordination degree between ecological protection and high-quality economic development in the Yellow River Regions[J]. Regional Economic Review, 2020(3): 25−34.
[2] 徐勇, 王传胜. 黄河流域生态保护和高质量发展: 框架、路径与对策[J]. 中国科学院院刊, 2020, 35(7): 875−883. Xu Y, Wang C S. Ecological protection and high-quality development in the Yellow River Basin: framework, path, and countermeasure[J]. Bulletin of Chinese Academy of Sciences, 2020, 35(7): 875−883.
[3] 李婷, 吕一河, 任艳姣, 等. 黄土高原植被恢复成效及影响因素[J]. 生态学报, 2020, 40(23): 8593−8605. Li T, Lü Y H, Ren Y J, et al. Gauging the effectiveness of vegetation restoration and the influence factors in the Loess Plateau[J]. Acta Ecologica Sinica, 2020, 40(23): 8593−8605.
[4] 王艳芬, 陈怡平, 王厚杰, 等. 黄河流域生态系统变化及其生态水文效应[J]. 中国科学基金, 2021, 35(4): 520−528. Wang Y F, Chen Y P, Wang H J, et al. Ecosystem change and its ecohydrological effect in the Yellow River Basin[J]. Bulletin of National Natural Science Foundation of China, 2021, 35(4): 520−528.
[5] Groot R, Brander L, van der Ploeg S, et al. Global estimates of the value of ecosystems and their services in monetary units[J]. Ecosystem Services, 2012, 1(1): 50−61. doi: 10.1016/j.ecoser.2012.07.005
[6] Ouyang Z, Zheng H, Xiao Y, et al. Improvements in ecosystem services from investments in natural capital[J]. Science, 2016, 352: 1455−1459. doi: 10.1126/science.aaf2295
[7] 陈玉福, 董鸣. 生态学系统的空间异质性[J]. 生态学报, 2003, 23(2): 346−352. Chen Y F, Dong M. Spatial heterogeneity in ecological systems[J]. Acta Ecologica Sinica, 2003, 23(2): 346−352.
[8] Rey S J, Janikas M V. STARS: space-time analysis of regional systems[J]. Geographical Analysis, 2006, 38(1): 67−86. doi: 10.1111/j.0016-7363.2005.00675.x
[9] 景海超, 刘颖慧, 贺佩, 等. 青藏高原典型区生态系统服务空间异质性及其影响因素: 以那曲市为例[J]. 生态学报, 2022, 42(7): 2657−2673. Jing H C, Liu Y H, He P, et al. Spatial heterogeneity of ecosystem services and it’s influencing factors in typical areas of the Qinghai-Tibet Plateau: a case study of Nagqu City[J]. Acta Ecologica Sinica, 2022, 42(7): 2657−2673.
[10] 王耕, 张芙榕. 近30年辽河三角洲生态系统服务价值时空演变及影响因素分析[J]. 环境科学, 2024, 45(1): 228−238. Wang G, Zhang F R. Spatial and temporal evolution and impact factors analysis of ecosystem service value in the Liaohe River Delta over the past 30 years[J]. Environmental Science, 2024, 45(1): 228−238.
[11] 黄木易, 方斌, 岳文泽, 等. 近20 a来巢湖流域生态服务价值空间分异机制的地理探测[J]. 地理研究, 2019, 38(11): 2790−2803. Huang M Y, Fang B, Yue W Z, et al. Spatial differentiation of ecosystem service values and its geographical detection in Chaohu Basin during 1995−2017[J]. Geographical Research, 2019, 38(11): 2790−2803.
[12] 夏楚瑜, 国淏, 赵晶, 等. 京津冀地区生态系统服务对城镇化的多空间尺度动态响应[J]. 生态学报, 2023, 43(7): 2756−2769. Xia C Y, Guo H, Zhao J, et al. Dynamic responses of ecosystem services to urbanization at multi-spatial scales in the Beijing-Tianjin-Hebei region[J]. Acta Ecologica Sinica, 2023, 43(7): 2756−2769.
[13] 赵泳春, 苏方林. 经济差异化增长下生态系统服务价值的时空演变特征: 以珠江–西江经济带为例[J]. 自然资源学报, 2022, 37(7): 1782−1798. doi: 10.31497/zrzyxb.20220709 Zhao Y C, Su F L. Spatio-temporal dynamic characteristics of the ecosystem service values under differential economic growth: a case study of the Pearl River-West River Economic Belt[J]. Journal of Natural Resources, 2022, 37(7): 1782−1798. doi: 10.31497/zrzyxb.20220709
[14] Huang B, Wu B, Barry M. Geographically and temporally weighted regression for modeling spatio-temporal variation in house prices[J]. International Journal of Geographical Information Science, 2010, 24(3): 383−401. doi: 10.1080/13658810802672469
[15] Fotheringham A S, Charlton M E, Brunsdon C. Geographically weighted regression: a natural evolution of the expansion method for spatial data analysis[J]. Environment and Planning A, 1998, 30(11): 1905−1927. doi: 10.1068/a301905
[16] 吴小影, 杨山, 尹上岗, 等. 基于GTWR模型的长三角地区城市建设用地时空动态特征及其驱动机理[J]. 长江流域资源与环境, 2021, 30(11): 2594−2606. Wu X Y, Yang S, Yin S G, et al. Spatial-temporal dynamic characteristics and its driving mechanism of urban built-up area in Yangtze River Delta based on GTWR model[J]. Resources and Environment in the Yangtze Basin, 2021, 30(11): 2594−2606.
[17] 王海军, 张彬, 刘耀林, 等. 基于重心-GTWR模型的京津冀城市群城镇扩展格局与驱动力多维解析[J]. 地理学报, 2018, 73(6): 1076−1092. Wang H J, Zhang B, Liu Y L, et al. Multi-dimensional analysis of urban expansion patterns and their driving forces based on the center of gravity-GTWR model: a case study of the Beijing-Tianjin-Hebei urban agglomeration[J]. Acta Geographica Sinica, 2018, 73(6): 1076−1092.
[18] 田晓华, 卢飞, 徐翠蓉. 基于超效率SBM模型的黄河中下游生态效率测度与空间相关性分析[J]. 信阳师范学院学报(自然科学版), 2023, 36(1): 89−97. Tian X H, Lu F, Xu C R. Ecological efficiency measurement and spatial correlation analysis in the middle and lower reaches of the Yellow River based on the super-SBM model[J]. Journal of Xinyang Normal University (Natural Science Edition), 2023, 36(1): 89−97.
[19] 谢高地, 张彩霞, 张雷明, 等. 基于单位面积价值当量因子的生态系统服务价值化方法改进[J]. 自然资源学报, 2015, 30(8): 1243−1254. Xie G D, Zhang C X, Zhang L M, et al. Improvement of the evaluation method for ecosystem service value based on per unit area[J]. Journal of Natural Resources, 2015, 30(8): 1243−1254.
[20] 郑晓豪, 陈颖彪, 郑子豪, 等. 湖北省生态系统服务价值动态变化及其影响因素演变[J]. 生态环境学报, 2023, 32(1): 195−206. Zheng X H, Chen Y B, Zheng Z H, et al. Dynamic changes of ecosystem service value and evolution of its influencing factors in Hubei Province[J]. Ecology and Environmental Sciences, 2023, 32(1): 195−206.
[21] 封建民, 郭玲霞, 李晓华. 汉中市土地利用时空变化及其对生态系统服务价值的影响[J]. 水土保持研究, 2020, 27(1): 275−282. Feng J M, Guo L X, Li X H. Temporal and spatial variations of land uses and their influences on ecosystem service function values in Hanzhong City[J]. Research of Soil and Water Conservation, 2020, 27(1): 275−282.
[22] 乔斌, 祝存兄, 曹晓云, 等. 格网尺度下青海玛多县土地利用及生态系统服务价值空间自相关分析[J]. 应用生态学报, 2020, 31(5): 1660−1672. Qiao B, Zhu C X, Cao X Y, et al. Spatial autocorrelation analysis of land use and ecosystem service value in Maduo County, Qinghai Province, China at the grid scale[J]. Chinese Journal of Applied Ecology, 2020, 31(5): 1660−1672.
[23] 王蓓, 赵军, 胡秀芳. 基于 InVEST 模型的黑河流域生态系统服务空间格局分析[J]. 生态学杂志, 2016, 35(10): 2783−2792. Wang B, Zhao J, Hu X F. Spatial pattern analysis of ecosystem services based on InVEST in Heihe River Basin[J]. Chinese Journal of Ecology, 2016, 35(10): 2783−2792.
[24] Ye X, Rey S. A framework for exploratory space-time analysis of economic data[J]. The Annals of Regional Science, 2013, 50(1): 315−339. doi: 10.1007/s00168-011-0470-4
[25] 张中浩, 聂甜甜, 高阳, 等. 长三角城市群生态安全评价与时空跃迁特征分析[J]. 地理科学, 2022, 42(11): 1923−1931. Zhang Z H, Nie T T, Gao Y, et al. Ecological security assessment and spatio-temporal transition characteristics in the Yangtze River Delta urban agglomeration[J]. Scientia Geographica Sinica, 2022, 42(11): 1923−1931.
[26] Rey S J, Murray A T, Anselin L. Visualizing regional income distribution dynamics[J]. Letters in Spatial and Resource Sciences, 2011, 4(1): 81−90. doi: 10.1007/s12076-010-0048-2
[27] Rey S J. Spatial empirics for economic growth and convergence[J]. Geographical Analysis, 2001, 33(3): 195−214. doi: 10.1111/j.1538-4632.2001.tb00444.x
[28] Liu Z F, He C Y, Zhang Q F, et al. Extracting the dynamics of urban expansion in China using DMSP-OLS nighttime light data from 1992 to 2008[J]. Landscape and Urban Planning, 2012, 106(1): 62−72. doi: 10.1016/j.landurbplan.2012.02.013
[29] 罗盛锋, 闫文德. 广西北部湾沿岸地区生态系统服务价值变化及其驱动力[J]. 生态学报, 2018, 38(9): 3248−3259. Luo S F, Yan W D. Evolution and driving force analysis of ecosystem service values in Guangxi Beibu Gulf coastal areas, China[J]. Acta Ecologica Sinica, 2018, 38(9): 3248−3259.
[30] 李英利. 财政压力、土地财政与区域房价的时空演化: 基于GTWR模型的实证研究[J]. 财政研究, 2020(5): 78−89. Li Y L. Fiscal pressure, land finance and the spatio-temporal evolution of regional house price: empirical research based on GTWR model[J]. Public Finance Research, 2020(5): 78−89.
[31] 张琨, 林乃峰, 徐德琳, 等. 中国生态安全研究进展: 评估模型与管理措施[J]. 生态与农村环境学报, 2018, 34(12): 1057−1063. Zhang K, Lin N F, Xu D L, et al. Research advance on ecological security in China: assessment models and management measures[J]. Journal of Ecology and Rural Environment, 2018, 34(12): 1057−1063.
[32] 王聪, 伍星, 傅伯杰, 等. 重点脆弱生态区生态恢复模式现状与发展方向[J]. 生态学报, 2019, 39(20): 7333−7343. Wang C, Wu X, Fu B J, et al. Ecological restoration in the key ecologically vulnerable regions: current situation and development direction[J]. Acta Ecologica Sinica, 2019, 39(20): 7333−7343.
[33] 李宗善, 杨磊, 王国梁, 等. 黄土高原水土流失治理现状、问题及对策[J]. 生态学报, 2019, 39(20): 7398−7409. Li Z S, Yang L, Wang G L, et al. The management of soil and water conservation in the Loess Plateau of China: present situations, problems, and counter-solutions[J]. Acta Ecologica Sinica, 2019, 39(20): 7398−7409.
[34] 孙梦华, 牛文浩, 张蚌蚌, 等. 黄河流域土地利用变化下生态系统服务价值时空演变及其响应: 以陕甘宁地区为例[J]. 应用生态学报, 2021, 32(11): 3913−3922. Sun M H, Niu W H, Zhang B B, et al. Spatial-temporal evolution and responses of ecosystem service value under land use change in the Yellow River Basin: a case study of Shaanxi-Gansu-Ningxia region, Northwest China[J]. Chinese Journal of Applied Ecology, 2021, 32(11): 3913−3922.
[35] Austin A T, Vitousek P M. Nutrient dynamics on a precipitation gradient in Hawaii[J]. Oecologia, 1998, 113(4): 519−529. doi: 10.1007/s004420050405
[36] 张明阳, 王克林, 刘会玉, 等. 桂西北典型喀斯特区生态系统服务价值对景观格局变化的响应[J]. 应用生态学报, 2010, 21(5): 1174−1179. Zhang M Y, Wang K L, Liu H Y, et al. Responses of ecosystem service values to landscape pattern change in typical karst area of northwest Guangxi, China[J]. Chinese Journal of Applied Ecology, 2010, 21(5): 1174−1179.
[37] 王婧, 康荣华, 方媛, 等. 宁夏清水河平原生态脆弱性评价及驱动因子分析[J]. 水利水电快报, 2023, 44(8): 98−105. Wang J, Kang R H, Fang Y, et al. Eco-environmental vulnerability evaluation in Qingshui River Plain of Ningxia and driving force analysis[J]. Express Water Resources & Hydropower Information, 2023, 44(8): 98−105.
[38] 王殿茹, 薛少冰. 中国人口密度对生态环境的影响: 基于双向固定效应模型的实证分析[J]. 河北地质大学学报, 2022, 45(5): 97−103. Wang D R, Xue S B. The impact of population density on ecosystem in China based on the two-way fixed effects model[J]. Journal of Hebei Geo University, 2022, 45(5): 97−103.
[39] 桑海云. 财政转移支付、就业机会与相对贫困缓解[J]. 技术经济与管理研究, 2022(9): 124−128. Sang H Y. Fiscal transfer payments, employment opportunities and relative poverty alleviation[J]. Journal of Technical Economics & Management, 2022(9): 124−128.
[40] 安彬, 肖薇薇, 朱妮, 等. 近60 a黄土高原地区降水集中度与集中期时空变化特征[J]. 干旱区研究, 2022, 39(5): 1333−1344. An B, Xiao W W, Zhu N, et al. Temporal and spatial variations of precipitation concentration degree and precipitation concentration period on the Loess Plateau from 1960 to 2019[J]. Arid Zone Research, 2022, 39(5): 1333−1344.
[41] 程永生, 张德元, 赵梦婵. 黄河流域生态保护和高质量发展的时空演变与驱动因素[J]. 经济体制改革, 2021(5): 61−69. Cheng Y S, Zhang D Y, Zhao M C. Spatial-temporal evolution and driving factors of ecological protection and high-quality development in the Yellow River Basin[J]. Reform of Economic System, 2021(5): 61−69.
[42] 杨洁, 谢保鹏, 张德罡. 基于InVEST和CA-Markov模型的黄河流域碳储量时空变化研究[J]. 中国生态农业学报(中英文), 2021, 29(6): 1018−1029. Yang J, Xie B P, Zhang D G. Spatio-temporal evolution of carbon stocks in the Yellow River Basin based on InVEST and CA-Markov models[J]. Chinese Journal of Eco-Agriculture, 2021, 29(6): 1018−1029.
[43] 张世举, 王娟. 黄河河南段滩涂水资源利用研究[J]. 节水灌溉, 2013(5): 37−39. Zhang S J, Wang J. Comprehensive utilization of water resources in the Yellow River intertidal zone of Henan Province[J]. Water Saving Irrigation, 2013(5): 37−39.
[44] Li F Z, Yin X X, Shao M. Natural and anthropogenic factors on China’s ecosystem services: comparison and spillover effect perspective[J/OL]. Journal of Environmental Management, 2022, 324: 116064[2023−05−30]. https://doi.org/10.1016/j.jenvman.2022.116064.
[45] 危小建, 辛思怡, 张颖艺, 等. 不同格网尺度下生态服务价值空间分异及其影响因素差异探析[J]. 生态学报, 2023, 43(18): 7585−7597. Wei X J, Xin S Y, Zhang Y Y, et al. Spatial difference of ecological services and its influencing factors under different scales: taking the Nanchang Urban Agglomeration as an example[J]. Acta Ecologica Sinica, 2023, 43(18): 7585−7597.
-
期刊类型引用(3)
1. 张欢欢，范挺秀，高双成，范丙友，史国安. 2个牡丹切花品种核型、花药形态与花粉特性的比较. 浙江农业学报. 2025(02): 329-337 . 
百度学术
2. 孙榕泽，宋焕芝，蒋珈琦，于晓南. 三倍体芍药品种染色体制片优化和核型分析. 河北农业大学学报. 2024(03): 46-55 . 百度学术
3. 段英姿，客绍英，王晓英，张胜珍，马艳芝. 15个油用紫斑牡丹品种核型及亲缘关系分析. 种子. 2023(07): 96-104 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 359
- HTML全文浏览量: 85
- PDF下载量: 63
- 被引次数: 3
