基于EALCO模型对中亚热带人工针叶林CO2通量季节变异的模拟

 EALCO模型是一个基于生理生态学过程,模拟生态系统下垫面与大气之间水、热和碳通量交换的综合模型。将该模型应用在亚热带常绿针叶林, 对其生态系统过程进行了模拟,以深入探讨季节性干旱对生态系统过程的影响。对EALCO模型进行了参数化与初始化并对模型的光合作用时段和 落叶机制进行了改进,以更好地模拟亚热带人工针叶林生态系统。千烟洲通量观测站自2002年底开始应用涡度相关技术对中亚热带人工针叶林 生态系统进行通量观测,该站点2003年经历了一次较严重的季节性干旱(由高温与少雨综合作用造成),降水量仅为多年平均值的65%,而2004年 的年降水量与多年平均值较为接近,利用该站点2003和2004年特殊的气候条件,使用其通量观测数据对模型的模拟效果进行检验。从模拟结果 的总体趋势来看,模型能较好地从半小时、日及年尺度上反映两年内土壤-植被-大气之间的碳交换状况。总初级生产力(Gross primary production, GPP)在一年中呈现单峰型变化,遇高温及干旱胁迫GPP值下降。由于受到干旱胁迫的影响,2003年GPP值比2004年偏低12.9%。模拟 结果显示,2003年GPP值比2004年偏低11.2%。观测数据与模拟结果均显示,水分胁迫期间净碳交换量(Net ecosystem production, NEP)模拟值 与实测值的日变化均呈现一种“偏态",即一天中生态系统碳交换量最大值出现在上午某一时刻,之后逐渐降低。 模拟结果显示,水分匮缺对 光合能力的影响比对生态系统呼吸作用的影响更为强烈,因而导致了净生态系统生产力的降低。进一步分析表明,水分匮缺期间,晴天正午之 前,深层土壤( >20 cm) 水分的匮缺抑制了光合作用能力,正午之后,高温与深层土壤水分匮缺共同削弱光合作用能力,影响各占一半。     

基于FORCCHN模型的中国典型森林生态系统碳通量模拟

森林生态系统是陆地碳循环的重要组成部分,其固碳能力显著高于其他陆地生态系统,研究森林生态系统碳通量是认识和理解全球变化对碳循环影响的关键。碳循环模型是研究森林生态系统碳通量有效工具。以长白山温带落叶阔叶林、千烟洲亚热带常绿针叶林、鼎湖山亚热带常绿阔叶林和西双版纳热带雨林等4种中国典型森林生态系统为研究对象,利用涡度相关2003-2012年观测数据,评估FORCCHN模型对生态系统呼吸（ER）,总初级生产力（GPP）,净生态系统生产力（NEP）的模型效果。结果表明：（1） FORCCHN模型能够较好的模拟中国4种典型森林生态系统不同时间尺度的碳通量。落叶阔叶林和常绿针叶林ER和GPP的逐日变化模拟效果较好（ER的相关系数分别为0.94和0.92,GPP的相关系数分别为0.86和0.74）;（2）4种森林生态系统碳通量季节动态模拟值和观测值显著相关（P<0.01）,ER、GPP、NEP的观测值和模拟值的R²分别为0.77-0.93、0.54-0.88和0.15-0.38;模型可以很好地模拟森林生态系统不同季节碳汇（NEP>0）,碳源（NEP<0）的变化规律;（3）4种森林生态系统碳通量模拟值与观测值的年际变化有很好的吻合度,但在数值大小上存在差异,模型高估了常绿阔叶林的ER和GPP,略微低估了其他3种森林生态系统ER和GPP。     

基于模型数据融合的千烟洲亚热带人工林碳水通量模拟

