Evaluation of freeze-thaw wind erosion and analysis of influencing factors on the Tibetan Plateau based on the improved RWEQ
-
摘要:
冻融风蚀是指气温变化导致地表结构发生变化的过程，对地表地貌的形成和土壤的性质产生重要影响。然而，青藏高原地区的特殊性和复杂性使得冻融风蚀研究尚未得到广泛关注。针对青藏高原冻融风蚀定量化分析研究相对较少，本文采用引入冻结N因子的修正风蚀方程模型（Revised Wind Erosion Equation，RWEQ），分析2022年青藏高原冻融风蚀量的分布特征，并对青藏高原进行冻融风蚀强度分级。结果表明：①青藏高原冻融风蚀模数为8.90×108～4.95×105 （t·km−2·a−1），冻融风蚀总量为2.87×1013 t，冻融风蚀总面积2.41×106 km2；②青藏高原轻度以上冻融风蚀占全部冻融风蚀面积的84.26%，说明青藏高原冻融风蚀以轻微和轻度为主；③5个因子的显著性检验结果均为显著，对青藏高原冻融风蚀模数影响程度排序为气象因子>土壤结皮因子>土壤可蚀性因子>地表粗糙度因子>植被覆盖度因子。本研究不仅为青藏高原地区的环境保护和可持续发展提供了基础资料和实践参考，而且通过模型的应用，还为其他具有类似气候和地理特征的地区在冻融风蚀研究方面提供了多元方法和经验。
Abstract:Freeze-thaw wind erosion refers to the process in which temperature fluctuations cause changes in the surface structure, significantly impacting on the formation of surface landform and soil properties. However, due to the particularity and complexity of the Tibetan Plateau, the study of freeze-thaw wind erosion has not received widespread attention. Due to the relative scarcity of quantitative studies on the freeze-thaw wind erosion on the Tibetan Plateau, this study utilizes a Revised Wind Erosion Equation (RWEQ) that incorporates a freezing N-factor to analyze the distribution characteristics of freeze-thaw wind erosion on the Tibetan Plateau in 2022. The study also categorizes the intensity of freeze-thaw wind erosion across the Tibetan Plateau. The results indicate: (1) The freeze-thaw wind erosion modulus on the Qinghai- Tibet Plateau ranges from 8.90×108 to 4.95×105 (t·km−2·a−1), with a total erosion amount of 2.87×1013 t and a total affected area of 2.41×10 km2; (2) Moderate or greater freeze-thaw wind erosion accounts for 84.26% of the total affected area, indicating that light and mild erosion are predominant. (3) The significance tests of the five factors were all significant, and the influence of the freezing-thawing wind erosion modulus on the Qinghai-Tibet Plateau was ranked as meteorological factors > soil crust factors > soil erodibility factors > surface roughness factors>vegetation coverage factors. This study not only provides foundational data and practical references for environmental protection and sustainable development in the Qinghai-Tibet Plateau, but also offers a diverse range of methods and experiences for studying freeze-thaw wind erosion in other regions with similar climatic and geographical characteristics through the application of the model.
-
0. 引言
地裂缝为常见的地质灾害之一，一般指在内外应力协同作用下，岩土体的连续性被破坏，并发生变形破裂，进而在浅表层形成的具一定规模的变形带或裂缝的灾害现象，其变形破坏常呈现为水平剪切、垂直剪切与拉张破坏[1 − 3]。地裂缝在许多国家都有广泛发育，在美国、中国、墨西哥、澳大利亚以及埃塞俄比亚、肯尼亚等地均出现了大规模的地裂缝灾害，对当地居民生产生活造成严重影响，产生了难以计量的破坏和损失[4 − 10]。
调查发现，自上世纪60年代起，北京平原区勘察到地裂缝超过40条，至今已发育有7个地裂缝带，其中以高丽营、顺义和通州地区所发育的地裂缝最为典型，对周围建筑和公共设施造成了严重破坏，引起了人们的高度关注[11]。在北京地裂缝的发育特征和成因机理研究方面：贾三满等认为高丽营地裂缝主要受断裂的控制，为基底断裂活动延伸至地表导致的，是断裂构造和地面沉降共同作用的结果[12]；彭建兵等通过对顺义区土沟村地裂缝进行研究，认为该地裂缝与下伏断裂相连，并呈现出同沉积断层的活动特征[13]；刘德成等认为通州地裂缝是非构造性地裂缝，主要受地貌形态所控制，并分布于河流两岸，顺河堤方向展布[14]；姜媛等在对大量资料分析的基础上，得到了高丽营地裂缝和当地地面沉降在时间、空间和活动方式上的关联，认为二者在加剧活动时间上具有同步性，在活动特征上具有相似性，并指出地面沉降在断裂两侧引起的不均匀沉降是高丽营地裂缝形成的主要因素[15]；卢全中等对北七家至高丽营段的地裂缝进行研究，确定其上、下盘影响带总宽度为43 m，地裂缝对房屋与道路的破坏形式有3种[16]；刘方翠等认为是区域构造活动及差异沉降变形机制的共同作用导致了北京平原区顺义地裂缝、高丽营地裂缝、羊房地裂缝和北小营地裂缝的快速形成发育[17]；周永恒等通过地球物理探测技术与联合钻孔剖面方法揭示了黄庄-高丽营断裂（房山-涞水段）具典型的正断活动特征[18]；关金环等发现北京首都国际机场发生的差异性沉降主要是由顺义-良乡断裂所发育的顺义地裂缝造成的[19]；孟振江等以北京地裂缝为原型研究了耦合型地裂缝的发育活动特征和成因机制，揭示了不同位错量和水位下降量引发的地层位移场和应力场的变化特征[20]；任雅哲等依据断层位错理论定量分析了顺义断裂蠕滑活动对机场地裂缝形成的影响，认为顺义断裂上盘中、深层地下水抽采引起的地面差异沉降是导致机场地裂缝加剧的主要因素[21]。以上研究在一定程度上揭示了北京地裂缝的成因与影响因素，为后续的深入研究提供了重要的研究基础和理论参考。
作为特殊的工程场地，地裂缝场地对动力的响应明显且产生的影响较大。相关研究表明，傅里叶频率与场地条件、土层结构完整性及力学性质等紧密相关，而地脉动频率及卓越周期都可以反映出试验场地土体结构的变化[22 − 23]。在利用地脉动测试方法研究方面：薛捷等通过地脉动测试的方法对西安地裂缝场地进行了研究，从而为西安市抗震设防工作提供理论依据[24]；慕焕东等使用FLAC3D数值分析软件对西安地裂缝场地开展研究，发现地裂缝场地具有明显的放大效应，且上盘效应显著[25]；张磊刚等应用地脉动测试手段对西安F6地裂缝场地进行研究，探明了场地因地裂缝而产生的放大效应和峰值变化规律[26]；王晗等通过地脉动现场测试方法，得到了西安地区典型地裂缝的场地放大效应和变化规律[27]；崔思颖等选取河北平原具有代表性的地裂缝进行地脉动测试揭示了动力响应特征[28]；Chang J和Xuan Y等利用地脉动测试技术获得了山西临汾和大同两处盆地的典型地裂缝场地的数据并开展研究，研究结果表明可以通过地脉动信号来分析地裂缝场地的动力响应特征，并揭示了其变化规律与放大效应[29 − 32]。
在对北京地区的地裂缝研究方面，前人所做研究大多集中于对地裂缝破坏现象的描述和产状勘察，以及从发育场地地貌和地质构造角度对其成因机理进行研究，在利用地脉动测试北京地裂缝发育场地的动力响应研究方面还属于空白，因该方法的应用还未被普及，加上北京多个地裂缝段现场的干扰较大，故尚未开展专门的地脉动现场测试[33]。鉴于此，本文的研究对象选取为北京通州区宋庄地裂缝场地，采用现场踏勘、地脉动测试和动力分析等方法，开展针对性的宋庄地裂缝场地的动力响应研究，以探究地裂缝的存在对场地动力响应的“峰值”及其“卓越频率”所造成的影响，并在此基础上揭示地裂缝发育场地的动力响应规律，以期为宋庄地区场地工程建设和防震减灾工程的设计和规划提供理论依据，同时为其他地裂缝发育区开展此类研究提供借鉴与参考。
1. 北京地裂缝分布及活动特征
北京平原区地处太行山隆起与燕山隆起和华北断陷盆地相交汇的地带，自新生代以来，该地区主要处于伸展变形的构造环境，地质构造复杂且独特，断裂构造发育。近年来，随着人类工程活动的日益加剧，引发了多起地裂缝灾害[33]。
通州区地裂缝于上世纪60年代开始零星出现，唐山大地震发生后，郎府、西集、宋庄等地都出现了大量的地裂缝，其中以宋庄地裂缝最为典型。宋庄镇地处北京迭断陷与大兴迭隆起相交界的区域，其下伏具有活动迹象的早、中更新世南苑-通县断裂，该断裂带长约110 km，走向NE35～50°，倾向NW，倾角60°～75°，NW盘下降，SE盘上升，为正断层（图1）。断层西北侧震旦亚界埋深400～
1000 余米，上覆有巨厚的新近系、第四系沉积，而东南侧寒武系及中上元古界埋深仅60～300米，上覆有第四系及少量新近系沉积；对断层附近钻孔岩芯观察，具有明显挤压现象，应属压型断层。宋庄地裂缝于2010年后活动尤为剧烈[34]。经野外调查发现，宋庄地裂缝沿南苑-通县断裂展布，南起小中河在地表呈近直线延伸，经双埠头村、沟渠村、大庞村至平家疃村后向北东方向延伸，长度约8.7 km，整体走向NE31°，倾向NW。该地裂缝以水平拉张破坏为主，地表水平拉张量在2.15～56.7 cm之间，垂直位错量较小，地裂缝影响变形带最宽处达到400 m（图2）。![]()
图 1 宋庄地裂缝与周边断裂分布图Figure 1. Distribution map of mainly fractures and location of Songzhuang ground fissures![]()
图 2 宋庄地裂缝致灾示意图Figure 2. Schematic diagram of disaster caused by songzhuang ground fissures宋庄地裂缝的地层剖面如图3示，地表以下由粉砂、细砂、粉土、粉质黏土互层组成[35]。在浅表层，地裂缝近直立，地裂缝以偏拉张开裂为主，张开量为1～15 cm；在地下埋深50 m以下的地层开始逐渐出现位错现象，埋深53.2 m处的钻孔岩芯显示，在浅黑灰色砂质粘土层中，发育不连续构造面，构造面倾角64°，且两侧地层颜色不同。总体上分析发现地层上下盘位错量随埋深逐渐增大，地裂缝NW盘地层相对下降，西北侧第四系地层沉积厚度大于东南侧。在宋庄地裂缝影响范围内，致灾最为严重的是双埠头村至大庞村段，据统计，2016年、2018年和2022年地裂缝造成的直接经济损失总值分别为
4623.24 万元、7033.49 万元和10751.14 万元，可见地裂缝的影响和破坏在不断加剧[36]。![]()
图 3 宋庄地裂缝浅表部地层剖面图（据赵龙，2018修改）Figure 3. Stratigraphic section of the shallow strata of Songzhuang ground fissure (modified from Zhao Long, 2018)2. 地脉动测试
2.1 地脉动基本原理
地脉动通常是指地表面振动周期在0.05～10s、振幅在千分之几至几微米的微幅振动，是由环境振动在地球表面所产生的复杂随机振动，包括地球内部应力、地震等自然因素及交通等人类活动所引起的振动，并经过不同传波介质和不同场地环境后综合作用于地球表面的结果。目前地脉动研究大致分为两类：一类是常时微动，一般是指周期小于1s的微动，侧重其振幅、周期性与观测点场地分类和振动特性研究；另一类是长波微动，一般是指周期大于1s的微动，研究的重点是由微动提取面波，进而解释推断观测点地下横波波速1s的常时微动，可用以研究地裂缝场地动力响应特性。
2.2 地脉动测试
地脉动测试要求的频段范围在0.1～15Hz之间，本次进行地脉动测试所使用的设备是高灵敏度伺服型速度网络地震仪，可满足测试精度要求，主要规格参数见表1。
表 1 选取地震仪的主要技术参数Table 1. Main technical parameters of the Seismic Monitor仪器型号 CV-374AV 动态范围 136 db 灵敏度 1000 mv/（cm/s）通频带 0.1～100 Hz 测量范围 ±0.02 m/s 线性 0.03% 采样频率 100 Hz AD分辨频率 24 bit 高灵敏度伺服型速度网络地震仪尺寸为180*120*100 mm，灵敏度较高，在进行地脉动测试时，严格按照操作规程测试，确保测试数据的精确性。为避开不利环境因素和人类活动对地脉动监测结果的影响，选择天气条件良好情况下进行地脉动测试，测试工作选择在凌晨1点至5点期间。
选择在北京通州区宋庄镇发育较为典型的双埠头村地裂缝段进行地脉动测试，现场布设测点具体原则为：在与主裂缝走向垂直的方向布设2条测线，结合地裂缝的影响带宽度范围，布设长度为60 m的测线，其中，上、下盘长度各设30 m。将地裂缝作为测线的中点，向远离地裂缝的方向进行测点的布设。以上盘为例，首先，在距地裂缝1.5 m处布设第一个测点，再沿测线以3 m为间距布设测点，共布设5个；接着按每隔5 m布设一个测点，共布设3个。因此，所布设测点与地裂缝间的水平距离由近到远设成1.5、3、6、9、12、15、20、25和30 m。下盘测点与上盘对称布设。
针对地脉动观测各个测点的数据，在稳定波段内以10秒为单位进行波段的截取并通过Pwave32软件将其转换导出为加速度时程曲线，之后应用SeismoSignal软件对数据进行预处理，以得到反应谱曲线、傅里叶谱曲线以及Arias烈度曲线，通过频谱特征分析，进而获得地裂缝场地地脉动响应的加速度幅值、卓越频率、影响范围等场地动力响应特征；通过利用3种从不同角度出发的动力分析方法，对采集并经处理的地脉动数据进行分析，进一步提炼地脉动产生的场地动力特性。
3. 测试结果分析
3.1 傅里叶分析
通过傅里叶谱分析，可以将地脉动加速度时程曲线转换成傅里叶谱曲线，从而得到场地卓越频率和傅里叶谱峰值，通过对不同测点卓越频率和峰值的对比，可得到场地不同区域地脉动响应的分布特征，从而获得地裂缝场地动力响应特征。
3.1.1 双埠头村地裂缝傅里叶分析
图4、图5为双埠头村地裂缝测线上盘和下盘在三个方向上的傅里叶谱图（其中X、Y、Z分别为垂直地裂缝走向方向、平行地裂缝走向方向及和垂直地面方向，下文同），图中各测点的傅里叶谱峰值显著，呈现出多峰的特点。主频带宽度较大，而各测点的卓越频率相对集中，因此其波速的变化没有明显的规律。傅里叶谱图反映出的多峰、谱面积大以及谐波丰富等特征，这与双埠头村地裂缝发育场地中较为松软的粘性土层结构特征相一致。通过对双埠头村地裂缝上下盘的傅里叶谱统计分析，得到该区域地裂缝场地的卓越频率在2.10～4.30 Hz之间浮动。而场地卓越频率的大小和测点距地裂缝的水平距离间的关系并不显著，均值约为3.30 Hz。
![]()
图 4 双埠头村地裂缝上盘不同测点XYZ方向上傅里叶谱Figure 4. Fourier spectra in the XYZ direction at different measurement points on the hanging wall of the ground fissures in Shuangbutou village![]()
图 5 双埠头村地裂缝下盘不同测点XYZ方向上傅里叶谱Figure 5. Fourier spectra in the XYZ direction at different measurement points on the footwall of the ground fissures in Shuangbutou village图6、图7为双埠头村地裂缝各测点的傅里叶卓越峰值图。从图中可以看出，在距离地裂缝区域较近的位置，地脉动响应由于地裂缝的存在而出现了显著的放大效应，傅里叶谱峰值与该距离呈负相关关系，距地裂缝的距离越大，峰值越低，在距离较远处曲线趋于平缓，其临界值约在距地裂缝12 m处，傅里叶谱峰值大致保持稳定不变。在距离地裂缝12 m以内的区域，振动的幅值和能量密度相对较高，而在距离地裂缝12米以外的区域，地表振动的幅值和能量密度则相对较低，这种现象表明靠近地裂缝的区域地表振动会比远离地裂缝的区域更强烈。
![]()
图 6 双埠头村地裂缝上盘XYZ方向上傅里叶卓越峰值Figure 6. The peak value of Fourier spectrum of the hanging wall in XYZ directions in Shuangbutou village![]()
图 7 双埠头村地裂缝下盘XYZ方向上傅里叶卓越峰值Figure 7. The peak value of Fourier spectrum of the footwall in XYZ directions in Shuangbutou village3.1.2 沟渠村地裂缝傅里叶分析
图8、图9分别为沟渠村地裂缝各测点X、Y、Z三个方向上傅里叶谱峰值与距地裂缝的距离变化关系的散点及曲线图，发现距离地裂缝越远的位置则傅里叶谱峰值越小，直至在较远的位置，曲线逐渐接近平缓；最终在距地裂缝约12 m处，曲线基本保持不变，傅里叶谱峰值趋于一致。通过统计分析，得出沟渠村地裂缝的场地卓越频率在2.30～4.30 Hz范围内变化，卓越频率的高低与测点所在的位置之间没有明显关系，其均值约为3.40 Hz。
![]()
图 8 沟渠村地裂缝上盘XYZ方向上傅里叶卓越峰值Figure 8. The peak value of Fourier spectrum of the hanging wall in XYZ directions in Gouqu village![]()
图 9 沟渠村地裂缝下盘XYZ方向上傅里叶卓越峰值Figure 9. The peak value of Fourier spectrum of the footwall in XYZ directions in Gouqu village3.2 反应谱分析
通过反应谱分析，可以得到特征周期和相应的反应谱峰值，通过对不同测点特征周期和反应谱峰值的对比，可以得到场地不同区域地脉动响应的分布特征，从而获得地裂缝场地动力响应特征。
3.2.1 双埠头村地裂缝反应谱分析
经数据分析得到，双埠头村地裂缝各测点在同一方向上的反应谱谱型均为“单峰型”，且具有主峰突出、频带窄和谱面积小的特点；反应谱卓越周期集中在0.1 s～0.2 s之间，部分卓越周期可达0.3s，此外卓越周期值与地裂缝和测点间的水平距离无明显关系（图10、11）。
![]()
图 10 双埠头村地裂缝上盘不同测点XYZ方向上反应谱Figure 10. Response spectrum in the XYZ direction at different measurement points on the hanging wall of the ground fissures in Shuangbutou Village![]()
图 11 双埠头村地裂缝下盘不同测点XYZ方向上反应谱Figure 11. Response spectrum in the XYZ direction at different measurement points on the footwall of the ground fissures in Shuangbutou Village由图12、图13可得，在距离地裂缝较近处，因地裂缝的存在地脉动响应表现出显著的放大效应，随着与地裂缝距离的增大，地脉动的响应加速度基本不变，在距地裂缝12m处趋于平稳，反应谱峰值基本一致。
![]()
图 12 双埠头村地裂缝上盘XYZ方向上反应谱加速度峰值变化曲线Figure 12. Response spectrum acceleration peak value variation curves in XYZ direction on the hanging wall of the ground fissures in Shuangboutou village![]()
图 13 双埠头村地裂缝下盘XYZ方向上反应谱加速度峰值变化曲线Figure 13. Variation curves of peak values of reaction spectral acceleration in XYZ direction on the footwall of the ground fissures in Shuangbutou Village3.2.2 沟渠村地裂缝反应谱分析
由图14、图15可以看出，沟渠村地裂缝各测点在X、Y、Z三个方向的上反应谱加速峰值随与地裂缝间的垂直距离变化，距地裂缝较近的区域的地脉动响应表现出显著的放大效应，且距离地裂缝越远，反应谱峰值越小，直至曲线接近稳定不变，反应谱峰值基本一致。
![]()
图 14 沟渠村地裂缝上盘XYZ方向上反应谱加速度峰值变化曲线Figure 14. Variation Curve of Peak Acceleration in the XYZ Directions of the Response Spectrum on the Hanging Wall of the Ground Fissure in Gouqu Village![]()
图 15 沟渠村地裂缝下盘XYZ方向上反应谱加速度峰值变化曲线Figure 15. Variation Curve of Peak Acceleration in the XYZ Directions of the Response Spectrum on the Footwall of the Ground Fissure in Gouqu Village通过对比双埠头村地裂缝和沟渠村地裂缝反应谱曲线，发现布设在宋庄地裂缝上下盘的测点，在与地裂缝同等距离时，上盘反应谱的峰值均大于下盘，这表明地裂缝的存在使得地裂缝两侧场地的地脉动响应发生了差异变化，存在显著“上盘效应”，且上盘峰值大于下盘。
3.3 Arias烈度分析结果
Arias烈度是一种目前较为成熟的评判场地地震动响应强弱、反映场地抗震能力和设防标准的方法，通过分析Arias烈度曲线的变化特征进而了解地裂缝场地不同位置的抗震特性，可从场地抗震烈度方面揭示地裂缝对场地的放大效应。
3.3.1 双埠头村地裂缝Arias烈度分析
图16、图17分别为双埠头村地裂缝上盘与下盘不同测点的Arias烈度曲线图，图中注释栏“MP”表示监测点。结果表明，Arias烈度呈现出随时间逐渐增长的趋势，但各个测点的增长幅度有所不同，且有着不同的最大值。
![]()
图 16 双埠头村地裂缝上盘不同测点XYZ方向上Arias烈度曲线Figure 16. Arias intensity profiles in the XYZ direction at different measurement points on the hanging wall of the ground fissure in Shuangbutou village![]()
图 17 双埠头村地裂缝下盘不同测点XYZ方向上Arias烈度曲线Figure 17. Arias intensity profiles in the XYZ direction at different measurement points on the footwall of the ground fissure in Shuangbutou village为反映地裂缝对场地不同位置烈度值产生的影响，将双埠头村地裂缝各测点Arias烈度的峰值沿垂直地裂缝走向方向展开，得到上下盘各测点在X、Y、Z三个方向上的Arias烈度峰值随与地裂缝距离变化的散点数据曲线图（图18、19）。可见，场地不同位置的烈度值受地裂缝的影响程度也不同。Arias烈度在距离地裂缝较近的区域，其峰值较大，且该烈度值随测点与地裂缝间的距离增大而减小并逐渐趋于稳定，其临界距离为12 m，而临界区域以外的范围受放大效应的影响较小。
![]()
图 18 双埠头村地裂缝上盘XYZ方向上Arias烈度峰值随距离变化曲线Figure 18. Variation curves of peak Arias intensity with distance in the XYZ direction on the hanging wall of the ground fissures in Shuangbutou village![]()
图 19 双埠头村地裂缝下盘XYZ方向上Arias烈度峰值随距离变化曲线Figure 19. Variation curves of peak Arias intensity with distance in the XYZ direction on the footwall of the ground fissures in Shuangbutou village4. 结论
（1）宋庄地裂缝延伸8.7 km，地表破坏以拉张开裂为主，水平拉张量在2.15～56.7 cm之间；其中双埠头村地裂缝的傅里叶卓越频率主要集中在2.10～4.30 Hz之间，反应谱的卓越周期集中在0.1～0.2 s之间；沟渠村地裂缝的傅里叶卓越频率集中在2.30～4.30 Hz之间，反映了地裂缝不同区段活动性的变化。
（2）通过反应谱分析法、傅里叶谱分析法分析得到，测点的卓越频率无明显的规律性变化；Arias烈度峰值在靠近地裂缝处较大，且随着距离的增大而逐渐减小，在距地裂缝较远的区域趋于平缓。
（3）宋庄地裂缝场地动力响应在反应谱、傅里叶谱和Arias烈度上都存在明显的“放大效应”，即靠近地裂缝的区域大于远离地裂缝的区域，且该临界距离为距地裂缝12 m。
（4）宋庄地裂缝场地有着较为显著的“上盘效应”，即距地裂缝相同距离情况下，位于上盘的傅里叶卓越频率、反应谱卓越周期和Arias烈度峰值总大于下盘。
-
表 1 数据引用来源
Table 1 Data citation sources
数据 数据来源 数据精度 单位 青藏高原矢量边界 中国科学院资源环境科学与数据中心 青藏高原气象站点 中国科学院资源环境科学与数据中心 降雨数据 中国科学院资源环境科学与数据中心 日 mm 风速数据 中国科学院资源环境科学与数据中心 日 m/s 蒸散发量数据 时空三极环境大数据平台 月 mm 雪盖数据 中国科学院资源环境科学与数据中心 日 mm 地表温度数据 欧洲中期天气预报中心 日 °C 土壤数据 国家冰川冻土沙漠科学中心 土地利用数据 Zenodo平台 1km 数字高程模型数据 中国科学院资源环境科学与数据中心 1km 归一化植被指数 地球资源数据云平台 年 
下载: 导出CSV
表 2 青藏高原冻融风蚀强度分级
Table 2 Classification of freeze-thaw wind erosion intensity on the Tibetan Plateau
等级 轻微 轻度 中度 强烈 极强烈 剧烈 分级标准
/（t·km−2·a−1）<2 2～25 25～50 50～80 80～150 >150 分布面积/km2 1388638 645200 149399 111901 74017 44472 面积占比/% 57.53% 26.73% 6.19% 4.64% 3.07% 1.84%
下载: 导出CSV
表 3 137Cs核素示踪实验值与模型模拟值对比
Table 3 Comparison of experimental values from the 137Cs nuclide tracer and the model simulation values
样点编码 S1 S3 S4 S5 S6 S7 实验值/（t·km−2·a−1） 8414 6943 2262 3068 2020 4759 改进模型计算结果
/（t·km−2·a−1）9447 4725 2387 3083 2103 4819 传统模型计算结果
/（t·km−2·a−1）10167 5896 1823 2110 1531 3290
下载: 导出CSV
表 4 个变量间的相关性
Table 4 correlations between variables
$ SL $ $ WF $ $ EF $ $ SCF $ $ {K}' $ $ C $ 冻融风蚀模数$ SL $ 1 −0.027 −0.013 0.020 0.007 0.007 气象因子$ WF $ −0.027 1 0.007 0.052 0.053 −0.072 土壤可蚀性因子$ EF $ −0.013 0.007 1 0.188 0.024 −0.033 土壤结皮因子$ SCF $ 0.020 0.052 0.188 1 0.045 0 地表粗糙度因子$ {K}' $ 0.007 0.053 0.024 0.045 1 −0.062 植被覆盖度因子$ C $ 0.007 −0.072 −0.033 0.000 −0.062 1
下载: 导出CSV
表 5 变量的显著性及共线性
Table 5 Significance and collinearity of variables
皮尔逊相关性 显著性 95%置信区间 共线性统计 下限 上限 容差 VIF $ SL-WF $ 0.027 0.001 −0.029 −0.025 0.981 1.019 $ SL-EF $ 0.013 0.001 −0.011 −0.015 0.309 3.233 $ SL-SCF $ 0.020 0.001 0.018 0.022 0.309 3.231 $ SL-K' $ 0.007 0.001 0.005 0.009 0.987 1.014 $ SL-C $ 0.007 0.001 0.005 0.009 0.989 1.011
下载: 导出CSV
表 6 回归系数
Table 6 Regression coefficients
模型 未标准化系数 标准化系数 $ t $ 显著性 β 标准错误 β 常量 3.246 0.029 113.648 0.000 $ WF $ −0.012 0.000 −0.028 −28.883 0.001 $ EF $ −0.762 0.043 −0.017 −17.528 0.001 $ SCF $ 0.499 0.020 0.024 24.415 0.001 $ {K}' $ 0.152 0.018 0.008 8.527 0.001 $ C $ 0.094 0.018 0.005 5.120 0.001
下载: 导出CSV
-
[1] 杨珍珍,倪万魁,李兰,等. 基于GIS的多因子土壤冻融侵蚀强度分级评价——以长江上游为例[J]. 水利与建筑工程学报,2022,20(6):89 − 95. [YANG Zhenzhen,NI Wankui,LI Lan,et al. Grading Evaluation of Multi Factor Soil Freeze-thaw Erosion Intensity Based on GIS:A case study of the upper reache of the Yangtze River[J]. Journal of Water Resources and Architectural Engineering,2022,20(6):89 − 95. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1672-1144.2022.06.014 YANG Zhenzhen, NI Wankui, LI Lan, et al. Grading Evaluation of Multi Factor Soil Freeze-thaw Erosion Intensity Based on GIS: A case study of the upper reache of the Yangtze River［J］. Journal of Water Resources and Architectural Engineering, 2022, 20（6）: 89 − 95. （in Chinese with English abstract） DOI: 10.3969/j.issn.1672-1144.2022.06.014
[2] 魏霞,丁永建,李勋贵. 冻融侵蚀研究的回顾与展望[J]. 水土保持研究,2012,19(2):271 − 275. [WEI Xia,DING Yongjian,LI Xungui. Review and prospect of freeze-thaw-induced erosion research[J]. Research of Soil and Water Conservation,2012,19(2):271 − 275. (in Chinese with English abstract)] WEI Xia, DING Yongjian, LI Xungui. Review and prospect of freeze-thaw-induced erosion research［J］. Research of Soil and Water Conservation, 2012, 19（2）: 271 − 275. （in Chinese with English abstract）
[3] 景国臣. 冻融侵蚀的类型及其特征研究[J]. 中国水土保持,2003(10):17 − 18. [JING Guochen. Study on types of freeze-thaw erosion and its characteristics[J]. Soil and Water Conservation in China,2003(10):17 − 18. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-0941.2003.10.010 JING Guochen. Study on types of freeze-thaw erosion and its characteristics［J］. Soil and Water Conservation in China, 2003（10）: 17 − 18. （in Chinese with English abstract） DOI: 10.3969/j.issn.1000-0941.2003.10.010
[4] 谢胜波,屈建军,韩庆杰. 青藏高原冻融风蚀形成机理的实验研究[J]. 水土保持通报,2012,32(2):64 − 68. [XIE Shengbo,QU Jianjun,HAN Qingjie. Mechanisms of freezing-thawing induced wind erosion in Qinghai-Tibet Plateau[J]. Bulletin of Soil and Water Conservation,2012,32(2):64 − 68. (in Chinese with English abstract)] XIE Shengbo, QU Jianjun, HAN Qingjie. Mechanisms of freezing-thawing induced wind erosion in Qinghai-Tibet Plateau［J］. Bulletin of Soil and Water Conservation, 2012, 32（2）: 64 − 68. （in Chinese with English abstract）
[5] BAGNOLD R A. The Physics of Blown Sand and Desert Dunes[M]. DordrechtSpringer Netherlands,1971.
[6] CHEPIL W S,Milne R A. Comparative study of soil drifting in the field and in a wind tunnel[J]. Scientific Agriculture,1939,19(5):249 − 257.
[7] CHEPIL W S. Measurement of wind erosiveness of soils by dry sieving procedure[J]. Scientific Agriculture,1942,23(3):154 − 160.
[8] CHEPIL W S. Factors that influence clod structure and erodibility of soil by wind[J]. Soil Science,1953,75(6):473 − 484.
[9] CHEPIL W S. Factors that influence clod structure and erodibility of soil by wind[J]. Soil Science,1954,77(6):473 − 480.
[10] CHEPIL W S. Influence of moisture on erodibility of soil by wind[J]. Soil Science Society of America Journal,1956,20(2):288 − 292.
[11] 江凌,肖燚,欧阳志云,等. 基于RWEQ模型的青海省土壤风蚀模数估算[J]. 水土保持研究,2015,22(1):21 − 25. [JIANG Ling,XIAO Yi,OUYANG Zhiyun,et al. Estimate of the wind erosion modules in Qinghai Province based on RWEQ model[J]. Research of Soil and Water Conservation,2015,22(1):21 − 25. (in Chinese with English abstract)] JIANG Ling, XIAO Yi, OUYANG Zhiyun, et al. Estimate of the wind erosion modules in Qinghai Province based on RWEQ model［J］. Research of Soil and Water Conservation, 2015, 22（1）: 21 − 25. （in Chinese with English abstract）
[12] 殷海军. RWEQ模型的改进及其在风沙区的应用研究[D]. 北京:中国农业大学,2007. [YIN Haijun. Research on improvement and application of revised wind erosion model in the windy and sandy regions [D]. Beijing:China Agricultural University,2007. (in Chinese with English abstract)] YIN Haijun. Research on improvement and application of revised wind erosion model in the windy and sandy regions ［D］. Beijing: China Agricultural University, 2007. （in Chinese with English abstract）
[13] GREGORY J M,BORRELLI J,FEDLER C B. TEAM:Texas erosion analysis model[J]. American Society of Agricultural Engineers (Microfiche collection)(USA),1988.
[14] SHAO Y P,RAUPACH M R,LEYS J F. A model for predicting aeolian sand drift and dust entrainment on scales from paddock to region[J]. Soil Research,1996,34(3):309.
[15] 董治宝. 建立小流域风蚀量统计模型初探[J]. 水土保持通报,1998,18(5):55 − 62. [DONG Zhibao. Establishing statistic model of wind erosion on small watershed basis[J]. Bulletin of Soil and Water Conservation,1998,18(5):55 − 62. (in Chinese)] DOI: 10.3969/j.issn.1000-288X.1998.05.013 DONG Zhibao. Establishing statistic model of wind erosion on small watershed basis［J］. Bulletin of Soil and Water Conservation, 1998, 18（5）: 55 − 62. （in Chinese） DOI: 10.3969/j.issn.1000-288X.1998.05.013
[16] 严平,董光荣,张信宝等. 137Cs 法测定青藏高原土壤风 蚀 的 初 步 结 果[J]. 科 学 通 报,2000(2):199 − 204. (in Chinese). [17] 张国平. 基于遥感和 GIS 的中国土壤风力侵蚀研究[D]. 中国科学院研究生院(遥感应用研究所),2002. [ZHANG Remote sensing and GIS based study on the soil wind erosion in China[D]. Graduate University of Chinese Academy of Sciences (Institute of Remote Sensing Digital Earth),2002. (in Chinese with English abstract)] ZHANG Remote sensing and GIS based study on the soil wind erosion in China［D］. Graduate University of Chinese Academy of Sciences （Institute of Remote Sensing Digital Earth）, 2002. （in Chinese with English abstract）
[18] 程江浩. 1990-2020年青藏高原生态系统调节服务时空变化及其驱动力研究[D]. 大理大学,2023. [Cheng Jianghao. Spatial and temporal changes in typical ecosystem services and the driving forces,Tibetan Plateau,1990-2020[D]. Dali University,2023. (in Chinese with English abstract)] Cheng Jianghao. Spatial and temporal changes in typical ecosystem services and the driving forces, Tibetan Plateau, 1990-2020［D］. Dali University, 2023. （in Chinese with English abstract）
[19] 赖锋,乔占明连. 青海省风蚀量及防风固沙量时空特征分析[J]. 测绘学,2023,48(1):148 − 156. [LAI Feng,QIAO Zhanming,XIONG Zenglian. Spatiotemporal characteristics of wind erosion and wind prevention and sand fixation in Qinghai Province[J]. Science of Surveying and Mapping,2023,48(1):148 − 156. (in Chinese with English abstract)] LAI Feng, QIAO Zhanming, XIONG Zenglian. Spatiotemporal characteristics of wind erosion and wind prevention and sand fixation in Qinghai Province［J］. Science of Surveying and Mapping, 2023, 48（1）: 148 − 156. （in Chinese with English abstract）
[20] 巩俐,王发科,李积芳,等. 五道梁地区土壤风蚀变化特征及气象影响因子分析[J]. 青海草业,2020,29(3):59 − 63. [GONG Li,WANG Fake,LI Jifang,et al. Characteristics of soil wind erosion and analysis of meteorological influencing factors in Wudaoliang area[J]. Qinghai Prataculture,2020,29(3):59 − 63. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1008-1445.2020.03.015 GONG Li, WANG Fake, LI Jifang, et al. Characteristics of soil wind erosion and analysis of meteorological influencing factors in Wudaoliang area［J］. Qinghai Prataculture, 2020, 29（3）: 59 − 63. （in Chinese with English abstract） DOI: 10.3969/j.issn.1008-1445.2020.03.015
[21] 杨栋,王全成,姜昭群. 季节性冻土区滑坡防治工程的冻融效应分析[J]. 中国地质灾害与防治学报,2021,32(6):82 − 89. [YANG Dong,WANG Quancheng,JIANG Zhaoqun. Analysis on the effect of freeze-thaw on landslide prevention projects in seasonal frozen soil area[J]. The Chinese Journal of Geological Hazard and Control,2021,32(6):82 − 89. (in Chinese with English abstract)] YANG Dong, WANG Quancheng, JIANG Zhaoqun. Analysis on the effect of freeze-thaw on landslide prevention projects in seasonal frozen soil area［J］. The Chinese Journal of Geological Hazard and Control, 2021, 32（6）: 82 − 89. （in Chinese with English abstract）
[22] 沈凌铠,周保,魏刚,等. 气温变化对多年冻土斜坡稳定性的影响——以青海省浅层冻土滑坡为例[J]. 中国地质灾害与防治学报,2023,34(1):8 − 16. [SHEN Lingkai,ZHOU Bao,WEI Gang,et al. Influence of air temperature change on stability of permafrost slope:A case study of shallow permafrost landslide in Qinghai Province[J]. The Chinese Journal of Geological Hazard and Control,2023,34(1):8 − 16. (in Chinese with English abstract)] SHEN Lingkai, ZHOU Bao, WEI Gang, et al. Influence of air temperature change on stability of permafrost slope: A case study of shallow permafrost landslide in Qinghai Province［J］. The Chinese Journal of Geological Hazard and Control, 2023, 34（1）: 8 − 16. （in Chinese with English abstract）
[23] ZHAO Dongsheng,WU Shaohong. Projected changes in permafrost active layer thickness over the Qinghai-Tibet Plateau under climate change[J]. Water Resources Research,2019,55(9):7860 − 7875.
[24] PANG Qiangqiang,ZHAO Lin,LI Shuxun,et al. Active layer thickness variations on the Qinghai–Tibet Plateau under the scenarios of climate change[J]. Environmental Earth Sciences,2012,66(3):849 − 857.
[25] 庞强强,李述训,吴通华,等. 青藏高原冻土区活动层厚度分布模拟[J]. 冰川冻土,2006,28(3):390 − 395. [PANG Qiangqiang,LI Shuxun,WU Tonghua,et al. Simulated distribution of active layer depths in the frozen ground regions of Tibetan Plateau[J]. Journal of Glaciology and Geocryology,2006,28(3):390 − 395. (in Chinese with English abstract)] PANG Qiangqiang, LI Shuxun, WU Tonghua, et al. Simulated distribution of active layer depths in the frozen ground regions of Tibetan Plateau［J］. Journal of Glaciology and Geocryology, 2006, 28（3）: 390 − 395. （in Chinese with English abstract）
[26] 徐晓明,吴青柏,张中琼. 青藏高原多年冻土活动层厚度对气候变化的响应[J]. 冰川冻土,2017,39(1):1 − 8. [XU Xiaoming,WU Qingbai,ZHANG Zhongqiong. Responses of active layer thickness on the Qinghai-Tibet Plateau to climate change[J]. Journal of Glaciology and Geocryology,2017,39(1):1 − 8. (in Chinese with English abstract)] XU Xiaoming, WU Qingbai, ZHANG Zhongqiong. Responses of active layer thickness on the Qinghai-Tibet Plateau to climate change［J］. Journal of Glaciology and Geocryology, 2017, 39（1）: 1 − 8. （in Chinese with English abstract）
[27] 徐晓明,吴青柏. 三江源多年冻土活动层厚度变化特征研究[J/OL]. 冰川冻土:1 − 15. [XU Xiaoming,WU Qingbai. Research on the variation characteristics of active layer thickness of permafrost in the Three River Source Region[J/OL]. Glacial permafrost:1 − 15. (in Chinese with English abstract)] XU Xiaoming, WU Qingbai. Research on the variation characteristics of active layer thickness of permafrost in the Three River Source Region［J/OL］. Glacial permafrost: 1 − 15. （in Chinese with English abstract）
[28] 张中琼,吴青柏. 气候变化情景下青藏高原多年冻土活动层厚度变化预测[J]. 冰川冻土,2012,34(3):505 − 511. [ZHANG Zhongqiong,WU Qingbai. Predicting changes of active layer thickness on the Qinghai-Tibet Plateau as climate warming[J]. Journal of Glaciology and Geocryology,2012,34(3):505 − 511. (in Chinese with English abstract)] ZHANG Zhongqiong, WU Qingbai. Predicting changes of active layer thickness on the Qinghai-Tibet Plateau as climate warming［J］. Journal of Glaciology and Geocryology, 2012, 34（3）: 505 − 511. （in Chinese with English abstract）
[29] 刘文惠,谢昌卫,刘海瑞,等. Stefan方程在土壤冻融过程模拟中的应用[J]. 冰川冻土,2022,44(1):327 − 339. [LIU Wenhui,XIE Changwei,LIU Hairui,et al. Application of Stefan equation on simulating freezing-thawing process of permafrost[J]. Journal of Glaciology and Geocryology,2022,44(1):327 − 339. (in Chinese with English abstract)] LIU Wenhui, XIE Changwei, LIU Hairui, et al. Application of Stefan equation on simulating freezing-thawing process of permafrost［J］. Journal of Glaciology and Geocryology, 2022, 44（1）: 327 − 339. （in Chinese with English abstract）
[30] 张云霞. 2000-2020年青藏高原生态质量时空变化及其影响因素[D]. 兰州:兰州大学,2023. [ZHANG Yunxia. Spatiotemporal changes of ecological quality and itsinfluencing factors of the Qinghai-Tibetan Plateaufrom 2000 to 2020 [D]. Lanzhou:Lanzhou University,2023. (in Chinese with English abstract)] ZHANG Yunxia. Spatiotemporal changes of ecological quality and itsinfluencing factors of the Qinghai-Tibetan Plateaufrom 2000 to 2020 ［D］. Lanzhou: Lanzhou University, 2023. （in Chinese with English abstract）
[31] NI Jie,WU Tonghua,ZHU Xiaofan,et al. Simulation of the present and future projection of permafrost on the Qinghai-Tibet Plateau with statistical and machine learning models[J]. Journal of Geophysical Research:Atmospheres,2021,126(2):e2020jd033402.
[32] 高荣,韦志刚,董文杰. 青藏高原土壤冻结始日和终日的年际变化[J]. 冰川冻土,2003,25(1):49 − 54. [GAO Rong,WEI Zhigang,DONG Wenjie. Interannual variation of the beginning date and the ending date of soil freezing in the Tibetan Plateau[J]. Journal of Glaciology and Geocryology,2003,25(1):49 − 54. (in Chinese with English abstract)] GAO Rong, WEI Zhigang, DONG Wenjie. Interannual variation of the beginning date and the ending date of soil freezing in the Tibetan Plateau［J］. Journal of Glaciology and Geocryology, 2003, 25（1）: 49 − 54. （in Chinese with English abstract）
[33] 王一菲,郑粉莉,张加琼,等. 冻融作用对典型黑土土壤风蚀的影响[J]. 水土保持学报,2020,34(5):34 − 41. [WANG Yifei,ZHENG Fenli,ZHANG Jiaqiong,et al. Effects of freeze-thaw on wind erosion in the typical mollisol region[J]. Journal of Soil and Water Conservation,2020,34(5):34 − 41. (in Chinese with English abstract)] WANG Yifei, ZHENG Fenli, ZHANG Jiaqiong, et al. Effects of freeze-thaw on wind erosion in the typical mollisol region［J］. Journal of Soil and Water Conservation, 2020, 34（5）: 34 − 41. （in Chinese with English abstract）
[34] 董瑞琨,许兆义,杨成永. 青藏高原冻融侵蚀动力特征研究[J]. 水土保持学报,2000,14(4):12 − 16. [DONG Ruikun,XU Zhaoyi,YANG Chengyong. Dynamic and characteristic of freezing thawing erosion on Qinghai Tibet Plateau[J]. Journal of Soil Water Conservation,2000,14(4):12 − 16. (in Chinese with English abstract)] DOI: 10.3321/j.issn:1009-2242.2000.04.003 DONG Ruikun, XU Zhaoyi, YANG Chengyong. Dynamic and characteristic of freezing thawing erosion on Qinghai Tibet Plateau［J］. Journal of Soil Water Conservation, 2000, 14（4）: 12 − 16. （in Chinese with English abstract） DOI: 10.3321/j.issn:1009-2242.2000.04.003
[35] 吴雅琼,刘国华,傅伯杰,等. 青藏高原土壤有机碳密度垂直分布研究[J]. 环境科学学报,2008,28(2):362 − 367. [WU Yaqiong,LIU Guohua,FU Bojie,et al. Study on the vertical distribution of soil organic carbon density in the Tibetan Plateau[J]. Acta Scientiae Circumstantiae,2008,28(2):362 − 367. (in Chinese with English abstract)] WU Yaqiong, LIU Guohua, FU Bojie, et al. Study on the vertical distribution of soil organic carbon density in the Tibetan Plateau［J］. Acta Scientiae Circumstantiae, 2008, 28（2）: 362 − 367. （in Chinese with English abstract）
[36] 刘勇,魏良帅,黄安邦,等. 气候变化下长江源土壤水时空演化及其环境响应[J]. 水文地质工程地质,2023,50(5):39 − 52. [LIU Yong,WEI Liangshuai,HUANG Anbang,et al. Spatial and temporal evolution of soil water and its response to the environment in the Yangtze River source area under climate change[J]. Hydrogeology & Engineering Geology,2023,50(5):39 − 52. (in Chinese with English abstract)] LIU Yong, WEI Liangshuai, HUANG Anbang, et al. Spatial and temporal evolution of soil water and its response to the environment in the Yangtze River source area under climate change［J］. Hydrogeology & Engineering Geology, 2023, 50（5）: 39 − 52. （in Chinese with English abstract）
[37] 张永双,郭长宝,李向全,等. 川藏铁路廊道关键水工环地质问题:现状与发展方向[J]. 水文地质工程地质,2021,48(5):1 − 12. [ZHANG Yongshuang,GUO Changbao,LI Xiangquan,et al. Key problems on hydro-engineering-environmental geology along the Sichuan-Tibet Railway corridor:Current status and development direction[J]. Hydrogeology & Engineering Geology,2021,48(5):1 − 12. (in Chinese with English abstract)] ZHANG Yongshuang, GUO Changbao, LI Xiangquan, et al. Key problems on hydro-engineering-environmental geology along the Sichuan-Tibet Railway corridor: Current status and development direction［J］. Hydrogeology & Engineering Geology, 2021, 48（5）: 1 − 12. （in Chinese with English abstract）
[38] 支泽民,刘峰贵,周强,等. 基于流域单元的地质灾害易发性评价——以西藏昌都市为例[J]. 中国地质灾害与防治学报,2023,34(1):139 − 150. [ZHI Zemin,LIU Fenggui,ZHOU Qiang,et al. Evaluation of geological hazards susceptibility based on watershed units:A case study of the Changdu City,Tibet[J]. The Chinese Journal of Geological Hazard and Control,2023,34(1):139 − 150. (in Chinese with English abstract)] ZHI Zemin, LIU Fenggui, ZHOU Qiang, et al. Evaluation of geological hazards susceptibility based on watershed units: A case study of the Changdu City, Tibet［J］. The Chinese Journal of Geological Hazard and Control, 2023, 34（1）: 139 − 150. （in Chinese with English abstract）
计量
- 文章访问数: 84
- HTML全文浏览量: 57
- PDF下载量: 0
邮件订阅
RSS