全 文 ：作物学报 ACTA AGRONOMICA SINICA 2015, 41(4): 601612 http://zwxb.chinacrops.org/
ISSN 0496-3490; CODEN TSHPA9 E-mail: xbzw@chinajournal.net.cn

本研究由国家重点基础研究发展计划(973 计划)项目(2012CB417103), 国家自然科学基金青年项目(31400418), 国家公益性行业(农业)
科研专项(201303022)和中国博士后科学基金面上项目(2014M561044)资助。
 通讯作者(Corresponding author): 徐明, E-mail: mingxu@igsnrr.ac.cn
第一作者联系方式: E-mail: zhengyunpu_000@sina.com
Received(收稿日期): 2014-11-03; Accepted(接受日期): 2015-02-06; Published online(网络出版日期): 2015-03-03.
URL: http://www.cnki.net/kcms/detail/11.1809.S.20150303.1649.001.html
DOI: 10.3724/SP.J.1006.2015.00601
玉米叶片气孔特征及气体交换过程对气候变暖的响应
郑云普 1,2,4 徐 明 2,* 王建书 3 邱 帅 2 王贺新 4
1河北工程大学水电学院, 河北邯郸 056038; 2中国科学院地理科学与资源研究所生态系统网络观测与模拟重点实验室, 北京 100101;
3河北工程大学农学院, 河北邯郸 056038; 4大连大学现代农业研究院, 辽宁大连 116622
摘 要: 气孔是植物叶片表面控制大气与植物间气体交换的孔状结构, 对于生态系统碳、水循环过程的调节起着非
常重要的作用。本文利用典型农田生态系统实验增温平台, 研究了未来气候变暖对玉米叶片的气孔特征(包括气孔频
度、气孔开口大小和形状以及气孔分布格局)和气体交换过程的影响。结果表明: (1)尽管增温并没有改变气孔密度
(P>0.05), 但却由于表皮细胞数目的减少导致气孔指数显著增加 12% (P<0.05); (2)增温使气孔开口的长度显著减小
18% (P<0.01), 宽度增加 26% (P<0.01), 面积和周长分别增加 31% (P <0.01)和 13% (P<0.05); (3)实验增温还使单个气
孔之间最近邻域的平均距离显著增加, 表明气孔在玉米叶片上的分布变得更加均匀; (4)增温导致玉米叶片的净光合
反应速率(Pn)、气孔导度(Gs)和蒸腾速率(Tr)分别增加 52% (P<0.05)、163% (P<0.001)和 81% (P<0.05); 与此相反, 玉
米叶片的暗呼吸速率(Rd)却显著降低 24% (P<0.01)。增温没有对细胞间 CO2浓度(Ci)和水分利用效率(WUE)产生显著
的影响(P>0.05)。本研究结果表明, 未来全球气候变暖可能通过改变玉米叶片的气孔频度、气孔开口大小和形状及其
在叶片上的空间分布格局来改变其气体交换过程。
关键词: 全球变暖; 玉米; 气孔特征; 气体交换; 华北平原
Responses of the Stomatal Traits and Gas Exchange of Maize Leaves to Climate
Warming
ZHENG Yun-Pu1,2,4, XU Ming2,*, WANG Jian-Shu3, QIU Shuai2, and WANG He-Xin4
1 School of Water Conservancy and Hydropower, Hebei University of Engineering, Handan 056038, China; 2 Key Laboratory of Ecosystem Network
Observation and Modeling, Institute of Geographical Sciences and Natural Resources, Chinese Academy of Sciences, Beijing 100101, China; 3 School
of Agriculture, Hebei University of Engineering, Handan 056038, China; 4 Institute of Modern Agricultural Research, Dalian University, Dalian
116622, China
Abstract: Stomata are the pores on leaf surfaces controlling gas exchanges, mainly CO2 and water vapor, between the atmosphere
and plants, and thus regulate carbon and water cycles in various ecosystems. This study investigated the effects of experimental
warming on the stomatal frequency, stomatal aperture size and shape, and stomatal distribution pattern, and their relationships
with the leaf gas exchange rates of maize (Zea may L.) leaves through a field manipulative warming experiment with infrared
heaters in a typical agriculture ecosystem in the North China Plain. Our results showed that experimental warming had little effect
on stomatal density, but increased stomatal index by 12% (P<0.05) due to the reduction in the number of epidermal cells under the
warming treatment. Warming also decreased stomatal aperture length by 18% (P<0.01) and increased stomatal aperture width
26% (P<0.01). As a result, experimental warming increased the average stomatal aperture area by 31% (P<0.01) and stomatal
aperture circumference by 13% (P<0.05), and resulted in a more regular stomatal distribution on both the adaxial and abaxial sur-
faces in leaves with an increased average nearest neighbor distance between stomata. In addition, experimental warming also af-
facted the gas exchange of maize leaves. Experimental warming significantly increased net photosynthetic rate (Pn), stomatal
conductance (Gs), and transpiration rate (Tr) by 52% (P<0.05), 163% (P<0.001), and 81% (P<0.05), respectively. Meanwhile,
experimental warming decreased the leaf dark respiration (Rd) by 24% (P<0.01), but had no significant effects on intercellular
602 作 物 学 报 第 41卷


