Señal frente a ruido: cuándo un cambio en el carbono del suelo es real

Un comprador pregunta al asesor de un agricultor de tierras de cultivo si el programa de cultivos de cobertura está secuestrando carbono. El asesor informa de +0.4 toneladas de C por hectárea y año durante cinco años, con 60 muestras pareadas por momento temporal. El verificador pide el intervalo de confianza. Cruza el cero. La afirmación se descarta en silencio, el presupuesto de MRV de la temporada se da por perdido, y el comprador vuelve al equipo de cadena de suministro sin nada que mostrar.

El patrón está ahora bien documentado. Bradford y colaboradores (2023) pusieron a prueba la viabilidad de detectar empíricamente el cambio de COS en 45 campos comerciales de cultivo y hallaron que la detección de una acumulación significativa (en torno a 3 Mg C/ha en diez años) a escala de campo individual no es fiable ni siquiera con densidades de muestreo de 1.2 hectáreas por muestra1. La señal existe. Solo que no es mayor que el ruido a la escala que suponen la mayoría de los diseños de campo. El colapso en 2024 de la sonada emisión australiana de créditos de pastoreo, en la que se emitieron unos 250,000 ACCU en gran parte por un año húmedo de decil 10 y no por un cambio de manejo, dejó claro lo mismo a escala de política pública2.

La detección del cambio de COS es un problema de relación señal-ruido. De un lado, tres fuentes de ruido. Del otro, una señal de manejo pequeña. El tamaño de muestra, el horizonte temporal y las decisiones de diseño se derivan de esa desigualdad; no son el punto de partida.

Anatomía del umbral de ruido

El ruido de la medición del COS es la varianza del cambio en las reservas estimado en torno a su valor verdadero. Se acumula a partir de tres fuentes, y la mayor rara vez es la que preocupa a los profesionales.

Varianza espacial dentro de la unidad de inferencia

El carbono del suelo es heterogéneo a todas las escalas. Dentro de una sola parcela de cultivo de aspecto homogéneo, Poeplau y colaboradores (2022) documentan desviaciones estándar a escala de parcela de 5 a 8 Mg C por hectárea en las reservas de 0–30 cm, incluso al volver a muestrear el mismo perfil con un desplazamiento posicional mínimo, un coeficiente de variación en torno al 10 por ciento que es en su mayor parte irreducible3. En paisajes heterogéneos, donde varían la posición en la ladera, el material parental y el historial de manejo, son habituales CV del 30 al 60 por ciento14. En perennes tropicales, las desviaciones estándar publicadas para las reservas de 0–30 cm en agroforestería de cacao se agrupan en torno a 10 a 14 Mg C por hectárea, determinadas por la posición en la ladera, la arquitectura de sombra y la edad de la cronosecuencia5.

La varianza espacial es la fuente de ruido que la mayoría de los proyectos presupuestan por debajo de lo debido, porque la pasada piloto suele muestrear unas pocas hectáreas de apariencia uniforme. La varianza se despliega cuando el diseño se extrapola a toda la unidad de manejo. Si una nueva zona está siquiera dentro del dominio del modelo es una cuestión relacionada que abordamos en ¿Representativo de qué?

Varianza temporal y estacional

Incluso en el mismo punto de la misma finca, muestrear de nuevo a la misma profundidad en fechas distintas arroja valores distintos. Dominan dos efectos. Primero, la humedad del suelo y la reserva de materia orgánica recién mineralizada fluctúan con las precipitaciones, motivo por el cual la emisión de créditos australiana desapareció el año siguiente al fin de las condiciones de triple La Niña2. Segundo, la profundidad de muestreo en el campo es una operación estocástica, no un plano fijo: unos pocos centímetros de diferencia entre la línea base y el remuestreo pueden alterar las reservas de 0–30 cm en varios Mg C por hectárea cuando la densidad aparente cambia con la labranza reducida o una enmienda orgánica. Wendt y Hauser (2013) introdujeron la corrección por masa de suelo equivalente (ESM) precisamente para absorber este artefacto, y Raffeld et al. (2024) muestran que la elección entre profundidad fija y ESM puede alterar el cambio inferido más que la propia señal de manejo con los mismos datos brutos67.

