A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

A ground hydrologic model with inclusion of a layer of vegetation canopy that can interact with general-circulation model

In this paper, a ground hydrologic model(GHM) is presented in which the vapor, heat and momentum exchanges between ground surface covers (including vegetation canopy) and atmosphere is described more realistically. The model is used to simulate three sets of field data and results from the numerical simulation agree with the field data well. GHM has been tested using input data generated by general circulation model (GCM) runs for both the North American regions and the Chinese regions, The results from GHM are quite different from those of GHMs in GCMs. It shows that a more active concerted effort on the land surface process study to provide a physically realistic GHM for predicting the exchange between land and atmosphere is important and necessary.

Radiation-use efficiency and the harvest index of winter wheat at different nitrogen levels and their relationships to canopy spectral reflectance

Radiation-use efficiency (RUE, g/MJ) and the harvest index (HI, unitless) are two helpful characteristics in interpreting crop response to environmental and climatic changes. They are also increasingly important for accurate crop yield simulation, but they are affected by various environmental factors. In this study, the RUE and HI of winter wheat and their relationships to canopy spectral reflectance were investigated based on the massive field measurements of five nitrogen (N) treatments. Crop production can be separated into light interception and RUE. The results indicated that during a long period of slow growth from emergence to regreening, the effect of N on crop production mainly showed up in an increased light interception by the canopy. During the period of rapid growth from regreening to maturity, it was present in both light interception and RUE. The temporal variations of RUEAPAR (aboveground biomass produced per unit of photosynthetically active radiation absorbed by the canopy) during the period from regreening to maturity had different patterns corresponding to the N deficiency, N adequacy and N-excess conditions. Moreover, significant relationships were found between the RUEAPAR and the accumulative normalised difference vegetation index (NDVI) in the integrated season (R-2 = 0.68), between the HI and the accumulative NDVI after anthesis (R-2 = 0.89), and between the RUEgrain (ratio of grain yield to the total amount of photosynthetically active radiation absorbed by the canopy) and the accumulative NDVI of the whole season (R-2 = 0.89) and that after anthesis (R-2 = 0.94). It suggested that canopy spectral reflectance has the potential to reveal the spatial information of the RUEAPAR, HI and RUEgrain. It is hoped that this information will be useful in improving the accuracy of crop yield simulation in large areas.

沙坡头人工植被区陆气耦合模式及生物结皮与植被演变的机理研究

论述了研究陆气相互作用的意义和现状。在以前工作基础上，针对腾格里沙漠人工植被区陆气水热传输过程，提出了一个多层陆气的耦合模式。特别给出了导水率的计算模型和修正后的根系吸水模式，考虑了结皮层对于土壤水分入渗的影响以及植被演变的过程。陆气耦合模式分别对大气、植被、土壤作用多层划分，以助于细致了解沿高度分布的各物理量。同时，介绍了当地气候根况和野外观测情况。利用本模式对中国科学院沙坡头沙漠站人工植被区陆气水热效换过程进行了数值模拟，模拟结果与实测值吻合较好，可为当地合理利用水热资源和治理沙漠提供科学依据。

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

东灵山森林生态系统若干水分问题的研究