人工林生态系统是我国森林生态系统的重要组成部分,在全球碳平衡中的作用越来越受到重视.利用千烟洲亚热带人工针叶林通量观测站的碳水通量和气象观测数据,通过模型数据融合方法对碳水循环过程模型——SIPNET模型关键参数进行反演,模拟了2004-2009年千烟洲人工林生态系统的碳水通量.结果表明:仅用碳通量观测数据优化模型参数时,净生态系统碳交换量(NEE)模拟效果较好(R2=0.934),而生态系统蒸散(ET)模拟效果较差(R2=0.188);同时用碳水通量观测数据优化时,NEE模拟效果稍差(R2=0.929),但ET模拟效果显著提升(R2=0.824),说明利用碳水通量观测数据同时优化,SIPNET模型才能较好地模拟试验站点碳水通量.在此基础上,开展了人工林生态系统碳通量对降水变化响应的敏感性分析,发现降水量减少对光合作用的影响比对呼吸作用的影响更为强烈,且碳水通量同时参与优化时模型才能较好地模拟碳通量随降水减少而快速降低的趋势,表明如果不能同时利用碳水通量进行参数优化,模型无法正确揭示生态系统碳循环对降水变异的响应.     

DROUGHT EFFECTS ON CARBON EXCHANGE IN A SUBTROPICAL CONIFEROUS PLANTATION IN CHINA

 干旱对陆地生态系统的影响已成为全球变化研究的焦点问题之一。该研究基于生态系统过程模型——CEVSA2, 结合涡度相关通量观测, 分析了不同程度干旱对亚热带人工针叶林碳交换的影响及其关键控制因素。结果表明: 1)干旱使生态系统碳交换显著下降, 2003和2004年的干旱使得年净生态系统生产力(Net ecosystem production, NEP)相比无干旱影响情景的模拟结果分别减少了63%和47%; 2)光合和呼吸对干旱具有不同的响应, 干旱时光合的下降比呼吸更为显著, 这导致了NEP的显著下降; 3)当饱和水气压差(Vapor pressure deficit, VPD)达到1.5 kPa以上时, 生态系统的光合、呼吸和净碳吸收均开始下降, 当VPD大于2.5 kPa、土壤相对含水量(土壤含水量/土壤饱和含水量)(Relative soil water content, RSW)低于40%时, 生态系统的碳收支由碳汇转为碳源; 4)土壤干旱是造成碳交换下降的主要驱动因素, 对年NEP下降的平均贡献率为46%, 而大气干旱的贡献率仅为4%。     

中国陆地生态系统通量观测站点空间代表性

涡度相关技术是测定大气与陆地生态系统之间CO₂交换、水分和能量通量最直接的方法,可用于研究土壤、植被与大气间的CO₂交换及其调控机制。收集了11个影响净碳交换量的主要变量信息,包括气象因素、土壤因素和地形因素的非生物因子、实际植被状态以及植被生产力,采用多元地理变量空间聚类分析方法,绘制出不同聚类数（25、50、75、85、100、150和200类）的通量生态区。结合中国现有通量观测站点的空间分布格局,与新生成的通量生态区和已有的自然地理区划进行对比分析,发现由于中国地形复杂,生态系统类型多样,现有85个涡度相关通量观测站点仅能刻画部分中国生态系统类型的净碳交换量时空特征,通量生态区划分为100-150类比较合适。考虑到涡度相关通量观测运行成本,通量站点可增加至150个,从而使得优化后的通量观测网络能够代表中国主要类型的生态系统,并且有利于通量观测数据与遥感资料的有效结合,提高碳水通量观测从站点扩展到区域尺度的精度,从而更好地检验过程机理模型的模拟结果。     

基于模型数据融合的中国温带和亚热带典型森林生态系统碳通量模拟

