Los modelos de Markov: Un acercamiento a la dinámica reciente de los huracanes en México

Por: Pedro Andrés Sánchez Gutiérrez

INTRODUCCIÓN

Los modelos de Markov son una conceptualización del camino al equilibrio que las cosas o ciertos entes transitan (Page, 2012). Estos modelos tienen amplia aplicación en evaluación de la confiabilidad de sistemas instrumentados de seguridad (Guo & Yang, 2008), y la determinación de niveles de integridad de seguridad en la industria de procesos (Knegtering & Brombacher, 1999).  Los modelos de Markov consisten en dos aspectos principalmente:

  • Estados. Los entes se mueven entre dichos estados.
  • Probabilidades de transición. Son la probabilidad de un ente para moverse de un estado a otro.

En el presente estudio se analiza la evolución, de acuerdo con la máxima categoría alcanzada, de los huracanes que han tocado territorio mexicano, del año 1989 al año 2018; considerando tres “estados”: las categorías -en la escala Saffir-Simpson- H1 y H2 como el estado de menor magnitud, la categoría H3 como estado intermedio, y las categorías H4 y H5 juntas como el estado de mayor magnitud.

No se trata de un estudio cuantitativo en cuanto a la incidencia del número de huracanes en México; es más bien una forma de observar el cambio, en el tiempo, en la distribución de los huracanes que año con año afectan al país en relación con su fuerza o magnitud. En este caso, huracanes que impactaron territorio mexicano provenientes tanto del Océano Pacífico como del Atlántico.  

CONSIDERACIONES GENERALES

En general, los procesos naturales y artificiales pueden tender al equilibrio, ser cíclicos, aleatorios o complejos. Bajo la suposición de que las emisiones humanas de gases efecto invernadero (GEI) a la atmósfera llegarán a un límite1, se puede considerar que la dinámica de la distribución de las magnitudes de los huracanes seguirá un comportamiento convergente al equilibrio.

De esta manera, dicha dinámica, puede describirse mediante modelos de Markov, los cuales operan bajo las siguientes condiciones (Page, 2012):

  1. Un número finito de estados.
  2. Las probabilidades de transición deben ser fijas. De un periodo a otro, la probabilidad de un ente de moverse de un estado a otro permanece constante.
  3. Los entes pueden eventualmente transitar de un estado a cualquier otro.
  4. No se trata de un ciclo simple.

Con base en la información disponible oficial (CONAGUA, 2021), para el periodo 1989–2018, se consideró agrupar los eventos -huracanes- en intervalos de 5 años contabilizando los eventos en cada uno de ellos y agrupándolos en los 3 estados mencionados -bajo la lógica de modelos de Markov- debido a su intensidad en la escala Saffir-Simpson.

Tabla 1. Escala de huracanes Saffir-Simpson (Comisión Nacional del Agua, 2002).

CategoríaVelocidad de los vientos (km/h)
H1119-153
H2154-177
H3178-209
H4210-249
H5mayor de 250

RESULTADOS

Tras la agrupación mencionada, en la Figura 1 puede observarse la dinámica de la intensidad de los huracanes que han impactado en México durante este periodo. Puede advertirse el efecto del cambio climático al no presentarse ningún huracán de categoría H4 o H5 al inicio del periodo de análisis hasta representar el 46% del total de dichos eventos en el último quinquenio del periodo de análisis. Asimismo, si se considera el número de eventos que tocaron territorio mexicano en el primer quinquenio del periodo de análisis -9 eventos- vs. el número de huracanes que se presentaron en el último quinquenio considerado -13 eventos- el aumento en el número de eventos fue del orden del 44%. Es decir, se presentan más eventos y más potentes cada vez.

Lo anterior, habla de una energización de los fenómenos hidrometeorológicos, especialmente los que emanan en buena medida de las condiciones termodinámicas del océano (Sánchez-Gutiérrez, 2020).

Aquí, cabe señalar que los océanos acumulan la mayor parte de la energía que queda atrapada en la atmósfera por acción de los GEI. Tan sólo en el periodo 1971-2010, se estima que los océanos en el planeta absorbieron más del 90% de la energía calorífica extra acumulada resultante de los niveles crecientes de GEI en la atmósfera (Pachauri, 2014).

Otro dato destacable es que la mitad del incremento de la energía calorífica acumulada en los océanos ha ocurrido en los últimos 20 años (Cheng, y otros, 2020).

