植物气孔导度的环境响应模拟及其尺度扩展

气孔导度是衡量植物和大气间水分、能量及CO2平衡和循环的重要指标,探讨气孔导度与环境因子的关系及其模拟,以及气孔导度在叶片、冠层及区域尺度间的尺度转换及累积效应,对更好地认识植被与大气间的水热运移过程,合理评价植被在陆面过程中的地位和作用都具有重要意义。从植物气孔导度与环境因子的关系、气孔导度模拟以及尺度扩展三个方面,对前人的研究成果进行了概括总结。从叶片和冠层两个尺度出发,归纳总结了前人对于不同植物气孔导度与环境因子关系的研究成果,发现由于不同植物的遗传特性、测定时的环境、时间尺度的不同,以及未考虑各个环境因子的相互作用对气孔导度的影响,由此得到的气孔导度与环境因子之间的关系也不尽一致。对各单一环境因子与气孔导度的关系,给出了生理学解释,从根本上说明了环境因子变化对气孔导度的影响,而研究环境因子对气孔导度的综合影响时,应对各环境因子进行系统控制与同步观测。模拟计算植物气孔导度的模型主要有Jarvis模型和BWB模型两类,这些模型的模拟能力随着研究对象、试验区域、环境条件的改变而存在一定的差异,在具体使用时应结合实际情况选择最优模型进行模拟。除上述常用模型外,还总结了其他学者分别从不同角度提出的新的模型,对现有气孔导度模型进行了全面的总结。从叶片-冠层、冠层-区域两个方面归纳总结了前人关于气孔导度尺度扩展的研究成果,发现叶片-冠层的尺度扩展研究较成熟而冠层-区域的尺度扩展在模拟精度的验证方面存在困难。针对以下几个方面提出了今后气孔导度的研究重点:(1)结合研究对象所在的区域及环境条件,选择最优模型进行模拟;(2)综合考虑环境因子之间的相互作用及其对气孔导度的累积影响;(3)BWB模型与光合模型的耦合;(4)提高大尺度范围内的气孔导度模拟精度。     

芦苇湿地多时空尺度蒸散模拟研究进展

湿地蒸散是湿地水分损失的主要途径,由于其对气候变化的高度敏感性和重要的反馈作用而倍受关注。湿地蒸散过程可以分为叶片、植株、冠层、景观至区域等多种尺度,如何实现不同时空尺度的过程和参数耦合成为湿地蒸散研究的重点和难点。本文回顾了湿地蒸散模拟研究的发展历程,着重总结了芦苇湿地蒸散监测及模拟方面的技术和方法,并从不同时间尺度(日、月、年)和不同空间尺度(叶片、植株、冠层、区域)以及实现不同时空尺度拓展的关键技术研究等方面进行归纳和总结,探讨了今后的研究重点。实现不同时空尺度的芦苇湿地蒸散过程的准确模拟,关键在于大气-植被-土壤界面的参数耦合,而目前针对不同尺度的芦苇湿地蒸散参数化研究很少,因此,有必要探讨适合芦苇湿地多尺度蒸散模拟及参数化方法,为芦苇湿地蒸散评估模型的选择提供理论支持。     

黄土高原半干旱区柠条(Caragana korshinskii)树干液流动态及其影响因子

应用热脉冲技术在黄土高原神木县六道沟小流域于2006年6月13至25日测定了两种不同密度柠条(Caragana korshinskii)群落的树干液流动态.同时测量了土壤水分、太阳辐射、大气温度、相对湿度、风速、水汽压亏缺和作物参考蒸散等环境因子,并根据植物蒸腾的P-M公式,反推计算冠层导度.结果表明,除风速外,柠条树木液流与太阳辐射、大气温度、相对湿度、水汽压亏缺、作物参考蒸散均显著相关,且可用太阳辐射的线性表达式来估测.不同密度群落的日蒸腾量随叶面积指数增大而增加,叶面积指数为2.3的群落平均日蒸腾为3.83mm d-1m-2,而叶面积指数为1.1的林分平均日蒸腾1.64mm d-1m-2.冠层导度与气象因子关系复杂,当土壤水分不存在亏缺时,冠层导度与太阳辐射、大气温度、作物参考蒸散因子显著相关,与水汽亏缺和相对湿度因子无相关性;当土壤水分存在亏缺时,冠层导度与太阳辐射、大气温度、作物参考蒸散因子无相关关系,而与水汽亏缺和相对湿度因子显著相关.     

