青藏高原草地地上生物量和理论载畜量

准确模拟和预测草地地上生物量(Aboveground biomass,AGB)和理论载畜量对于维持草地生态系统平衡、优化放牧管理至关重要。当前很多研究以围栏外草地AGB为基础,估算了青藏高原草地AGB的现存量。但是,牛羊啃食后的草地AGB现存量无法准确评估草地理论载畜量。围栏内草地不受家畜采食影响,其年际变率由环境因子驱动,可视为草地潜在AGB (potential AGB,AGB_p),更适用于草地理论载畜量的评估。本研究以青藏高原345个围栏内AGB观测数据为基础,结合气候、土壤和地形数据,利用随机森林算法构建草地潜在地上生物量估算模型,并对当前气候条件(2000-2018年)和未来20年(2021-2040年)4种气候变化情景(SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5)下的草地AGB_p和高寒草地理论载畜量进行模拟与预测。结果表明:(1)随机森林算法可准确模拟当前气候条件下的青藏高寒草地AGB_p(R²=0.76,P<0.001);2000-2018年青藏高寒草地AGB_p平均值为102.4 g m^-2,时间上增加趋势不明显(P>0.05);AGB_p年际波动和生长季降水显著正相关(R²=0.57,P<0.001),和生长季温度日较差显著负相关(R²=0.51,P<0.001)。(2)当前气候条件下,青藏高寒草地平均理论载畜量为0.94 SSU ha^-1(standardized sheep unit ha^-1);在过去20年约有54.1%草地理论载畜量呈提升状态。(3)和当前相比,未来20年青藏高原中部和北部草地AGB_p和理论载畜量呈下降态势。因此,建议未来在厘清气候变化影响下草畜关系的基础上进行有针对性的草牧业规划和管理,以缓解区域气候变化引起的草畜矛盾。     

基于EOS/MODIS数据的NDVI与EVI比较研究

作为NOAA/AVHRR 归一化植被指数(NDVI) 的延续和发展, EOS/MODIS 归一化植被指 数(NDVI) 和增强植被指数( EVI) 在许多领域得到广泛应用。应用数理统计和地统计学方法对二 者进行的对比研究表明: NDVI 在植被生长旺盛期容易达到饱和, 而EVI 则能克服这一现象, 比 较真实地反映植被的生长变化过程; 相同空间分辨率下, EVI 取值范围、标准差与变异系数均高 于NDVI, NDVI 数据比较均一, 其空间相关性高于EVI, EVI 更能反映研究区域内植被空间差异。 关键词:MODIS; 归一化植被指数(NDVI) ; 增强植被指数( EVI) ; 对比     

利用MODIS数据进行旱情动态监测研究

MODIS-EVI植被指数与不同覆盖程度植被的线性关系已得到明显改善,可以更有效地反映地表植被的生长状态。利用MODIS合成数据MYD11A2和MYD13A2获取的增强型植被指数(EVI)和陆地表面温度(Ts)构建EVI-Ts特征空间,并以该特征空间计算的温度植被干旱指数(TVDI)作为干旱监测指标,分析广西2006年秋旱分布。结果表明:以地表温度和MODIS-EVI为基础的温度植被干旱指数能较好地反映区域旱情分布和旱情发展过程,2006年9月中旬—11月中旬广西受旱区域不断扩大,旱情持续加重。     

Estimation of Daily Vapor Pressure Deficit Using MODIS Potential Evapotranspiration on the Tibetan Plateau