CO2 concentration (Ci) and water use efficiency (WUE; P>0.05). In conclusion, the experimental warming may affect the gas
exchange of maize leaves through the changes of the stomatal traits including stomatal frequency, stomatal aperture size and shape,
and stomatal distribution on leaves.
Keywords: Global warming; Maize; Stomatal traits; Gas exchange; The North China Plain
气孔是植物叶片表面控制大气与植物间进行
CO2、水蒸气等气体交换的孔状结构[1-2], 对于调节
生态系统碳、水循环过程起着非常重要的作用[3-5]。
就全球范围而言, 通过这种叶片与大气之间的气体
交换过程植物每年大约在光合作用中固定440×1015
g CO2和蒸发32×1018 g H2O[6-7]。植物叶片通过改变
气孔频度(气孔密度及气孔指数)、气孔开度、气孔大
小和气孔空间分布格局等特征来最优化其气体交换
效率。植物叶片的气体交换效率量化的气孔导度来
表征, 这主要由气孔大小、形状、频度以及空间分
布格局所决定[8]。然而, 目前有关增温对叶片气孔频
度及大小等方面特征的影响尚无定论。有些研究发
现, 增温对植物叶片气孔密度和气孔指数均没有影
响 [9-12]; 然而 , 还有一些研究却发现增温可能导致
气孔密度和气孔指数的减少或增加[13-17]。另据报道,
增温还能改变单个气孔的大小和形状 [14,18]。例如 ,
Ferris等 [14]的研究结果显示 , 增温显著增加了多年
生黑麦草(Lolium perenne)叶片的气孔长度 ; 然而 ,
张立荣等 [18]的研究却发现增温使我国青藏高原4种
亚高山草甸物种的叶片气孔长度减小。除了叶片气
孔的数量、大小和形状之外, 增温还可能通过对细
胞分裂和分化过程的影响来改变气孔在叶片上的空
间分布格局[19-21]。同时, 细胞的分裂及分化过程也
会受到遗传信号[22-23]以及环境因素[24-25]的调控和影
响。另外, 气孔在植物叶片上的空间分布格局在物
种间也呈现出较高的变异性[14-15]。气孔在植物叶片
上空间分布格局的变异性通常表现在多个尺度上 ,
例如近/远轴面的比例、同一叶面不同部位的变异、
同一叶片单个气孔的变异[19-21]。早期的研究报道叶片
近轴面和远轴面的气孔密度显著不同[14-15,19,21,26-28]; 气
孔的空间分布格局在同一叶片的不同部位也呈现出
较高的变异性[14,16,29-33]。此外, 叶片近轴面和远轴面
的气孔密度也存在着显著的差异[14-15,19,26-28]。
华北平原是我国最重要的粮食产区之一, 种植
制度以玉米和小麦轮作系统为主 , 每年全国大约
50%的小麦和30%的玉米产自该区[34]。值得关注的
是, 许多气候模型已经证实我国的气候正在逐渐呈
现变暖的趋势, 特别是自从20世纪70年代以后我国
北部地区(如华北平原)的气候呈现出年平均气温升
高的特征[35]。这种区域甚至全球范围内的气候变暖
势必改变粮食作物的蒸散和水分利用效率[36], 并进
一步改变农作物的生长速率和粮食产量 [37]。此外 ,
近年来的许多研究结果还表明, 全球变暖也会通过
改变农作物的生长季而对其物候期产生深远的影
响。例如, Tao等[38]的研究表明, 全球气候变暖导致
作物产量形成的重要生长阶段(从开花到成熟)呈现
延长的趋势[38]。Hou等[39]在华北平原的增温实验研
究发现, 气候变暖将会缩短小麦的生长季, 但并不
改变其生育期, 最终导致小麦产量减少。同样, 以往
的研究结果已经发现气候变暖导致全球许多地区玉
米产量的减少[40]。由此看来, 气候变暖对农作物产
量的影响将会直接关系到未来区域 /全球粮食供应
的安全问题。然而, 目前全球气候变暖对作物产量
产生影响的机制还不清楚。此外, 由于植物叶片表
面分布的气孔是植物与大气之间进行气体交换的重
要器官[1-2], 它不仅控制着植物的光合、呼吸及蒸腾
等关键的生理过程[3-5], 而且还最终决定着农作物的
产量[37]。因此, 深入探讨作物叶片气孔结构特征及
其气体交换过程对实验增温的响应有助于全面理解
未来全球变暖对作物产量影响的潜在机制。
本研究通过开展野外增温控制实验, 利用光学
显微镜、扫描电镜、空间分析技术等手段, 研究实
验增温对华北平原重要作物玉米叶片的气孔特征
(气孔频度、气孔开口大小和形状以及气孔的空间分
布格局)及其气体交换过程的影响。这将对深入理解
实验增温对玉米气孔特征影响机制具有重要的理论
意义, 为预测未来气候变暖对农田生态系统结构和
功能的影响提供重要的数据和理论支持。
1 材料与方法
1.1 研究区概况
中国科学院禹城综合试验站(36°57′N, 116°38′E,
海拔 23.4 m)位于山东省禹城市, 属大陆性季风气候
暖温带地区。年平均气温 13.1℃, 1月平均温度3℃,
7 月平均温度 26.9℃; 年平均降水量 610 mm, 降水
季节分配不均匀, 3 月至 5 月平均降水量 75.7 mm,
占年的 12.4%, 6 月至 8 月降水量平均为 419.7 mm,
占年的 68.8%。该实验站土壤质地以粉沙和轻壤为
第 4期 郑云普等: 玉米叶片气孔特征及气体交换过程对气候变暖的响应 603


