Characterization of Vegetation Phenology in Response to Drought in the Middle and High Latitudes of China
-
摘要:
植被物候直接影响其生物量，调控生态系统碳循环过程。目前，气候变化(尤其干旱)对中国中高纬度植被物候的影响依然不清楚。因此，文章基于GIMMS NDVI3g数据集，提取中国30° N以北地区中多种植被类型生长季的开始日期(Start of the Season，SOS)与结束日期(End of the Season，EOS)两物候参数。然后结合野外观测数据，验证提取物候参数结果可靠性，并结合饱和水汽压差(Vapor Pressure Deficit, VPD)与改进后的标准化降水蒸散指数(Standardized Precipitation Evapotranspiration Index, SPEI)探究植被物候对干旱的响应特征规律。结果表明：(1)不同地区的植被物候变化呈现明显的差异性，单季植被与双季植被第1个生长季的SOS集中在每年的第30~180天，而双季植被第2个生长季的SOS集中在每年的第200~220天。单季植被与双季植被第1个生长季的EOS主要集中在每年的第180~300天，双季植被第2个生长季的EOS主要集中在每年的第260~300天。(2)森林季前VPD的上升导致植被的SOS提前及EOS延迟；草地季前VPD上升导致植被的SOS滞后以及EOS提前。(3)研究区内大部分地区的SPEI与植被的SOS、EOS均呈正相关，即干旱促使该地区植被的SOS、EOS提前。
-
关键词:
- 植被物候 /
- 植物生长季开始日期 /
- 植物生长季结束日期 /
- 饱和水汽压差 /
- 标准化降水蒸散发指数
Abstract:Vegetation phenology directly affects its biomass and regulates ecosystem carbon cycling processes. Currently, the impacts of climate change (especially drought) on vegetation phenology at mid- to high-latitudes in China are still unclear. Therefore, based on the GIMMS NDVI3g data, the Start of the Season (SOS) and End of the Season (EOS) parameters for the growing seasons of several vegetation types in the region north of 30° N were extracted. The extracted data were then combined with field observations to verify the effects of climate change on ecosystem carbon cycling. The results were then combined with field observations to verify the reliability of the parameters, and the saturated vapor pressure difference (VPD) and the improved standardized precipitation evapotranspiration index (SPEI) were used to investigate the characteristic patterns of vegetation phenology in response to drought. The results showed that: (1) the vegetation phenology varied significantly in different regions, with the SOS of the first growing season of mono- and bi-coastal vegetation concentrated between the 30th and 180th days of the year, and the SOS of the second growing season of bi-coastal vegetation concentrated between the 200th and 220th days of the year. The EOS spanning the 1st growing season of single-season and double-season vegetation was mainly concentrated between the 180th and 300th day of the year, and the EOS of the 2nd growing season of double-season crops was mainly concentrated between the 260th and 300th day of the year. (2) The rise of forest pre-season VPD leads to advancement of SOS and delay of EOS of vegetation; the rise of grass pre-season VPD leads to lagging of SOS and advancement of EOS of vegetation. (3) SPEI was positively correlated with SOS and EOS of vegetation in most parts of the study area, drought tended to contribute to the advancement of SOS and EOS of vegetation in the area.
-
植被物候是植被生长发育的周期性过程，在调节陆地生态系统碳平衡方面具有重要作用[1-2]，通常可通过观察植被生长、发育、衰老过程中的关键节点(如幼芽萌发时间、叶片展开时间、叶片颜色变化时间)和卫星物候指标(如生长季节的开始日期、高峰期和结束日期)来确定[3]。由于气候变化的加剧将导致全球和地区干旱发生显著变化，干旱会对地区植被物候造成明显影响[4]，因此，必须继续研究干旱与植被物候间的复杂作用机制，以了解未来生态系统物候的动态[5]。
随着卫星数据的广泛使用，全球大尺度物候学崭露头角。遥感方法可以更宏观、便捷地识别植被物候。一般通过定义生长季的开始日期(Start of the Season，SOS)和结束日期(End of the Season，EOS)来研究地区植被物候生长季长度及其对气候变化的响应特征[6-7]。植被物候遥感识别往往是基于植被指数或其他地表参量的季节变化而得。最常用的植被指数包括：归一化植被指数(Normalized Difference Vegetation Index，NDVI)[8]、增强型植被指数(Enhanced Vegetation Index，EVI)[9]、宽动态植被指数(Wide Dynamic Range Vegetation Index，WDRVI)[10]和陆地叶绿素指数(MERIS Terrestrial Chlorophyll Index，MTCI)[11]等。其中，NDVI是植被物候研究中首选的植被指数，其简便有效的计算方法在全球范围内得到了广泛认可[12]。
植被关键物候参数一般是对去噪后的植被指数时间序列进行提取。目前常用的关键物候参数确定的方法主要包括阈值法和植被指数变化检测法，其中，植被指数变化检测法通过直接检测植被指数时间序列曲线的变化特征来提取物候参数[12]，能有效确定植被物候的关键参数[13]。如TATEISHI和EBATA[14]将植被指数时间序列上升幅度最大、下降幅度最大的时间点分别确定为春季物候的开始日期、秋季物候的结束日期。然而，在过去的植被物候与气候影响研究中，通常忽略了双季植被，并且缺乏地面数据校验，物候识别精度有待提高[15-16]。
气候变化导致自然灾害越发频繁。干旱是全球最严重的自然灾害之一，具有复杂、随机、多因素引发的特点[17]，与植被物候变化关系密切。干旱指标是研究物候响应的基础，也是评估干旱对植被物候影响程度的重要指标。在干旱监测和研究物候对干旱的响应中，其定量化标准具有关键地位[18]。据统计，全球有上百种可供研究的干旱指标，各有优缺点，目前还没有一种适用于所有地区和所有时间尺度的指标[19]。标准化降水蒸散指数(Standardized Precipitation Evapotranspiration Index，SPEI)继承了标准化降水指数(Standardized Precipitation Index，SPI)多时间尺度的优点，并考虑了气温对蒸发的影响，能够较好满足多种类型的干旱评估需求[20-21]。同时作为区域干燥程度的气候指标，饱和水汽压差(Vapor Pressure Deficit，VPD)也是植被物候变化的前兆和标志[20-22]。
综上，本研究以中国中高纬度地区为研究区，综合利用GIMMS NDVI3g、气象站点数据和农气站点观测等多源数据，结合干旱指标SPEI和VPD，从更宏观、更全面的角度探讨干旱对植被物候的影响机制，以期为促进区域植被恢复、确保粮食产量与粮食安全提供科学参考。
1. 数据与方法
1.1 研究区概况
鉴于中国30° N以北地区(图 1)的农作物主要以一年一季或一年两季轮种为主，种植方式相对简单。因此，本研究重点探究中国30° N以北地区的植被物候特征及其与VPD、SPEI的响应特征。为了排除非植被土地利用类型变化的影响，研究利用MODIS卫星的土地覆盖气候模拟网格产品(https://modis.gsfc.nasa.gov/)，剔除了2001—2014年间发生土地利用变化的区域、非植被用地类型(如荒地、水体、建筑用地、永久冰雪和永久湿地)以及小于100个像元的植被用地类型。
![]()
图 1 研究区及其植被覆被类型分布图注：该图基于审图号为GS(2023)2767号的标准地图制作，底图无修改。Figure 1. Distribution of the study area and its vegetation cover types1.2 数据来源及预处理
本研究使用GIMMS NDVI3g数据集(半月合成最大值数据，空间分辨率为0.083 33°，来源：https://iridl.ldeo.columbia.edu/SOURCES)，计算了1982—2014年植被的关键物候参数(SOS、EOS)。使用样条插值法，将中国气象数据网(CMDN)提供的1982—2014年786个监测站的逐日降雨量和气温数据插值为栅格数据，以匹配研究区的物候参数(空间分辨率为0.833 33°，来源：http://data.cma.cn/)。地面植被物候数据来自中国气象数据网的农作物生长发育和农田土壤湿度旬值数据集(http://data.cma.cn/)，包含了1992—2013年在778个地点观测的不同植被的发育时期、株高等信息[23]。鉴于农作物种类繁多，为提高提取精度，本研究只选择研究区内种植面积较为广泛的小麦、玉米和水稻作为研究对象，剔除其他不满足条件的站点数据，最终选定了1992—2013年671个站点的小麦、玉米和水稻的生长记录。将地面数据中玉米的三叶期、小麦的出苗期和水稻的移栽期定义为SOS，而玉米、小麦和水稻的成熟期定义为EOS[13]。植被分类数据采用2001—2014年MODIS MCD12C1的IGBP分类数据(空间分辨率为0.833 33°，来源：https://modis.gsfc.nasa.gov/)。SPEI数据来源于中国区域日尺度SPEI数据集[21]。
1.3 植被物候提取与精度验证
1.3.1 Savizky-Golay滤波器
由于大气和云层等因素可能对植被物候提取精度产生影响, 为了提高准确性，本研究采用SG滤波器对每个像素的NDVI时间序列曲线进行重建。SG滤波的公式[24]如下：
Fi=12m+1m∑j=−mCjfi+j, (1) 其中：i是数据序列中数据点的索引, 表示当前正在计算的是哪一个数据点的平滑值；j是与i相关的一个位移索引, 用于在平滑窗口2m+1内对数据点进行遍历；fi、Fi分别为第i个原始数据、第i个过滤后的数据；Cj是滤波窗口的第j个权重值；m是滤波窗口大小的一半。最后, 使用二次样条法插值数据[25]。SG滤波器通过调整滤波窗口大小和多项式拟合阶数这2个参数，可以自由控制曲线的拟合程度[26]。
1.3.2 物候参数提取
鉴于植物生长具有阶段性，本研究通过关注植被生长季的开始日期(SOS)和结束日期(EOS)来探究研究区物候的时空特征。为了识别双季物候，本研究利用重构后的NDVI曲线，通过计算每个像素点的NDVI曲线的二阶导数，找到每个生长季的前半季的局部最大值点，将其作为SOS，而后半季的局部最大值点则被定义为EOS[27]。本研究使用地面实测的物候数据，计算了各省不同植被类型的平均SOS和EOS。为了方便处理数据，文中所有物候参数指标以时间序日(day of the year, DOY)表示，例如1月10日记为一年中的第10天。时间窗口的定义是各省不同植被的SOS和EOS平均值的前后30天[13]。
1.3.3 物候参数精度验证
为确保提取物候参数的准确性，本研究使用地面物候数据集进行验证；鉴于数据获取和目标研究区的植被生长状况，选择了3种植被类型：小麦、玉米和水稻；采用纳什决定系数(Nash-coefficient of determination，R2)[28]、纳什效率系数(Nash-Sutcliffe Efficiency，NSE)[29]和百分比偏差(Percentage Bias，PBIAS)[30]3个统计指标进行精度验证。3个统计指标的计算公式如下：
R2=[∑ni=1(Oi−ˉO)(Pi−ˉP)]2∑ni=1(Oi−ˉO)2∑ni=1(Pi−ˉP)2, (2) NSE=∑ni=1(Oi−ˉO)2−∑ni=1(Pi−ˉP)2∑ni=1(Oi−ˉO)2, (3) PBIAS =∑ni=1(Pi−Oi)∑ni=1Oi, (4) 其中，Pi是被评估成分的第i个预测值，Oi是第i个观察值，P是预测值的平均值，Q是观察值的平均值，n是观察总数。已有研究[31-32]表明: 当R2＞0.6、NSE＞0.5、PBIAS＜0.2时，精度结果可以满足本研究要求。
1.3.4 VPD计算
饱和蒸汽压差(Vapor Pressure Deficit，VPD)是饱和水汽压与实际水汽压之间的差值，能够有效反映区域的空气干燥程度，值越大代表空气越干燥[33]。计算公式如下：
VPD=es−ea, (5) 其中，es、ea分别是饱和水汽压、实际水汽压，计算公式如下：
es=0.61078e17.27∗TT+237.3, (6) ea=RH100∗es, (7) 其中，RH为空气相对湿度。
1.3.5 相关分析
本研究采用Pearson相关分析法，探索不同季前长度SOS、EOS与每个栅格的气候因素之间的相关关系。相关系数(r)的计算公式[34]如下：
r=∑ni=1(xi−ˉx)(yi−ˉy)√∑ni=1(xi−ˉx)2√∑ni=1(yi−ˉy)2, (8) 其中，x、y分别为VPD或SPEI的时间序列、植被物候时间序列，x、y分别为其平均值，n为时间序列长度。当数据服从正态分布时, 可使用如下公式检验相关系数r的显著性:
t=r∗√n−2√1−r2, (9) 其中，t是检验统计量，n是数据的个数。在给定的自由度n-2和显著性水平α下, 通过查t分布临界值表, 查找P值大小，本研究选择P＜0.05作为显著性水平标准，判断是否通过显著性检验。
2. 研究结果与分析
2.1 植被物候参数精度验证
由1982—2014年的地面观测物候与遥感植被物候的验证结果(表 1)可知：3种植被综合物候提取精度验证结果中R2、NSE和PBAIS分别为0.97、0.95和-5.3%。其中水稻的拟合效果最好，其R2、NSE和PBAIS分别为0.99、0.95和-0.6%。结果表明遥感提取的物候参数符合本研究对数据精度的要求。
表 1 植被物候提取精度Table 1. Precision of climate extract candidates作物类型 R2 NSE PBAIS/% 玉米 0.98 0.88 -6.5 小麦 0.97 0.90 -6.3 水稻 0.99 0.95 -0.6 综合 0.97 0.95 -5.3 2.2 植被物候空间分布特征
1982—2014年，中国30° N以北地区中植被的多年平均SOS的空间分布如图 2所示。由单季植被与双季植被第1个生长季的SOS空间分布(图 2A)可知：(1)单季植被与双季植被第1个生长季的SOS并无明显纬度分布差异，在不同植被类型之间差异较大。华北平原地区植被的SOS最早，内蒙古自治区植被的SOS较晚。单季植被与双季植被第1个生长季的SOS集中在每年的第30~180天。(2)内蒙古自治区西部、吉林省和辽宁省等地植被的SOS出现较晚，集中在每年的第110~185天，该地区的植被类型以草地为主。而内蒙古自治区西北部与黑龙江省交界处的植被的SOS主要集中在每年的第100~110天。(3)东三省南部地区植被的SOS相对较早，集中在每年的第80~100天。青海省、西藏自治区与四川省交界处的大部分地区植被的SOS晚于每年的第110天。由双季植被第2个生长季的SOS的空间分布(图 2B)可知：(1)从空间上看，双季植被主要集中在黄淮海平原，双季植被第2个生长季的SOS集中在每年的第200~220天。(2)除河南省西部、安徽省西部的部分地区外，其他地区双季植被第2个生长季的SOS均早于每年的第200天。
![]()
图 2 1982—2014年研究区植被的平均SOS空间分布图Figure 2. Spatial distribution of mean SOS of vegetation in the study area from 1982 to 20141982—2014年，中国30° N以北地区中植被的平均EOS的空间分布如图 3所示。由单季植被与双季植被第1个生长季的EOS的空间分布(图 3A)可知: (1)从内蒙古自治区中部逐渐向西至青海省、西藏自治区等地，植被主要以草地为主，植被的EOS大多出现在每年的第280天之后。(2)在江苏省北部和黄淮海平原等地以作物为主的区域，植被第1个生长季的EOS大多在每年的第220天之前。由双季植被第2个生长季的EOS的空间分布(图 3B)可知：双季植被以作物为主，主要分布在黄淮海平原，其植被的EOS从西南到东北逐渐推迟；双季植被的第2个生长季的EOS主要集中在每年的第260~300天。
![]()
图 3 1982—2014年研究区植被的平均EOS空间分布图Figure 3. Spatial distribution of mean EOS of vegetation in the study area from 1982 to 2014为了进一步探究不同植被类型物候期的差异性，本研究分析了不同植被类型SOS、EOS的差异。由植物生长季的SOS的箱型分布图(图 4A)可知: (1)森林、草地区域的SOS的时间整体差异不大，均集中在每年的第90~110天。(2)由于耕地种植的农作物不同，单双季种植植被间的SOS差异明显。其中早稻和冬小麦的SOS最早平均为每年的第53.64天和第59.43天，而夏玉米和晚稻作为第2季种植植被，它们的SOS也是最晚的，平均为每年的第194.04天和第193.54天。由植物生长季的EOS的箱型分布图(图 4B)可知: (1)早稻和冬小麦的EOS是最早的，平均为每年的第147.63天和第154.01天。(2)第2季植被与单季植被的EOS平均时间非常接近，集中在每年的第285~295天。
![]()
图 4 不同植被类型的SOS与EOS的箱型分布图注：图中数字为对应植被类型的SOS/EOS的平均值。Figure 4. Box distribution diagram of SOS and EOS in different vegetation types2.3 植被物候对干旱的响应特征
2.3.1 SOS对季前VPD的响应特征
1982—2014年，中国30° N以北地区植被的SOS与季前VPD间相关性的空间分布(图 5A)和统计结果(图 5B)表明: (1)研究区内，有64.48%的地区植被的SOS与季前VPD呈负相关，其中38.1%的地区的相关性结果通过了显著性检验。这些区域主要位于吉林省、黑龙江省和内蒙古自治区的东北部。(2)研究区内，有13.23%的地区植被的SOS与季前VPD呈显著正相关，主要分布在辽宁省与内蒙古自治区的交界地区。由于全球变暖，全球饱和水汽压差会呈现增加的趋势，将使得我国东北北部植被物候的SOS明显推迟，而我国的30° N~45° N的中纬度地区植被的SOS将提前。
![]()
图 5 植被的SOS与季前VPD相关性及显著性结果注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。Figure 5. Correlation and significance result diagram between SOS of vegetation and pre-season VPD2.3.2 EOS对季前VPD的响应特征
中国30° N以北地区植被的EOS与季前VPD间相关性、显著性的空间分布(图 6A)和统计结果(图 6B)表明: (1)研究区内，有66.05%的地区植被的EOS与季前VPD呈正相关，其中37%的地区的相关性结果通过了显著性检验，主要分布在吉林省、黑龙江省以及内蒙古自治区的东北部。(2)13.46%的地区植被的EOS与季前VPD呈显著负相关，主要集中在黑龙江省东部以及辽宁省与内蒙古自治区之间的交界地区。随着季前VPD的增加，我国大部分地区植被的EOS将提前，并且高纬度地区植被的EOS的提前情况更明显。
![]()
图 6 植被的EOS与季前VPD相关性及显著性结果注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。Figure 6. Correlation and significance results between EOS of vegetation and pre-season VPD2.3.3 SOS对SPEI的响应特征
由中国30° N以北地区植被的SOS与SPEI相关性、显著性的空间分布(图 7A)和统计结果(图 7B)可知: (1)研究区内，有58.88%的地区植被的SOS与SPEI呈正相关，而41.12%的地区为负相关。(2)研究区内，有17.50%的地区的相关性结果通过了显著性检验，这些区域主要分布在东北部的大兴安岭地区，呈现显著正相关的栅格占比为14.25%。总体而言，东北大兴安岭部分地区植被的SOS的响应相对较为敏感, 干旱多促使该地区植被的SOS提前。在吉林省、辽宁省与内蒙古自治区交界处，SPEI与植被的SOS之间则以负相关为主，干旱增加使得该地区植被的SOS被推迟。
![]()
图 7 植被的SOS与SPEI相关性及显著性结果注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。Figure 7. Correlation and significance result diagram between SOS of vegetation and SPEI2.3.4 EOS对SPEI的响应特征
由中国30° N以北地区植被的EOS与SPEI相关性、显著性的空间分布(图 8A)和统计结果(图 8B)可知: (1)SPEI与植被的EOS主要呈正相关，其面积占比约为67.49%，但只有6.06%的地区的相关性结果通过了显著性检验，主要集中在呼伦贝尔东部。(2)研究区内1.72%的地区植被的EOS与SPEI呈现显著负相关。总体而言，研究区内大部分地区，干旱发生将促使研究区内大部分地区植被的EOS提前。而在东北北部大部分地区植被的EOS随着干旱的发生将逐渐推迟。
![]()
图 8 植被的EOS与SPEI相关性及显著性结果注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。Figure 8. Correlation and significance result diagram between EOS of vegetation and SPEI3. 结论
本文以中国中高纬度区域为研究区，利用GIMMS NDVI3g长时间序列遥感数据集，采用Savizky-Golay滤波器与变点监测法，在栅格水平上提取不同植被类型的物候参数，并结合地面物候数据，验证遥感提取物候的精度，分析不同植被物候的时空演变特征。此外，基于气象数据计算VPD指数，结合改进的日SPEI数据集，表征研究区范围内植被物候对干旱的响应机制, 解析植被物候对干旱的响应特征。主要结论如下:
(1) 不同地区的植被物候变化呈现明显的差异性。单季植被与双季植被第1个生长季的SOS集中在每年的第30~180天，而双季植被第2个生长季的SOS集中在每年的第200~220天。单季植被与双季植被第1个生长季的EOS主要集中在每年的第180~300天，双季植被的第2个生长季的EOS主要集中在每年的第260~300天。
(2) VPD对植被物候的影响在不同地区呈现明显的空间差异。在以森林为主的吉林省、黑龙江省和内蒙古自治区东北部等地，季前VPD上升导致植被的SOS提前和EOS延迟；而在以草地为主的辽宁省与内蒙古自治区交界处，季前VPD上升导致植被的SOS滞后和EOS提前。
(3) 研究区内大部分地区的SPEI与植被的SOS和EOS呈正相关，即干旱促使植被的SOS和EOS提前。然而，这种变化并不显著，只有山西省北部地区的变化趋势通过显著性检验。
本研究提出基于遥感数据识别物候参数，结合VPD和SPEI较好地揭示了植被物候对干旱的响应特征，可以为区域遥感监测植被状况的变化、预警干旱可能对地区植被、农业造成的危害提供技术支持和理论依据。然而，植被生长物候是一个复杂而缓慢的过程，与气候因子之间关系复杂。本研究仅局限于干旱与植被物候间的响应关系，后续将从气温、降雨量和土壤湿度等其他可能影响植被物候的因素综合考虑，结合物理过程模型，从机理上探究植被物候对气候的响应机制。
-
![]()
图 1 研究区及其植被覆被类型分布图
注：该图基于审图号为GS(2023)2767号的标准地图制作，底图无修改。
Figure 1. Distribution of the study area and its vegetation cover types
![]()
图 2 1982—2014年研究区植被的平均SOS空间分布图
Figure 2. Spatial distribution of mean SOS of vegetation in the study area from 1982 to 2014
![]()
图 3 1982—2014年研究区植被的平均EOS空间分布图
Figure 3. Spatial distribution of mean EOS of vegetation in the study area from 1982 to 2014
![]()
图 4 不同植被类型的SOS与EOS的箱型分布图
注：图中数字为对应植被类型的SOS/EOS的平均值。
Figure 4. Box distribution diagram of SOS and EOS in different vegetation types
![]()
图 5 植被的SOS与季前VPD相关性及显著性结果
注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。
Figure 5. Correlation and significance result diagram between SOS of vegetation and pre-season VPD
![]()
图 6 植被的EOS与季前VPD相关性及显著性结果
注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。
Figure 6. Correlation and significance results between EOS of vegetation and pre-season VPD
![]()
图 7 植被的SOS与SPEI相关性及显著性结果
注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。
Figure 7. Correlation and significance result diagram between SOS of vegetation and SPEI
![]()
图 8 植被的EOS与SPEI相关性及显著性结果
注：A图中的黑点区域为相关性结果通过显著性检验的区域，其他区域为相关性结果未通过显著性检验的区域；B图中的百分数为统计的不同相关性结果面积占比。
Figure 8. Correlation and significance result diagram between EOS of vegetation and SPEI
表 1 植被物候提取精度
Table 1 Precision of climate extract candidates
作物类型 R2 NSE PBAIS/% 玉米 0.98 0.88 -6.5 小麦 0.97 0.90 -6.3 水稻 0.99 0.95 -0.6 综合 0.97 0.95 -5.3 
下载: 导出CSV
-
[1] CHEN M M, LOPEZ L, BHAVSAR S P, et al. What's hot about mercury? Examining the influence of climate on mercury levels in Ontario top predator fishes[J]. Environmental Research, 2018, 162(1): 63-73.
[2] HILL A C, VAZQUEZ L A, VARGAS R. Linking vegetation spectral reflectance with ecosystem carbon phenology in a temperate salt marsh[J]. Agricultural and Forest Meteorology, 2021, 307: 108481/1-14.
[3] 王敏钰, 罗毅, 张正阳, 等. 植被物候参数遥感提取与验证方法研究进展[J]. 遥感学报, 2022, 26(3): 431-455. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB202203001.htm WANG M Y, LUO Y, ZHANG Z Y, et al. Recent advances in remote sensing of vegetation phenology: retrieval algorithm and validation strateg[J]. National Remote Sensing Bulletin, 2022, 26(3): 431-455. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB202203001.htm
[4] ZHANG J Z, GYO X H, ALEXANDER M, et al. Drought limits wood production of Juniperus przewalskii even as growing seasons lengthens in a cold and arid environment[J]. Catena, 2021, 196: 104936/1-11.
[5] BADECK F W, BONDEAU A, BOTTCHER K, et al. Responses of spring phenology to climate change[J]. New Phytologist, 2004, 162(2): 295-309.
[6] 赵珍, 贾文雄, 张禹舜, 等. 祁连山区植被物候遥感监测与变化趋势[J]. 中国沙漠, 2015, 35(5): 1388-1395. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSS201505037.htm ZHAO Z, JIA W X, ZHANG Y S, et al. Monitoring and change trends of phenological phase in the Qilian Mountains based on remote sensing[J]. Journal of Desert Research, 2015, 35(5): 1388-1395. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSS201505037.htm
[7] 俎佳星, 杨健. 东北地区植被物候时序变化[J]. 生态学报, 2016, 36(7): 2015-2023. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201607025.htm ZU J X, YANG J. Temporal variation of vegetation phenology in northeastern China[J]. Acta Ecologica Sinica, 2016, 36(7): 2015-2023. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201607025.htm
[8] 马晓芳, 陈思宇, 邓婕, 等. 青藏高原植被物候监测及其对气候变化的响应[J]. 草业学报, 2016, 25(1): 13-21. https://www.cnki.com.cn/Article/CJFDTOTAL-CYXB201601003.htm MA X F, CHEN S Y, DENG J, et al. Vegetation phenology dynamics and its response to climate change on the Tibetan Plateau[J]. Acta Prataculturae Sinica, 2016, 25(1): 13-21. https://www.cnki.com.cn/Article/CJFDTOTAL-CYXB201601003.htm
[9] DESCALS A, VERGER A, YIN G, et al. Improved estimates of arctic land surface phenology using Sentinel-2 time series[J]. Remote Sensing, 2020, 12(22): 3738/1-13.
[10] SAKAMOTO T, WARDLOW B D, GITELSON A A, et al. A two-step filtering approach for detecting maize and soybean phenology with time-series MODIS data[J]. Remote Sensing of Environment, 2010, 114(10): 2146-2159. doi: 10.1016/j.rse.2010.04.019
[11] BOYD D S, ALMOND S, DASH J, et al. Phenology of vegetation in Southern England from Envisat MERIS terrestrial chlorophyll index (MTCI) data[J]. International Journal of Remote Sensing, 2011, 32(23): 8421-8447.
[12] TANG J, ZENG J Y, ZHANG Q, et al. Self-adapting extraction of cropland phenological transitions of rotation agroecosystems using dynamically fused NDVI images[J]. International Journal of Biometeorology, 2020, 64(8): 1273-1283. doi: 10.1007/s00484-020-01904-1
[13] LUO Y C, ZHANG Z, CHEN Y, et al. ChinaCropPhen1km: a high-resolution crop phenological dataset for three staple crops in China during 2000-2015 based on leaf area index (LAI) products[J]. Earth System Science Data, 2020, 12(1): 197-214. doi: 10.5194/essd-12-197-2020
[14] TATEISHI R, EBATA M. Analysis of phenological change patterns using 1982-2000 Advanced Very High Resolution Radiometer(AVHRR) data[J]. International Journal of Remote Sensing, 2004, 25(12): 2287-2300. doi: 10.1080/01431160310001618455
[15] REINERS P, ASAM S, FREY C, et al. Validation of AVHRR land surface temperature with MODIS and in situ LST—a timeline thematic processor[J]. Remote Sensing, 2021, 13(17): 3473/1-30.
[16] ZHAO J J, WANG Y Y, ZHANG Z X, et al. The variations of land surface phenology in Northeast China and its responses to climate change from 1982 to 2013[J]. Remote Sensing, 2016, 8(5): 400/1-23.
[17] HAO Z C, SINGH V P, XIA Y L. Seasonal drought prediction: advances, challenges, and future prospects[J]. Reviews of Geophysics, 2018, 56(1): 108-141. doi: 10.1002/2016RG000549
[18] 王永利, 侯琼, 苗百岭, 等. 内蒙古马铃薯干旱风险区划[J]. 应用气象学报, 2017, 28(4): 504-512. https://www.cnki.com.cn/Article/CJFDTOTAL-YYQX201704011.htm WANG Y L, HOU Q, MIAO B L, et al. Drought risk regionalization of potatoes in Inner Mongolia[J]. Journal of Applied Meteorological Science, 2017, 28(4): 504-512. https://www.cnki.com.cn/Article/CJFDTOTAL-YYQX201704011.htm
[19] FUNG K F, HUANG Y F, KOO C H, et al. Drought forecasting: a review of modelling approaches 2007-2017[J]. Journal of Water and Climate Change, 2020, 11(3): 771-799. doi: 10.2166/wcc.2019.236
[20] VICENTE-SERRANO S M, BEGUERIA S, LOPEZ-MORENO J I. A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index[J]. Journal of Climate, 2010, 23(7): 1696-1718. doi: 10.1175/2009JCLI2909.1
[21] WANG Q F, ZENG J Y, QI J Y, et al. A multi-scale daily SPEI dataset for drought characterization at observation stations over mainland China from 1961 to 2018[J]. Earth System Science Data, 2021, 13(2): 331-341. doi: 10.5194/essd-13-331-2021
[22] MARQUES-GALVEZ J E, MORTE A, NAVARRO-RODENAS A. Spring stomatal response to vapor pressure deficit as a marker for desert truffle fruiting[J]. Mycorrhiza, 2020, 30(4): 503-512. doi: 10.1007/s00572-020-00966-8
[23] TAN Q H, LIU Y J, DAI L, et al. Shortened key growth periods of soybean observed in China under climate change[J]. Scientific Reports, 2021, 11(1): 8197/1-12.
[24] 蔡天净, 唐瀚. Savitzky-Golay平滑滤波器的最小二乘拟合原理综述[J]. 数字通信, 2011, 38(1): 63-68. https://www.cnki.com.cn/Article/CJFDTOTAL-SZTX201101022.htm CAI T J, TANG H. A review of least squares fitting principles for Savitzky-Golay smoothing filters[J]. Digital Communication, 2011, 38(1): 63-68. https://www.cnki.com.cn/Article/CJFDTOTAL-SZTX201101022.htm
[25] KOIDE D, IDE R, OGUMA H. Detection of autumn leaf phenology and color brightness from repeat photography: accurate, robust, and sensitive indexes and modeling under unstable field observations[J]. Ecological Indicators, 2019, 106: 105482/1-11.
[26] 梁守真, 施平, 邢前国. MODIS NDVI时间序列数据的去云算法比较[J]. 国土资源遥感, 2011(1): 33-36. https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201101007.htm LIANG S Z, SHI P, XING Q G. A comparison between the algorithms for removing cloud pixel from MODIS NDVI time series data[J]. Remote Sensing for Natural Resources, 2011(1): 33-36. https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201101007.htm
[27] WANG Q F, TANG J, ZENG J Y, et al. Regional detection of multiple change points and workable application for precipitation by maximum likelihood approach[J]. Arabian Journal of Geosciences, 2019, 12(23): 745/1-16.
[28] PARAJULI P B, NELSON N O, FREES L D, et al. Comparison of AnnAGNPS and SWAT model simulation results in USDA-CEAP agricultural watersheds in south-central Kansas[J]. Hydrological Processes, 2008, 23(5): 748-763.
[29] NASH J E, SUTCLIFFE J V. River flow forecasting through conceptual models part Ⅰ-A discussion of principles[J]. Journal of Hydrology, 1970, 10(3): 282-290.
[30] MORADI M. Trend analysis and variations of sea surface temperature and chlorophyll-a in the Persian Gulf[J]. Marine Pollution Bulletin, 2020, 156: 111267/1-14.
[31] 张晓萱, 崔耀平, 刘素洁, 等. 自然植被物候遥感提取精度评估[J]. 生态学杂志, 2019, 38(5): 1589-1599. https://www.cnki.com.cn/Article/CJFDTOTAL-STXZ201905038.htm ZHANG X X, CUI Y P, LIU S J, et al. Evaluation of the accuracy of phenology extraction methods for natural vegetation based on remote sensing[J]. Chinese Journal of Ecology, 2019, 38(5): 1589-1599. https://www.cnki.com.cn/Article/CJFDTOTAL-STXZ201905038.htm
[32] 赵一静, 王晓利, 侯西勇, 等. 2003—2019年山东省冬小麦关键物候期时空特征[J]. 生态学报, 2021, 41(19): 7785-7795. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB202119027.htm ZHAO Y J, WANG X L, HOU X Y, et al. Spatio-temporal characteristics of key phenology of winter wheat in Shandong Province from 2003 to 2019[J]. Acta Ecologica Sinica, 2021, 41(19): 7785-7795. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB202119027.htm
[33] 张红梅, 吴炳方, 闫娜娜. 饱和水汽压差的卫星遥感研究综述[J]. 地球科学进展, 2014, 29(5): 559-568. https://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ201405005.htm ZHANG H M, WU B F, YAN N N. Remote sensing estimates of vapor pressure deficit: an overview[J]. Advance in Earth Science, 2014, 29(5): 559-568. https://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ201405005.htm
[34] PEARSON K. Note on regression and inheritance in the case of two parents[J]. Proceedings of the Royal Society of London, 1895, 58(1): 240-242.
计量
- 文章访问数: 55
- HTML全文浏览量: 22
- PDF下载量: 36
