大气CO2浓度升高及气候变化对作物冠层光合影响的数值模拟

利用美国Licor-6200光合作用测定仪,对黄淮海地区代表性冬小麦品种鲁麦23号叶片光合作用速率进行了较为全面的测定,分别确定了冬小麦叶片光-光合作用响应曲线和CO2-光合作用响应曲线,在此基础上,建立了叶片光合作用模式,并进而建立了一个具有瞬时时间尺度,空间积分为叶片尺度的冬小麦冠层模式,利用模式分别分析了大气中CO2浓度升高和温度变化对冠层光合作用的不同影响,并在此基础上进一步进行了综合数值分析.单因子分析表明晴天状况下,冠层光合速率随CO2浓度升高而上升,当CO2浓度由330×10-6上升至660×10-6时,冠层光合日总量可增加19.7 %;冠层光合速率随辐射增加而增大,辐射量增加10.0 %,冠层光合日总量可增加6.7 %;冠层光合速率随温度升高而下降,温度升高1 ℃,冠层光合日总量减少2.9 %.多因子综合数值分析表明在辐射量较大的气候背景下,冠层光合日总量对温度和CO2变化响应更加敏感.本文的实测数据为研究气候变化对中国农业影响提供了最基本的可靠模型参数,冠层光合模型为未来改进作物模型提供了理论基础.     

O3,CO2浓度及太阳光谱变化对作物光合作用影响的数值模拟研究

对光合作用 蒸腾作用 气孔调节进行耦合 ,从生物化学尺度扩展至冠层尺度 ,发展了一个冬小麦冠层光合作用生态动力模式 ,模式考虑了O3,CO2 和光谱变化对作物光合的综合影响。利用美国光合作用实测资料对模式进行验证 ,叶片模式通过了相关显著性检验并具有较高的准确度。数值分析表明 :当O3 浓度由 0× 10 -9V/V上升至2 0 0× 10 -9V/V时 ,冠层光合速率下降 2 9%左右 ;当CO2 浓度由 330× 10 -6V/V上升至 6 6 0× 10 -6V/V时 ,冠层光合速率增加大约 37% ;当光谱比例系数由目前的 0 .5下降至 0 .4时 ,冠层光合速率将下降 2 7%左右。对于污染严重、易发生光化学烟雾的城郊附近 ,在阳光强烈的典型晴天 ,中午O3 浓度达到 2 0 0× 10 -9V/V时 ,即使气候条件不发生改变 ,CO2 浓度对作物光合作用的正效应也不足以弥补O3 浓度升高所造成的负效应 ,冠层光合速率将比目前干洁地区略有下降 ,如果进一步考虑光合作用有效辐射光谱成分下降至 0 .4左右 ,冠层光合作用将比目前的BASE值下降 35 %左右。     

大气CO2浓度升高及气候变化对作物冠层光合影响的数值模拟

利用美国Licor-6200光合作用测定仪,对黄淮海地区代表性冬小麦品种鲁麦23号叶片光合作用速率进行了较为全面的测定,分别确定了冬小麦叶片光—光合作用响应曲线和CO₂—光合作用响应曲线,在此基础上,建立了叶片光合作用模式,并进而建立了一个具有瞬时时间尺度,空间积分为叶片尺度的冬小麦冠层模式,利用模式分别分析了大气中CO₂浓度升高和温度变化对冠层光合作用的不同影响,并在此基础上进一步进行了综合数值分析。单因子分析表明:晴天状况下,冠层光合速率随CO₂浓度升高而上升,当CO₂浓度由330×10^-6上升至660×10^-6时,冠层光合日总量可增加19.7%;冠层光合速率随辐射增加而增大,辐射量增加10.0%,冠层光合日总量可增加6.7%;冠层光合速率随温度升高而下降,温度升高1℃,冠层光合日总量减少2.9%。多因子综合数值分析表明:在辐射量较大的气候背景下,冠层光合日总量对温度和CO₂变化响应更加敏感。本文的实测数据为研究气候变化对中国农业影响提供了最基本的可靠模型参数,冠层光合模型为未来改进作物模型提供了理论基础。     

冠层内动量和物质吸收特征的数值研究

本文设计了一冠层(CL)和大气边界层(ABL)之间物质和动量交换的耦合模式,并对CL内风速、物质随高度分布和日变化作了数值模拟.结果表明,由大尺度扩散引起冠层低层的第二个风速极大和多极值的浓度分布,CL内湍流通量和物质浓度随高度减小而迅速降低,以及CL动量减小对浓度分布的重要影响,模式都能很好地描述,模拟结果与观测事实有好的一致性.利用浓度和温度廓线相似假设,导出了质量汇的经验关系.     

农田风状况及其模式

