Introducción

Los cambios en el clima de la Tierra a lo largo de su existencia pueden deberse a variaciones en la dinámica interna del sistema del clima o bien por forzamientos externos. Los forzamientos pueden ser de origen natural, como variaciones en la actividad solar o en la actividad volcánica, o variaciones de origen humano, a través del aumento de emisiones de gases de efecto invernadero y cambios en el uso del suelo, entre otros factores.

Cuando no existe variación en los forzamientos externos, los cambios del clima son influenciados sólo por los procesos de la dinámica interna del sistema. De este modo se tiene conocimiento de cambios del clima de corto o largo plazos relacionados con sistemas de interacción océano-atmósfera, como en el caso de la variabilidad climática asociada con el fenómeno de El Niño-Oscilación del Sur (ENSO, por sus siglas en inglés), o de las variaciones en periodos más largos, como las vinculadas con la Oscilación Decadal del Pacífico (PDO, por sus siglas en inglés) y la Oscilación Multidecadal del Atlántico (AMO, por sus siglas en inglés), o incluso con grandes erupciones volcánicas.

A fin de identificar tendencias de cambios del clima se realizan evaluaciones de “detección” para demostrar que el clima o el sistema afectado por el clima han cambiado en un sentido estadístico sin explicar el origen o la “atribución” de dichos cambios. Es necesario realizar investigaciones detalladas adicionales para determinar la atribución a causa alguna. Un cambio identificado es detectado en observaciones instrumentales si se determina que es pequeña su probabilidad de ocurrencia debido sólo a la variabilidad interna (^{Hegerl et al., 2010}).

De esta forma se tiene que el calentamiento global detectado en la tendencia de incremento de la temperatura media del planeta es uno de los indicadores de un cambio climático en nuestro tiempo. En informes recientes de evaluación, el Panel Intergubernamental de Expertos sobre el Cambio Climático (IPCC, por sus siglas en inglés) concluyó que “El calentamiento del sistema climático es inequívoco” (^{Trenberth et al., 2007}). Además de registrar un calentamiento multidecadal notable, la temperatura media global en superficie muestra una variabilidad decadal e interanual considerable (^{IPCC, 2013}, ²⁰¹⁴).

Así, la influencia antropogénica en el clima ha sido detectada y atribuida robustamente a escala global, pero para muchas aplicaciones es más útil una estimación de la contribución antropogénica en las tendencias de temperatura en una región en particular (^{Hegerl et al., 2007}), en este caso para México. La detección y atribución del cambio climático a escala continental

o más pequeña es más difícil de detectar que a escala global debido a la mayor contribución de la variabilidad interna en escalas pequeñas en relación con la respuesta al forzante antropogénico, la mayor dificultad para distinguir entre diferentes factores que causan el cambio a escala pequeña, y errores incrementales en la representación de detalles regionales por parte de modelos del clima que omiten forzantes que pueden ser importantes a menor escala, como el cambio de uso de suelo o el carbono negro (^{Bindoff et al., 2013}).

En el presente trabajo se detectan y analizan las tendencias positivas de la temperatura media a nivel local en diferentes regiones de México, a partir del procesamiento de series históricas de la temperatura registrada en 112 estaciones climáticas selectas para el periodo 1950-2013, con respecto al periodo base 1961-1990. En el presente estudio se consideró tomar en cuenta estudios climáticos en México previos, donde en mayor medida se realizaron esfuerzos sobre la calidad de la información (^{Sigro, 2012}; ^{Bravo, Aspra, Zarraluqui, Gay, & Estrada, 2012}; ^{Méndez, Ramírez, Cornejo, Zarate, & Cavazos, 2010}; ^{Pavia, Graef, & Reyes, 2006}; ^{Englehart & Douglas, 2002}, ²⁰⁰⁴; ^{Pieyns, 2003}; ^{Mantua & Hare, 2002}).

Las estaciones elegidas cumplieron criterios objetivos de control de calidad aplicados a 1 099 estaciones identificadas como las que se han empleado en estudios climáticos para México (^{Bravo et al., 2012}; ^{Sigro, 2012}; ^{Méndez et al., 2010}; ^{Pieyns, 2003}), a partir, la mayoría de ellas, de la base de datos climática nacional del Climate Computing Project (CLICOM), iniciativa promovida internacionalmente por la Organización Meteorológica Mundial (OMM) desde 1985 (^{OMM, 2015}), que en México es administrada por el Servicio Meteorológico Nacional (SMN).

Datos y métodos

Bases de datos y metodología