生态系统碳循环过程对水分响应的研究已成为全球变化关注的焦点问题之一。基于长白山温带针阔混交林与千烟洲亚热带人工针叶林观测站2003—2009年生长季的碳通量(NEE)和气象观测数据,综合考虑水分对光合、呼吸作用的影响,构建不同的NEE模型,并应用模型数据融合方法优化模型参数、遴选最适模型,系统分析了水分因子对不同森林生态系统碳循环的影响。结果表明:(1)优化后的模型参数均能被NEE实测数据较好约束。长白山生长季的光合、呼吸参数值均高于千烟洲,未考虑空气饱和水汽压差(VPD)的模型高估了千烟洲温度敏感性参数(Q10)值、低估了千烟洲基础呼吸速率参数(BR)值;(2)仅考虑VPD对光合作用影响的模型是长白山生长季碳通量模拟的最优模型,但模拟精度提高不显著。不同模型间碳通量组分模拟结果差异较小;(3)考虑VPD和土壤含水量对光合、呼吸作用共同影响的模型是千烟洲生长季碳通量模拟的最优模型,并且显著提高了模拟精度。未考虑水分的模型在生长季高估了总生态系统生产力(GEP)总量2.0%(21.85 g C/m~2),同时更大幅度地高估了生态系统呼吸(RE)总量4.4%(38.02 g C/m~2),从而导致NEE总量低估于实测值7.8%(18.55 g C/m~2)。     

遥感GPP模型在中亚干旱区4个典型生态系统的适用性评价

总初级生产力（GPP）是全球生态系统碳循环的重要组成部分,对全球气候变化有重要影响。目前有多种遥感模型可以模拟总初级生产力,比较不同遥感模型在中亚干旱区上的适用性对推进全球干旱区碳收支估算具有重要意义。基于涡度协相关技术观测的四个地面站数据验证MOD17、VODCA2、VPM、TG、SANIRv五种模型的模拟精度。结果表明：（1）基于光能利用率理论的MOD17、VPM模型模拟咸海荒漠植被和阜康荒漠植被GPP的精度最高（R²分别为0.52和0.80）,但在模拟草地、农田生态系统生产力时存在较明显的低估（RE>20%）;基于植被指数的遥感模型TG模型、SANIRv模型模拟巴尔喀什湖草地生态系统和乌兰乌苏农田生态系统GPP的精度最高（R²分别为0.91和0.81）,同时模拟值与实测值的相对误差也较低;基于微波的VODCA2模型模拟各生态系统生产力的效果最差。（2）水分亏缺是限制植被GPP的主要因素,因此是否合理考虑水分胁迫是影响GPP模型在中亚干旱区适用性的重要因素。研究揭示了遥感GPP模型在中亚干旱区的应用潜力,为推进全球植被碳通量的准确估算提供参考。     

山黧豆幼苗对干旱胁迫的生理响应

用PEG6000模拟干旱对山黧豆进行胁迫。结果发现,胁迫初期（0~48 h）,突然的水分亏缺使气孔导度（Sc）和蒸腾速率（Tr）迅速下降,而净光合速率（Pn）和水分有效利用率（WUE）基本维持稳定。胁迫后期（48~108 h）,上述四个光合指标均不同程度地向对照水平恢复;脯氨酸和相对电导率（REC）受胁迫强度影响较小,但与胁迫时间呈明显正相关,分别在48和60 h增幅最大,而且当REC上升至对照的50%左右时,脯氨酸含量已达对照的14倍之多。由此认为,山黧豆可能主要通过迅速减小Sc和Tr以及大量而急剧的积累脯氨酸来减轻干旱所造成的伤害。     