东灵山地区年均降水量659.7mm，单次降水以雨量小、雨强低的降水为主。水汽压（年均17.7mb）、相对湿度（年均66％）的季节变化呈现生长季高、冬季低的趋势。年均蒸发量1019.5mm;气温、风速、日照时间和水压与月蒸发量和日蒸量相关显著;气温、日照时间和水压分别在11-6月、7－8月和9－10月为决定蒸发量的首要因子。枯枝落叶层、土壤层湿度主要受前十日降水量和坡向影响。 植物体含水量生长季节较高，冬季较低;含水量随径级的增大而降低。六个灌木树种的平均含水量大小顺序为：毛榛（48.62%）最高荆条（36.32%）最低;七个乔木树种水分含量为油松，56.14%;蒙椴，54.19%;华北落叶松，52.91%;五角枫，43.64%;辽东栎，41.87%;棘皮桦，41.13%;大叶白腊，37.79%。几种植被类型的储水量为：辽东栎林，82.08mm;华北落叶松林，47.35mm;混交林，34.60mm;油松林，31.33mm;灌丛，12.40mm。各树种叶片日最低水势的季节均值为：辽东栎，-16.1bar;五角枫，-15.8bar;大叶白腊，-15.1bar;糠椴，-13.4bar;棘皮桦，-12.3bar;蒙椴，-12.2bar。叶片水势的日间变化均呈一“V”形曲线;光照在叶片水势的日间变化中起着决定性作用。 96年各树种平均单株树干茎流量为辽东栎，30.3mm(4.19%);华北落叶松，16.1mm(2.22%);油松，8.9mm(1.23%);棘皮桦，2.9mm(0.40%)。两个生长季各林分冠层的水量平衡为：辽东栎林，树干流茎量101.87mm(9.18%)，穿透降水量823.08mm(74.15%)，截留量185.05mm(16.67%);华北落叶松林，树干径流量66.88mm(6.03%)，穿透降水量836.92mm(75.40%)，截留量206.20mm(18.58);混交林，树干径流量50.13(4.52%)，穿透降水量846.78mm(76.29%)，截留量212.20mm(19.12%);油松林，树干径流量33.90mm(3.05%)，穿透降水量934.88mm(84.22%)，截留量141.22mm(12.72%)。多元回归分析表明，树干流茎量S与降水量P和前24小时降水量P_1呈显著正相关关系;穿透降水量T与降水量P和最大雨强M正相关显著。附加截留量与降水时间成正比。 枯枝落叶层的生物量为：油松林，25.56t/hm~2;华北落叶松林20.01t/hm~2;辽东栎林，8.31t/hm~2;混交林，7.98t/hm~2。枯枝落叶层的平均实际持水量和有效持水量均以油松林最大，其次是华北落叶松林，而混交林和辽东栎林较低;枯枝落叶层的实际持水量和有效持水量的季节变化分别与前十日降水量P10成正相关和负相关关系。枯枝落叶层的截留量为油松林＞华北落叶松林＞辽东栎林＞混交林;油松林（145.632mm和90.800mm）混交林(61.816mm和54.504mm)。油松林、辽东栎林、混交林和华北落叶松林去除枯枝落叶层后，土壤入渗量比对照平均降低100mm以上;表层土壤含水量分别比对照土壤下降了6.26、18.26、15.06和15.07个百分点。地表径流量分别增加了，辽东栎林34.299mm(603%)和15.816mm(525%);油松林14.593mm(732%)和10.584mm(1321%);混交林12.004mm(181%)和7.275mm(364%);华北落叶松林3.555mm(118%)，3.275mm(229%)。96年生长季，各土壤流失量分别增加了：油松林172.751t/hm~2(124倍);辽东栎林836.500t/hm~2(119倍);混交林172.499t/hm~2(47倍);华北落叶松林11.557t/hm~2(11倍)。表层土壤容重分别增加了：油松林15.0%和20.6%，辽东栎林18.4%和28.2%，混交林11.5%和38.5%，华北落叶松林4.3%和17.1%。 0-60cm深度土壤容重平均值的大小顺序为：草地>灌丛>辽东栎林>油松林>混交林>华北落叶松林;而土壤孔隙度的大小顺序为华北落叶松林>混交林>油松林>辽东栎林>灌丛>草地。两个生长季为土壤实际储水量的均值：油松林，124.45mm，78.62mm;辽东栎林，131.23mm，87.72mm;混交林，180.41mm，113.90mm;华北落叶松林，165.53mm，127.95mm;灌丛，172.50mm，89.81mm;草地，152.92mm，89.59 mm分别比干旱年份97年高出45.83mm、43.51mm、51.63mm、37.58mm、82.69mm和63.33mm。两个生长季的地表径流量为草地，30.930mm(2.79%);灌丛，16.321mm(147%);油松林，2.911mm(0.26%);辽东栎林，8.703mm(0.78%);混交林，8.625mm(0.78%);华北落叶松林，4.447mm(0.40%)。油松林、混交林和华北落叶松林地表径流量与降水量P(mm)和最大雨强(mm/h)正相关显著;而辽东栎林、灌丛和草地的地表径流量则与降水量P(mm)、平均雨强Q(mm/hr)和最大雨强M(mm/hr)三者之间呈显著正相关关系。与草地相比(1220.093kg/hm~2，100％)，灌丛、辽东栎林、混交林、油松林和华北落叶松林96年生长季的土壤流失量分别降低了85.05%、94.26%、96.99%、98.86和99.14%。 降水量是影响小流域径流量时间变化的主要因素;南沟和马牙石沟96年的径流量分别是97年的8.19倍和7.87倍，而径流深(46.25mm，52.75mm)分别比97年(5.65mm,6.70mm)高出40.60mm和46.05mm。两个小流域由于面积的差异而使南沟两年的径流量分别比马牙石沟高出2773.136m~3(13.15%)和235.434m~3(8.79%)。96年和97年马牙石沟径流深比南沟高出6.5mm(14.05%)和1.05mm(18.58%)。在地处大陆性季风气候区的东灵山地区，用0.010m~3/min/km~2/hr能较好地分割小流域的洪峰和基流。在五次暴雨水文曲线中，马牙石沟的快速径流量分别比南沟高出25.00%到143.33%。五次洪水水文响应R的平均值南沟为0.218%，马牙石沟为0.404%;与海洋性气候地区相比，东灵山地区小流域的R值要低一到两个数量级。马牙石沟洪峰流量Qp的平均值为418.772L/min要比南沟(281.191L/min)大48.9%。东灵山地区小流域的洪水径流过程可分为三种类型。

