Foram utilizados conjuntos de dados geoespaciais, de sensoriamento remoto, de planejamento e socioeconômicos agregados, disponíveis publicamente, neste estudo. Não foram envolvidos participantes humanos, sujeitos animais, materiais clínicos ou informações pessoais identificáveis. A aprovação ética institucional, portanto, não foi necessária. Os ambientes de software específicos, pacotes computacionais (incluindo bibliotecas de econometria espacial e de visualização de alta resolução), identificadores digitais (DOIs) dos conjuntos de dados e URLs de origem de todos os recursos necessários para reproduzir este fluxo de trabalho estão detalhados na abrangente Tabela de Materiais.
1. Área de estudo e estrutura analítica
A análise foi conduzida na Região Autônoma de Xinjiang Uigur, no noroeste da China, com foco analítico específico nas extensões urbanas delimitadas de cidades-chave (Figura 1) para capturar com precisão a dinâmica de expansão localizada. Esse limite foi utilizado de forma consistente para o recorte de matrizes raster, simulação de uso da terra, avaliação de estoques de carbono e agregação de unidades administrativas. O tamanho da amostra analítica (n) compreendeu 105 unidades administrativas de nível municipal analisadas em três nós temporais distintos (2000, 2010 e 2020), com simulações de uso da terra incorporando 50 réplicas estocásticas para capturar a incerteza espacial.
Utilizou-se uma estrutura analítica multiescalar que integra análises em nível de grade e em unidades administrativas. Na escala da grade, utilizaram-se dados de uso e cobertura da terra, intensidade da luz noturna, cobertura de superfícies impermeáveis, índices de vegetação, variáveis topográficas e restrições de planejamento ecológico para identificar padrões de urbanização da terra, simular cenários futuros de uso da terra e avaliar a distribuição de estoques de carbono. Na escala da unidade administrativa, variáveis socioeconômicas, de transporte e ambientais foram incorporadas a análises econométricas espaciais para avaliar os efeitos diretos e de derramamento da urbanização da terra e do planejamento ecológico sobre a dinâmica dos estoques de carbono.
A estrutura analítica compreendeu cinco componentes. Primeiro, dados multi-temporais de sensoriamento remoto e de uso da terra foram utilizados para caracterizar os padrões espaço-temporais da urbanização da terra, incluindo a expansão de áreas edificadas, a intensidade de desenvolvimento e a reestruturação do uso da terra. Segundo, foram desenvolcidos cenários de proteção ecológica, proteção de terras agrícolas e planejamento de baixo carbono para simular padrões futuros de uso da terra. Terceiro, o módulo de armazenamento de carbono do InVEST foi utilizado para quantificar a distribuição e as alterações nos estoques de carbono sob condições históricas e futuras de uso da terra. Quarto, foram realizadas análises de autocorrelação espacial global e local para identificar padrões de agrupamento e regiões de foco associadas à variação dos estoques de carbono. Quinto, modelos econométricos espaciais foram utilizados para avaliar a magnitude, a direção e os efeitos de derramamento da urbanização da terra, da intensidade do planejamento ecológico e dos fatores ambientais e socioeconômicos sobre a dinâmica dos estoques de carbono.
Figura 1 apresenta o fluxo analítico geral. O fluxo integra a definição da área de estudo, a integração de dados multissource, a identificação da urbanização da terra, o desenvolvimento de cenários de planejamento ecológico, a simulação de uso da terra, a avaliação de estoques de carbono, a análise de autocorrelação espacial e a modelagem econométrica espacial. Dados espaciais multissource serviram como camada de entrada. A identificação da urbanização da terra e a análise histórica do uso da terra formaram a camada de caracterização de padrões. Os cenários de planejamento ecológico e a avaliação dos estoques de carbono formaram a camada de previsão e avaliação de impactos. As análises de autocorrelação espacial e econométrica foram utilizadas para identificar mecanismos espaciais e apoiar a interpretação de políticas.