主 ; pH 值为 7.9~8.0; 土壤有机质含量较低 , 为
0.6%~1.0%; 全氮含量为 0.050%~0.065%。从土壤自
然条件看, 该站在华北平原具有典型代表性, 当地
种植制度以冬小麦-夏玉米轮作为主。
1.2 野外原位增温实验
在2009年9月份, 随机选取试验区内6个2 m ×
3 m的样方用以布置增温实验。随机选取其中的3个
样方作为增温组, 另3个样方为对照组。选取的3个
增温样方从2009年11月18日开始利用悬挂于样方正
上方距离地表2.25 m高度的MSR-2420型红外发射器
(165.0 cm × 1.5 cm, Kalglo Electronics Inc., Bethle-
hem, PA)进行全天24 h不间断的连续增温处理(图1)。
每个对照样方内同样在距离地面2.25 m高度处悬挂
一个相同大小和形状的“伪”红外增温装置模拟红外
发射器对下方土壤产生阴影的效果。为了避免红外
增温装置影响到对照样方, 增温样方和对照样方相
距约5 m。于2011年6月24日将玉米种子(Zea may L.
cv. Zhengdan 958)分别播种到增温和对照样地, 同
时进行连续不间断的增温处理(全天24 h增温处理),
利用PT100热电偶系统 (Unism Technologies Incor-
porated, Beijing, China)自动记录每小时2.25 m高度
处的空气温度和5 cm深处的土壤温度。另外用红外
温度计(FLUKE 574, Fluke Inc., USA)测量玉米冠层
叶片表面的温度。在整个玉米生长季从2011年6月24
日到10月7日, 相对于对照样地, 增温样地的空气温
度和土壤温度在日间分别增加了 (1.42±0.18)℃和
(1.68±0.90)℃; 同时在夜间分别增加了(1.77±0.24)℃
和(2.04±0.16)℃。增温样地的灌层温度相对于对照
样地在日间增加(2.08±0.72)℃。为避免干旱胁迫对
玉米植株生长的影响, 在整个生长季利用地下水对
增温和对照的各个样地内植株定期灌溉。同时, 利
用FDS100土壤湿度传感器(Unism Technologies In-
corporated, Beijing, China)监测样方内0~10 cm深度
土壤湿度的变化状况。在整个监测时期, 对照样地
10 cm深处土壤湿度为26%, 而增温样地10 cm深处
土壤湿度为25%[39]。
1.3 印迹法取样
为了确定玉米叶片最大的气孔开度, 气孔印迹
样品的采集工作在气孔开张的最适环境条件下开
展。2011年8月24日上午10:30—11:00的天气晴朗 ,
空气温度约为32℃, 其自然环境条件对于气孔的开
张是最理想的, 有理由认为此刻气孔具有最大的开
张度。从3个对照样地和3个增温样地中各随机选取5
棵玉米植株采样。由于玉米的穗位叶主要决定着玉
米的产量, 本研究从上述选定的10株玉米中随机选
取1片穗位叶(10个叶片)利用透明无色的指甲油分别
从近轴面和远轴面的3个不同部位(基部、中部和尖
部)采集面积为5 mm × 15 mm的气孔印迹样品, 共
60个样品(10个叶片 × 3个部位 × 2个不同叶面 =
60个)用于气孔特征的观察和分析。
1.4 光学显微镜观察
将上述叶片印迹玻片在装备有照相机 (DFC
300-FX, Leica Corp, Germany)的莱卡光学显微镜
(DM2500, Leica Corp, Germany)下观察并拍照。随机
选择 5个视野, 每个视野拍 3张照片, 即在玉米叶片
远轴面和近轴面的不同部位分别得到 15张气孔的显
微照片(面积为 0.30 mm × 0.45 mm)。然后从其中随
机选取 5 张计算气孔密度、表皮细胞密度及气孔指
数, 即分别从玉米叶片的近轴面和远轴面得到 75个
气孔印迹显微照片(面积为 0.30 mm × 0.45 mm)。再
将照片综合到一起, 探讨增温对整个叶片气孔特征
的影响。气孔和表皮细胞的密度用单位叶面积内气
孔和表皮细胞的个数来表示。气孔指数＝气孔数/(气
孔数+表皮细胞数)×100%。另外, 分别在上述每张照
片上选取 6 个气孔和 6 个表皮细胞(共计为 75 个照
片×6 个气孔/表皮细胞=450 个气孔或表皮细胞), 利
用 Image J 分析软件(NIH, USA)分别测量气孔和表
皮细胞的长度、宽度、周长及面积[50-51]。分别将这
些数据综合在一起来研究叶片近轴面和远轴面气孔
和表皮细胞特征对增温响应的差异。此外, 利用上
述叶片近轴面和远轴面气孔和表皮细胞数据的平均
值来代表整个叶片气孔和表皮细胞的特征。
1.5 扫描电镜观察
从上述每个处理叶子的中部采集3个2 mm ×
2 mm的叶片用于气孔及表皮细胞特征的扫描电镜
观察。野外采集下来的玉米叶片立即利用2.5% (v/v)
的戊二醛(0.1 mol L–1磷酸缓冲液, pH 7.0)固定样品,
放到4℃环境下冷藏保存 , 并迅速将样品运回实验
室, 再用相同浓度的磷酸缓冲液冲洗6次, 并在室温
环境下用1% (v/v)锇酸固定3 h, 用同样的磷酸缓冲
液冲洗干净。将该叶片组织经一系列浓度的酒精梯
度脱水 , 进行临界点干燥 , 固定在观察台上 , 利用
高压涂膜装置喷金处理。最后, 用Quanta 200扫描电
子显微镜(FEI Corp, USA)观察和照相。
1.6 气孔的空间格局分析
分别从增温及对照叶片的扫描电镜照片中各随
604 作 物 学 报 第 41卷



图 1 华北平原典型农田生态系统实验增温样地
Fig. 1 Experimental warming site in a typical agricultural ecosystem in the China Northern Plain

机选取3张(放大100倍)用于研究增温对叶片气孔空
间格局的影响。在本项分析中, 认为每一个气孔都
是玉米表面上分布的单点, 气孔开口的最中间位置
是该单点的位置。首先利用空间分析软件ArcGIS
10.0 (ESRI Inc. USA)将所选的显微照片在相同的坐
标系下进行数字化处理 , 得到这3个叶片近轴面和
离轴面上每一个气孔的坐标。然后, 利用空间统计
分析方法Ripley’s K-function对数字化后表征气孔分
布状况的点进行空间分析[41]。Ripley’s K-function是
一个分布累加函数, 该函数利用所有单点距离的二
阶矩阵探究这些点在不同尺度上的两维分布格局。
分析结果由L(t)值来表达。
L(t) = ( ) / πK t −t (1)
当该分布格局为泊松随机分布时, 所有的 t 值
到 L(t)的距离均相等。为了确定 95%的可信任区间,
采用蒙特卡洛算法模拟随机分布点 1000次。假如叶
片表面的气孔在给定尺度 t 下为随机分布, 则计算
出来的 L(t)值应该位于 95%可信任区间之内。假如
L(t)值大于 95%可信任区间, 则气孔在该尺度下为
簇状分布。若 L(t)值小于 95%可信任区间时, 气孔在
该尺度下为规则分布。
1.7 玉米叶片气体交换参数测量
利用 Li-6400便携式光合测定系统(LI-COR Inc.
Lincoln, Nebraska, USA)测定玉米叶片的净光合速
率、气孔导度和蒸腾速率等气体交换参数。该测定
系统配套使用的 2 cm × 3 cm标准气室能独立控制光
合光量子通量密度 (Photosynthetic Photon Flux
Density, PPFD)、CO2浓度、叶片温度及气室湿度。
测量时设定叶室内的光照强度为 1500 μmol m–2 s–1
光合有效辐射 (Photosynthetically Active Radiation,
PAR), 水蒸汽压亏缺(Vapor Pressure Deficit, VPD)控
制为 2.0 kPa, CO2浓度为 380 μmol mol–1, 叶片温度
为 25℃。
1.8 数据统计
利用单因素和多因素方差分析方法分析增温对
气孔特征的影响, 再用 Duncan’s multiple range test
比较处理间的显著性差异(P<0.05)。本研究中的统计
分析均利用 SPSS 13.0 (Chicago, IL)统计软件完成,
利用软件 Sigmaplot作图。
2 结果与分析
2.1 实验增温对气孔密度及气孔指数的影响
增温显著增加了叶片的气孔指数(P<0.05), 但
并没有改变气孔密度(表1)。增温增加叶片气孔数约
5%, 即由单位叶面积上67个增加为70个。然而, 增
温环境下的气孔指数显著高于对照12% (P<0.05)。
第 4期 郑云普等: 玉米叶片气孔特征及气体交换过程对气候变暖的响应 605