Con los estados de las categorías de la escala de Saffir-Simpson agrupadas en el periodo de análisis es posible construir una matriz de transición de Markov, la cual integra un conjunto de estados y proporciona la información sobre qué porcentaje de determinado estado en el tiempo t va a transitar a un determinado estado en el tiempo t+1. Ver Tabla 2.

Tabla 2. Matriz de Transición (Modelo de Markov) propuesta para la dinámica observada en la distribución de los huracanes en México, 1989-2018; con 3 estados que agrupan las categorías en la escala Saffir-Simpson.

[H1 – H2]t[H3]t[H4 – H5]t
[H1 – H2]t+10.90.10.03
[H3]t+10.020.750.07
[H4 – H5]t+10.080.150.9

Bajo una aproximación frecuentista, el 90% de los huracanes permanecerá en categorías 1 y 2 en el tiempo t+1, 2% evolucionará a categoría 3 y 8% se convertirá en huracanes de categoría 4 y 5; por otra parte, del subconjunto de huracanes en el estado intermedio correspondiente a la categoría 3, 10% transitará hacia el estado de las intensidades 1 y 2 en el tiempo t+1, 75% permanecerá en la categoría 3 y 15% se convertirá en huracanes de intensidad 4 ó 5. Finalmente, de los huracanes pertenecientes al estado que agrupa las categorías 4 y 5, sólo 3% de ellos se convertirá en huracanes de intensidad 1 ó 2 en el tiempo t+1, 7% se convertirá en huracanes de categoría 3, y 90% permanecerá en el estado que agrupa a las categorías 4 y 5 de la escala Saffir-Simpson.  Donde t y t+1 son quinquenios para el presente caso.

De esta manera, aplicando las probabilidades de la matriz de transición a la distribución de huracanes correspondientes al quinquenio 1989-1993, se obtiene la dinámica mostrada en la Figura 2. En ella, los valores de las categorías se estabilizan, tienden a un límite o asíntota; es decir, se aproximan al equilibrio.

Figura 2. Distribución de las magnitudes de los huracanes en México bajo un modelo de Markov quinquenalmente, 1989-2018, que representa la convergencia al equilibrio de los 3 estados que agrupan las categorías H1 a H5 (escala Saffir-Simpson) bajo un modelo de Markov propuesto para el caso de los huracanes que impactaron territorio mexicano en años recientes.

CONCLUSIONES

Se analizaron las magnitudes de los huracanes que impactaron territorio mexicano en el periodo 1989–2018. Por razones de simplificación el periodo se dividió en 6 quinquenios y las 5 categorías de huracanes en la escala Saffir-Simpson se agruparon en 3 grupos, considerados estados en la lógica de los modelos de Markov.

Asumiendo que las magnitudes de los huracanes eventualmente alcanzarán el equilibrio en el tiempo –bajo el supuesto verosímil de que las emisiones de GEI producto de las actividades humanas se estabilizarán en un futuro próximo- se llega a una aproximación relativamente satisfactoria de la dinámica real observada. Lo anterior, mediante la construcción de una matriz de transición que precisamente representa una narrativa del comportamiento de la distribución de las magnitudes de los huracanes en el tiempo, durante el periodo de análisis.

Un análisis más minucioso habrá que realizarse respecto a la pertinencia de esta metodología tomando en cuenta la climatología, que considera oscilaciones tanto oceánicas como atmosféricas. Por lo pronto, se propone como un punto de partida para el análisis de este tipo de fenómenos hidrometeorológicos en México a la luz de las causas y efectos del cambio climático, así como de las medidas de la mitigación de emisiones de GEI, las cuales muestran una tendencia a la mejora, sobre todo en los principales países emisores de GEI.

REFERENCIAS

Buckley , T. (27 de Febrero de 2019). Institute for Energy Economics and Financial Analysis. Obtenido de Over 100 Global Financial Institutions Are Exiting Coal, With More to Come | Every Two Weeks a Bank, Insurer or Lender Announces New Restrictions on Coal: https://ieefa.org/wp-content/uploads/2019/02/IEEFA-Report_100-and-counting_Coal-Exit_Feb-2019.pdf

Chandler, D. L. (20 de Noviembre de 2018). Massachusetts Institute of Technology. Obtenido de MIT News on Campus and Around the World: http://news.mit.edu/2018/explaining-dropping-solar-cost-1120

