В данном исследовании использовались общедоступные геопространственные данные, данные дистанционного зондирования, данные планирования и агрегированные социально-экономические наборы данных. В работе не участвовали люди, животные, не использовались клинические материалы или идентифицируемая персональная информация. В связи с этим получение одобрения комитета по этике не требовалось. Конкретные программные среды, вычислительные пакеты (включая библиотеки пространственной эконометрики и визуализации высокого разрешения), цифровые идентификаторы объектов (DOI) наборов данных и исходные URL-адреса всех ресурсов, необходимых для воспроизведения данного рабочего процесса, подробно описаны в исчерпывающей Таблице материалов.
1. Район исследования и аналитическая база
Анализ проводился в Синьцзян-Уйгурском автономном районе на северо-западе Китая, при этом особое внимание уделялось определенным границам городской застройки ключевых городов (Рисунок 1) для точного фиксирования динамики локального расширения. Данная граница последовательно использовалась для обрезки растров, моделирования землепользования, оценки запасов углерода и агрегации административных единиц. Аналитическая выборка (n) включала 105 административных единиц уровня уезда, проанализированных в трех дискретных временных точках (20, 2010 и 2020 годы), при этом моделирование землепользования включало 50 стохастических повторностей для учета пространственной неопределенности.
Был применен многомасштабный аналитический подход, объединяющий анализ на уровне сетки и на уровне административных единиц. На уровне сетки данные о землепользовании и почвенно-растительном покрове, интенсивность ночного освещения, степень покрытия водонепроницаемыми поверхностями, вегетационные индексы, топографические переменные и ограничения экологического планирования использовались для выявления закономерностей урбанизации земель, моделирования будущих сценариев землепользования и оценки распределения запасов углерода. На уровне административных единиц социально-экономические, транспортные и экологические переменные были включены в пространственный эконометрический анализ для оценки прямых и косвенных эффектов урбанизации земель и экологического планирования на динамику запасов углерода.
Аналитическая база состояла из пяти компонентов. Во-первых, многовременные данные дистанционного зондирования и землепользования использовались для характеристики пространственно-временных закономерностей урбанизации земель, включая расширение застроенных территорий, интенсивность застройки и реструктуризацию землепользования. Во-вторых, были разработаны сценарии охраны экологии, защиты сельскохозяйственных угодий и низкоуглеродного планирования для моделирования будущих структур землепользования. В-третьих, модуль хранения углерода InVEST использовался для количественной оценки распределения и изменения запасов углерода при исторических и будущих условиях землепользования. В-четвертых, был проведен анализ глобальной и локальной пространственной автокорреляции для выявления паттернов кластеризации и регионов «горячих точек», связанных с вариациями запасов углерода. В-пятых, пространственные эконометрические модели использовались для оценки величины, направления и эффектов перелива урбанизации земель, интенсивности экологического планирования, а также экологических и социально-экономических факторов влияния на динамику запасов углерода.
Рисунок 1 представляет общую аналитическую схему рабочего процесса. Этот процесс объединяет определение области исследования, интеграцию данных из нескольких источников, идентификацию урбанизации земель, разработку сценариев экологического планирования, симуляцию землепользования, оценку запасов углерода, анализ пространственной автокорреляции и пространственное эконометрическое моделирование. Пространственные данные из различных источников служили входным слоем. Идентификация урбанизации земель и анализ исторического землепользования образовали слой характеристики структуры. Сценарии экологического планирования и оценка запасов углерода составили слой прогнозирования и оценки воздействия. Пространственная автокорреляция и эконометрический анализ использовались для выявления пространственных механизмов и обоснования интерпретации политических решений.