2.2 实验增温对气孔开度和形状的影响
实验增温不但影响气孔密度和指数, 还影响气
孔的张开度。增温显著减小气孔长度, 但增加气孔
宽度。由于温度的升高叶片气孔长度显著减小了
18% (P<0.01; 表 1); 而气孔宽度却显著增加了 26%
(P<0.01; 表 1)。上述结果同直接观察的结果相一致。
通过光学显微镜和扫描电镜观察到增温条件下的玉
米叶片与对照叶片相比具有更短、更宽的气孔开口
(图 2和图 3)。另外, 增温分别增加了气孔的面积和周
长 31% (P<0.01)和 13% (P<0.05), 还导致气孔面积指数
(单位叶面积上的总气孔面积)显著升高 40%, 但是对于
气孔的形状指数没有显著的影响(P>0.05; 表 1)。

表 1 增温对玉米叶片气孔结构特征的影响
Table 1 Effects of experimental warming on stomatal traits of maize leaves
对照 Control 增温 Warming 气孔参数
Stomatal parameter 近轴面 Adaxial 远轴面 Abaxial 近轴面 Adaxial 远轴面 Abaxial
增量
Increase (%)
P值
P-value
56±5 b 77±3 a 58±9 b 81±1 a — — 气孔密度 SD
(number mm–2) 67 70 4.5 P>0.05
13.8±0.1 c 19.4±0.5 ab 16.5±0.3 b 20.6±0.8 a — — 气孔指数 SI
(%) 16.6 18.6 11.7 P<0.05
36.8±5.1 a 35.5±2.1 a 30.8±3.8 b 28.5±3.7 b — — 气孔长度 SAL
(μm) 36.2 29.7 18.0 P<0.01
3.7±1.0 bc 3.2±0.3 c 4.2±1.0 ab 4.5±0.9 a — — 气孔宽度 SAW
(μm) 3.5 4.4 25.7 P<0.01
118±43 bc 100±15 c 150±45 a 135±31 ab — — 气孔面积 SAA
(μm2) 109 143 31 P<0.01
74±10 b 70±3 b 85±10 a 77±7 b — — 气孔周长 SAC
(μm) 72 81 12.5 P<0.05
0.66±0.24 c 0.77±0.12 bc 0.87±0.26 b 1.10±0.25 a — — 气孔面积指数 SAI
(%) 0.71 0.99 39.9 P<0.01
13.8±1.0 b 14.3±0.7 b 14.4±0.9 b 15.2±0.9 a — — 气孔形状指数 SASI
(%) 14.0 14.8 5.7 P>0.05
— — 远/近轴面 SD比率
Adaxial/abaxial
0.73±0.09 0.72±0.13
–1.4 P>0.05
— — 远/近轴面 SI比率
Adaxial/abaxial
0.71±0.09 0.80±0.15
12.3 P<0.05
所有数据为平均值±标准偏差, 气孔密度、气孔指数、气孔面积指数、气孔形状指数为 n = 75, 气孔长度、气孔宽度、气孔面积、气孔
周长为 n = 450。所有数据的平均值利用单因素方差分析在 0.05水平上进行比较。同列标以不同字母的值在 0.05水平上差异显著。SD: 气
孔密度; SI: 气孔指数; SAL: 气孔长度; SAW: 气孔宽度; SAA: 气孔面积; SAC: 气孔周长; SAI: 气孔面积指数; SASI: 气孔形状指数。
Values are means ± standard deviation for SD, SI, SAAI, and SASI (75 subsamples and 3 replicates), and for SAL, SAW, SAA, and
SAC (450 subsamples and 3 replicates). Means are compared by the one-way analysis of variance (ANOVA) at P<0.05. Values followed by
different letters are significantly different at P<0.05 and the same letters indicate P>0.05. SD: stomatal density; SI: stomatal index; SAL:
stomatal aperture length; SAW: stomatal aperture width; SAA: stomatal aperture area; SAC: stomatal aperture circumference; SAI: stomatal
area index; SASI: stomatal shape index.

图 2 玉米叶片气孔和表皮细胞的光学显微照片
Fig. 2 Optical microscopic photographs of stomata and epidermal cells in maize leaves
a: 对照; b: 增温; 标尺: 100 μm。a: control; b: warming; Bar: 100 μm.
606 作 物 学 报 第 41卷



图 3 玉米叶片气孔和表皮细胞的扫描电子显微照片
Fig. 3 Scanning electron microscopic photographs of stomata
and epidermal cells in maize leaves
a、b、c为对照环境下的玉米叶片气孔特征, d、e、f为增温环境
下的气孔特征。与对照相比, 增温使近轴面(b和 e)和远轴面(c
和 f)的气孔均变得更短、更宽。同时, 增温还增加了表皮细胞的
长度和宽度。标尺: 10 μm (a和 d)、40 μm (b、c、e和 f)。
Shorter and wider stomatal pores were observed on both the ad-
axial surface ( band e) and abaxial surface (c and f) of maize
leaves grown at elevated temperature than those grown at ambient
temperature. In addition, elevated temperature also increased the
width of epidermal cells. Bars, 10 μm (a and d) and 40 μm
(b, c, e, and f).

