带横隔板圆柱绕流特性数值模拟

为揭示尾迹区添加横隔板对圆柱绕流流场特性影响,把多步格式引入到特征线算子分裂有限元法中,建立了基于多步格式的特征线算子分裂有限元法:在每个时间步内将Navier-Stokes方程分裂成对流项和扩散项,对流项时间离散采用多步格式,在每一子时间步内沿特征线展开并显式求解。方腔流数值模拟结果表明该算法既可降低对整体时间步长的要求又可提高计算精度。对比有无横隔板圆柱绕流流场和圆柱表面压力变化表明,横隔板可以有效地抑制绕流尾迹区涡旋脱落,提高圆柱背流面压力,减少圆柱上下表面的压力差。     

坡面薄层流水动力学特性的实验研究

采用变坡实验水槽研究了坡面薄层水流水动力学特性——流态、流速、水深及阻力系数随流量和坡度的变化规律,实验坡度5°～25°、单宽流量为0.625×10-3～12.5×10-3m3/(s·m).研究结果表明,坡面薄层水流流态与水深密切相关,坡面流的流态基本上呈过渡流和紊流;薄层水流平均流速、水深和Darcy-Weisbach阻力系数主要受流量控制,坡度的影响并不显著,其关系均可用流量和坡度的幂函数形式模拟,相关系数(R2)分别为0.98、0.97和0.95,去掉坡度后相关系数分别下降7%、16%和3%.随着雷诺数的增大,Darcy-Weisbach阻力系数下降.当单宽流量小于0.008m3/(s·m)时,坡度对阻力系数的影响较为显著;当流量大于0.008m3/(s·m)时,阻力系数基本受流量控制.随着流量增大,阻力系数呈幂函数形式递减.     

地震荷载下危岩运动特征的模拟研究

西南地区地形复杂,位于地震高发区域。为了进一步研究该区域地震作用下落石的运动特征,以我国九寨沟景区内荷叶寨后山危岩带2为例。首先研究危岩带2中单一危岩体在自重荷载下的失稳情况,分别采用RocFall和非连续变形分析方法(DDA)进行模拟,结果表明:DDA和RocFall的模拟结果相近。DDA中落石水平位移模拟结果为210.6 m,RocFall中大部分落石水平位移分布在205 m附近。鉴于RocFall的局限性,采用DDA方法继续研究地震作用下危岩体的密度、高度、体积和地震波峰值加速度等因素对落石运动特征如运动轨迹、水平位移和速度的影响规律,并监测落石在坡脚附近的最大弹跳高度。研究结果表明:(1)危岩带2所处边坡自然坡角范围为65°~75°,落石在坡面上加速,在高程约2292 m(与居民区相对高差约为132 m)处发生抛射,运动方式以滚动、碰撞、弹跳为主;(2)岩体软弱结构面的贯通时间随落石密度、体积以及地震波峰值加速度增大而减小;(3)峰值速度随边坡高度、体积及地震波峰值加速度增大而增大;(4)水平位移随落石密度、边坡高度及地震波峰值加速度增大而增大;(5)坡脚处落石最大弹跳高度随各因素变量增大而增大。继续研究内部存在节理裂隙的岩体(3条平行节理将危岩分为K1~K4部分)在地震荷载下的运动特征:(1)在地震荷载作用下,危岩沿节理面发生倾倒破坏,落石水平位移模拟结果为191.5~243.8 m;(2)落石在坡脚附近的最大弹跳高度范围为0~13.8 m。     

用Gao-Yong湍流方程组数值模拟高雷诺数顶盖驱动方腔流

基于有限元方法中的虚位移原理,写出了不可压流Gao-Yong湍流方程组弱形式的微分方程表达式,利用有限元程序自动生成系统(FEPG)得到了全部有限元计算程序。对二维顶盖驱动方腔流进行了数值模拟,得到了方腔内的涡量分布和中心线上的速度变化曲线,并且成功地模拟出了高雷诺数时方腔内的双涡结构,计算结果与文献提供基准解符合良好。进一步证实了Gao-Yong湍流方程组能够准确模拟复杂湍流粘性场,并且为Gao-Yong湍流方程组的工程应用奠定了基础。     

四川茂县周场坪滑坡发育特征与变形监测分析

