CoLM模式对塔克拉玛干沙漠北缘陆面过程模拟评估及修正

采用2011年3月22日至7月26日肖塘陆气相互作用观测资料测试了不同陆面参数对公共陆面模式（CoLM）模拟效果的影响。地表参数包括地表反照率（α）、地表比辐射率（ε）、空气动力学和热力学粗糙度（z_0m、z_0h）、零置位移（d）、热传输附加阻尼（kB^-1）。结果表明：感热通量H和地表温度T_g对地表反照率、动力学粗糙度和地表比辐射率比较敏感,对零置位移不敏感。通过观测资料获得的这些参数及kB^-1参数化方案均被用来替换原CoLM模式中相应值及参数化方案。CoLM模式基本上能较好地模拟净辐射R_n、感热通量H、地表土壤热通量G₀和地表温度T_g日变化特征,只是在日峰值及其峰值出现时间的模拟上不理想;而CoLM模式对该地区土壤湿度M的模拟效果非常不好。误差统计值Bias、SEE、NSEE表明优化参数后的CoLM模式使得R_n的模拟误差被降低4.21%（SEE）,H的模拟误差被降低25.19%（Bias）,T_g的模拟误差被降低33.33%（Bias）、10.45%（SEE）和25%（NSEE）。     

固定沙质床面的空气动力学粗糙度

通过固定沙质床面的空气动力学粗糙度z₀的风洞实验研究,结果表明,z₀与沙粒粒径的关系为:z₀=a₁+b₁D²ln(D),并非Bagnold的1/30定律,也不是其他研究者的z₀与沙粒粒径D的固定关系。Bagnold的1/30定律仅对有限风速情况和非常狭窄的颗粒粒径范围有效。风也是影响空气动力学粗糙度的活跃因子。z₀与沙粒粒径,z₀与摩阻速度,z₀与摩阻雷诺数之间的关系都与风速有关。类似于水利学粗糙度,将空气动力学粗糙度划分为空气动力学完全光滑区、过渡区和完全粗糙区三个区域。     

砾石床面的空气动力学粗糙度

通过砾石床面的空气动力学粗糙度(Z₀)的风洞实验研究,结果表明,砾石床面的空气动力学粗糙度(Z₀)与砾石粒径、砾石覆盖度和自由风速有关,Bagnold的1/30定律和其它有关空气动力学粗糙度与粗糙元高度的固定比例对砾石床面都不适用。在各种一定砾石覆盖度条件下,砾石床面的Z₀随自由风速的增加而呈指数衰减。在各种一定自由风速条件下,Z₀随砾石覆盖度C的变化遵循二次曲线:Z₀=F₁+F₂C+F₃C^1.5+F₄C²,砾石覆盖度为40%~75%时,Z₀达到最大值。建立了包含自由风速和砾石覆盖度两个因子的Z₀的双因子综合模型。     

数值模拟气流特征和砾石几何参数对床面空气动力学粗糙度的影响

本文应用流体计算软件FLUENT6.3,采用非结构化网格划分技术模拟了气流特征和砾石几何参数对床面空气动力学粗糙度的影响。结果表明：空气动力学粗糙度（z₀）与风速（u）、摩阻速度（u_*）之间存在定量关系：z₀=a exp（bu/u_*）。砾石高度对空气动力学粗糙度的影响显著优于砾石直径,空气动力学粗糙度随砾石密度的变化比较复杂,先增加后减小。FLUENT在模拟风洞砾石床面动力学过程中的成功应用,是我们在研究方法上的一次有益尝试。     

粗糙度研究的现状及展望

略述边界层剪切湍流各区域速度廓线规律和地面粗糙分类后,对各种地面粗糙度(包括沙质粗糙度、动力粗糙度、植被粗糙度、复杂地面粗糙度、有效粗糙度)的研究现状和计算方法作了简要综述,认为动力粗糙度算式中系数A和植被或有效粗糙度算式中位移高度d将会是进一步研究的重点。     

科尔沁沙地风沙土地温特征及其分析

通过对科尔沁沙地3种类型风沙土中地温的变化特征分析,有植被覆盖和裸露地表条件下地温的关系探讨。结果表明:①地温随时间进程有明显的日变化规律,随深度加深地温日变化振幅按指数规律减小;②半固定风沙土的浅层(5~20 cm)地温略高于流动风沙土和固定风沙土,尤其表现在地温达到最大值时;③流动风沙土土壤热扩散率λ较大,固定风沙土土壤热扩散率λ较小;4有植被覆盖的风沙土中平均地温略高于裸露地表下的平均地温。     

植被恢复方式对黄土丘陵区土壤理化性质及微生物特性的影响