樱桃冠层导度特征及模拟

为了揭示樱桃冠层蒸腾、冠层导度对环境因子的响应规律,评价Jarvis模型在樱桃冠层尺度上应用的适用性,利用Granier热消散式探针连续监测了北京四季青果林所试验地3年生盆栽樱桃(Prunus avium L.)4-8月份蒸腾动态变化,同步监测了气象与土壤水分数据。以实测液流为基础,利用Penman-Monteith方程反推方法获取了长期连续冠层导度,在分析樱桃冠层蒸腾、冠层导度的动态变化规律的基础上,采用十字交叉法对多元回归模型与Jarvis模型进行参数率与误差分析,结果显示盆栽樱桃冠层蒸腾规律性强、时滞效应小,不同辐射条件下,冠层导度随水汽压亏缺增加呈负指数函数下降趋势,采用水汽压亏缺、光合有效辐射、气温的不同组合方式构建了多元回归和Jarvis冠层模型,模拟结果显示Jarvis模型精度高于多元回归模型,环境因子对模型精度的影响程度依次为:水汽压亏缺光合有效辐射气温,考虑了水汽压亏缺和太阳辐射的Jarvis模型精度最高,最低相对误差仅为12.12%,均方根误差为0.271。     

华北平原冬小麦冠层导度的环境响应及模拟

通过引入叶面积指数，将叶片水平的气孔导度组合模型扩展到冠层水平，建立了冠层导度环境响应组合模型，组合模型所需参数较少，且均可在冠层水平直接测量，便于应用；模型由潜在气孔导度（PSC）和相对气孔开度（RDO）组成，二者分别由环境变量的日际（inter-day）和日间（intra-day）的值决定。分析表明，冠层导度在日际尺度和日间尺度上对环境变量具有多尺度响应特性，在日际尺度上，温度是影响冠层导度的主要因子，在日间尺度，光是影响气孔开闭的主要因素。利用以温度和光合有效辐射为输入变量构建的组合模型，模拟了华北平原冬小麦农田生态系统的冠层导度，并用Penman-Monteith方程估算的表面导度进行验证。结果显示，在不同天气情况下，二者的日变化均具有较好的一致性；将组合模型与电学类比模型结合，进一步估算了拔节．灌浆期的冠层潜热通量，利用涡度相关系统观测的潜热通量数据进行验证，结果表明对冬小麦冠层潜热通量模拟精度较高，直线回归斜率为0．7054，R^2=0．7894。     

林木蒸腾耗水量测算方法的比较与应用

本文简述了林木蒸腾耗水的概念、研究价值和研究进展,从实测法和估测法的角度,综述了林木蒸腾耗水量的主要测算方法,对比了各方法的优缺点、适用性、局限性、应用现状以及适用尺度。本研究认为: 实测法可应用于多种空间尺度耗水量的测定,估测法常应用于大空间尺度蒸散量的测算。实测法是估测法的基础,因此应对实测法的测定结果进行合理的质量控制与评价,为校正估测法的估算结果提供数据基础。不管是同一大空间尺度还是不同空间尺度,实测法和估测法的结合能够提高蒸散量测算结果的准确性。提升复杂下垫面和恶劣气候环境下大空间尺度林木蒸散量的测算精度将成为未来的研究热点与难点。随着科学技术的不断进步,现有的测定装置和测算方法将被改进,精确测定林木蒸腾耗水量的新方法也将随之诞生。     

基于改进SW模型的千烟洲人工林蒸散组分拆分及其特征