基于模型数据融合的长白山阔叶红松林碳循环模拟

 充分、有效地利用各种陆地生态系统碳观测数据改善陆地生态系统模型, 是当前我国陆地生态系统碳循环研究领域亟待解决的重要问题之一。该研究以2003~2005年长白山阔叶红松林的6组生物计量观测数据和涡度相关技术测定的碳通量数据为基础, 利用马尔可夫链-蒙特卡罗方法对陆地生态系统模型的关键参数(即碳滞留时间)进行了反演, 进而预测了长白山阔叶红松林生态系统碳库、碳通量及其不确定性。反演结果表明, 长白山阔叶红松林叶凋落物和微生物碳的平均滞留时间最短, 为2~6个月; 其次是叶和细根生物量碳, 二者的平均滞留时间为1~2 a; 慢性土壤有机碳的平均滞留时间为8~16 a; 碳在木质生物量和惰性土壤有机质库中的滞留时间最长, 平均滞留时间分别为77~109 a和409~1 879 a。模拟结果显示, 碳库和累积碳通量模拟值的不确定性将随着模拟时间的延长而增大。当气温升高10%和20%时, 长白山阔叶红松林总初级生产力年总量将分别增加6.5%和9.9%, 净生态系统生产力(NEP)年总量的变化取决于土壤温度的变化。若土壤温度保持不变, NEP年总量将分别增加11.4%~21.9%和17.6%~33.1%; 若土壤温度也相应升高10%和20%, NEP年总量的增幅反而下降甚至低于原来的水平。假设气候和植被保持在2003~2005年的状态, 2020年长白山阔叶红松林NEP年总量为(163±12) g C·m^–2·a^–1, 土壤呼吸年总量为(721±14) g C·m^–2·a^–1。马尔可夫链-蒙特卡罗方法是反演模型参数、优化模拟结果和评估模拟结果不确定性的有效方法, 但今后仍需在惰性土壤碳滞留时间的估计、驱动数据和模型结构的不确定性分析、模型数据融合方法方面进行深入研究, 以进一步提高碳循环模拟的准确性。     

水分胁迫对云杉光合特性的影响

比较研究了3年生红皮云杉和嫩江云杉苗木在不同水分条件下的光合特征,并探讨其对水分胁迫的适应机制。研究结果表明,不同光强下两种云杉的最大光合速率（P_max）随土壤含水量的降低而逐渐下降,且嫩江云杉的降幅较小;表观量子效率（AQY）、蒸腾速率（Tr）、气孔导度（G_s）也都呈下降的趋势。暗呼吸速率（R_d）和水分利用率（WUE）在逐渐干旱的过程中先上升后下降,而胞间CO₂浓度先下降,后随上升。这表明在水分胁迫较为严重时,光合速率的降低主要是非气孔限制造成的。比较上述指标的变化情况可以看出,嫩江云杉对干旱的适应性要强于红皮云杉。     

数据处理方法不确定性对CO_2通量组分估算的影响

基于中国陆地生态系统通量观测研究网络(ChinaFLUX)4个站点(2个森林站和2个草地站)的涡度相关通量观测资料,分析了CO2通量数据处理过程中异常值剔除参数设置、夜间摩擦风速(u*)临界值(u*c)确定及数据插补模型选择对CO2通量组分估算的影响.结果表明:3种数据处理方法均对净生态系统碳交换量(NEE)年总量估算有显著影响,其中u*c确定是影响NEE估算的重要因子;异常值剔除、u*c确定及数据插补模型选择导致NEE年总量估算偏差分别为0.62～21.31 g C.m-2.a-1(0.84%～65.31%)、4.06～30.28 g C.m-2.a-1(3.76%～21.58%)和0.69～27.73 g C.m-2.a-1(0.23%～55.62%),草地生态系统NEE估算对数据处理方法参数设置更敏感;数据处理方法不确定性引起的总生态系统碳交换量和生态系统呼吸年总量估算相对偏差分别为3.88%～11.41%和6.45%～24.91%.     

NEE观测误差分布类型对陆地生态系统机理模型参数估计的影响——以长白山温带阔叶红松林为例