本文讨论了农田中风速的水平和铅直分布的基本规律性。指出农田植被中的风速廓线具有上下不对称的“S”型分布。分析了农田植被中风速水平和铅直分布廓线常见模式的缺点,在理论分析和经验拟合的基础上提出了新的描述农田植被中风速的通用模式。按此模式计算的平均风速,无论水平分布或铅直廓线都与实际情况比较一致。     

用K理论研究南京大厂地区污染物扩散

本文用三维K理论扩散方程对南京大厂地区污染物浓度分布特征进行了讨论。在K理论扩散方程中,输入的风速廓线、扩散系数廓线采用Businger通量廓线关系式。用塔层风资料求取了南京地区近地面特征量。用具有四阶精度的Chapeau函数方法对三维K理论扩散方程进行数值求解,计算了大厂地区SO_2年平均浓度分布,和实测值基本上吻合。本文还研究了忽略水平扩散系数、不考虑风速随高度变化,对污染物地面浓度分布的影响。结果表明:在点源的下风方1—5千米范围内,不管是较强的还是较弱的风,忽略水平扩散系数都会使计算的地面浓度估计过高;忽略风速随高度的变化,会给计算带来一定的误差。     

包含CO2因子的冠层光合作用简化模型及其关键参数数值敏感分析

利用美国Licor公司生产的Li-6200便携式光合作用测定仪对冬小麦叶片光合作用进行了较为系统的测定研究,在实测资料支持下建立了包含CO2因子的叶片和冠层光合作用模型,验证结果表明,模型具有较高准确性.利用模型对CO2影响参数μ进行数值敏感分析,结果表明:对叶片光合而言,μ=0.83时叶片光合作用模型拟合最好,随μ减小或增大,拟合精度均下降;对冠层光合而言,μ=0.83时模拟光合日总量达到最大值,μ减小或增大,光合日总量均下降.多因子数值分析表明:在辐射值较大的状况下,μ对冠层光合的影响更为显著.     

黑云杉林大气CO2稳定同位素组成δ13 C,δ18 O的分析

研究CO2稳定同位素特征可以揭示光合、呼吸作用等众多信息，从而有助于了解生态系统与环境之间的碳循环过程。利用大气CO2浓度及其稳定同位素的测定资料，分析稳定同位素比δ^13Cδ^18O，发现两者具有相似的时空分布特征。主要表现在δ^13C和δ^18O在冠层内具有明显的垂直变化趋势，冠层上部重同位素含量较高，而底部含量较低。从时间变化看，δ^13Cδ^18O在午后到日落具有较高的水平，而凌晨，δ^13Cδ^18O较低。运用Keeling图法分析δ^13Cδ^18O和CO2浓度的关系，发现δ^13C，值主要受光合和呼吸作用影响，与CO2浓度的高低密切相关；而δ^18O的变化比较复杂，因为δ^18O除受光合、呼吸强度影响外，还取决于环境空气湿度。     

多层城市冠层模式的建立及数值试验研究

为在城市气象数值模拟中更好地体现由城市发展引起的下垫面土地利用改变及人为活动对大气过程的影响,建立了基于建筑物三维分布的多层城市冠层模式,冠层内动力方程组考虑了建筑物冠层拖曳力的作用及雷诺应力的影响,通过引入建筑物宽度、间距以及垂直分布密度指数等建筑物形态特征参数,以更好地体现城市复杂地表对大气温度、湿度及动量方程的影响.同时,该模式分屋顶、4个侧壤及地面分别考虑辐射及能量平衡求解表面温度,计算各表面与大气的通量交换,并考虑辐射阴影效应、冠层内部各个面之间的可视因子、以及与冠层内建筑物密度指数、可视因子等相关的多重反射辐射导致的辐射截陷作用.模式的离线检验结果表明:(1)冠层模式计算风廓线与风洞实验测量数据吻合良好;(2)离线冠层模式能够模拟实际小区的风速、温度垂直廓线,并能够较好地体现小区内气温日变化.冠层模式与区域边界层模式耦合检验结果表明:(1)耦合模拟的近地面(2 m处)气温及地表温度的结果明显优于传统的水泥平板方案,尤其是在夜间,水泥平板方案与实测气温最大偏差4 K左右,耦合模拟方案为1-2 K;(2)耦合模拟方案考虑了建筑物对冠层之上的拖曳力影响以及建筑物形态结构对雷诺应力的影响,风速(10 m处)计算结果与观测值相差约在1 m/s,水泥平板方案偏差3 m/s左右.     

丘陵浅山地区风速幂指数及其在大气环境质量预测中的应用

根据丘陵浅山地区8个测点低空风速实测资料进行风速廓线拟合，求得风速幂指数，将大气环境质量预测中的风速模式化，用来预测大气污染物的浓度分布。     

城市街道峡谷汽车尾气污染的数值模拟