Varianza analítica y de procesamiento

El laboratorio también añade varianza. Incluso en un único laboratorio acreditado, las mediciones replicadas sobre la misma muestra tienen un CV del 1 al 3 por ciento para el carbono total por combustión seca. La varianza entre laboratorios y entre métodos (combustión seca frente a Walkley-Black, deriva de calibración del infrarrojo medio (MIR), diferencias en la preparación de la muestra) es mayor y rara vez se rastrea hasta la cifra final. Heuvelink y colaboradores (2021) sostienen que la varianza de predicción de los modelos de aprendizaje automático puede dominar la estimación final del cambio en las reservas a escala espacial fina, y que sin una propagación explícita de los tres términos de varianza, una cifra publicada de cambio de COS es estructuralmente ininterpretable8.

En una línea base típica de tierras de cultivo, las tres aportan en torno al 80, 15 y 5 por ciento del total. En un diseño pareado de línea base y remuestreo en el que se revisitan las mismas parcelas, σ²_spatial se desploma (el componente entre parcelas se cancela), y σ²_temporal más σ²_analytical pasan a ser el ruido dominante. La proporción entre ambos es la palanca central del diseño, y la razón por la que los diseños pareados son ahora el estándar en los regímenes de crédito serios.

La señal que intentamos detectar

El otro lado de la desigualdad es la señal impulsada por el manejo. La literatura es por fin lo bastante rica como para fundamentar rangos realistas, y son menores de lo que reconoce la mayoría de las presentaciones comerciales.

• Cultivos de cobertura. El metaanálisis de 37 sitios de Poeplau y Don halla 0.32 ± 0.08 Mg C por hectárea y año en el suelo superficial9. Honesto, real, pequeño.
• Siembra directa y labranza reducida. Las ganancias aparentes de 0.3 a 0.7 Mg C por hectárea y año en 0–15 cm desaparecen en gran medida al integrar a 0–40 cm. Powlson y colaboradores (2014) mostraron que esto es principalmente un artefacto de redistribución en profundidad, no un secuestro neto10.
• Agroforestería templada. Mayer et al. (2022) reportan una media en el suelo superficial de 0.21 Mg C por hectárea y año en 51 estudios, con una desviación estándar entre estudios de 0.79; los setos vivos alcanzan 0.32, el cultivo en callejones 0.26 y el silvopastoreo un valor ligeramente negativo11. La DE entre estudios es mayor que la media, que es el problema central.
• Agroforestería de cacao. Adiyah et al. (2022, 2023) documentan una gran variación según la clase de edad y la posición en la ladera, con tasas de acumulación en el suelo superficial de 0.2 a 0.5 Mg C por hectárea y año y ganancias sustanciales en el subsuelo en rodales de más de ~18 años5.
• Biochar y enmiendas. Eventos discretos y puntuales de 5 a 15 Mg C por hectárea seguidos de una atenuación lenta; en la práctica, un escalón más que una pendiente.
• Pastoreo adaptativo multipotrero. Dependiente del sitio y del equilibrio; diferencias medias del orden de 0.5 a 1.5 Mg C por hectárea y año donde existen, pero con una gran varianza entre sitios.

De ello se derivan dos implicaciones. La señal de manejo realista se sitúa en el rango de 0.1 a 0.5 Mg C por hectárea y año para las prácticas de tipo pendiente, menor en términos absolutos que el error de laboratorio en una sola muestra. Y la DE entre sitios del efecto, en los metaanálisis, es comparable o mayor que la media, lo que limita cuánta compresión de varianza puede extraer un solo proyecto de la literatura previa.

La desigualdad de detectabilidad

Para un diseño pareado de línea base y remuestreo entre dos momentos temporales, el cambio medio más pequeño que una campaña puede distinguir estadísticamente de cero es la diferencia mínima detectable (DMD):