基于观测数据的陆地生态系统模型参数估计有助于提高模型的模拟和预测能力,降低模拟不确定性.在已有参数估计研究中,涡度相关技术测定的净生态系统碳交换量(NEE)数据的随机误差通常被假设为服从零均值的正态分布.然而近年来已有研究表明NEE数据的随机误差更服从双指数分布.为探讨NEE观测误差分布类型的不同选择对陆地生态系统机理模型参数估计以及碳通量模拟结果造成的差异,以长白山温带阔叶红松林为研究区域,采用马尔可夫链-蒙特卡罗方法,利用2003～2005年测定的NEE数据对陆地生态系统机理模型CEVSA2的敏感参数进行估计,对比分析了两种误差分布类型(正态分布和双指数分布)的参数估计结果以及碳通量模拟的差异.结果表明,基于正态观测误差模拟的总初级生产力和生态系统呼吸的年总量分别比基于双指数观测误差的模拟结果高61～86 g C m-2 a-1和107～116 g C m-2 a-1,导致前者模拟的NEE年总量较后者低29～47 g C m-2 a-1,特别在生长旺季期间有明显低估.在参数估计研究中,不能忽略观测误差的分布类型以及相应的目标函数的选择,它们的不合理设置可能对参数估计以及模拟结果产生较大影响.     

干旱胁迫下藜的光合特性研究

通过人工控制水分模拟干旱来研究生长期的藜对干旱胁迫的生理生化反应,以期望为干旱农业的高效生产提供理论依据。以盆栽的藜为材料,用称重控制浇水的方法,研究了干旱胁迫对藜叶片的光合特性。结果表明：干旱胁迫下藜的光合日变化呈双峰型,有“午休”现象（13:00）且受气孔限制;最大净光合速率出现在上午8:00。与正常条件下生长的藜相比,干旱胁迫下藜的光饱和点(LSP)、最大净光合速率(Pn)、表观量子效率(AQY)、二氧化碳饱和点(CSP)和羧化效率(CE)均降低,分别为1 200 μmolphoton·m^-2·s^-1、8.01 μmol CO₂·m^-2·s^-1、0.016 1 μmol CO₂·mol^-1 photons、1 200 μmol CO₂·mol^-1、0.017 6 μmol CO₂·m^-2·s^-1;光补偿点(LCP)、二氧化碳补偿点(CCP)升高,分别达到44.88 μmol photon·m^-2·s^-1、和46 μmol CO₂·mol^-1,干旱使藜的光合能力下降。干旱胁迫下藜的光合能力虽有所下降,但与其它C3植物相比仍具有较强的CO₂同化能力。藜是一种耐旱力较强的植物。     

中亚热带人工针叶林生态系统碳通量拆分差异分析

涡度通量观测可直接获取陆地生态系统与大气之间CO2净交换量(NEE),但深入认识碳循环过程和校验生态系统模型需要不同时间尺度总初级生产力(GPP)和生态系统呼吸(Re)等碳通量数据。利用中国陆地生态系统通量观测与研究网络(ChinaFLUX)中亚热带人工针叶林生态系统2003—2009年的涡度通量和气象观测数据,分析了两种NEE拆分方法对不同时间尺度GPP和Re评估的影响,结果表明:(1)两种拆分方法得到的生态系统碳通量组分(GPP和Re)的季节动态变化一致,都在生长季7、8月份达到峰值;(2)非线性回归模型拆分得到的全年Re和GPP相较于光响应曲线模型分别高出2%—28.6%和1.6%—23%,最大高出317.6 gC·m-2·a-1(2006年),逐月最大差值主要发生在8、9月份;(3)不同时间尺度上,两种方法拆分得到的GPP和Re之间差值的环境响应因子不同。在广泛采用非线性回归模型进行拆分时,如果当月光合有效辐射接近到905mol·m-2·月-1,月平均空气饱和水汽压差接近1.18 kPa时,需要考虑使用光响应曲线模型拆分该月通量,结合两种拆分方法以减小全年的误差。     

三江并流核心区亚高山森林非生长季净生态系统CO2交换量及其影响因素