位于四川茂县南新镇的周场坪滑坡是一大型古滑坡,曾于1982年发生大规模快速复活,目前滑体半堰塞岷江。野外地质调查表明,周场坪滑坡在平面上呈不规则长舌形,长约850 m,滑体前后缘高差约350 m;在剖面上发育3级滑动,钻探揭露滑带埋深以50~70 m为主,推测潜在失稳滑坡体积为1 500×10 ⁴~2 000×10 ⁴ m ³。周场坪滑坡在平面上分为4个变形区,在滑体中后部和前缘坡脚发育大量拉裂缝与下错陡坎,拉裂缝宽度以0.2~3.0 m为主,陡坎下错高度2~10 m。在野外调查和钻探分析的基础上,对该滑坡开展了地表位移(GNSS)、深部蠕滑变形(钻孔测斜仪)和雨量等监测。监测分析表明,目前周场坪滑坡后缘变形速率达0.80 m/a,中部和前缘分别为0.69 m/a和0.51 m/a,呈推移式滑动变形,整体向NW310°方向滑动;地表位移速率在监测期内基本在1~3 mm/d之间波动,波动主要受降雨量影响,且略滞后于降雨量;滑移加速度基本在0~6 mm/d ²范围波动。ZK2钻孔测斜仪监测数据表明,滑坡在80 m深度内主要沿2层滑带蠕滑,其中浅层滑带埋深在22 m左右,深层滑带埋深在46 m左右,滑移速率在0~5 mm/d范围波动。综合研究认为,周场坪滑坡目前处于缓慢变形的深层蠕滑中,其变形速率受降雨和河流侵蚀等因素的影响,在极端内、外动力条件下,可能会加速滑动,并再次整体复活,造成堵塞岷江等危害。     

渗透率对干热岩开采过程储层变化规律的影响

热储层由基质系统和裂隙系统共同构成,二者热量传递的方式存在很大区别。仅考虑基质渗透率或裂隙渗透率,与实际采热过程并不相符。只有明确基质-裂隙双重渗透率下热储层的变化规律,才能更为合理、有效地开发地热资源。因此以青海共和盆地地热田GR1井为研究对象,基于热流固耦合理论,运用COMSOL数值模拟软件,建立双重孔隙介质渗透率水流传热模型。通过考虑不同基质渗透率(0,1×10-18,1×10-16m2)、裂隙渗透率(5×10-11,1×10-10,2×10-10m2),得到了储层温度场、应变场、应力场、位移场变化规律。研究发现:（1）仅考虑裂隙渗透率,会高估储层的开采寿命和产出温度;会低估采热过程中储层产生的压应变和沉降量,表明基质渗透率不能忽略。（2）最优裂隙渗透率为1×10-10m2,此时最适宜进行热开采;裂隙渗透率为2×10-10m2,储层寿命低于50 a。（3）采热初期,相比裂隙渗透率5×10-11m2时的最大压应变,裂隙渗透率为2×10-10m2时最大压应变提高了2. 74倍;采热40 a,相比裂隙渗透率为5×10-11m2,裂隙渗透率为2×10-10m2时,储层沉降量增加0. 164 05 m,沉降区域扩大3倍左右。所得结论对青海共和盆地干热岩开采过程中渗透率与储层变化规律的研究提供了一定参考。     

水射流自驱钻头自进能力影响因素研究

管道内流场结构特性是影响自驱钻头自进力和解堵除垢能力的关键因素。针对多股射流共同作用下管道内流场变化影响自进力和返水阻力的问题,基于FLUENT数值分析不同条件下管道内流场变化,得到后置喷嘴角度、转速、环空比和系统压力对自进力的影响规律,揭示综合摩擦因素理论值存在误差的原因,通过正交数值模拟明确自进力影响因素的主次关系。基于自进力测试实验装置,开展不同条件下自进力测试实验,验证数值模拟结果的正确性。结果表明：涡旋挤压碰撞作用和射流卷吸作用导致自驱钻头前后产生压力梯度,从而提高自进力。后置喷嘴倾角由20°增加至45°,壁面射流逆流速度由6.2 m/s增加至14.5 m/s,涡旋碰撞挤压后压力峰值由68.47 kPa增加至80.79 kPa,导致压差力增大。提高转速能够减小一次涡旋范围,当涡旋在旋转体区域发生碰撞挤压时,旋转体所受压差力增加。提高转速、系统压力或减小环空比均会增强后置射流的封隔作用,从而提高自进力。由于理论计算返水阻力时忽略了涡旋区的影响,因此综合摩擦因数理论值需要修正。环空比、系统压力、转速和后置喷嘴倾角的极差分别为300.07、111.87、60.42、36.32,影响...     

