共查询到20条相似文献，搜索用时 281 毫秒
1.
基于时空不确定性的对流尺度集合预报效果评估检验 总被引：3，自引：0，他引：3
针对对流尺度天气系统的高度非线性特征和高分辨率模式预报结果存在时、空不确定性现象,以及当前邻域概率法主要考虑高分辨率预报结果的空间位移误差,而不能有效解决预报结果存在时间超前与滞后问题,将时间因素引入到邻域概率法中,结合一次强飑线过程进行对流尺度集合预报试验,并基于改进后的新型邻域概率法与分数技巧评分,对降水预报进行了不同时、空尺度的效果评估检验。结果表明:（1）邻域集合概率法和概率匹配平均法在极端降水的分数技巧评分远高于传统集合平均,弥补了集合平均对极端降水预报能力偏低的缺陷。（2）对于此类飑线过程的对流尺度天气系统而言,邻域半径为15—45 km的空间尺度能够改善降水位移误差的空间不确定性,并使其预报效果达到最优,其中15—30 km的邻域半径对于尺度更小的大量级降水事件预报能力更强。（3）对流尺度降水预报考虑时间尺度与降水强度存在着对应关系,不同时间尺度可以捕获到不同量级降水的时间不确定性。同时,时间尺度与空间尺度对于降水预报效果的影响是相互关联的。（4）改进的邻域概率法能够同时体现高分辨率模式预报结果在对流尺度降水事件上存在的时、空不确定性,实现了对流尺度降水在时、空尺度上的综合评估,并能为不同量级降水提供与其时、空尺度相匹配的概率预报结果。 相似文献
2.
How to accurately address model uncertainties with consideration of the rapid nonlinear error growth characteristics in a convection-allowing system is a crucial issue for performing convection-scale ensemble forecasts. In this study, a new nonlinear model perturbation technique for convective-scale ensemble forecasts is developed to consider a nonlinear representation of model errors in the Global and Regional Assimilation and Prediction Enhanced System(GRAPES)Convection-Allowing Ensemble Predi... 相似文献
3.
Summary Selected small domain LAM forecasts modulated by highly corrugated underlying topography, and driven by different state-of-science outer models suggest that uncertain outer model guidance for LAMs produces large, domain averaged sensitivity. A further literature survey indicates that many LAM forecasts are relatively insensitive to details of the local initial state, and that mesoscales show slight error growth, in contradiction to classical predictability theory. A series of global predictability experiments is presented in order to reconcile the contradiction. The experiments imply that, even in baroclinically unstable atmospheres, the most common sources of local error growth are associated with small uncertainties of the larger spatial scales rather than small uncertainties of the smaller spatial scales. Variable resolution, real-data experiments of barotropic versions of the global model display substantial mesoscale error growth, due principally to the effect of larger scales. The uncertainties possessing largest spatial scale appear as boundary uncertainties in LAMs, and explain the strong boundary sensitivity and weak local initial data sensitivity observed in many LAMs. We infer that accurate depiction of the largest spatial scales is a first order priority for accurate local prediction, and that for the advective portion of the dynamics, errors of the outer model that provides lateral boundary conditions may impose the largest current practical limitation for many LAM predictions.With 10 Figures 相似文献
4.
The limits of predictability of El Niño and the Southern Oscillation (ENSO) in coupled models are investigated based on retrospective forecasts of sea surface temperature (SST) made with the National Centers for Environmental Prediction (NCEP) coupled forecast system (CFS). The influence of initial uncertainties and model errors associated with coupled ENSO dynamics on forecast error growth are discussed. The total forecast error has maximum values in the equatorial Pacific and its growth is a strong function of season irrespective of lead time. The largest growth of systematic error of SST occurs mainly over the equatorial central and eastern Pacific and near the southeastern coast of the Americas associated with ENSO events. After subtracting the systematic error, the root-mean-square error of the retrospective forecast SST anomaly also shows a clear seasonal dependency associated with what is called spring barrier. The predictability with respect to ENSO phase shows that the phase locking of ENSO to the mean annual cycle has an influence on the seasonal dependence of skill, since the growth phase of ENSO events is more predictable than the decay phase. The overall characteristics of predictability in the coupled system are assessed by comparing the forecast error growth and the error growth between two model forecasts whose initial conditions are 1 month apart. For the ensemble mean, there is fast growth of error associated with initial uncertainties, becoming saturated within 2 months. The subsequent error growth follows the slow coupled mode related the model’s incorrect ENSO dynamics. As a result, the Lorenz curve of the ensemble mean NINO3 index does not grow, because the systematic error is identical to the same target month. In contrast, the errors of individual members grow as fast as forecast error due to the large instability of the coupled system. Because the model errors are so systematic, their influence on the forecast skill is investigated by analyzing the erroneous features in a long simulation. For the ENSO forecasts in CFS, a constant phase shift with respect to lead month is clear, using monthly forecast composite data. This feature is related to the typical ENSO behavior produced by the model that, unlike the observations, has a long life cycle with a JJA peak. Therefore, the systematic errors in the long run are reflected in the forecast skill as a major factor limiting predictability after the impact of initial uncertainties fades out. 相似文献
5.
不同海温强迫的月动力延伸集合预报试验 总被引：1，自引：0，他引：1
利用全球谱模式T106L19和增长模繁殖法（BGM）建立了月动力延伸集合预报系统,基于气候海表面温度（SST）和预测海表面温度,设计了三组集合预报试验,一组为气候SST作为模式下边界条件的集合预报试验（CSST试验）,另一组为预测SST作为模式的下边界条件的集合预报试验（FSST试验）,第三组为前两组试验的集合预报结果之和（AVE30试验）,对两种海温强迫分别进行了48个月的试验,并对预报结果进行了检验和分析。结果表明：相对于单一的控制预报,不管是CSST试验还是FSST试验,利用BGM方法制作的初值集合预报能显著提高月平均环流的预报技巧,集合预报对PNA区域的预报技巧改进显著,特别是预测SST强迫有正的贡献;同时考虑初值和边值不确定性影响的集合预报试验（AVE30试验）,其全球预报技巧不仅高于控制预报,也分别高于FSST试验和CSST试验,这说明要提高月延伸预报技巧,必须同时考虑初值和边值的影响;大气对SST强迫的响应在模式积分10天开始显著,SST对第二旬和第三旬的作用直接影响月平均环流的预报效果,而SST对第二旬和第三旬预报的影响不仅与SST本身变化有关,还与初值有关,不同的初值其作用不同;集合预报对我国夏季月平均温度分布具有较强预报能力,采用预报海温强迫的预报结果,总体上优于气候海温强迫的结果。 相似文献
6.
Mesoscale predictability of mei-yu heavy rainfall 总被引：1，自引：0，他引：1
Recently reported results indicate that small amplitude and small scale initial errors grow rapidly and subsequently contaminate short-term deterministic mesoscale forecasts. This rapid error growth is dependent on not only moist convection but also the flow regime. In this study, the mesoscale predictability and error growth of mei-yu heavy rainfall is investigated by simulating a particular precipitation event along the mei-yu front on 4-6 July 2003 in eastern China. Due to the multi-scale character of th... 相似文献
7.
The quantitative precipitation forecast (QPF) in very-short range (0-12 hours) has been investigated in this paper by using a convective-scale (3km) GRAPES_Meso model. At first, a latent heat nudging (LHN) scheme to assimilate the hourly intensified surface precipitation data was set up to enhance the initialization of GRAPES_Meso integration. And then based on the LHN scheme, a convective-scale prediction system was built up in considering the initial “triggering” uncertainties by means of multi-scale initial analysis (MSIA), such as the three-dimensional variational data assimilation (3DVAR), the traditional LHN method (VAR0LHN3), the cycling LHN method (CYCLING), the spatial filtering (SS) and the temporal filtering (DFI) LHN methods. Furthermore, the probability matching (PM) method was used to generate the QPF in very-short range by combining the precipitation forecasts of the five runs. The experiments for one month were carried out to validate the MSIA and PM method for QPF in very-short range.
The numerical simulation results showed that: (1) in comparison with the control run, the CYCLING run could generate the smaller-scale initial moist increments and was better for reducing the spin-up time and triggering the convection in a very-short time; (2) the DFI runs could generate the initial analysis fields with relatively larger-scale initial increments and trigger the weaker convections at the beginning time (0-3h) of integration, but enhance them at latter time (6-12h); (3) by combining the five runs with different convection triggering features, the PM method could significantly improve the QPF in very-short range in comparison to any single run. Therefore, the QPF with a small size of combining members proposed here is quite prospective in operation for its lower computation cost and better performance. 相似文献
8.
北京“7.21”特大暴雨不同集合预报方案的对比试验 总被引：11，自引：0，他引：11
采用6套扰动方案(初值、多物理、3组随机物理和初值与随机物理的混合)对2012年7月21日(“7.21”)北京特大暴雨过程进行了集合降水预报试验,检验了不同方案的集合平均预报、集合区间预报和概率预报较控制预报改进的相对程度,分析了它们对该过程时空不确定性的预报能力、不同扰动方法的离散度贡献以及不同尺度扰动对预报误差的贡献等。结果表明:(1)所有集合方案特别是初值扰动、多物理和混合扰动的集合预报相对控制预报在暴雨强度和位置上都有较显著的改进,并为用户决策提供了包括预报不确定性在内的更多预报信息。(2)3组随机物理产生的集合预报离散度很相似, 都远小于初值扰动和多物理方案产生的离散度, 并且主要集中在强降水中心附近, 因此在初值扰动的基础上加入随机扰动,可以提高强降水中心的离散度, 但对强降水中心以外的地区作用甚微;尺度分析进一步表明随机物理产生的离散度贡献主要集中在较小尺度上(<320 km),在更小的尺度上(<160 km)它甚至可以与初值和多物理扰动的贡献相当,而初值扰动和多物理过程的贡献可以比随机物理过程多延伸400—500 km直到较大的尺度(如>1000 km), 其中多物理过程在较小尺度上(<100 km)可比初值扰动贡献更大, 并且能部分消除预报系统偏差。(3) 所有集合扰动方案所产生的离散度尺度谱都与实际预报误差尺度谱分布一致, 即随空间尺度增大而减小,但在幅度上都小于预报误差(离散度不够大),并且这种差异随着空间尺度的减小而加速增大,在小尺度上相差甚大。 相似文献
9.
The spatial propagation of meso- and small-scale errors in a Meiyu frontal heavy rainfall event,which occurred in eastern China during 4 -6 July 2003,is investigated by using the mesoscale numerical mo... 相似文献
10.
不同大尺度强迫条件下考虑初始场与侧边界条件不确定性的对流尺度集合预报试验 总被引：5，自引：0，他引：5
研究了不同大尺度强迫条件下的暴雨个例中,考虑不同尺度特征的初始扰动与侧边界扰动相互作用构造对流尺度集合预报的可行性,为进一步构建“自适应”于不同强对流天气的对流尺度集合预报系统提供依据。结果表明,在大尺度强迫显著的个例1中,以大尺度扰动信息为主的动力降尺度的增长趋势较集合转换卡尔曼滤波(ETKF)更为显著,且总扰动能量在预报中后期超过集合转换卡尔曼滤波,而在大尺度强迫较弱的个例2中,集合转换卡尔曼滤波扰动能量始终高于动力降尺度。此外,当大尺度强迫显著时,初始扰动与侧边界扰动相匹配会产生相互促进的作用,而不匹配时初始扰动会在预报中后期抑制侧边界扰动的发展,当大尺度强迫较弱时,即使是互相间不匹配的初始扰动与侧边界扰动也能在大部分预报时段起到相互促进的作用,说明初始扰动与侧边界扰动的相互作用机制不仅与天气形势相关,也与二者是否匹配挂钩,另外,扰动的发展特征同样依赖于天气形势;从集合离散度的角度来看,当大尺度强迫明显时,侧边界扰动的作用会在更短的时间内取代初始扰动,从而对离散度起到主导地位;两种初始扰动方法相比,集合转换卡尔曼滤波在多数情况下对总离散度的贡献均大于动力降尺度;从降水量预报及概率预报情况来看,大尺度强迫明显的个例可预报性更高,且各集合成员间的差异较小,大尺度强迫较弱的个例则相反,且当两种初始扰动方案与侧边界扰动相结合时,较仅侧边界扰动均有一定提高。 相似文献
11.
The impact of applying three-dimensional variational data assimilation (3D-Var DA) on convective-scale forecasts is investigated by using two mesoscale models, the Weather Research and Forecasting model (WRF-ARW) and the Hirlam and Aladin Research Model On Non-hydrostatic-forecast Inside Europe (HARMONIE-AROME). One month (1 to 30 December 2013) of numerical experiments were conducted with these two models at 2.5 km horizontal resolution, in order to partly resolve convective phenomena, on the same domain over a mountainous area in Iran and neighboring areas. Furthermore, in order to estimate the domain specific background error statistics (BES) in convective scales, two months (1 November to 30 December 2017) of numerical experiments were carried out with both models by downscaling operational ECMWF forecasts. For setting the numerical experiments in an operational scenario, ECMWF operational forecast data were used as initial and lateral boundary conditions (ICs/LBCs). In order to examine the impact of data assimilation, the 3D-Var method in cycling mode was adopted and the forecasts were verified every 6 hours up to 36 hours for selected meteorological variables. In addition, 24 h accumulated precipitation forecasts were verified separately. Generally, the WRF and HARMONIE-AROME exhibit similar verification statistics for the selected forecast variables. The impact of DA on the numerical forecast shows some evidence of improvement in both models, and this effect decreases severely at longer lead times. Results from verifying the 24 h convective-scale precipitation forecasts from both models with and without DA suggest the superiority of the WRF model in forecasting more accurately the occurred precipitation over the simulation domain, even for the downscaling run. 相似文献
12.
The article is to report some results of numerical experiments on the error growth and the atmosphericpredictability Experiments with two-level global baroclinic primitive equation spectral model have mainresults as follows.The magnitude of initial errors directly affects the error growth,but its distributionform has little effect on the growth.The loss of predictability resulting from small-scale error is much greaterthan that from large-scale error.The small-scale error rapidly grows and is transferred to the large-scaleerror by interaction between different scale waves,which stimulates the growth of error for the whole systemOrographic forcing restrains planetary-scale error(wavenumbers 0—3)but enhances the small-scale error(wavenumbers 8 or greater).Hence,orographic effects on the error growth closely depend on the characteris-tic scale of initial errors,and there may be a critical wavenumber between 4 and 7.The error growth is great-er in Northern Hemisphere than in Southern Hemisphere if initial errors are the same.In the end we givesome discussions about model,initialization scheme,etc.,to improve model prediction. 相似文献
13.
针对2018年4月22日发生在湖北西部山地的一次特大暴雨过程,采用降尺度方案和显式对流参数化方案模式,开展了高分辨率对流许可尺度(3 km)的集合预报试验,并对全球集合预报(GEFS)和对流尺度集合预报(SSEF)的降水预报进行了对比评估,结果表明:(1)SSEF集合平均的雨量和落区预报均优于GEFS。(2)SSEF各成员的降水离散度分布更合理,因而具有更优的降水区间预报,其“离散度-误差关系”更优,能更好地给出预报误差的分布及其可能的大小。(3)SSEF的概率预报在所有空间尺度上均优于GEFS,且在短历时强降水上的优势更加明显。由此可见,针对此类山地暴雨过程,对流尺度集合预报相对于全球集合预报具有巨大的改进潜力。 相似文献
14.
2013年欧洲中心台风集合预报的检验 总被引：1，自引：0，他引：1
广州中心气象台利用中国气象局下发的欧洲中心台风集合预报数据,制作了台风集合预报产品,供业务参考应用。利用欧洲中心台风集合预报数据,对2013年1307—1331号热带气旋的集合预报路径和强度进行检验,通过对比集合平均、模式高分辨率确定性预报和预报员主观预报,发现路径集合平均在24~120 h预报误差最小;在有限的预报样本数中,从热带风暴到台风级别的热带气旋,各预报时效路径集合平均的误差随强度增强而减小;强引导气流背景下的热带气旋预报误差小于弱引导气流的误差。对比强度集合平均和模式高分辨率确定性预报,发现各时效集合平均的误差比确定性预报大,随着预报时效的延长误差没有明显增大或减小的趋势,而且强度集合平均预报,在中心最低气压、中心最大风速、热带气旋等级都表现出明显的系统性偏弱特征;对不同级别的热带气旋强度预报,集合平均的误差随强度增强而增大,即强度集合预报对强度较弱的热带气旋有更高的准确率;对比受强、弱引导气流影响的两类热带气旋,集合平均对受弱引导气流影响的一类预报误差更小。 相似文献
15.
采用2004年Casati提出的强度-尺度检验技术,选取2008年汛期代表不同类型降水(对流云降水、层状云降水、混合云降水)的4个降水过程,从尺度分解角度入手,对"世界气象组织天气研究计划——北京奥运会预报示范项目"(WWRP B08FDP)项目中4个I临近预报参加系统(BJANC,GRAPES-SWIFT,STEPS,CARDS)的1h定量降水预报进行时空尺度分解检验,研究降水预报技巧与降水时空尺度和强度之间的关系。结果表明:尽管目前国际先进的临近预报系统的水平分辨率已高达1～2km,但其有技巧的临近预报能力主要集中于空间尺度大于32km、时间尺度大于1h的降水系统,而对小于这些尺度的降水系统预报能力仍非常有限;在不同时空尺度的临近预报降水误差中,60%以上的误差来自于空间尺度小于8km的降水,85%以上的误差来自于时间尺度小于1h的降水,传统的外推技术不能满足这些较小时空尺度降水预报的需求,要发展有效的预报方法来提高较小时空尺度降水的预报能力。将基于外推的临近预报和基于稠密观测资料、快速更新的数值预报的潜势预报相结合可能是一条有效的解决途径。 相似文献
16.
Can the uncertainties of Madden-Jullian Oscillation cause a significant “Spring Predictability Barrier” for ENSO events? 
下载免费PDF全文
下载免费PDF全文 With the Zebiak-Cane model and a parameterized stochastic representation of intraseasonal forcing, the impact of the uncertainties of Madden-Jullian Oscillation (MJO) on the ??Spring Predictability Barrier (SPB)?? for El Ni?o-Southern Oscillation (ENSO) prediction is studied. The parameterized form of MJO forcing is added physically to the Zebiak-Cane model to obtain the so-called Zebiak-Cane-MJO model and then the effects of initial error, stochastic model error, and their joint error mode on the SPB associated with El Ni?o prediction are estimated. The results show that the model errors caused by stochastic MJO forcing could hardly lead to a significant SPB while initial errors can do; furthermore, the joint error mode of initial error and model error associated with the stochastic MJO forcing can also lead to a significant SPB. These demonstrate that the initial error is probably the main error source of the SPB, which may provide a theoretical foundation of data assimilation for ENSO forecasts. 相似文献
17.
X. Li 《Meteorology and Atmospheric Physics》2006,94(1-4):87-102
Summary Cloud microphysical and precipitation responses to a large-scale forcing in the tropical deep convective regime are investigated
based on hourly zonally-averaged, vertically-integrated simulation data from a two-dimensional coupled ocean-cloud resolving
atmosphere model. The model is forced by the large-scale vertical velocity and zonal wind observed and derived from TOGA COARE
for a 50-day period. The accretion of cloud water by graupel induces growth of graupel that enhances raindrops through its
melting during a weak-forcing period, whereas the large deposition rate of vapor associated with a large upper-tropospheric
upward motion causes growth of snow from the conversion of cloud ice and enhancement of graupel from the accretion of snow
during a strong-forcing period. The local changes of raindrops and graupel switch from the negative to positive values as
the forcing strengthens in the weak-forcing case, whereas the variations of cloud hydrometeors are not sensitive to the strength
of the forcing in the strong-forcing case. Phase analysis indicates that cloud water leads the surface rain rate by 1 hour.
The surface rain rate can be calculated based on the conservation of vapor and cloud hydrometeors and the budget of raindrops.
The vapor source and local changes of cloud hydrometeors could have impacts in the calculation of the surface rain rate. The
vapor source determines the surface rain rate in the strong-forcing case whereas the cloud variations could become important
in the weak-forcing case. In the budget of raindrops, the sum of the collection of cloud water by raindrops, the melting of
graupel, and the evaporation of raindrops determines the surface rain rate in the strong-forcing case whereas the other rain-related
microphysical processes become important in the weak-forcing case. 相似文献
18.
基于北京“7.21”特大暴雨个例,设计了一种考虑地形不确定性对降水影响的集合预报方案,在对该方案进行初步评估的基础上,重点通过计算相关系数、扰动总能量和尺度分解,对包括地形扰动方案在内的4种集合预报方案（初值、多物理、地形、初值-多物理混合）中离散度演变的异同性进行了分析。结果表明：（1）考虑模式地形不确定性的扰动方案,在不影响集合平均降水预报质量的基础上,对集合降水预报的离散度和概率预报略有正贡献。（2）离散度空间结构的演变与天气形势的演变密切相关。不同扰动方案产生的离散度在初始时刻的空间分布各不相同,但随模式向前积分其离散度的相似度快速增大,其中0-6 h内增长速度最快,离散度场之间的相关系数可以超过0.6。混合方案与单一扰动方案相比,对离散度空间结构的贡献不大。（3）虽然不同方案的离散度空间结构相似,但其幅度却存在明显的差异,如地形扰动方案的离散度幅度明显小于初值扰动和物理过程扰动方案。混合方案可以增加原单一扰动方案的离散度振幅,但这种增加在高层明显,而在近地面层并不明显,因而增加降水和其他近地面大气变量的离散度要比增加上层大气变量的离散度更困难。（4）尺度分离的结果表明,随着空间尺度的增大和积分时间的延长,不同扰动方法产生的离散度结构会逐渐变得相似,但在积分早期（<12 h）和较小的空间尺度（<448 km）上离散度结构的差异仍明显,并且在较小的空间尺度（<448 km）上,不同扰动方法产生的离散度幅度有明显的差异。所以对于小空间尺度系统或甚短期预报,选择扰动方案比大尺度和较长期的预报更重要。以上研究可为集合预报如何合理采用不同扰动方案或不同方案的组合提供科学依据。 相似文献
19.
本文基于气候系统模式CESM4.0长期积分试验,分析评估了工业革命前（1850年）及当前（2000年）两种辐射强迫作用（分别为太阳活动和温室气体）对中国北方干旱半干旱区降水的影响。结果表明,模拟结果与观测之间尽管存在一定的偏差,但仍能再现降水气候态的空间分布以及季节变化特征;两种辐射强迫下的降水长期变化均无明显趋势,但二者的差异却呈现出70~100年的准周期振荡;由人类活动引起的当前辐射强迫作用对降水的多年际变率幅度有一定影响,造成极端强降水事件出现的概率增多,而由太阳活动引起的辐射强迫作用主要对降水多年代际周期具有一定的调制作用。进一步分析表明,两种辐射强迫下中国北方干旱半干旱区降水年多年代际变率的主要模态基本一致,但人类活动引起的辐射强迫作用会影响降水多年代际变率与热带海温异常的相互作用的强度,从而改变降水多年代际变率的幅度。 相似文献
20.
Annalisa Cherchi Andrea F. Carril Claudio G. Menéndez Laura Zamboni 《Climate Dynamics》2014,42(1-2):219-236
An ensemble of nine experiments with the same interannually varying sea surface temperature (SST), as boundary forcing, and different initial conditions is used to investigate the role of tropical oceans in modulating precipitation variability in the region of La Plata Basin (LPB). The results from the ensemble are compared with a twentieth-century experiment performed with a coupled ocean-atmosphere model, sharing the same atmospheric component. A rotated empirical orthogonal functions analysis of South America precipitation shows that the dominant mode of variability in spring is realistically captured in both experiments. Its principal component (RPC1) correlated with global SST and atmospheric fields identifies the pattern related to El Niño Southern Oscillation and its large-scale teleconnections. Overall the pattern is well simulated in the tropical southern Pacific Ocean, mainly in the ensemble, but it is absent or too weak in other oceanic areas. The coupled model experiment shows a more realistic correlation in the subtropical South Atlantic where air-sea interactions contribute to the relationship between LPB precipitation and SST. The correspondence between model and data is much improved when the composite analysis of SST and atmospheric fields is done over the ensemble members having an RPC1 in agreement with the observations: the improvement relies on avoiding climate noise by averaging only over members that are statistically similar. Furthermore, the result suggests the presence of a high level of uncertainty due to internal atmospheric variability. The analysis of some individual years selected from the model and data RPC1 comparison reveals interesting differences among rainy springs in LPB. For example, 1982, which corresponds to a strong El Niño year, represents a clean case with a distinct wave train propagating from the central Pacific and merging with another one from the eastern tropical south Indian Ocean. The year 2003 is an example of a rainy spring in LPB not directly driven by remote SST forcing. In this case the internal variability has a dominant role, as the model is not able to reproduce the correct local precipitation pattern. 相似文献