采用开路式涡度相关系统,针对三江并流核心区西藏红拉山滇金丝猴国家自然保护区,通过测量和分析非生长季亚高山常绿针叶林净生态系统碳交换量(NEE),探讨了亚高山森林非生长季CO₂通量特征及其主要影响因子。保护区常绿针叶林NEE值在非生长季具有明显“U”型变化曲线,白天表现为碳吸收,夜间表现为碳释放,日间CO₂吸收高峰介于12:00到15:00之间,平均每天碳汇时间在10 h左右。非生长季各月NEE大小依次为：4月> 3月>2月>11月>1月>12月。研究期内气温(T)、相对湿度(RH)、饱和水汽压差(VPD)和光合有效辐射(PAR)等气象因子对净生态系统CO₂交换量影响显著。此外,森林碳吸收对温度响应敏感,光合作用在整个非生长季较为明显。各影响因子中光合有效辐射对碳交换影响最大;夜间NEE与5 cm土壤温度呈极显著相关性(P<0.01)且NEE随着土壤温度升高而增大;整个非生长季NEE、生态系统呼吸量(Re)和总生态系统CO₂交换量(GEE)分别为-596.759 g...     

内蒙古克氏针茅草原生态系统-大气通量交换特征

基于内蒙古克氏针茅草原生态系统的涡动相关观测资料和小气候梯度系统观测资料,对研究区碳通量和水热通量的日、季动态进行了研究.结果表明：克氏针茅草原生态系统生长季的碳通量日动态呈U型曲线,即日出前释放CO₂,日出后开始吸收CO₂,正午前后达峰值,午后吸收CO₂减弱,日落后重新转为释放CO₂;9月白天CO₂吸收最为强烈,8月次之,10月最低.克氏针茅草原的感热和潜热通量的日动态均呈倒U型曲线,与碳通量日动态相反,即白天感热和潜热通量多为正值,夜间感热为负值,潜热接近于零;感热通量以5月最高,潜热通量以9月最高.冬季草地为弱碳源,CO₂通量较小,夏季表现为明显的碳汇.     

Fisher discriminant analysis method applied in drought detection: an instance in an alpine meadow ecosystem

Fisher判别法在识别干旱中的应用——以高寒草甸生态系统为例 Fisher判别法可综合考虑事物的多重属性特征,进而分辨事物类型。若能将其应用于对干旱的识别,或将成为一种准确而有效的干旱识别方法。本研究以高寒草甸生态系统为研究对象,基于9年碳通量和小气候观测数据,以土壤含水量(soil water content, SWC)和饱和水汽压差(vapor pressure deficit, VPD)作为判别因子,利用Fisher判别法识别干旱。Fisher判别法可用于对高寒草甸生态系统干旱的识别。因其综合考虑了土壤水分匮缺和大气水分匮缺,故可以更合理准确地判断干旱的开始和结束时间。基于干旱样本和非干旱样本的特征,构建判别方程为：y = 24.46SWC − 4.60VPD。当y > 1时,样本点位于临界线上方,若持续10天以上,则判定为发生干旱。两次干旱过程间隔2天以内,可认为是一次干旱过程。随着研究年限的增加和观测数据的积累,该临界线方程尚有优化空间,其对干旱识别的准确度可进一步提高。     

华北平原玉米田生态系统光合作用特征及影响因素

采用涡度相关法对华北平原夏玉米田进行了连续4a(2003-2006年)的碳通量观测,结果表明:夏玉米田生态系统初始量子效率(α)、最大光合速率(P_max)、暗呼吸速率(R_d)和总初级生产力(GPP)随作物生长发育而变化。在夏玉米生育前期和后期,α、P_max、R_d和GPP都比较小,其最大值出现在抽穗期/灌浆期。2003-2006年,夏玉米生长季平均α、P_max、R_d的范围分别为0.054-0.124 μmol/μmol、1.72-2.93 mg CO₂ · m^-2 · s^-1、0.23-0.38 mg CO₂ · m^-2 · s^-1。α、P_max和R_d均随叶面积指数(LAI)增加呈指数增长。2003-2006年夏玉米生长季GPP总量分别为806.2、741.5、703.0、817.4 g C/m²,年际差异较大。玉米田生态系统GPP随温度升高呈指数增长。在玉米营养生长阶段,GPP随LAI增加而增大,两者之间的关系可用直角双曲线方程来表示;生殖生长阶段,GPP随LAI降低而下降.相同LAI下,生殖生长阶段的GPP明显低于营养生长阶段。     