con z_α/2 = 1.96 para α = 0.05 (bilateral) y z_β = 0.84 con una potencia del 80 por ciento, de modo que (z_α/2 + z_β)² ≈ 7.85. σ_paired es la desviación estándar del cambio por parcela entre momentos temporales (mucho menor que la σ absoluta de las reservas si se revisitan las mismas parcelas y la autocorrelación intraparcelaria es alta). n es el número de parcelas pareadas.

Tres propiedades de esta desigualdad importan para el diseño.

La DMD es simétrica. Reducirla a la mitad cuadruplica la n requerida. Duplicar la σ también cuadruplica la n. No hay una vía gradual desde una campaña que detecta 1 Mg C por hectárea y año hasta otra que detecta 0.25; cuesta dieciséis veces la muestra.

σ_paired se desploma donde σ no lo hace. Por eso los diseños pareados son ahora la estructura creíble para la verificación. Un diseño pareado de 60 parcelas con σ_paired = 4 Mg C por hectárea puede detectar un cambio acumulado a 5 años de unos 2 Mg C por hectárea con una potencia del 80 por ciento. Las mismas 60 parcelas, muestreadas de forma independiente en la línea base y en el remuestreo, necesitarían una σ más cercana a 2 Mg C por hectárea para alcanzar la misma DMD, algo que no se logra en tierras de cultivo reales.

La señal se acumula con el tiempo, no solo con el número de muestras. Una acumulación de 0.3 Mg C por hectárea y año llega a 1.5 Mg C por hectárea en cinco años y a 3 Mg C por hectárea en diez. La misma n que no detecta el cambio en el año cinco puede detectarlo con holgura en el año diez sin añadir una sola muestra. Esta es la paradoja del horizonte temporal: la forma más barata de comprar poder de detección es esperar.

Simulador de la envolvente de la señal

El primer diagrama concreta la desigualdad. Elija una tasa de acumulación anual, una desviación estándar pareada y un número de parcelas; el gráfico traza la trayectoria media verdadera y el corredor de ruido en ±ε, y marca el año en que la trayectoria supera por primera vez el corredor. Baje σ con un diseño pareado y el año de detección se reduce a la mitad; suba la tasa de acumulación y la trayectoria sale antes; recorte n y el corredor se ensancha hasta que la trayectoria no puede escapar de él en 15 años.

El simulador permite ver por qué la mayoría de las afirmaciones publicadas del tipo «secuestramos X» son inverificables con sus propios datos: los controles que resuelven una señal creíble quedan fuera de lo que el diseño de campo realmente pagó.

La paradoja del horizonte temporal

La implicación de la desigualdad de la DMD es incómoda. Para la mayoría de las prácticas de tipo regenerativo, la ventana de detección honesta es de 8 a 12 años, no los 3 a 5 años sobre los que suelen construirse las hojas de ruta corporativas y los ciclos de crédito112. Smith (2004) modeló esto hace dos décadas: con un incremento del 20 al 25 por ciento en los aportes de carbono, la detección con un 90 por ciento de confianza requiere aproximadamente de 6 a 10 años si el régimen tolera un desplazamiento del 3 por ciento en el COS de fondo, y es esencialmente indetectable en ciclos más cortos si el régimen solo tolera un 15 por ciento del cambio de fondo12.

De ello se derivan tres consecuencias para el diseño de proyectos. Primero, aceptar una contabilidad híbrida de medición y modelización para los primeros años en lugar de sobreafirmar. La consulta de VM0042 v3.0 de Verra (el Soil Sampling and Analysis Handbook publicado en marzo de 2026) y el argumento económico de medición y remedición de Potash et al. (2025) apuntan ambos en esta dirección1314. Segundo, estructurar el contrato en torno a la ventana de verificación, no a la de marketing: un informe a 3 años que dice «por debajo del límite de detección, trayectoria de la señal coherente con el modelo» es más duradero que un informe a 3 años que dice «+1.4 Mg C por hectárea» con un IC que cruza el cero. Tercero, tomarse en serio la cuestión de la profundidad. Las ganancias de COS en perennes de raíces profundas de más de ~18 años suelen superar las del suelo superficial; un diseño de 0–30 cm sencillamente no las ve57.