对甘肃省定西市黄土丘陵区不同植被恢复方式下的土壤水分、pH、有机碳（SOC）、全氮（TN）、碱解氮（AN）、全磷（TP）、速效磷（AP）、速效钾（AK）及微生物生物量、呼吸速率、呼吸熵（qCO₂）以及微生物熵（q_MB）进行了测定与分析。结果表明：不同植被恢复方式下,土壤pH值变化差异不显著（p>0.05）;0~20 cm土SOC、TN、AN、TP、AP、微生物量碳（SMBC）、微生物量氮（SMBN）含量及微生物呼吸速率（MR）均为天然植被（CL）下最高,AK为自然恢复（Q）方式下最高,qCO₂与q_MB在人工重建+自然恢复（RZ）方式下最高;20~40 cm土壤的SOC、TN、AN、TP、AP含量变化同0~20 cm,AK、SMBC、SMBN、 q_MB及MR在RZ方式下最高;40~60 cm土壤的SOC、TN、AN、AP及TP含量在CL方式下最高,SMBC、SMBN、 q_MB及MR在RZ方式下最高,AK、 qCO₂在Q方式下最高。相同植被恢复方式下,随土层深度增加,SOC、TN、AN、TP、AP、AK、SMBC、SMBN含量及q_MB、MR均逐渐降低,但qCO₂逐渐升高（除RZ方式）。不同植被恢复方式,SOC、TN、AN、SMBC、SMBN、q_MB、MR及qCO₂之间存在不同程度的相关性。     

准噶尔中南部沙尘暴源区地表沉积物粒度特征

对准噶尔中南部沙尘暴源区——艾比湖地区、克拉玛依农业开发区和古尔班通古特沙漠地表沉积物粒度组成、传统粒度参数及分维值进行了分析计算,结果表明:艾比湖流域及克拉玛依开发区广泛分布的沙质、粘土质、盐土荒漠区,特别是春季翻耕的大面积裸露农田地表物质含有大量小于0.063mm的粉砂和粘土,且粒度分维值较高,是沙尘暴的丰富物源。古尔班通古特沙漠风成沙粒径较粗,分维值较低,加之较好的植被状况,沙尘暴发生频数相对较低。粒度参数C-M图、M_z-S₀散点图和粒度分维值可以有效地表征地表松散物质的搬运-沉积特征,是描述沙尘暴发生地表条件的重要指标。     

青海喇家遗址地层划分及齐家文化废墟覆盖层成因分析

对喇家遗址进行了全面的野外考察,在其北部发现完整连续的全新世黄土-古土壤地层剖面。通过系统采样、粒度和磁化率等多种指标测试,结合OSL技术测年断代,建立了完整的土壤地层序列。深入分析土壤沉积物性质与地表过程,揭示了覆盖齐家文化废墟的沉积物的成因。研究结果表明：喇家遗址土壤沉积物地层层序自地表向下依次为：现代土壤（MS,0~1.50 ka B.P.）→全新世晚期风成黄土（L₀,1.50~3.10 ka B.P.）→全新世中期古土壤（S₀,3.10~8.50 ka B.P.）→全新世早期过渡性风成黄土（L_t,8.50~11.50 ka B.P.）→晚更新世晚期马兰黄土（L₁,>11.50 ka B.P.）。其中全新世中期黑垆土古土壤（S₀）被3组红色粘土质泥流（RC₃、RC₂和RC₁）和山洪沙土（FFD）沉积层穿插分隔为4个亚层（S_0上、S_0中上、S_0中下和S_0下）。这证明黑垆土发育过程中,在OSL年龄3.96~3.65 ka,对应于¹⁴C年龄3.85~3.60 ka B.P.出现一个山洪泥流盛行期,来自于盆地北侧沟谷的大规模暴雨山洪泥流,在大红山前古洪积扇前沿溢出沟槽、扩散覆盖了黄河第二级阶地面的土地。其中与大地震相伴随的第一期暴雨山洪泥流过程,摧毁了喇家遗址齐家文化聚落,造成居室内人群的死亡。     

沙地春玉米孕穗期和开花期土壤水分与产量的关系研究

采用回归旋转正交组合设计,对盆栽春玉米进行了孕穗期和开花期不同土壤水分处理试验。建立了孕穗期(X₁)和开花期(X₂)土壤水分与产量(Y)的关系模型:Y=90.5+7.98X₁+6.97X₂-2.25X₁X₂+6.38X₁²-7.38X₂²。孕穗期的土壤相对含水量不应低于54%。开花期的理想土壤相对含水量为65%。     

南极威德尔海西北区域的冰-气拖曳系数和脊帆形拖曳力分析