Figura 1: Área de estudo, integração de dados e fluxo técnico de trabalho. (A) Contexto espacial multiescalar e limites analíticos urbanos precisos para cidades-chave da Região Autônoma Uigur de Xinjiang, no noroeste da China. (B) Fluxo técnico integrado que conecta dados espaciais multissensoriais, identificação da urbanização do solo baseada em sensoriamento remoto, projeto de cenários de planejamento ecológico, simulação de uso da terra com o modelo PLUS, avaliação de estoque de carbono com o InVEST, autocorrelação espacial, análise de hotspots e modelagem econométrica espacial. ND = desenvolvimento natural; EP = proteção ecológica; CP = proteção de áreas agrícolas; LC = otimização de baixo carbono; OLS = mínimos quadrados ordinários; SAR = modelo autorregressivo espacial; SEM = modelo de erro espacial; SDM = modelo de Durbin espacial. Clique aqui para visualizar uma versão maior desta figura.
2. Fontes de dados e sistema de variáveis
Um conjunto de dados multissensorial foi compilado para apoiar a identificação da urbanização do solo, a simulação de cenários de planejamento ecológico, a avaliação de estoques de carbono e a análise espacial econométrica. O conjunto de dados incluiu informações sobre uso e cobertura da terra, fração de superfícies impermeáveis, intensidade da luz noturna, índice de vegetação por diferença normalizada (NDVI), variáveis topográficas, indicadores de acessibilidade ao transporte, variáveis hidrológicas, densidade populacional, densidade do produto interno bruto (PIB), restrições de planejamento, variáveis climáticas e parâmetros dos reservatórios de carbono. Tabela 1 resume a categoria do conjunto de dados, a descrição das variáveis, unidade, resolução espacial e temporal, fonte e aplicação analítica. Os principais conjuntos de dados foram obtidos para os anos 2000, 2010 e 2020. Os dados de uso da terra e os principais produtos de sensoriamento remoto foram harmonizados para resoluções que variam de 30 m a 1 km, conforme as exigências analíticas.
Tabela 1: Fontes de dados e sistema de variáveis. A tabela apresenta a categoria do conjunto de dados, variável, descrição, unidade, resolução espacial/temporal, fonte e uso analítico para a identificação da urbanização da terra, simulação de cenários, contabilização de estoque de carbono no InVEST e modelagem econométrica espacial. NDVI = índice de vegetação por diferença normalizada; DEM = modelo digital de elevação; LULC = uso da terra/cobertura do solo. Clique aqui para baixar esta tabela.
Todas as camadas espaciais foram projetadas em um sistema de coordenadas comum, recortadas até o limite da área de estudo mostrado na Figura 1 e reamostradas ou agregadas à escala de grade ou unidade administrativa requerida. Valores ausentes foram verificados antes da análise de sobreposição espacial. Os conjuntos de dados matriciais utilizados nas análises do PLUS e do InVEST foram alinhados ao nível do pixel, enquanto as variáveis socioeconômicas foram agregadas a unidades administrativas para regressão espacial. Este fluxo de trabalho de pré-processamento garantiu unidades espaciais consistentes e pontos temporais comparáveis nas análises temporais, simulações de cenários, cálculos de estoque de carbono e modelagem econométrica.
Os dados de uso/ocupação da terra serviram como entrada principal para a detecção de mudanças históricas e simulação futura do uso da terra. As categorias originais de uso da terra foram reclassificadas em seis classes: terra arável, terra florestal, pastagem, corpos d'água, área urbanizada e terra não utilizada. Este esquema de classificação foi utilizado de forma consistente na análise de mudanças do uso da terra, no desenvolvimento de cenários e nos cálculos de estoque de carbono no InVEST. A proporção de área urbanizada, a cobertura de superfícies impermeáveis e a intensidade da luz noturna foram utilizadas como indicadores de urbanização da terra. A cobertura de superfícies impermeáveis e a intensidade da luz noturna funcionaram como proxies de sensoriamento remoto para intensidade de desenvolvimento e atividade humana. O NDVI foi incluído como um indicador ecológico nas análises econométricas espaciais subsequentes.
As variáveis topográficas e de localização incluíram elevação, declividade, distância às principais estradas, densidade de estradas e proximidade com rios principais. Essas variáveis foram utilizadas como fatores impulsionadores na simulação de uso da terra e como variáveis de controle em modelos econométricos espaciais. Em conjunto, representaram restrições do terreno, acessibilidade ao transporte e conectividade hidrológica.
As variáveis socioeconômicas incluíram densidade populacional e densidade de PIB, complementadas pela densidade da rede viária quando disponível. Essas variáveis foram utilizadas para caracterizar a intensidade do desenvolvimento e as pressões da atividade humana associadas às mudanças nos estoques de carbono. A estrutura de variáveis integrou indicadores de desenvolvimento territorial, pressão socioeconômica e contexto ecológico para representar os múltiplos fatores que influenciam a dinâmica dos estoques de carbono.
Um índice de intensidade de planejamento ecológico foi desenvolvido utilizando linhas vermelhas ecológicas, reservas naturais, zonas tampão de corpos d'água, áreas com restrição de declividade e outras regiões ambientalmente sensíveis. Este índice serviu tanto como uma camada de restrição em simulações futuras de uso da terra quanto como uma variável explicativa em análises econométricas espaciais.
Para a avaliação de estoque de carbono no InVEST, foram compilados quatro parâmetros de compartimentos de carbono para cada classe de uso da terra: carbono da biomassa acima do solo, carbono da biomassa abaixo do solo, carbono orgânico do solo e carbono da matéria orgânica morta. Os valores dos parâmetros foram obtidos a partir de estudos regionais publicados, documentos orientadores do InVEST e características locais de cobertura do solo, e foram mapeados para o sistema unificado de classificação de uso da terra. Estoque de carbono refere-se ao carbono total estimado contido nesses quatro compartimentos, enquanto densidade de estoque de carbono refere-se ao estoque de carbono por unidade de área.
As variáveis foram agrupadas em três categorias. A primeira categoria compreendia indicadores de urbanização da terra, incluindo a proporção de área construída, cobertura de superfícies impermeáveis e intensidade da luz noturna. A segunda categoria incluía variáveis ecológicas e de planejamento, tais como NDVI, intensidade de planejamento ecológico, elevação, declividade e proximidade hidrológica. A terceira categoria compreendia fatores socioeconômicos, incluindo densidade populacional, densidade de PIB e densidade de vias. Essas variáveis foram utilizadas para avaliar as relações entre urbanização da terra, planejamento ecológico e dinâmicas de estoque de carbono.
3. Identificação baseada em sensoriamento remoto da urbanização do solo
Os padrões espaço-temporais da urbanização da terra foram identificados utilizando três indicadores complementares: expansão da área construída, cobertura de superfícies impermeáveis e intensidade da luz noturna. A expansão da área construída foi extraída de mapas sequenciais de uso da terra para delimitar a extensão física e a progressão do desenvolvimento urbano a partir das áreas centrais construídas até as regiões circunvizinhas. A cobertura de superfícies impermeáveis foi calculada para quantificar a intensidade do desenvolvimento em nível de grade, enquanto os dados matriciais de intensidade da luz noturna foram processados e normalizados para representar a atividade humana e a concentração funcional. Essas camadas espaciais foram posteriormente integradas para construir o conjunto de dados histórico de urbanização da terra. Os indicadores combinados foram utilizados para identificar áreas de intensificação contínua do desenvolvimento, zonas de transição e regiões relativamente estáveis, bem como para avaliar a correspondência espacial entre o desenvolvimento físico da terra e a urbanização funcional.
Foi construída uma matriz de transição do uso da terra para quantificar a magnitude e a direção da conversão do uso da terra entre terras aráveis, florestas, pastagens, corpos d'água, áreas de construção e terras não utilizadas. Uma atenção especial foi dada às transições das classes de terra ecológica e agrícola para áreas de construção. A matriz foi utilizada para identificar os principais caminhos de conversão e as classes de terra de origem que mais contribuem para a expansão urbana.
A análise de padrões da paisagem foi realizada para avaliar as alterações estruturais associadas à urbanização da terra. As métricas incluíram densidade de fragmentos, densidade de bordas, índice de forma da paisagem e indicadores relacionados à fragmentação. Essas métricas foram calculadas para cada período de estudo e utilizadas para quantificar as mudanças na configuração da paisagem, na continuidade espacial e na fragmentação associadas à expansão das áreas edificadas.
4. Projeto de cenários de planejamento ecológico e simulação de uso da terra
Padrões futuros de uso da terra foram simulados utilizando o modelo de Simulação de Uso da Terra Gerador de Fragmentos (PLUS) em quatro cenários de planejamento: desenvolvimento natural (ND), proteção ecológica (EP), proteção de terras agrícolas (CP) e otimização de baixo carbono (LC). As premissas dos cenários, regras de conversão de uso da terra, tipos de terra restritos e os resultados esperados de estoque de carbono são resumidos na Tabela 2.
Tabela 2: Regras de controle de cenário e restrições de transição. A tabela define as probabilidades de conversão de uso da terra e as restrições espaciais para os cenários de desenvolvimento natural (ND), proteção ecológica (EP), proteção de áreas agrícolas (CP) e otimização de baixo carbono (LC). Observações: as abreviações dos cenários correspondem rigorosamente às utilizadas na modelagem PLUS e nas seções de Resultados. As configurações dos cenários definem regras de conversão orientadas por políticas e tipos de terra restritos utilizados nas simulações; todas as zonas de exclusão estatutárias permaneceram não conversíveis nas camadas espaciais finais. Clique aqui para baixar esta tabela.
Mapas históricos de uso da terra e variáveis espaciais de condução foram integrados ao modelo PLUS para estimar probabilidades de transição de terra de referência utilizando o módulo de autômatos celulares (CA). Para os cenários de proteção ecológica (EP), proteção de terras agrícolas (CP) e otimização de baixo carbono (LC), linhas vermelhas ecológicas, terras agrícolas básicas permanentes, zonas tampão ao redor de corpos d'água e outras restrições de planejamento foram incorporadas como camadas espaciais de restrição para limitar a conversão de terras de acordo com as regras predefinidas dos cenários. Quantidades de demanda de terra-alvo para 2030 foram posteriormente definidas para cada cenário, e o módulo CA foi executado utilizando 50 réplicas estocásticas para gerar as projeções finais de uso da terra.
A simulação do uso da terra foi realizada utilizando o modelo PLUS e sua estrutura de autômatos celulares (CA) com geração estocástica de fragmentos de múltiplos tipos. Mapas históricos de uso da terra e variáveis ambientais e socioeconômicas condutoras foram utilizados para estimar as probabilidades de expansão da terra para cada classe de uso da terra. As quantidades futuras de demanda de terra foram então especificadas de acordo com os requisitos de cada cenário de planejamento, e os mapas correspondentes de uso da terra foram gerados para o período-alvo da simulação.
O desempenho do modelo foi avaliado por meio de retroprojeção histórica antes da simulação futura. Mapas anteriores de uso da terra e variáveis associadas de drivers foram utilizados para simular um mapa observado posterior de uso da terra. A concordância entre as distribuições simuladas e observadas de uso da terra foi avaliada utilizando a acurácia geral (AG), o coeficiente Kappa e a Figura de Mérito (FoM). A validação foi realizada tanto no nível geral quanto por classes principais de uso da terra. A retroprojeção histórica de 2010 a 2020 apresentou uma acurácia geral (AG) de 93,4%, um coeficiente Kappa de 0,89 e uma Figura de Mérito (FoM) de 0,26, indicando uma capacidade altamente confiável de projeção espacial nas simulações subsequentes de múltiplos cenários.
Análises de incerteza e sensibilidade foram conduzidas para avaliar a robustez dos resultados da simulação. As análises de sensibilidade do PLUS examinaram os efeitos de diferentes configurações de resistência à transição e parâmetros de peso de vizinhança para as principais classes de uso da terra. As análises de sensibilidade do InVEST avaliaram a influência da variação dos coeficientes dos estoques de carbono entre os diferentes tipos de uso da terra. As análises de sensibilidade econométricas espaciais compararam especificações alternativas da matriz de pesos espaciais. Essas análises foram utilizadas para verificar se as classificações dos cenários e a direção dos principais efeitos da urbanização e do planejamento ecológico permaneceram consistentes sob diferentes configurações de parâmetros. Especificamente, a estabilidade dos resultados foi rigorosamente confirmada: as classificações dos cenários e os efeitos negativos de derrapagem espacial da urbanização da terra permaneceram inalterados quando os parâmetros de resistência à transição e os coeficientes dos estoques de carbono foram perturbados em ±15%.
5. Avaliação do estoque de carbono
O estoque de carbono foi avaliado utilizando o modelo de carbono InVEST. Parâmetros de referência para os quatro reservatórios de carbono (biomassa acima do solo, biomassa abaixo do solo, carbono orgânico do solo e matéria orgânica morta) foram atribuídos a cada tipo de uso da terra reclassificado, utilizando os valores biofísicos resumidos na Tabela 3. Mapas históricos de uso da terra (2000–2020) e conjuntos de dados matriciais de uso da terra futuros simulados pelo PLUS foram então importados para o modelo e integrados com os respectivos parâmetros de densidade de carbono. O modelo foi posteriormente executado para estimar o estoque total de carbono regional (Tg C), a densidade de estoque de carbono por grade (Mg C/ha) e mapas espaciais de mudança no estoque de carbono (ΔC) para os cenários histórico e futuro.
Os mapas de uso da terra reclassificados foram vinculados aos valores paramétricos correspondentes das reservas de carbono de referência relatados na Tabela 3. A sensibilidade das estimativas de estoque de carbono às incertezas dos parâmetros foi testada ajustando esses valores de referência em ±15%, com resultados detalhados de sensibilidade. O modelo foi utilizado para calcular o estoque total de carbono, a densidade de estoque de carbono e a alteração no estoque de carbono para cada classe de uso da terra sob condições históricas e cenários futuros de planejamento.
Três categorias de resultados foram avaliadas. Primeiro, o estoque total de carbono regional e as tendências temporais foram calculados para quantificar a magnitude e a direção da mudança no estoque de carbono ao longo do tempo. Segundo, as distribuições espaciais do estoque de carbono e da alteração no estoque de carbono foram mapeadas para identificar áreas de retenção e perda de carbono. Terceiro, as estimativas de estoque de carbono foram comparadas entre cenários de planejamento para avaliar os efeitos relativos das estratégias de proteção ecológica, proteção de terras agrícolas e otimização de baixo carbono na conservação do estoque de carbono.
Tabela 3: Parâmetros de reservatórios de carbono na linha de base para diferentes tipos de uso da terra. A tabela apresenta o carbono da biomassa acima do solo, o carbono da biomassa abaixo do solo, o carbono orgânico do solo, o carbono da matéria orgânica morta e a densidade total de carbono utilizados no módulo de estoque de carbono do InVEST. As unidades são Mg C/ha. Notas: Os valores representam os parâmetros de linha de base utilizados no modelo InVEST. A densidade total de carbono equivale à soma dos quatro reservatórios de carbono. Análises de sensibilidade que alteram esses valores de linha de base em ±15% são apresentadas na Tabela 6. Clique aqui para baixar esta tabela.
6. Autocorrelação espacial e análise econométrica espacial
A análise de autocorrelação espacial foi realizada para determinar se o estoque de carbono e a alteração no estoque de carbono exibiam dependência espacial significativa. O índice global de Moran foi calculado para avaliar o grau geral de agrupamento espacial na distribuição do estoque de carbono e na alteração do estoque de carbono em toda a área de estudo. O índice de Moran local foi então utilizado para identificar padrões locais de associação espacial, incluindo agrupamentos alto-alto, baixo-baixo, alto-baixo e baixo-alto. A análise de hotspots foi realizada para identificar áreas com perda concentrada de carbono e retenção do estoque de carbono.
Modelos econométricos espaciais foram utilizados para examinar as relações entre urbanização da terra, planejamento ecológico, condições ambientais, fatores socioeconômicos e dinâmicas do estoque de carbono. A densidade do estoque de carbono ou a alteração no estoque de carbono serviram como variável dependente. As variáveis explicativas incluíram o índice de urbanização da terra, proporção de área construída, intensidade de planejamento ecológico, índice de vegetação por diferença normalizada (NDVI), densidade populacional, densidade de PIB, densidade de estradas, elevação, declividade, precipitação média anual e temperatura média anual.
A regressão de mínimos quadrados ordinários (OLS) foi utilizada como modelo de referência. A dependência espacial dos resíduos foi avaliada antes da estimação dos modelos autorregressivos espaciais (SAR), de erro espacial (SEM) e de Durbin espacial (SDM). O desempenho dos modelos e as estimativas dos coeficientes foram comparados entre as especificações.
Uma matriz espacial padronizada por linha foi construída para representar as relações de vizinhança entre unidades administrativas. A especificação principal baseou-se na contiguidade espacial, e análises de robustez compararam matrizes espaciais alternativas baseadas em distância e nos vizinhos mais próximos, quando disponíveis.
Os efeitos diretos, indiretos e totais foram calculados a partir do modelo de equações estruturais espaciais (SDM) para avaliar as relações locais e de derramamento espacial. Variáveis de terreno, acessibilidade, vegetação, clima e variáveis socioeconômicas foram incluídas como controles para reduzir o viés de variável omitida. Os coeficientes do modelo foram interpretados como associações espaciais condicionais, e não como efeitos causais definitivos. Variáveis com valores-p acima dos limiares convencionais de significância foram interpretadas como evidência fraca ou sugestiva e não foram consideradas efeitos estatisticamente robustos.