黄河流域旅游经济系统韧性的时空差异特征及其影响因素.pdf

第 卷第期 年月地 理 与 地 理 信 息 科 学 收稿日期：；修回日期：基金项目：国家自然科学基金项目（）作者简介：王兆峰（），男，博士，教授，博士生导师，主要从事旅游管理等研究。：黄河流域旅游经济系统韧性的时空差异特征及其影响因素王 兆 峰，张 先 甜（湖南师范大学旅游学院，湖南 长沙 ）摘要：探析旅游经济系统韧性状态对推进旅游高质量发展具有科学参考价值。该文以黄河流域 个地级市为研究单元，采用核密度估计、变差函数和空间计量模型，揭示 年黄河流域旅游经济系统韧性水平的时空差异特征及影响因素。结果表明：黄河流域旅游经济系统韧性水平整体偏低，时序上呈现出先升后降的阶段性特征，空间上存在显著的正相关性，呈现出“下游中游上游”的块状分布格局。旅游经济系统韧性空间变异幅度不断缩小，由空间自相关引起的结构化分异不断增强；韧性水平的空间差异主要体现在局部微观尺度，东西方向空间差异最大，南北方向空间差异最小。黄河流域旅游经济系统韧性存在显著的空间溢出效应，财政自给水平、市场规模、旅游专业化、信息化水平对黄河上、中、下游韧性的提升具有直接促进作用，信息化水平、市场规模对韧性水平的提高存在正向溢出效应。关键词：旅游经济系统；韧性；时空差异；驱动机制；黄河流域中图分类号：文献标识码：文章编号：（）引言旅游业是跨行业、跨地区的综合性产业，具有辐射面宽、融合性强、影响力大的特点，诸多突发性事件和社会经济环境的衰退均会对旅游业的稳定发展产生直接冲击。已有研究通过探讨旅游业与气候变化、自然灾害的关系，以及政治不稳定、冲突、恐怖主义等对旅游业的影响，论证了旅游业隐含的脆弱性特征。当前，面对世界大变局和日益复杂的外部环境，旅游业自身的敏感性、脆弱性等特征加速显现，“韧性（）”作为表述或刻画事物应对不确定性因素的重要属性，为深刻理解和剖析旅游业如何应对内外部冲击提供了良好的研究视角。韧性理论最早应用于物理学科，而基于旅游学角度的旅游韧性定义为：“旅游系统应对压力、冲击的抵御能力及恢复其增长路径的能力”。经济属性是旅游业的本质属性，旅游经济系统运转需要面对区域内外和行业之间的多重干扰，并最终生存于系统对于外界干扰的敏感脆弱性与应对能力间的制衡博弈之中，即表现出旅游经济系统韧性 由旅游经济系统受到干扰后恢复到初始或更理想状态的能力，可衡量旅游经济可持续发展水平，。目前国内外关于旅游经济韧性已有一定成果，从研究方法上看，国内学者采用本底趋势线法、集成分析法、模型 等分析旅游经济抗周期性和后危机反弹能力，国外学者通过构建韧性模型、恢复力模型衡量旅游经济韧性 ；从测度方法上看，早期学者多采用单一指数法，如根据旅游人数、旅游收入、旅游就业 等要素的变化情况衡量区域旅游经济韧性水平，目前越来越多的学者开始采用多维指标法 ，通过构建综合性指标体系量化旅游经济韧性。总体而言，国内外旅游经济系统韧性研究尚处于探索阶段，有关旅游经济韧性时空演化的量化研究依然薄弱，在理论评述和实证剖析过程中也缺乏较系统的驱动机制探究；此外，既有定量研究以全国省域尺度居多，关于流域尺度的研究颇少，而旅游业是区域经济的重要组成部分，旅游经济韧性更是关乎流域内外甚至是中国旅游业的可持续发展。黄河流域作为我国区域旅游经济开发的重要单元，已成为国家“十四五”旅游发展规划的重要篇章，然而，目前黄河流域旅游经济面临着整体发展能级偏低、区域内部发展差距拉大、产业转型压力加剧、发展内生动力不足等问题，阻碍了黄河“黄金旅游线”价值的最大化释放，制约了黄河沿线旅游经济健康可持续发展。因此，本文从流域视角切入，对 年黄河流域旅游经济系统韧性的时空差异特征及驱动机制进行探析，以期为黄河流域旅游经济运行与政策仿真、旅游危机管理与影响评估提供指导依据，从而提高黄河流域旅游经济系统韧性，实现黄河流域旅游高质量发展。研究方法与数据来源 研究区域黄河流域级旅游景区占全国的 ，是我国旅游资源最富集的区域之一，但 年黄河流域全年旅游总收入仅占全国的 ，旅游资源赋存未能有效转化为旅游经济优势，导致区域旅游经济发展水平整体不高；受新冠疫情影响，年黄河流域旅游总收入下滑近 ，流域内旅游经济系统抵御危机、应对风险显化的能力问题日渐凸出，因此，强化区域旅游经济系统韧性能力至关重要。综合考虑沿黄省市行政区划的完整性和区域经济社会高质量发展与黄河流域的关联性，本文以黄河自然流经地区为主体，排除已纳入长江经济带的四川省和广义上属于东北地区的内蒙古蒙东地区，兼顾数据的可得性和准确性，最终从青海、甘肃、陕西、宁夏、内蒙古、山西、河南和山东省域共 个市域中选取 个市级行政区作为研究区域，同时，为便于进行对比和细化分析，借鉴文献 并参考 黄河年鉴，以内蒙古河口镇和河南省桃花峪为上、中游和中、下游的分界线，将研究区域划分为黄河流域上、中、下游（图）。注：审图号为 （），底图无修改，下同。图 黄河流域研究区域示意犉 犻 犵 犛 犽 犲 狋 犮 犺犿 犪 狆狅 犳 狋 犺 犲 狊 狋 狌 犱 狔犪 狉 犲 犪 犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀 评价指标体系构建目前，关于旅游经济系统韧性尚无统一的评价模型，提出了一种包含抵抗力、恢复力、调整力、创造力个维度的经济韧性分析框架，已在韧性研究领域获得一定共识。因此，本文基于 的理论框架，充分借鉴文献 ，的研究成果，并结合黄河流域的旅游经济实际，从敏感与脆弱性、抵抗与恢复力、适应与调整力、创新与转型力个维度选取 个指标，构建黄河流域旅游经济系统韧性评价指标体系（表）。敏感与脆弱性即系统受自身特征的制约和对内外扰动的敏感性导致结构、功能容易受损的一种内在属性，旅游收入增长弹性、旅游外汇依存度、旅游地承载力等均是旅游经济系统所表现出的敏感与脆弱性。其中，游客密度指数采用旅游人数与地区常住人口总数的比值表征，旅游收入增长弹性用旅游总收入增长率与 增长率的比值表征。抵抗与恢复力即系统抵御冲击、吸收冲击以恢复到原始动态平衡的能力，旅游市场及旅游经济水平、地方经济及产业结构基础等是影响抵抗与恢复力的关键因素。适应与调整力即系统通过重新整合内外部资源并调整自身结构适应冲击后的环境来维持自身稳定发展的能力，旅游业内外部资源包括旅游人力资源条件、旅游交通条件、旅游设施条件、旅游接待条件等方面。其中，旅游资源丰度指数采用对 级以上旅游景区赋分表征，游客空间指数用游客数量与地区总面积的比值表征。创新与转型力即区域旅游经济系统可以恢复到破坏前状态甚至达到更高状态的能力，人表 黄河流域旅游经济系统韧性评价指标体系犜 犪 犫 犾 犲 犈 狏 犪 犾 狌 犪 狋 犻 狅 狀 犻 狀 犱 犲 狓 狊 狔 狊 狋 犲 犿狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 狅 犳 狋 犺 犲 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀目标层准则层指标层属性旅游经济系统韧性敏感与脆弱性旅游收入增长弹性 负旅游外汇收入 旅游总收入 负城镇失业登记率 负游客密度指数 负旅游总收入 正抵抗与恢复力国内旅游收入 万元正旅游外汇收入 百万美元正国内旅游人次 千人正入境旅游人次 千人正住宿餐饮营业收入 千元正地区生产总值 亿元正居民消费水平 元正产业结构高级化 正适应与调整力第三产业从业人口数量 人正旅游资源丰度 个正星级饭店数量 家正路网密度 正人均汽车占有量 辆正游客空间指数 （人 ）正对外开放水平 正建成区公共绿地率 正旅游固定资产投资 万元正创新与转型力城镇化率 正旅游院校数量 所正教育支出占 比重 正专利授权数 项正页第第期 王兆峰，张先甜：黄河流域旅游经济系统韧性的时空差异特征及其影响因素才、教育、创新、城镇化等是系统创新与转型不可或缺的因素。为消除旅游经济系统韧性数据量纲的不一致性与极端值问题，对指标进行极差标准化处理；在利用熵值法计算各指标权重的基础上，采用综合指数法测算黄河流域旅游经济系统韧性水平。研究方法 核密度估计核密度估计作为一种非参数估计方法，通过考虑空间因素能更全面地衡量黄河流域旅游经济系统韧性水平的时序变化，本文选择常用的 核函数进行估计，计算公式为：犵（狓）狀 犿狀犻 狕（狓犻珚狓犿）（）狕（狌）（狌）槡（）式中：犵（狓）为随机变量狓的密度函数，狀为观测数量，珚狓为观测值均值，犿为带宽，狕（狌）为核函数（即权重函数），狌为样本均值。分析采用 方法中的全局 测定黄河流域旅游经济系统韧性水平的空间聚集程度，采用 犌指数进一步探究系统韧性在局部空间上的依赖性及异质性。空间变差函数为描述区域化变量变异性和结构性，本文使用空间变差函数（犺）对各地理单元空间格局的变异程度进行测度，计算公式为：（犺）犖（犺）犖（犺）犻 犣（狓犻）犣（狓犻犺）（）（犺）犺 犇（）式中：犣（狓犻）和犣（狓犻犺）分别为犣（狓）在狓犻和（狓犻犺）空间位置上的属性值（犻，犖（犺），分维数犇表示变异函数曲线率，其值越接近，空间均衡性特征越显著。数据来源本文原始数据来源于 年 中国城市统计年鉴、国家统计局数据库网站（：）、年黄河流域各地级市统计年鉴以及国民经济和社会发展统计公报、文化与旅游厅官方网站等，需要说明的是，由于 年莱芜市并入济南市，为保证前后核算口径的一致性，将 年莱芜市相关数据汇总至济南市，其他缺失数据均通过插值法补充。黄河流域旅游经济系统韧性时空格局 黄河流域旅游经济系统韧性水平时序格局由图 可知，黄河流域旅游经济系统韧性水平不高，但整体呈缓慢上升趋势，由期初的 上升至期末的 ；从地带层面看，黄河流域旅游经济系统韧性存在异质性特征，年流域内呈现“下游中游上游”的分异格局，其原因在于，此阶段支持旅游发展的各项政策密集出台，推动旅游经济发展措施落地实施，致使旅游产业规模持续扩大，旅游投入要素相应增加，旅游产品体系不断完善，旅游经济成效凸显，旅游经济韧性水平攀升；年下游和整体水平表现出下降趋势，中、上游地 区 变 动 较 平 稳，各 区 域 涨 幅 呈 现“上 游（）中游（）下游（）”趋势，这可能是由于下游城市的资源体量不大、旅游承载力有限、旅游资源禀赋不足的现实情况与大量的旅游固定资产投资不匹配，致使旅游经济系统内外部资源、环境、需求等各种矛盾逐渐凸显，整体旅游经济韧性水平下降。虽然平均值能体现黄河流域旅游经济韧性整体状况及变化趋势，但会掩盖市域间差距，因此，进一图 年黄河流域旅游经济系统韧性水平及核密度分布犉 犻 犵 犚 犲 狊 犻 犾 犻 犲 狀 犮 犲 犾 犲 狏 犲 犾 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犪 狀 犱 犻 狋 狊 犽 犲 狉 狀 犲 犾 犱 犲 狀 狊 犻 狋 狔犱 犻 狊 狋 狉 犻 犫 狌 狋 犻 狅 狀 犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀 犳 狉 狅 犿 狋 狅 页第地 理 与 地 理 信 息 科 学第 卷步借助核密度估计图（图）判断韧性水平的动态演进轨迹及趋势特征。由图 可以看出：在位置变化上，核密度曲线波峰所处位置均在左侧且曲线右侧“长尾”现象明显，表明黄河流域旅游经济系统韧性水平大部分处于中低状态，但整体水平在缓慢抬升；在分布形态上，密度曲线主峰高度不断下降，由“尖峰型”向“宽峰型”小幅转变，说明离散化程度逐渐缩小，地区绝对差异有所减小，未来要对黄河流域旅游经济系统韧性建设能力进一步优化，缩小地区之间的发展差距。黄河流域旅游经济系统韧性水平空间格局从图可知，黄河流域旅游经济系统韧性高水平和较高水平区域主要围绕西安、郑州、济南等省会城市以及青岛、洛阳、烟台、威海等省重点发展城市分布，在一定程度上说明黄河流域旅游经济韧性建设可以这类城市为增长极带动周边地区乃至全省域发展；在空间分布上，黄河中游和下游城市韧性水平整体高于上游地区，总体分布格局变动较小，济南、郑州、西安等省会城市的韧性水平指数在个时间剖面上分别增至 、以上，位于所有城市前列；上游高水平地区主要包括兰州和银川，而甘肃的平凉、定西和宁夏的中卫、固原等城市韧性水平低于 且增幅较小；中游地区整体增幅较大，如陕西的宝鸡、汉中、咸阳及山西的临汾、吕梁等。研究期间黄河流域旅游经济系统韧性水平“中上游高、下游低”的整体格局变化甚微。图 黄河流域旅游经济系统韧性水平空间格局犉 犻 犵 犛 狆 犪 狋 犻 犪 犾 狆 犪 狋 狋 犲 狉 狀狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 犾 犲 狏 犲 犾 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀 黄河流域旅游经济系统韧性空间关联格局 全局空间关联格局借助 软件计算 年黄河流域旅游经济系统韧性水平的 值（表），分析可知：年 值均大于，且通过 水平的显著性检验，表明各城市旅游经济系统韧性水平之间存在空间正相关关系；值伴有“增大减小”的波动趋势，表明韧性水平的空间差异具有阶段性特征，其中，年 值呈增大态势，空间聚集特征不断增强，年 值递减，表明韧性水平的空间差异性有所减小，黄河流域各城市韧性水平之间的正向空间聚集态势由强转弱。表 黄河流域旅游经济系统韧性水平的全局 犜 犪 犫 犾 犲 犌 犾 狅 犫 犪 犾犕 狅 狉 犪 狀 狊 犐 狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 犾 犲 狏 犲 犾 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀年份 犘值犣得分 页第第期 王兆峰，张先甜：黄河流域旅游经济系统韧性的时空差异特征及其影响因素 局域空间关联格局借助 软件分析 犌犻统计量，得到黄河流域旅游经济系统韧性水平冷热点分布（图），发现黄河流域旅游经济系统韧性水平在个时间点均表现出“上游偏冷，下游偏热，中游不显著”的特点，具有显著的“块状”分布特征，且上游冷点区数量逐年减少，上游冷点区和次冷点区逐渐向四周扩散，呈现“中心离散，周围聚集”的特点。图 黄河流域旅游经济系统韧性水平冷热点分布格局犉 犻 犵 犇 犻 狊 狋 狉 犻 犫 狌 狋 犻 狅 狀狆 犪 狋 狋 犲 狉 狀狅 犳 犮 狅 犾 犱犪 狀 犱犺 狅 狋 狊 狆 狅 狋 狊 犳 狅 狉 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 犾 犲 狏 犲 犾 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀 黄河流域旅游经济系统韧性空间分异格局为进一步解析旅游经济系统韧性的空间差异演化格局，将黄河流域旅游经济系统韧性水平值作为空间变量赋予各单元的几何中心点，并设置采样步长为 ，计算变差函数及个方向的分维数，进一步分析空间分异情况。）由表可知，块金值和基台值均呈下降趋势，块金系数先增至 后减至 ，说明黄河流域旅游经济系统韧性空间差异不断缩小，且空间分异的随机成分开始减弱，而由空间自相关造成的空间分异态势趋于显著；模型拟合系数犚均大于 ，拟合效果较好，说明韧性水平空间格局具有较好的连续性和稳定性，空间自组织性整体相对较好；年旅游经济系统韧性水平变程参数在既定步长下先减后增，说明在结构化空间梯度引起的空间关联范围趋于缩小，高韧性水平区逐渐表 黄河流域旅游经济系统韧性变差函数拟合参数犜 犪 犫 犾 犲 犘 犪 狉 犪 犿 犲 狋 犲 狉 狊 狅 犳 狏 犪 狉 犻 狅 犵 狉 犪 犿狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀年份变程块金值基台值块金系数 拟合模型犚 靠近热点极核，各区域空间关联作用范围不断扩大，高韧性水平区域的涓滴作用和辐射作用开始显现。）由表可知，旅游经济系统韧性水平全方向分维值呈先降后升变化，说明韧性水平空间差异在研究期末具有缩小态势。其中，东西方向上的分维数始终较小，偏离理想均值较远，表明该方向上的韧性水平空间差异始终最大，但 年该方向差异性出现缩小趋势并具有明显的迁移现象；南北、东南西北方向的分维数值相对较高且整体稳定，南北方向的韧性水平空间差异最小；东北西南方向分维数由 减至 ，犚由 增至 ，因此，该方向可能会成为黄河流域旅游经济系统韧性水平空间差异的主要方向。）由图可知，年黄河流域旅游经济系统韧性水平的空间差异格局存在一定规律性和延续性特征，韧性水平整体表现为东高西低的分布格局，空间差异形态由显著的梯度变化态势逐渐向较为均衡的结构性分布转变。黄河流域旅游经济系统韧性水平均衡性的提升一方面得益于上游西部区域经济、文化环境的进步，突破了旅游发展瓶颈，旅游经济发展后劲增强，旅游经济韧性建设取得重大成果，另一方面是由于“沿黄九省（区）黄河之旅旅游联盟”“中国黄河沿线城市与旅游产业联盟”等组织页第地 理 与 地 理 信 息 科 学第 卷表 黄河流域旅游经济系统韧性变差函数分维数犜 犪 犫 犾 犲 犉 狉 犪 犮 狋 犪 犾 犱 犻 犿 犲 狀 狊 犻 狅 狀狅 犳 狏 犪 狉 犻 狅 犵 狉 犪 犿狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀年份全方向南北东北西南东西东南西北犇犚犇犚犇犚犇犚犇犚 图 黄河流域旅游经济系统韧性变差函数犉 犻 犵 犞 犪 狉 犻 犪 狋 犻 狅 狀 犳 狌 狀 犮 狋 犻 狅 狀犱 犻 犪 犵 狉 犪 犿狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀的建立，构建了黄河流域旅游经济合作共赢的联盟新机制，充分发挥了韧性高值区对整体流域旅游经济系统韧性水平的带动作用。黄河流域旅游经济系统韧性的影响因素 指标选取为厘清黄河流域旅游经济系统韧性的影响因素及可能存在的溢出效应，参考王倩等的研究成果 ，基于数据的可得性，在单位根检验（均拒绝了变量不平稳的原假设）和影响因子共线性检验（均小于）的基础上进行指标遴选，具体如下：旅游专业化（）：旅游专业化能反映区域旅游经济运行的区位优势，区域旅游专业化水平越高，资源消耗及产业成本越低，旅游经济系统产出效率越高，从而影响旅游经济系统韧性水平，用旅游总收入的区位熵衡量。财政自给水平（）：财政自给水平的提高可以增加城市经济发展的内生活力，并产生强大的集聚经济优势，为旅游经济的平稳运行提供物质财富支持，从而提高旅游经济系统的抗风险能力，用地方财政收入与财政支出的比值表征。环境规制（）：环境规制可以反映政府对环境的管控力度，当环境规制的成本比重增加时，一方面可能会增加区域经济负担，另一方面可能会倒逼旅游企业采用先进技术和管理方法以提高旅游资源利用效率，因此，环境规制对旅游经济系统韧性的影响具有不确定性，还需进一步检验，参考叶琴等的研究，计算城市各类污染物排放量的综合指数以表征城市环境规制强度。市场规模（）：旅游市场作为旅游业发展的重要保证，直接影响旅游资源的分配及设施的使用，区内市场规模越大，旅游经济运行所需的人力资源越丰富，旅游经济发展的潜力越大，用地方常住人口数量取对数表征。信息化水平（）：旅游信息化通过旅游大数据的整合，激活智慧旅游、虚拟旅游等旅游活动要素，促进技术外溢和先进管理理念的传播，推动旅游经济增长方式转变，因此也与区域旅游经济系统韧性息息相关，本文采用互联网宽带接入用户数衡量信息化服务程度，用邮电业务量衡量信息传递服务能力，用页第第期 王兆峰，张先甜：黄河流域旅游经济系统韧性的时空差异特征及其影响因素二者综合指数表征信息化水平。模型构建空间计量模型主要包括空间滞后、空间误差和空间杜宾模型种，均可通过设定为一般化的空间模型的不同参数推导而来，其形式为：狔犻 狋 狔犻，狋 犻狔狋犡 犻 狋 犻犡狋狌犻狋犻 狋（）式中：犻狔狋和 犻犡狋为空间滞后项，、为回归系数，犻为空间权重行向量，狌犻为个体效应，狋为时间效应，犻 狋为随机扰动项。首先利用 检验判断建立空间计量模型的合理性和对模型的选择性问题（表），及 均通过了 水平的显著性检验，验证了运用空间面板数据模型的合理性；同时，检验值拒绝随机效应的原假设，采用固定效应模型；再次，运用、检验以及犚和 等指标对 、和 进行检验，结果证明采用空间面板杜宾模型更合理（表）。模型具体表达式为：犾犻 狋 犠犻 犼犾犻 狋犡犻 狋 犽犽犠犻 犼犡犻 狋 犽犽狌犻狋犻 狋（）式中：犾为韧性水平，犡分别为 、等解释变量，犠为空间邻接权重矩阵，为空间自相关系数，犽为样本量，犠犻 犼犾犻 狋为被解释变量滞后项，犠犻 犼犡犻 狋 犽为解释变量滞后项。表 空间杜宾模型估计结果犜 犪 犫 犾 犲 犈 狊 狋 犻 犿 犪 狋 犻 狅 狀 狉 犲 狊 狌 犾 狋 狊 狅 犳 犛 犇 犕犿 狅 犱 犲 犾变量（）空间固定效应（）直接效应（）溢出效应（）总效应 （）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）犠 （）犠 （）犠 （）犠 （）犠 （）（）（）犚 空间 ；空间 ；犖 注：、分别表示在 、水平下显著（括号内为标准误），下同。结果分析运用 软件对空间杜宾模型进行回归分析，根据对比结果，选择拟合效果最好的个体固定效应模型进行分析（表）。综合表、表可知，黄河流域旅游经济系统韧性水平的变化是在多因素共同作用下的复杂过程（图），具体而言：表 分流域空间杜宾模型估计结果犜 犪 犫 犾 犲 犈 狊 狋 犻 犿 犪 狋 犻 狅 狀狉 犲 狊 狌 犾 狋 狊 狅 犳 犛 犇犕犿 狅 犱 犲 犾 犻 狀 狊 狌 犫 狑 犪 狋 犲 狉 狊 犺 犲 犱变量（）黄河上游（）黄河中游（）黄河下游 （）（）（）（）（）（）（）（）（）（）（）（）（）（）（）犠 （）（）（）犠 （）（）（）犠 （）（）（）犠 （）（）（）犠 （）（）（）（）（）（）（）（）（）犖 犚 ）旅游专业化（）。根据黄河流域总体空间杜宾模型测算结果，旅游专业化对旅游经济系统韧性水平具有显著的正向影响，影响效应为 ，表明旅游专业化水平的提升是培育旅游经济系统韧性能力的重要一环。旅游专业化的发展易产生规模效应，从而减少旅游经济流通成本，促进旅游产业结构合理化、高级化，使旅游经济系统本身起到“冲击吸收器”的作用，增强分散冲击带来的负面影响并调整恢复的能力。从各流域城市旅游经济韧性的影响效应看，旅游专业化对黄河中、下游韧性水平产生积极的促进作用，对上游的促进作用尚不显著，可能的原因是：上游存在着旅游布局不合理、旅游服务设施不完善、资源利用水平低等问题，旅游经济运行未形成完整的旅游产业链体系，因此难以充分发挥旅游专业化对旅游经济系统韧性的保障作用。）财政自给水平（）和环境规制（）。从黄河流域总体空间杜宾模型测算结果看，财政自给水平的估计系数为正但不显著，原因可能在于，多数页第地 理 与 地 理 信 息 科 学第 卷图 黄河流域旅游经济系统韧性驱动机制犉 犻 犵 犇 狉 犻 狏 犻 狀 犵犿 犲 犮 犺 犪 狀 犻 狊 犿狅 犳 狉 犲 狊 犻 犾 犻 犲 狀 犮 犲 狅 犳 狋 狅 狌 狉 犻 狊 犿犲 犮 狅 狀 狅 犿 犻 犮 狊 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犫 犪 狊 犻 狀黄河流域所在城市是发展资源“缺口”地区，旅游发展财政资金存在支出结构不均衡、使用效率不高、资金项目管理困难等现象，因此，从全流域看，该因素尚未有效发挥拉动效应；但从分流域影响情况看，财政自给水平对黄河流域上游和下游韧性水平的促进作用已经开始显现，说明探索建立财政自给机制是培育黄河流域建设旅游经济系统韧性的重要手段，未来更应该强化全流域财政支持政策。无论是从全流域还是从分流域空间杜宾模型测算结果看，环境规制对黄河全流域及上、中游的影响均不显著，这被认为是由环境规制的技术效应和成本效应共同作用造成的，表明短期内环境规制难以促进黄河流域旅游经济系统韧性水平提升。）市场规模（）。人口规模对旅游经济系统韧性水平具有显著的正向影响，人是旅游活动的主体，旅游业是存在明显人口规模效益的产业，人口数量是能直接影响潜在旅游市场规模的要素，因此与旅游经济系统韧性直接相关；同时，市场规模还具有正向空间溢出效应，随着地区经济发展水平的不断提升，熟悉本地区旅游目的地的旅游者增多，旅游休闲边际效益不断降低，往往会向周围邻近地区扩散，形成旅游流“外溢”，从而导致周边地区承接的旅游者人数和旅游消费增多，促进旅游经济系统韧性水平的提高。）信息化水平（）。信息化水平对全流域旅游经济系统韧性具有显著的正向影响及空间溢出效应，旅游业是关联性极强的产业，在智慧旅游如火如荼的背景下，旅游信息化水平的提高有利于不同目的地分流，促进区域客流的空间置换、旅游资源设施的跨区共享，同时作为信息密集型产业，信息化拓宽了旅游营销、旅游业管理、旅游市场化发展的内涵，从而提高整个流域旅游经济系统对危机风险的预测和感知能力，有利于建立完善的旅游危机或风险的预警机制，优化对抗逆境和风险的韧性能力，虽然黄河流域信息化建设已经取得很大进步，但在信息化浪潮的冲击中找到旅游经济新的增长点，也是增强黄河流域旅游经济系统韧性能力的当务之急。进一步采用反距离权重矩阵代替空间邻接权重矩阵进行稳健性检验，在同一模型、两种权重矩阵下，模型参数值虽有所不同，但其方向及显著性较一致，表明模型估计结果较为稳健，支持上述分析。结论与建议本文基于韧性视角，综合运用核密度估计、空间变差函数、空间计量等方法对黄河流域旅游经济系统韧性水平的时空差异及其影响因素进行探索，结论如下：从时空格局看，年黄河流域旅游经济系统韧性水平缓慢提高，各城市之间的差距具有缩小趋势；旅游经济系统韧性整体呈现出“下游中游上游”的空间分异格局，且高值区集中于省会城市或省重点发展城市；流域内韧性水平的空间相关性显著为正，局部冷热点分布不均衡，具有“上游偏冷，下游偏热，中游不显著”的块状分布特征。从区域差异看，研究期内黄河流域旅游经济系统韧性空间变异幅度不断缩小，由空间自相关引起的结构化分异不断增强。韧性水平在全方向上的差异较小，空间差异主要体现在局部微观尺度，其中，东西方向空间差异最大，南北方向空间差异最小。从影响因素看，黄河流域旅游经济系统韧性具有显著的空间溢出效应，韧性水平高的城市将发挥涓滴和辐射作用，带动周围地区韧性水平提高。当前，财政自给水平、市场规模、旅游专业化、信息化水平对黄河上、中、下游旅游经济韧性的提高具有直接促进作用，信息化水平、市场规模对韧性水平的提页第第期 王兆峰，张先甜：黄河流域旅游经济系统韧性的时空差异特征及其影响因素高存在显著的空间溢出效应。黄河流域上、中、下游旅游经济系统韧性水平表现出显著的区段性、差异性，应根据区域旅游经济发展阶段与旅游经济韧性实际情况，制定有针对性的提升策略：上游要积极转变旅游业增长方式，以平凉、定西等低水平城市为试点，培养高素质、创新型旅游人才，以独特的传统文旅资源开发为导向，使之成为黄河上游旅游经济的重要增长点，补足城市短板，推进区域协调发展，同时还要制定科学合理且操作性强的地方资助实施路径，为上游旅游经济系统韧性建设提供支撑；中游在推动黄河流域旅游可持续发展与韧性能力建设中起到承上启下的作用，必须优化旅游危机预警与响应机制，同时根据后疫情时代旅游市场的消费模式与偏好，调整旅游产品供给结构，增加安全、健康旅游产品和服务供给，缓解新冠疫情对旅游经济的冲击，推进旅游经济恢复重振；对下游流域而言，一方面要因地制宜地用好沿海的区位优势，推动旅游投融资项目，为黄河旅游经济发展奠定坚实的物质基础，另一方面要合理规避“虹吸效应”的消极影响，重组黄河流域旅游经济系统合作的内部结构、创新旅游经济系统合作的实施功能，强化旅游经济系统组团优势和集聚效应，充分发挥下游地区旅游经济韧性的正向溢出作用。由于旅游经济系统韧性理论的丰富性和外延复杂性，本文尚存在如下不足之处：受统计数据资源限制，指标体系构建有待进一步完善，未来可通过实地调研获得一手数据或遥感数据以完善指标体系；此外，旅游经济系统韧性是一个动态演化的概念，受内外部主客观因素驱动，黄河流域与周边区域旅游经济系统韧性时空关联的路径和驱动机制等议题也将是后续研究关注的重点方向。参考文献：杜文蠧，施益军，徐丽华，等风险扰动下的城市经济韧性多维测度与分析：以长三角地区为例地理科学进展，（）：，（）：，：，：，：，（）：孙业红，宋雨新从韧性看可持续旅游发展旅游学刊，（）：，：，：珦，（）：，：，（）：杨友宝，王荣成边境民族地区旅游经济系统脆弱性评价及其优化路径：以延边朝鲜族自治州为例地域研究与开发，（）：白丽明旅游经济系统弹性的研究综述中国市场，（）：冯凌，郭嘉欣，刘乙疫情后期旅游目的地恢复力及其时空动态分析的研究方法（英文），（）：张铁生，孙根年，马丽君危机事件对张家界旅游影响评价：基于本底线模型的高分辨率分析经济地理，（）：李锋基于集对分析法（）的中国旅游经济系统脆弱性测度研究旅游科学，（）：徐凯舟，吴宇冉，严晓琳，等疫情冲击下旅游经济韧性评价指数分析当代旅游，（）：，：，：，：，：，：，：，（）：，：，（）：王倩，赵林，于伟，等中国旅游经济系统韧性的时空变化特征与影响因素分析地理与地理信息科学，（）：崔文静，薛涛粤港澳大湾区旅游产业经济韧性影响因素研究特区经济，（）：庞冬彦，赵林，于伟，等山东省旅游经济系统的韧性测度与障页第地 理 与 地 理 信 息 科 学第 卷碍因素分析资源与产业，（）：田磊，孙凤芝，张宗斌，等黄河流域旅游业与资源环境协调发展实证研究干旱区资源与环境，（）：中华人民共和国文化和旅游部中国旅游统计年鉴 北京：中国旅游出版社，张梦梦，张丽君，秦耀辰，等黄河流域乡镇人口：经济增长与收缩格局及影响因素地理科学进展，（）：，（）：马慧强，廉倩文，论宇超，等基于 神经网络的旅游经济系统脆弱性省际空间分异资源科学，（）：王兆峰，陈青青 年以来长江经济带旅游生态安全时空格局演化及趋势预测生态学报，（）：孙根年，冯茂娥西部入境旅游市场竞争态与资源区位的关系西北大学学报（自然科学版），（）：曹乃刚，赵林，贾建琦基于 框架的中国海洋经济发展质量测度与影响因素分析地理与地理信息科学，（）：许瑶，纪建悦，许玉洁中国养殖海域利用效率空间非均衡格局及成因资源科学，（）：靳诚，陆玉麒基于空间变差函数的长江三角洲经济发展差异演变研究地理科学，（）：王兆峰，霍菲菲湖南武陵山区旅游产业集聚与区域经济发展关系测度地域研究与开发，（）：朱媛媛，汪紫薇，乔花芳，等大别山革命老区旅游地“乡土生态”系统韧性演化规律及影响机制自然资源学报，（）：叶琴，曾刚，戴劭?，等不同环境规制工具对中国节能减排技术创新的影响：基于 个地级市面板数据中国人口资源与环境，（）：张辉，耿成轩，谢申祥中国旅游产业发展的人口因素分析数学的实践与认识，（）：谢涓，李玉双，韩峰环境规制与经济增长：基于中国省际面板联立方程的分析经济经纬，（）：犛 狆 犪 狋 犻 狅 狋 犲 犿 狆 狅 狉 犪 犾犇 犻 犳 犳 犲 狉 犲 狀 犮 犲 狊 犪 狀 犱犐 狀 犳 犾 狌 犲 狀 犮 犻 狀 犵犉 犪 犮 狋 狅 狉 狊 狅 犳犚 犲 狊 犻 犾 犻 犲 狀 犮 犲狅 犳犜 狅 狌 狉 犻 狊 犿犈 犮 狅 狀 狅 犿 犻 犮犛 狔 狊 狋 犲 犿犻 狀 狋 犺 犲犢 犲 犾 犾 狅 狑犚 犻 狏 犲 狉犅 犪 狊 犻 狀 ，（犆 狅 犾 犾 犲 犵 犲 狅 犳犜 狅 狌 狉 犻 狊 犿，犎 狌 狀 犪 狀犖 狅 狉 犿 犪 犾犝 狀 犻 狏 犲 狉 狊 犻 狋 狔，犆 犺 犪 狀 犵 狊 犺 犪 ，犆 犺 犻 狀 犪）犃 犫 狊 狋 狉 犪 犮 狋：，犓 犲 狔狑 狅 狉 犱 狊：；页第第期 王兆峰，张先甜：黄河流域旅游经济系统韧性的时空差异特征及其影响因素