Vapor pressure deficit (VPD) is an important parameter in modelling hydrologic cycles and vegetation productivity. Meteorological stations are scarce in remote areas, which often results in imprecise estimations of VPD on the Tibetan Plateau. Moderate Resolution Imaging Spectroradiometer (MODIS) provides evapotranspiration data, which may offer the possibility of scaling up VPD estimations on the Tibetan Plateau. However, no studies thus far have estimated VPD using MODIS evapotranspiration data on the Tibetan Plateau. Therefore, this study used MODIS potential evapotranspiration (PET) to estimate VPD in alpine meadows, alpine steppes, croplands, forests and shrublands for the year, spring, summer, autumn and winter in 2000-2012. A series of root-mean- squared-error (RMSE) and mean-absolute-error (MAE) values were obtained for correlating measured VPD and estimated VPD using MODIS PET data for each listed time period and vegetation type: whole year (0.98-2.15 hPa and 0.68-1.44 hPa), spring (0.95-2.34 hPa and 0.72-1.54 hPa), summer (1.39-2.60 hPa and 0.89-1.96 hPa), autumn (0.78-1.93 hPa and 0.56-1.36 hPa), winter (0.48-1.40 hPa and 0.36-0.98 hPa), alpine steppes (0.48- 1.39 hPa and 0.36-1.00 hPa), alpine meadows (0.58-1.39 hPa and 0.44-0.90 hPa), croplands (1.10-2.55 hPa and 0.82-1.74 hPa), shrublands (0.98-1.90 hPa and 0.78-1.37 hPa), and forests (1.40-2.60 hPa and 0.98-1.96 hPa), respectively. Therefore, MODIS PET may be used to estimate VPD, and better results may be obtained if future studies incorporate vegetation types and seasons when the VPD data are estimated using MODIS PET on the Tibetan Plateau.     

基于MODIS数据的青藏高原旱情监测研究

本文利用温度植被旱情指数（TVDI）和植被供水指数（VSWI）分别对2009、2010年3—10月青藏高原土壤湿度状况进行监测分析,同时利用气象台站实测地面降水资料进行了验证。利用MODIS资料提取的归一化植被指数（NDVI）和地表温度（T_S）,构建NDVI-T_S特征空间,依据该特征空间计算出的反映青藏高原土壤湿度的TVDI与同期累积降水相关性显著;VSWI计算过程简单,但所反映的土壤湿度与同期累积降水的相关性较差。因此,对青藏高原这种范围广、下垫面多变复杂区域而言,TVDI能够更好地反映土壤湿度状况,对干旱监测具有一定的科学意义。     

Multifunctionality and Thresholds of Alpine Grassland on the Tibetan Plateau

Ecosystems can simultaneously provide multiple functions and services. Knowledge on the combinations of such multi-dimensional functions is critical for accurately assessing the carrying capacity and implementing sustainable management. However, accurately quantify the multifunctionality of ecosystems remains challenging due to the dependence and close association among individual functions. Here, we quantified spatial patterns in the multifunctionality of alpine grassland on the Tibetan Plateau by integrating four important individual functions based on data collected from a field survey and remote sensing NDVI. After mapping the spatial pattern of multifunctionality, we extracted multifunctionality values across four types of grassland along the northern Tibet Plateau transect. Effects of climate and grazing intensity on the multifunctionality were differentiated. Our results showed that the highest values of multifunctionality occurred in the alpine meadow. Low values of multifunctionality were comparable in different types of grassland. Annual precipitation explained the large variation of multifunctionality across the different types of grassland in the transect, which showed a significantly positive effect on the multifunctionality. Grazing intensity further explained the rest of the variation in the multifunctionality (residuals), which showed a shift from neutral or positive to negative effects on multifunctionality across the different types of grassland. The consistently rapid declines of belowground biomass, SOC, and species richness resulted in the collapse of the multifunctionality as bare ground cover amounted to 75%, which corresponded to a multifunctionality value of 0.233. Our results are the first to show the spatial pattern of grassland multifunctionality. The rapid decline of the multifunctionality suggests that a collapse in the multifunctionality can occur after the vegetation cover decreases to 25%, which is also accompanied by rapid losses of species and other individual functions. Our results are expected to provide evidence and direction for the sustainable development of alpine grassland and restoration management.     

基于植被盖度和高度的不同退化程度高寒草地地上生物量估算