建立了一个简单的三维街道峡谷空气污染模式,并用实测资料进行了验证。利用建立的模式设计了7种试验方案,对街道峡谷内的污染状况进行了模拟。以CO为模拟对象的数值试验结果表明,街道峡谷上空的风速风向条件是决定街道峡谷内的污染状况的重要因素。峡谷上空风向与街道轴线的夹角越大、风速越小,则街道地面CO浓度越高。以现有的兰州典型车流量和排放因子,兰州街道地面CO浓度容易超标;若不控制车流量,到2008年,即使兰州上路汽车排放达标,但街道地面CO浓度仍然容易超标。     

An Analytical Model for Mean Wind Profiles in Sparse Canopies

Existing analytical models for mean wind profiles within canopies are applicable only in dense canopy scenarios, where all momentum is absorbed by canopy elements and, hence, the effect of the ground on turbulent mixing is not important. Here, we propose a new analytical model that can simulate mean wind profiles within sparse canopies under neutral conditions. The model adopts a linearized canopy-drag parametrization and a first-order turbulence closure scheme taking into account the effects of both the ground and canopy elements on turbulent mixing. The resulting wind profile within a sparser canopy appears to be more like a logarithmic form, with the no-slip condition at the ground being satisfied. The analytical solution converges exactly to the standard surface-layer logarithmic wind profile in the case of zero canopy density (i.e., no-canopy scenario) and tends to be an exponential wind profile for a dense canopy; this feature is unique compared with existing analytical models for canopy wind profiles. Results from the new model are in good agreement with those from laboratory experiments and numerical simulations.     

Line-source carbon dioxide release

Model predictions of CO₂ concentrations downwind from a line source were calibrated using experimental data. Agreement between the model and experimental data was improved by adjusting for wind direction meander and cup anemometer overshoot. The model predictions showed that by using a negative exponential wind speed profile within the crop canopy, predictions were closer to observed CO₂ concentration profiles than when experimentally-observed wind speed profiles, which were constant with height in the lower canopy, were used. This finding suggests that much of the lower canopy airflow was not direct mass flow in the downwind direction. Eddy diffusivity profiles which showed a within-canopy local minimum resulted in arestriction in the predicted loss of CO₂ out of the canopy system. Two-dimensional plots of predicted null vertical flux and CO₂ concentration portrayed vividly the turbulent diffusion and mass flow transport of CO₂ from the line source.     

Are Urban-Canopy Velocity Profiles Exponential?

Using analyses of data from extant direct numerical simulations and large-eddy simulations of boundary-layer and channel flows over and within urban-type canopies, sectional drag forces, Reynolds and dispersive shear stresses are examined for a range of roughness densities. Using the spatially-averaged mean velocity profiles these quantities allow deduction of the canopy mixing length and sectional drag coefficient. It is shown that the common assumptions about the behaviour of these quantities, needed to produce an analytical model for the canopy velocity profile, are usually invalid, in contrast to what is found in typical vegetative (e.g. forest) canopies. The consequence is that an exponential shape of the spatially-averaged mean velocity profile within the canopy cannot normally be expected, as indeed the data demonstrate. Nonetheless, recent canopy models that allow prediction of the roughness length appropriate for the inertial layer’s logarithmic profile above the canopy do not seem to depend crucially on their (invalid) assumption of an exponential profile within the canopy.     

An Inverse Method for Determining Source Characteristics for Emergency Response Applications

Turbulence above and within canopies has characteristics distinct from that over rough surfaces. The vertical transport of momentum and scalars is dominated by coherent structures whose origin is now thought to be the result of the unstable inflexion in the profile of the mean wind speed established by the application of canopy drag. This distinctive property leads to the failure of the standard Monin–Obukhov flux–profile relationships over homogeneous canopies, relationships that are assumed in many surface exchange schemes within numerical weather prediction and general circulation models. A modification of the flux–profile relationships is presented that incorporates the effects of the canopy turbulence. The subsequent impacts on the evolution of the surface energy balance and boundary-layer state are investigated within a simple numerical model for the evolution of the boundary layer and canopy state. By comparing cases with and without the modification it is shown that canopy-generated turbulence can lead, not only to the alteration of the flux–profile relationships above the canopy, but also to a different evolution of the surface energy balance and differences in near-surface conditions that would be significant in numerical weather prediction. More fundamentally, the modifications to the flux–profile relationships imply that parameters such as the roughness length and displacement height for canopies should not be considered as invariant properties, but rather as properties that depend on the flow and hence vary systematically with the diabatic stability of the boundary layer.     

Numerical Simulation Experiment of Land Surface Physical Processes and Local Climate Effect in Forest Underlying Surface