基于机载激光高度计测得的南极威德尔海西北区域冬季海冰表面起伏数据,结合脊帆形态参数和空间分布对中性条件(大气层处于中性平衡态)下对应10m高度处风速的冰-气拖曳系数Cdn(10)和脊帆形拖曳力Fd的参数化方案进行了改进,并探讨了冰-气拖曳系数Cdn(10)和脊帆形拖曳力对总拖曳力的贡献随脊帆强度(脊帆高度与间距的比值,其中脊帆高度是指脊帆顶点到平整冰面的垂直距离,脊帆间距是指相邻两个脊帆顶点之间的距离)和冰面粗糙长度的变化情况。结果表明,冰-气拖曳系数Cdn(10)随脊帆强度增大呈递增趋势,对较小的脊帆强度,Cdn(10)随粗糙长度的增大而增大,但脊帆强度较大时,Cdn(10)随粗糙长度减小而增大,即存在某一脊帆强度阈值,使Cdn(10)在该阈值两侧随粗糙长度的变化趋势相反;脊帆形拖曳力对总拖曳力的贡献随脊帆强度减小而减小,随着粗糙长度增大而减小。通过分析发现,造成以上不同变化趋势的主要原因是:随着脊帆强度的增大,摩拖曳力在总拖曳力中的优势地位逐渐由形拖曳力代替。     

孤立冰块流拖曳系数的优化辨识

流拖曳系数是海冰动力学模型中的重要参数,实际海冰运动时受到水流的底面剪切作用和侧面正压力作用,相应的拖曳系数分别为摩拖曳系数和形拖曳系数。利用实验室波流水槽设备,对长方平底、长方粗糙底和圆饼型淡水冰试样进行拖曳运动试验,获得26组冰样运动数据。根据这些试验数据,利用动量法原理和数学优化方法构造了摩拖曳系数和形拖曳系数的辨识模型,将摩拖曳系数和形拖曳系数两项分开。在此基础上,分析了冰样粗糙度指标(支承长度率和冰厚均方根偏差)对这两个系数的影响;发现了摩拖曳系数同冰样水平尺寸和垂直尺寸之比的关系,该关系能将冰山和浮冰的摩拖曳系数归结到同一条曲线上。     

风沙流中沙粒拖曳力系数研究

 拖曳力系数是计算风沙流中沙粒受空气阻力的重要参数。考虑实际风沙流中沙粒浓度及沙粒形状,利用FLUENT软件首先计算了不同风速下距地表不同高度处两沙粒的拖曳力系数,给出了影响沙粒拖曳力系数的间距范围,然后计算了真实风沙流中不同高度沙粒拖曳力系数。结果表明,给定风速下拖曳力系数随距地面高度的增加先减小后增加,并将沙粒拖曳力系数拟合成距床面高度的函数,该函数与风速有关。     

戈壁风蚀面与植被覆盖面地表性质粗糙度长度的确定

以Monin-Obukhov相似性理论为基础, 利用量纲分析法分别推导出不同层结稳定度下确定戈壁风蚀面与植被覆盖面空气动力学参数的物理模型, 并利用该模型研究了粗糙度长度与粗糙元性质, 流经近地层流体特征以及大气层结稳定度之间的关系。得出以下结论: 戈壁风蚀面上空气动力学粗糙度长度与砾石粒径、高度、覆盖度、自由风速、摩擦速度以及大气层结稳定度有关; 植被覆盖面空气动力学粗糙度长度取决于植被类型、植被高度、覆盖度、零平面位移高度、自由风速、摩擦速度以及大气层结稳定度。     

The Bowen Ratio of an Alpine Grassland in Three-River Headwaters,Qinghai-Tibet Plateau,from 2001 to 2018

The Bowen ratio (β) is used to quantify heat transfer from the land surface into the air, which is becoming a hot topic in research on the biogeophysical effects of land use and cover changes. The Three-River Headwaters (TRH), as a sensitive and fragile region, was selected as the study area. The β for 2001-2018 was estimated from the evapotranspiration product (ET_MOD16) of MODIS and the net radiation of the land surface through the albedo from GLASS. The ET_MOD16 data were evaluated against the observation data (ET_OBS) at two alpine grassland flux towers obtained from ChinaFLUX. The interannual trend of the β was analyzed by multiple linear regression (MLR) and structure model (SEM) with the multiple factors of precipitation, temperature, humidity, albedo, and normalized difference vegetation index (NDVI, MOD09Q1). The results show that the ET_MOD16 values were significantly correlated with ET_OBS, with a correlation coefficient above 0.70 (P < 0.01) for the two sites. In 2001-2018, the regional mean β was 2.52 ± 0.77 for the whole grassland, and its spatial distribution gradually increased from the eastern to western region. The interannual β showed a downward trend with a slope of -0.025 and a multiple regression coefficient (R ²) of 0.21 (P = 0.056). Most of the variability (51%) in the interannual β can be explained by the linear regression of the above multiple factors, and the temperature plays a dominant role for the whole region. The SEM analysis further shows that an increasing NDVI results in a decreasing albedo with a path coefficient of -0.57, because the albedo was negatively correlated with NDVI (R ² = 0.52, P < 0.01), which indicates a negative and indirect effect on β from vegetation restoration. An obvious warming climate was found to prompt more evapotranspiration, and restoring vegetation makes the land surface receive more radiation, which both resulted in a decreasing trend in the annual β. This study revealed the biogeophysical mechanisms of vegetation restoration under a changing climate, and demonstrated the Bowen ratio can be applied as an indicator of climate-regulating functions in ecosystem assessments.     

