共查询到20条相似文献，搜索用时 127 毫秒
1.
三江并流河源区植被覆盖度对气候要素的响应 总被引：2，自引：0，他引：2
基于AVHRR/NDVI和MODIS/NDVI遥感数据,通过拟合两种数据源延长NDVI时间序列来反演了三江并流河源区1982—2012年植被覆盖度空间格局的变化规律,并结合气温和降水量数据,从不同时空尺度上分析了植被覆盖变化趋势及其对气候变化的响应。结果表明:1.研究区植被覆盖度整体上随着热量梯度呈南高北低、东高西低的态势分布,且海拔在4 000 m以下的区域植被覆盖度较高,平均值在0.68~0.75;2.研究期间,植被覆盖度以0.02/(10 a)的变化率呈增加趋势,海拔在4 600 m左右和低覆盖度的区域植被增加趋势显著,显著增加区域面积占研究区总面积的36.4%,主要分布于沱沱河流域、怒江流域(除源头外)和研究区东南部等地区,这些区域植被覆盖度平均值仅为0.36;3.相对于降水量,研究区的植被覆盖度与气温要素相关性更为显著,特别是高原面上的植被覆盖度整体上与气温呈显著正相关关系,显示了强烈的热量限制性生态类型的特点。 相似文献
2.
基于RS的近30年滹沱河流域植被覆盖度动态变化研究 总被引：1，自引：0，他引：1
以滹沱河流域1984年、2001年和2014年3期遥感影像为数据源,基于像元二分模型,利用归一化植被指数(NDVI)反演流域植被覆盖度信息,分析了滹沱河流域近30年来植被覆盖度时空演变特征及其变化原因。滹沱河流域近30年来植被覆盖度总体呈增加趋势,其中高植被覆盖度区面积比例增幅最明显,增加了1.1倍。空间分布上,受气候和人类活动双重因素影响,流域不同区域植被覆盖度变化差异较显著,上中游植被总体以改善为主,海拔相对较低的区域则表现为退化,尤其是下游植被覆盖度持续下降,植被退化面积占该区的71.52%。 相似文献
3.
2000—2020年渭河流域植被覆盖度及产水量时空变化 总被引：1，自引：0，他引：1
渭河流域是黄土高原最早实施生态治理工程的区域,研究其植被覆盖度和产水量的变化及相互关系,对于黄河流域以及黄土高原来说具有重要意义。基于InVEST产水模型和Landsat遥感影像分析了流域产水量和植被覆盖度的变化特征。结果表明：渭河流域2000—2020年植被覆盖度整体呈显著增加趋势,空间上渭河流域的NDVI从东南至西北呈现减小—增大—减小的趋势;极显著改善区域的面积占整个流域的面积69.2%,植被覆盖度等级整体向高植被覆盖度等级转换,流域生态环境正在逐渐改善。产水量整体呈上升趋势,坡度较大、植被覆盖度较高的区域产水量较低,而坡度较缓、人口活动较为强烈的区域产水量高。除夏季外,各季节产水量整体呈增加趋势;受生态工程的影响,夏季产水量呈先增加后降低的趋势,冬季产水量呈先减少后增加的趋势。植被类型多为林地,坡度较大的区域,产水量和植被覆盖度呈显著负相关,在坡度较缓、植被类型多为草地的区域,产水量和植被覆盖度呈正相关。随着植被覆盖度增加,产水量呈增加趋势,当植被覆盖度达到阈值时,流域产水量不再增加,甚至出现降低的趋势。 相似文献
4.
针对中亚地区的强生态脆弱性、高敏感性特征,有必要开展广域、长期的植被覆盖监测以匹配“绿色丝绸之路”的可持续发展目标。鉴于此,联合Landsat 5和Landsat 8卫星数据集,利用Google Earth Engine（GEE）地理空间数据云计算平台,估算了中亚地区1993—2018年间共12期的植被覆盖度。结果表明：（1）中亚地区植被覆盖总体水平较低,但也具有较为显著的空间异质性。（2）中亚地区1993—2018年间多数区域植被覆盖趋势较为稳定,哈萨克斯坦丘陵、费尔干纳盆地等区域植被覆盖度呈增加趋势,乌拉尔河流域和锡尔河流域等区域植被覆盖趋势为负。（3）植被覆盖度时序特征上,中亚地区1993—2018年间总体植被覆盖度累积增加3%,其中吉尔吉斯斯坦和塔吉克斯坦植被覆盖分别增加3.96%和5.86%。（4）裸土区呈退缩趋势,面积总计减少25.9×104 km2,低植被覆盖区、中植被覆盖区和高植被覆盖区范围在呈现出的振荡式增加。研究结合遥感大数据和地理云计算对中亚地区进行区域尺度的植被覆盖动态监测,能对中亚地区生态评估和演替分析提供技术支持和定量数据。 相似文献
5.
2000-2015年石羊河流域植被覆盖度及其对气候变化的响应 总被引：1，自引：1，他引：0
内陆河流域植被覆盖度的敏感性是预测未来生物多样性变化的重要指标,是植被应对气候变暖的重要反馈。分析了2000-2015年MODIS-NDVI数据反演的植被覆盖度时空动态变化趋势,结合平均气温、降水量、日照时数、相对湿度、地面温度和蒸发量数据,研究了流域及各生态功能区的植被覆盖度与气候因子的相关性,探讨植被覆盖度变化过程中的气候因素制约方式,了解不同时空尺度下内陆河流域植被覆盖度在全球暖湿化过程中对气候的响应。结果表明：（1）石羊河流域平均植被覆盖度较低,上游的植被覆盖度59.4%,下游13.6%;2000-2015年,流域植被呈现改善趋势的面积远远大于退化的面积,盆地绿洲区植被覆盖度增加趋势最明显。流域植被总体恢复较好,但高海拔地区、城市和民勤绿洲的周边地区植被有不同程度的退化。（2）2000-2015年,石羊河流域各气候因子对植被覆盖度表现为不显著的相关关系,其中与降水量呈正相关的面积最大,与蒸发量呈负相关的面积最大;从上游到下游,植被生长与热量的相关程度逐渐变弱,与水分的相关程度则逐渐增强。（3）石羊河流域的植被覆盖度与气候因子的样条函数存在极显著的线性相关,水原涵养区和荒漠区的植被覆盖度对气候因子的响应较高;绿洲区的植被覆盖度对气候因子的响应相对较低。地面温度的变化是影响石羊河流域植被覆盖度空间格局变化的主要气候制约因素。 相似文献
6.
艾比湖流域NDVI垂直梯度变化特征 总被引：2，自引：0，他引：2
利用2000—2011年MODIS-NDVI数据分析了艾比湖流域植被变化特征,以MVC法合成年最大NDVI,研究不同尺度下海拔与NDVI的关系。建立了12年的NDVI变化趋势的一元线性回归方程,结合DEM数据,分析了24等分下海拔间的NDVI变化趋势。结果显示:(1)总体上,NDVI与海拔的拟合程度较好,判定系数均大于0.88,随着等分数目减少,判定系数增大,10等分时,海拔与NDVI的相关性最大,为0.97;2 500m为NDVI的分水岭,该海拔以下区域植被覆盖逐渐增加,以上区域植被覆盖逐渐减少;(2)研究区12年NDVI趋势变化平均值为1.03,总体属轻度改善;1 500m以上区域植被呈减少态势,1 200m以下植被呈增长态势;(3)各海拔带内NDVI趋势平均值为-42~54,研究区大部分像元NDVI趋势值较小,植被覆盖情况属于轻度改善和轻度退化。 相似文献
7.
《云南地理环境研究》2016,(2)
2003年来随着普洱市"边三县"(孟连、西盟、澜沧)大面积人工桉树林的引种,该地区植被覆盖发生了较大变化。以桉树引种较多的普洱市"边三县"为研究区、TM影像为数据源,在RS和GIS的支持下,通过计算并处理各期影像NDVI(归一化植被指数)值,并用MODIS数据处理结果进行验证的方法,对三县范围内的植被覆盖时空变化进行分析。结果表明:研究区在2000~2010年植被覆盖情况得到很大改善,2010~2014年整体植被覆盖度有下降的趋势;2000~2014年高植被覆盖区和低植被覆盖区相对不变,各乡镇植被覆盖变化趋势相同。 相似文献
8.
《干旱区地理》2021,44(3):700-707
为了补给塔里木河下游断流的河道和挽救濒死的荒漠群落,自2001年塔里木河实施了生态输水工程,流域生态环境得到初步改善,实现了流域经济社会发展与生态环境保护的双赢。在生态输水20 a之际,为了评估输水的生态效应,基于多源遥感数据MOD13Q1和MCD12Q1数据,通过数学统计方法调查了塔里木河下游归一化植被指数(NDVI)和植被覆盖度的时空变化。结果表明:经过近20 a的生态输水,塔里木河下游的NDVI从2000年的0.14增大到2020年的0.21,增幅达33.3%,植被面积从2000年492 km~2扩大到2020年的1423 km~2;其中,低、中、高植被覆盖度的面积相较于2000年分别增多了277 km~2,537 km~2和132 km~2,增幅分别达到20.8%,448.0%和190.0%,中段的植被面积和植被覆盖度的增幅均高于上、下段区域的增幅。NDVI和植被覆盖度均表现为距离河道2 km范围内的值较大,且增幅也较大,2 km以外NDVI和植被覆盖度逐渐变小,且增幅也较小。空间上来看,约有57.1%的区域NDVI和植被覆盖度呈现显著增加趋势,显著减小的区域都不超过2.0%。Hurst指数表明超过75.0%以上的区域这种变绿的趋势会在未来继续保持,仅有6.1%的区域生态会继续恶化。综合评估塔里木河输水工程的生态效应,对于干旱区荒漠的保育恢复治理、生态水资源管理等方面具有重要的科学意义。 相似文献
9.
植被NDVI对城市扩展及气候变化的响应——以西安及其附近区域为例 总被引：1，自引：0，他引：1
《干旱区地理》2017,(2)
植被NDVI对城市扩展及气候变化响应研究,对于科学评估区域生态环境变化及调整与约束人类活动具有重要理论和现实意义。以西安及其附近区域为例,基于区域土地利用、MODIS NDVI、气温和降水数据,分析了植被NDVI对城市扩展及气候变化的响应,结果表明:(1)2000-2014年研究区植被NDVI变化过程划分为2000-2007年的显著增加阶段和2007-2014年的显著减少阶段,前者主要分布于区域北部黄土高原、南部秦岭北坡,后者主要分布于区域中部关中平原尤其是西安及其附近区域。(2)2000-2015年研究区建设用地增加1 428.27 km2,建设用地增加区域植被NDVI呈显著减少趋势。(3)研究区植被NDVI与年降水量的相关性高于年平均气温,同时西安及其以南区域植被NDVI与年平均气温、年降水量均呈负相关关系,反映出城市扩展等人类活动对植被NDVI变化的影响超过了气候变化的影响。研究结果表明植被NDVI总体受气候变化控制,但局部受人类活动影响更为严重,并且植被NDVI对气候变化的响应表现出波动性,而对城市扩展表现出线性减少趋势性,为通过植被NDVI变化区分自然因素与人为因素对环境影响提供了可能。 相似文献
10.
2000-2015年毛乌素沙区植被覆盖度变化趋势 总被引：1，自引：1，他引：0
毛乌素沙区是中国生态安全屏障的重要组成部分,土地沙漠化严重,近年来国家和地方政府实施了一系列生态保护与建设工程。利用2000-2015年MODIS13Q1 NDVI产品、年平均气温、年降水量数据,采用回归分析方法和显著性检验,对毛乌素沙区植被覆盖度变化趋势及其对气温、降水变化的响应进行了分析研究。结果表明：（1）毛乌素沙区植被覆盖度总体由东向西呈减少趋势,大部分区域植被覆盖度在30%以下,沙地腹地依然明显存在极低植被覆盖的流沙区。（2）21世纪以来毛乌素沙区植被覆盖度总体呈增加趋势,但空间差异明显。中东部大部分地区及腹地流动沙带之间植被覆盖度呈显著、极显著增加趋势;西部西鄂尔多斯荒漠草原区植被覆盖度变化不显著;植被覆盖度变化显著或极显著减小的地区极少,呈点状零散分布。（3）毛乌素沙区植被覆盖度对气温和降水响应的敏感性存在空间差异,且存在时滞性差异。东部黄土高原过渡区和西部西鄂尔多斯荒漠草原区植被对降水和气温的响应敏感,植被覆盖度的变化与气温和降水因子呈显著相关关系;毛乌素沙地主体区植被覆盖度变化与当年总降水量和年平均气温相关性不强,但与时滞降水的相关系数显著增大,可能与毛乌素沙地土壤质地和植被类型对降水的分配与利用方式有关。 相似文献
11.
This paper used five years (2001-2006) time series of MODIS NDVI images with a 1-km spatial resolution to produce a land cover map of Qinghai Province in China. A classi-fication approach for different land cover types with special emphasis on vegetation, espe-cially on sparse vegetation, was developed which synthesized Decision Tree Classification, Supervised Classification and Unsupervised Classification. The spatial distribution and dy-namic change of vegetation cover in Qinghai from 2001 to 2006 were analyzed based on the land cover classification map and five grade elevation belts derived from Qinghai DEM. The result shows that vegetation cover in Qinghai in recent five years has been some improved and the area of vegetation was increased from 370,047 km2 in 2001 to 374,576 km2 in 2006. Meanwhile, vegetation cover ratio was increased by 0.63%. Vegetation cover ratio in high mountain belt is the largest (67.92%) among the five grade elevation belts in Qinghai Prov-ince. The second largest vegetation cover ratio is in middle mountain belt (61.80%). Next, in the order of the decreasing vegetation cover ratio, the remaining grades are extreme high mountain belt (38.98%), low mountain belt (25.55%) and flat region belt (15.46%). The area of middle density grassland in high mountain belt is the biggest (94,003 km2), and vegetation cover ratio of dense grassland in middle mountain belt is the highest (32.62%), and the in-creased area of dense grassland in high mountain belt is the greatest (1280 km2). In recent five years the conversion from sparse grass to middle density grass in high mountain belt has been the largest vegetation cover variation and the converted area is 15931 km2. 相似文献
12.
利用 2000—2015 年的 EOS／MODIS 数据，采用趋势分析、Hurst 指数、变异系数法对伊犁河 谷植被时空变化及未来趋势进行分析，结果显示：空间分布上，伊犁河谷植被覆盖度呈北部、南部、 东部偏高，西部、中部偏低的分布特征；时间变化上，2000—2015 年，伊犁河谷植被覆盖度波动减 小，减速为 6.25%·(10 a)-1；区域分布上，伊犁河谷植被表现为低波动变化，波动程度中等以及下占 73.16%，波动程度高的区域占 26.84%。未来预测表明，伊犁河谷植被覆盖呈退化趋势，其中，持续 退化的面积占 57.55%，持续改善的面积占 13.51%。 相似文献
13.
GIMMS NDVI database and geo-statistics were used to depict the spatial distribution and temporal stability of NDVI on the Mongolian Plateau.The results demonstrated that:(1) Regions of interest with high NDVI indices were distributed primarily in forested mountainous regions of the east and the north,areas with low NDVI indices were primarily distributed in the Gobi desert regions of the west and the southwest,and areas with moderate NDVI values were mainly distributed in a middle steppe strap from northwest to southeast.(2) The maximum NDVI values maintained for the past 22 years showed little variation.The average NDVI variance coefficient for the 22-year period was 15.2%.(3) NDVI distribution and vegetation cover showed spatial autocorrelations on a global scale.NDVI patterns from the vegetation cover also demonstrated anisotropy;a higher positive spatial correlation was indicated in a NW-SE direction,which suggested that vegetation cover in a NW-SE direction maintained increased integrity,and vegetation assemblage was mainly distributed in the same specific direction.(4) The NDVI spatial distribution was mainly controlled by structural factors,88.7% of the total spatial variation was influenced by structural and 11.3% by random factors.And the global autocorrelation distance was 1178 km,and the average vegetation patch length(NW-SE) to width(NE-SW) ratio was approximately 2.4:1.0. 相似文献
14.
基于MOD13Q1数据的宁夏生长季植被动态监测 总被引：1，自引：0，他引：1
宁夏自2000年后实施退耕还林以来,局部地区的生态环境得到明显改善,为探求近15来年宁夏地区植被的动态变化及其影响因子,本文以MOD13Q1为数据源,结合DEM数据、土地利用分类图,采用Sen+Mann-Kendall非参数检验方法和Hurst模型,分析了宁夏不同行政区、不同海拔、不同坡度、不同坡向及不同植被类型生长季NDVI的空间变化特征及未来变化趋势;并利用重心迁移模型和转移矩阵分析宁夏2000-2014年间植被覆盖的时空演变特征。结果表明：①从空间分布看,宁夏南部六盘山、北部贺兰山及引黄灌溉区植被长势较好,中部干旱地区植被长势较差;且植被NDVI与海拔高程和坡度呈显著正相关。②从植被覆盖的转移矩阵看,较高植被覆盖的面积占比从2000年的17.29%增长到2014年的31.55%,主要是由较低植被覆盖转化而来的。③从重心迁移方向看,中度植被覆盖和较高植被覆盖的重心迁移最为明显,分别向东北方向偏移了129.49 km和向东南方向偏移了89.49 km。④从变化趋势看,生长季植被NDVI整体呈上升趋势,明显改善的面积占总面积的59.63%,轻微改善区域占31.72%;林地和水田显著改善的面积分别占总面积的71.50%和70.80%;显著改善的面积比例随海拔高程和坡度的增加均先增加后减少,且南部各行政区植被改善的面积均高于北部。⑤从可持续性看,植被恢复的持续性较强,89.24%的植被NDVI呈现持续改善的趋势;南部地区的持续改善的面积大于北部地区。 相似文献
15.
蒙古高原NDVI的空间格局及分异 总被引：5，自引：2，他引：3
GIMMS NDVI database and geo-statistics were used to depict the spatial distribu-tion and temporal stability of NDVI on the Mongolian Plateau. The results demonstrated that: (1) Regions of interest with high NDVI indices were distributed primarily in forested moun-tainous regions of the east and the north, areas with low NDVI indices were primarily distrib-uted in the Gobi desert regions of the west and the southwest, and areas with moderate NDVI values were mainly distributed in a middle steppe strap from northwest to southeast. (2) The maximum NDVI values maintained for the past 22 years showed little variation. The average NDVI variance coefficient for the 22-year period was 15.2%. (3) NDVI distribution and vege-tation cover showed spatial autocorrelations on a global scale. NDVI patterns from the vegetation cover also demonstrated anisotropy; a higher positive spatial correlation was in-dicated in a NW-SE direction, which suggested that vegetation cover in a NW-SE direction maintained increased integrity, and vegetation assemblage was mainly distributed in the same specific direction. (4) The NDVl spatial distribution was mainly controlled by structural factors, 88.7% of the total spatial variation was influenced by structural and 11.3% by random factors. And the global autocorrelation distance was 1178 km, and the average vegetation patch length (NW-SE) to width (NE-SW) ratio was approximately 2.4:1.0. 相似文献
16.
为科学认识喀斯特山区植被变化及其地形效应,基于MODIS NDVI数据,采用统计学方法,系统分析2000-2016年喀斯特山区植被变化的时空特征及其与海拔、地形起伏度、坡度、坡向的关系。研究表明,黔桂喀斯特山区植被绿度中部高,西北及东南较低,年均NDVI随海拔和地形起伏度的增加呈单峰曲线变化,峰值位于400~600 m,NDVI随坡度和坡向的变化不明显;2000-2016年大部分地区NDVI呈增长趋势,其中超过20%的地区呈显著增长（P<0.05）,年均增长率约0.0018。西部和东南部绿化趋势最为显著,仅在东北和中东部,NDVI呈下降趋势;NDVI呈增长趋势的比例随海拔的增加而增加,说明该喀斯特山区近年来植被恢复向着良性化方向发展,高海拔植被恢复速率更快,低海拔缓坡处的植被生态建设需要进一步加强。 相似文献
17.
近5年青海省植被覆盖变化的遥感监测 总被引：3，自引：2，他引：1
This paper used five years （2001-2006） time series of MODIS NDVI images with a 1-km spatial resolution to produce a land cover map of Qinghai Province in China. A classification approach for different land cover types with special emphasis on vegetation, especially on sparse vegetation, was developed which synthesized Decision Tree Classification, Supervised Classification and Unsupervised Classification. The spatial distribution and dynamic change of vegetation cover in Qinghai from 2001 to 2006 were analyzed based on the land cover classification map and five grade elevation belts derived from Qinghai DEM. The result shows that vegetation cover in Qinghai in recent five years has been some improved and the area of vegetation was increased from 370,047 km^2 in 2001 to 374,576 km^2 in 2006. Meanwhile, vegetation cover ratio was increased by 0.63%. Vegetation cover ratio in high mountain belt is the largest （67.92%） among the five grade elevation belts in Qinghai Province. The second largest vegetation cover ratio is in middle mountain belt （61.80%）. Next, in the order of the decreasing vegetation cover ratio, the remaining grades are extreme high mountain belt （38.98%）, low mountain belt （25.55%） and flat region belt （15.46%）. The area of middle density grassland in high mountain belt is the biggest （94,003 km^2）, and vegetation cover ratio of dense grassland in middle mountain belt is the highest （32.62%）, and the increased area of dense grassland in high mountain belt is the greatest （1280 km^2）. In recent five years the conversion from sparse grass to middle density grass in high mountain belt has been the largest vegetation cover variation and the converted area is 15931 km^2. 相似文献
18.
2000~2014年黑龙江流域（中国）植被覆盖时空变化及其对气候变化的响应 总被引：2，自引：0，他引：2
采用MODIS/NDVI数据,利用Theil-Sen Median 趋势分析、Mann-Kendall 以及Hurst指数方法分析2000~2014年黑龙江流域（中国）植被的时空变化特征、植被变化发展趋势及可持续性特征;应用相关分析法研究了气候变化对植被生长的影响。结果表明,2000~2014年黑龙江流域（中国）植被NDVI指数呈缓慢增加趋势,山区植被覆盖增加显著,东北部平原区植被覆盖持续退化,总体上植被覆盖持续改善能力较弱。植被NDVI对气候响应的季节差异显著,且不同类型植被对气候因子的响应不一致：春季植被NDVI主要受温度影响,夏季植被NDVI主要受降水量影响,秋季林地NDVI与温度正相关、草地NDVI与降雨量正相关。 相似文献
19.
长白山区植被生长季NDVI时空变化及其对气候因子敏感性 总被引：7，自引：1，他引：6
本文利用长白山区SPOT/VGT NDVI数据和气象数据,分析该区不同植被类型NDVI时空变化特征以及与气候因子的相关关系,并探讨了植被对气候变化响应的滞后性。结果表明：①2000-2009年,长白山区植被NDVI逐年变化总体呈增长趋势,增长区域的面积占全区面积的83.91%,在空间上主要集中在北坡和西坡,NDVI减少区域集中在南坡;②NDVI变化率随季节和植被类型变化而不同,NDVI增长主要集中在5月和9月,而7月NDVI变化较小,甚至出现下降趋势;③植被NDVI与温度和降水存在着显著的正相关性（p〈0.01）,且NDVI与温度的相关性高于与降水的相关性,且随海拔升高,NDVI与温度相关性增强;④NDVI对气温和降水变化的响应存在滞后期,不同植被类型,滞后期存在差异。苔原NDVI对温度和降水响应的滞后期大约10天,而针阔混交林和针叶林ND-VI对温度和降水响应的滞后期约为20天。 相似文献
20.
科尔沁沙地沙漠化风险评价 总被引：2，自引：0，他引：2
基于气象、土地利用类型与植被覆盖指数等数据,采用层次分析法与加权综合评分法评价了科尔沁沙地沙漠化风险。结果表明：1995、2000、2010年科尔沁沙地沙漠化风险指数平均值分别为0.420、0.428、0.437,16年来科尔沁沙地沙漠化处于高度风险水平。科尔沁沙地沙漠化风险等级在空间分布上表现为北部与西部沙漠化风险等级低,南部与东部为沙漠化高风险区。1995-2010年科尔沁沙地各旗县沙漠化风险等级存在差异,开鲁县的沙漠化风险由高度风险转变为极高度风险等级,扎鲁特旗的沙漠化风险由中度风险转变为低度风险等级,敖汉旗、巴林右旗、科尔沁区与科尔沁左翼中旗的沙漠化风险等级波动变化,其余旗县的沙漠化风险等级未变化。1995-2010年科尔沁沙地沙漠化中度,高度风险等级的面积分别减少了16 544.76、4 223.25 km2,低度,极高度风险等级的面积分别增加了16 544.76、4 223.25 km2。区域沙漠化风险的主要影响因素为年降水量、年均风速、植被覆盖指数、沙漠化土地面积与载畜量,植被覆盖指数增加有助于区域沙漠化风险的降低。 相似文献