水分胁迫对脂松幼苗叶绿素荧光特征的影响

以2年生盆栽脂松苗木为材料,研究了不同土壤含水量（干旱和水渍胁迫及解除胁迫后）对脂松幼苗的叶绿素荧光参数的影响,旨在为脂松的扩大栽培和管理提供理论依据。结果表明：（1）土壤水分胁迫对脂松苗木针叶初始荧光(F_o)、最大荧光(F_m)、可变荧光（F_v）、PSⅡ最大光化学效率(F_v/F_m)、PSⅡ电子传递量子产率（Φ_PSⅡ）有显著影响（p<0.05）;（2）干旱和淹水胁迫中脂松针叶光化学效率（F_v/F_m）均表现为随胁迫强度增大和时间的延长而降低,总体表现为FC60>FC80>FC>FC40>FC20>FL,解除胁迫后有不同程度的升高,且干旱处理比水淹处理恢复快。（3）60%土壤含水量叶绿素荧光参数水平最佳,据此认为中等轻度干旱立地较适宜脂松苗木生长。     

动态及对气候变化响应的模拟分析

水分条件不仅影响半干旱区群落的组成, 而且在一定程度上决定了群落的功能。处于不同水分条件生境下群落的优势物种在水分利用和同化物利用效率方面的功能特征会存在差异, 这些差异将导致群落对于气候变化产生不同的响应, 进而影响到景观和区域尺度上对于全球变化下碳动态和格局的分析。该文选取了锡林河流域典型草原区沿水分梯度的4个代表群落, 在野外实验测定并结合长期定位研究成果基础上, 利用BIOME-BGC模型对代表群落的长期净初级生产力(Net primary productivity,NPP)动态进行了模拟和模型验证。 通过分析该地区1953～2005年气候变化趋势, 推测了未来可能的气候变化情景, 进而模拟了气候变化下4个群落长期NPP动态的响应。结果表明,当前气候条件下, 羊草(Leymus chinensis)群落NPP平均值为197.76 gC·m^-2 (SE=7.11), 大针茅(Stipa grandis)群落NPP平均值为198.95 gC·m^-2 (SE=6.41), 贝加尔针茅(Stipa baicalensis)群落NPP平均值为210.41 gC·m^-2(SE=7.87), 克氏针茅(Stipa krylovii)群落NPP平均值为144.92 gC·m^{- 2} (SE=4.64), 4个群落NPP平均值为188.01 gC·m^-2 (SE=3.72); 气候变化情景下, 温度增加下(P0T1),NPP平均下降14.2%,降水增加下(P1T0), NPP平均增加13.2%,温度与降水都增加情景下(P1T1), NPP平均下降2 .7%, 但由于生境水分条件差别和优势物种功能特征差异, 4个群落表现出了增减幅度不同的趋势。对气候因子的敏感性分析及回归分析表明, 降水是该地区NPP最主要的决定因子, 而温度决定作用相对较小,主要通过影响植物的呼吸和水分蒸散等过程影响NPP。在最有可能代表未来气候 变化的温度增加的两种情景下(P0T1、P1T1), NPP均呈下降趋势。群落NPP对气候变化的响应趋势与水分胁迫系数(Water stress index, WSI)、碳胁迫系数(Carbon stress index, CSI)变化密切相关。克氏针茅群落由于所处生境水分条件差,WSI高,对降水的依赖程度最大;贝加尔针茅群落一方面处于较好的水分生境,具有较小的WSI,另一方面,由于具有高碳氮比,维持呼吸消耗的光合产物比例低,CSI远低于其它3个群落, 未来气候变化下, NPP较其它3个群落仍较高。