FeCl3膨胀石墨层间化合物制备及电导性能研究

在N2保护下用熔盐法合成了FeCl_3膨胀石墨层间化合物(简称FeCl_3-EGICs)。研究表明,在FeCl3过量、360℃、12h的条件下可得到纯一阶产物;460℃、12h为纯二阶。293K时,一阶FeCl3_EGICs的电导率为7.513×106S·m-1,二阶为6.87×106S·m-1,分别是原料膨胀石墨的3.7倍和3.4倍。在293K至分解温度范围内,FeCl3_EGICs具有金属电导特性,其电阻与温度间呈线性关系,随温度升高电阻增大,其机理为声子散射。     

三江平原地下水流场演化趋势及影响因素

近年来,三江平原部分区域地下水位呈持续下降的态势,引起广大学者和相关部门的高度关注,为了查明三江平原地下水流场时空演化规律,揭示其主要影响因素,以三江平原1980年72组及2019、2020年同期1 092组地下水位统测数据和44组国家地下水监测数据为基础,应用ArcGIS插值分析、栅格代数运算、对比分析等方法,查明了三江平原地下水流场时空演化特征,阐明了不同影响因素对地下水流场演化的控制作用。结果表明:与1980年相比,三江平原地下水位整体呈下降趋势。西部平原区地下水位累计降幅1～5 m的区域面积2.6×104 km2;东部建三江垦区累计降幅大于5 m的区域面积为1.17×104 km2,其中,累计降幅大于10 m的区域面积为3 400 km2,地下水位年均降幅约0.29 m。地下水开采引起地下水流场变化,“西砂、东黏”的水文地质条件促使区域地下水降幅的时空演化差异;水田种植规模的不断扩大引起地下水超采,浅地表黏土层阻挡降水入渗补给,地下水无法实现以丰补欠自平衡,造成了建三江垦区地下水位的持续下降。本研究成果为进一步查清三江平原地下水变化规律和现状特征奠定了基础,为科学指导地下水资源合理开发利用和管理提供了技术支持。     

四川茂县新磨村高位滑坡铲刮作用分析

2017年6月24日,四川省茂县叠溪镇新磨村发生高位顺层山体滑坡,滑动高差达1 160 m,滑动平距约2 200 m。该滑坡的滑动方量巨大,与其滑动过程中产生的铲刮效应有关。为分析其铲刮效应,文章通过现场调查、遥感影像解译和无人机航拍图像,确定该滑坡的滑动全过程为:多次历史地震造成滑坡源区岩体结构破碎,降雨沿顶部裂隙入渗导致水压力增大及石英砂岩中的薄层板岩软化,在长期疲劳效应下斜坡上部岩体最终发生滑动;上部滑体在运移过程中,对斜坡中部浅表风化层、部分基岩及下部老滑坡堆积体进行铲刮并重新堆积。采用Rockfall软件模拟源区滑体的运动路径、速度与能量,结果表明:在碎屑流区和老滑坡堆积区都存在明显的集中铲刮作用,整个滑坡的高危险区也主要位于该区域,所以危险性分区可代表不同滑坡区域的铲刮程度。计算得两个区域的铲刮方量分别为4.9×106,4.38×106 m3,滑坡总方量为13.35×106 m3。该模拟和计算方法迅速有效,可为以后类似滑坡的应急、救灾和铲刮方量计算提供参考。     

西北天山莱历斯高尔岩体年代学、地球化学及其成因

西北天山莱历斯高尔岩体主要由花岗闪长斑岩和二长花岗斑岩组成,二者具有相似的地球化学特征,可能为同源岩浆分异演化的产物。锆石SHRIMP U-Pb年龄显示岩体形成于晚泥盆世(374±4 Ma)。莱历斯高尔岩体具有较高的SiO₂(67.75% ～74.71%)含量和K₂O/Na₂O(1.24～2.20)比值,属于高钾钙碱性系列。稀土元素配分曲线为轻稀土富集型(La/Yb)_N(5.25～10.57),具有中等-弱的Eu负异常δEu(0.55～0.92)。微量元素特征显示富集大离子亲石元素(LILE),而亏损高场强元素(HFSE),具有较低的Sr(122×10^-6～356×10^-6)含量、较高的Y(17.19×10^-6～21.82×10^-6)和Yb(1.78×10^-6～2.57×10^-6)含量。岩体的Sr-Nd同位素特征为I_Sr=0.707 867～0.709 654,ε_Nd(t)=-2.79～-1.46。铅同位素特征为²⁰⁶Pb/²⁰⁴ Pb=18.423～19.915,²⁰⁷ Pb/²⁰⁴Pb=15.576～15.685,²⁰⁸ Pb/²⁰⁴Pb=38.344～39.305。元素和同位素地球化学特征表明,莱历斯高尔岩体岩浆源区主要为中元古代下地壳的部分熔融,并有部分地幔物质的加入。岩体形成于晚泥盆世准噶尔洋向伊犁-中天山微板块陡角度俯冲的大陆弧构造环境。     

