CN103971011A - 全球陆地植被生物量变化监测方法 - Google Patents
全球陆地植被生物量变化监测方法 Download PDFInfo
- Publication number
- CN103971011A CN103971011A CN201410212907.4A CN201410212907A CN103971011A CN 103971011 A CN103971011 A CN 103971011A CN 201410212907 A CN201410212907 A CN 201410212907A CN 103971011 A CN103971011 A CN 103971011A
- Authority
- CN
- China
- Prior art keywords
- biomass
- vegetation
- global
- gpp
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种全球陆地植被生物量变化监测方法,能够应用在气候变化研究、环境监测、土地管理、农情监测、以及灾害监测等遥感应用部门。该方法包含三个步骤:第一步骤是对原始数据进行处理和空间插值处理,然后用等式1对Terra卫星上MODIS的GPP/NPP产品计算空间平均值。第二个步骤是计算全球或者区域植被平均生物量状况。第三步骤是利用第二步计算的全球或者区域植被平均生物量值进一步计算全球或者区域植被生物量变化趋势,对全球或者区域植被生物量变化进行时空变化分析。全球分析表明,全球生物量整体下降趋势。
Description
技术领域
本发明涉及一种利用对地观测卫星上MODIS Gross/Net Primary Production (GPP/NPP)产品监测全球生物量变化的方法,从而有利于人类了解全球陆地不同区域的生物量变化情况,从而采取相应措施,该方法突破了传统方法利用站点对局部观测的局限,直观地展示全球陆地不同区域生物量变化情况。能够应用在气象、农业、环境监测等部门。
背景技术
陆地植被生产力是人类主要食物,纤维和燃料的最原始材料。定量计算全球植被的生物量变化,对预测气候变化、减缓和适应气候变化,研究全球碳循环与气候变化具有非常重要的意义。传统方法是通过直接或间接测定全球植被的生产力与生物现存量,再乘以不同植被生物量中碳元素的含量推算而得。由于全球植被观测受时间和空间,以及地表情况的影响,全球地面观测站点数量有限,而且不是均匀分布,特别是在山区和极低地区,地面观测和分类插值得到的结果不是非常好,精度不是很高[Aber, J.D., Melillo, J.M., Nadelhoffer, K.J., Pastor, J. and Boone, R.D., 1991. Factors controlling nitrogen cycling and nitrogen saturation in northern temperate forest ecosystems. Ecological Applications, 1(3): 303-315. Baldocchi, D., Hicks, B.B. and Meyers, T.P., 1988. Measuring biosphere-atmosphere exchanges of biologically related gases with micrometeorological methods. Ecology, 69: 1331-1340. Baldocchi, D., Valentini, R., Running, S., Oechel, W. and Dahlman, R., 1996. Strategies for measuring and modeling carbon dioxide and water vapour fluxes over terrestrial ecosystems. Global Change Biology, 2: 159-168.]。这种传统估算方法是基于植被(森林、草地和农作物等)调查数据,相关报道比较多,优点是直接、明确、技术简单,但其耗时耗力,且只能应用于小尺度的研究,难以及时反映全球陆地生态系统碳储量的动态变化[Asrar, G., Fuchs, M., Kanemasu, E.T. and Hatfield, J.H., 1984. Estimating absorbed photosynthetic radiation and leaf area index from spectral reflectance in wheat.Agronomy Journal, 76: 300-306.]。
卫星遥感技术已经被用来监测全球陆地表面,对利用卫星遥感技术监测地表生物量的变化已经有很多研究[Asrar, G., Myneni, R.B. and Choudhury, B.J., 1992. Spatial heterogeneity in vegetation canopies and remote sensing of absorbed photosynthetically active radiation: a modeling study. Remote Sensing of Environment, 41: 85-103.],已经提出了很多估算算法,美国对地观测卫星提供了GPP/NPP产品[Reeves MR, Zhao M, Heinsch FA and Running SW. (in prep). Characterizing moisture driven rangeland biomass fluctuations using MODIS primary productivity estimates.]。MODIS 遥感器于1999和2002年搭载对地观测卫星发射成功,为全球和区域资源环境动态监测开辟了又一新的途径。MODIS 是一个拥有36个波段的中分辨率遥感系统,每天可获得4次全球观测数据(1:30,10:30,13:30, 22:3),其飞行与太阳同步,并且是免费接收,因此非常适合全球生物量监测。
发明内容
本发明的目的在于提供一种从遥感数据MODIS GPP/NPP产品对全球生物量时空变化进行的方法,以克服现有大区域生物量时空变化监测的实际困难,直观地向人们展示全球不同区域生物量的变化,为气候变化研究提供植被相应的高精度数据,以便人类更好地适用气候变化和应对粮食安全。
为实现上述目的,本发明提供的从遥感数据MODIS GPP/NPP产品估算全球和区域生物量时空变化方法为:
第一步:对原始数据进行处理和空间插值处理,然后用等式1对Terra卫星上MODIS的GPP/NPP产品计算空间平均值
(1)
式中GPP mj 是像元j植被覆盖的平均状况,i是一年中的天数。
第二步:等式2计算全球或者区域植被平均生物量状况
(2)
式中是植被的平均状况,i是每天天数,j是像元数,S(j)像元j面积权重函数,是第i天j象元的植被状态。
第三步:利用第二步计算的全球或者区域植被平均生物量值和式3进一步计算全球或者区域植被生物量变化趋势,对全球或者区域植被生物量变化进行时空变化分析。
(3)
式中GPP_Rate表示植被生物量变化率,k表示第几年,GPPmk 是第k年的植被生物量,n表示这个年份。在MODIS产品中,GPP是以8天为单位进行输出,NPP是以年为单位进行输出。我们在区域估算中,对生物量的估算可以采用以天为计算单位,也可以直接利用8天产品进行计算。NPP直接采用公式3进行变化率计算。
本发明的有益效果是,利用MODIS GPP/NPP产品对全球陆地植被生物量变化进行变化监测分析,有效地克服了以往传统方法站点观测数量不足,分布不均,且校正不一致的缺点。为气候变化研究,粮食安全等监测等提供了有效手段和技术支撑。其操作实用性比传统的利用地面观察站点插值得要简单,面上精度要高。事实上,地面气象观测站也是本方法进一步提高精度的数据重要补充来源,二者结合将大大提高全球植被生物量的变化监测精度。
附图说明
下面结合附图和实施例对本发明进一步说明。
附图说明
图1是2003-2010年(a)全球陆地植生物量平均分布图,(b)全球陆地植被生物量变化率图
图2是2003至2011年全球陆地植被生物量时间变化率图
图3是2003至2011年世界主要大洲植被生物量变化图(a)北美洲,(b) 南美洲,(c) 亚洲,(d) 非洲
图4是2003至2011年世界主要国家植被生物量变化图(a)中国,(b)俄罗斯,(c)美国,(d)加拿大
具体实施方式
这里提供一个计算2001年至2011年全球平均水汽计算的一个实例:
本实例实现(方法)主要包括三个步骤:
第一步:对原始数据进行处理和空间插值处理,然后用等式1对Terra卫星上MODIS上午(10:30)的GPP/NPP产品计算空间平均值
第二步:等式2计算全球或者区域植被生物量平均生物量状况
第三步:利用第二步计算的全球或者区域植被生物量平均值和式3进一步计算和分析全球或者区域植被生物量变化趋势,对全球或者区域生物量变化进行时空变化分析。
图1是2003-2010年(a)全球陆地植生物量平均分布图,(b)全球陆地植被生物量变化率图,从图1a可以看出,全球生物量最高的地区分布在南美洲北部;从图1b看出,生物增加最明显的地方在南美洲和北高纬度地区,减少最明显的地方是非洲中部。
图2是2003至2011年全球陆地植被生物量时间变化率图,图中可以看出,全球生物量整体呈下降趋势。
图3是2003至2011年世界主要大洲植被生物量变化图(a)北美洲,(b) 南美洲,(c) 亚洲,(d) 非洲,非洲呈明显下降趋势。
图4是2003至2011年世界主要国家植被生物量变化图(a)中国,(b)俄罗斯,(c)美国,(d)加拿大,俄罗斯和加拿大植被生物量呈平稳增加,中国和美国植被生物量波动比较大。
Claims (1)
1. 全球陆地植被生物量变化监测方法,其步骤为:
第一步:对原始数据进行处理和空间插值处理,然后用等式1对Terra卫星上MODIS上午(10:30)的GPP/NPP产品计算空间平均值
(1)
式中GPP mj 是像元j植被覆盖的平均状况,i是一年中的天数;
第二步:等式2计算全球或者区域植被平均生物量状况
(2)
式中是植被的平均状况,i是每天天数,j是像元数,S(j)像元j面积权重函数,是第i天j象元的植被状态;
第三步: 利用第二步计算的全球或者区域植被平均生物量值和式3进一步计算全球或者区域植被生物量变化趋势,对全球或者区域植被生物量变化进行时空变化分析
(3)
式中GPP_Rate表示植被覆盖变化率,k表示第几年,GPPmk 是第k年的植被覆盖状态,n表示这个年份,在MODIS产品中,GPP是以8天为单位进行输出,NPP是以年为单位进行输出, 在区域估算中,对生物量的估算可以采用以天为计算单位,也可以直接利用8天产品进行计算,NPP直接采用公式3进行变化率计算,分析表明,全球生物量整体呈下降趋势。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410212907.4A CN103971011A (zh) | 2014-05-20 | 2014-05-20 | 全球陆地植被生物量变化监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410212907.4A CN103971011A (zh) | 2014-05-20 | 2014-05-20 | 全球陆地植被生物量变化监测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103971011A true CN103971011A (zh) | 2014-08-06 |
Family
ID=51240499
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410212907.4A Pending CN103971011A (zh) | 2014-05-20 | 2014-05-20 | 全球陆地植被生物量变化监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103971011A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285196A (zh) * | 2018-11-01 | 2019-01-29 | 云南省农业科学院质量标准与检测技术研究所 | 一种花毛茛花色变化特征确定方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103616490A (zh) * | 2013-12-03 | 2014-03-05 | 中国科学院南京地理与湖泊研究所 | 一种大型浅水湖泊水华蓝藻总存量估算的方法 |
CN103678914A (zh) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | 基于卫星遥感数据的高寒草地土壤呼吸估算方法 |
CN103761446A (zh) * | 2014-02-18 | 2014-04-30 | 中国农业科学院农业资源与农业区划研究所 | 从modis温度产品估算全球和区域平均温度方法 |
-
2014
- 2014-05-20 CN CN201410212907.4A patent/CN103971011A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103616490A (zh) * | 2013-12-03 | 2014-03-05 | 中国科学院南京地理与湖泊研究所 | 一种大型浅水湖泊水华蓝藻总存量估算的方法 |
CN103678914A (zh) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | 基于卫星遥感数据的高寒草地土壤呼吸估算方法 |
CN103761446A (zh) * | 2014-02-18 | 2014-04-30 | 中国农业科学院农业资源与农业区划研究所 | 从modis温度产品估算全球和区域平均温度方法 |
Non-Patent Citations (3)
Title |
---|
冯金朝等: "《生态学实验》", 31 October 2011, 中央名族大学出版社 * |
周珺: ""基于遥感数据的重庆市净初级生产力(NPP)时空特征研究"", 《中国硕士论文全文数据库工程科技Ⅰ辑》 * |
邢旗等: ""应用MODIS-NDVI对草地植被变化监测研究-以锡林郭勒盟为例"", 《草地学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285196A (zh) * | 2018-11-01 | 2019-01-29 | 云南省农业科学院质量标准与检测技术研究所 | 一种花毛茛花色变化特征确定方法及系统 |
CN109285196B (zh) * | 2018-11-01 | 2020-11-13 | 云南省农业科学院质量标准与检测技术研究所 | 一种花毛茛花色变化特征确定方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gerhards et al. | Challenges and future perspectives of multi-/Hyperspectral thermal infrared remote sensing for crop water-stress detection: A review | |
Gao et al. | Assessing the variability of corn and soybean yields in central Iowa using high spatiotemporal resolution multi-satellite imagery | |
Siebert et al. | Impact of heat stress on crop yield—on the importance of considering canopy temperature | |
Zhao et al. | Assimilating remote sensing information with crop model using Ensemble Kalman Filter for improving LAI monitoring and yield estimation | |
Pan et al. | Joint assimilation of leaf area index and soil moisture from sentinel-1 and sentinel-2 data into the WOFOST model for winter wheat yield estimation | |
Morel et al. | Coupling a sugarcane crop model with the remotely sensed time series of fIPAR to optimise the yield estimation | |
An et al. | Precipitation and minimum temperature are primary climatic controls of alpine grassland autumn phenology on the Qinghai-Tibet Plateau | |
Knipper et al. | Using high-spatiotemporal thermal satellite ET retrievals for operational water use and stress monitoring in a California vineyard | |
Yan et al. | A drought monitoring operational system for China using satellite data: design and evaluation | |
Cui et al. | Classification management for grassland using MODIS data: a case study in the Gannan region, China | |
Monteiro et al. | Sugarcane yield gap: can it be determined at national level with a simple agrometeorological model? | |
Conesa et al. | Terrestrial and remote indexes to assess moderate deficit irrigation in early-maturing nectarine trees | |
Fernández-Novales et al. | Monitoring and mapping vineyard water status using non-invasive technologies by a ground robot | |
Senay | Modeling landscape evapotranspiration by integrating land surface phenology and a water balance algorithm | |
Wang et al. | Grassland productivity response to droughts in northern China monitored by satellite-based solar-induced chlorophyll fluorescence | |
Peereman et al. | Assessing typhoon-induced canopy damage using vegetation indices in the Fushan Experimental Forest, Taiwan | |
Huang et al. | Comparison of remote sensing yield estimation methods for winter wheat based on assimilating time-sequence LAI and ET | |
He et al. | Satellite-observed energy budget change of deforestation in northeastern China and its climate implications | |
Alordzinu et al. | Rapid estimation of crop water stress index on tomato growth | |
García-Martín et al. | Spatiotemporal analysis of the frost regime in the Iberian Peninsula in the context of climate change (1975–2018) | |
Gu et al. | Crop water stress index as a proxy of phenotyping maize performance under combined water and salt stress | |
Gerardo et al. | Sentinel-2 Satellite Imagery-Based Assessment of Soil Salinity in Irrigated Rice Fields in Portugal | |
Li et al. | Surface temperature correction in TVDI to evaluate soil moisture over a large area | |
Balbontín et al. | Irrigation performance assessment in table grape using the reflectance-based crop coefficient | |
Liu et al. | Rapid extraction of regional-scale agricultural disasters by the standardized monitoring model based on google earth engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140806 |
|
RJ01 | Rejection of invention patent application after publication |