Please wait a minute...

当期目录

    2018年, 第73卷, 第9期 刊出日期:2018-09-25 上一期    下一期
    土地利用
    基于地貌分区的1990-2015年中国耕地时空特征变化分析
    程维明,高晓雨,马廷,徐新良,陈印军,周成虎
    2018, 73 (9):  1613-1629.  doi: 10.11821/dlxb201809001
    摘要 ( 1052 )   HTML ( 54 )   PDF(8195KB) ( 2301 )  

    地形地貌通过分配地表水分和热量制约耕地的利用形式和成效,对耕地质量具有非常重要的决定意义,以地貌分区的视角研究中国耕地变化具有重要意义。根据1990-2015年中国6期土地利用空间数据和地貌分区数据,运用GIS空间叠加分析方法,分析中国耕地1990-2015年的面积变化和空间分布格局,进一步探讨不同地貌分区下的耕地新增与流失方向。结果表明,中国耕地面积略有增加,但总体变化不大,耕地面积从1990年的17715万hm2增加到2015年的17851万hm2,平均每年增加5.44万hm2,每年增幅仅为0.03%。耕地主要分布在平原地区,台地、丘陵次之;新增耕地主要来源为草地、林地和未利用地。东部平原低山丘陵区(I)耕地面积最大,而西北高中山盆地高原区(IV)耕地动态度明显高于其他地貌区。空间上呈现“南减北增,新增耕地的重心向西北移动”的特征。1990-2015年间,西北高中山盆地高原区(IV)和东部平原低山丘陵区(I)为耕地面积增长区,平均每年增加耕地面积分别为8.9万hm2和5.4万hm2;东南低山丘陵平原区(II)和西南中低山高原盆地区(V)为耕地面积减少区,平均每年减少耕地面积分别为5.9万hm2和2.8万hm2;而华北—内蒙东中山高原区(III)和青藏高原高山极高山盆地谷地区(VI)耕地面积几乎没有变化,平均每年变化仅为0.15 万hm2和0.06万 hm2。耕地流失主要发生在东部平原低山丘陵区(I)和东南低山丘陵平原区(II),主要原因为城镇化进程加快带来的建设用地对优质耕地的大量占用;而新增耕地主要发生在西北高中山盆地高原区(IV),多来自于对草地和未利用地的开垦。

    参考文献 | 相关文章 | 计量指标
    三峡库区规模农地时空变化特征及其驱动机制
    梁鑫源,李阳兵
    2018, 73 (9):  1630-1646.  doi: 10.11821/dlxb201809002
    摘要 ( 640 )   HTML ( 34 )   PDF(6350KB) ( 1042 )  

    探讨库区农地利用的转型及其机制,可反映中西部结合带、山区、库区的农村变化和人与环境关系的变化,为典型山区的农业多功能转型发展方向提供科学依据。选取三峡库区腹地5个区县为研究区,通过定义“规模农地”相关概念,利用ArcGIS 10.2软件并结合SPSS、地理空间探测器等技术手段对规模农地的变化与机制进行分析。研究发现:规模农地的转移来源主要由耕地和灌木林地组成,除部分高山反季节蔬菜,规模农地的增长一般发生在低海拔地区,且在坡度30°范围内均有分布;不同类型的规模农地空间变化有所不同,但均与道路和乡镇级行政中心的关系明显。斑块水平上,高程、坡度等自然因素对规模农地的贡献相对较大;乡镇水平上以道路交通、农业人口密度等社会经济与人文因素驱动力为主;区域水平上,各区县的驱动因子影响力较研究区整体尺度更为显著。规模农地的变化与发展受多种因素驱动,但农业政策始终是各驱动因素的主导,不同水平的行为主体会促使规模农地发展趋向不同的方向,但原则上均受国家农业发展现状与政府政策导向影响。

    参考文献 | 相关文章 | 计量指标
    城市化时期郊区土地利用结构信息熵上升的原因
    王玉明,王瑞康
    2018, 73 (9):  1647-1657.  doi: 10.11821/dlxb201809003
    摘要 ( 875 )   HTML ( 37 )   PDF(1614KB) ( 1218 )  

    本文认为,此前一直被当作城市化地区的土地利用结构信息熵快速上升,应是指城市化地区郊区的土地利用结构信息熵上升。城市化时期郊区土地利用结构向均衡化方向的明显变动,引起了郊区土地利用结构信息熵的上升;城市建成区的土地利用结构信息熵却未有上升,其各种土地利用类型的面积占比在城市化时期基本稳定,或向居民点及工矿用地、交通用地类型有所集中,导致其土地利用结构信息熵基本不变或有所下降。郊区的局部城市化区域对郊区土地利用结构信息熵的上升有着主要贡献。数据分析显示,城市化时期郊区的农用地较多较快地转化成了城市建设用地,且后者的面积占比与土地利用结构信息熵之间具有高度正相关关系,说明城市建设用地的快速增加是引起郊区土地利用结构信息熵上升的主要原因。熵变化分析表明,城市化时期郊区系统的熵变化与郊区土地利用结构信息熵的变动方向一致,均大于零,表征了郊区的自然地理环境演化在此期间是向着有序性降低的方向发展的。

    参考文献 | 相关文章 | 计量指标
    撂荒耕地的提取与分析——以山东省庆云县和无棣县为例
    肖国峰,朱秀芳,侯陈瑶,夏兴生
    2018, 73 (9):  1658-1673.  doi: 10.11821/dlxb201809004
    摘要 ( 753 )   HTML ( 24 )   PDF(6472KB) ( 1092 )  

    由于城镇化的快速发展和农村劳动力不断流失,中国部分地区出现大面积的撂荒现象。利用遥感技术可以间接提取撂荒耕地的规模和数量,对耕地的保护和粮食安全有重要意义。以山东省庆云县和无棣县为研究区,基于Landsat数据和HJ1A数据,采用CART决策树分类方法,制作了1990-2017年的土地利用图,制定了撂荒地的识别规则,在此基础之上提取了研究区的撂荒地空间分布、持续撂荒时间分布和撂荒地复垦区域。结果显示:1990-2017年基准期影像的CART决策树分类精度高于85%;1992-2017年间,研究区撂荒地面积最大值为5503.86 hm2,最大撂荒率为5.37%,其中1996-1998年撂荒率最高,2006-2017年撂荒地面积的整体趋势逐年降低;1992-2017年间最大持续撂荒时间为15年,大部分耕地持续撂荒时间在4年之内,少数耕地持续撂荒时间超过10年;1993-2017年撂荒耕地复垦面积最大为2022.3 hm2,最小复垦面积约为20 hm2,复垦率最大值为67.44%,年均复垦率为31.83%。研究结果不仅能够为研究区撂荒驱动因素分析提供数据支撑,而且也可以为其他地区的撂荒耕地识别提供参考。

    参考文献 | 相关文章 | 计量指标
    地表过程
    基于地理探测器的喀斯特不同地貌形态类型区土壤侵蚀定量归因
    王欢,高江波,侯文娟
    2018, 73 (9):  1674-1686.  doi: 10.11821/dlxb201809005
    摘要 ( 708 )   HTML ( 30 )   PDF(2849KB) ( 953 )  

    土壤侵蚀形成机制与影响因素识别是当前研究的核心与前沿议题,然而从多因素综合作用的角度进行定量归因仍需加强。以喀斯特典型峰丛洼地流域为研究区,基于GIS手段和RUSLE模型模拟土壤侵蚀,综合土地利用、坡度、降雨、高程、岩性、植被覆盖度等影响因子,应用地理探测器方法针对喀斯特不同地貌形态类型区进行土壤侵蚀的定量归因研究。结果表明,各影响因子对土壤侵蚀的解释力及因子间耦合作用程度在不同地貌形态类型区差异显著,其中土地利用和坡度是决定土壤侵蚀空间异质的主导因子,但在山地丘陵区,随着地形起伏度的升高,坡度的控制作用下降,即地理探测器q值表现为中海拔丘陵>小起伏中山>中起伏中山;生态探测器显示土地利用对土壤侵蚀的影响相比于其他因子有显著差异;双因子交互作用有助于增强对土壤侵蚀的解释力,土地利用与坡度的协同作用对土壤侵蚀的解释力达到70%以上;对于土壤侵蚀空间分布的差异性检验,风险探测器显示在小起伏中山、中起伏中山等地貌形态类型中,具有显著差异的影响因子分层组合数占比至少55%。因而,喀斯特地区土壤侵蚀的治理应综合考虑不同地貌形态类型区土壤侵蚀影响机制的空间异质性。

    参考文献 | 相关文章 | 计量指标
    短时间高强度旅游活动下洞穴CO2的变化特征及对滴水水文地球化学的响应
    张结,周忠发,汪炎林,潘艳喜,薛冰清,张昊天,田衷珲
    2018, 73 (9):  1687-1701.  doi: 10.11821/dlxb201809006
    摘要 ( 409 )   HTML ( 8 )   PDF(4652KB) ( 587 )  

    洞穴高强度旅游活动产生的CO2对洞穴滴水水文地球化学及洞穴沉积物沉积具有重要影响。本文于2017年9月30日-10月9日对贵州绥阳大风洞洞穴CO2、温度、相对湿度、游客数量及洞穴滴水水化学指标等进行连续监测,运用系统分析方法对各要素进行综合分析。结果发现:① 受游客数量和洞穴通风效应等因素影响,洞穴空气CO2分压(PCO2(A))在时间序列上呈现明显的昼夜变化和日际变化,表现为白昼高、夜间低,游客多的天数高,游客少的天数低。在空间变化上,由于通风程度和洞腔体积不同,不同监测点PCO2(A)存在明显差异,由洞内深处至洞口分别为3#(神泉玉露)>1#(时光隧道)>2#(夜明珠);② 通过比较PCO2(A)和滴水水温,前者对CO2溶解度影响比后者更为显著,表明PCO2(A)是洞穴沉积过程中最重要的驱动因素之一;③ 洞穴滴水水温、滴水PCO2分压(PCO2(W))与PCO2(A)变化趋势大体一致,也呈现出明显的昼夜变化和日际变化。pH、SIc和HCO3-变化趋势大体上与PCO2(A)相反,EC和Ca2+则无明显的昼夜变化,但存在一定的日际变化。随着旅游活动强度的增加,滴水水化学变化幅度逐渐增大。此外,不同滴水点所在洞腔结构、大小、封闭程度等不同,使PCO2(A)的扩散和通风程度存在差异,进而影响洞穴滴水水化学组分和洞穴沉积物沉积状况。因此,本研究对于洞穴环境保护和管理及其岩溶洞穴碳循环的研究具有重要意义。

    参考文献 | 相关文章 | 计量指标
    走廊南山河流纵剖面高海拔裂点的成因
    陈苗,胡小飞,王维
    2018, 73 (9):  1702-1713.  doi: 10.11821/dlxb201809007
    摘要 ( 800 )   HTML ( 19 )   PDF(7112KB) ( 657 )  

    河流水力侵蚀物理模型表明基岩河道纵剖面在均衡状态时表现为平滑上凹的形态,其特征反映了构造、基岩抗侵蚀能力和气候的作用;然而自然界河道纵剖面多呈现以裂点为特征的不均衡形态,不均衡的剖面形态以及裂点的研究同样可以对外力作用的变化起到很好的指示作用。位于北祁连的走廊南山高海拔河道纵剖面普遍呈现不均衡形式且发育海拔较高的裂点。通过对裂点成因分析发现,这些裂点并不主要受控于岩性、气候、构造等因素,而反映了冰川作用遗留地形与河流地形的分界。这一结果说明在对河道纵剖面高海拔裂点进行分析时要考虑到古冰川遗留地形也会对现代河道纵剖面产生重要影响,为进一步认识和理解造山带地貌演化以及控制因素提供了思路。

    参考文献 | 相关文章 | 计量指标
    三峡水库蓄水后荆江洲滩变化特征
    薛兴华,常胜,宋鄂平
    2018, 73 (9):  1714-1727.  doi: 10.11821/dlxb201809008
    摘要 ( 571 )   HTML ( 18 )   PDF(5741KB) ( 772 )  

    目前对三峡水库蓄水后荆江河段的洲滩演变还缺乏完整认识。基于三峡水库蓄水前后枯水期遥感影像,分析了荆江洲滩的冲淤变化与分布及形态演变。结果表明,蓄水后荆江洲滩总面积持续冲刷减小,累计冲刷4.56 km2,大部分发生在蓄水后前6年(冲刷速率0.55 km2/a)。上、下荆江洲滩的冲淤演变存在差异性。上荆江洲滩总面积一直处于冲刷萎缩中,且其强度明显大于下荆江,累计冲刷6.46 km2;下荆江前期(2002-2009年)冲刷、后期(2009-2015年)淤增,累计淤1.90 km2。在冲淤分布上,上荆江凸岸滩持续冲刷萎缩,凹岸滩前期冲刷、后期略有淤增,心滩(洲)前期淤积增长而后期冲刷萎缩;下荆江主要是凸岸滩冲刷,凹岸滩和心滩(洲)有所发展。根据滩体位置活动和冲淤动态性,荆江心滩(洲)演变被划分为8种典型类型。在形态演变上,上荆江以凸岸突出滩体和边滩发育的凹岸滩冲刷显著,形态变化不大。下荆江凸岸滩上游弯侧冲蚀后退、湾顶退缩、下游弯侧淤积伸长,形态趋向低弯扁平化,在高弯曲特定河湾平面形态格局下凸岸冲刷—淤积过程延伸到相邻河湾凹岸,成为下荆江凹岸滩和心滩淤积发展的重要因素,但淤积一般不越过凹岸湾顶。

    参考文献 | 相关文章 | 计量指标
    金沙江金江街段河流阶地年代及对河谷水系演化历史的启示
    董铭,苏怀,史正涛,明庆忠,何回丽
    2018, 73 (9):  1728-1736.  doi: 10.11821/dlxb201809009
    摘要 ( 665 )   HTML ( 20 )   PDF(3097KB) ( 945 )  

    100多年来,关于金沙江独特水系格局的形成历史一直是地学界争论的重要话题之一。多数学者认为,现代金沙江水系是古长江袭夺古红河上游发展过来的。红河海底扇5.5 Ma泥沙供给中断被认为与这一袭夺事件有关。然而,长期以来人们一直没有找到与这一时代相匹配的地貌证据。最近在金沙江金江街段找到了多达8级的河流阶地序列,ESR测年结果显示这些阶地的形成年代为1.07 Ma、0.70 Ma、0.65 Ma、0.51 Ma、0.47 Ma、0.44 Ma、0.30 Ma和0.18 Ma,结合GPS高程测量数据,推算最近1.0 Ma以来的河谷平均下切速率为147 mm/ka。以填充河谷地形为主要手段的古地形恢复结果(基于DEM数据)显示,古长江袭夺古红河上游形成现代金沙江水系发生在这一区域内海拔2000 m左右的古地形面解体之后,依照河谷平均下切速率外推,古地形面解体时代为5.5 Ma,即现代金沙江水系形成于5.5 Ma之后。我们的研究结果与红河海底扇的资料形成一个相互呼应的证据链,为重建现代金沙江水系格局形成历史提供重要依据。

    参考文献 | 相关文章 | 计量指标
    基于洪涝情景模拟的城市公共服务灾害应急响应空间可达性评价——以医疗急救为例
    殷杰,许世远,经雅梦,尹占娥,廖邦固
    2018, 73 (9):  1737-1747.  doi: 10.11821/dlxb201809010
    摘要 ( 749 )   HTML ( 24 )   PDF(9899KB) ( 638 )  

    在气候变化和城市化背景下,日益频发的洪涝灾害业已成为当前中国城市面临的重大挑战,也是城市灾害管理和应急响应研究的热点问题之一。本文旨在构建一套适用于城市尺度的洪涝灾害应急响应能力评估范式,以上海市外环以内中心城区为研究区,采用洪涝数值模拟与GIS网络分析相结合的研究方法,评估了当前以及未来正常条件和不同重现期洪涝情景下,城市关键公共服务部门(120急救)应急响应的空间可达性。结果显示:洪涝淹没强度(范围和水深)、道路交通状况(车流速度)以及应急服务机构的数量和位置共同决定了城市洪涝灾害医疗急救的服务范围及响应时间;由于上海市中心城区洪涝影响范围主要分布在黄浦江两岸2~3 km以内区域,因此洪涝对整个中心城区应急医疗服务的影响有限,主要是位于滨江地区部分医院的应急响应范围较正常状态下明显减少,120急救车辆无法或延迟达到部分救援点。研究表明基于洪涝情景模拟的城市公共服务应急响应空间可达性定量化评估方法,具有重要的科学价值和实践意义,可为中国城市洪涝灾害应急管理提供决策依据。

    参考文献 | 相关文章 | 计量指标
    气候变化与生态环境
    中国气温变化对全球变暖停滞的响应
    杜勤勤,张明军,王圣杰,车存伟,邱雪,马转转
    2018, 73 (9):  1748-1764.  doi: 10.11821/dlxb201809011
    摘要 ( 934 )   HTML ( 45 )   PDF(8243KB) ( 1402 )  

    1998-2012年出现的全球变暖停滞(global warming hiatus)现象,近年来受到各界的广泛关注。基于中国622个气象站的气温数据,研究了全国及三大自然区气温变化对全球变暖停滞的响应。结果表明:① 1998-2012年间,中国气温变化率为-0.221 ℃/10 a,较1960-1998年增温率下降0.427 ℃/10 a,存在同全球变暖停滞类似的增温减缓现象,且减缓程度更明显,其中冬季对中国增温减缓的贡献最大,贡献率为74.13%,夏季最小;② 中国气温变化对全球变暖停滞的响应存在显著的区域差异,从不同自然区看,1998-2012年东部季风区和西北干旱区降温显著,其中东部季风区为中国最强降温区,为全国增温减缓贡献了53.79%,并且具有显著的季节依赖性,减缓期冬季气温下降了0.896 ℃/10 a,而夏季上升了0.134 ℃/10 a。青藏高寒区1998-2012年增温率达0.204 ℃/10 a,对全球变暖停滞的响应并不显著;③ 中国增温减缓可能受太平洋年代际振荡(PDO)负相位、太阳黑子数与太阳总辐照减小等因素的影响;④ 1998-2012年中国虽出现增温减缓现象,但2012年之后气温快速升高,且从周期变化看,未来几年可能持续升温。

    参考文献 | 相关文章 | 计量指标
    五华河流域非点源污染风险区和风险路径识别
    陈裕婵,张正栋,万露文,张杰,杨传训,叶晨,李青圃
    2018, 73 (9):  1765-1777.  doi: 10.11821/dlxb201809012
    摘要 ( 599 )   HTML ( 12 )   PDF(4914KB) ( 639 )  

    非点源污染是亟待解决的水环境问题之一,确定非点源污染过程中的潜在风险区和风险路径是解决非点源污染问题的关键。引入景观生态学中的最小累积阻力模型,以高污染负荷的耕地、建设用地为“源”,运用表示下垫面产流性质的地形湿度指数和CN值构建阻力面,对五华河流域的非点源污染高风险区和风险路径进行可视化识别和分析。结果显示:① 在下垫面产流作用下,五华河流域低产流区主要分布在流域西南部,高产流区呈“人”形贯穿全流域。② 最小累积阻力模型能够有效识别流域内非点源污染风险区和风险路径,五华河流域非点源污染高风险区分布于迥龙、田心、龙母、铁场、登云、通衢、鹤市、紫市、岐岭、华城、转水、潭下、水寨等地的河道两岸,以耕地为“源”的非点源污染风险路径与建设用地为“源”的非点源污染风险路径在空间分布上差异较大。③ 耕地对五华河水质的影响大于建设用地对五华河水质的影响,耕地中的富营养物质和沉积物更容易随地表径流进入受纳水体。④ 流域尺度上治理以耕地为“源”的非点源污染应在邻近耕地的河流两岸建立一定宽度的植被缓冲区,治理以建设用地为“源”的非点源污染宜围绕关键源区进行治理。本研究为非点源污染风险区和风险路径的识别提供一种新的思路,为进一步开展非点源污染治理提供理论依据。

    参考文献 | 相关文章 | 计量指标
    基于卫星遥感和再分析数据的青藏高原土壤湿度数据评估
    范科科,张强,史培军,孙鹏,余慧倩
    2018, 73 (9):  1778-1791.  doi: 10.11821/dlxb201809013
    摘要 ( 831 )   HTML ( 16 )   PDF(4967KB) ( 926 )  

    土壤水是地表与大气在水热交换方面的关键纽带,是关键水循环要素,更是地表产汇流过程的关键控制因子。青藏高原是地球第三极,也是亚洲水塔,探讨青藏高原土壤水变化对于探讨青藏高原热力学特征变化及其对东亚乃至全球气候变化的影响具有重要意义,而获取高精度长序列大尺度土壤水数据集则是其关键。本文利用青藏高原100个土壤水站点观测数据,从多空间尺度(0.25°×0.25°,0.5°×0.5°,1°×1°)、多时间段(冻结和融化期)等角度,采用多评价指标(R、RMSE、Bias),对多套遥感反演和同化数据(ECV、ERA-Interim、MERRA、Noah)进行全面评估。结果表明:① 除ERA外,其他数据均能反映青藏高原土壤水变化,且与降水量变化一致。而在那曲地区,遥感反演和同化数据均明显低估实测土壤水含量。从空间分布来看,MERRA和Noah与植被指数最为一致,可很好地反映土壤水空间变化特征;② 青藏高原大部分地区土壤水变化主要受降水影响,其中青藏高原西部边缘与喜马拉雅地区土壤水变化则受冰雪融水和降水的共同影响;③ 除阿里地区外,大部分遥感反演和同化数据在融化期与实测土壤水相关性高于冻结期,其中在那曲地区,遥感反演和同化数据均高估冻结期土壤含水量,却低估融化期土壤含水量。另外,遥感反演和同化数据对中大空间尺度土壤水的估计要好于对小空间尺度土壤水的估计。本研究为青藏高原土壤水研究的数据集选择提供重要理论依据。

    参考文献 | 相关文章 | 计量指标
    湖南省地形因素对水稻生产的影响
    王琛智,张朝,张静,陶福禄,陈一,丁浒
    2018, 73 (9):  1792-1808.  doi: 10.11821/dlxb201809014
    摘要 ( 719 )   HTML ( 28 )   PDF(8543KB) ( 1019 )  

    作为一种人工生态系统,水稻的生产受自然条件和人为管理的共同影响。其中,地形因素不仅决定了水稻生长必须的光热等资源的再分配,而且还影响人类的农业管理活动。以中国水稻主产区中地形较为复杂的湖南省为研究对象,基于多种数据,从生态学的“格局—过程”思想出发,讨论了地形因素对水田分布、水稻生长过程及水稻产量的影响。研究发现,地形因素对水田分布的空间格局影响最为明显,对水稻生长过程影响次之,而对水稻产量影响较弱。在各地形因素中,高程对于水稻生产的影响最为突出。研究结果还反映出地形因素造成了湖南省水田分布和热量分配在北部平原地区的不匹配,政府部门应根据区域特征,因地制宜采取措施,提高该地区水稻生产效益。

    参考文献 | 相关文章 | 计量指标
    基于贝叶斯网络的水源涵养服务空间格局优化
    曾莉,李晶,李婷,杨晓楠,王彦泽
    2018, 73 (9):  1809-1822.  doi: 10.11821/dlxb201809015
    摘要 ( 630 )   HTML ( 26 )   PDF(3325KB) ( 841 )  

    以渭河流域关中—天水经济区段(简称“渭河流域关天段”)为研究区,基于贝叶斯网络和水量平衡原理建立了水源涵养服务网络模型;将CA-Marcov模型与贝叶斯网络模型相结合,预测了2050年不同土地利用情景及其水源涵养服务分布概率;提出了关键变量关键状态子集方法,对研究区水源涵养服务空间格局进行优化。结果表明:① 保护情景下,林地面积增长了18.12%,其主要来源为耕地;草地和城市面积增长缓慢,分别增加了0.73%和0.38%;水体和未利用地分别减少了5.08%和0.92%,该情景下的水源涵养量值偏高的概率在3种情景中最大,保护情景的设计对未来的土地利用政策制定具有一定参考价值。② 水源涵养服务的关键影响因子是降水、蒸散发和土地利用,水源涵养量最高状态对应的关键变量关键状态子集是:﹛降水= 1,蒸散发= 2,土地利用= 2﹜,该子集主要分布在年平均降雨量和蒸散发量较大,植被覆盖率高的地区。③ 研究区适宜优化水源涵养的区域主要分布在天水市麦积区南部、宝鸡市陇县西南部和渭滨区南部、咸阳市旬邑县东北部和永寿县西北部,以及铜川市耀州区西部。结合贝叶斯网络模型研究水源涵养服务的优化区域,不仅有助于提升对生态系统水源涵养服务过程的直观认识,而且增加了情景设计和格局优化的合理性。在此基础上提出的关键状态关键因子方法,对研究区水源涵养生态环境建设和政策制定都具有重要的意义。

    参考文献 | 相关文章 | 计量指标