黑潮与毗邻陆架海域的碳交换

作为沟通陆架边缘海与开放大洋的北太平洋西边界流,黑潮对毗邻中国陆架海域碳循环有着不可忽视的影响。分析总结了黑潮与毗邻陆架海域海水中溶解有机碳(DOC)、颗粒有机碳(POC)、溶解无机碳(DIC)和颗粒无机碳(PIC)分布与交换的现有研究成果。分析表明,黑潮通过表层水涡旋、入侵以及次表层以下水体上涌等形式携碳进入毗邻陆架海域。在黑潮与毗邻陆架海域,DOC和POC的分布受生物活动、水团混合以及沉积物—水界面相互作用的共同影响,表层与陆架水中含量较高,DIC受海水碳酸盐体系平衡影响,表层含量低于底层,在陆架海域与上升流区有高值,PIC含量较低。在台湾东北海域,DOC,POC和DIC由东海陆架海域向黑潮有净输出,输出量分别为2 431×109,1 051×109和21 733×109mol C/a,PIC通过陆架源悬浮颗粒物向黑潮主流输送,通量为1 852×109mol C/a。DOC,POC和PIC由南海经吕宋海峡向黑潮有净输出,输出量分别为2 652×109,1 009×109和230×109mol C/a,但DIC由黑潮向南海有净输出,通量为13 005×109mol C/a。台湾海峡碳循环受黑潮输运的影响,但影响机制及过程至今并不清楚。因此,黑潮与毗邻陆架海域的碳交换过程研究还有若干重大科学问题亟待解决,尚需深入系统的研究。     

青藏高原北部白榴碧玄岩年代学及地球化学研究

本文在可可西里北部边缘银顶山地区首次发现含橄榄石的高镁高钾质超基性火山岩。火山岩在阿尼玛卿-昆仑-木孜塔格缝合带南侧呈熔岩残丘形式产出,分布面积约为0.2 km²,其形成时代为上新世(5～4 Ma),是岩浆的快速侵位-喷溢的结果。斑晶主要为橄榄石(10%～12%)+白榴石(21%～27%)+霞石(13%～16%)±透辉石(29%～31%),基质主要为透辉石+云母+Fe-Ti氧化物+玻璃,定名为白榴碧玄岩。对3个样品51个点的橄榄石电子探针数据计算显示,橄榄石牌号为Fo_74-87Fa_13-26,属贵橄榄石,岩石中的橄榄石均不存在扭折带,表明这些橄榄石很可能是岩浆结晶的产物,而非岩浆捕获体。通过计算,银顶山地区橄榄石的结晶温度大约在1 226 ℃～1 234 ℃左右。白榴碧玄岩的化学组分:SiO₂<45%,MgO>10%, Na₂O+K₂O>8%,Mg^#(68.84～70.80),而Cr(149×10^-6～156×10^-6)和Ni(138×10^-6～151×10^-6),⁸⁷Sr/⁸⁶Sr(约0.708 6)和¹⁴³Nd/¹⁴⁴ Nd(约0.512 4),处于EMII源区范围内。初步认为藏北地区上新世超基性火山岩来自于较深的源区,最有可能是软流圈顶部地幔的局部熔融,并受到俯冲地壳物质的交代。     

黄渤海有机碳的分布特征及收支评估研究

