Scandal: BoM thermometer records fiddled “by month” — mysterious square wave pattern discovered
What would you say if you knew that the official Perth thermometer was accurate at recording minimums for most of time in October in the eighties, but 0.7°C too warm all of December and 1.2°C too cool in January. Bizarrely that same thermometer was back to being too warm in February? Try to imagine what situation could affect that thermometer and require post hoc corrections of this “monthly” nature. Then imagine what could make that same pattern happen year after year.
Lo que se sugiere es que el proceso de homogeneización introduce desviaciones importantes en las tendencias de temperatura. Para ello trata algunos ejemplos, entre ellos el de la serie de temperatura de la ciudad de Perth. El caso es que solo se sugiere que hay desviaciones importantes, pero no se demuestra.
El escándalo del título se refiere al BoM, el Bureau of Meteorology del Gobierno de Australia. Como en la entrada del blog se señalan donde se pueden obtener los datos de temperatura medidos directamente (raw ) y los homogeneizados, he descargado algunos para ver qué salía.
He descargado datos de 3 estaciones (en total hay datos de 112, pero no era cuestión de descarglas todas): Perth (codigo 09021) que es de la que más hablan, y otras dos al azar, simplemente porque me sonaban: Sydney (codigo 66062) y Darwin (14015).
Los datos de temperaturas maximas y mínimas diarias, una vez homogeneizadas se pueden obtener de la dirección:
http://www.bom.gov.au/climate/change/acorn-sat/#tabs=Data-and-networks . Hay un mapa con las 112 estaciones y pinchando en cada una de ellas se puede acceder a los datos y descargarlos en formato texto.
Para los datos sin homogeneizar se puede utilizar el formulario de la dirección:
http://www.bom.gov.au/climate/data/ , que los enseña en pantalla, pero tiene una opción de descarga en un fichero comprimido. Los códigos que identifican las estaciones son los mismos en los dos casos.
Lo que he hecho es simplemente calcular el ajuste lineal con la temperatura media diaria (media entre maxima y minima) con los datos homogeneizados y sin homogeneizar para los últimos 5, 10, 15, 20, 25, 30, 35, 40, 45 y 50 años. Lo que he obtenido es:
(el valor de Diferencia es la tendencia homogeneizada - tendencia raw, así una diferencia positiva es que hay mayor calentamiento en la serie homogeneizada que en la raw, y negativa mayor calentamiento en la raw.
[1] "Estacion 09021. Perth"
[1] "Ultimos 05 años. Tendencia ºC/Año Raw: -0.24003 . Homogeneizada: -0.24003. Diferencia 0.00000 "
[1] "Ultimos 10 años. Tendencia ºC/Año Raw: 0.08335 . Homogeneizada: 0.08242. Diferencia -0.00093 "
[1] "Ultimos 15 años. Tendencia ºC/Año Raw: 0.04411 . Homogeneizada: 0.04495. Diferencia 0.00084 "
[1] "Ultimos 20 años. Tendencia ºC/Año Raw: 0.00361 . Homogeneizada: 0.00908. Diferencia 0.00547 "
[1] "Ultimos 25 años. Tendencia ºC/Año Raw: 0.01761 . Homogeneizada: 0.02550. Diferencia 0.00789 "
[1] "Ultimos 30 años. Tendencia ºC/Año Raw: 0.01839 . Homogeneizada: 0.02160. Diferencia 0.00320 "
[1] "Ultimos 35 años. Tendencia ºC/Año Raw: 0.01834 . Homogeneizada: 0.01825. Diferencia -0.00009 "
[1] "Ultimos 40 años. Tendencia ºC/Año Raw: 0.01445 . Homogeneizada: 0.01310. Diferencia -0.00135 "
[1] "Ultimos 45 años. Tendencia ºC/Año Raw: 0.01738 . Homogeneizada: 0.01611. Diferencia -0.00127 "
[1] "Ultimos 50 años. Tendencia ºC/Año Raw: 0.01887 . Homogeneizada: 0.01901. Diferencia 0.00014 "
[1] "Estacion 66062 Sydney "
[1] "Ultimos 05 años. Tendencia ºC/Año Raw: 0.02291 . Homogeneizada: 0.02291. Diferencia 0.00000 "
[1] "Ultimos 10 años. Tendencia ºC/Año Raw: -0.00265 . Homogeneizada: -0.00320. Diferencia -0.00055 "
[1] "Ultimos 15 años. Tendencia ºC/Año Raw: 0.00303 . Homogeneizada: 0.00325. Diferencia 0.00022 "
[1] "Ultimos 20 años. Tendencia ºC/Año Raw: 0.03161 . Homogeneizada: 0.03160. Diferencia -0.00001 "
[1] "Ultimos 25 años. Tendencia ºC/Año Raw: 0.03081 . Homogeneizada: 0.03082. Diferencia 0.00000 "
[1] "Ultimos 30 años. Tendencia ºC/Año Raw: 0.02458 . Homogeneizada: 0.02481. Diferencia 0.00022 "
[1] "Ultimos 35 años. Tendencia ºC/Año Raw: 0.01954 . Homogeneizada: 0.01711. Diferencia -0.00243 "
[1] "Ultimos 40 años. Tendencia ºC/Año Raw: 0.01740 . Homogeneizada: 0.01476. Diferencia -0.00264 "
[1] "Ultimos 45 años. Tendencia ºC/Año Raw: 0.01838 . Homogeneizada: 0.01653. Diferencia -0.00185 "
[1] "Ultimos 50 años. Tendencia ºC/Año Raw: 0.01889 . Homogeneizada: 0.01786. Diferencia -0.00102 "
[1] "Estacion 14015 Darwin"
[1] "Ultimos 05 años. Tendencia ºC/Año Raw: 0.04210 . Homogeneizada: 0.04210. Diferencia 0.00000 "
[1] "Ultimos 10 años. Tendencia ºC/Año Raw: 0.00932 . Homogeneizada: 0.00931. Diferencia -0.00001 "
[1] "Ultimos 15 años. Tendencia ºC/Año Raw: 0.01214 . Homogeneizada: 0.01198. Diferencia -0.00015 "
[1] "Ultimos 20 años. Tendencia ºC/Año Raw: 0.00288 . Homogeneizada: 0.00322. Diferencia 0.00034 "
[1] "Ultimos 25 años. Tendencia ºC/Año Raw: -0.00078 . Homogeneizada: -0.00069. Diferencia 0.00009 "
[1] "Ultimos 30 años. Tendencia ºC/Año Raw: -0.00726 . Homogeneizada: -0.00721. Diferencia 0.00005 "
[1] "Ultimos 35 años. Tendencia ºC/Año Raw: -0.00412 . Homogeneizada: -0.00410. Diferencia 0.00001 "
[1] "Ultimos 40 años. Tendencia ºC/Año Raw: 0.00026 . Homogeneizada: 0.00059. Diferencia 0.00033 "
[1] "Ultimos 45 años. Tendencia ºC/Año Raw: 0.00033 . Homogeneizada: 0.00177. Diferencia 0.00144 "
[1] "Ultimos 50 años. Tendencia ºC/Año Raw: 0.00428 . Homogeneizada: 0.00625. Diferencia 0.00197 "
Es decir, la máxima diferencia es de 0.0078 ºC/año ( 0.07ºC/década) para algunos datos de Perth, aunque la gran mayoría de las diferencias son del orden de 0.001 ºC/año, es decir una centésima por década, dentro del margen de error que se puede esperar. Y lo que es más llamativo, algunas diferencias salen negativas, es decir, en algunos casos la homogeneización supone menos calentamiento que los datos sin homogeneizar.
Así que no veo que la homogeneización esté introduciendo desviaciones importantes en los datos. Naturalmente, para poder sacar una conclusión fiable sería necesario un estudio más amplio (y desde luego mucho más serio que el que hace el blog).
En la propia web del Bureau of Meteorology australiano (
http://www.bom.gov.au/climate/change/acorn-sat/#tabs=Adjustments) hay una discusión sobre este tema, y un gráfico en el que se muestran los resultado de las series homogeneizadas y sin homogeneizar: las diferencias en los últimos 50 años son casi inapreciables.