La base de datos climáticos del CLICOM (^{OMM, 2015}) utilizada, actualizada a julio de 2015, cuenta con registros de 5 031 estaciones que reportaron temperatura mínima diaria y 5 034 que reportaron temperatura máxima diaria, aunque el número de estaciones con información en una misma fecha es mucho menor, similar a lo encontrado por ^{Bravo, Aspra, Zarraluqui, Gay y Estrada (2012)}.

Se conjuntó una relación de 2 252 estaciones climáticas a partir de las utilizadas en siete investigaciones identificadas para México (cuadro 1). Los estudios considerados han evaluado regiones climáticas con características similares entre las estaciones climáticas que las representan y al mismo tiempo características diferentes entre regiones (^{Bravo et al., 2012}; ^{Sigro, 2012}; ^{Pieyns, 2003}), así como la teleconexión de patrones climáticos en el país con oscilaciones estacionales, interanuales, decadales y multidecadales, como El Niño o la PDO, entre otros (^{Méndez et al., 2010}).

Se cotejaron en forma cruzada los grupos de estaciones de los siete estudios, pues éstas no son excluyentes y algunas están duplicadas entre las investigaciones (cuadro 1). A partir del análisis se identifica que:

• 1 099 estaciones son independientes.

• 1 153 se duplican entre los grupos, lo cual implica que están contenidas como parte de las 1 099 estaciones independientes.

Cuadro 1 Resultados de la comparación de las estaciones utilizadas en estudios climáticos. 

Estudio climático	Estaciones duplicadas	Estaciones independientes	Total
1. Monitor de Sequías de México (MSM)	0	358	358
2. Sigro, 2012	276	10	286
3. Dr. Douglas	74	139	213
4. Bravo et al., 2012	131	218	349
5. Méndez et al., 2010	310	271	581
6. Pieyns, 2003	116	59	175
7. Análisis propio	246	44	290
Total	1 153	1 099	2 252

Después de identificar el conjunto de 1 099 estaciones climáticas independientes se obtuvieron y procesaron los datos diarios de las mismas para las variables de temperatura máxima (Tmx) en °C y temperatura mínima (Tmn) en °C de la base de datos del CLICOM para México.

La temperatura media diaria se calculó como el promedio de la temperatura máxima y mínima diarias, tal como lo recomienda la Organización Meteorológica Mundial (^{OMM, 2011}):

Tmed=Tmx+Tmn2 (1)

A partir de los datos diarios del CLICOM se calcularon valores mensuales y anuales; estos últimos corresponden a la escala temporal de análisis considerada en el artículo. El criterio considerado es que el dato mensual debe contener más del 80% de datos diarios y para calcular el dato anual se debe tener más del 80% de los meses con datos en cada año; en caso contrario se considera como dato perdido. En la figura 1 se presenta el total de años con que disponen de información de la temperatura media las 1 099 estaciones climáticas. En varias partes del país, las estaciones tienen 50 o más años de información, excepto algunas zonas de Guerrero, Chihuahua, Coahuila, Durango, Sonora, Sinaloa y norte de la Península de Yucatán, donde las estaciones cuentan con menos de 50 años con información.

Figura 1 Representación de la cobertura temporal de información de la temperatura media de 1 099 estaciones climáticas. Los puntos indican la ubicación de las estaciones. 

Se procesó el periodo 1950 a 2013 como el rango temporal común para el análisis de tendencias debido a que la mayoría de las series de tiempo en las estaciones climáticas contaron con información continua y con los menores vacíos a partir de 1950.

Se aplicaron criterios de control de calidad a las 1 099 estaciones climáticas, con el fin de garantizar que la información posea buena calidad, distribución espacial y suficientes observaciones como se sugiere para estudios climáticos (^{Bravo et al., 2012}; ^{Méndez et al., 2010}; ^{Alexander et al., 2006}). Los criterios aplicados para considerar las estaciones en la mayoría de los siete estudios analizados se enfocan en la disponibilidad mínima de información diaria, mensual y anual en la serie.

En tanto que los criterios de control de calidad aplicados en el presente estudio (cuadro 2) son similares a los sugeridos por el Grupo de Expertos en Detección e Índices de Cambio Climático (ETCCDI, por sus siglas en inglés) coordinado por la Comisión de Climatología de la OMM (CCl/OMM) (^{Alexander et al., 2006}), los considerados por ^{Sigro (2012)}, y ^{Méndez, Ramírez, Cornejo, Zarate y Cavazos (2010)}. Asimismo, si se presentan cinco o más días consecutivos con datos perdidos, se consideró el mes como dato perdido (^{OMM, 2011}).

