Dynamics and simulation of leaf area index for Artemisia ordosica community in the Mu Us Desert of northwestern China
-
摘要:目的本研究以毛乌素沙地灌木油蒿群落为研究对象,研究如何快速、准确、连续地获取油蒿群落LAI,获取长期连续的LAI数据对研究生态系统过程与环境变化间的相互关系具有重要意义。方法本文采用原位连续监测和野外试验的方法,于2013和2014年每年4—10月使用LAI-2000冠层分析仪定期测定油蒿群落LAI,同时利用太阳辐射传感器和光量子传感器测得的辐射数据计算地面观测NDVI(NDVIground)并获取同期遥感MODIS NDVI(NDVIMODIS)。用LAI数据分析油蒿叶面积群落指数的季节动态,并与归一化积温、NDVIground和NDVIMODIS分别建立经验模型,得到最优油蒿群落LAI模拟模型。结果油蒿群落LAI在整个生长季内的变化与物候具有一致性,整体上随时间呈现单峰型变化趋势,4—8月末油蒿持续生长,LAI稳定增大,8月份达到峰值,2013年为1.09m2/m2,2014年为1.33m2/m2。9月初至9月中旬油蒿冠层结构趋于稳定,LAI变化较小;9月下旬油蒿叶片开始枯黄掉落,LAI迅速下降。经验证,利用NDVIground模拟荒漠地区油蒿群落LAI效果最优(R2=0.76, P < 0.01)。结论在荒漠地区简单测量NDVI并使用本文研究模型,可以快速、无破坏地获取长期连续可靠的LAI数据,研究结果对理解和预测荒漠生态系统对全球气候变化的响应具有重要意义。Abstract:ObjectiveOur objective here was to rapidly and accurately obtain long-term continuous LAI for Artemisia ordosica communities in the Mu Us Desert of northwestern China, which is of great importance for understanding the relationship between ecosystem processes and environmental changes.MethodWe measured LAI in a typical A. ordosica community using the LAI-2000 canopy analyzer (LI-COR, USA) from April to October in the growing season of 2013 and 2014, and obtained MODIS NDVI (NDVIMODIS) data for the same periods. Ground-based NDVI (NDVIground) was calculated from incident and reflected solar radiation and photosynthetically active radiation, which were measured by radiation sensors mounted on a tower in the center of the community. LAI measurements were used to examine seasonal changes and construct LAI models. Normalized effective accumulated temperature, NDVIground and NDVIMODIS were used to construct optimal LAI model for the community.ResultThe seasonal dynamics was consistent with the phenophases of A. ordosica and showed a hump-shaped pattern. LAI increased from April to late August, reached a peak in August (1.09 and 1.33m2/m2 in 2013 and 2014, respectively), remained relatively stable from early September to mid-September, and rapidly declined from late September due to the defoliation of A. ordosica. The optimal LAI model was NDVIground (R2=0.76, P < 0.01).ConclusionOur results indicate that long-term continuous LAI estimates for A. ordosica communities in the Mu Us Desert can be easily obtained using NDVI measurements and the model. The model can be used to understand and predict ecological impacts of climate change.
-
Keywords:
- Artemisia ordosica /
- leaf area index /
- accumulated temperature /
- NDVI
-
随着天然林的过度采伐,优质木材资源短缺,开发人工林速生材成为了木材工业的重要途径。但人工林速生材的密度低、力学性能差,进一步限制了其有效利用。木材化学改性是通过将化学物质引入到木材内部,占据木材细胞壁中的孔隙并与木材组分中的亲水性羟基反应,降低其亲水性的同时增强木材的物理力学性能,从而延长其使用寿命,扩大其应用范围[1-3]。
纳米技术的发展为木材化学改性提供了新的思路。木材是多孔性材料,具备宏观孔隙﹑介观孔隙和微观孔隙在内的多级孔隙结构[4]。纳米孔隙的存在为纳米颗粒的进入创造了条件。以蒙脱土(montmorillonite,MMT)为代表的纳米黏土被广泛应用于木材的化学改性中,在添加量较少(3% ~ 5%)的情况下就可以大幅度提高木材的力学、防水、阻燃等性能[5-7]。通常天然MMT亲水且易团聚,因此对木材改性效果不佳,需要采用有机改性剂改性为有机蒙脱土(organo-montmorillonite,OMMT)。而OMMT难以均匀分散到水中,通常需要用水溶性树脂作为中间介质,先将OMMT分散进入树脂中再浸渍处理木材[8]。但由于承载OMMT的树脂乳液粒径较大,黏度较高,对木材的渗透效果较差,因而固化后蒙脱土多填充于木材细胞腔内,仅部分纳米片层进入到了细胞壁的无定形区[9]。虽然这种传统MMT改性方法可以提高木材的阻燃性能,但对木材物理力学性能的提高有限。
超支化聚合物是一种具有高度分支结构的树枝状聚合物。将水性超支化聚合物接枝到纳米OMMT上,由于其球形以及多支化结构,可以使OMMT在水中稳定分散[10-12]。Li等[13]成功制备了一种水性超支化聚丙烯酸酯(hyperbranched polyacrylate,HBPA)分散蒙脱土乳液,其平均粒径为100 nm。木材经浸渍处理后,在细胞壁中发现了剥离的OMMT片层以及固化的HBPA,因而处理材的防水、硬度等性能有所提升,但OMMT是首先分散到丙烯酸酯单体中再聚合制备成乳液,因此受乳液pH、黏度和离子等因素的影响,分散进入HBPA乳液中的OMMT含量较少,只能小幅度提高木材的物理力学性能。为进一步提高OMMT的含量,本研究首先合成了HBPA乳液,再将OMMT直接添加到HBPA乳液中。另外,Xu等[14]通过侧链丙烯酸控制还可以合成一种能与聚乙二醇(poly(ethyl glycol), PEG)羟基发生酯化或氢键作用形成物理网络结构的水性树脂,使复合的PEG/HBPA能与木材羟基进一步发生氢键作用,与OMMT产生协同效应。因此,本研究尝试在HBPA中加入低分子量聚乙二醇200(PEG-200),形成物理交联的网络结构从而进一步提高木材的性能,并比较了不同层间离子的OMMT对复合乳液改性材力学性能的影响。
1. 材料与方法
1.1 材 料
青杨(Populus cathayana)取自河南省漯河市,选择无明显节子、腐朽等自然缺陷的边材作为试材,根据相应测试标准锯切成规定尺寸,自然气干备用。4种不同层间离子的有机蒙脱土编号分别为OMMT-1、OMMT-2、OMMT-3、OMMT-4,其有机长链分别为[CH3(CH2)17]2(CH3)2N+、HCH3(CH2)16CH2NH2+、CH3(CH2)17N(CH3)[(CH2CH2OH)2]+、HOOC(CH2)17NH3+,均过200目筛,购买自北京怡蔚特化科技发展有限公司。水性超支化聚丙烯酸乳液,平均粒径100 nm,黏度120 mPa·s,固体含量10%,实验室自制,其化学结构式如图1所示。PEG-200,购买自国药集团化学试剂有限公司。
1.2 蒙脱土/聚乙二醇/超支化聚丙烯酸酯乳液的制备
将PEG-200加入到HBPA溶液中混合均匀,PEG的添加质量为HBPA固含质量的1/3。用去离子水将混合液稀释至4%后加入质量分数2%的OMMT,之后在1 000 r/min的速度下搅拌10 min,即得到OMMT/PEG/HBPA乳液。同时,制备一组未添加OMMT的PEG/HBPA乳液作为对照组。
1.3 木材浸渍处理
处理前,先将青杨试件在103 ℃下干燥至恒重,之后将试件置于浸渍罐中用处理液进行真空–加压浸渍处理。真空度为0.01 MPa,时间30 min;加压压力为2 MPa,加压时间120 min。浸渍完成后,将试件从浸渍罐中取出,擦去表面多余水分,先在40 ℃恒温干燥箱中干燥24 h,随后在103 ℃下干燥至恒定。
1.4 表征与测试
将复配改性剂乳液在室温下放置24 h,观测改性剂乳液是否出现了分层、沉淀等现象,并采用激光粒度仪测试乳液的粒径,旋转黏度计测试乳液黏度,考察改性剂乳液的稳定性。
参照GB/T 1935—2009《木材顺纹抗压强度试验方法》、GB/T 1936.1—2009《木材抗弯强度试验方法》、GB/T 1941—2009《木材硬度试验方法》分别测试处理材的顺纹抗压强度、抗弯强度和端面硬度,每组重复试件数为10个。
采用美国尼高力公司生产的IS10型傅里叶红外交换光谱仪对改性剂乳液以及处理材进行傅里叶红外交换光谱(FTIR)测试。将处理材粉碎过100目筛子,按绝干质量1∶100与溴化钾混合均匀后压片,改性剂乳液则稀释100倍后取一滴与溴化钾压片混合后烘干。测试波长范围为400 ~ 4 000 cm−1,分辨率为4 cm−1,扫描次数为32次。
采用德国布鲁克公司生产的D8 Advance型X射线衍射仪对OMMT以及处理材进行X射线衍射(XRD)测试。将处理材粉碎过100目筛子,采用Cu-Kα进行测试,扫描范围为1.5° ~ 40°,扫描速度为每步0.5 s。
采用美国FEI公司生产的Tecnai G2 F30型场发射透射电子显微镜对处理材进行透射电子显微镜(TEM)观测。处理材首先用环氧树脂包埋,然后用超薄切片刀切片后置于铜网上观测,并采用X射线能谱仪(EDX)点扫描方式对选定区域进行Si和Al元素确定,测试电压为10 kV。
2. 结果与分析
2.1 改性剂稳定性
图2为OMMT/PEG/HBPA改性剂放置24 h后的外观状态。改性剂乳液呈浅黄色,无明显的分层和沉淀出现。从表1中也可以发现改性剂乳液的平均粒径和乳液黏度基本无变化,说明了这4种OMMT均能稳定分散进入到PEG/HBPA乳液中,可用于木材浸渍改性处理。
表 1 OMMT/PEG/HBPA改性剂乳液静置24 h后的平均粒径与黏度变化Table 1. Average particle size and viscosity changes of OMMT/PEG/HBPA modifier emulsions after standing for 24 hours编号 No. 平均粒径 Average particle size/nm 黏度 Viscosity/(mPa·s) 静置前 Before standing 静置后 After standing 静置前 Before standing 静置后 After standing PEG/HBPA 108 110 118 113 OMMT-1/PEG/HBPA 106 107 120 122 OMMT-2/PEG/HBPA 109 112 123 121 OMMT-3/PEG/HBPA 110 106 130 127 OMMT-4/PEG/HBPA 108 117 124 129 注:PEG/HBPA是聚乙二醇/超支化聚丙烯酸酯乳液;OMMT-1/PEG/HBPA是层间离子[CH3(CH2)17]2(CH3)2N+的有机蒙脱土/聚乙二醇/超支化聚丙烯酸酯乳液;OMMT-2/PEG/HBPA是层间离子为的HCH3(CH2)16CH2NH2+的有机蒙脱土/聚乙二醇/超支化聚丙烯酸酯乳液;OMMT-3/PEG/HBPA是层间离子为CH3(CH2)17N(CH3)[(CH2CH2OH)2]+的有机蒙脱土/聚乙二醇/超支化聚丙烯酸酯乳液;OMMT-4/PEG/HBPA是层积钠离子为HOOC(CH2)17NH3+的有机蒙脱土/聚乙二醇/超支化聚丙烯酸酯乳液。Notes: PEG/HBPA is poly(ethyl glycol)/hyperbranched polyacrylate emulsion. OMMT-1/PEG/HBPA is OMMT with interlayer ions of [CH3(CH2)17]2(CH3)2N+/poly(ethyl glycol)/hyperbranched polyacrylate emulsion. OMMT-2/PEG/HBPA is OMMT with interlayer ions of HCH3(CH2)16CH2NH2+/poly(ethyl glycol)/hyperbranched polyacrylate emulsion. OMMT-3/PEG/HBPA is OMMT with interlayer ions of CH3(CH2)17N(CH3)[(CH2CH2OH)2]+/poly(ethyl glycol)/hyperbranched polyacrylate emulsion. OMMT-4/PEG/HBPA is OMMT with interlayer ions of HOOC(CH2)17NH3+/poly(ethyl glycol)/hyperbranched polyacrylate emulsion. 2.2 力学性能
OMMT/PEG/HBPA改性材的力学性能如表2所示。未改性木材的顺纹抗压强度、抗弯强度和端面硬度分别为45.6 MPa、62.3 MPa和5 010 N。经过PEG/HBPA改性后,木材的顺纹抗压强度和抗弯强度均有所提升。Xu等[14]发现:PEG能和HBPA侧链中的羧基发生酯交换反应或形成氢键结合,同时这两者较长的分子链互相缠绕,可形成网络的结构。而当PEG/HBPA浸渍进入木材后,由于PEG和HBPA干燥过程形成的交联以及与木材的羟基作用,使木材抵抗外力的能力增加,因而力学强度明显上升。而PEG/HBPA对端面硬度的提升不明显,这也是由于PEG和HBPA均为柔性的高分子树脂,不能有效提升木材的表面硬度。添加OMMT后,各组的顺纹抗压强度、抗弯强度和端面硬度数值进一步增加,其中OMMT的添加对提升木材端面硬度值效果明显。这是由于OMMT自身是一种强度较高的刚性纳米填料,不论是凝聚在木材外表面还是进入到细胞壁内,都能够有效提高其力学强度[15]。而当PEG/HBPA进入到OMMT层间后,其分子的网络结构能够使OMMT片层剥离,从而可能进入到木材细胞壁内[13]。由表2的结果也可以看出:对于不同层间离子的OMMT,其顺纹抗压强度、抗弯强度和端面硬度这3项力学性能大小排序皆为OMMT-2 > OMMT-3 > OMMT-4 > OMMT-1。这说明不同层间离子的OMMT对木材的改性效果有所不同,具体原因将在FTIR和TEM结果中进行分析。本试验中优化配方为OMMT-2添加组,其顺纹抗压强度达到82.2 MPa,抗弯强度98.2 MPa,端面硬度8 920 N。
表 2 OMMT/PEG/HBPA改性材的力学性能Table 2. Mechanical properties of wood modified with OMMT/PEG/HBPA编号 No. 顺纹抗压强度 Radial compressive strength/MPa 抗弯强度 Flexural strength/MPa 端面硬度 End hardness/N 未改性 Unmodified 45.6 (6.1) 62.3 (8.2) 5 010 (110) PEG/HBPA 59.1 (6.3) 77.4 (5.5) 5 220 (310) OMMT-1/PEG/HBPA 70.2 (8.3) 85.6 (6.3) 7 740 (430) OMMT-2/PEG/HBPA 82.2 (6.6) 98.2 (7.3) 8 920 (320) OMMT-3/PEG/HBPA 75.3 (2.8) 92.1 (6.3) 8 660 (440) OMMT-4/PEG/HBPA 73.3 (2.5) 88.1 (2.5) 8 110 (280) 注:括号内数值为10个重复试件的标准偏差。Note: values in the parentheses are standard deviations of 10 replicates. 2.3 FTIR分析
OMMT/PEG/HBPA改性材的FTIR结果如图3a所示。未处理材在3 333 cm−1处的羟基吸收峰经PEG/HBPA改性后向高峰3 378 cm−1处偏移,并且峰值增加。这是由于PEG/HBPA中存在大量羟基,能和木材中的羟基形成氢键结合[16]。未添加OMMT和OMMT-1组的羟基吸收峰强度较高,而OMMT-2、OMMT-3和OMMT-4组的羟基吸收峰强度相对低一些,说明加入OMMT-1后,OMMT-1与木材基本上为物理吸附,而其余几组的OMMT可能与木材之间形成了氢键,使羟基吸收峰下降。另外,各组在1 736 cm−1处的C=O吸收峰也有所增强,这是由于木材半纤维素中的C=O与HBPA中的C=O发生了重叠。经PEG/HBPA改性后,在1 050 cm−1处出现了新的吸收峰,此处是PEG/HBPA中C—O—C的伸缩振动。对PEG/HBPA乳液进行FTIR测试(图3b),乳液的C—O—C伸缩振动出现在1 110 cm−1处,进入木材后,此峰发生了明显的偏移,说明了改性剂与木材纤维素形成了氢键结合。相比未添加OMMT组,添加OMMT后,木材在2 918 cm−1处出现了较尖锐的吸收峰。此处归属于—CH3/—CH2的伸缩振动,说明了OMMT中有机长链进入到了木材中。此外,OMMT-2/PEG/HBPA在1 667 cm−1处的峰有所增强,同时1 736和1 586 cm−1处的峰均有所下降,这可能是由于OMMT-2与其余3组不同,为非季铵盐改性OMMT,其特有的官能团—NH2能够与木材以及改性剂中的羟基、羧基等反应[17]。
2.4 XRD分析
OMMT/PEG/HBPA改性材的XRD结果如图4所示。2θ在2° ~ 10°范围内,OMMT有一个明显的(001)晶面衍射峰,根据布拉格公式(1),可以计算得到OMMT的层间距。
d=λ2sinθ (1) 式中:d为层间距(nm);
λ 为X射线的波长,λ = 0.154 nm;θ为衍射角(°)。由图4可以看出,OMMT-1的层间距最小,为1.24 nm,而OMMT-2的层间距最大,为4.48 nm。因此,PEG/HBPA分子更容易进入到OMMT-2的层间,形成剥离结构。在浸渍处理木材后,在2θ = 15°和2θ = 22.5°的纤维素衍射峰没有发生明显变化,说明本研究所使用的改性方法没有破坏纤维素的结晶结构。同时,改性材的XRD图在2θ = 2° ~ 10°中未发现明显的OMMT晶面衍射峰,说明各组OMMT进入木材后均可能呈剥离状态。
2.5 TEM-EDX分析
通过TEM可以进一步确定OMMT在木材细胞壁中的状态,如图5所示。从图5a中看出未处理木材细胞壁清晰可见。经过PEG/HBPA处理后,PEG/HBPA固化后形成一层薄膜黏附在细胞壁上,部分进入到了细胞壁内(图5b)。加入OMMT后,可发现有片层状物质进入到了木材的细胞壁内,为确定其为OMMT,对片层区域(如框所示)进行了EDX分析,探测Si和Al两种特征元素,图6为OMMT-1/PEG/HBPA所选区域的EDX结果,结果表明所选区域出现了Si和Al元素,因此这种片层状物质为OMMT。其余OMMT-2、OMMT-3和OMMT-4组的EDX结果基本和图6一致,也出现了Si和Al元素。对比图5c、图5d、图5f可以发现OMMT-2组对木材细胞壁的渗透性最好,在细胞壁内发现了较多的层状片层结构(图5d),可能是由于自身OMMT的层间距较大,使用PEG/HBPA作为分散剂时能更加容易地进入到OMMT层间,同时,OMMT-2中存在的氨基官能团使之与木材、改性剂的结合更强。OMMT-1组对木材的渗透性最差,图5c中发现有OMMT片层黏附在木材细胞壁上。OMMT-3和OMMT-4组有少量的OMMT存在于细胞壁中(图5e和图5f),这可能是由于OMMT-3和OMMT-4的自身的层间距相对OMMT-1较大,且有羟基、羧基等官能团,能与木材和改性剂之间形成氢键。
3. 结 论
本研究采用PEG/HBPA分散OMMT,用于改性木材,并比较了不同层间离子OMMT对改性效果的影响。4种OMMT均能稳定分散进入到PEG/HBPA乳液中,经过24 h后无明显的分层和沉淀现象,且粒径和黏度无明显变化。处理材经过PEG/HBPA处理后,PEG/HBPA能固化并黏附在细胞壁表面或进入细胞壁内,使木材的抗压强度和抗弯强度有所提高,但对端面硬度增强作用不大。添加OMMT能进一步提高处理材的力学性能,并增加端面硬度。OMMT层间离子中含有氨基、羟基、羧基等官能团能使OMMT更好地进入到木材细胞壁中,其中OMMT-2处理组的增强效果较优,这可能是由于OMMT-2自身层间距较大,同时含有的氨基官能团与改性剂、木材形成较强的氢键结合,为本试验的优化处理组。经过处理后,改性材顺纹抗压强度为82.2 MPa,抗弯强度为98.2 MPa,端面硬度为8 920 N。
-
图 2 毛乌素沙地油蒿叶面积指数及归一化植被指数的季节动态
NDVIground.地面观测NDVI Ground-based NDVI; NDVIMODIS.遥感MODIS NDVI Remote sensing MODIS NDVI;LE.展叶期Leaf expansion stage; FB.开花结果期Flowering and bearing stage; LC.变色期Leaf coloration stage; DP.凋落期Defoliation stage
Figure 2. Seasonal variations of LAI and NDVI in Artemisia ordosica in Mu Us Desert
图 4 2014年不同方法油蒿群落叶面积指数模拟值与实测值的比较和相对误差
a、b为积温模拟法;c、d为地面观测NDVI模拟法;e、f为遥感反演法。
Figure 4. Comparison and relative error between simulated and measured LAI of Artemisia ordosica community in 2014
a and b, normalized effective accumulated temperature; c and d, ground-based NDVI; e and f, remote sensing retrieval.
表 1 2013和2014年油蒿群落物候期
Table 1 Phenophases of Artemisia ordosica community in 2013 and 2014
年份
Year展叶期
Leaf expansion period开花结果期
Flowering and bearing period变色期
Leaf coloration period凋落期
Defoliation period2013 94~145 170~248 235~287 242~287 2014 102~147 171~247 232~291 238~291 数据为一年中的第几天。Note: data is the ordinal date in a year. 表 2 不同方法模拟油蒿LAI的结果
Table 2 Simulated results of Artemisia ordosica LAI by different methods
模拟方法Simulation method 拟合方程Model expression R2 RMSE AIC 积温模拟法Simulation by NEAT LAI=0.63NEAT+0.47 0.83 0.09 -103.90 LAI=0.54exp(0.72NEAT) 0.77 0.10 -97.56 LAI=0.27ln(NEAT)+1.04 0.93 0.07 -123.77 LAI=1.07NEAT0.36 0.91 0.07 -118.80 LAI=1.01/(1+exp(0.97-6.03NEAT)) 0.95 0.05 -132.36 LAI=1.00/(1+exp(0.82-4.76NEAT-2.82NEAT2)) 0.95 0.05 -132.65 地面观测NDVI模拟法
Simulation by NDVIgroundLAI=4.32NDVIground-0.92 0.74 0.11 -90.07 LAI=0.10exp(5.10NDVIground) 0.69 0.12 -86.02 LAI=1.69ln(NDVIground)+2.37 0.76 0.11 -91.37 LAI=5.38NDVIground2.08 0.72 0.12 -88.03 遥感反演法
Simulation by remote sensing retrievalLAI=5.87NDVIMODIS-0.34 0.82 0.09 -103.05 LAI=0.20exp(7.07NDVIMODIS) 0.77 0.11 -97.35 LAI=1.10ln(NDVIMODIS)+2.62 0.84 0.09 -106.04 LAI=7.78NDVIMODIS1.39 0.81 0.10 -101.45 注:所有模型均通过0.01水平的显著性检验。Note: all models were examined by a significance test of 0.01 level. 表 3 不同植被基于NDVI的LAI估算模型
Table 3 LAI estimation models based on NDVI for different vegetation species
植被
Vegetation species模型方程
Model expressionR2 小麦Triticum aestivum LAI=exp(2.77-0.52/NDVI) 0.96 玉米Zea mays LAI=3.20NDVI+0.53 0.65 水稻Oryza sativa LAI=1.06NDVI+1.05 0.99 杨树Populus L. LAI=0.48exp(2.89NDVI) 0.51 草地Grassland LAI=1.42NDVI-0.09 0.62 LAI=0.13exp(NDVI/0.31) 0.59 -
[1] Chen J M, Black T A. Defining leaf area index for non-flat leaves[J]. Plant, Cell & Environment, 1992, 15(4):421-429. doi: 10.1111-j.1365-3040.1992.tb00992.x/
[2] 王龑, 田庆久, 王琦, 等.杨树林全生长期LAI遥感估算模型适用性[J].生态学报, 2016, 36(8):2210-2216. http://d.old.wanfangdata.com.cn/Periodical/stxb201608012 Wang Y, Tian Q J, Wang Q, et al. Study on the applicability of leaf area index estimation model in an aspen forest over a growth period[J]. Acta Ecologica Sinica, 2016, 36(8):2210-2216. http://d.old.wanfangdata.com.cn/Periodical/stxb201608012
[3] Bonan G B. Land-atmosphere interactions for climate system models:coupling biophysical, biogeochemical, and ecosystem dynamical processes[J]. Remote Sensing of Environment, 1995, 51(1):57-73. doi: 10.1016/0034-4257(94)00065-U
[4] Veroustraete F, Patyn J, Myneni R B. Estimating net ecosystem exchange of carbon using the normalized difference vegetation index and an ecosystem model[J]. Remote Sensing of Environment, 1996, 58(1):115-130. doi: 10.1016/0034-4257(95)00258-8
[5] 赵传燕, 沈卫华, 彭焕华.祁连山区青海云杉林冠层叶面积指数的反演方法[J].植物生态学报, 2009, 33(5):860-869. doi: 10.3773/j.issn.1005-264x.2009.05.004 Zhao C Y, Shen W H, Peng H H. Methods for determining canopy leaf area index of Picea crassifolia forest in Qilian Mountains, China[J]. Chinese Journal of Plant Ecology, 2009, 33(5):860-869. doi: 10.3773/j.issn.1005-264x.2009.05.004
[6] 曾小平, 赵平, 饶兴权, 等.鹤山丘陵3种人工林叶面积指数的测定及其季节变化[J].北京林业大学学报, 2008, 30(5):33-38. http://j.bjfu.edu.cn/article/id/8743 Zeng X P, Zhao P, Rao X Q, et al. Measurement of leaf area index of three plantations and their seasonal changes in Heshan hilly land[J]. Journal of Beijing Forestry University, 2008, 30(5):33-38. http://j.bjfu.edu.cn/article/id/8743
[7] Xiao Z Q, Liang S L, Jiang B. Evaluation of four long time-series global leaf area index products[J]. Agricultural and Forest Meteorology, 2017, 246:218-230. doi: 10.1016/j.agrformet.2017.06.016
[8] Ryu Y, Sonnentag O, Nilson T, et al. How to quantify tree leaf area index in an open savanna ecosystem:a multi-instrument and multi-model approach[J]. Agricultural and Forest Meteorology, 2010, 150(1):63-76. doi: 10.1016/j.agrformet.2009.08.007
[9] Bréda N J J. Ground-based measurements of leaf area index:a review of methods, instruments and current controversies[J]. Journal of Experimental Botany, 2003, 54(392):2403-2417. doi: 10.1093/jxb/erg263
[10] 刘志理, 金光泽.光学仪器法测定叶面积指数季节变化的误差分析[J].林业科学, 2016, 52(9):11-21. http://d.old.wanfangdata.com.cn/Periodical/lykx201609002 Liu Z L, Jin G Z. Bias analysis of seasonal changes of leaf area index derived from optical methods[J]. Scientia Silvae Sinicae, 2016, 52(9):11-21. http://d.old.wanfangdata.com.cn/Periodical/lykx201609002
[11] 温一博, 常颖, 范文义.基于MISR数据大兴安岭地区叶面积指数反演及尺度转换验证研究[J].北京林业大学学报, 2016, 38(5):1-10. doi: 10.13332/j.1000-1522.20150204 Wen Y B, Chang Y, Fan W Y. Algorithm for leaf area index inversion in the Great Xing' an Mountains using MISR data and spatial scaling for the validation[J]. Journal of Beijing Forestry University, 2016, 38(5):1-10. doi: 10.13332/j.1000-1522.20150204
[12] She W W, Zhang Y Q, Qin S G, et al. Habitat effect on allometry of a xeric shrub (Artemisia ordosica Krasch) in the Mu Us Desert of northern China[J]. Forests, 2015, 6(12):4529-4539. doi: 10.3390/f6124385
[13] 阎伟, 宗世祥, 骆有庆, 等.逐步回归模型在油蒿钻蛀性害虫预测中的应用[J].北京林业大学学报, 2009, 31 (3): 140-144. http://j.bjfu.edu.cn/article/id/8580 Yan W, Zong S X, Luo Y Q, et al. Application of stepwise regression model in predicting the movement of Artemisia ordosica boring insects[J]. Journal of Beijing Forestry University, 2009, 31(3): 140-144. http://j.bjfu.edu.cn/article/id/8580
[14] 张军, 黄永梅, 焦会景, 等.毛乌素沙地油蒿群落演替的生理生态学机制[J].中国沙漠, 2007, 27(6):977-983. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsm200706012 Zhang J, Huang Y M, Jiao H J, et al. Ecophysiological characteristics of Artemisia ordosica community succession in Mu Us Sandy grassland[J]. Journal of Desert Research, 2007, 27(6):977-983. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgsm200706012
[15] 张德魁, 王继和, 马全林, 等.油蒿研究综述[J].草业科学, 2007, 24(8):30-35. doi: 10.3969/j.issn.1001-0629.2007.08.006 Zhang D K, Wang J H, Ma Q L, et al. Summary of Artemisia ordosica studies[J]. Pratacultural Science, 2007, 24(8):30-35. doi: 10.3969/j.issn.1001-0629.2007.08.006
[16] 张军红, 吴波, 雷雅凯, 等.毛乌素沙地油蒿植株形态与结构特征分析[J].西南林业大学学报, 2011, 31(5):6-9. doi: 10.3969/j.issn.2095-1914.2011.05.002 Zhang J H, Wu B, Lei Y K, et al. Analysis of Artemisia ordosica plant morphology and structure characteristics in Mu Us Sandland[J]. Journal of Southwest Forestry University, 2011, 31(5):6-9. doi: 10.3969/j.issn.2095-1914.2011.05.002
[17] Xie J, Zha T S, Jia X, et al. Irregular precipitation events in control of seasonal variations in CO2 exchange in a cold desert-shrub ecosystem in northwest China[J]. Journal of Arid Environments, 2015, 120:33-41. doi: 10.1016/j.jaridenv.2015.04.009
[18] 唐思凌, 贾昕, 郭建斌, 等.沙蒿(Artemisia ordosica)叶面积指数的测定及模拟[J].生态学杂志, 2014, 33(2):547-554. http://d.old.wanfangdata.com.cn/Periodical/stxzz201402040 Tang S L, Jia X, Guo J B, et al. Measuring and modeling leaf area index for Artemisia ordosica[J]. Chinese Journal of Ecology, 2014, 33(2):547-554. http://d.old.wanfangdata.com.cn/Periodical/stxzz201402040
[19] 孔德胤, 杨松, 黄淑琴, 等.河套地区玉米叶面积指数的动态模拟[J].中国农业气象, 2014, 35(3): 281-286. doi: 10.3969/j.issn.1000-6362.2014.03.007 Kong D Y, Yang S, Huang S Q, et al. Dynamic simulation of maize leaf area index in Hetao region[J]. Chinese Journal of Agrometeorology, 2014, 35(3):281-286. doi: 10.3969/j.issn.1000-6362.2014.03.007
[20] Huemmrich K F, Black T A, Jarvis P G, et al. High temporal resolution NDVI phenology from micrometeorological radiation sensors[J/OL]. Journal of Geophysical Research:Atmospheres, 1999, 104(D22):27935-27944[2017-06-13]. http://onlinelibrary.wiley.com/doi/10.1029/1999JD900164/citedby.
[21] Wilson T B, Meyers T P. Determining vegetation indices from solar and photosynthetically active radiation fluxes[J]. Agricultural and Forest Meteorology, 2007, 144(3):160-179. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ021944193/
[22] 杜占池, 杨宗贵, 崔骁勇.内蒙古典型草原地区5类植物群落叶面积指数的比较研究[J].中国草地学报, 2001, 23(5):13-18. http://d.old.wanfangdata.com.cn/Periodical/zgcd200105003 Du Z C, Yang Z G, Cui X Y. A comparative study on leaf area index of five plant communities in typical steppe region of Inner Mongolia[J]. Chinese Journal of Grassland, 2001, 23(5):13-18. http://d.old.wanfangdata.com.cn/Periodical/zgcd200105003
[23] Huang L, Zhang Z S, Li X R. Sap flow of Artemisia ordosica and the influence of environmental factors in a revegetated desert area:Tengger Desert, China[J]. Hydrological Processes, 2010, 24(10):1248-1253. https://www.researchgate.net/publication/229927469_Sap_flow_of_Artemisia_ordosica_and_the_influence_of_environmental_factors_in_a_revegetated_desert_area_Tengger_Desert_China
[24] Xiao B, Hu K L. Moss-dominated biocrusts decrease soil moisture and result in the degradation of artificially planted shrubs under semiarid climate[J]. Geoderma, 2017, 291:47-54. doi: 10.1016/j.geoderma.2017.01.009
[25] Asner G P, Scurlock J M O, Hicke J A. Global synthesis of leaf area index observations:implications for ecological and remote sensing studies[J]. Global Ecology and Biogeography, 2003, 12(3):191-205. doi: 10.1046/j.1466-822X.2003.00026.x
[26] 邵璞, 曾晓东. CLM3.0-DGVM中植物叶面积指数与气候因子的时空关系[J].生态学报, 2011, 31(16):4725-4731. http://d.old.wanfangdata.com.cn/Periodical/gyqx200803001 Shao P, Zeng X D. Spationtemporal relationship of leaf area index simulated by CLM3.0-DGVM and climatic factors[J]. Acta Ecologica Sinica, 2011, 31(16):4725-4731. http://d.old.wanfangdata.com.cn/Periodical/gyqx200803001
[27] 麻雪艳, 周广胜.春玉米最大叶面积指数的确定方法及其应用[J].生态学报, 2013, 33(8):2596-2603. http://d.old.wanfangdata.com.cn/Periodical/stxb201308031 Ma X Y, Zhou G S. Method of determining the maximum leaf area index of spring maize and its application[J]. Acta Ecologica Sinica, 2013, 33(8):2596-2603. http://d.old.wanfangdata.com.cn/Periodical/stxb201308031
[28] Wang Q, Adiku S, Tenhunen J, et al. On the relationship of NDVI with leaf area index in a deciduous forest site[J]. Remote Sensing of Environment, 2005, 94(2):244-255. doi: 10.1016/j.rse.2004.10.006
[29] Bertram A, Clasen A, Kleinschmit B, et al. Estimation of the seasonal leaf area index in an alluvial forest using high-resolution satellite-based vegetation indices[J/OL]. Remote Sensing of Environment, 2014, 141:52-63[2017-06-21]. http://www.sciencedirect.com/science/article/pii/S0034425713003866.
[30] Street L E, Shaver G R, Williams M, et al. What is the relationship between changes in canopy leaf area and changes in photosynthetic CO2 flux in arctic ecosystems?[J]. Journal of Ecology, 2007, 95(1):139-150. doi: 10.1111/jec.2007.95.issue-1
[31] Maire G L, François C, Soudani K, et al. Forest leaf area index determination:a multiyear satellite-independent method based on within-stand normalized difference vegetation index spatial variability[J]. Journal of Geophysical Research:Biogeosciences, 2006, 111(G2): 1967-1970. http://cn.bing.com/academic/profile?id=f8035fbb83829f900cae2b5bb00b54d9&encoded=0&v=paper_preview&mkt=zh-cn
[32] Schleppi P, Thimonier A, Walthert L. Estimating leaf area index of mature temperate forests using regressions on site and vegetation data[J]. Forest Ecology and Management, 2011, 261(3):601-610. doi: 10.1016/j.foreco.2010.11.013
[33] 张学艺, 李剑萍, 官景得, 等.两种叶面积指数动态模拟方法的对比研究[J].国土资源遥感, 2011, 23(3):43-47. http://d.old.wanfangdata.com.cn/Periodical/gtzyyg201103008 Zhang X Y, Li J P, Guan J D, et al. A comparative study of two dynamic simulated methods for spring wheat leaf area index in Ningxia irrigation area[J]. Remote Sensing for Land & Resources, 2011, 23(3):43-47. http://d.old.wanfangdata.com.cn/Periodical/gtzyyg201103008
[34] 刘珺, 庞鑫, 杜灵通.夏玉米叶面积指数遥感反演研究[J].农业机械学报, 2016, 47(9):309-317. http://d.old.wanfangdata.com.cn/Periodical/nyjxxb201609042 Liu J, Pang X, Du L T. Study on leaf area index inversion of summer maize using remote sensing[J]. Transaction of the Chinese Society for Agricultural Machinery, 2016, 47(9):309-317. http://d.old.wanfangdata.com.cn/Periodical/nyjxxb201609042
[35] 彭虓, 张树文.基于NDVI与LAI的水稻生长状况研究[J].东北测绘, 2002, 25(4):16-19. doi: 10.3969/j.issn.1672-5867.2002.04.010 Peng X, Zhang S W. Research on rice growth status based on NDVI and LAI[J]. Northeast Surveying and Mapping, 2002, 25(4):16-19. doi: 10.3969/j.issn.1672-5867.2002.04.010
[36] Fan L, Gao Y, Brück H, et al. Investigating the relationship between NDVI and LAI in semi-arid grassland in Inner Mongolia using in-situ measurements[J]. Theoretical and Applied Climatology, 2009, 95(1-2):151-156. doi: 10.1007/s00704-007-0369-2
[37] López-Serrano F R, García-Morote A, Andrés-Abellán M, et al. Site and weather effects in allometries:a simple approach to climate change effect on pines[J]. Forest Ecology and Management, 2005, 215(1-3):251-270. doi: 10.1016/j.foreco.2005.05.014
[38] Wang X P, Fang J Y, Tang Z Y, et al. Climatic control of primary forest structure and DBH-height allometry in Northeast China[J]. Forest Ecology and Management, 2006, 234(1):264-274. doi: 10.1016-j.foreco.2006.07.007/
[39] Maestre F T, Salguerogómez R, Quero J L. It is getting hotter in here:determining and projecting the impacts of global environmental change on drylands[J]. Philosophical Transactions of the Royal Society of London, 2012, 367:3062-3075. doi: 10.1098/rstb.2011.0323