古尔班通古特沙漠南缘土壤：空间异质性与尺度特征的深度剖析

古尔班通古特沙漠南缘土壤：空间异质性与尺度特征的深度剖析一、引言1.1研究背景古尔班通古特沙漠作为中国第二大沙漠，位于新疆准噶尔盆地中央，其南缘处于天山北麓，地理位置独特且重要，处于干旱区与绿洲的过渡地带，是生态环境脆弱敏感区域。在该区域内，土壤是生态系统的关键组成部分，不仅承载着植被生长，还在物质循环和能量流动中扮演着重要角色。然而，由于特殊的地理位置和复杂的气候条件，古尔班通古特沙漠南缘土壤的空间异质性明显，不同区域土壤的物理、化学和生物学性质存在显著差异。土壤空间异质性是指土壤特性在空间分布上的不均匀性和变异性。在古尔班通古特沙漠南缘，这种异质性体现在多个方面。从土壤质地来看，不同地段的土壤颗粒组成差异较大，有的区域以砂粒为主，土壤通气性良好但保水保肥能力差；有的区域则含有一定比例的黏粒和粉粒，土壤的结构和肥力状况有所不同。在土壤养分方面，有机质、氮、磷、钾等含量在空间上分布不均，呈现出斑块状或梯度变化。这些土壤特性的空间变化并非随机，而是受到多种因素的综合影响。地形地貌是重要的影响因素之一，沙丘、丘间低地等不同地貌部位，土壤的堆积和侵蚀过程不同，导致土壤性质产生差异。植被类型和覆盖度也与土壤空间异质性密切相关，不同植被通过根系活动、凋落物归还等方式，对土壤的物质组成和生态过程产生不同作用。尺度特征是理解土壤空间异质性的重要维度。在不同尺度下，土壤空间异质性的表现和影响因素各不相同。在较小尺度上，如植物个体根际附近，根系分泌物、微生物活动等生物因素对土壤性质的影响较为突出，可能导致土壤养分、酸碱度等在短距离内发生明显变化。随着尺度增大到种群、群落乃至区域尺度，地形、气候等宏观因素的作用逐渐增强，它们在更大范围内控制着土壤的形成和演化过程，使得土壤空间异质性呈现出更为复杂的格局。例如，在区域尺度上，降水、温度等气候因子的空间差异，会导致土壤水分和热量条件不同，进而影响土壤的风化程度、有机质分解速率等，最终造成土壤性质在区域范围内的分异。土壤的空间异质性及其尺度特征对古尔班通古特沙漠南缘的生态系统和土地利用具有深远影响。在生态系统方面，土壤异质性影响着植被的分布、生长和多样性。土壤养分丰富、水分条件适宜的区域，往往植被生长茂盛，物种丰富度较高；而土壤贫瘠、干旱的区域，植被则较为稀疏，物种组成相对单一。这种植被分布格局又反过来影响土壤的侵蚀、水分涵养和生物地球化学循环等生态过程，形成土壤与植被之间的相互作用和反馈机制。从土地利用角度来看，了解土壤空间异质性及其尺度特征，对于合理规划农业生产、生态保护和土地整治至关重要。在农业种植中，根据土壤肥力的空间分布差异，可以精准地进行施肥、灌溉等农事操作，提高资源利用效率，减少农业面源污染。在生态保护方面，明确土壤异质性的关键区域和尺度，可以有针对性地制定生态修复和保护策略，保护脆弱的生态环境。1.2研究目的与意义本研究旨在深入剖析古尔班通古特沙漠南缘土壤的空间异质性及其尺度特征，具体目的如下：其一，精确测定该区域土壤的多种属性，包括土壤质地、养分含量、酸碱度、电导率等，全面掌握土壤的基础信息。其二，运用地理信息系统（GIS）和地统计学等方法，分析土壤空间异质性特征，构建空间变异结构模型，直观展现土壤特性在空间上的分布和变化规律。其三，深入研究不同尺度下土壤空间异质性的变化规律，明确各尺度上影响土壤异质性的主导因素，以及不同尺度间土壤空间异质性的内在联系。其四，探究土壤空间异质性与地形、气候、植被、土地利用等环境因素的相互关系，揭示土壤空间异质性形成的内在机制。本研究具有重要的理论与实践意义。从理论层面来看，有助于深化对干旱区土壤形成与演化过程的理解，丰富土壤地理学和生态学的理论体系。通过研究不同尺度下土壤空间异质性，进一步完善生态系统等级结构理论，为理解生态过程在不同尺度上的表现提供依据。从实践意义出发，对古尔班通古特沙漠南缘的土地利用规划和生态环境管理具有重要的指导作用。在土地利用规划方面，依据土壤空间异质性特征，可以合理划分土地利用类型，优化农业布局。对于土壤肥力较高、水分条件较好的区域，可规划为优质农田，进行高效农业生产；而对于土壤贫瘠、生态脆弱的区域，则适宜发展生态林业或畜牧业，避免过度开垦导致土地退化。在生态环境管理方面，明确土壤空间异质性与生态环境的关系，有助于制定精准的生态保护和修复策略。针对土壤侵蚀严重、生态功能薄弱的区域，采取有效的水土保持措施，如植树造林、种草固沙等，促进生态系统的恢复和稳定。此外，研究成果还能为区域水资源合理利用、农业面源污染防治等提供科学依据，推动古尔班通古特沙漠南缘地区的可持续发展。1.3国内外研究现状土壤空间异质性及其尺度特征一直是土壤学、生态学等领域的研究热点。国外在该领域的研究起步较早，20世纪中叶，随着地统计学的兴起，为土壤空间异质性研究提供了有力工具。如法国著名土壤学家Matheron首次将地统计学引入土壤科学研究，提出了区域化变量理论和半方差函数分析方法，开启了定量研究土壤空间异质性的先河。此后，众多学者运用地统计学方法对不同区域、不同类型土壤的空间异质性展开深入研究。在尺度特征研究方面，国外学者强调多尺度研究的重要性，通过构建不同尺度的土壤空间变异模型，分析尺度效应及其对生态过程的影响。如在北美草原地区的研究中，学者们发现土壤养分的空间异质性在小尺度上主要受植物根系和微生物活动影响，而在大尺度上则与气候、地形等因素密切相关。在森林生态系统中，研究表明土壤碳氮等化学性质的空间异质性在不同尺度上表现出不同的主导因素，且尺度间存在复杂的相互作用。国内对土壤空间异质性及其尺度特征的研究始于20世纪80年代，随着地理信息系统（GIS）和全球定位系统（GPS）技术的广泛应用，研究水平得到了快速提升。在干旱半干旱地区，众多学者针对不同沙漠和沙地开展研究。例如，在腾格里沙漠的研究中，通过对土壤水分、养分等指标的测定和地统计学分析，揭示了土壤空间异质性与沙丘地貌、植被分布之间的关系。在毛乌素沙地，研究发现土壤理化性质的空间异质性在不同植被恢复阶段存在明显差异，植被恢复过程对土壤空间异质性的形成和演变具有重要影响。在尺度特征研究方面，国内学者结合不同生态系统类型，探讨了土壤空间异质性在不同尺度上的变化规律。如在黄土高原地区，研究了不同尺度下土壤侵蚀与土壤空间异质性的关系，发现随着尺度增大，地形因素对土壤侵蚀和土壤空间异质性的影响逐渐增强。在湿地生态系统中，研究了土壤碳、氮、磷等养分在不同尺度上的空间分布特征，以及尺度效应与湿地生态功能的关系。然而，针对古尔班通古特沙漠南缘土壤的空间异质性及其尺度特征的研究相对较少。虽然已有一些关于古尔班通古特沙漠土壤的研究，但多集中在土壤化学性质、微生物群落等方面的整体分析，对空间异质性及其尺度特征的系统研究尚显不足。在空间异质性研究中，缺乏对土壤多种属性在不同地貌类型、植被覆盖条件下的全面分析，难以准确揭示土壤空间异质性的形成机制和分布规律。在尺度特征研究方面，对不同尺度下土壤空间异质性的变化规律和尺度转换关系研究不够深入，尚未建立起完善的尺度效应模型。此外，将土壤空间异质性及其尺度特征与区域生态系统功能、土地利用变化等相结合的研究也较为薄弱，限制了研究成果在区域生态环境保护和土地资源管理中的应用。二、研究区域与方法2.1研究区概况古尔班通古特沙漠南缘位于新疆维吾尔自治区准噶尔盆地南部，处于北纬43°50′-45°20′，东经84°50′-89°10′之间，南邻天山北麓，北接古尔班通古特沙漠主体。该区域是典型的干旱区与绿洲过渡地带，地理位置独特，生态环境脆弱且敏感，在区域生态系统和土地利用中具有重要地位。古尔班通古特沙漠南缘属于温带大陆性干旱气候，其显著特点是降水稀少、蒸发量大。年降水量通常在100-200毫米之间，且降水分布极为不均，主要集中在春季和夏季，其中夏季降水相对较多，但多以暴雨形式出现，难以有效被土壤储存和植被利用。年蒸发量高达2000-3000毫米，远远超过降水量，导致土壤水分长期处于亏缺状态。该地区光照充足，年日照时数可达2800-3200小时，为植物的光合作用提供了有利条件。气温年较差和日较差都很大，年平均气温在6-8℃之间，夏季炎热，极端最高气温可达40℃以上；冬季寒冷，极端最低气温可降至-30℃以下，日较差一般在10-15℃之间，这种较大的温差有利于植物体内糖分的积累。研究区内地形地貌复杂多样，主要包括山前冲积扇、洪积平原、沙漠边缘沙丘和丘间低地等。山前冲积扇是由河流出山口后，因坡度骤减，流速降低，携带的大量碎屑物质堆积而成，地势较高，地形较为平坦，土壤颗粒较粗，多为砾石和粗砂。洪积平原是洪水携带的泥沙等物质在山前堆积形成，地形相对平坦，土壤质地以砂壤土和壤土为主。沙漠边缘沙丘形态各异，主要有新月形沙丘、纵向沙垄等，沙丘高度一般在5-20米之间，沙丘上的土壤以风沙土为主，质地疏松，保水保肥能力差。丘间低地地势较低，是沙丘之间的低洼区域，土壤水分条件相对较好，部分区域有季节性积水，土壤类型多为盐土和草甸土。植被类型以荒漠植被为主，植物种类相对较少，群落结构简单。主要植被物种包括梭梭（Haloxylonammodendron）、红柳（Tamarixramosissima）、琵琶柴（Reaumuriasoongorica）、芨芨草（Achnatherumsplendens）等。梭梭是该区域的优势物种之一，具有耐旱、抗风沙的特性，常形成大面积的梭梭林，对于固定沙丘、保持水土具有重要作用。红柳多生长在丘间低地和水分条件相对较好的区域，其根系发达，能够适应干旱和盐碱环境，具有良好的防风固沙能力。琵琶柴分布较为广泛，对干旱和贫瘠的土壤有较强的适应性。芨芨草是一种多年生草本植物，主要生长在盐渍化较轻的低地和河滩地带，是牲畜的重要饲料。在绿洲边缘，有人工种植的杨树（Populusspp.）、柳树（Salixspp.）等防护林树种，以及小麦（Triticumaestivum）、玉米（Zeamays）、棉花（Gossypiumhirsutum）等农作物。土地利用现状主要包括绿洲农业用地、荒漠生态用地和少量的城镇建设用地。绿洲农业用地主要分布在山前冲积扇和洪积平原的水源充足区域，土壤肥沃，灌溉条件良好，是该区域的主要农业生产基地，种植作物以小麦、玉米、棉花等为主。随着农业的发展，部分地区存在过度开垦和不合理灌溉的现象，导致土壤盐渍化和土地退化问题日益严重。荒漠生态用地面积广阔，主要包括沙漠、戈壁和自然植被覆盖区域，是维持区域生态平衡的重要组成部分，但由于长期受到风沙侵蚀和人类活动的干扰，生态系统较为脆弱。城镇建设用地集中在交通便利、水源充足的区域，近年来，随着区域经济的发展，城镇规模不断扩大，对周边土地资源的占用和生态环境的影响也逐渐加剧。2.2研究方法2.2.1土壤采样在古尔班通古特沙漠南缘研究区域内，采用系统抽样法确定土壤采样点位置。首先，利用ArcGIS软件将研究区域划分为500m×500m的网格，以每个网格的中心点作为潜在采样点。借助GPS定位仪，精确记录每个采样点的经纬度信息，确保采样点位置的准确性。在实际采样过程中，充分考虑地形地貌、植被覆盖等因素对土壤特性的影响，对部分采样点位置进行适当调整，以保证所采集的土壤样品能够全面、准确地代表研究区域内不同类型的土壤。例如，在沙丘顶部、丘间低地、绿洲边缘等具有明显地形和植被差异的区域，增加采样点的密度；而在地形相对平坦、植被类型较为单一的区域，适当减少采样点数量。最终，共确定了200个土壤采样点。在每个采样点，使用不锈钢土壤采样器采集不同层次的土壤样品。对于表层土壤（0-20cm），将采样器垂直插入土壤，然后轻轻旋转，使土壤进入采样器内，取出采样器后，将土壤样品放入已编号的密封塑料袋中。对于深层土壤（20-50cm），采用分层采样的方法，每隔10cm采集一个样品。在采集深层土壤样品时，先用土钻钻出一个垂直的孔，然后将采样器插入孔中，按照上述方法采集土壤样品。为了保证样品的代表性，在每个采样点的周围1m范围内，随机采集3个次样品，将这3个次样品混合均匀后，作为该采样点的土壤样品。在完成土壤样品采集后，立即用记号笔在密封塑料袋上标记采样点编号、采样日期、采样深度等信息，并将样品放入冷藏箱中，尽快带回实验室进行处理。回到实验室后，将土壤样品平铺在干净的塑料布上，置于通风良好的室内自然风干，风干过程中，定期翻动土壤样品，以加速风干速度，并防止土壤样品发霉变质。待土壤样品完全风干后，用木棍将其轻轻碾碎，去除其中的植物根系、石块等杂物，然后过2mm筛，将过筛后的土壤样品装入密封袋中，保存备用。2.2.2土壤理化性质分析采用电位法测定土壤pH值，称取10g过2mm筛的风干土样置于250mL烧杯中，按照土水比1:2.5的比例加入去离子水，用玻璃棒搅拌均匀，使土样充分分散，静置30min，待土壤悬浊液稳定后，使用校正后的pH计测定悬液的pH值，测定时将玻璃电极球部浸入悬液泥层中，甘汞电极浸在悬液上部清液中，读取pH值。土壤电导率（EC）使用DDS-307A型电导率仪测定，同样称取10g风干土样，按照土水比1:5的比例加入去离子水，搅拌均匀后，静置1h，然后将电导率仪的电极插入上清液中，测定电导率，测定结果以mS/cm表示。土壤有机质含量采用重铬酸钾容量法-稀释热法测定。准确称取0.2-0.5g过0.25mm筛的风干土样于硬质试管中，加入10mL0.8mol/L（1/6K₂Cr₂O₇）溶液，然后缓慢加入20mL浓硫酸，在试管口加一小漏斗，将试管放入已预热至170-180℃的油浴锅中，煮沸5min，取出试管，稍冷后，将试管内容物小心转移至250mL三角瓶中，用去离子水冲洗试管和漏斗，使三角瓶内溶液总体积约为150mL，冷却后，加入3-4滴邻菲罗啉指示剂，用0.2mol/L硫酸亚铁标准溶液滴定，溶液由黄色经过绿色、淡绿色突变至棕红色即为终点。同时做两个空白试验，取其平均值。有机质含量计算公式为：有机质（g/kg）=[((V₀-V)×N×0.003×1.1×1.724)/样品重]×1000，其中V₀为滴定空白液时所用去的硫酸亚铁毫升数，V为滴定样品液时所用去的硫酸亚铁毫升数，N为标准硫酸亚铁的浓度。土壤全氮含量测定采用凯氏定氮法。称取0.5-1.0g风干土样于消化管中，加入1g混合加速剂（K₂SO₄:CuSO₄:Se=100:10:1）和5mL浓硫酸，将消化管置于消化炉上，先低温加热10-15min，待管内反应缓和后，加强火力使消煮的土液保持微沸，消煮至溶液和土粒全部变为灰白稍带绿色后，再继续消煮1h。消煮完毕，待消化液稍冷，加入约50mL去离子水，冷却后转移至100mL容量瓶中，定容摇匀。取10mL消化液于蒸馏装置中，加入20mL10mol/LNaOH溶液，蒸馏出的氨用2%硼酸溶液吸收，以0.02mol/L硫酸标准溶液滴定，溶液由蓝绿色变为紫红色即为终点。土壤全氮含量计算公式为：土壤全氮N（g/kg）=[(V-V₀)×c×14.0×10⁻³]/m，其中V为滴定样品时所用的酸标准液的体积，V₀为滴定空白时所用的酸标准液的体积，c为酸标准液的浓度，m为烘干土的质量。土壤速效氮采用碱解扩散法测定。称取2.00g风干土样于扩散皿外室，加入0.20gFeSO₄・7H₂O粉末，轻轻转动扩散皿，使土样与FeSO₄・7H₂O粉末均匀铺平，在扩散皿外室边缘涂抹碱性胶液，内室加入2mL2%硼酸-指示剂溶液，然后迅速加入10mL1.07mol/LNaOH溶液，立即盖严毛玻璃，用橡皮筋固定，将扩散皿置于40±1℃恒温箱中，碱解扩散24±0.5h后取出，用0.01mol/L硫酸标准溶液滴定内室吸收液中的氨，溶液由蓝绿色变为紫红色即为终点。碱解氮含量计算公式为：碱解氮N（mg/kg）=[(V-V₀)×c×14.0×10³]/m，其中V为样品滴定时用去硫酸标准液的体积，V₀为空白试验滴定时用去硫酸标准液的体积，c为硫酸标准液的浓度，m为烘干土的质量。土壤全磷含量采用氢氧化钠熔融-钼锑抗比色法测定。称取0.5g风干土样于镍坩埚中，加入2gNaOH，将坩埚放入高温炉中，从低温升至720℃，保持15min，取出冷却后，将坩埚放入250mL烧杯中，用热水浸取，洗出坩埚，将溶液转移至100mL容量瓶中，定容摇匀。吸取5mL上清液于50mL容量瓶中，加入5mL钼锑抗显色剂，定容摇匀，放置30min后，用分光光度计在波长700nm处比色测定。土壤全磷含量计算公式为：土壤全磷（g/kg）=（显色液磷浓度×显色液体积×分取倍数×10⁻⁶）/样品质量。土壤速效磷采用0.5mol/LNaHCO₃浸提-钼锑抗比色法测定。称取1.0g风干土样于50mL三角瓶中，加入20mL0.5mol/LNaHCO₃溶液，振荡30min后过滤，吸取5mL滤液于50mL容量瓶中，加入5mL钼锑抗显色剂，定容摇匀，放置30min后，用分光光度计在波长700nm处比色测定。土壤速效磷含量计算公式为：土壤速效磷P（mg/kg）=（显色液磷浓度×显色液体积×分取倍数×10⁻³）/样品质量。土壤颗粒组成采用激光粒度分析仪测定。称取10g风干土样于500mL烧杯中，加入50mL0.5mol/LNa₂CO₃溶液，用玻璃棒搅拌均匀，浸泡过夜，使土壤颗粒充分分散。将分散后的土壤悬浊液转移至激光粒度分析仪的样品池中，按照仪器操作规程进行测定，仪器自动分析并输出土壤颗粒组成数据，包括砂粒（2-0.05mm）、粉粒（0.05-0.002mm）和黏粒（<0.002mm）的含量。2.2.3空间分析方法运用地统计学方法对土壤理化性质的空间异质性特征进行分析。地统计学是以区域化变量理论为基础，以变异函数为主要工具，研究自然现象空间分布特征和变异规律的一门学科。在本研究中，将土壤的各项理化性质视为区域化变量，利用半方差函数来描述其空间变异特征。半方差函数的计算公式为：\gamma(h)=\frac{1}{2N(h)}\sum_{i=1}^{N(h)}[Z(x_i)-Z(x_i+h)]^2其中，\gamma(h)为半方差值，h为空间滞后距离（步长），N(h)为间距为h的样本对数，Z(x_i)和Z(x_i+h)分别为区域化变量Z(x)在空间位置x_i和x_i+h处的观测值。通过计算不同步长h下的半方差值\gamma(h)，绘制半方差函数曲线图。根据半方差函数曲线的形状和特征，选择合适的理论模型进行拟合，常用的理论模型包括球状模型、指数模型、高斯模型等。拟合得到的理论模型参数包括块金值（Nugget）、基台值（Sill）和变程（Range）。块金值表示随机因素引起的空间变异，反映了小尺度上的空间异质性；基台值是半方差函数在变程处的取值，代表了区域化变量的总空间变异程度；变程表示区域化变量在空间上的自相关距离，在变程范围内，变量具有空间相关性，超过变程，空间相关性消失。利用克里格插值法对土壤理化性质进行空间插值，绘制空间分布图。克里格插值是一种基于半方差函数模型的最优内插法，它在考虑观测点与待估点之间空间位置关系和变量空间相关性的基础上，对待估点的值进行无偏最优估计。在ArcGIS软件中，选择地统计分析模块，根据拟合得到的半方差函数理论模型，进行克里格插值操作，生成土壤pH值、电导率、有机质、全氮、速效氮、全磷、速效磷等理化性质的空间分布图。通过这些空间分布图，可以直观地了解土壤各项理化性质在研究区域内的空间分布格局和变化趋势。2.2.4尺度分析方法利用分形分析法研究土壤空间异质性的尺度特征。分形理论是研究自然界中复杂现象和不规则图形的有力工具，它能够定量描述事物的自相似性和标度不变性。在土壤空间异质性研究中，分形维数可以作为衡量土壤空间异质性程度的指标，分形维数越大，表明土壤空间异质性越强。采用容量维数法计算土壤理化性质的分形维数。容量维数的计算方法如下：D=\frac{\lnN(r)}{\ln(1/r)}其中，D为分形维数，N(r)为尺度r下的土壤性质变异系数，r为尺度因子，可以是采样间距、研究区域面积等。在本研究中，通过改变采样间距，分别计算不同采样间距下土壤各项理化性质的变异系数N(r)，然后以\ln(1/r)为横坐标，\lnN(r)为纵坐标，进行线性回归分析，得到回归直线的斜率，即为分形维数D。通过分析不同尺度下土壤理化性质的分形维数变化规律，探究土壤空间异质性在不同尺度上的表现和变化特征。例如，研究随着采样间距的增大，土壤分形维数的变化趋势，以及不同土壤理化性质分形维数对尺度变化的响应差异。同时，结合地统计学分析结果，探讨不同尺度下土壤空间异质性的主导因素和形成机制，以及尺度效应与土壤空间变异结构之间的关系。三、古尔班通古特沙漠南缘土壤空间异质性分析3.1土壤理化性质的描述性统计对采集的200个土壤样品进行各项理化性质分析后，得到的描述性统计结果如表1所示。土壤pH值范围在7.85-8.96之间，平均值为8.32，表明研究区土壤整体呈碱性。最大值与最小值之间的差值为1.11，标准差为0.24，变异系数为2.89%，说明土壤pH值在研究区域内的离散程度较小，变异性较弱。这种相对稳定的pH值可能与该地区的母质类型、气候条件以及长期的成土过程有关。研究区气候干旱，降水稀少，土壤中的碱性物质难以被淋溶，使得土壤pH值维持在相对较高且稳定的水平。土壤电导率（EC）平均值为1.25mS/cm，最大值达3.56mS/cm，最小值为0.32mS/cm，标准差为0.56，变异系数为44.80%，属于中等变异性。电导率反映了土壤中可溶性盐分的含量，其较大的变异性表明研究区域内土壤盐分含量存在明显的空间差异。在沙漠边缘的一些低洼区域，由于地下水位较高，盐分随水分蒸发向上运移并积累，导致土壤电导率较高；而在沙丘顶部等地势较高、排水良好的区域，土壤盐分含量相对较低，电导率也较低。土壤有机质含量平均值为12.56g/kg，最大值为35.68g/kg，最小值仅为2.34g/kg，标准差为6.54，变异系数高达52.10%，表现出较强的变异性。土壤有机质主要来源于植物残体、根系分泌物以及土壤微生物的活动等。在植被覆盖较好的区域，如丘间低地生长有梭梭、红柳等植物，植物残体归还量大，土壤有机质含量相对较高；而在植被稀疏的沙丘区域，植物生长受限，土壤有机质来源少，含量较低。土壤全氮含量平均值为0.86g/kg，最大值为1.89g/kg，最小值为0.23g/kg，标准差为0.35，变异系数为40.70%，呈现中等变异性。全氮含量的空间差异与土壤有机质含量密切相关，同时也受到植被类型和土地利用方式的影响。绿洲农业用地中，由于施肥等农事活动，土壤全氮含量相对较高；而在荒漠生态用地中，自然植被生长条件相对较差，土壤全氮含量较低。土壤速效氮平均值为35.68mg/kg，最大值为89.56mg/kg，最小值为10.23mg/kg，标准差为18.56，变异系数为52.00%，变异性较强。速效氮是土壤中可被植物直接吸收利用的氮素形态，其含量受到土壤微生物活动、氮肥施用以及土壤水分等因素的影响。在水分条件较好、微生物活性较高的区域，土壤中有机氮的矿化作用较强，速效氮含量相对较高。土壤全磷含量平均值为0.65g/kg，最大值为1.23g/kg，最小值为0.25g/kg，标准差为0.21，变异系数为32.30%，属于中等变异性。土壤全磷含量主要取决于成土母质，同时也受到土壤侵蚀、淋溶等因素的影响。在山前冲积扇等母质富含磷素的区域，土壤全磷含量相对较高；而在长期遭受风沙侵蚀的沙丘区域，土壤中细颗粒物质被吹蚀，全磷含量相对较低。土壤速效磷平均值为10.56mg/kg，最大值为35.68mg/kg，最小值为2.34mg/kg，标准差为7.56，变异系数为71.60%，变异性很强。速效磷的含量除了与土壤全磷含量有关外，还受到土壤酸碱度、铁铝氧化物含量等因素的影响。在碱性土壤中，磷素容易与钙结合形成难溶性磷酸盐，降低速效磷含量；而在一些富含铁铝氧化物的区域，磷素会被吸附固定，也会导致速效磷含量较低。土壤砂粒含量平均值为75.68%，最大值为92.34%，最小值为52.10%，标准差为8.56，变异系数为11.30%，变异性较弱。研究区地处沙漠边缘，风沙活动频繁，土壤以风沙土为主，砂粒含量较高且相对稳定。土壤粉粒含量平均值为15.68%，最大值为32.10%，最小值为5.68%，标准差为5.56，变异系数为35.50%，属于中等变异性。粉粒含量的变化与土壤的搬运和沉积过程有关，在不同地貌部位和植被覆盖条件下，粉粒的沉积和流失情况不同，导致其含量存在一定差异。土壤黏粒含量平均值为8.64%，最大值为20.10%，最小值为2.34%，标准差为3.56，变异系数为41.20%，呈现中等变异性。黏粒含量的空间分布受到成土过程中物理、化学和生物作用的综合影响，在植被根系活动频繁、土壤团聚体稳定性较好的区域，黏粒含量相对较高。3.2土壤空间自相关性分析运用空间自相关分析方法，深入探究古尔班通古特沙漠南缘土壤理化性质在空间上的分布特征和相关性。空间自相关分析能够揭示土壤属性在不同位置之间的依赖关系，判断其分布是呈现随机状态还是具有某种空间格局。本研究采用全局Moran'sI指数来衡量土壤属性的空间自相关程度，其计算公式为：I=\frac{n}{\sum_{i=1}^{n}\sum_{j=1}^{n}w_{ij}}\frac{\sum_{i=1}^{n}\sum_{j=1}^{n}w_{ij}(x_i-\overline{x})(x_j-\overline{x})}{\sum_{i=1}^{n}(x_i-\overline{x})^2}其中，n为样本数量，x_i和x_j分别为第i个和第j个样本的属性值，\overline{x}为样本属性值的平均值，w_{ij}为空间权重矩阵，表示样本i和样本j之间的空间关系，若i和j相邻，则w_{ij}=1，否则w_{ij}=0。Moran'sI指数的取值范围为[-1,1]，当I>0时，表示正空间自相关，即相似的属性值在空间上趋于聚集；当I<0时，表示负空间自相关，即相似的属性值在空间上趋于分散；当I=0时，表示属性值在空间上呈随机分布。计算得到的土壤各项理化性质的全局Moran'sI指数及相关统计结果如表2所示。土壤pH值的Moran'sI指数为0.125，P值为0.086，在0.1的显著性水平下，表现出较弱的正空间自相关。这表明土壤pH值在一定程度上存在空间聚集现象，但聚集程度相对较弱，可能是由于研究区域内成土母质、地形等因素在较大范围内的相对均一性，使得土壤pH值的空间变化较为平缓。土壤电导率的Moran'sI指数为0.356，P值小于0.01，呈现显著的正空间自相关。说明土壤电导率在空间上具有较强的聚集性，高电导率区域和低电导率区域分别相对集中分布。如前所述，这主要与该区域的地形地貌和水分条件密切相关，在沙漠边缘的低洼区域，地下水位较高，盐分随水分蒸发向上运移并积累，导致这些区域土壤电导率较高且在空间上聚集分布；而在沙丘顶部等地势较高、排水良好的区域，土壤盐分含量低，电导率也低，同样呈现出空间聚集的特征。土壤有机质的Moran'sI指数为0.456，P值小于0.01，表现出显著的正空间自相关。土壤有机质主要来源于植物残体和根系分泌物等，在植被覆盖较好的丘间低地，植物生长繁茂，残体归还量大，土壤有机质含量高且在空间上聚集；而在植被稀疏的沙丘区域，土壤有机质来源少，含量低，也呈现出空间聚集的分布格局。土壤全氮的Moran'sI指数为0.389，P值小于0.01，具有显著的正空间自相关。全氮含量与土壤有机质含量密切相关，同时受植被类型和土地利用方式影响。绿洲农业用地施肥等农事活动使土壤全氮含量高且聚集分布；荒漠生态用地植被生长条件差，土壤全氮含量低且相对聚集。土壤速效氮的Moran'sI指数为0.423，P值小于0.01，呈现显著正空间自相关。速效氮受土壤微生物活动、氮肥施用和土壤水分等因素影响，在水分条件好、微生物活性高的区域，速效氮含量高且聚集；反之则低且聚集。土壤全磷的Moran'sI指数为0.256，P值为0.023，在0.05的显著性水平下表现出正空间自相关。土壤全磷含量主要取决于成土母质，在山前冲积扇等母质富含磷素的区域，土壤全磷含量高且相对聚集；在长期遭受风沙侵蚀的沙丘区域，土壤全磷含量低且聚集。土壤速效磷的Moran'sI指数为0.523，P值小于0.01，具有显著的正空间自相关。速效磷含量受土壤酸碱度、铁铝氧化物含量等因素影响，在碱性土壤和富含铁铝氧化物的区域，速效磷含量低且聚集；在其他区域则相对较高且聚集。土壤砂粒含量的Moran'sI指数为0.189，P值为0.056，在0.1的显著性水平下表现出较弱的正空间自相关。研究区地处沙漠边缘，风沙活动频繁，土壤以风沙土为主，砂粒含量较高且在一定程度上存在空间聚集现象，但由于风沙活动的随机性和局部地形的微变化，这种聚集程度相对较弱。土壤粉粒含量的Moran'sI指数为0.321，P值小于0.01，呈现显著的正空间自相关。粉粒含量的空间分布受到土壤搬运和沉积过程的影响，在不同地貌部位和植被覆盖条件下，粉粒的沉积和流失情况不同，导致高粉粒含量区域和低粉粒含量区域在空间上分别相对集中。土壤黏粒含量的Moran'sI指数为0.367，P值小于0.01，表现出显著的正空间自相关。黏粒含量受成土过程中物理、化学和生物作用的综合影响，在植被根系活动频繁、土壤团聚体稳定性较好的区域，黏粒含量较高且在空间上聚集；在其他区域则相对较低且聚集。为了更直观地展示土壤理化性质的空间自相关关系，绘制了Moran'sI散点图（图1）。以土壤有机质为例，Moran'sI散点图的横坐标为土壤有机质的标准化属性值，纵坐标为其空间滞后值。散点主要分布在第一象限（高-高聚集）和第三象限（低-低聚集），进一步验证了土壤有机质在空间上呈现显著的正空间自相关，即高值区域周围多为高值区域，低值区域周围多为低值区域。其他土壤理化性质的Moran'sI散点图也呈现出类似的分布特征，与全局Moran'sI指数的分析结果一致。通过空间自相关分析，明确了古尔班通古特沙漠南缘土壤各项理化性质在空间上的分布并非随机，而是存在一定的相关性和聚集特征，这为深入理解该区域土壤的空间异质性提供了重要依据。3.3土壤空间异质性特征运用地统计学中的半方差函数对古尔班通古特沙漠南缘土壤理化性质的空间变异结构进行深入分析。半方差函数能够有效描述土壤属性在空间上的变异程度和自相关特征。通过计算不同步长下土壤各理化性质的半方差值，并拟合相应的理论模型，得到了土壤各项理化性质的半方差函数参数，结果如表3所示。土壤pH值的半方差函数最佳拟合模型为球状模型，块金值为0.025，基台值为0.068，变程为3500m。块金值反映了由随机因素引起的小尺度空间变异，其值相对较小，说明土壤pH值的随机变异程度较低。基台值代表了区域化变量的总空间变异程度，变程表示在该距离范围内土壤属性具有空间自相关性。土壤pH值的变程较大，表明在3500m的范围内，土壤pH值具有一定的空间相关性，在该距离之外，空间相关性减弱。这可能是由于研究区域内成土母质、地形等因素在较大范围内相对均一，使得土壤pH值的变化较为平缓。土壤电导率的半方差函数拟合为指数模型，块金值为0.086，基台值为0.325，变程为1800m。块金值相对较大，说明土壤电导率受到随机因素的影响较为明显。电导率主要反映土壤中可溶性盐分的含量，其空间变异受到地形、水分条件等多种因素的综合影响。在沙漠边缘的低洼区域，地下水位较高，盐分随水分蒸发向上运移并积累，导致土壤电导率较高；而在沙丘顶部等地势较高、排水良好的区域，土壤盐分含量低，电导率也低。这种地形和水分条件的差异使得土壤电导率的空间分布较为复杂，随机变异程度较大。变程为1800m，表明在该距离内土壤电导率具有空间自相关性，超过此距离，空间相关性逐渐减弱。土壤有机质的半方差函数符合球状模型，块金值为0.865，基台值为2.568，变程为2200m。土壤有机质主要来源于植物残体和根系分泌物等，其空间分布与植被覆盖状况密切相关。在植被覆盖较好的丘间低地，植物生长繁茂，残体归还量大，土壤有机质含量高；而在植被稀疏的沙丘区域，土壤有机质来源少，含量低。这种植被分布的差异导致土壤有机质的空间变异较大，块金值和基台值都相对较高。变程为2200m，说明在该距离范围内土壤有机质的空间相关性较强，超过此距离，空间相关性逐渐减弱。土壤全氮的半方差函数拟合为球状模型，块金值为0.065，基台值为0.189，变程为2000m。全氮含量与土壤有机质含量密切相关，同时受植被类型和土地利用方式的影响。在绿洲农业用地中，由于施肥等农事活动，土壤全氮含量相对较高；而在荒漠生态用地中，自然植被生长条件相对较差，土壤全氮含量较低。这种土地利用方式和植被类型的差异使得土壤全氮含量存在一定的空间变异，块金值和基台值适中。变程为2000m，表明在该距离内土壤全氮含量具有空间自相关性。土壤速效氮的半方差函数符合指数模型，块金值为2.568，基台值为6.568，变程为1500m。速效氮是土壤中可被植物直接吸收利用的氮素形态，其含量受到土壤微生物活动、氮肥施用以及土壤水分等因素的影响。在水分条件较好、微生物活性较高的区域，土壤中有机氮的矿化作用较强，速效氮含量相对较高。这些因素的空间变异性导致土壤速效氮含量的空间变异较大，块金值和基台值较高。变程为1500m，说明在该距离范围内土壤速效氮具有空间自相关性，超过此距离，空间相关性逐渐减弱。土壤全磷的半方差函数拟合为球状模型，块金值为0.023，基台值为0.078，变程为2500m。土壤全磷含量主要取决于成土母质，同时也受到土壤侵蚀、淋溶等因素的影响。在山前冲积扇等母质富含磷素的区域，土壤全磷含量相对较高；而在长期遭受风沙侵蚀的沙丘区域，土壤中细颗粒物质被吹蚀，全磷含量相对较低。由于成土母质在较大范围内相对稳定，所以土壤全磷含量的随机变异较小，块金值较低。变程为2500m，表明在该距离内土壤全磷含量具有空间自相关性。土壤速效磷的半方差函数符合指数模型，块金值为1.568，基台值为4.568，变程为1200m。速效磷的含量除了与土壤全磷含量有关外，还受到土壤酸碱度、铁铝氧化物含量等因素的影响。在碱性土壤中，磷素容易与钙结合形成难溶性磷酸盐，降低速效磷含量；而在一些富含铁铝氧化物的区域，磷素会被吸附固定，也会导致速效磷含量较低。这些因素的复杂性使得土壤速效磷含量的空间变异较大，块金值和基台值较高。变程为1200m，说明在该距离范围内土壤速效磷具有空间自相关性，超过此距离，空间相关性逐渐减弱。土壤砂粒含量的半方差函数拟合为球状模型，块金值为0.015，基台值为0.045，变程为3000m。研究区地处沙漠边缘，风沙活动频繁，土壤以风沙土为主，砂粒含量较高且相对稳定。因此，土壤砂粒含量的随机变异较小，块金值较低。变程为3000m，表明在该距离内土壤砂粒含量具有空间自相关性。土壤粉粒含量的半方差函数符合指数模型，块金值为0.045，基台值为0.125，变程为1600m。粉粒含量的空间分布受到土壤搬运和沉积过程的影响，在不同地貌部位和植被覆盖条件下，粉粒的沉积和流失情况不同，导致其含量存在一定差异。这种差异使得土壤粉粒含量具有一定的空间变异，块金值和基台值适中。变程为1600m，说明在该距离范围内土壤粉粒含量具有空间自相关性，超过此距离，空间相关性逐渐减弱。土壤黏粒含量的半方差函数拟合为球状模型，块金值为0.035，基台值为0.105，变程为1800m。黏粒含量受成土过程中物理、化学和生物作用的综合影响，在植被根系活动频繁、土壤团聚体稳定性较好的区域，黏粒含量相对较高。这些因素的空间变化导致土壤黏粒含量存在一定的空间变异，块金值和基台值适中。变程为1800m，表明在该距离内土壤黏粒含量具有空间自相关性。为了更直观地展示土壤理化性质的空间分布特征，利用克里格插值法对土壤各项理化性质进行空间插值，并绘制了空间分布图（图2）。从土壤pH值空间分布图可以看出，研究区域内土壤pH值整体呈碱性，且空间分布相对较为均匀，在山前冲积扇和洪积平原等区域，土壤pH值略有升高，可能与这些区域的母质特性有关。土壤电导率空间分布图显示，电导率高值区域主要集中在沙漠边缘的低洼地带和部分绿洲内部的灌溉区，这些区域地下水位较高，盐分容易积累；而沙丘顶部等地势较高的区域，电导率较低。土壤有机质空间分布图表明，有机质含量高值区主要分布在丘间低地和植被覆盖较好的绿洲边缘，低值区主要位于沙丘区域，与植被分布格局密切相关。土壤全氮和速效氮的空间分布趋势与有机质相似，在绿洲农业用地和植被生长良好的区域含量较高。土壤全磷空间分布图显示，全磷含量高值区主要分布在山前冲积扇等母质富含磷素的区域，低值区分布在沙丘区域。土壤速效磷空间分布图表明，速效磷含量在不同区域差异较大，在碱性土壤和富含铁铝氧化物的区域含量较低。土壤砂粒含量空间分布图显示，研究区域内砂粒含量普遍较高，且空间分布相对稳定，在沙丘区域砂粒含量略有增加。土壤粉粒和黏粒含量空间分布图表明，它们在不同地貌部位和植被覆盖条件下存在一定差异，在丘间低地和绿洲边缘等区域，粉粒和黏粒含量相对较高。通过半方差函数分析和克里格插值法绘制的空间分布图，全面、直观地呈现了古尔班通古特沙漠南缘土壤理化性质的空间异质性特征。四、古尔班通古特沙漠南缘土壤尺度特征分析4.1不同尺度下土壤空间异质性变化规律为深入探究古尔班通古特沙漠南缘土壤空间异质性在不同尺度下的变化规律，本研究从根际、个体、种群、地貌和区域五个尺度对土壤的pH、电导率、有机质、全氮、速效氮、全磷和速效磷等理化性质进行分析。在根际尺度（10^{-3}m~10^{-2}m），土壤理化性质的异质性相对较小。以土壤pH为例，根际微生物的活动和根系分泌物虽然会对土壤微环境产生一定影响，但这种影响范围有限，导致土壤pH在根际尺度的变异程度较低。土壤电导率在根际尺度同样变化较小，这是因为根际附近盐分的迁移和积累主要受植物根系主动吸收和分泌过程的调控，且根际范围小，盐分来源相对单一，使得电导率的空间差异不明显。对于土壤养分，如有机质、全氮、速效氮、全磷和速效磷，根际微生物的分解作用和根系对养分的吸收利用虽有一定作用，但总体上根际尺度的养分异质性相对其他尺度并不突出，因为根际范围有限，土壤养分的来源和转化过程相对简单。在个体尺度（10^{-1}m~10^{0}m），土壤空间异质性明显增大。植物个体对土壤性质产生显著影响，以土壤pH值为例，不同植物种类的根系分泌物性质和数量不同，一些植物根系分泌的有机酸等物质会改变根际土壤的酸碱度。如梭梭根系分泌物可能会使根际土壤pH值略有降低，导致在个体尺度上土壤pH值围绕植物个体出现明显的空间变异。土壤电导率在个体尺度的异质性也较为突出，植物个体通过蒸腾作用影响土壤水分运动，进而影响盐分的分布。在干旱环境下，植物根系吸水导致根际周围土壤水分减少，盐分相对浓缩，使得根际附近土壤电导率升高，而远离植物个体的区域电导率相对较低，形成明显的空间差异。土壤有机质、全氮、速效氮和速效磷在个体尺度的异质性达到最大，这是因为植物个体的凋落物、根系分泌物以及根系周转等过程直接向土壤中输入有机物质和养分，在植物个体周围形成“肥岛”效应。例如，红柳个体周围，其凋落物丰富，土壤有机质和养分含量明显高于周边区域，导致土壤养分在个体尺度呈现出强烈的空间变异。种群尺度（10^{0}m~10^{1}m）下，土壤空间异质性相对个体尺度有所减小。以土壤pH为例，虽然种群内植物个体对土壤pH有影响，但种群整体的一致性使得土壤pH的变化相对个体尺度更为缓和。同一植物种群的根系分泌物和对土壤的影响方式具有一定相似性，使得土壤pH在种群尺度的空间变异程度降低。土壤电导率在种群尺度的异质性也减小，这是因为种群内植物个体的分布相对均匀，对土壤水分和盐分的影响在一定程度上相互抵消，使得土壤盐分分布在种群尺度上相对均匀。对于土壤有机质、全氮、速效氮和速效磷，虽然种群内植物个体共同向土壤输入有机物质和养分，但由于种群分布范围较大，养分的扩散和均化作用使得土壤养分在种群尺度的空间异质性小于个体尺度。然而，种群内不同个体之间仍存在一定差异，导致土壤养分在种群尺度仍具有一定程度的空间变异。在地貌尺度（10^{1}m~10^{2}m），土壤空间异质性又有所增大。不同地貌部位的土壤理化性质存在明显差异，土壤pH在沙丘顶部和丘间低地表现出不同的值，这是由于地貌部位影响土壤的水分状况和盐分运移。沙丘顶部地势高，水分容易流失，盐分不易积累，土壤pH值相对较低；而丘间低地地势低洼，水分相对较多，盐分容易积累，土壤pH值相对较高，从而在地貌尺度上形成明显的空间变异。土壤电导率在地貌尺度也呈现出显著差异，如在沙漠边缘的低洼地带，地下水位高，盐分随水分蒸发向上运移并积累，土壤电导率高；而在沙丘顶部等地势较高、排水良好的区域，土壤电导率低。对于土壤有机质、全氮、速效氮和速效磷，地貌部位影响植被的生长和分布，进而影响土壤养分的积累和分布。丘间低地植被生长较好，土壤有机质和养分含量相对较高；而沙丘顶部植被稀疏，土壤有机质和养分含量较低，导致土壤养分在地貌尺度的空间异质性较大。土壤全磷在地貌尺度的异质性最大，这是因为地貌特征对土壤母质的分布和侵蚀过程影响显著。在山前冲积扇等母质富含磷素的区域，土壤全磷含量高；而在长期遭受风沙侵蚀的沙丘区域，土壤中细颗粒物质被吹蚀，全磷含量低，使得土壤全磷在地貌尺度的空间变异最为明显。区域尺度（10^{3}m~10^{4}m）下，土壤空间异质性相对复杂。土壤pH在区域尺度上的变化受到气候、母质等大尺度因素的综合影响。研究区气候干旱，降水稀少，土壤中的碱性物质难以被淋溶，使得区域内土壤整体呈碱性，但在不同区域，由于母质类型的差异，土壤pH仍存在一定的空间变异。土壤电导率在区域尺度的异质性也受到气候和地形的影响，在区域内不同的地形地貌和水分条件下，土壤盐分含量不同，导致电导率存在空间差异。对于土壤有机质、全氮、速效氮和速效磷，区域尺度上的空间异质性受到植被类型、土地利用方式等因素的影响。绿洲农业用地中，由于施肥等农事活动，土壤养分含量相对较高；而在荒漠生态用地中，自然植被生长条件相对较差，土壤养分含量较低，从而在区域尺度上形成土壤养分的空间变异。综上所述，土壤pH在各尺度异质性的大小顺序为：个体>种群>地貌>区域>根际；土壤电导率在各尺度异质性的大小顺序为：个体>区域>种群>地貌>根际；土壤有机质、全氮、速效氮和速效磷的空间异质性均在个体尺度上最大，地貌尺度上次之，而种群尺度上最小；土壤全磷的异质性则在地貌尺度最大，个体尺度次之，根际尺度最小。这些变化规律表明，不同尺度下影响土壤空间异质性的主导因素不同，在小尺度上，植物个体的生物作用对土壤性质影响显著；随着尺度增大，地貌、气候等物理因素以及植被类型、土地利用方式等综合因素的作用逐渐增强。4.2土壤空间异质性与尺度的关系在根际尺度，土壤微生物活动和根系分泌物是影响土壤空间异质性的关键因素。土壤微生物的代谢活动会改变根际土壤的酸碱度、养分含量和氧化还原电位等性质。一些微生物能够分解有机物质，释放出氮、磷、钾等养分，增加土壤中速效养分的含量；同时，微生物在生长过程中会分泌有机酸、多糖等物质，这些物质可以与土壤中的金属离子结合，影响土壤的理化性质。根系分泌物中含有糖类、氨基酸、有机酸等成分，它们不仅为根际微生物提供了碳源和能源，还能直接参与土壤中物质的转化和迁移过程。例如，根系分泌的有机酸可以溶解土壤中的难溶性磷，提高磷的有效性，从而导致根际土壤中速效磷含量与非根际土壤存在差异。随着尺度增大到个体尺度，植物个体对土壤空间异质性的影响逐渐凸显。植物的生长、发育和代谢过程会对周围土壤环境产生显著影响。不同植物种类具有不同的根系形态、分布和生理特性，从而导致其对土壤养分的吸收和利用方式存在差异。深根系植物能够吸收深层土壤中的水分和养分，而浅根系植物主要利用表层土壤资源。这种根系分布的差异使得不同植物个体周围土壤的养分含量和分布呈现出明显的空间异质性。植物的凋落物也是影响土壤空间异质性的重要因素。凋落物分解后会向土壤中输入有机物质和养分，在植物个体周围形成“肥岛”效应。研究表明，在古尔班通古特沙漠南缘，梭梭等植物的凋落物富含氮、磷等养分，分解后显著提高了周围土壤的肥力，使得土壤有机质、全氮、速效氮等含量在个体尺度上呈现出明显的空间变异。在种群尺度，植物种群的结构和动态对土壤空间异质性产生重要影响。种群内植物个体之间的竞争和相互作用会改变土壤资源的分配格局。当种群密度较高时，植物个体之间对水分、养分等资源的竞争加剧，导致土壤资源在种群内的分布更加均匀；而当种群密度较低时，个体之间的竞争相对较弱，土壤资源的分布则相对不均匀。植物种群的更新和演替过程也会影响土壤空间异质性。在种群演替的不同阶段，植物种类和群落结构发生变化，从而导致土壤性质的改变。在古尔班通古特沙漠南缘，随着植被的恢复和演替，土壤有机质含量逐渐增加，土壤结构得到改善，土壤空间异质性也相应发生变化。地貌尺度上，地形地貌特征成为影响土壤空间异质性的主导因素。不同地貌部位的土壤形成过程和物质迁移规律存在显著差异。在沙丘顶部，由于风力侵蚀作用强烈，土壤颗粒较粗，保水保肥能力差，土壤养分含量较低；而丘间低地则容易积水，土壤颗粒较细，保水保肥能力较强，土壤养分含量相对较高。地形起伏还会影响土壤水分和盐分的运移。在地势较高的区域，土壤水分容易流失，盐分不易积累；而在低洼区域，水分汇聚，盐分随水分蒸发向上运移并积累，导致土壤盐分含量升高，电导率增大。在古尔班通古特沙漠南缘，沙漠边缘的低洼地带土壤电导率明显高于沙丘顶部，这与地貌对土壤水分和盐分运移的影响密切相关。区域尺度上，气候、母质、植被类型和土地利用方式等多种因素综合作用，导致土壤空间异质性更加复杂。气候因素中的降水和温度对土壤的风化、淋溶和有机质分解等过程产生重要影响。在降水较多的区域，土壤淋溶作用强烈，盐分含量较低，土壤pH值相对较低；而在干旱区域，土壤淋溶作用弱，盐分积累，土壤pH值较高。母质是土壤形成的物质基础，不同母质类型的土壤在矿物组成、质地和养分含量等方面存在差异。研究区山前冲积扇的母质富含矿物质，土壤颗粒较粗，而洪积平原的母质相对较细，这导致不同区域土壤的质地和养分含量存在明显的空间变异。植被类型和土地利用方式也对土壤空间异质性产生重要影响。绿洲农业用地中，由于施肥、灌溉等农事活动，土壤养分含量较高，且空间分布相对均匀；而荒漠生态用地中，自然植被生长条件相对较差，土壤养分含量较低，空间异质性较大。综上所述，随着尺度的增大，影响土壤空间异质性的主导因素从微观的生物因素逐渐转变为宏观的物理和地理因素，以及多种因素的综合作用。不同尺度下土壤空间异质性的变化规律和主导影响因子的差异，对于深入理解土壤形成和演化过程，以及合理进行土地利用和生态环境保护具有重要意义。4.3尺度效应分析运用分形理论对古尔班通古特沙漠南缘土壤的空间异质性进行深入研究，通过计算土壤属性的分形维数来揭示其尺度效应。分形维数能够定量描述土壤空间分布的复杂性和不规则性，分形维数越大，表明土壤空间异质性越强。采用容量维数法计算土壤各项理化性质的分形维数。以土壤有机质为例，通过改变采样间距，分别计算不同采样间距下土壤有机质的变异系数，以\ln(1/r)为横坐标，\lnN(r)为纵坐标进行线性回归分析，得到回归直线的斜率即为分形维数。计算结果显示，土壤有机质的分形维数为2.76，这表明土壤有机质在空间分布上具有较强的异质性。与其他土壤理化性质相比，土壤有机质的分形维数相对较高，进一步说明其空间分布的复杂性和不规则性更为突出。通过分析不同尺度下土壤理化性质的分形维数变化规律，深入探究尺度效应。结果表明，随着尺度的增大，土壤理化性质的分形维数呈现出不同的变化趋势。对于土壤pH值，分形维数在小尺度下相对较小，随着尺度的增大，分形维数逐渐增大，当尺度达到一定程度后，分形维数趋于稳定。这说明在小尺度上，土壤pH值的空间变异相对较小，主要受局部因素如根际微生物活动和根系分泌物的影响；随着尺度的增大，地形、母质等大尺度因素的影响逐渐增强，导致土壤pH值的空间异质性增大。当尺度足够大时，这些大尺度因素的作用趋于稳定，使得分形维数不再明显变化。土壤电导率的分形维数在不同尺度下的变化较为复杂。在小尺度上，由于土壤盐分的局部聚集和迁移，分形维数较大；随着尺度的增大，盐分的分布逐渐受到地形、水分等因素的综合影响，分形维数先减小后增大。在中等尺度上，盐分的分布相对均匀，分形维数达到最小值；当尺度继续增大时，区域内不同地形地貌和水分条件导致盐分分布的差异增大，分形维数再次增大。土壤养分（如全氮、速效氮、全磷、速效磷）的分形维数在小尺度上主要受植物个体和种群的影响，分形维数较大。随着尺度的增大，地貌、气候等因素的作用逐渐增强，分形维数呈现出先减小后增大的趋势。在中等尺度上，植物种群的均化作用和地貌因素的综合影响使得土壤养分的空间分布相对均匀，分形维数达到最小值；当尺度增大到区域尺度时，气候、母质、植被类型和土地利用方式等多种因素的综合作用导致土壤养分的空间异质性增大，分形维数再次增大。尺度效应在土壤研究中具有重要作用。在土壤采样设计方面，充分考虑尺度效应可以提高采样的代表性和准确性。例如，在研究土壤有机质时，由于其在小尺度上的空间异质性较大，为了准确反映其空间分布特征，在小尺度范围内应适当增加采样点的密度；而在大尺度上，采样点的间距可以适当增大。在土壤质量评价中，尺度效应也不容忽视。不同尺度下土壤质量的评价指标和标准可能不同，考虑尺度效应可以更加客观、准确地评价土壤质量。在制定土壤保护和管理策略时，根据尺度效应分析结果，可以针对不同尺度上的主导影响因素，采取相应的措施。在小尺度上，注重改善植物生长环境，促进植物对土壤养分的吸收和利用；在大尺度上，加强对地形、气候等因素的调控，保护土壤生态环境。五、影响土壤空间异质性及其尺度特征的因素5.1自然因素5.1.1地形地貌古尔班通古特沙漠南缘地形地貌复杂，涵盖了沙丘、丘间低地、山前冲积扇、洪积平原等多种类型，这些地形地貌的差异深刻影响着土壤空间异质性。地形起伏对土壤水分分布起着关键作用。在地势较高的区域，如沙丘顶部，重力作用使水分迅速下渗和流失，导致土壤含水量较低。相关研究表明，沙丘顶部的土壤水分含量可比丘间低地低30%-50%。而在丘间低地等地势低洼处，水分汇聚，土壤水分条件相对较好。这种水分分布的差异进一步影响土壤养分的迁移和积累。水分是养分运输的载体，在水分含量高的丘间低地，养分随水分的运移和淋溶作用更为活跃，使得土壤中可溶性养分含量相对较高；而在沙丘顶部，由于水分缺乏，养分的迁移和转化受到限制，土壤养分含量较低。坡度对土壤侵蚀和物质迁移也有重要影响。在坡度较大的区域，如沙丘的迎风坡，风力侵蚀作用强烈，土壤颗粒被吹蚀搬运，导致土壤质地变粗，细颗粒物质减少。研究发现，迎风坡土壤中的黏粒和粉粒含量比背风坡低20%-30%。同时，坡度影响降水的再分配，坡面径流会携带土壤中的养分和细颗粒物质向下坡方向移动，造成上坡土壤养分流失，下坡土壤养分富集。在古尔班通古特沙漠南缘，这种因坡度导致的土壤养分空间差异在植被生长状况上也有明显体现，下坡植被生长相对茂盛，而上坡植被则较为稀疏。坡向主要通过影响光照和热量条件来改变土壤水分蒸发和植被生长，进而影响土壤空间异质性。阳坡接受的太阳辐射多，温度较高，土壤水分蒸发量大，导致土壤相对干燥。阴坡则相反，光照较弱，温度较低，土壤水分蒸发量小，土壤含水量相对较高。不同的水分条件又影响植被的分布和生长，阳坡植被多为耐旱性较强的物种，而阴坡则更适合对水分需求较高的植物生长。植被的差异进一步影响土壤的理化性质，如植被凋落物的数量和质量不同，会导致土壤有机质含量和养分组成在阳坡和阴坡存在差异。在一些研究中发现，阴坡土壤的有机质含量可比阳坡高10%-20%。5.1.2气候条件古尔班通古特沙漠南缘属于温带大陆性干旱气候，降水、气温、风力等气候因素对土壤空间异质性产生重要影响。降水是影响土壤水分和物质迁移的关键因素。该地区年降水量稀少且分布不均，年降水量通常在100-200毫米之间，且多集中在春季和夏季。降水的空间差异导致土壤水分含量在不同区域存在明显变化。在降水较多的区域，土壤水分含量较高，有利于土壤中矿物质的溶解和淋溶，使土壤中的盐分和养分被带走，从而影响土壤的化学性质。而在降水较少的区域，土壤水分不足，植被生长受到限制，土壤有机质积累较少，土壤肥力较低。降水的季节性变化也对土壤空间异质性产生影响。在降水集中的季节，大量降水可能导致土壤侵蚀加剧，使土壤中的细颗粒物质和养分流失，改变土壤的质地和养分分布。气温对土壤水分蒸发和土壤中生物化学反应速率有重要影响。在高温季节，土壤水分蒸发强烈，导致土壤水分含量降低。同时，气温升高会加速土壤中有机质的分解和微生物的活动，影响土壤养分的转化和循环。在古尔班通古特沙漠南缘，夏季气温较高，土壤有机质的分解速率明显加快，导致土壤中速效养分含量在夏季相对较低。而在冬季，气温较低，土壤微生物活动受到抑制，土壤中生物化学反应速率减慢，土壤养分的转化和迁移也相应减缓。风力是该地区气候的重要特征之一，对土壤空间异质性的影响主要体现在土壤侵蚀和风沙堆积方面。强劲的风力会将沙漠表面的细颗粒物质吹起并搬运到其他地方，导致土壤质地发生变化。在风力侵蚀严重的区域，土壤中的细颗粒物质被吹走，留下粗颗粒的砂粒，土壤质地变粗，保水保肥能力下降。而在风力较弱的区域或有植被阻挡的地方，风沙物质会堆积下来，使土壤厚度增加，土壤质地也会发生改变。风力还会影响土壤中盐分的分布，在风沙活动过程中，盐分可能会随着风沙物质的搬运而重新分布，导致土壤盐分含量在空间上出现差异。5.1.3植被类型古尔班通古特沙漠南缘植被类型多样，主要包括梭梭、红柳、琵琶柴、芨芨草等荒漠植被以及绿洲边缘的人工植被。不同植被类型通过根系分布、凋落物归还等对土壤理化性质和空间异质性产生显著影响。根系是植物与土壤相互作用的重要界面，不同植被的根系分布特征各异。深根系植物如梭梭，其根系可深入地下数米，能够吸收深层土壤中的水分和养分，从而改变土壤水分和养分在垂直方向上的分布。研究表明，梭梭根系周围土壤在深层（20-50cm）的水分含量和养分含量相对较高，形成独特的土壤水分和养分分布格局。浅根系植物如芨芨草，根系主要分布在表层土壤（0-20cm），对表层土壤的物理结构和化学性质影响较大。芨芨草根系的生长和活动可以增加土壤孔隙度，改善土壤通气性和透水性，同时，根系分泌物还能影响土壤微生物的群落结构和活性，进而影响土壤养分的转化和循环。植被凋落物是土壤有机质的重要来源，不同植被的凋落物数量和质量存在差异。梭梭和红柳等灌木的凋落物富含木质素和纤维素，分解速度较慢，但能为土壤提供长期稳定的有机质输入。研究发现，在梭梭和红柳群落中，土壤有机质含量随着凋落物积累时间的增加而逐渐升高。而草本植物如芨芨草的凋落物分解速度相对较快，在短期内能够为土壤提供较多的速效养分。凋落物的分解过程还会释放出二氧化碳、氮氧化物等气体，影响土壤的气体组成和酸碱度。此外，凋落物在土壤表面的覆盖还能减少土壤水分蒸发，降低土壤侵蚀风险，对土壤空间异质性产生间接影响。五、影响土壤空间异质性及其尺度特征的因素5.2人为因素5.2.1土地利用方式在古尔班通古特沙漠南缘，不同土地利用方式对土壤结构、养分含量和空间异质性产生显著影响。耕地作为主要的农业用地类型，长期的耕作活动深刻改变了土壤结构。频繁的翻耕、耙地等农事操作，使得土壤团聚体遭到破坏，大团聚体减少，小团聚体增加。相关研究表明，与未开垦的自然土壤相比，耕地土壤的团聚体稳定性降低了20%-30%。这种结构变化影响土壤的通气性和透水性，进而改变土壤水分和养分的运移和储存。耕地的养分含量与施肥、灌溉等农事管理密切相关。为满足农作物生长需求，农民通常会大量施用化肥和有机肥，导致土壤中氮、磷、钾等养分含量明显高于自然土壤。然而，不合理的施肥方式，如过量施肥或施肥时间不当，会造成养分的流失和浪费，使土壤养分在空间分布上更加不均。长期不合理灌溉会导致土壤盐分积累，尤其是在干旱地区，灌溉水的蒸发会使盐分在土壤表层富集，增加土壤电导率，影响土壤的化学性质和空间异质性。林地对土壤具有重要的保护和改良作用。树木的根系能够深入土壤，增强土壤的稳定性，减少土壤侵蚀。研究表明，林地土壤的抗侵蚀能力比裸地提高了50%以上。树木的凋落物为土壤提供了丰富的有机质来源，分解后增加土壤有机质含量，改善土壤结构和肥力。不同树种的林地对土壤性质的影响存在差异，阔叶树的凋落物富含氮、磷等养分，分解后能显著提高土壤肥力；针叶树的凋落物则含有较多的木质素，分解速度较慢，但能长期维持土壤的有机质含量。草地在维持土壤生态系统平衡方面发挥着重要作用。草本植物的根系密集且浅，能够有效固定土壤颗粒，防止土壤侵蚀。草地植被的覆盖还能减少土壤水分蒸发，保持土壤水分。草地土壤的有机质含量相对较高，这得益于草本植物的凋落物和根系分泌物。然而，过度放牧会破坏草地植被，导致土壤裸露，土壤结构被破坏，土壤肥力下降。研究发现，过度放牧区域的草地土壤有机质含量比适度放牧区域降低了15%-25%。建设用地的扩张导致大量自然土壤被覆盖和压实，土壤的原有结构和功能遭到严重破坏。建筑物、道路等基础设施的建设改变了土壤的物理性质，使土壤通气性和透水性变差，土壤微生物活动受到抑制。建设用地周边的土壤往往受到人为活动的污染，如重金属污染、有机物污染等，导致土壤质量下降，空间异质性增加。5.2.2农业活动灌溉是农业生产中的关键环节，对古尔班通古特沙漠南缘土壤水分、养分分布和空间异质性产生重要影响。在该地区，由于降水稀少，灌溉成为维持农作物生长的主要水源补给方式。不同的灌溉方式会导致土壤水分分布的差异，漫灌是一种较为粗放的灌溉方式，容易造成水分在田面的不均匀分布，使得土壤水分含量在水平方向上存在较大差异。研究表明，漫灌条件下，田块边缘和中心的土壤水分含量差值可达10%-15%。滴灌和喷灌等精准灌溉方式能够更均匀地将水分输送到作物根系周围，减少水分的浪费和流失，使土壤水分在空间上的分布相对均匀。然而，长期采用单一的灌溉方式，如过度依赖滴灌，可能导致土壤盐分在滴头附近积累，形成盐分高值区，增加土壤空间异质性。施肥是调节土壤养分含量和供应的重要农业措施，但不合理的施肥会对土壤空间异质性产生负面影响。在古尔班通古特沙漠南缘，农民在施肥过程中往往存在盲目性，过量施用氮肥和磷肥，而忽视钾肥和微量元素肥料的使用。这种不平衡的施肥方式导致土壤中氮、磷等养分含量过高，而其他养分相对缺乏，使得土壤养分在空间分布上呈现出斑块状差异。研究发现，过量施肥区域的土壤中氮、磷含量比正常施肥区域高出30%-50%，而钾和微量元素含量则相对较低。此外，施肥方法不当，如撒施不均匀，也会造成土壤养分在局部区域的富集或缺乏，进一步加剧土壤空间异质性。耕作是改变土壤物理结构和化学性质的重要农事活动。传统的深耕和旋耕方式能够打破土壤板结，增加土壤通气性和透水性，但过度耕作会破坏土壤团聚体结构，导致土壤肥力下降。在古尔班通古特沙漠南缘，长期的高强度耕作使得土壤表层的有机质含量降低，土壤颗粒变得松散，容易受到风力侵蚀。免耕和少耕等保护性耕作方式近年来得到推广，这些方式能够减少对土壤的扰动，保持土壤结构的稳定性，有利于土壤有机质的积累和微生物的活动。研究表明，采用免耕和少耕的农田，土壤有机质含量比传统耕作农田提高了10%-15%，土壤空间异质性相对较小。然而，保护性耕作方式在实施过程中也面临一些挑战，如杂草控制和病虫害防治等问题，需要进一步研究和解决。六、结论与展望6.1研究主要结论本研究对古尔班通古特沙漠南缘土壤的空间异质性及其尺度特征进行了系统研究，取得以下主要结论：古尔班通古特沙漠南缘土壤理化性质呈现出明显的空间异质性。土壤pH值整体呈碱性，且空间分布相对较为均匀，变异系数为2.89%，变异性较弱。土壤电导率、有机质、全氮、速效氮、全磷、速效磷等含量在空间上存在显著差异，变异系数分别为44.80%、52.10%、40.70%、52.00%、32.30%、71.60%，表现出中等至强变异性。土壤砂粒含量较高且相对稳定，变异系数为11.30%，变异性较弱；粉粒和黏粒含量存在一定空间差异，变异系数分别为35.50%和41.20%，呈现中等变异性。通过空间自相关分析发现，土壤各项理化性质在空间上存在显著的正空间自相关，高值区域和低值区域分别相对聚集。利用半方差函数分析表明，不同土壤理化性质具有不同的空间变异结构，变程范围在1200-3500m之间，块金值与基台值的比值反映了随机因素和结构性因素对土壤空间变异的相对贡献。通过克里格插值绘制的空间分布图直观展示了土壤理化性质的空间分布格局。不同尺度下土壤空间异质性变化规律明显。在根际尺度，土壤理化性质的异质性相对较小；在个体尺度，土壤空间异质性明显增大，尤其是土壤有机质、全氮、速效氮和速效磷的空间异质性达到最大；在种群尺度，土壤空间异质性相对个体尺度有所减小；在地貌尺度，土壤空间异质性又有所增大，其中土壤全磷的异质性在地貌尺度最大；在区域尺度，土壤空间异质性相对复杂，受多种因素综合影响。土壤pH在各尺度异质性的大小顺序为：个体>种群>地貌>区域>根际；土壤电导率在各尺度异质性的大小顺序为：个体>区域>种群>地貌>根际；土壤有机质、全氮、速效氮和速效磷的空间异质性均在个体尺度上最大，地貌尺度上次之，而种群尺度上最小；土壤全磷的异质性则在地貌尺度最大，个体尺度次之，根际尺度最小。土壤空间异质性与尺度密切相关，随着尺度的增大，影响土壤空间异质性的主导因素从微观的生物因素逐渐转变为宏观的物理和地理因素，以及多种因素的综合作用。在根际尺度，土壤微生物活动和根系分泌物是主要影响因素；在个体尺度，植物个体的生长、发育和代谢过程对土壤空间异质性影响显著；在种群尺度，植物种群的结构和动态起重要作用；在地貌尺度，地形