Cuadro 2 Criterios de control de calidad considerados y resultados de su aplicación. 

Criterio	Datos sospechosos identificados	Datos sustituidos
Temperatura mínima es mayor o igual que la máxima	2 524	2 524
Diferencia mayor que 20 °C entre la temperatura mínima y máxima, en combinación con una:
1. Diferencia interdiaria mayor o igual que 20 °C con respecto a los dos días adyacentes	170	140
2. Diferencia interdiaria mayor o igual que 25 °C con respecto a uno de los dos días adyacentes	126	12
Total	2 820	2 676

De la aplicación de los criterios de control de calidad a las 1 099 estaciones climáticas se sustituyeron como datos perdidos 2 676 datos diarios de 2 820 datos sospechosos identificados (cuadro 2).

En la selección de las estaciones para la detección de la tendencia de la temperatura media local en regiones de México se aplicaron criterios objetivos presentados en el cuadro 3. A cada criterio se le ponderó una calificación en función del grado de rigurosidad de su aplicación, siendo 3 (gris claro, denotado como gris 3) para el más riguroso y 1 (gris más oscuro, denotado como gris 1) para el más flexible. Como parte de los criterios se considera el periodo climatológico base 1961-1990 y se excluyen de los análisis anomalías para un margen de cambio de ± 0.25 °C, como una aproximación que considera indirectamente potenciales errores instrumentales.

Cuadro 3 Criterios aplicados para la selección de estaciones a considerar en el análisis de tendencias. 

Criterio/semáforo	1. Longitud temporal de la serie	2. Año final de la serie	3. Disponibilidad de información anual en el periodo 1961-1990	4. Anomalía con respecto al periodo 1961-1990
3. Gris 3	Igual o mayor que 60 años	Igual o mayor a 2010	100% de los años	Mayor que 0.25 °C
2. Gris 2	Entre 50 y 60 años	Entre 2000 y 2010	Entre 90% y menor a 100%	Menor que -0.25 °C
1. Gris 1	Menor que 50 años	Menor a 2000	Menor a 90% de los años	Entre -0.25 y 0.25 °C

El número de estaciones que cumplen la secuencia de cada uno de los cuatro criterios aplicados, agrupadas por la calificación obtenida, se muestra en la figura 2. En la aplicación secuencial de los cuatro criterios se consideraron las estaciones con las dos calificaciones más altas (gris 3 y 2).

Figura 2 Total de estaciones agrupadas por calificación (o semáforo) obtenida durante la aplicación secuencial de los cuatro criterios de selección y de la revisión de las gráficas de las series de anomalías. 

Se realizó un análisis cualitativo de las series de anomalía de la temperatura media de las 383 estaciones que cumplieron los cuatro criterios, con el fin de realizar una evaluación más profunda de las estaciones que no presentaran lo siguiente: secuencias de tres o más años vacíos en la parte intermedia de las series, ausencia de ciclo anual y situaciones explícitas de inhomogeneidad. Posterior al análisis, se identificaron 112 estaciones climáticas para evaluar la tendencia de la temperatura local en México.

Resultados

De los registros históricos analizados se tiene que 112 estaciones climáticas cumplen con la calidad y estándares requeridos para el presente estudio; de éstas, 84 corresponden a anomalías positivas y 28 a negativas (figura 3); en esta proporción, 3 de cada 4 estaciones señalan un calentamiento, en comparación con 1 de cada 4, que señala un enfriamiento. Este hallazgo coincide con los publicados por el ^{IPCC (2013)} en el sentido que la Tierra experimenta un proceso de calentamiento global donde algunas regiones muestran diferentes condiciones o niveles de calentamiento, o incluso condiciones de enfriamiento en otras.

Figura 3 Distribución de estaciones de acuerdo con la intensidad de la anomalía (a) positiva (cuadros rellenos) y (b) negativa (círculos sin rellenar), con respecto al periodo base 1961-1990. Se incluye un gráfico de la serie anual con baja influencia continental (estación 12099) y alta influencia (estación 08044). 

De la figura 3 se observa que las estaciones con señal de calentamiento y enfriamiento están intercaladas, resultados consistentes con lo presentado por ^{Hartmann et al. (2013)} para estaciones climáticas de no sólo los países vecinos, sino para todo el orbe. Lo anterior probablemente esté relacionado con factores como el entorno urbano/rural de las estaciones, diferencias en altitud, patrones físicos locales del clima o el mantenimiento a la infraestructura de medición.