Рисунок 1Район исследования, интеграция данных и технический рабочий процесс. (A) Многомасштабный пространственный контекст и точные границы городского анализа для ключевых городов Синьцзян-Уйгурского автономного района, северо-запад Китая. (B) Интегрированный технический рабочий процесс, объединяющий пространственные данные из нескольких источников, идентификацию урбанизации земель на основе дистанционного зондирования, проектирование сценариев экологического планирования, моделирование землепользования PLUS, оценку запасов углерода InVEST, пространственную автокорреляцию, анализ «горячих точек» и пространственное эконометрическое моделирование. ND = естественное развитие; EP = экологическая защита; CP = защита сельскохозяйственных угодий; LC = низкоуглеродная оптимизация; OLS = обычный метод наименьших квадратов; SAR = пространственная авторегрессионная модель; SEM = пространственная модель ошибки; SDM = пространственная модель Дурбина. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
2. Источники данных и система переменных
Для обеспечения идентификации урбанизации земель, моделирования сценариев экологического планирования, оценки запасов углерода и пространственного эконометрического анализа был собран многоисточниковый набор данных. Набор данных включал данные о землепользовании и земельном покрове, долю водонепроницаемых поверхностей, интенсивность ночного освещения, нормализованный разностный вегетационный индекс (NDVI), топографические переменные, показатели транспортной доступности, гидрологические переменные, плотность населения, плотность валового внутреннего продукта (ВВП), планировочные ограничения, климатические переменные и параметры углеродного пула. В Таблице 1 обобщены категории наборов данных, описания переменных, единицы измерения, пространственное и временное разрешение, источники и варианты аналитического применения. Основные наборы данных были получены для 200, 2010 и 2020 годов. Данные о землепользовании и основные продукты дистанционного зондирования были приведены к разрешению от 30 m до 1 km в зависимости от аналитических требований.
Таблица 1: Источники данных и система переменных. В таблице представлены категория набора данных, переменная, описание, единица измерения, пространственное/временное разрешение, источник и аналитическое применение для идентификации урбанизации земель, моделирования сценариев, учета запасов углерода InVEST и пространственного эконометрического моделирования. NDVI = нормализованный разностный вегетационный индекс; DEM = цифровая модель рельефа; LULC = землепользование/земельный покров. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Все пространственные слои были спроецированы в единую систему координат, обрезаны по границам исследуемой области, показанной на Рисунке 1, и подвергнуты ресемплингу или агрегации до требуемого масштаба сетки или административных единиц. Перед проведением анализа пространственного наложения была выполнена проверка на наличие пропущенных значений. Растровые наборы данных, использованные в анализах PLUS и InVEST, были выровнены на уровне пикселей, в то время как социально-экономические переменные были агрегированы по административным единицам для пространственной регрессии. Данный рабочий процесс предварительной обработки обеспечил единообразие пространственных единиц и сопоставимость временных точек при проведении временных анализов, имитации сценариев, расчета запасов углерода и эконометрического моделирования.
Данные о землепользовании и почвенном покрове служили основным входным материалом для выявления исторических изменений и моделирования землепользования в будущем. Исходные категории землепользования были переклассифицированы на шесть классов: пашни, лесные земли, grassland, водные объекты, застроенные земли и неиспользуемые земли. Данная схема классификации последовательно использовалась для анализа изменений землепользования, разработки сценариев и расчета запасов углерода с помощью InVEST. Доля застроенных земель, покрытие водонепроницаемыми поверхностями и интенсивность ночного освещения использовались в качестве индикаторов урбанизации земель. Покрытие водонепроницаемыми поверхностями и интенсивность ночного освещения служили прокси-показателями дистанционного зондирования для оценки интенсивности застройки и человеческой деятельности. NDVI был включен в качестве экологического индикатора в последующий пространственный эконометрический анализ.
Топографические переменные и переменные местоположения включали высоту над уровнем моря, уклон, расстояние до основных дорог, плотность дорожной сети и близость к крупным рекам. Эти переменные использовались в качестве факторов в моделировании землепользования и в качестве контрольных переменных в пространственных эконометрических моделях. В совокупности они представляли собой ограничения рельефа, транспортную доступность и гидрологическую связность.
Социально-экономические переменные включали плотность населения и плотность ВВП, дополненные плотностью дорожной сети, где эта информация была доступна. Данные переменные использовались для характеристики интенсивности освоения и давления антропогенной деятельности, связанных с изменениями запасов углерода. Структура переменных объединяла показатели освоения земель, социально-экономического давления и экологического контекста для представления совокупности факторов, влияющих на динамику запасов углерода.
Индекс интенсивности экологического планирования был разработан с использованием экологических «красных линий», природных заповедников, буферных зон водных объектов, районов с ограничениями по уклону и других экологически чувствительных регионов. Этот индекс служил как ограничивающим слоем в симуляциях землепользования в будущем, так и объясняющей переменной в пространственном эконометрическом анализе.
Для оценки запасов углерода с помощью InVEST для каждого класса землепользования были определены четыре параметра пулов углерода: углерод надземной биомассы, углерод подземной биомассы, органический углерод почвы и углерод мертвого органического вещества. Значения параметров были получены из опубликованных региональных исследований, руководств по InVEST, а также на основе характеристик местного растительного покрова и соотнесены с единой системой классификации землепользования. Под запасом углерода понимается общее расчетное количество углерода, содержащегося в этих четырех пулах, в то время как плотность запасов углерода означает запас углерода на единицу площади.
Переменные были разделены на три категории. Первая категория включала показатели урбанизации земель, в том числе долю застроенных земель, покрытие водонепроницаемыми поверхностями и интенсивность ночного освещения. Вторая категория включала экологические и планировочные переменные, такие как NDVI, интенсивность экологического планирования, высота, уклон и близость к гидрологическим объектам. Третья категория включала социально-экономические факторы, в том числе плотность населения, плотность ВВП и плотность дорожной сети. Эти переменные использовались для оценки взаимосвязей между урбанизацией земель, экологическим планированием и динамикой запасов углерода.
3. Дистанционное зондирование для идентификации урбанизации земель
Пространственно-временные закономерности урбанизации земель были определены с помощью трех взаимодополняющих показателей: расширения застроенных земель, покрытия водонепроницаемыми поверхностями и интенсивности ночного освещения. Расширение застроенных земель было извлечено из последовательных карт землепользования для определения физических границ и прогрессии городского развития от центральных застроенных районов к прилегающим регионам. Покрытие водонепроницаемыми поверхностями рассчитывалось для количественной оценки интенсивности застройки на уровне сетки, в то время как растровые данные об интенсивности ночного освещения были обработаны и нормализованы для отображения человеческой деятельности и функциональной концентрации. Эти пространственные слои были впоследствии интегрированы для создания исторического набора данных по урбанизации земель. Комбинированные показатели использовались для выявления зон устойчивого усиления развития, переходных зон и относительно стабильных регионов, а также для оценки пространственного соответствия между физическим освоением земель и функциональной урбанизацией.
Для количественной оценки масштабов и направления преобразования земель между категориями пахотных земель, лесных массивов, grassland, водных объектов, застроенных территорий и неиспользуемых земель была построена матрица перераспределения землепользования. Особое внимание было уделено переходам из категорий экологических и сельскохозяйственных земель в категорию застроенных территорий. Матрица использовалась для выявления доминирующих путей преобразования и основных типов земель, способствующих расширению городов.
Был проведен анализ паттернов ландшафта для оценки структурных изменений, связанных с урбанизацией земель. Метрики включали плотность пятен, плотность границ, индекс формы ландшафта и показатели фрагментации. Эти метрики рассчитывались для каждого периода исследования и использовались для количественной оценки изменений в конфигурации ландшафта, пространственной непрерывности и фрагментации, связанных с расширением застроенных территорий.
4. Разработка сценариев экологического планирования и моделирование землепользования
Моделирование будущих структур землепользования проводилось с помощью модели Patch-generating Land Use Simulation (PLUS) в рамках четырех сценариев планирования: естественное развитие (ND), экологическая защита (EP), защита пахотных земель (CP) и низкоуглеродная оптимизация (LC). Допущения для сценариев, правила преобразования земель, ограниченные типы земель и ожидаемые результаты по запасам углерода обобщены в Таблице 2.
Таблица 2: Правила управления сценариями и ограничения переходов. В таблице определены вероятности преобразования земель и пространственные ограничения для сценариев естественного развития (ND), экологической защиты (EP), защиты пахотных земель (CP) и низкоуглеродной оптимизации (LC). Примечания: Сокращения сценариев строго соответствуют тем, что используются в разделах моделирования PLUS и «Результаты». Настройки сценариев определяют ориентированные на политику правила преобразования и ограниченные типы земель, используемые в симуляциях; все установленные законом зоны отчуждения остались неизменяемыми на итоговых пространственных слоях. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Исторические карты землепользования и пространственные переменные-драйверы были интегрированы в модель PLUS для оценки базовых вероятностей перехода земель с использованием модуля клеточных автоматов (CA). Для сценариев экологической защиты (EP), защиты пахотных земель (CP) и низкоуглеродной оптимизации (LC) экологические «красные линии», основные permanently basic cropland, буферные зоны водных объектов и другие планировочные ограничения были включены в качестве слоев пространственных ограничений для лимитирования конверсии земель в соответствии с предопределенными правилами сценариев. Впоследствии для каждого сценария были определены целевые количественные показатели спроса на землю на 2030 год, и модуль CA был запущен с использованием 50 стохастических повторностей для создания окончательных прогнозов землепользования.
Моделирование землепользования проводилось с использованием модели PLUS и ее структуры клеточных автоматов (CA) с многотипной стохастической генерацией участков. Для оценки вероятности расширения земель для каждого класса землепользования использовались исторические карты землепользования, а также переменные факторов окружающей среды и социально-экономические показатели. Затем были определены объемы потребности в земле в будущем в соответствии с требованиями каждого сценария планирования, и для целевого периода моделирования были созданы соответствующие карты землепользования.
Эффективность модели оценивали с помощью исторического обратного прогнозирования (back-casting) перед проведением симуляции на будущие периоды. Более ранние карты землепользования и соответствующие переменные-факторы использовались для симуляции более поздней наблюдаемой карты землепользования. Соответствие между симулированным и наблюдаемым распределением землепользования оценивали с помощью общей точности (OA), коэффициента Каппа и показателя качества (FoM). Валидация проводилась как на общем уровне, так и по основным классам землепользования. Историческое обратное прогнозирование с 2010 по 2020 год дало общую точность (OA) 93,4%, коэффициент Каппа 0,89 и показатель качества (FoM) 0,26, что указывает на высокую надежность способности к пространственному прогнозированию в последующих многосценарных симуляциях.
Для оценки устойчивости результатов моделирования были проведены анализ неопределенности и анализ чувствительности. Анализ чувствительности в модели PLUS позволил изучить влияние альтернативных настроек сопротивления переходу и параметров весов соседства для основных классов землепользования. Анализ чувствительности в InVEST был направлен на оценку влияния вариаций коэффициентов пулов углерода для различных типов землепользования. В рамках пространственного эконометрического анализа чувствительности сравнивались альтернативные спецификации матрицы пространственных весов. Эти анализы использовались для проверки того, сохраняются ли ранжирование сценариев и направление основных эффектов урбанизации и экологического планирования при альтернативных настройках параметров. В частности, была строго подтверждена стабильность результатов: ранжирование сценариев и отрицательные пространственные эффекты перелива от урбанизации земель оставались неизменными при изменении параметров сопротивления переходу и коэффициентов пулов углерода в пределах ±15%.
5. Оценка запасов углерода
Оценка запасов углерода проводилась с использованием модели углерода InVEST. Базовые параметры для четырех пулов углерода (надземная биомасса, подземная биомасса, органический углерод почвы и мертвая органическая материя) были присвоены каждому переклассифицированному типу землепользования с использованием биофизических значений, приведенных в Таблице 3. Затем в модель были импортированы исторические карты землепользования (20–2020 гг.) и растровые наборы данных о будущем землепользовании, смоделированные с помощью PLUS, и интегрированы с соответствующими параметрами плотности углерода. Впоследствии была запущена модель для оценки общего регионального запаса углерода (Tg C), плотности запаса углерода на уровне сетки (Mg C/ha) и пространственных карт изменения запасов углерода (ΔC) для исторических и будущих сценариев.
Переклассифицированные карты землепользования были сопоставлены с соответствующими базовыми значениями параметров пула углерода, представленными в Таблице 3. Чувствительность оценок запасов углерода к неопределенности параметров была протестирована путем изменения этих базовых значений на ±15% с подробным анализом результатов чувствительности. Модель использовалась для расчета общего запаса углерода, плотности запаса углерода и изменения запаса углерода для каждого класса землепользования в исторических условиях и при сценариях будущего планирования.
Были оценены три категории выходных данных. Во-первых, были рассчитаны общий региональный запас углерода и временные тренды для количественной оценки величины и направления изменения запасов углерода с течением времени. Во-вторых, было построено картирование пространственного распределения запасов углерода и их изменения для выявления зон удержания и потери углерода. В-третьих, оценки запасов углерода сравнивались между различными сценариями планирования для анализа относительного влияния стратегий экологической защиты, защиты сельскохозяйственных угодий и низкоуглеродной оптимизации на сохранение запасов углерода.
Таблица 3: Базовые параметры пулов углерода для различных типов землепользования. В таблице представлены показатели углерода надземной биомассы, углерода подземной биомассы, органического углерода почвы, углерода мертвой органики и общей плотности углерода, использованные в модуле запасов углерода InVEST. Единицы измерения — Mg C/ha. Примечания: Значения представляют собой базовые параметры, используемые в модели InVEST. Общая плотность углерода равна сумме четырех пулов углерода. Результаты анализа чувствительности при изменении этих базовых значений на ±15% представлены в Таблице 6. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
6. Пространственная автокорреляция и пространственный эконометрический анализ
Для определения того, проявляют ли запасы углерода и их изменение значимую пространственную зависимость, был проведен анализ пространственной автокорреляции. Для оценки общей степени пространственной кластеризации распределения запасов углерода и их изменения на всей исследуемой области рассчитывался глобальный индекс Морана (Global Moran's I). Затем использовался локальный индекс Морана (Local Moran's I) для выявления локальных закономерностей пространственной ассоциации, включая кластеры «высокий-высокий», «низкий-низкий», «высокий-низкий» и «низкий-высокий». Анализ «горячих точек» проводился для определения областей концентрированной потери углерода и сохранения запасов углерода.
Для изучения взаимосвязей между урбанизацией земель, экологическим планированием, состоянием окружающей среды, социально-экономическими факторами и динамикой запасов углерода использовались модели пространственной эконометрики. В качестве зависимой переменной выступали плотность запасов углерода или изменение запасов углерода. Объясняющие переменные включали индекс урбанизации земель, долю застроенных земель, интенсивность экологического планирования, нормализованный разностный вегетационный индекс (NDVI), плотность населения, плотность ВВП, плотность дорожной сети, высоту над уровнем моря, уклон, среднегодовое количество осадков и среднегодовую температуру.
В качестве базовой модели использовалась регрессия с обычным наименьшим квадратом (OLS). Перед оценкой моделей пространственной авторегрессии (SAR), пространственной ошибки (SEM) и пространственной модели Дурбина (SDM) оценивалась остаточная пространственная зависимость. Эффективность моделей и оценки коэффициентов сравнивались для различных спецификаций.
Для представления соседских отношений между административными единицами была построена матрица пространственных весов со стандартизацией по строкам. Основная спецификация была основана на пространственной смежности; в рамках анализа устойчивости сравнивались альтернативные матрицы пространственных весов на основе расстояния и ближайших соседей, где это было применимо.
Для оценки локальных взаимосвязей и пространственного перелива были рассчитаны прямые, косвенные и общие эффекты на основе пространственной модели зависимости (SDM). Переменные рельефа, доступности, растительности, климата и социально-экономические показатели были включены в качестве контрольных переменных для снижения систематической ошибки, вызванной пропуском значимых переменных. Коэффициенты модели интерпретировались как условные пространственные ассоциации, а не как окончательные причинно-следственные связи. Переменные с p-значениями, превышающими общепринятые пороги значимости, интерпретировались как слабые или косвенные доказательства и не рассматривались как статистически достоверные эффекты.