陆架边缘海是陆海相互作用研究中最为关键的区域,也是全球重要的碳储库,在区域物质循环过程中发挥着重要的作用。基于2012年5月和11月对黄渤海海域的综合调查,对该海域水体和沉积物中有机碳的含量与分布进行了分析,并结合相关文献资料对黄渤海有机碳的收支进行了估算。主要结论为:黄渤海溶解有机碳和颗粒有机碳均呈近岸河口区域高、离岸低的分布趋势;有机碳的组成以溶解有机碳为主,颗粒有机碳由海洋自生的有机碳和陆地来源的有机碳组成;黄渤海沉积物有机碳高值区主要分布在河口和泥质区,其组成也是由海洋自生和陆源混合而成,其中渤海以陆源为主,而黄海以海源为主。黄渤海有机碳收支评估表明,有机碳的主要来源为初级生产力产生的有机物,其贡献为(6 760±971)×104t/a,占有机碳输入总量的(74±10)%,沉积物再悬浮的通量为(884±200)×104t/a,东海向黄海输入的通量为(679±107)×104t/a,河流及陆源输入的通量为(643±63)×104t/a,大气干湿沉降的通量为(141±39)×104t/a,其贡献分别占有机碳输入总量的(10±2.2)%,(7.5±1.2)%,(7.0±0.7)%和(1.5±0.4)%;黄渤海有机碳的主要支出为呼吸消耗,其贡献为(5 190±746)×104t/a,占有机碳输出总量的(57±8.2)%,黄海向东海输出的通量为(2 150±370)×104t/a,有机碳沉积通量为(1 030±225)×104t/a,有机碳降解通量为(737±191)×104t/a,其贡献分别占有机碳输出总量的(24±4.1)%,(11±2.5)%和(8.0±2.1)%。有机碳收支评估表明黄渤海有机碳以海洋自生来源为主,且具有潜在碳的"汇"的特性,水体中外源输入和海洋自生有机碳的(1.6±0.3)%埋藏于该海域内。     

新疆开都-孔雀河流域绿洲需水量与稳定性分析

水是绿洲存在和发展的核心, 干旱区绿洲稳定性与水密切相关. 根据2000-2009年资料, 采用蒸发系数法和定额法估算开都-孔雀河流域绿洲自然生态系统和社会经济系统综合需水量, 并对水资源约束条件下的绿洲稳定性进行初步探讨. 结果表明: 2000-2009年, 绿洲年均总需水量理论值约为54.80×10⁸ m³, 其中开都河绿洲总需水量约为20.55×10⁸ m³, 孔雀河绿洲总需水量约为21.90×10⁸ m³, 博斯腾湖区耗水量约为12.35×10⁸ m³, 与绿洲10 a平均供水量相比, 供需表现出极大地不平衡性. 水资源可承载绿洲面积(不含博斯腾湖)约为3139.66 km², 其中可承载灌溉地面积约为1395.41 km², 与绿洲10 a平均面积5 248 km²相比, 差别较大, 绿洲处于不稳定状态, 现状绿洲面积应适当收缩. 最后, 对博斯腾湖最低生态水位进行讨论, 初步把大湖最低水位定为海拔1 045 m, 小湖最低生态水位定为海拔1 046.5 m.     

2003-2011年青藏高原佩枯错相对水量变化及其对气候变化的响应

湖泊的退缩与扩张是全球气候变化的指示器.利用2003-2011年Landsat ETM数据和2003-2009年ICESat激光测高数据, 分别对青藏高原佩枯错湖泊的面积和高程变化进行了分析, 并进一步估算了湖泊2003-2009年相对水量变化.结果表明: 佩枯错面积年内变化明显, 湖泊面积冬季最小, 春季出现小峰值, 秋季达到最大; 面积年内波动明显(1.18%), 但在冬季、 春季和秋季相对稳定, 波动范围分别为0.26%、 0.1%和0.29%. 2003-2011年湖泊呈退缩趋势, 冬季、 春季和秋季面积年际变化率分别为-0.52 km²·a^-1、-0.35 km²·a^-1和-0.61 km²·a^-1; 2003-2009年间湖泊水位下降了1.17 m, 变化率为-0.05 m·a^-1; 2003-2010年, 冬季总水量减少了2.51×10⁸ m³, 春季总水量减少了1.74×10⁸m³, 秋季总水量减少了2.80×10⁸ m³, 平均相对水量变化率分别为-0.35×10⁸ m³·a^-1、-0.21×10⁸ m³·a^-1、-0.37×10⁸ m³·a^-1. 从空间上看, 湖泊退缩主要发生在东北角、 东南角和西南角.气候因素分析表明, 佩枯错湖泊退缩秋季主要是因为夏半年平均气温的升高, 冬季和春季则主要是因为冬半年降水量的减少.     

油田地热资源评价研究新进展