Cheng, L., Abraham, J., Zhu, J., Trenberth, K. E., Fasullo, J., Boyer, T., . . . Mann, M. E. (2020). Record-setting ocean warmth continued in 2019. Adv. Atmos. Sci., 37(2), 137-142. Obtenido de https://doi.org/10.1007/s00376-020-9283-7

Comisión Nacional del Agua. (2002). Compendio Básico del Agua en México 2002. Subdirección General de Programación. México, D. F.: CONAGUA.

Comisión Nacional del Agua. (2012). Análisis de las temporadas de huracanes de los años 2009, 2010 y 2011 en México. México, D. F.: SEMARNAT.

CONAGUA. (15 de Agosto de 2021). Sistema Nacional de Información del Agua. Obtenido de https://www.gob.mx/conagua/acciones-y-programas/sistema-nacional-de-informacion-del-agua-sina

Di Paola, A. (6 de Febrero de 2020). Bloomberg. Obtenido de Markets | Economics: https://www.bloomberg.com/news/articles/2020-02-06/middle-east-s-2-trillion-wealth-could-be-gone-by-2034-imf-says

Guo, H., & Yang, X. (2008). Automatic creation of Markov models for reliability assessment of safety instrumented systems. Reliability Engineering and Safety System, 93, 807-815.

Harrington, R. (29 de Septiembre de 2015). Business Insider. Obtenido de Home | Tech Contributors: https://www.businessinsider.com/this-is-the-potential-of-solar-power-2015-9

Houser, T., & Pitt, H. (7 de Enero de 2020). Rhodium Group. Obtenido de Preliminary US Emissions Estimates for 2019: https://rhg.com/research/preliminary-us-emissions-2019/

IRENA. (2020). International Renewable Energy Agency. Obtenido de Renewable Energy and Jobs – Annual Review 2020: https://www.irena.org/publications/2020/Sep/Renewable-Energy-and-Jobs-Annual-Review-2020

Knegtering, B., & Brombacher, A. C. (1999). Application of micro Markov models for quantitative safety assessment to determine safety integrity levels as defined by the IEC 61508 standard for functional safety. Reliability Engineering and System Safety, 66, 171-175.

McKerracher, C. (15 de Abril de 2021). BloombergNEF. Obtenido de EV Charging Data Shows Widely Divergent Global Path: https://about.bnef.com/blog/ev-charging-data-shows-widely-divergent-global-path/

Myllyvirta, L. (1 de Marzo de 2021). CarbonBrief. Obtenido de China Policy: https://www.carbonbrief.org/analysis-chinas-co2-emissions-surged-4-in-second-half-of-2020/

Pachauri, R. K. (2014). Climate Change 2014: Synthesis Report. Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press. Obtenido de https://archive.ipcc.ch/pdf/assessment-report/ar5/syr/SYR_AR5_FINAL_full_wcover.pdf

Page, S. E. (2012). Markov Models. Model Thinking. NY: Coursera.

Sánchez-Gutiérrez, P. A. (30 de Abril de 2020). Blog de la Cátedra UNESCO-UDLAP en Riesgos Hidrometeorológicos. Obtenido de Datos energéticos y de salud de los océanos Pt. 1: https://unesco.udlap.xyz//datos-energeticos-y-de-salud-de-los-oceanos-pt-i/

Solar.com. (Abril de 2021). Solar Learning Center. Obtenido de Solar Panel Efficiency – Pick the Most Efficient Solar Panels: https://www.solar.com/learn/solar-panel-efficiency/

Varadhan, S. (25 de Junio de 2019). Reuters. Obtenido de India plans to add 500 GW renewable energy by 2030: government: https://www.reuters.com/article/us-india-renewables/india-plans-to-add-500-gw-renewable-energy-by-2030-government-idUSKCN1TQ1R9

Vogt-Schilb, A. (28 de Abril de 2021). Inter-American Development Bank. Obtenido de The path to net zero emissions can create millions of jobs in Latin America: https://blogs.iadb.org/sostenibilidad/en/the-path-to-net-zero-emissions-can-create-millions-of-jobs-in-latin-america/

Whiteman, A., Akande, D., Elhassan, N., Escamilla, G., Lebedys, A., & Arkhipova, I. (2021). International Renewable Energy Agency. Obtenido de Renewable Capacity Statistics 2021: https://www.irena.org/-/media/Files/IRENA/Agency/Publication/2021/Apr/IRENA_RE_Capacity_Statistics_2021.pdf