蒸散组分拆分是准确评估陆地生态系统生产力以及估算水分利用效率的重要基础。利用改进后的Shuttleworth-Wallace模型,将蒸散拆分为植被蒸腾、土壤蒸发和冠层截留蒸发,并采用Monte Carlo随机参数化方案对模型参数进行优化。将模型与千烟洲亚热带人工针叶林站点的2011年涡度相关及小气候观测资料结合,对千烟洲人工林蒸散及其组分进行模拟。研究结果表明:半小时尺度上蒸散量模拟值与实测值的一致性在晴天和雨天都较高。半小时尺度上全年蒸散模拟值与实测值的决定系数、均方根误差和平均偏差为0.73、1.55 mmol m~(-2)s~(-1)和0.21 mmol m~(-2)s~(-1)。蒸散是该生态系统水分输出的最主要贡献项,占全年降水的80%。在蒸散中,植被蒸腾约占总蒸散量的85%,可推测2011年千烟洲人工林生态系统有较高的水分利用效率。该生态系统的蒸腾量季节变化明显,主要受饱和水汽压差和气温两种环境因素以及植被的叶面积指数影响且与三者均呈正相关;土壤蒸发约占总蒸散量的5%,季节变化平缓;模拟的冠层截留蒸发量约占总蒸散量的10%,季节变化大,与降水量呈正相关,与暴雨频次呈负相关,说明冠层无法有效截留强降水。该模型参数较少、时间分辨率高且可以有效模拟蒸散及其组分特征,是陆地生态系统水分循环过程研究有力的模型工具。     

华北低丘山地人工林蒸散的控制因子

人工林蒸散的影响机制研究对指导我国林业生态工程建设有重要意义。基于涡度相关技术,对华北低丘山区30年生栓皮栎-刺槐-侧柏人工混交林进行了连续2a的观测,以探讨蒸散的控制因子。结果表明:退耦系数(Ω)与冠层导度的季节动态有很好的一致性。2007、2008年的快速生长季Ω变化范围分别为0.12—0.62、0.08—0.54,平均值分别为0.37、0.23,快速生长季蒸散主要受气孔控制。2008年比2007年偏旱,该年快速生长季的退耦系数、冠层导度和蒸散低于2007年,发生干旱时蒸散主要受气孔限制。气孔对蒸散的限制作用在大气湿度较低时较高,并且这种限制作用在发生干旱时会进一步加剧。快速生长季内辐射是影响蒸散的主要环境因子,辐射资源丰富的年份气孔对蒸散的控制程度更高。     

玉米农田蒸散过程及其对气候变化的响应模拟

应用基于生理生态学过程的EALCO模型,对玉米农田生态系统的蒸散(ET)过程进行了模拟,在模型检验基础上,使用该模型模拟了玉米农田生态系统ET过程对未来气候变化的响应。结果表明,EALCO模型中能量与水过程的动态耦合机制使模型能够较好地模拟农田蒸散过程,基于涡度相关法的观测值与模型模拟值在小时、日尺度上均吻合较好,模型可以解释67%的日蒸散的变化特征。对土壤蒸发与冠层蒸腾的分别模拟显示,生长季土壤蒸发约占ET的36%。温度的升高会引起ET与冠层蒸腾的增加,同时土壤蒸发减少;ET对降水减少的响应较为敏感,主要表现在土壤蒸发的下降。大气CO2浓度升高对冠层蒸腾影响显著,该情景下冠层蒸腾下降幅度最大。研究所假设的2100年气候情景下,该农田生态系统生长季蒸散将减少,然而相对于降水的减少而言,蒸散的减少量较小,即水分支出项相对增加,因此,发生土壤水分匮乏的可能性加大,这可能会加剧该地区的暖干化趋势,给作物产量及生态环境带来威胁。     

近似贝叶斯法在光合模型参数估计中的应用

长期以来,光合作用机理模型中参数的确定都是一个难点。该文提出一种参数反演的方法,称为近似贝叶斯法(APMC),用来确定Farquhar光合模型的生理参数。通过将整个冠层抽象为一片大叶的思维抽象,笔者进一步将APMC应用到冠层尺度的生理参数求解,使直接求算冠层尺度生理参数成为可能。该文详细介绍了使用APMC估算光合模型参数的具体算法,并用实测数据进行了验证。结果表明,APMC可以很好地应用于冠层光合模型参数的估计,估计所得的参数落在参数生理上下限值之间,应用1 948个实测数据进行检验,得到决定系数0.75。模拟值和实测值的线性回归曲线斜率为1.04,与理论上的1.0非常接近。这个方法对光合模型参数的获取或许有积极的意义。     