由于气候变化和不合理的人类活动,20世纪80年代以来青藏高原高寒草地发生严重退化。地上生物量是评价草地退化的直观指标。通常采用植被盖度和高度来估算草地地上生物量,但草地退化后,植被盖度和高度与地上生物量之间的关系是否会发生变化目前还不清楚,这影响着退化草地生物量估算的精度。通过多元回归分析研究了青藏高原中部和东北部高寒草甸、高寒草原在不同退化程度下植被盖度和高度与地上生物量的关系。结果表明：（1）高寒草甸与高寒草原地上生物量整体上及不同退化阶段都没有显著差异（P>0.05）。（2）随着退化程度的加剧植被盖度和高度对地上生物量的影响也发生改变,体现在未退化阶段地上生物量主要受植被高度影响,退化后主要受植被盖度影响。（3）无论是高寒草甸还是高寒草原分退化程度的回归模型估算结果都较不分退化程度模型估算的生物量更接近实测值。我们建议在退化高寒草地研究中采用盖度和高度估算生物量时,根据退化阶段采用不同的估算模型。     

三江源区MODIS植被指数时空分布特征

利用16天合成的250m分辨率的MODIS数据,对三江源区NDVI和EVI植被指数时空分布特征进行了系统性的研究。结果表明：NDVI和EVI空间分布呈现出由东南向西北逐步下降的趋势。植被指数随季节呈现出规律变化,最高值出现在8月。不同海拔高度两种植被指数对植被覆盖情况的指示作用比较一致。植被类型的NDVI与EVI的差异...     

Driving Mechanism of Gross Primary Production Changes and Implications for Grassland Management on the Tibetan Plateau

The contribution of climatic change and anthropogenic activities to vegetation productivity are not fully understood. In this study, we determined potential climate-driven gross primary production (GPPp) using a process-based terrestrial ecosystem model, and actual gross primary production (GPPa) using MODIS Approach in alpine grasslands on the Tibetan Plateau from 2000 to 2015. The GPPa was influenced by both climatic change and anthropogenic activities. Gross primary production caused by anthropogenic activities (GPPh) was calculated as the difference between GPPp and GPPa. Approximately 75.63% and 24.37% of the area percentages of GPPa showed increasing and decreasing trends, respectively. Climatic change and anthropogenic activities were dominant factors responsible for approximately 42.90% and 32.72% of the increasing area percentage of GPPa, respectively. In contrast, climatic change and anthropogenic activities were responsible for approximately 16.88% and 7.49% of the decreasing area percentages of GPPa, respectively. The absolute values of the change trends of GPPp and GPPh of meadows were greater than those of steppes. The GPPp change values were greater than those of GPPh at all elevations, whereas both GPPp and GPPh showed decreasing trends when elevations were greater than or equal to 5000 m, 4600 m and 4800 m in meadows, steppes and all grasslands, respectively. Climatic change had stronger effects on the GPPa changes when elevations were lower than 5000 m, 4600 m and 4800 m in meadows, steppes and all grasslands, respectively. In contrast, anthropogenic activities had stronger effects on the GPPa changes when elevations were greater than or equal to 5000 m, 4600 m and 4800 m in meadows, steppes and all grasslands, respectively. Therefore, the causes of actual gross primary production changes varied with elevations, regions and grassland types, and grassland classification management should be considered on the Tibetan Plateau.     

CO2 processes in an alpine grassland ecosystem on the Tibetan Plateau

In this paper, the CO2 concentrations profile from 1.5 m depth in soil to 32 m height in atmosphere were measured from July 2000 to July 2001 in an alpine grassland ecosystem located in the permafrost area on the Tibetan Plateau, which revealed that CO2 concentrations varied greatly during this study period. Mean concentrations during the whole experiment in the atmosphere were absolutely lower than the CO2 concentrations in soil, which resulted in CO2 emissions from the alpine steppe soil to the atmosphere. The highest CO2 concentration was found at a depth of 1.5 m in soil while the lowest CO2 concentration occurred in the atmosphere. Mean CO2 concentrations in soil generally increased with depth. This was the compositive influence of the increasing soil moistures and decreasing soil pH, which induced the increasing biological activities with depth. Temporally, the CO2 concentrations at different layers in air remained a more steady state because of the atmospheric turbulent milking. During the seasonal variations, CO2 concentrations at surface soil interface showed symmetrical patterns, with the lowest accumulation of CO2 occurring in the late winter and the highest CO2 concentration in the growing seasons.     

青藏高原高寒草原生态系统CO2过程研究