El 48.8% de las estaciones con señal de calentamiento presenta un incremento moderado menor a 1 °C; 41.7% de incremento alto, con valores entre 1 y 2 °C; y el 9.5% de incremento muy alto, con valores mayores a 2 °C. De las estaciones que indican enfriamiento: 53.6% presenta un decremento moderado mayor que -1 °C y 46.4% con un decremento alto, con valores entre -1 y -2 °C. Destaca que no se identificaron estaciones con anomalías de muy bajo decremento, menor que -2 °C.

Con el objeto de confirmar la hipótesis inicial, en el sentido de que la tendencia positiva de cambio climático está presente tanto en zonas urbanas como rurales, a partir de las 84 estaciones climáticas identificadas con tendencia positiva se consideraron 20 estaciones selectas para el análisis de detección de tendencias de calentamiento. Las estaciones se seleccionaron en función de la distribución espacial, la elevación en el país, la longitud temporal de las series y el entorno ambiental actual en el que se encuentran ubicadas: a) urbano, con población mayor a 100 000 habitantes o b) rural con población menor a 100 000 habitantes en el año 2010 (figura 4).

Figura 4 Representación de 20 estaciones selectas para el análisis de tendencia positiva de temperatura. Los números se asocian con los de la primera columna del recuadro. 

Se consideró la población para el año 1980 y 2010 en los municipios urbanos (recuadro en la figura 4), con el fin de inferir la potencial contribución de la isla de calor urbana a la señal del calentamiento local. El IPCC señala que es improbable que el efecto de isla de calor urbana y el cambio de uso de suelo hayan contribuido en más del 10% de la tendencia de la temperatura media global de la superficie. Sin embargo, en algunas regiones con rápido desarrollo, los impactos de la isla de calor urbana y el cambio de uso de suelo pueden ser sustancialmente elevados en las tendencias regionales de temperatura (^{Bindoff et al., 2013}).

A las 20 estaciones climáticas selectas se les aplicó un análisis de homogeneidad, complementario al control de calidad, mediante la herramienta “R-Standard Normal Homogeneity Test–rsnht.r” (^{Aguilar, 2010}) para la detección y corrección de “brincos” en la serie de datos mensuales, relacionados con el cambio de instrumento de medición por desperfecto, cambio de observador o de ubicación de la estación (y con ello el entorno), e inconsistencia en las series, validados con respecto a los registros en papel (metadatos) del acervo histórico del SMN.

Asimismo, y con el fin de identificar la potencial influencia de las características continentales sobre la tendencia de la temperatura de las 20 estaciones selectas se aplicó un índice de continentalidad termo-pluviométrica de Gorezyñski-Valle Melendo (^{del Valle, 1991}):

Donde A se refiere a la amplitud térmica (diferencia entre la temperatura media del mes más cálido y la del mes más frío, en °C); %P.V. es el porcentaje de la precipitación acumulada de junio, julio y agosto, con respecto al acumulado anual; %P.I. es el porcentaje de la precipitación acumulada de diciembre, enero y febrero con respecto al acumulado anual; Lat es la latitud geográfica de la estación climática en radianes.

A partir de los resultados, se encontró que es reducida la variabilidad interanual de la temperatura en estaciones con valores bajos del índice (cuadro 4) y ésta aumenta en estaciones que presentan mayor influencia continental, no obstante que algunas se encuentren en zonas costeras. Un rasgo interesante de continentalidad obedece al tipo de sistemas atmosféricos que inciden en la región y al patrón estacional de la precipitación donde está la estación climática. Por ejemplo, los valores del índice de continentalidad termo-pluviométrica en las estaciones ubicadas en la vertiente del Golfo de México muestran valores más altos que aquellos de la vertiente del Pacífico, sin que esto sea una regla. La interpretación de estos valores es que la vertiente del Golfo de México está influenciada por sistemas meteorológicos de tipo continental, como frentes fríos, y se da una reducción de precipitación en la temporada de lluvia, conocida como canícula, situaciones que ocurren con menos regularidad en el Pacífico.

La señal de calentamiento por influencia antropogénica está presente en la tendencia positiva de las series de la temperatura de las 20 estaciones climáticas, con independencia del grado de influencia continental en las estaciones; no obstante, la intensidad del calentamiento es mayor en estaciones del centro y norte del país (recuadros en la figura 3), con alta influencia continental. Estos hallazgos son congruentes con las proyecciones de cambio climático para las próximas décadas, los cuales señalan que los mayores cambios de temperatura en México se presentarán hacia el norte del país, siendo de alta intensidad en el interior del continente (Tercera y Cuarta Comunicación Nacional de México ante la Convención Marco de las Naciones Unidas sobre el Cambio Climático) (^{Semarnat-INE, 2006} y ²⁰⁰⁹).