油田地热资源评价和开发利用一直是我国陆上沉积盆地内中-低温水热型地热资源研究的重点工作。本文介绍了油田地热资源评价的基本思路和评价方法,并以大庆油田、华北油田和辽河油田为例,通过地温场研究、热储特征分析、岩石热物性测试,按照新建立的油田地热资源评价分级体系,重新评价了油田区的地热资源量。3大油田区11个层系地质资源总量为10 934×10¹⁸ J,地热水资源量为86 607×10⁸ m³,其中,可采水资源量为19 322×10⁸ m³,可采地热能资源为425×10¹⁸ J。根据油田采出水和综合利用等5种开发利用方式,评估了油田地热能开发潜力,为我国地热资源评价和油田地热规模化开发奠定了基础。     

黑河下游额济纳绿洲土地生产潜力的动态变化及影响因素分析

在通过Visual MODFLOW软件模拟地下水位的基础上, 采用GIS技术进行栅格运算, 把影响土地生产潜力的各个因子作为"修正"系数, 运用逆向因子修正法分析黑河下游土地生产潜力的变化.在计算不同影响因素组合条件下土地生产潜力的同时, 结合不同的分水方案, 对各影响因素进行分析.结果表明: 在狼心山来水量为7.5×10⁸m³·a^-1时, 区域的水分有效系数整体增加, 植被覆盖度增大, 土地生产潜力大幅度提高, 能有效改善当地的生态环境.与现状相比, 在来水量为2.5×10⁸m³·a^-1时, 东河下游的水分有效系数出现降低的情形; 在来水量为5.0×10⁸m³·a^-1和7.5×10⁸m³·a^-1 时, 其水分有效系数均显著增加. 通过对影响因素的剖析有利于预警当地荒漠化发展的趋势, 进一步建立土地沙漠化的判别指标体系, 从而为黑河流域的土地生产服务.     

近20a来黑河中游张掖盆地地下水动态变化

在研究黑河中游张掖盆地水文地质条件及对区内地下水进行动态分区的基础上, 总结了区内1981-2001年地下水水位、储量的时空变化规律. 结果表明: 1)盆地内部沿冲洪积扇自上而下的 4 个区域内由于赋存条件的差异, 地下水多年动态类型分异显著. 其中, 冲洪积扇中上部 20 世纪 80 年代的10 a内地下水储量累积减少56.47×10⁶m³, 90年代地下水储量累积减少150.83×10⁶m³; 冲洪积扇下部80年代储量减少11.80×10⁶m³, 90年代减少77.50×10⁶m³; 细土平原 20 a 内储量累积减少0.56×10⁶m³; 临高河谷平原近20 a内储量增加了 1.78×10⁶m³; 进入 90 年代后, 各区域变化幅度均有不同程度增加. 2)各区域内地下水年内动态类型复杂, 存在显著的季节差异, 在进入 90 年代后年内动态变化也均有不同程度的增强.     

Advances in Study of Phytolith Carbon Sequestration in Terrestrial Ecosystems of China

Phytolith-Occluded Carbon (PhytOC), a relatively stable carbon fraction, plays an important role in biogeochemical carbon cycle and mitigation of global warming. The formation mechanisms of PhytOC, the influence factors of phytolith carbon sequestration, the advances in study of phytolith carbon sequestration, and the management for enhancing the potential of phytolith carbon sequestration in terrestrial ecosystems of China were summarized in this review. Finally, future researches on phytolith carbon sequestration in terrestrial ecosystems of China were prospected. Climates, vegetation types, soil circumstances, the chemical compositions of the phytoliths and other factors will directly or indirectly affect the potential of phytolith carbon sequestration. In China, the PhytOC production quantity in grassland, cropland, forest, wetland and shrub ecosystems is (0.6±0.1)×10⁶,(4.9±1.7)×10⁶,(1.9±0.4)×10⁶,(0.6±0.5)×10⁶ and (1.3±0.3)×10⁶ t CO₂/a, respectively. Application of silicon-containing fertilizer, cultivation of high-silicon plant, and traditional enhancement of the plant net primary productivity can significantly improve the potential of phytolith carbon sequestration in terrestrial ecosystems of China. Future studies should focus more on ①the mechanisms of phytolith formation in different plants, ②the phytolith carbon sequestration in underground parts of plants from different terrestrial ecosystems, ③the quantification of soil PhytOC in different terrestrial ecosystems. Furthermore, more comprehensive, economical and reasonable management practices of improving the potential of phytolith carbon sequestration should also be further studied in different terrestrial ecosystems.