In this paper, the CO2 concentrations profile from 1.5 m depth in soil to 32 m height in atmosphere were measured from July 2000 to July 2001 in an alpine grassland ecosystem located in the permafrost area on the Tibetan Plateau, which revealed that CO2 concentrations varied greatly during this study period. Mean concentrations during the whole experiment in the atmosphere were absolutely lower than the CO2 concentrations in soil, which resulted in CO2 emissions from the alpine steppe soil to the atmosphere. The highest CO2 concentration was found at a depth of 1.5 m in soil while the lowest CO2 concentration occurred in the atmosphere. Mean CO2 concentrations in soil generally increased with depth. This was the compositive influence of the increasing soil moistures and decreasing soil pH, which induced the increasing biological activities with depth. Temporally, the CO2 concentrations at different layers in air remained a more steady state because of the atmospheric turbulent milking. During the seasonal variations, CO2 concentrations at surface soil interface showed symmetrical patterns, with the lowest accumulation of CO2 occurring in the late winter and the highest CO2 concentration in the growing seasons.     

Temporal Variations in Growth and Aboveground Biomass of Phragmites australis and EVI Analysis in Jiaozhou Bay Coastal Salt Marshes,China

To better ascertain leaf, stem and flower traits, and analyze aboveground allocation during the vegetation period, three sampling vegetation transects were settled and reed samples were collected at intervals to determine morphological parameters and dry and wet weights in Jiaozhou Bay wetland. Remote sensing data were also combined to evaluate regional biomass through EVI regression models. Results show that growth dynamics of leaves and stems changed during the vegetation period. Stem length increased rapidly and peaked in September (194.40 ± 23.89 cm), whereas leaf width peaked in July. There was a significantly negative correlation between stem length and stem diameter with a value of -0.785. Stem biomass was higher than that of leaves, and the maximum value of aboveground biomass was 27.17 ± 3.56 g. F/C exhibited a tendency to increase and values ranged from 0.37-0.76. The aboveground biomass of sample plots reached a peak of 2356 ± 457 g/m² in September. EVI was 0.05-0.5; EVI and biomass had a better fitting effect using the power-exponent model compared with other models and its function was y = 4219.30 x^0.88 (R² = 0.7810). R² of the other three models ranked as linear function > polynomial function > exponent function, with the values being 0.7769, 0.7623 and 0.6963, respectively. EVI can be used to estimate vegetation biomass, and effectively solved the problems of the destructive effect to sample plots resulting from traditional harvest methods.     

青藏高原和阿尔卑斯山山体效应的对比研究

山体效应不仅对气候产生重大影响,也对区域地理生态格局有深远影响,尤其是它对山地垂直带分布和结构类型等的影响已经为地理学家和地植物学家所认识。目前相关研究主要集中在山体效应定量化方面,缺少不同山地山体效应的对比研究,因此对山体效应的区域差异性了解不足。本文选择欧亚大陆上具有明显山体效应的两个山地青藏高原和阿尔卑斯山为研究对象,利用收集到的气象台站观测数据、林线和DEM数据以及基于MODIS地表温度估算的青藏高原和阿尔卑斯山气温数据等,通过对比分析青藏高原与阿尔卑斯山相同海拔高度上的气温以及林线分布高度等来探讨两个山地的山体效应差异性。分析结果表明青藏高原的山体效应比阿尔卑斯山更为强烈,表现为：① 由于山体效应影响,在相同海拔高度上（4500 m）,青藏高原内部气温远高于阿尔卑斯山的气温,尤其是在最热月高原内部气温比阿尔卑斯山内部气温高10~15℃,在最冷月高原内部气温比阿尔卑斯山内部气温高5~10℃。② 由于山体效应影响,青藏高原内部林线也远高于阿尔卑斯山内部林线,约高2000~3000 m。本研究将为山体效应的影响因素分析奠定基础,同时对于揭示欧亚大陆山地生态系统格局具有一定的科学意义。     

Comparison of Methods for Evaluating the Forage-livestock Balance of Alpine Grasslands on the Northern Tibetan Plateau

