本通报所用的NDVI数据主要是NASA提供的MODIS NDVI数据,通报所用NDVI数据覆盖的时间从2002年1月至2013年9月底。通报只选择MODIS Terra星陆地 3/4级产品(54),空间分辨率为1公里。在原始数据下载之后,用全球耕地分布数据对NDVI进行掩膜处理,剔除非耕地数据,确保NDVI数据集适于农作物长势监测及估产等研究。
此外,本通报还使用了VITO(55)提供的,基于SPOT-VEGETATION传感器,分辨率为0.185°的长时间序列(1999年-2012年)NDVI平均数据。
本通报生产的气温产品为覆盖全球(0.25×0.25度)的月产品,产品时间范围为2000年1月至2013年9月。该产品是基于全球变化总目录资料库中的全球地表日数据集(GSOD)生产的(56)。GSOD数据集是由美国国家气候中心(NCDC)生产,由地表小时数据集(ISH)DSI-3505 (C00532)插值得到。可下载的在线数据最早可追溯至1929年,最新数据在日观测数据获取1-2日内更新。该数据集可提供来自全球9000多站点的气温、露点温度、海面气压、风速、降雨、雪深等观测参量。
本通报使用了GSOD数据集中2000年1月至2013年9月的日平均气温计算出月平均气温,考虑高程对温度的影响,结合STRM_DEM数据(57)使用克里金插值法得到了全球0.25×0.25度的月气温产品。
光合有效辐射(PAR)是影响作物产量的一个重要参数,PAR指的是波长在400-700纳米之间的太阳辐射,在不同的光周期内PAR和太阳短波辐射有一定的转换关系(58)。所以光合有效辐射可以由太阳短波辐射计算得到。2001-2012年全球月平均PAR数据来自NASA(美国航空航天局)的CERES下行短波辐射产品(59)。2013年全球月平均PAR数据由Geoland-2项目的全球一小时的下行短波辐射产品累加得到,在Geoland-2产品不能覆盖的区域,引入了风云卫星(FY-2D)的地表辐射产品进行补充(60-62)。所有上述产品都被转换到了WGS-84投影坐标系下,空间分辨率统一为0.25度。
CropWatch通过集成不同来源的降水产品,获取了2000年10月至2013年9月,覆盖整个陆地表面的,空间分辨率为0.25°×0.25°,时间分辨率为月的降水时间序列。该降水数据集使用的两大数据源包括:(i) TRMM降水数据集:由NASA提供的第7版空间分辨率为0.25°×0.25°,覆盖南北纬50°之间地区的TRMM降水产品(63);(ii) GPCC降水数据集(64):通过重采样技术将其由0.5°与1°空间分辨率采样成0.25°。2000年10月-2010年12月之间的CropWatch降水时间序列中南北纬50°范围的降水量采用TRMM月降水产品数据(TRMM 3B43),其他区域的采用的是GPCC数据。2011年1月-2013年6月之间,南北纬50°之外的地区采用的是GPCC监测分析数据(65)。2013年7月-2013年9月,南北纬50°之间的地区,通过合成时间分辨率为3小时,空间分辨率为0.25°×0.25°的TRMM 3B42近实时降水产品数据而来,其他区域采用的是GPCC First_guess产品(66) 。将CropWatch降水产品与实测站点的月降水量进行比较,结果表明二者具有较高的一致性。由于采用不同的数据源合成降水时间序列数据,在一些月份二者的拼接处有明显的(切割线)现象。
植被健康指数(VHI)可以有效地指示作物生长状况。本通报采用温度状态指数(TCI)和植被状态指数(VCI)加权的方法计算植被健康指数(67) (68) (69) 。温度状态指数和植被状态指数数据均可以通过NOAA数据中心的卫星数据应用和研究数据库下载 (70).该方法的计算方程如下:
但是,本文使用的经验系数(a)在全球不同区域是变化的,而不是前面研究者使用的常数0.5。
通报根据Lieth的“迈阿密”模型(71) (72)计算了净初级生产力(NPP),以此作为全球潜在生物量的度量。在迈阿密模型中,NPP可表达为温度和降雨两个环境因子的函数,如下式所示:
其中,TC 为年平均气温,单位为摄氏度(°C),Prec表示年平均降雨,单位为毫米(mm).