El radar Popocatépetl está instalado dentro del Parque Nacional Iztaccíhuatl-Popocatépetl, a 4000 m s.n.m. en el cerro Altzomoni y a 11 km de distancia del cráter. Su función principal es monitorear las cenizas del volcán. El proyecto fue llevado a cabo con financiamiento de la Secretaría de Educación, Ciencia, Tecnología e Innovación del Gobierno de la Ciudad de México, mediante el convenio SECITI/90/2017. Para la instalación se contó con el apoyo técnico de Ingeniería Geofísica y Sistemas S.A. de C.V.

Integrantes del proyecto

Dr. Hugo Delgado Granados

Líder del proyecto. Experto en vulcanología.

Instituto de Geofísica, UNAM

Dr. José Agustín García Reynoso

Experto en contaminación ambiental.

Instituto de Ciencias de la Atmósfera y Cambio Climático, UNAM

Fís. Eric B. Téllez

Técnico encargado de la operación.

Instituto de Geofísica, UNAM

Ing. Gerardo Agustín Galguera Rosas

Enlace administrativo SECTEI-UNAM.

Dr. Adolfo V. Magaldi Hermosillo

Experto en radares.

Escuela Nacional de Estudios Superiores Unidad Juriquilla, UNAM

Dr. Jorge Zavala Hidalgo

Experto en radares.

Instituto de Ciencias de la Atmósfera y Cambio Climático, UNAM

M. en C. Dulce Rosario Herrera Moro

Visibilidad web de productos.

Instituto de Ciencias de la Atmósfera y Cambio Climático, UNAM

Ing. Miguel Ángel Vela Rosas

Desarrollador de página web.

Instituto de Geofísica, UNAM

¿Cómo opera?

El radar opera en cinco formas diferentes:

¿Qué nos permite obtener o visualizar?

Productos primarios