Livestock grazing is one of primary way to use grasslands throughout the world, and the forage-livestock balance of grasslands is a core issue determining animal husbandry sustainability. However, there are few methods for assessing the forage-livestock balance and none of those consider the dynamics of external abiotic factors that influence forage yields. In this study, we combine long-term field observations with remote sensing data and meteorological records of temperature and precipitation to quantify the impacts of climate change and human activities on the forage-livestock balance of alpine grasslands on the northern Tibetan Plateau for the years 2000 to 2016. We developed two methods: one is statical method based on equilibrium theory and the other is dynamic method based on non-equilibrium theory. We also examined the uncertainties and shortcomings of using these two methods as a basis for formulating policies for sustainable grassland management. Our results from the statical method showed severe overgrazing in the grasslands of all counties observed except Nyima (including Shuanghu) for the entire period from 2000 to 2016. In contrast, the results from the dynamic method showed overgrazing in only eight years of the study period 2000-2016, while in the other nine years alpine grasslands throughout the northern Tibetan Plateau were less grazed and had forage surpluses. Additionally, the dynamic method found that the alpine grasslands of counties in the northeastern and southwestern areas of the northern Tibetan Plateau were overgrazed, and that alpine grasslands in the central area of the plateau were less grazed with forage surpluses. The latter finding is consistent with field surveys. Therefore, we suggest that the dynamic method is more appropriate for assessment of forage-livestock management efforts in alpine grasslands on the northern Tibetan Plateau. However, the statical method is still recommended for assessments of alpine grasslands profoundly disturbed by irrational human activities.     

Mapping Air Temperature in the Lancang River Basin Using the Reconstructed MODIS LST Data

Air temperature is an important climatological variable and is usually measured in meteorological stations. Accurate mapping of its spatial and temporal distribution is of great interest for various scientific disciplines, but low station density and complexity of the terrain usually lead to significant errors and unrepresentative spatial patterns over large areas. Fortunately the current studies have shown that the regression models can help overcome the problem with the help of time series remote sensing data. However, noise induced by cloud contamination and other atmospheric disturbances variability impedes the application of LST data. An improved Savizky-Golay (SG) algorithm based on the LST background library is used in this paper to reconstruct MODIS LST product. Data statistical analysis included 12 meteorological stations and 120 reconstructed MODIS LST images of the period from 2001 to 2010. The coeffificient of correlations (R²) for 80% of the stations was higher than 0.5 (below 0.5 for only 2 stations) which illustrated that there is a considerably close agreement between monthly mean TA (air temperature) and the reconstructed LST in the Lancang River basin. Comparing to the regression model for every month with only LST data, the regression model with LST and NDVI had higher R² and RMSE. Finally, the LST-NDVI regression method was applied as an estimate model to produce distributed maps of air temperature with month intervals and 1 km spatial in the Lancang River basin of 2010.     

归一化水汽指数与地表温度的关系

城市热岛的遥感研究,传统上主要利用归一化植被指数作为指标来评估地表温度-植被之间的关系。本研究探讨了以归一化水汽指数作为评估其与地表温度关系的适应性,比较了归一化水汽指数和归一化植被指数在表达城市热岛效应方面的效力。利用"单窗算法"和3个不同时段的TM/ETM+数据获取了珠江三角洲核心区域的地表温度分布图,分析了两种指数与地表温度之间的关系。对局部区域逐像元和对区域总体分析的结果表明:不同时段的地表温度与归一化水汽指数之间均存在较强的线性关系,而地表温度与归一化植被指数之间的线性关系相对较弱。说明在城市化的环境下,与传统上常用归一化植被指数作为指标来定量分析地表温度相比,归一化水汽指数提供了一种互补的、甚至更加有效的衡量指标。     

CO<Subscript>2</Subscript> processes in an alpine grassland ecosystem on the Tibetan Plateau

In this paper, the CO2 concentrations profile from 1.5 m depth in soil to 32 m height in atmosphere were measured from July 2000 to July 2001 in an alpine grassland ecosystem located in the permafrost area on the Tibetan Plateau, which revealed that CO2 concentrations varied greatly during this study period. Mean concentrations during the whole experiment in the atmosphere were absolutely lower than the CO2 concentrations in soil, which resulted in CO2 emissions from the alpine steppe soil to the atmosphere. The highest CO2 concentration was found at a depth of 1.5 m in soil while the lowest CO2 concentration occurred in the atmosphere. Mean CO2 concentrations in soil generally increased with depth. This was the compositive influence of the increasing soil moistures and decreasing soil pH, which induced the increasing biological activities with depth. Temporally, the CO2 concentrations at different layers in air remained a more steady state because of the atmospheric turbulent milking. During the seasonal variations, CO2 concentrations at surface soil interface showed symmetrical patterns, with the lowest accumulation of CO2 occurring in the late winter and the highest CO2 concentration in the growine seasons.     