森林水源涵养功能的多尺度内涵、过程及计量方法

近年来国内森林生态系统服务功能研究较多,但应用价值不高,其科学性受到诸多质疑。从森林水源涵养功能的多尺度内涵、过程以及其计量方法出发,对国内外研究动态和发展趋势进行了总结分析,重新审视森林水源涵养功能的研究意义,探讨森林水源涵养功能的多尺度特征。从森林水源涵养功能作用的空间尺度上看,其拦蓄洪水削减洪峰的功能仅在较小尺度上有效,而调节径流的功能只有当森林土壤的入渗量超过蒸散量时,才可能有更多地下水补给河道径流,进而增加旱季河道流量。同时森林水源涵养功能也具有明显的时间尺度特征,具体表现为在次降水事件中,由于蒸散发量较小,森林水源涵养功能的物质量等于森林不同层次的截留量,在功能上表现为拦蓄降水;在长时间尺度上,由于林地蒸散要耗去大量水分,森林水源涵养功能的物质量等于森林不同层次的截留量减去林地蒸散发量,在功能上表现为净化水质和调节径流。大多数研究仅对单一林分的个别层次蓄水功能进行研究,缺乏流域尺度或者更大空间尺度的森林水源涵养功能研究。建议从区域降尺度到流域或将坡面尺度上推到流域,集中在流域尺度解决森林水源涵养空间异质性将是解决森林水源涵养功能尺度外推的有效办法。就目前国内流行的森林水源涵养功能计量方法而言,其与尺度及研究目的有较大相关性,在研究中应根据研究目的、研究尺度和可获取的数据情况选择合适的计量方法。研究突出了不同尺度作用下森林水文过程的复杂性及重要性,并结合森林与水关系的争论问题,分析目前国内对森林水源涵养功能研究的一些误区,提出森林水源涵养功能研究的一些关键科学问题及未来可能的发展方向,主要包括:1)明确界定森林水源涵养功能的边界,探索森林水源涵养功能计量的新方法;2)加强不同时空尺度关联的森林水源涵养功能研究,包括正确评价森林水源涵养功能的时空变异规律,森林生态系统水源涵养功能的尺度效应,森林水源涵养与下游水生态安全,森林水源涵养研究范式转变等核心问题。     

利用半球图像法提取植被冠层结构特征参数

植被冠层结构深刻地影响着植物群落与环境的相互作用,对植被冠层结构的研究是深入理解植被生态系统格局、过程及其运作机制的重要基础。冠层结构特征参数的快速测量方法是植被冠层结构研究的前提,目前测量方法主要是基于实际测量的地面法,地面法一般费时费力,受人为因素影响较大,因此本文探索利用半球图像法获取植被冠层结构特征参数。通过对半球图像进行几何纠正并建立参数图层,与分类后的植被冠层图层进行运算提取植被冠层结构特征参数。将该方法应用于祁连山旺腰沟流域青海云杉冠层结构特征参数的提取,包括植被冠幅、冠层面积、冠层周长等,结果显示：半球图像法能够较好的提取植被冠层结构特征参数,该方法具有简单、客观、可重复等优点,也可作为植被冠层结构变化的监测方法。     

Water inputs across a tropical montane landscape in Veracruz,Mexico: synergistic effects of land cover,rain and fog seasonality,and interannual precipitation variability