2.3 实验增温对气孔空间分布格局的影响
增温几乎没有对气孔密度的近轴面 /远轴面比
例产生影响(对照与处理的近轴面/远轴面气孔密度
比例分别为 0.72和 0.73), 但却导致气孔指数的近轴
面/远轴面比例约 12%由 0.71显著增加到 0.80 (表 1)。
另外, 增温对玉米叶片不同部位(尖部、中部和基部)
的气孔特征也产生了影响 ,增加了近轴面基部和远
轴面尖部的气孔密度(P<0.05; 表 2)。同时, 增温却
显著减少了近轴面基部和中部以及远轴面中部的表
皮细胞密度(P<0.05)。因此, 相对于对照, 增温使近
轴面基部和中部气孔指数分别增加约 40%和 25%,
远轴面中部气孔指数约增加 10% (P<0.05)。另外, 增
温还分别使近轴面尖部和远轴面中部的气孔长度减
小约 30%和 35% (P<0.05)。相对于对照组, 增温分
别导致叶片近轴面尖部和中部气孔宽度增加 20%和
40% (P<0.05)。相似地, 增温也同时分别使远轴面尖
部和中部的气孔宽度增加约 30% (P<0.05)。对于近
轴面, 增温使叶片中部和基部的气孔面积显著增加
34%和 33%及气孔面积指数增加了 28%和 61%
(P<0.05)。增温分别使远轴面尖部、中部及基部的气
孔面积增加 43%、40%、30%以及使气孔面积指数增
加 61%、44%、32% (P<0.05)。另外, 多因素方差分
析的结果显示, 增温和叶片位置显著改变了玉米叶
片的气孔密度、气孔指数、气孔长度、气孔宽度、
气孔面积及气孔面积指数(P<0.05)。
无论是对照还是增温环境条件下玉米叶片近轴
面和远轴面气孔均在小尺度范围内为规则分布
(<140 μm), 而在大尺度范围内呈随机分布特征(图 3
和图 4)。然而, 远轴面的气孔分布比近轴面的气孔
分布更加规则, 因为在相同尺度下远轴面的 L(t)值
更小, 尤其是在对照条件下(图 4和图 5)。对照温度
下玉米叶片气孔最规则的分布格局发生在尺度约
25 μm处, 此时近轴面平均的最小 L(t)为1.81, 而远
轴面平均的最小 L(t)为3.25。增温使玉米叶片近轴
面和远轴面气孔的分布格局更加规则, 这主要体现
在增温环境下玉米叶片的 L(t)值比对照环境变得更
小, 即近轴面平均的最小 L(t)从1.81 减小为4.43
(图 4), 远轴面平均的最小 L(t)从3.25 减小为4.80
(图 5)。另外, 增温还增加了玉米叶片气孔规则分布
的尺度范围, 最规则分布的尺度约为 60 μm, 而对照
条件下仅为 25 μm。
2.4 实验增温对叶片气体交换过程的影响
实验增温不仅改变了玉米叶片的气孔特征, 还
对其气体交换过程产生了显著的影响(图 6)。增温导
致玉米叶片的净光合速率(Pn)、气孔导度(Gs)和蒸腾
速率(Tr)增加, 同时使玉米叶片的暗呼吸速率(Rd)降
低, 但却并没有改变叶片的细胞间CO2浓度(Ci)和水
分利用效率(WUE)。实验增温使玉米叶片的净光合
速率(Pn)增加约 50% (P<0.05), 气孔导度(Gs)和蒸腾
速率(Tr)分别增加 163% (P<0.001)和 81% (P<0.05),
而使暗呼吸速率(Rd)降低约 25% (P<0.01)。此外, 实
验增温增加玉米叶片的细胞间 CO2 浓度(Ci) 25%
(P>0.05)和降低叶片水分利用效率 (WUE) 15%
(P>0.05), 但是该影响并没有达到显著水平。
第 4期 郑云普等: 玉米叶片气孔特征及气体交换过程对气候变暖的响应 607


表 2 增温对玉米叶片上不同位置气孔特征的影响
Table 2 Effects of warming on the stomatal traits among different sections of maize leaves
气孔参数
Stomatal parameter
气孔密度 SD
(Number mm–2)
气孔指数
SI (%)
气孔长度
SAL (mm)
气孔宽度
SAW (mm)
气孔面积
SAA (mm2)
气孔面积指数
SAAI (%)
上部 Tip 51 13.4 38.9 6.9 112 0.57
中部 Middle 61 14.0 49.8 6.8 117 0.71 近轴面
Adaxal
下部 Base 56 14.0 48.9 7.2 126 0.71
上部 Tip 73 19.2 38.7 5.8 84 0.61
中部 Middle 79 19.7 49.0 6.6 97 0.77
对照
Control
远轴面
Abaxial
下部 Base 79 19.4 46.0 8.2 120 0.95
上部 Tip 50 13.2 35.2 8.3 121 0.60
中部 Middle 58 17.3 34.8 9.6 157 0.91 近轴面
Adaxal
下部 Base 67 19.3 34.8 8.2 168 1.14
上部 Tip 81 19.7 36.3 7.6 120 0.98
中部 Middle 82 21.7 31.8 8.8 136 1.11
增温
Warming
远轴面
Abaxial
下部 Base 80 20.6 31.1 9.6 156 1.25
温度 Temperature <0.050 <0.001 <0.001 <0.001 <0.001 <0.001
叶面 Leaf surfaces <0.001 <0.001 0.087 0.971 <0.001 <0.001
叶片部位 Leaf sections <0.001 <0.001 <0.010 <0.010 <0.001 <0.001
温度×叶面 Temperature × Leaf surfaces 0.536 0.059 0.587 0.714 0.330 <0.050
温度×部位 Temperature×Leaf sections 0.252 <0.010 <0.001 0.159 0.064 <0.050
叶面×部位 Leaf surfaces×Leaf sections <0.050 <0.050 0.252 <0.010 0.364 0.363
温度×叶面×部位
Temperature×surfaces×sections
<0.050 0.068 0.630 0.661 0.086 <0.001
所有数据为平均值, 气孔密度和气孔指数均为 n = 25, 气孔长度、气孔宽度及气孔面积为 n = 150。所有数据的平均值利用多因
素方差分析在 0.05水平上进行比较, 对差异显著的组再进行多重比较(Duncan’s multiple range test)。缩写同表 1。
Values are means ± standard deviation for SD, SI (25 subsamples and 3 replicates) and for SAL, SAW, and SAA (150 subsamples and 3
replicates). Means are compared by the ANOVA followed by Duncan’s multiple range test. Values followed by different letters are signifi-
cantly different at P<0.05. Abbreviations are the same as those given in Table 1.

3 讨论
气孔密度和气孔指数是确定叶片气体交换有效
面积关键的生理生态参数, 直接决定着叶片潜在的
最大气孔导度[1]。许多研究发现增温显著改变了气
孔的密度, 但不影响气孔指数[13,15-17]。然而, 另一些
研究的结果显示增温并没有影响气孔密度, 却由于
表皮细胞密度的变化而改变了气孔指数 [14]。另外 ,
还有一些研究发现增温既不影响气孔密度也不改变
气孔指数[9-12]。本研究的结果表明, 增温主要是通过
减少表皮细胞的数目而增加玉米叶片的气孔指数 ,
并不是直接影响单位叶面积内的气孔数目。另外 ,
增温导致玉米叶片表皮细胞密度减少 9% (P<0.05),
而表皮细胞面积增加 35% (P<0.01)直接支持了上述
结论(表 3)。此外, 扫描电子显微镜观察的结果也显
示增温条件下玉米叶片表皮细胞与对照相比体积更
大、数目更少(图 2)。综上所述, 尽管全球变暖不会
显著影响气孔数目, 但可能会通过直接影响表皮细
胞的分化和扩展而进一步影响气孔的指数。
以往的研究发现, 拟南芥保卫细胞的长度由其
基因组的大小决定, 并且不受CO2浓度、干旱胁迫、
相对湿度、辐射强度、紫外线以及病原体侵染等环
境因素的影响[42]; 但不幸的是该研究并没有涉及到
温度对保卫细胞长度的影响。然而, 有趣的是本研
究的结果显示 , 增温显著减小了保卫细胞的长度 ,
从而导致气孔长度的减小。本研究结果也直接支持
了张立荣等[18]的研究结果, 该研究在对青藏高原上
亚高山草本植物进行增温实验时也同样发现了增温
显著降低该区域内4种草本植物叶片的气孔长度。然
而, 值得注意的是增温使气孔长度减小的同时也导
致了气孔宽度的增加, 以至于增温显著增大了气孔
的面积和气孔面积的指数(P<0.01; 表1)。由于叶片
气孔导度和蒸腾速率通常与气孔指数、气孔密度及
气孔开度大小呈正相关[2-3], 本研究的结果表明, 增
温增加了玉米叶片的气孔导度和蒸腾速率。同时 ,
608 作 物 学 报 第 41卷