Based on the basic principles of atmospheric boundary layer and plant canopy micrometeorology, a forest underlying surface land surface physical process model and a two-dimensional atmospheric boundary layer numerical model are developed and numerical simulation experiments of biosphere and physiological processes of vegetation and soil volumetric water content have been done on land surface processes with local climate effect. The numerical simulation results are in good agreement with realistic observations, which can be used to obtain reasonable simulations for diurnal variations of canopy temperature, air temperature in canopy, ground surface temperature, and temporal and spatial distributions of potential temperature and vertical wind velocity as well as relative humidity and turbulence exchange coefficient over non-homogeneous underlying surfaces. It indicates that the model developed can be used to study the interaction between land surface process and atmospheric boundary layer over various underlying surfaces and can be extended to local climate studies. This work will settle a solid foundation for coupling climate models with the biosphere.     

Numerical Investigations of Mean Winds Within Canopies of Regularly Arrayed Cubical Buildings Under Neutral Stability Conditions

Recently, several attempts have been made to model the wind velocity in an urban canopy in order to accurately predict the mixing and transport of momentum, heat, and pollutants within and above the canopy on an urban scale. For this purpose, unverified assumptions made by Macdonald (Boundary-Layer Meteorol 97:25–45, 2000) to develop a model for the profile of the mean wind velocity within an urban canopy have been used. In the present study, in order to provide foundations for improving the urban canopy models, the properties of the spatially-averaged mean quantities used to make these assumptions have been investigated by performing large-eddy simulations (LES) of the airflow around square and staggered arrays of cubical blocks with the following plan area densities: λ _p = 0.05, 0.11, 0.16, 0.20, 0.25, and 0.33. The LES results confirm that the discrepancy between the spatial average of wind velocity and Macdonald’s five-point average of wind velocity can be large in both types of arrays for large λ _p . It is also confirmed that Prandtl’s mixing length varies significantly with height within the canopy, contrary to Macdonald’s assumption for both types of arrays and for both small and large λ _p . On the other hand, in accordance with Macdonald’s assumption, the sectional drag coefficient is found to be almost constant with height except in the case of staggered arrays with high λ _p .     

一个包含土壤和植被的区域气候模式及其性能检验

本文研制了一个三维地气耦合的区域气候模式，该模式由大气、土壤和植被三个子模式组成，各子模式间通过动量、能量和水份交换过程进行相互作用。在大气子模式中提出用z-z*混合坐标系处理地形的新方法，并把诊断分析中的散度修订方案引入数值模式，使其在物理依据上更加合理。此外，模式中较为详细地考虑了影响区域气候形成和变化的物理过程。简单的数值试验表明，区域气候模式比较真实地模拟了土壤和植被冠层的温度日变化，对区域气候平均态也有一定的模拟能力，可较好地处理大地形，具有良好的稳定性能，并对下垫面热力非均匀性的响应也是敏感的。     

从单Doppler速度场反演风矢量场的VAP方法

从单Doppler速度场反演风矢量场的VAP(Velocity Azimuth Processing)方法是在假定相邻方位角的风矢量相等的前提下,从单部Doppler雷达所测得的各个距离圈上的径向速度随方位角分布的廓线推算风向和风速的方法。这个假定可以保证反演所得的风矢量场中仍保留次天气尺度的风场信息。 对VAP方法反演公式的误差分析表明,反演出的风向误差为20—30度,与常规测风的误差相近。风速误差一般为10~0m/s,但在某些特定的情况下可达到与风速相同的量级(10~1m/s)。 用VAP方法反演风矢量场前必须首先消除原始Doppler速度资料中的脉动。它可以通过简单的一维线性平滑实现。 文中给出了一个反演实例,从中可以分析出次天气尺度的切变线和辐合线。常规天气图和UHF雷达的风廓线观测表明,这些次天气尺度系统是存在的。它们与反射率因子分布的配置还表明这些次天气尺度系统是和降水相联系的。     

积云团对大气层结影响的模拟研究

文中设计了一个二维完全弹性的积云团数值模式 .模式中考虑了水平均一和时间定常的大尺度抬升速度的影响 .大度抬升速度只对热力学方程和水汽守恒方程有影响 .积云团的启动方式采用冷池和随机扰动相结合的方法 .文中对一个给定温湿层结和大尺度抬升速度情况下 ,就不同扰动结构、水平环境切变、模式区域大小和冷池强度等进行了 1 1个模拟试验 .结果表明 ,在不同扰动结构、是否有水平环境风切变、模式区域不同以及不同云型谱时 ,给定温湿层结和大尺度抬升速度下的积云团活动使得温湿层结向着同一温湿层结调整 .这说明大气由于积云团活动达到准平衡状态 ,这个准平衡状态主要依赖于初始温湿层结和大尺度抬升速度 .这与 Betts积云参数化方案的假设相一致 .