Las series de tiempo de la temperatura media en las 20 estaciones presentan variaciones caracterizadas por periodos fríos y cálidos intercalados que responden a la influencia, entre otros, de oscilaciones como El ENSO, la PDO y la AMO, la topografía, las variaciones en tipos de vegetación y uso de suelo, cercanía con el mar, fenómenos físicos como el desplazamiento latitudinal de la ZCIT, la actividad interanual de los ciclones tropicales y del monzón de Norteamérica, así como la actividad de la corriente en chorro, los frentes fríos y los sistemas convectivos de mesoescala (^{Bravo et al., 2012}; ^{Méndez et al., 2010}; ^{Englehart & Douglas, 2002}, ²⁰⁰⁴).

En las figuras 5 y 6 se presentan las gráficas de las series de la temperatura y su tendencia para estaciones de caso, con respecto al periodo base 1961-1990. Del análisis de las series se identifica que:

Figura 5 (a) Series de la anomalía de temperatura media anual en entornos urbanos con respecto al promedio del periodo 1961-1990. Región noroeste: estación 02033 en Mexicali, Baja California; región norte: estación 08044 en Delicias, Chihuahua; región centro: estación 09043 en Gustavo A. Madero, Ciudad de México; región Golfo de México: estación 28111 en Tampico, Tamaulipas; región Sur: estación 07205 en Comitán de Domínguez, Chiapas; (b) serie temporal de anomalías de temperatura relacionadas con El Niño y representación de las fases de la AMO y PDO. Fuente: elaborada a partir de ^{Englehart y Douglas (2002)}, y de datos de NCEP disponibles en http://www.cpc.ncep.noaa.gov/data/indices//uri]. 

Figura 6 (a) Series de la anomalía de temperatura media anual en entornos rurales con respecto al promedio del periodo 1961-1990. Región noroeste: estación 02038 en Tijuana, Baja California; región norte: estación 10021 en Inde, Durango; región centro: estación 01018 en San José de Gracia, Aguascalientes; región Golfo de México: estación 27030 en Macuspana, Tabasco; región sur: estación 12099 en Benito Juárez, Guerrero; (b) serie temporal de anomalías de temperatura relacionadas con El Niño y representación de las fases de la AMO y PDO. Fuente: elaborada a partir de ^{Englehart y Douglas (2002)}, y de datos de NCEP disponibles en http://www.cpc.ncep.noaa.gov/data/indices//uri]. 

Del análisis de las tendencias de la temperatura media, con respecto al periodo base 1961-1990, en algunas de las estaciones de caso de las figuras 5 y 6, con respecto a las oscilaciones de la PDO, AMO y El Niño, se identifica que:

Como parte del análisis de la tendencia de calentamiento se generaron curvas de la función de distribución de probabilidad (FDP) para la temperatura media, ajustadas con la función de Gauss (^{Bidegain & Diaz, 2011}):

Donde x es la temperatura media, µ es el valor medio y σ es la desviación estándar.

A partir de la curva de Gauss, calculada para tres periodos climatológicos de 30 años, se compara y analiza la evolución de la temperatura media con respecto al periodo base 1961-1990. En la figura 7 se presentan, a modo de ejemplo, los resultados para la estación rural 01018 en Aguascalientes y la estación urbana 09043 en la Ciudad de México, como casos representativos de las 20 estaciones selectas.

Figura 7 Curvas de la FDP de Gauss para la temperatura media en a) la estación rural 01018 en San José de Gracia, Aguascalientes, y en b) la urbana 09043 en Gustavo A. Madero, Ciudad de México. 

De la comparación de los tres periodos climatológicos de 30 años se identifica que: a) la media de la temperatura promedio se ha desplazado a valores más cálidos en ambas estaciones, con mayor incremento en la estación con entorno urbano con respecto a la rural, lo que reafirma los resultados descritos en la presente sección e indica que muy probablemente es un indicio de un cambio de clima de cada periodo; b) los máximos de la temperatura media se han desplazado hacia valores más cálidos en las dos estaciones, en tanto que los mínimos se han desplazado con mayor amplitud en la estación urbana con respecto a la rural, la primera en el orden de 1.5 °C y la segunda en menos de 0.5 °C, y c) la variabilidad de la temperatura es mayor en los dos periodos recientes en la estación 01018, en comparación con el periodo base 1961-1990, y se ha reducido en el periodo 1981-2010 en la estación 09043, al estar más acotada la curva en torno a su media, con respecto a la curva del periodo base.