NDVI、NDMI与地表温度关系的对比研究

通过研究归一化植被指数(NDVI)、归一化水汽指数(NDMI)与地表温度(LST)的相关关系,对比NDVI和NDMI定量分析LST的适宜程度。以Landsat 8遥感影像为数据源,以郑汴都市区为例,反演LST,计算NDVI和NDMI。从整体、分区、像元等不同层面分析NDVI和NDMI与LST的相关关系,并利用GEO-Da软件,分析150 m、300 m、450 m三种不同采样间隔数据的NDVI、NDMI与LST的空间相关性。主要结论为：①NDVI与LST的线性拟合度较差,而NDMI与LST具有很强的线性关系,剖面分析显示NDMI与LST呈显著负相关关系。②缓冲区分析结果表明,随着距城市中心距离的增加,用地类型增多,LST和NDMI之间的相关性逐渐增强。③在Moran's I空间相关性分析中,不同采样间隔下两指数与LST的负相关关系均比较明显,但由于水体在2个指数中数值的差异,使NDMI与LST表现出更强的空间负相关性,而NDVI与LST空间相关性则相对较弱。因此,总的来说NDMI是热环境研究的有效指标,与NDVI相比,NDMI与LST的负相关关系更强、更稳定,也更适宜于对LST的定量分析。     

汉江流域NDVI与水热指数时空变化及相关性分析

文章提出水热指数的概念和计算模型,运用经验正交函数和奇异值分解法对汉江流域1999-2010年的归一化植被指数及38个气象站点的气温和降水资料进行分析。结果表明:①汉江流域植被覆盖具有明显的空间差异,地形和农耕活动是最主要原因。②水热指数的空间分布也具有明显的规律,气温的纬度差异和雨带季节推移是根本原因;大地貌单元和季风气候是重要原因。其时间系数具有明显的冬夏差异,冬半年时间系数的绝对值及变率均较夏半年小,即水热组合冬半年差异小,夏半年差异大。③NDVI与水热指数之间呈现出高度正相关,但人类通过各种方式改变着地表覆盖和水热组合,使得两者之间的高度正相关在局部地区受到影响。     

古尔班通古特沙漠南缘荒漠化过程演变的景观格局特征分析

荒漠化过程是荒漠化研究的核心问题。以1977年MSS,1990年、2010年TM和2001年ETM+四个时期为遥感影像数据源,研究了近40 a古尔班通古特沙漠南缘荒漠化过程演变的景观格局特征。结果表明:(1)近40 a古尔班通古特沙漠南缘荒漠化过程演变经历了微弱逆转-加剧-逆转变化,荒漠化处于正逆交替动态变化中,总体上呈逆转趋势。(2)不同时期荒漠化过程演变格局特征呈波动性变化,以2001年为分水岭。近40 a来轻度和中度始终为主导作用,且有增强趋势。(3)年降水量与归一化植被指数(NDVI)相关系数达0.805,有明显正相关关系。年降水量对荒漠化具有显著作用,且其影响明显大于年均气温。(4)荒漠化过程景观格局特征对年降水量有明显响应。不同时期荒漠化过程演变中斑块密度随着降水量的增加总体呈减小趋势,年降水量达到184.4 mm,斑块密度趋于稳定,景观破碎度和景观异质性趋于稳定。随着降水的增加,轻度荒漠化类型破碎度呈逐渐减小趋势;中度荒漠化土地破碎度呈波动式变化;重度荒漠化破碎度波动下降趋势。(5)近40 a来气候变化是古尔班通古特沙漠的荒漠化与景观格局演化主要因素,但人为因素不容忽视,预示着未来干旱区荒漠化处于继续逆转趋势中。