地形对包兰铁路沙坡头段防护体系的影响

选择沙坡头铁路防护体系典型断面,利用野外风场观测与风洞模拟,讨论地形在道路防沙体系中的作用。结果表明,防护体系前沿栅栏沙丘下风向10倍沙丘高度的水平范围内,地表摩阻风速通常小于临界值;风速递减系数随防护距离的增大而增大,地形倾斜和起伏对贴地层风力减弱的贡献占防护体系风速削弱程度的43%。自流沙区至铁路路基,潜在输沙率迅速波动降低,其中地形变化使得防护体系内潜在输沙率降低50%以上。这表明,地形是沙坡头铁路防护体系发挥卓越功能的重要因素,栅栏沙丘已成为防护体系的重要构成部分。     

Changes in intensity of wind erosion at different stages of degradation development in grasslands of Inner Mongolia, China

Overgrazing by increasing numbers of livestock in the Horqin Sandy Land of Inner Mongolia, China, has led to extensive degradation of the region's sandy steppes. Degraded grasslands are generally classified into four main types: fixed (least degradation), semi-fixed (light degradation), semi-shifting (moderate degradation) and shifting (severe degradation) sandy lands, representing four stages of degradation development. An experiment was conducted in the Horqin Sandy Land to investigate changes in intensity of wind erosion at different stages of degradation development in sandy grasslands and determine the extent to which surface wind erosion was affected by surface-related soil and vegetation factors through their effects on surface roughness length and wind regimes. Daily wind erosion rate was monitored at four sites of degraded grassland over an erosive period from 1 April to 10 June in 2001. Soil and vegetation properties for these sites were also measured twice: one in mid-April prior to the establishment of vegetation and again in mid-June after the establishment of vegetation. Relationships between surface roughness length and soil and vegetation variables were examined at each of the two stages of vegetation development. This study shows striking differences in the intensity of surface wind erosion among sites. The daily wind erosion rate in the fixed sandy land was, on average, only about 1/5 of the rate in the semi-fixed sandy land, 1/14 of the rate in the semi-shifting sandy land and 1/47 of the rate in the shifting sandy land suggesting a much higher resistance of the fixed sandy land to wind erosion compared to other sites. Differences in rate of wind erosion between sites were attributed to between-site differences in soil and vegetation properties that exerted significant effects on wind regimes by altering surface roughness length. At the pre-establishment stage of vegetation, surface roughness length was determined by a combination of litter amount on the ground, soil surface hardness and soil moisture content, with litter amount explaining the greatest proportion of the variation. At the post-establishment stage of vegetation, the development of the surface roughness effects was mainly governed by vegetation characteristics (vegetation cover in particular), while the effects of soil surface hardness and soil moisture on surface roughness length are likely to be masked by vegetation effects. The findings suggest that better management practices of restoring vegetation in degraded grasslands are required to reduce soil erosion losses and achieve a sustainable livestock production in the Horqin Sandy Land, an ecologically fragile sandy land ecosystem.     

确定坡面径流过程曼宁糙率系数的实验方法研究

曼宁糙率系数是用水动力学方法进行流速计算的关键参数。坡面流曼宁糙率系数与明渠流的不同。为确定坡面径流过程的曼宁糙率系数,自行研发了一种包括供水系统、实验水槽和数据观测记录系统的室内可变糙率坡面实验系统。通过87场预实验验证了供水系统的稳定性和准确性。以坡度、实测流量、实测水深、不同糙率板上河砂的平均直径和地表粗糙度为自变量,以曼宁糙率系数为因变量,选用均方根误差(RMSE)和决定系数(R ²)为评价指标,对166种实验场景进行了支持向量机(Support Vector Machines, SVM)训练与预测,发现：① 紊流的训练结果难以预测层流和过渡流的曼宁糙率系数,说明流态不同时,实验因素对水流的影响机制不同;② 若要较为准确地预测曼宁糙率系数,至少需要包括实测水深在内的3种因素;③ 当同时考虑4种及更多种因素时,紊流状态下均可对曼宁糙率系数进行较为准确的预测。