对于叶片气体交换参数直接测量的结果也支持了上
述结论。因此, 未来气候变暖不仅可能增加玉米叶
片气孔指数和气孔开度大小, 而且在一定程度上导
致叶片气孔导度和蒸腾速率的增加。
许多的研究结果显示, 玉米叶片的背腹面(近轴
面/远轴面)极性在其分生组织阶段就已经建立起来,
并在随后整个叶片发育的阶段始终维持这种极性[23]。
本研究发现增温没有改变气孔密度的近轴面 /远轴
面比例, 表明玉米叶片的背腹面极性确实是受基因
控制的, 并不受外界环境温度的影响。另外, Driscoll
等[28]也报道了相似的研究结果, 他们发现玉米气孔
的背腹面极性不受 CO2 浓度的影响。然而, 增温却
显著增加了玉米叶片气孔指数的近轴面 /远轴面比
例, 这种结果可能主要是由于叶片背腹面表皮细胞
减少的比例不同而造成的。上述的研究结果表明 ,
玉米响应未来全球变暖可能会通过改变近轴面 /远
轴面表皮细胞比例来影响气孔在叶片背腹面的分
布。除了植物叶片不同的叶面外, 气孔在同一叶片
的不同部位也呈现出较高的变异性[14,16,29-33]。本研究
的结果也显示, 增温不仅影响玉米近轴面和远轴面
的气孔频度, 而且还改变玉米同一叶面上不同部位
(即尖部、中部及下部)的气孔频度, 表明未来气候变
暖不仅导致玉米叶面间的高变异, 还会引起玉米叶
面内不同部位的变异。相似的研究结果也同样在两
个常绿草本物种黑麦草(Lolium perenne)[14]和羊草
(Leymus chinensis)[16]中发现。然而, 以往多数的研究
在探讨增温对植物叶片气孔特征的影响取样时仅考
虑在叶片远轴面的中间部位取样[10-11,13,17], 这种取
样方法不够全面。因此, 今后在对植物叶片气孔研
究的取样方法上应该改进, 至少在取样时应该考虑
整个叶片多点采样。另外还应该注意的是, 尽管本
研究中的样品叶是同一位置的 , 但增温还可能会

图 4 玉米叶片近轴面气孔空间分布格局
Fig. 4 Spatial distribution pattern in the adaxial surface of maize leaves
对照环境下玉米叶片 1、叶片 2、叶片 3近轴面(a~c)气孔的空间分布格局; 增温环境下玉米叶片 1、叶片 2、叶片 3近轴面(d~f)气孔
的空间分布格局。上下两条包迹线代表 95%的置信区间。
Point pattern analyses of stomata on the adaxial surface (ac) in leaf 1, leaf 2, and leaf 3 of maize plants grown at ambient temperature and on
those (df) at elevated temperature, respectively. The dotted lines give a 95% confidence envelope for complete spatial randomness. The data
were given for three leaves from three ambient or warmed plots.
第 4期 郑云普等: 玉米叶片气孔特征及气体交换过程对气候变暖的响应 609



图 5 玉米叶片远轴面气孔空间格局分析
Fig. 5 Spatial distribution pattern in the abaxial surface of maize leaves
对照环境下玉米叶片 1、叶片 2、叶片 3远轴面(图 4-a~c)气孔的空间格局; 增温环境下玉米叶片 1、叶片 2、叶片 3远轴面(d~f)气孔
的空间格局。上、下 2条包迹线代表 95%的置信区间。
Point pattern analyses of stomata on the abaxial surface (ac) in leaf 1, leaf 2, and leaf 3 of maize plants grown at ambient temperature and on
those (df) at elevated temperature, respectively. The dotted lines give a 95% confidence envelope for complete spatial randomness. The data
were given for three leaves from three ambient or warmed plots.

显著提前生育进程。所以即使是同一个位置的叶片,
其生育日龄也可能存在显著的差异, 进而其叶的气
孔特征也可能发生显著变化。因此, 增温导致叶片
生育进程的不同, 也可能对叶片气孔特征产生潜在
的影响。
以往的研究结果对于叶片光合对增温的响应还
存在较大的争议[9,43]。本研究发现: (1)增温显著增加
了叶片的气孔指数(P<0.05; 表 1), 使玉米叶片上气
孔和表皮细胞的比例改变, 即气孔在叶片上所占的
比例增大; (2)增温还增加了玉米叶片气孔开口的面
积(P<0.01; 表 1), 即气孔的开张度更大; (3)增温也
同时增加了气孔面积指数(P<0.01; 表 1), 即单位叶
面积上气孔的总面积。这些参数的增加可能在一定
程度上有利于大气 CO2更加顺利地通过气孔进入叶
片的内部, 并扩散到相应的光合作用位点, 从而提
高玉米叶片的净光合速率。更为有趣的是增温导致
玉米叶片气孔的分布更加规则(图 4和图 5)。气孔在
玉米叶片上的分布越规则, 叶片的气体交换效率就
越高, 因为假如气孔在叶片上随机或簇状分布的话,
CO2分子经过气孔扩散到叶片内部的有效面积减小,
即气体交换的有效面积减小。该结果表明增温可以
改变叶片气孔的分布格局使更多 CO2分子更容易通
过气孔扩散到叶片内部, 从而有利于玉米叶片净光
合速率的增加。
4 结论
未来全球气候变暖并不会对我国华北平原玉米
叶片的气孔密度产生直接的影响, 但却可能通过改
变表皮细胞的数目来影响气孔指数。同时, 未来气
候变暖还会显著减小气孔开口的长度, 增加气孔开
口的宽度, 从而增加气孔开口的面积和周长。另外,
全球气候变暖还可能使气孔在玉米叶片上的分布格
610 作 物 学 报 第 41卷