玉米高产优质生产中若干生理生态问题的研究

玉米（Zea mays L.）是我国十分重要粮食、饲料和工业原料作物，种植区域覆盖我国大部分农业区。随着玉米品种改良和新栽培技术的应用，我国玉米产量大幅度增加。自1950s以来，我国玉米产量递增幅度为126kg/hm2/yr。在玉米产量提高过程中，单叶光合作用与产量之间存在什么样的关系？当代玉米品种的品质和养分利用效率如何？高密度种植条件下是否存在“根系拥挤”及如何调控等。为探讨上述科学问题，本研究选择中国北方常见的大田玉米品种，在高肥力自然光照条件下，探讨玉米高产优质栽培过程中生理生态特征的变化趋势，以指导科学育种和栽培。主要研究结果如下： 　　1）光合与产量的演变我国 1950s、1970s、1990s等不同年代推广的玉米品种中，当代品种叶片光合速率高且高值持续期长，光合色素叶绿素a、叶绿素b、类胡萝卜素等的含量高且持续时间长，与光合有关的蒸腾速率（E.）、细胞间隙CO2浓度（Ci.）、气孔导度（gs）等也有较大改良，中下部叶片尤其明显;在生育后期，当代品种具有更高的光合优势。老品种饱和光合速率(Psat)在灌浆期下降，并非RuBPCase 和PEPCase的活性降低，而是由于叶绿素含量和可溶性蛋白含量的降低。在花后期间，由于PS2功能的下降，造成了光合能力下降，而现代品种的PS2 功能在衰老前一致保持旺盛状态。 　　老品种光合特征对缺氮的反应表现更敏感。花后缺氮光合作用下降是非气孔限制的，因为气孔导度和胞间CO2浓度没有发生明显的变化。其主要原因是缺素造成老品种叶片早衰，叶绿素含量、可溶性蛋白含量、PEP羧化酶活性下降。现代品种表现较强的抗衰老能力，其N素利用效用高于老品种。我国玉米产量的大幅度提高在很大程度上应归功于叶片光合性能的改良。 　　随玉米品种更替，群体光合速率增强，群体光合衰减率降低，呼吸消耗所占总光合的百分率下降。灌浆期当代品种中下部叶片的群体光合速率明显高于老品种。种植密度是影响玉米群体光合速率的主要因素，在高中低三种密度条件下，当代品种均有较高的群体光合速率，表现出耐密性强、适应性广、源足库大、产量高的特点。 　　2）高油玉米的产量受到叶源大小和叶源活力的双重限制在 1.5 株/m2密度下，与普通玉米相比较，高油玉米单株籽粒产量显著低于普通玉米，产量构成中穗粒数差异不显著，千粒重较低（P＜0.01）;两类型玉米的单株库容量相当，高油玉米籽粒灌浆速率小，籽粒充实度低，单粒重对叶源相对减少（剪叶）或相对增多（疏库）的反应比普通玉米更为敏感，其产量受到同化产物供应（叶源）相对不足的限制。高油玉米授粉后的叶面积、叶面积持续期小，叶片含氮量和光合速率较低，说明高油玉米的产量受到叶源活力（光合速率）小和叶源数量少的双重限制。 　　3）我国北方玉米品种的个体产量潜力、氮素利用效率及籽粒与秸秆粗蛋白质含量在充分发挥个体生产潜力的低密度条件下，我国北方1990s 以来大面积种植的50个玉米主栽品种中，个体产量潜力和氮素利用效率高度正相关(P=0.01)，而子粒千粒重与NUE 呈显著性负相关(P=0.002)。对玉米产量和氮素利用效率进行分层聚类，可将北方玉米品种划分为高产高NUE 型、低产低NUE 型和中间型，高产高NUE 型玉米品种相对较少，仅占24%。籽粒粗蛋白质含量（CPC）与秸秆CPC 相关性不显著（P＞0.05）。对籽粒和秸秆的CPC 进行分层聚类，将北方玉米品种划分为籽粒高秸秆低型、籽粒与秸秆双低型和籽粒与秸秆双高型，CPC 双高型品种相对较少，仅占20%。 　　4）玉米根系拥挤效应对产量影响的生理生态机制及其调控随玉米品种更替根系的空间分布呈“横向紧缩，纵向延伸”的特点。当代三类型玉米根系分布特性与株型、穗型相关。紧凑型品种根系分布深，下层根系所占比率大，适合密植，群体产量潜力大;平展大穗型品种根量多，分布较浅，在低密度下可获得较高的个体生产力，但不适合密植，群体产量潜力小。 　　“根系拥挤”显著影响玉米产量，减小根系横向伸展空间，下层土壤中的根系分配比率增多。在地上部充分生长条件下，紧凑型品种横向空间为30-50cm即可满足要求，平展型品种大于50cm;紧凑型品种对纵向空间受限制的反应更为敏感，平展型品种对横向空间受限制的反应更为敏感。“根系拥挤”影响根系活性、分布、氮素吸收利用和花后光合与14C同化物的分配。 　　在根系受限制条件下，增施肥料产量提高，根系总重增加，增加了根系在深层土壤（60-100cm）中的根系比率，显著增加了根系的TTC 还原量、SOD、CAT、POD活性。土壤加沙，根量减少，但根系TTC 还原量增加、产量提高，提高幅度以大穗型品种更为显著。 　　随种植密度增加耕层根系密度与群体产量同步增大，各类品种均在最高根系密度下获得最高产量。根系负荷的籽粒产量潜力三类型品种存在极大差异，在一定范围内增大种植密度，根系伸展空间减小，群体产量提高，紧凑大穗型品种产量最高，品种的耐密性是限制根系负荷籽粒产量潜力的主导因素。因此，培育株型紧凑、耐密性强、大穗玉米良种，采取有效的调控措施是玉米进一步高产的主攻方向。 　　5）我国夏玉米高产田的培创理论研究与实践相结合，2005 年在我国华北地区的山东莱州培创出籽粒实产21 042.9kg/hm2 （ 14% 含水量， 实收面积=45.7m×15.9m=726.63m2）的夏玉米高产纪录。主要采用以增加密度为保障的“群体结构性挖潜”和以提高整齐度为保障的“个体功能性挖潜”途径，生理生态指标包括：选用紧凑抗倒耐密植品种DH3719，种植密度102 030 株/hm2，收获密度98 610 株/hm2，花后具有较长的叶面积高值持续期，达60d以上，叶面积指数最大为6.53，收获2.59。上部叶片光合值对外界光强度变化敏感，其光合峰值出现时间提前，而后迅速衰减;中部叶片光合值的降低较慢，下部叶片变幅最小，可能是长期处于争光环境表现出的生态适应性。粒叶比0.32，经济系数0.542，单株产量216g，千粒重375.1g。