Global Warming Has Been Accelerating Since At Least 1990

Resumen Técnico: El calentamiento global se ha acelerado desde al menos 1990

Planteamiento del problema
La literatura reciente ha reportado que la Temperatura Media de la Superficie Global (GMST, por sus siglas en inglés) se ha acelerado, con algunos estudios sugiriendo que esta aceleración comenzó alrededor de 2015. Sin embargo, las metodologías existentes para detectar la aceleración poseen limitaciones. La especificación de tendencia cuadrática estándar ($xt = \alpha + \lambda t + \beta t^2 + \epsilon_t$) impone una forma funcional rígida donde la tasa de calentamiento aumenta a un ritmo constante, lo cual puede no reflejar la realidad. Además, los modelos cuadráticos a menudo requieren procedimientos de dos pasos (filtrar los datos para efectos de El Niño y solares antes de la estimación de la tendencia), lo que puede reducir la precisión de las estimaciones de los parámetros. Por el contrario, las tendencias lineales por tramos permiten rupturas estructurales pero asumen cambios discretos en fechas específicas, lo que podría malinterpretar una tasa que aumenta continuamente como una serie de rupturas. El artículo busca abordar estos problemas mediante el desarrollo de un marco estadístico más flexible para probar tendencias supralineales (mayores que lineales) y determinar el inicio de la aceleración de manera más robusta.

Metodología
Los autores proponen una especificación linearitmica como una alternativa flexible a la tendencia cuadrática. El modelo se define como:
$$x_t = \alpha + \beta \log(t)t + \epsilon_t$$
donde $x_t$ es la anomalía de la GMST en el tiempo $t$. La presencia de aceleración se prueba evaluando la hipótesis $\beta > 0$. Teóricamente, el término linearitmico $\log(t)t$ crece más rápido que una tendencia lineal pero más lento que cualquier tendencia polinómica con un grado $\delta > 1$, proporcionando un límite inferior para la aceleración supralineal.

El estudio utiliza cinco de los principales conjuntos de datos de temperatura (HadCRUT, GISTEMP, NOAA, Berkeley Earth y ERA5) que cubren el período desde enero de 1970 hasta diciembre de 2025. El análisis emplea un procedimiento de estimación de un solo paso que incluye variables de control para el efecto de El Niño ($N_t$) y la irradiancia solar ($S_t$) directamente dentro de la regresión, en lugar de utilizar un enfoque de filtrado de dos pasos. Los modelos se estiman mediante Mínimos Cuadrados Ordinarios (OLS) con errores estándar consistentes con heterocedasticidad y autocorrelación (estimador de Newey-West) para dar cuenta de la correlación serial y la heterocedasticidad.

Para determinar el momento de la aceleración, los autores emplean un marco de ventana expansiva, comenzando con una muestra de 1970–1990 y extendiendo incrementalmente la fecha final hasta 2025. Para abordar el problema de las pruebas múltiples inherente a este enfoque, utilizan un procedro de bootstrap silvestre de max-|t| para calcular los valores críticos y ajustar los p-valores, asegurando la robustez contra falsos positivos.

Contribuciones clave

1. Especificación estadística novedosa: El artículo introduce la especificación linearitmica ($\beta \log(t)t$) como una prueba robusta de aceleración. A diferencia del modelo cuadrático, no asume una tasa constante de aumento en la tasa de calentamiento y evita los problemas de multicolinealidad que suelen encontrarse entre los términos lineales y cuadráticos.
2. Perspectiva teórica sobre las rupturas estructurales: Los autores proporcionan una prueba teórica (Teorema 0.3) que demuestra que, si la tendencia subyacente real es supralineal, las pruebas de ruptura estructural estándar eventualmente detectarán cambios en la pendiente de un modelo de tendencia lineal. Esto sugiere que las "rupturas estructurales" reportadas en la GMST (por ejemplo, en la década de 1970 o en 2015) pueden ser artefactos estadísticos de ajustar un modelo lineal a una tendencia supralineal, en lugar de evidencia de cambios discretos en la tasa de calentamiento.
3. Refinamiento del tiempo de la aceleración: Al aplicar el marco de ventana expansiva con ajustes de bootstrap, el estudio reevalúa el inicio de la aceleración, desafiando la noción de que es un fenómeno puramente post-2015.

Resultados

• Especificación linearitmica: Los resultados indican evidencia estadísticamente significativa de aceleración ($\beta > 0$) en los cinco conjuntos de datos utilizando la especificación linearitmica. Esta significancia se mantiene tanto con como sin variables de control adicionales. Crucialmente, el análisis de ventana expansiva muestra que la evidencia de aceleración está presente desde al menos 1990, y la significancia estadística se fortalece a medida que se incluyen datos más recientes (hasta 2025).
• Especificación cuadrática: En contraste, la especificación cuadrática arroja resultados significativos solo en las ventanas de estimación más largas y principalmente cuando se incluyen variables de control adicionales. Sin controles, el coeficiente cuadrático es significativo en solo tres de los cinco conjuntos de datos.
• Estimación del exponente polinómico: Un ejercicio adicional de estimación del exponente $\delta$ en una tendencia polinómica general ($t^\delta$) sugiere que el exponente real es menor que dos (el caso cuadrático) pero mayor que uno. Esto respalda el uso del modelo linearitmico como un límite inferior más apropiado para la dinámica observada.
• Rupturas estructurales: Los resultados teóricos sugieren que la detección de rupturas estructurales en modelos de tendencia lineal es consistente con una tendencia supralineal subyacente, ofreciendo una explicación alternativa para los puntos de cambio previamente reportados en las tendencias de la GMST.

Significancia y afirmaciones
El artículo afirma que la aceleración del calentamiento global no es meramente un fenómeno reciente que comienza alrededor de 2015, sino que ha estado presente durante varias décadas, siendo detectable desde al menos 1990. Los autores argumentan que la especificación linearitmica proporciona un marco más robusto y flexible para probar la aceleración en comparación con el rígido modelo cuadrático.

El estudio destaca que la aceleración observada implica una tasa de calentamiento creciente en el tiempo, lo que tiene implicaciones para los ecosistemas, la salud humana y la economía. Los autores señalan que sus hallazgos son consistentes con otras evidencias físicas, como el aumento de la fracción aérea de las emisiones de carbono antropogénicas y la duplicación del desequilibrio energético de la Tierra en las últimas décadas.

Finalmente, el artículo enfatiza la importancia de distinguir entre tendencias supralineales y rupturas estructurales para la política y la modelización. Si la tendencia real es supralineal, los modelos lineales por tramos pueden subrepresentar la verdadera dinámica de la GMST al capturar únicamente cambios discretos en la pendiente en lugar de una tasa que aumenta continuamente. Los autores concluyen que la presencia de aceleración subraya la urgencia de las estrategias de mitigación y adaptación.