El diseño de muestreo como palanca de reducción de la varianza

Si σ está fijada por el suelo y esperar más tiempo está limitado por el contrato, la palanca restante es el diseño. Tres movimientos de diseño se combinan, y un proyecto serio usa los tres.

Estratificación. Asignar el esfuerzo de muestreo a estratos que capturan varianza (topografía, edad del manejo, tipo de suelo) reduce la σ efectiva y por tanto la n. de Gruijter et al. (2016) introdujeron el algoritmo Ospats para una asignación estratificada de coste óptimo usando un mapa previo de COS, y mostraron que la ganancia es sustancial cuando la información previa es buena15. Empíricamente, un trabajo posterior sobre los datos de Bradford et al. muestra que estratificar con covariables fácilmente disponibles rinde una reducción media del tamaño de muestra del 17 por ciento (IC del 95% del 11 al 23 por ciento) frente al muestreo aleatorio simple. Información previa sólida más un muestreo doblemente equilibrado puede llevar eso al 30 a 50 por ciento1. Obsérvese el mínimo: VM0042 v2.2 fija un mínimo de 15 muestras por estrato, lo que puede borrar la ganancia de Neyman en unidades de manejo pequeñas.

Línea base y remuestreo pareados. Volver a las mismas parcelas en el remuestreo desploma el componente de varianza entre parcelas. Las reducciones de ruido típicas de un diseño pareado son del 30 al 70 por ciento según la autocorrelación temporal del COS del micrositio. La Red de Monitoreo del Carbono del Suelo del USDA NRCS lo pone en práctica con seis muestras por campo de 4 hectáreas en un calendario de diseño pareado, aproximadamente una muestra por cada 0.7 hectáreas16.

Compuesta dentro del estrato, independiente entre estratos. Agrupe las muestras dentro de un estrato para absorber la varianza a microescala en una estimación de media del estrato a menor coste; mantenga los estratos como réplicas independientes para que la prueba conserve grados de libertad. Lo contrario (agrupar entre estratos en una sola muestra compuesta) es rápido y barato e inútil para la inferencia.

Selector de diseño de muestreo

El segundo diagrama compone las tres palancas de diseño. Un diseño base aleatorio simple con σ fijada por el usuario produce una estimación de años hasta la detección. La estratificación reduce la σ efectiva en un factor realista según la literatura. La línea base y el remuestreo pareados desploman aún más la varianza entre parcelas. El gráfico actualiza los años hasta la detección, el recuento por estrato tras aplicar el mínimo de 15 al estilo de VM0042 y el recuento total de muestras resultante.

El objetivo del selector no es optimizar hacia un único número. Es poner de manifiesto dónde rinde cada decisión de diseño, y dónde el mínimo de verificación (la restricción de 15 por estrato, el requisito de profundidad, la corrección ESM) cuesta más de lo que Neyman ahorra.

Qué cambia el giro hacia la fusión de modelo y datos

El cambio más trascendental en 2024–2026 ha sido el paso de un crédito basado solo en la medición o solo en el modelo hacia diseños híbridos. La arquitectura está ya bien descrita: modelos de Nivel 3 basados en procesos o empíricos calibrados localmente con redes de medición pareadas, con la varianza de predicción y la varianza de parámetros propagadas conjuntamente al intervalo de incertidumbre reportado1714. El documento de consulta de VM0042 v3.0 de Verra y el GHG Protocol Land Sector and Removals Standard (en vigor el 1 de enero de 2027) exigen ambos datos empíricos y específicos por reservorio para cualquier remoción de carbono del suelo que se reporte, poniendo fin al resquicio del «solo modelo» que usaron varios programas de primera generación1318. La revisión independiente de la guía para compradores de CarbonPlan concluyó en 2023 que ningún protocolo disponible comercialmente cumplía a la vez requisitos adecuados de rigor, durabilidad, adicionalidad y profundidad; el ciclo de normas de 2026 es, en parte, una respuesta explícita19.