Hay ocho productos primarios provenientes del radar:

  • 1. Reflectividad horizontal

    $$ \mathbf{Z_{hh} \; [dBz] }$$

    Definición

    El factor de reflectividad es el parámetro más importante para interpretar los datos del radar. Se obtiene a partir de un modelo de dispersión de Rayleigh y se define teóricamente como la suma de los diámetros de las partículas (gotas) contenidas en el volumen de muestreo, elevados a la sexta potencia:

    $$ Z = \sum_{vol} N (D) D^6 $$

    donde $ Z $ se expresa se expresa en milímetros elevados a la sexta potencia entre metro cúbico. En muchos casos, el número, la composición y la forma de las partículas no se conocen y se define una reflectividad equivalente o efectiva $ Z_e $. Por ejemplo, para las partículas de nieve y de hielo hay que tomar una reflectividad equivalente $ Z_e $ que representa a $ Z $, suponiendo que todas las partículas retrodispersoras son gotas esféricas de densidad $ \rho $.

    Con el fin de abarcar el intervalo de valores, con frecuencia se trabaja en una escala logarítmica o con unidades de $ dBZ $ que se definen numéricamente como $ dBZ_e = 10log_{10}Z_e$.

    En el caso del radar polarimétrico, la cantidad $ Z_{hh} $ se refiere únicamente a la reflectividad equivalente en la dirección horizontal. Mientras que la cantidad $ Z_{vv} $ se refiere a la reflectividad en la dirección vertical.

    La escala empleada en las imágenes utiliza colores fríos (azules a verdes) para indicar zonas con intensidades débiles y colores cálidos (verdes a rojos) para zonas con intensidades fuertes.

    Uso

    Visualizar los objetivos dentro del volumen de análisis del radar.

  • 2. Reflectividad diferencial

    $$ \mathbf{ZDR \; [dB] }$$

    Definición

    La reflectividad diferencial es la razón logarítmica de la reflectividad polarizada horizontalmente con la reflectividad polarizada verticalmente. Valores positivos de $ Z_{DR} $ indican que los hidrometeoros predominantes son más grandes en la dirección horizontal que en la vertical. Valores negativos indican que indican que los hidrometeoros predominantes son más grandes en la dirección vertical que en la horizontal. Valores cercanos a cero indican que los hidrometeoros son similares en tamaño en la dirección horizontal y en la vertical.

    $$ Z_{DR} = 10\log_{10}(\frac{Z_{hh}}{Z_{vv}}) $$

    Uso

    Es un buen indicador de la forma de los objetos. Los objetos esféricos tendrán un $ Z_{DR} $ casi igual a cero, mientras que los objetos elipsoidales tendrán un valor positivo o negativo dependiendo de si son más grandes en la dirección horizontal o en la vertical respectivamente.

  • 3. Fase diferencial de la polarización cruzada

    $$ \mathbf{ \phi_{DP} } $$

    Definición

    Es la diferencia del cambio de fase entre las ondas polarizadas verticales y horizontales.

    $$ \phi_{DP} = \phi_{hh} - \phi_{vv} $$

    En el que $ \phi_{hh} $ y $ \phi_{vv} $ son el cambio de fase diferencial acumulado horizontal y vertical, respectivamente, en el trayecto de ida y vuelta entre el radar y el volumen de estudio.

    Uso

    La forma de los objetos afectan a $ \phi_{DP} $ igual que a $ Z_{DR} $. Los objetos esféricos tendrán una $ \phi_{DP} $ casi igual a cero, mientras que los objetos elipsoidales tendrán un valor positivo o negativo dependiendo de si son más grandes en la dirección horizontal o en la vertical respectivamente. Además, $ \phi_{DP} $ depende de la concentración de las partículas: mientras más partículas haya en el volumen del pulso, más grande será $ \phi_{DP} $.

  • 4. Fase diferencial específica

    $$ \mathbf{ K_{DP} } $$

    Definición

    Es el gradiente de la fase diferencial con respecto a la distancia.

    $$ K_{DP} = \frac{\phi_{DP}(r_2)-\phi_{DP}(r_1)}{2(r_2 - r_1)} $$

    Uso

    Para caracterizar la cantidad de lluvia en el volumen de datos del radar, ya que el aumento en $ K_{DP} $ indica el aumento en el tamaño y concentración de gotas de lluvia y, por lo tanto, aumento en la tasa de precipitación.

  • 5. Coeficiente de correlación copolar

    $$ \mathbf{ \rho_hv } $$

    Definición

    El coeficiente de correlación copolar $ \rho_{hv} $ o rho es una correlación estadística entre los ecos polarizados vertical y horizontalmente. Proporciona una medida de la consistencia de las formas y tamaños de los objetivos dentro del haz del radar. Un valor más alto indica una consistencia más alta en el tamaño y la forma de los objetivos, mientras que un valor menor indica una mayor variabilidad en formas y tamaños.

    $$ \rho_{hv} = \frac{\lt S_{vv}S^*_{hh}\gt}{\sqrt{\lt | S_{hh} |^2 \gt \cdot \lt | S_{vv} |^2 \gt }} $$

    donde $ S_{hh} $ y $ S_{vv} $ son las amplitudes dispersadas en la dirección horizontal y vertical respectivamente.

    Uso

    Se puede usar para ayudar a distinguir entre objetivos meteorológicos (agua, granizo) y no meteorológicos (aves, insectos, ceniza).

  • 6. Intensidad de la lluvia

    $$ \mathbf{ R \; [mm/h] } $$

    Definición

    Es una cantidad calculada a partir de la reflectividad $ (dBz) $ con la ecuación de Marshall-Palmer:

    $$ R[\frac{mm}{h}] = (\frac{10(\frac{dBz}{10})}{200})^\frac{5}{8} $$

    Uso

    Para estimar la cantidad de lluvia en $ mm/h $ que puede caer en la región de análisis.

  • 7. Velocidad Doppler

    $$ \mathbf{ V \; [m/s] } $$

    Definición

    Se calcula con la diferencia entre la frecuencia emitida y recibida por el radar debido a que la señal impacta con un objetivo móvil. Un número positivo indica que se está alejando del radar y un número positivo indica que va hacia el radar.

    Uso

    Para medir la velocidad de objetos dentro del volumen de observación.

  • 8. Ancho del espectro de la velocidad Doppler

    $$ \mathbf{ W \; [m/s] } $$

    Definición

    Es la desviación estándar de la velocidad Doppler. Se extiende por perturbaciones y ruido.

    Uso

    Para eliminar el ruido de la velocidad Doppler.

Productos secundarios

En el Instituto de Geofísica se desarrolló una presentación de los datos del radar a partir de los productos primarios, adaptando la biblioteca de software Wradlib:


Vol-CAPPI

Otro método para visualizar los datos se llama Indicador de Posición en Planta de Altitud Constante en Volumen. Esta visualización permite ver un eco en tres secciones transversales. Además de la vista de planta usual (al centro), se muestran dos elevaciones: una correspondiente a una elevación lateral de este a oeste (sección transversal horizontal) y otra correspondiente a una elevación lateral de norte a sur (sección transversal vertical). En general, este método se combina con los datos de reflectividad.

Información técnica

  • Radar Doppler meteorológico compacto de polarización doble banda X marca Furuno modelo WR-2100.
  • Alcance máximo: 70 km de radio.
  • Frecuencia de operación: 9.4 GHz
  • Ancho del haz: 2.7 grados (vertical y horizontal).
  • Velocidad máxima de rotación: 16 rpm.
  • Medición Doppler máxima: ±48 m/s.
  • Rango angular de escaneo vertical: -2 a 182 grados.
  • Rango angular de escaneo horizontal: 360 grados.
  • Resolución angular: 0.1 grados.
  • Precisión angular: 0.2 grados.