Tesis doctoral de Iñigo Errasti Arrieta
Los modelos climáticos acoplados atmósfera-océano son sofisticadas herramientas diseñadaspara simular el sistema climático terrestre y las complejas interacciones entre sus componentes.Hoy en día más de una docena de centros meteorológicos en el mundo desarrollan modelosde clima para mejorar nuestra comprensión del clima y el cambio climático, dando soporte alas actividades desarrolladas por el panel intergubernamental para el cambio climático (ipcc)de las naciones unidas.Sin embargo, estos modelos no son perfectos. Nuestra comprensión del clima de la tierrano es perfecta y por tanto en los modelos se debe asumir una serie de suposiciones ysimplificaciones, que introducen por ejemplo biases (desviaciones) en las simulaciones delclima, a veces difíciles de subsanar o eliminar. Aun así, en los últimos años se ha mejorado lacomprensión teórica del clima, ajustado la base física de los modelos climáticos que lo simulane incrementado igualmente el número y la calidad de las observaciones, aumentando muysignificativamente la potencia de cálculo necesaria para llevar a cabo las integraciones de losmodelos.En particular, los modelos climáticos de circulación general (gcms) están basados en lasecuaciones fundamentales del movimiento, la termodinámica y otras propiedadesmacroscópicas de la atmósfera, océanos y capa de hielo, obtenidas de las leyes básicas de lafísica. Estos principios básicos son particularmente la conservación de momento, masa yenergía, una ecuación de conservación de salinidad y humedad y ciertas ecuaciones de estadoque relacionan parámetros como presión, temperatura, densidad y salinidad.Los modelos gcms usan un conjunto diverso de esquemas numéricos y de parametrización,que parametrizan por ejemplo, los correspondientes términos de convección, flujos de capalímite, transferencia radiativa, condensación y flujos superficiales sobre tierra. Estasparametrizaciones pretenden relacionar el estado climático que representa el modelo en unaresolución espacial por debajo de la escala espacial de ciertos procesos físicos, como puedenser la microfísica, la humedad o el transporte de energía por convección en pequeños cúmulos,la influencia de las nubes en los flujos de radiación o los flujos de energía y humedad debidosa transporte turbulento en la capa límite. De hecho, estas parametrizaciones son la única formade representar estos procesos en código numérico, implementado en ordenadores que van adesarrollar las integraciones de los modelo climáticos, ya que no se pueden resolver de formaexplícita con el tamaño de celda del modelo.La idea general de esta tesis es que los modelos de clima son herramientas que se usanpara entender mejor las complejas interacciones del sistema climático. Los experimentos congcms ofrecen la aproximación más adecuada para mejorar nuestro conocimiento de losmecanismos físicos subyacentes en climas pasados y futuros. Van a permitir evaluarparcialmente el impacto del aumento de las concentraciones de gases invernaderos de origenxivprefacioxvantropogénico en el clima de la tierra, y sus estimaciones deben ser tenidas en cuenta en lafase de atribución del problema del cambio climático. La atmósfera terrestre está siendosometida a unos cambios sin precedentes en la historia, y aunque variaciones en el clima tanconsiderables como el actual sí han sido registradas en el pasado, pocas han acontecido a lavelocidad observada en el cambio climático actual. Aun así, los modelos climáticos sonherramientas imperfectas, y por tanto, para perfeccionarlos es fundamental comparar sussalidas (estimaciones) con datos (observacionales y de reanálisis) que reflejan el clima actual.Dado que las ecuaciones atmosféricas no son lineales, la evolución en el tiempo de lasolución o integración es muy sensible a las condiciones iniciales del sistema. Si dosrealizaciones (runs) de un mismo sistema son inicializadas desde dos condiciones inicialesmuy próximas pero distintas, las dos soluciones divergirán pronto. El crecimiento del error enlas condiciones iniciales depende claramente de las propias condiciones iniciales. Por tanto,la predictibilidad varía y hace la predicción no válida tras un número finito de tiempo. Estehecho puede ser explicado usando el modelo matemático original desarrollado por lorenz(1963) para explicar el caos, conocido popularmente como efecto mariposa. Lorenz desarrollóesta teoría a partir de una forma simplificada de las ecuaciones de movimiento de la atmósfera.La propiedad del límite de la predictibilidad meteorológica sugirió entonces la necesidad deconsiderar conjuntos (ensembles) que muestran la predictibilidad del sistema climático. Estaultima está a su vez relacionada con la predictibilidad de las funciones densidad deprobabilidad de los elementos climáticos y no con las realizaciones individuales de cadatrayectoria en el espacio de fases del sistema.Las predicciones de los modelos gcm no son sólo sensibles a las condiciones iniciales,sino también a las incertidumbres de la parametrización del modelo. De hecho, el sistemaclimático global no es completamente predecible en términos de las condiciones inicialesexistentes. El problema de la sensibilidad a las condiciones iniciales es normalmente evitadomediante el uso de un conjunto (ensemble) de predicciones. Esta técnica consiste en ejecutarun modelo de predicción varias veces a partir de condiciones iniciales que difieren encantidades consistentes con las incertidumbres en el estado inicial.Sin embargo, las incertidumbres también aparecen como consecuencia de una incompletarepresentación de la realidad en la formulación de los diferentes modelos gcms. No todoslos modelos usan los mismos esquemas numéricos de parametrización para representarmatemáticamente los mismos procesos físicos. Como consecuencia, modelos diferentesinicializados con las mismas condiciones iniciales formularán resultados diferentes tras ciertotiempo de predicción (palmer, 1999).Otra fuente de incertidumbre en los modelos y por tanto también en sus simulaciones, sonlas condiciones de contorno. Los modelos cgm dependen de entradas (forzamientos) querepresentan influencias externas naturales como son el ciclo solar o el vulcanismo terrestre yantropogénicas, principalmente emisiones de gases invernadero cuya evolución está asociadaa factores tecnológicos y socio-económicos no predecibles.Así, que un modelo climático simule adecuadamente el clima actual no garantiza que generepredicciones fiables para el clima futuro. A pesar de estas dificultades y limitaciones, laconcordancia de los modelos con las observaciones del clima actual es la única forma parapoder tener confianza en un modelo, asumiendo que un modelo que describe con precisión elclima actual hará probablemente una mejor proyección del clima futuro.Prefacioxvien esta tesis, se van a usar las simulaciones de los modelos gcms obtenidos del climatemodel intercomparison project (cmip), proyecto que archiva de una manera organizada losdatos proporcionados por los modelos gcms más recientes y desarrollados. Dichos datos sondistribuidos por el program for climate model diagnosis and intercomparison (pcdmi). Elprimer proyecto de intercomparación de modelos de clima acoplados fue llevado a cabo amediados de los años noventa del siglo pasado. A esta primera generación de modelos de climase le denominó cmip1, y a la siguiente disponible a inicios del siglo actual, cmip2. Los datosusados en este trabajo son los correspondientes a la tercera generación de modelos cmip3,coordinado por el world climate research program (wcrp). Estos modelos están incluidosen el cuarto informe de evaluación (ar4) del ipcc y por tanto, son denominados modeloscmip3 o modelos ar4 (ipcc ar4, 2007).Los datos obervacionales usados en este estudio de tesis han sido obtenidos a partir de dosproyectos de reanálisis y de una base de datos de precipitación. En los años noventa del siglopasado, el national center for environmental prediction (ncep), en cooperación con elnational center for atmospheric research (ncar) creó un registro de análisis globales dediferentes variables climáticas mediante el análisis de datos meteorológicos observados ymodernos métodos de asimilación de datos y modelización. Dicho proyecto es conocido conel nombre de reanálisis ncep/ncar (kalnay et al., 1996; kistler et al., 2001). El proyectode reanálisis europeo era40 (simmons and gibson, 2000; uppala et al., 2005) complementaal reanálisis ncep/ncar usando otros esquemas de asimilación de datos.Igualmente, hoy en día un gran numero de bases de datos de observaciones de variablesclimáticas como son la presión a nivel de mar, la temperatura superficial o la precipitaciónestán disponibles para la comunidad científica gracias al trabajo realizado por instituciones declima y meteorológicas. Ejemplo de estos centros, son el met office hadley centre encolaboración con climatic research unit (cru) de la universidad de east anglia (jones etal., 1999), el goddard institute for spaces studies (hansen et al., 2010) o por citar el caso dela base de datos de precipitación utilizada en esta tesis, el global precipitation climatologyproject (gpcp) descrito por huffman et al. (1995) y adler et al. (2003).Las simulaciones de modelos ar4 seleccionadas en este trabajo de tesis para lacomparación con el clima actual fueron llevadas a cabo con un conjunto realista deforzamientos externos observados en la atmósfera, que incluyen la historia estimada o conocidade un rango dado de emisiones antropogénicas y también naturales, como son las variacionesen la radiación solar, la actividad volcánica, los gases invernadero y aerosoles sulfato. Laformulación exacta de estos forzamientos varía de un modelo a otro con implicacionespotenciales en el comportamiento de cada modelo.Los datos de las simulaciones de clima, son para el siglo xx. El experimento o conjuntode integraciones 20c3m, usa condiciones iniciales que parten de una simulación de controlpreindustrial, y su integración se realiza con el forzamiento histórico observado de gases deefecto invernadero, aerosoles, actividad volcánica y radiación solar durante el siglo x . Lasconcentraciones de co2 varían desde 286 a 367 ppm. Cada simulación 20c3m se inicializaa partir de diferentes estados climáticos obtenidos en la simulación de control preindustrial yes forzada con las concentraciones históricas de gases invernadero y aerosoles del periodo1860-2000. En la simulación de control preindustrial, la concentración de co2 es la de losniveles previos a la era industrial (simulación de 500 años anteriores a 1860 con unaconcentración de co2 constante de 275 ppm).Prefacioxviipara el siglo xxi, las simulaciones de los modelos han sido clasificadas atendiendo a losdiferentes escenarios de emisión de gases invernadero (nakicenovic et al., 2000). Estosescenarios corresponden a concentraciones de gases invernadero relativamente altos (escenarioa2), moderado (a1b) y bajo (b1) a final de siglo, y pueden ser considerados extensiones deintegraciones previas (simulación de control preindustrial y experimento 20c3m).Los actuales modelos de clima acoplados han avanzado poco en la predicción a nivelregional de los cambios climáticos de una forma fiable y unánime. Esta falta de progreso estáíntimamente ligada con el papel que juega la dinámica atmosférica en la conformación deciertos aspectos de las respuesta climáticas, bien a pequeña escala que afectan a laestratificación atmosférica o a otros efectos diabáticos, o bien a la escala sinóptica, cuyacomplejidad inherente e interacción no lineal, también va a limitar la probabilidad depredicciones regionales más exactas. Es interesante, por tanto, profundizar en estudios depredicción a escala regional, por ejemplo mediante técnicas de downscaling o regionalizaciónen trabajos futuros.Sin embargo, para poder desarrollar dichos estudios de downscaling, es necesario evaluarinicialmente el comportamiento relativo de los diferentes modelos de clima sobre el área deestudio. Nosotros vamos a desarrollar un estudio de evaluación de modelos en la penínsulaibérica, para posteriormente ampliar la región de estudio al atlántico norte y finalmente a unaescala global. Aun así, se ha de tener en cuenta y tener presente, que ciertos estudios hanmostrado que modelos que reproducen bien el clima observado histórico sobre una regióngeográfica no reproducen bien el clima en otra región del globo. En otras palabras, unaevaluación o ranking de modelos ar4 hecha para una región puede no ser válida para otraregión.Estudios previos han mostrado que, para temperaturas globales, el incremento observadoen las últimas décadas del siglo xx, se debe mayoritariamente a la actividad humana sobre elplaneta basada en la combustión de combustibles fósiles y en el incremento de gasesinvernadero inducido por esta combustión en la atmósfera. Cambios en la temperatura yprecipitación en la tierra, al fin y al cabo pueden provocar un serio impacto sobre lasactividades del hombre en muchas regiones. Estos modelos gcms no son unánimes en señalarla distribución de los cambios regionales en precipitación.Por tanto, el área de investigación en la que queda enmarcada esta tesis, es hoy en día devital importancia debido a razones sociales, económicas, científicas, y puede permitir obtenerunos resultados para ser publicados en revistas científicas de renombre.El grupo de investigación eolo de meteorología, clima y medio ambiente donde eldoctorando ha desarrollado este trabajo de tesis, tiene experiencia reconocida enpredictibilidad de temperatura sobre la península ibérica usando modelos gcms (frías et al.,2005). Igualmente tiene amplia experiencia en el análisis de la variabilidad de varios elementosclimáticos a nivel regional (ezcurra et al., 2008) y en la comparación con los datos salida demodelos gcms (sáenz et al., 2001a; sáenz et al., 2001b y fernández et al., 2003).Particularmente, durante los últimos años, este grupo ha liderado uno de los proyectos másimportantes en la estrategia de investigación vasca sobre cambio climático (eklima21),dentro del programa etortek. Una parte importante de esta estrategia está dedicada al temacentral de esta tesis (evaluación de modelos climáticos ar4), y por tanto, este trabajo deinvestigación es parte esencial del proyecto eklima21.Prefacioforman parte de este grupo de investigación, agustín ezcurra, gabriel ibarra-berastegi,jon sáenz, javier díaz de argandoña, isabel herrero e iñigo errasti, profesores de losdepartamentos de física aplicada ii, ingeniería nuclear y mecánica de fluidos, y expresióngráfica y proyectos de ingeniería de las facultades de ciencia y tecnología, farmacia yescuelas de ingeniería de la universidad del país vasco (upv/ehu).
Datos académicos de la tesis doctoral «Evaluation and diagnosis of coupled climate models included in the ipcc fourth assessment report«
- Título de la tesis: Evaluation and diagnosis of coupled climate models included in the ipcc fourth assessment report
- Autor: Iñigo Errasti Arrieta
- Universidad: País vasco/euskal herriko unibertsitatea
- Fecha de lectura de la tesis: 13/12/2013
Dirección y tribunal
- Director de la tesis
- Jon Saenz Aguirre
- Tribunal
- Presidente del tribunal: eugenio Ruiz urrestarazu
- concepcion Rodriguez puebla (vocal)
- alejandro Oscoz abad (vocal)
- fabienne Lohou — (vocal)