Bovinos por km² e estabelecimentos com adubo orgânico nos municípios de São Paulo

Português · English

Prancha Nº 114, a quarta de cinco da coleção Bioenergia cruzada; a anterior é a prancha 113; a seguinte, a prancha 115.

A coleção cruza a bioenergia paulista com o que se vende nos postos, com o campo e com as prefeituras, nos 645 municípios de São Paulo. Os pares saíram de uma triagem de correlações entre o potencial de biometano do CP2B, as vendas de combustíveis do Anuário da SEMIL, a silvicultura, o rebanho, o Censo Agropecuário, a frota e a MUNIC. Cada par passou por correlação de postos bruta e parcial (descontadas urbanização, densidade e renda), amostra efetiva corrigida pela autocorrelação espacial, valores q de Benjamini–Hochberg e LISA bivariado. A última prancha é descritiva, sem teste, e o texto diz isso.

As chamadas do mapa apontam Ferraz de Vasconcelos, no cinturão de hortas da metrópole, onde 81 % dos estabelecimentos usaram adubo orgânico com 2,2 bois por km², e Presidente Bernardes, no Pontal, com 202 bois por km² e 16 %: um dos 136 municípios de pouco adubo cercados de muito boi.

O que vem nos arquivos

Camada (GPKG)	Geometria	Conteúdo
municipios	polígono	645 municípios com cana e etanol, biometano e diesel, frota, mata nativa e madeira para celulose, bovinos e adubo orgânico, distâncias ao gasoduto, gás canalizado e as respostas da MUNIC sobre resíduos; classes b111, b113, b114, arranjo115 e aglomerados lisa111, lisa113, lisa114
gasodutos	linha	gasodutos de transporte e rede de distribuição (tipo)
regioes_intermediarias, limite_estadual	polígono	referência

O GPKG é o mesmo nas cinco pranchas e traz um estilo do QGIS para cada uma. O CSV tem a tabela dos municípios, sem geometria. Campos vazios são falta de dado, não zero. Os aglomerados LISA vêm como AA (alto-alto), BB (baixo-baixo), AB, BA e ns (não significativo); a primeira letra é a variável do município, a segunda a dos vizinhos.

Método

1. bovinos_km2: efetivo de bovinos em 31/12/2024 (IBGE, PPM, tabela 3939) sobre a área do município.
2. estab_agro, estab_adubo_organico e pct_estab_adubo_organico: estabelecimentos agropecuários, os que fizeram adubação orgânica ou química e orgânica, e a parte, em % (IBGE, Censo Agropecuário 2017, SIDRA 6847).
3. Terços de bovinos: até 22,0, de 22,0 a 59,7 e mais de 59,7 cabeças por km²; de adubo orgânico: até 20,1 %, de 20,1 % a 33,3 % e mais de 33,3 %. b114: classe bivariada, dezena = terço de bovinos (1 a 3), unidade = terço de adubo orgânico (1 a 3). De baixo para cima e da esquerda para a direita, as nove casas têm 37, 54, 114; 70, 77, 58; 98, 75, 34 municípios. No estado, 28,7 % dos estabelecimentos usaram adubo orgânico; 28 municípios ficam sem dado, quase todos na Grande São Paulo e no litoral: a PPM não publica rebanho bovino neles, e em alguns o Censo Agropecuário também não tem dado. Falta de dado não é zero.
4. n_kg_ha_agro: nitrogênio dos dejetos dos rebanhos por hectare de área agropecuária (prancha 081). Pelos quartos de nitrogênio, a mediana do adubo orgânico cai de 30,8 % para 21,6 %.
5. Adubo orgânico × bovinos por km² (log): correlação de postos de Spearman −0,41 (617 municípios), parcial −0,32, amostra efetiva 351, p < 0,001, q < 0,001. lisa114: aglomerados LISA bivariados (adubo orgânico no município × bovinos nos vizinhos): 46 alto-alto, 62 baixo-baixo, 113 alto-baixo e 136 baixo-alto; I de Moran bivariado −0,39. Baixo-alto por região intermediária: São José do Rio Preto (50), Presidente Prudente (36) e Araçatuba (26).
6. Correlação parcial: correlação entre os resíduos dos postos de cada variável regredidos nos postos dos controles, que são urbanização (parte de moradores em situação urbana), log da densidade demográfica e log da renda por habitante (Censo 2022). Amostra efetiva sob autocorrelação espacial: n·(1 − Iₓ·Iᵧ)/(1 + Iₓ·Iᵧ), com Iₓ e Iᵧ o I de Moran dos resíduos numa vizinhança dos 6 vizinhos mais próximos; valores q de Benjamini–Hochberg sobre todos os pares da triagem. LISA bivariado: valor do município contra a média dos 6 vizinhos na outra variável, 999 permutações, semente 20261019, p < 0,05.

Limitações

• O Censo Agropecuário é de 2017 e o rebanho, de 2024: sete anos separam as duas medidas.
• Adubo orgânico não é só esterco: inclui composto e outros resíduos, e o Censo conta estabelecimentos, não toneladas nem hectares. O pasto quase não recebe adubo, e o esterco do gado solto fica no campo, sem ser recolhido; o mapa mostra onde o esterco existe e onde ele é usado como adubo, não quanto dele se perde.
• Leitura ecológica: o mapa compara municípios, não pessoas nem fazendas. Uma correlação entre municípios não diz que as mesmas propriedades ou os mesmos motoristas têm as duas características, nem que uma causa a outra.