Discusión y conclusiones

El presente estudio muestra que las tendencias de cambio en temperatura a largo plazo definen de manera objetiva que está presente en dicha variable en nuestro país. Reafirmando la influencia de continentalidad en las estaciones climáticas analizadas, se puede asegurar que independientemente de que la influencia sea mayormente marina o continental, la señal de calentamiento global está presente.

Las 1 099 estaciones, consideradas a partir de diversos conjuntos de estaciones utilizadas en estudios climáticos, conforman una base de datos climática depurada que se puede tomar como referencia, donde sólo quedaría realizar el complemento del análisis de calidad a través de la homogeneidad, proceso pendiente que representa un gran reto debido a la falta del metadato correspondiente para cada sitio. En este sentido, es recomendable que se realice en México una “antropología de metadatos”, para documentar in situ y con registros históricos los cambios en las estaciones que permitan sistematizar los metadatos.

Con este estudio también se muestra que a pesar de contar con un número mayor a 5 000 estaciones climáticas en el país, para estudios específicos de tendencias de temperatura este número se reduce sustancialmente a poco más de 100 estaciones irregularmente distribuidas en México, situación que conlleva a retomar el tema de la confiabilidad de la información climática y su distribución espacio-temporal.

De la distribución de 112 estaciones seleccionadas objetivamente se identifican regiones con baja cobertura en zonas del norte y sur del país, principalmente. Esta situación representa una ventana de oportunidad para fortalecer las capacidades institucionales, la infraestructura de la red de medición, el control de calidad y homogeneidad, la sistematización y difusión de la información meteorológica, climática e hidrométrica a lo largo del territorio, específicamente en las regiones que se identifican con escasa cobertura. Un beneficio directo de las mejoras será contar con insumos para una mejor gestión del riesgo de desastres ante fenómenos hidrometeorológicos y climáticos extremos, así como para la adaptación al cambio climático.

También es necesario considerar que este estudio muestra objetivamente tendencias tanto positivas como negativas, siendo las primeras mayores en una tasa de 3/1, denotando que la tendencia positiva predomina; además, bajo este análisis se confirman las declaraciones del IPCC en el sentido de que las tendencias pueden llegar a ser diferentes, dependiendo de la latitud, longitud, topografía (^{Solomon et al., 2007}) y ecoentorno, entre otros. En el caso de México, la interacción humana cercana a los sitios de las estaciones climáticas también es un factor por considerar. Otro factor en el que se debe realizar análisis se relaciona con los errores instrumentales en la medición, pues existe el riesgo de que las observaciones puedan parecer y comportarse como “buenos” datos, seguir rigurosamente teorías físicas válidas, y aun así conducir a conclusiones erróneas.

Asimismo, algunas estaciones con tendencia positiva y negativa están cercanas e intercaladas, lo cual tal vez esté relacionado con el entorno urbano/rural de las estaciones, diferencias en altitud, procesos físicos locales del clima y mantenimiento a la infraestructura de medición.

A partir del análisis de tendencias locales se identificó que el periodo más reciente de calentamiento intenso inicia generalmente a mediado de la década de 1970, siendo congruente con la señal identificada por el ^{IPCC (2013)} para la temperatura media global. En las zonas rurales se observa un desfase de hasta 10 años en el inicio del periodo reciente de calentamiento, en comparación con zonas urbanas para una misma región de México. El calentamiento y la diferencia entre periodos cálidos y fríos son mayores, regularmente, en áreas urbanas, en comparación con las rurales. El mayor cambio de la temperatura, con valores entre 1 y 3 °C, se presenta en el noroeste y en el norte, seguido por estaciones en el centro del país y la costa norte del Pacífico, con rangos de calentamiento entre 1 y 2 °C, y con los menores cambios, menor a 1 °C, en la costa del Golfo de México, con excepción de la ubicada en la ciudad de Tampico, Tamaulipas.

Queda pendiente el tema de la atribución. Si bien se proponen hipótesis de algunas de sus posibles causas, es necesario realizar estudios de investigación científica, para poder determinar ya sea en forma fraccional o completa algunos de los forzantes tanto naturales como inducidos de forma completamente objetiva.