Land‐cover change can alter the spatiotemporal distribution of water inputs to mountain ecosystems, an important control on land‐surface and land‐atmosphere hydrologic fluxes. In eastern Mexico, we examined the influence of three widespread land‐cover types, montane cloud forest, coffee agroforestry, and cleared areas, on total and net water inputs to soil. Stand structural characteristics, as well as rain, fog, stemflow, and throughfall (water that falls through the canopy) water fluxes were measured across 11 sites during wet and dry seasons from 2005 to 2008. Land‐cover type had a significant effect on annual and seasonal net throughfall (NTF <0=canopy water retention plus canopy evaporation; NTF >0=fog water deposition). Forest canopies retained and/or lost to evaporation (i.e. NTF<0) five‐ to 11‐fold more water than coffee agroforests. Moreover, stemflow was fourfold higher under coffee shade than forest trees. Precipitation seasonality and phenological patterns determined the magnitude of these land‐cover differences, as well as their implications for the hydrologic cycle. Significant negative relationships were found between NTF and tree leaf area index (R²=0.38, P<0.002), NTF and stand basal area (R²=0.664, P<0.002), and stemflow and epiphyte loading (R²=0.414, P<0.001). These findings indicate that leaf and epiphyte surface area reductions associated with forest conversion decrease canopy water retention/evaporation, thereby increasing throughfall and stemflow inputs to soil. Interannual precipitation variability also altered patterns of water redistribution across this landscape. Storms and hurricanes resulted in little difference in forest‐coffee wet season NTF, while El Niño Southern Oscillation was associated with a twofold increase in dry season rain and fog throughfall water deposition. In montane headwater regions, changes in water delivery to canopies and soils may affect infiltration, runoff, and evapotranspiration, with implications for provisioning (e.g. water supply) and regulating (e.g. flood mitigation) ecosystem services.     

Increased Drought Impacts on Temperate Rainforests from Southern South America: Results of a Process-Based,Dynamic Forest Model

Increased droughts due to regional shifts in temperature and rainfall regimes are likely to affect forests in temperate regions in the coming decades. To assess their consequences for forest dynamics, we need predictive tools that couple hydrologic processes, soil moisture dynamics and plant productivity. Here, we developed and tested a dynamic forest model that predicts the hydrologic balance of North Patagonian rainforests on Chiloé Island, in temperate South America (42°S). The model incorporates the dynamic linkages between changing rainfall regimes, soil moisture and individual tree growth. Declining rainfall, as predicted for the study area, should mean up to 50% less summer rain by year 2100. We analysed forest responses to increased drought using the model proposed focusing on changes in evapotranspiration, soil moisture and forest structure (above-ground biomass and basal area). We compared the responses of a young stand (YS, ca. 60 years-old) and an old-growth forest (OG, >500 years-old) in the same area. Based on detailed field measurements of water fluxes, the model provides a reliable account of the hydrologic balance of these evergreen, broad-leaved rainforests. We found higher evapotranspiration in OG than YS under current climate. Increasing drought predicted for this century can reduce evapotranspiration by 15% in the OG compared to current values. Drier climate will alter forest structure, leading to decreases in above ground biomass by 27% of the current value in OG. The model presented here can be used to assess the potential impacts of climate change on forest hydrology and other threats of global change on future forests such as fragmentation, introduction of exotic tree species, and changes in fire regimes. Our study expands the applicability of forest dynamics models in remote and hitherto overlooked regions of the world, such as southern temperate rainforests.     

六盘山辽东栎、少脉椴天然次生林夏季蒸散研究

2004年8～9月份,利用热扩散技术,结合微型蒸渗仪和水文学方法,研究了辽东栎、少脉椴次生林蒸散组成及其与林分结构的关系.结果表明,辽东栎和少脉椴树干的液流密度在"相对静止期"内比较稳定和微弱,其值在0.05μl·cm^-2·min^-1以下;在"活跃期"内树干液流密度上升较快,并呈单峰、双峰或多峰曲线,其值在0.25μl·cm^-2·min^-1以下;两树种单株蒸腾量有明显的种间差异,前者晴天和阴雨天单株蒸腾量分别为5.31和2.48 L·d^-1,为后者的2.3倍和3.75倍.林下灰_栒子和黄刺玫蒸腾速率日均值接近,分别为0.331和0.321 g·g^-1·h^-1.次生林日均蒸散量1.4 mm·d^-1,其中蒸腾量0.72 mm·d^-1、土壤蒸发量0.19 mm·d^-1、林冠截留量0.4 mm·d^-1,各占总量的49.6%、13.3%和37.1%.乔、灌木树种组成对次生林蒸腾量影响的表现不同,前者表现为个体蒸腾量的种间差异,而后者取决于单位林地面积上各树种的叶量.乔木层、灌木层和草本层(含土壤层)日均蒸散量分别为0.96、0.30和0.19 mm·d^-1,各占总量的65.8%、20.9%和13.3%,说明乔木层对林分日蒸散量大小起主要作用,灌木层次之,草本和土壤蒸发量的贡献最小.     