Esto importa para el marco de señal frente a ruido. Un modelo por sí solo tiene una predicción determinista con una cota de incertidumbre estructuralmente optimista. Una campaña de medición por sí sola tiene el umbral de ruido descrito arriba. Un modelo calibrado con mediciones puede comprimir el intervalo de incertidumbre por debajo de cualquiera de los dos componentes, pero solo si las mediciones son lo bastante densas como para restringir los parámetros que importan y el error estructural del modelo se incluye en la varianza propagada. Los híbridos de «modelo ligero, medición ligera» no son creíbles; la «medición como calibración del modelo» es el estándar operativo hacia el que convergen los regímenes de 2026.

Plantillas de reporte honesto

La disciplina que hace operativo un marco de señal frente a ruido está en el reporte. Cuatro números deberían viajar juntos con cualquier afirmación de cambio de COS, y un verificador de la era de 2026 los pedirá.

1. La estimación puntual del cambio, con unidades y profundidad (p. ej., «cambio de COS de 0–30 cm de +1.4 Mg C por hectárea en 5 años, corregido por ESM»).
2. El intervalo de confianza del 95 por ciento calculado a partir de la varianza real del diseño, no de un CV genérico de la literatura.
3. La DMD que el diseño puede resolver con una potencia del 80 por ciento, calculada antes de la campaña, y el resultado expresado como «afirmación dentro / fuera del límite de detección».
4. La taxonomía de incertidumbre propagada siguiendo las convenciones del Volumen 4 (AFOLU) del IPCC20: varianza de los datos de actividad, varianza de las reservas de referencia y varianza del cambio en las reservas o del error del modelo, combinadas de forma aditiva.

Una afirmación de +1.4 Mg C por hectárea con un IC de −0.2 a +3.0 y una DMD de 1.6 es honesta, útil e inverificable; se consigna así en el informe y el proyecto planifica la siguiente remedición. Una afirmación de +1.4 Mg C por hectárea con un IC de +0.6 a +2.2 y una DMD de 0.8 es honesta, útil y defendible; entra en el inventario.

Qué se sigue para el diseño de proyectos

El marco de señal frente a ruido simplifica varias decisiones de contratación que, a lo largo de los últimos cinco años, se han litigado repetidamente como si fueran una cuestión de opinión.

La unidad de inferencia es el proyecto, no el campo. La detección por campo de las prácticas de tipo pendiente no es estadísticamente fiable a una densidad asequible, mientras que la detección de la media del proyecto en 20 a 60 campos pareados es operativamente viable. Diseñe en consecuencia, y encamine la afirmación del comprador al nivel que los datos pueden defender.

La ventana de detección honesta es el horizonte de manejo, no la hoja de ruta corporativa. Para los cultivos de cobertura, la siembra directa y la mayoría de los tipos de agroforestería, eso son de 8 a 12 años desde la línea base. Los ciclos de reporte de tres años siguen siendo útiles, siempre que reporten frente a la trayectoria de la DMD y no frente a una señal presunta.

El perfil profundo ya no es opcional. Las ganancias en el subsuelo en sistemas perennes de más de ~18 años a menudo superan las del suelo superficial y son invisibles para un diseño de 0–30 cm5721. La corrección ESM es igualmente innegociable en las comparaciones entre labranza y siembra directa67.

El protocolo contra el que se audita al comprador fija el mínimo. VM0042 v2.2/v3.0, el CRCF de la UE (Reglamento 2024/3012), el GHG Protocol Land Sector and Removals Standard y SBTi FLAG convergen en requisitos de datos empíricos más estrictos a lo largo de 2026–202713182223. Diseñar hacia la norma más permisiva con la que el equipo está familiarizado, cuando el comprador será auditado contra la más estricta, es un modo de fallo a cámara lenta.

Ideas clave