图 6 实验增温对玉米叶片气体交换的影响
Fig. 6 Effects of experimental warming on the gas exchange of maize leaves
Pn: 净光合速率; Gs: 气孔导度; Tr: 蒸腾速率; Rd: 暗呼吸速率; Ci: 细胞间 CO2浓度; WUE: 水分利用效率。
Pn: net photosynthetic rates; Gs: stomatal conductance; Tr: transpiration rates; Rd: dark respiration rates; Ci: intercellular CO2 concentration;
WUE: water use efficiency.

表 3 增温对玉米叶片表皮细胞和保卫细胞特征的影响
Table 3 Effects of experimental warming on the characteristics of epidermal and guard cells of maize leaves
对照 Control 增温 Warming 参数
Parameter 近轴面 Adaxial 远轴面 Abaxial 近轴面 Adaxial 远轴面 Adaxial
增量
Increase (%)
P值
P-value
表皮细胞 Epidermal cells
349±27 a 321±11 ab 296±28 b 314±16 ab — — 表皮细胞密度 ECD
(Number mm–2) 335 305 –9.0 P<0.05
90.7±5.4 a 74.5±14.3 b 94.3±19.0 a 82.6±16.4 ab — — 表皮细胞长度
ECL (µm) 81.7 89.3 9.3 P>0.05
24.0±3.0 bc 21.9±4.4 c 29.1±4.4 a 25.6±3.0 b — — 表皮细胞宽度
ECW (µm) 22.9 27.6 20.6 P<0.01
2104±326 b 1496±309 c 2638±551 a 2051±392 b — — 表皮细胞面积
ECA (μm2) 1766 2384 35.0 P<0.01
247±24 a 183±26 b 293±82 a 249±54 a — — 表皮细胞周长
ECC (µm) 212 274 29.6 P<0.01
第 4期 郑云普等: 玉米叶片气孔特征及气体交换过程对气候变暖的响应 611


(续表 3)
对照 Control 增温 Warming 参数
Parameter 近轴面 Adaxial 远轴面 Abaxial 近轴面 Adaxial 远轴面 Adaxial
增量
Increase (%)
P值
P-value
表皮细胞密度比率
Adaxial/abaxial
1.09 0.94 –13.3 P<0.05
保卫细胞 Guard cells
52.4±4.3 a 49.8±4.0 b 47.1±3.3 c 42.2±2.6 d — — 保卫细胞长度
GCL (μm) # 51 45 –11.8 P<0.05
6.2±2.87 d 8.3±1.24 c 11.6±1.6 b 12.6±1.2 a — — 保卫细胞宽度
GCW (μm) # 7.3 12.1 66 P<0.01
453±78 b 468±77 ab 497±70 ab 504±72 a — — 保卫细胞面积
GCA (μm2) 461 501 8.0 P>0.05
126±11 a 129±9 a 130±15 a 134±22 a — — 保卫细胞周长
GCC (μm) 128 132 3.1 P>0.05
所有数据为平均值±标准偏差, 表皮细胞密度为 n = 75, 其他所有参数均为 n = 450。所有数据的平均值利用单因素方差分析在
0.05水平上进行比较, 对差异显著的组再进行多重比较(Duncan’s multiple range test)。标以不同字母的值。#保卫细胞长度定义为保卫
细胞长度方向上最长的距离, 保卫细胞宽度定义为保卫细胞宽度方向上最宽的距离。ECD: 表皮细胞密度; ECL: 表皮细胞长度; ECW:
表皮细胞宽度; ECA: 表皮细胞面积; ECC: 表皮细胞周长; GCL: 保卫细胞长度; GCW: 保卫细胞宽度; GCA: 保卫细胞面积; GCC: 保
卫细胞周长。
Values are means ± standard deviation for ECD (75 subsamples, 3 replicates) and for other parameters (450 subsamples, 3 replicates).
Means were compared by the ANOVA followed by Duncan’s multiple range tests. Values followed by the different letters are significantly
different at P<0.05. # Guard cell length is the longest dimension, and the guard cell width is the widest dimension. ECD: epidermal cell den-
sity; ECL: epidermal cell length; ECW: epidermal cell width; ECA: epidermal cell area; ECC: epidermal cell circumference; GCL: guard cell
length; GCW: guard cell width; GCA: guard cell area; GCC: guard cell circumference.

局变得更加均匀。这些特征参数的改变导致华北平
原玉米叶片气体交换过程显著变化, 从而可能最终
影响到玉米的品质和产量。
References
[1] Woodward F I. Stomatal numbers are sensitive to increases in
CO2 from preindustrial levels. Nature, 1987, 327: 617–618
[2] Hetherington A M, Woodward F I. The role of stomata in sensing
and driving environmental change. Nature, 2003, 424: 901–908
[3] Franks P J, Beerling D J. Maximum leaf conductance driven by
CO2 effects on stomatal size and density over geologic time. Proc
Natl Acad Sci USA, 2009, 106: 10343–10347
[4] Haworth M, Heath J, McElwain J C. Differences in the response
sensitivity of stomatal index to atmospheric CO2 among four
genera of Cupressaceae conifers. Ann Bot, 2010, 105: 411–418
[5] Taylor S H, Franks P J, Hulme S P. Photosynthetic pathway and
ecological adaptation explain stomatal trait diversity amongst
grasses. New Phytol, 2012, 193: 387–396
[6] Lake J A, Woodward F I. Response of stomatal numbers to CO2
and humidity: control by transpiration rate and abscisic acid. New
Phytol, 2008, 179: 397–404
[7] Ciais P, Denning A S, Tans P P. A three-dimensional synthesis of
vegetation feedbacks in doubled CO2 climate experiments. J
Geophys Res, 1997, 102: 5857–5872
[8] Buckley T N, Farquhar G D, Mott K A. Qualitative effects of
patchy stomatal conductance distribution features on gas-
exchange calculations. Plant Cell Environ, 1997, 20: 867–880
[9] Apple M E, Olszyk D M, Ormrod D P. Morphology and stomatal
function of douglas fir needles exposed to climate change: ele-
vated CO2 and temperature. Int J Plant Sci, 2000, 161: 127–132
[10] Hovenden M J. The influence of temperature and genotype on the
growth and stomatal morphology of southern beech, Nothofagus
cunninghamii (Nothofagaceae). Aust J Bot, 2001, 49: 427–434
[11] Kouwenberg L L R, Kürschner W M, McElwain J C. Stomatal
frequency change over altitudinal gradients: Prospects for pa-
leoaltimetry. Rev Mineral Geochem, 2007, 66: 215–241
[12] Fraser L H, Greenall A, Carlyle C. Adaptive phenotypic plasticity
of Pseudoroegneria spicata: response of stomatal density, leaf
area and biomass to changes in water supply and increased tem-
perature. Ann Bot, 2009, 103: 769–775
[13] Beerling D J, Chaloner W G. The impact of atmospheric CO2 and
temperature change on stomatal density: observations from
Quercus robur Lammas leaves. Ann Bot, 1993, 71: 231–235
[14] Ferris R, Nijs I, Behaeghe T. Elevated CO2 and temperature have
different effects on leaf anatomy of perennial ryegrass in spring
and summer. Ann Bot, 1996, 78: 489–497
[15] Reddy K R, Robana R R, Hodges H F. Interactions of CO2 en-
richment and temperature on cotton growth and leaf characteris-
tics. Environ Exp Bot, 1998, 39: 117–129
[16] Xu Z Z, Zhou G S, Shimizu H. Effects of soil drought with noc-
turnal warming on leaf stomatal traits and mesophyll cell ultra-
structure of a perennial grass. Crop Sci, 2009, 49: 1843–1851
[17] Xu Z Z, Zhou G S. Effects of water stress and high nocturnal
temperature on photosynthesis and nitrogen level of a perennial
grass Leymus chinensis. Plant Soil, 2005, 269: 131–139
[18] 张立荣, 牛海山, 汪诗平, 李英年, 赵新全. 增温与放牧对矮
嵩草草甸 4 种植物气孔密度和气孔长度的影响. 生态学报,
612 作 物 学 报 第 41卷