森林蒸散模型参数的确定

以长白山阔叶红松林为研究对象,利用长白山阔叶红松林气象观测塔安装的常规气象梯度观测系统、开路涡动相关系统的观测数据,依据空气动力学基本理论与能量平衡方程建立了森林蒸散机理模型,确定了模型参数,即风速廓线稳定度订正函数φ_m、温度廓线稳定度订正函数φ_h和零平面位移高度d.其中,研究地的零平面位移高度d为17.8m,是平均冠层高度(26m)的0.68.同时,给出了φ_m和φ_h随梯度理查逊数R_i变化的数学表达式.     

森林蒸散模型参数的确定

以长白山阔叶红松林为研究对象,利用长白山阔叶红松林气象观测塔安装的常规气象梯度观测系统、开路涡动相关系统的观测数据,依据空气动力学基本理论与能量平衡方程建立了森林蒸散机理模型,确定了模型参数,即风速廓线稳定度订正函数ψm、温度廓线稳定度订正函数ψh和零平面位移高度d。其中,研究地的零平面位移高度d为17．8m,是平均冠层高度(26m)的0．68。同时,给出了ψm和ψh随梯度理查逊数Ri变化的数学表达式。     

Do the energy fluxes and surface conductance of boreal coniferous forests in Europe scale with leaf area?

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (G_s), water‐ and light‐use efficiency and surface–atmosphere coupling of European boreal coniferous forests was explored using eddy‐covariance (EC) energy and CO₂ fluxes. The observed scaling relations were then explained using a biophysical multilayer soil–vegetation–atmosphere transfer model as well as by a bulk G_s representation. The LAI variations significantly alter radiation regime, within‐canopy microclimate, sink/source distributions of CO₂, H₂O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem‐scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry‐canopy evapotranspiration (ET) was reasonably ‘conservative’ over the studied LAI range 0.5–7.0 m² m^?2. Both ET and G_s experienced a minimum in the LAI range 1–2 m² m^?2 caused by opposing nonproportional response of stomatally controlled transpiration and ‘free’ forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m² m^?2) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m² m^?2). This finding emphasizes the significance of stand‐replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light‐saturated water‐use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests.     

用光合-蒸散耦合模型模拟冬小麦CO2通量的日变化

根据SPAC理论建立了一个冬小麦光合和蒸散的耦合模型.冬小麦CO2通量包括冠层光合、呼吸和土壤呼吸.冠层光合采用了Farquhar光合作用生化模型,并通过冠层阻力的参数化将光合作用与蒸腾作用耦合起来.用涡度相关方法观测了CO2通量,对模型进行了验证,结果显示模型可以较好地模拟CO2通量日变化过程.对模型的敏感性分析发现日间CO2通量最敏感的参数是初始量子效率.其次,CO2通量对光响应曲线凸度、CO2补偿点、凋萎点和叶面积指数的变化也有着较强的敏感性;夜间CO2通量敏感的参数是最适温度下Rubisco催化能力和暗呼吸参数.     

用光合-蒸散耦合模型模拟冬小麦CO2通量的日变化

根据 SPAC理论建立了一个冬小麦光合和蒸散的耦合模型。冬小麦 CO2 通量包括冠层光合、呼吸和土壤呼吸。冠层光合采用了 Farquhar光合作用生化模型 ,并通过冠层阻力的参数化将光合作用与蒸腾作用耦合起来。用涡度相关方法观测了 CO2通量 ,对模型进行了验证 ,结果显示模型可以较好地模拟 CO2 通量日变化过程。对模型的敏感性分析发现日间 CO2 通量最敏感的参数是初始量子效率。其次 ,CO2 通量对光响应曲线凸度、CO2 补偿点、凋萎点和叶面积指数的变化也有着较强的敏感性 ;夜间 CO2 通量敏感的参数是最适温度下 Rubisco催化能力和暗呼吸参数