2010, 30: 6961–6969
Zhang L R, Niu H S, Wang S P, Li Y N, Zhao X Q. Effects of tem-
perature increase and grazing on stomatal density and length of four
alpine Kobresia meadow species, Qinghai-Tibetan Plateau. Acta
Ecol Sin, 2010, 30: 6961–6969 (in Chinese with English abstract)
[19] Croxdale J L. Stomatal patterning in monocotyledons: Trades-
cantia as a model system. J Exp Bot, 1998, 49: 279–292
[20] Berger D, Altmann T. A subtilisin-like serine protease involved in
the regulation of stomatal density and distribution in Arabidopsis
thaliana. Genes Dev, 2000, 14: 1119–1131
[21] Croxdale J L. Stomatal patterning in angiosperms. Am J Bot,
2000, 87: 1069–1080
[22] Bergermann D C, Lukowitz W, Somerville C R. Stomatal deve-
lopment and pattern controlled by a MAPKK Kinase. Science,
2004, 304: 1494–1497
[23] Juarez M, Twigg R, Timmermans M. Specification of adaxial cell
fate during maize leaf development. Development, 2004, 131:
4533-4544
[24] Wang H, Ngwenyama N, Liu Y. Stomatal development and pat-
terning are regulated by environmentally responsive mito-
gen-actived protein kinases in Arabidopsis. Plant Cell, 2007, 19:
63–73
[25] Casson S A, Hetherington A M. Environmental regulation of
stomatal development. Curr Opin Plant Biol, 2010, 13: 90–95
[26] Ciha A J, Brun W A. Stomatal size and frequency in soybeans.
Crop Sci, 1975, 15: 309–313
[27] Green R L, Beard J B, Casnoff D M. Leaf blade stomatal charac-
terizations and evapotranspiration rates of 12 cool-season peren-
nial grasses. HortScience, 1990, 25: 760–761
[28] Driscoll S P, Prins A, Olmos E. Specification of adaxial and
abaxial stomata, epidermal structure and photosynthesis to CO2
enrichment in maize leaves. J Exp Bot, 2006, 57: 381–390
[29] Salisbury E J. On the causes and ecological significance of
stomatal frequency, with special reference to the woodland flora.
Philos Trans R Soc Lond B Biol Sci, 1927, 216: 1–65
[30] Sharma G K, Dunn D B. Environmental modifications of leaf surface
traits in Datura stramonium. Can J Bot, 1969, 47: 1211–1216
[31] Tichá I. Photosynthetic characteristics during ontogenesis of
leaves: 7. Stomata density and sizes. Photosynthetica, 1982, 16:
375–471
[32] Zacchini M, Morini S, Vitagliano C. Effect of photoperiod on
some stomatal characteristics of in vitro cultured fruit tree shoots.
Plant Cell Tissue Organ Cult, 1997, 49: 195–200
[33] Stancato G C, Mazzoni-Viveiros S C, Luchi A E. Stomatal char-
acteristics in different habitat forms of Brazilian species of
Epidendrum (Orchidaceae). Nord J Bot, 1999, 19: 271–275
[34] Tao F, Yokozawa M, Xu Y, Hayashi Y, Zhang Z. Climate changes
and trends in phenology and yields of field crops in China,
1981–2000. Agric For Meteorol, 2006, 138: 82–92
[35] Lin E. Agricultural vulnerability and adaptation to global warm-
ing in China. Water Air Soil Pollut, 1996, 92: 63–73
[36] Mo X G, Liu S X, Lin Z H, Guo R P. Regional crop yield, water
consumption and water use efficiency and their responses to cli-
mate change in the North China Plain. Agric Ecosyst Environ,
2009, 134: 67–78
[37] Liu S X, Mo X G, Lin Z H, Xu Y Q, Ji J J, Wen G, Richey J. Crop
yield responses to climate change in the Huang-Huai-Hai Plain of
China. Agric Water Manag, 2010, 97: 1195–1209
[38] Tao F L, Zhang S, Zhang Z. Spatiotemporal changes of wheat
phenology in China under the effects of temperature, day length
and cultivar thermal characteristics. Eur J Agron, 2012, 43:
201–212
[39] Hou R X, Ou-Yang Z, Li Y S. Is the change of winter wheat yield
under warming caused by shortened reproductive period? Ecol
Evol, 2012, 2: 2999–3008
[40] Lobell D B, Field C B. Global scale climate–crop yield relation-
ships and the impacts of recent warming. Environ Res Lett, 2007,
2: 1–7
[41] Ripley B D. The second-order analysis of stationary processes. J
Appl Prob, 1976, 13: 255–266
[42] Lomax B H, Woodward F I, Leitch I J. Genome size as a predic-
tor of guard cell length in Arabidopsis thaliana is independent of
environmental conditions. New Phytol, 2009, 181: 311–314
[43] Djanaguiraman M, Prasad P V V, Boyle D L. High-tem-
perature stress and soybean leaves: leaf anatomy and photo-
synthesis. Crop Sci, 2011, 51: 2125–2131