Site Quality Analysis for an Indian Submillimeter Telescope: A Reanalysis-Based Approach

Resumen Técnico: Evaluación de Sitios para un Observatorio Submilimétrico en Ladakh

1. Planteamiento del Problema

La astronomía submilimétrica (ondas milimétricas y submilimétricas) es fundamental para estudiar procesos astrofísicos en entornos fríos y difusos, como la formación de estrellas y planetas, la química del medio interestelar y el universo de alto desplazamiento al rojo. Sin embargo, la observación desde tierra está severamente limitada por la absorción atmosférica, principalmente causada por el vapor de agua.

La métrica clave para la viabilidad de un sitio es el Vapor de Agua Precipitable (PWV, por sus siglas en inglés). Para acceder a ventanas de frecuencia altas y mantener una baja temperatura de brillo del cielo (ruido), se requiere un PWV $\lesssim 1$ mm. Aunque existen observatorios de clase mundial en el hemisferio sur (como ALMA en Chile y sitios en la Antártida), la falta de un sitio seco en el hemisferio norte limita la cobertura de objetivos del norte. La región de Ladakh, en la meseta del Himalaya indio, ofrece un entorno desértico frío, alta altitud y baja actividad antropogénica, pero carecía de un análisis sistemático de su idoneidad para telescopios submilimétricos de clase científica.

2. Metodología

El estudio emplea un enfoque basado en datos de reanálisis atmosférico para realizar una exploración preliminar de gran escala antes de realizar mediciones in situ.

• Fuente de Datos: Se utilizó el conjunto de datos de reanálisis ERA5 (quinta generación del ECMWF), que cubre el período de enero de 2010 a abril de 2025 (184 meses). ERA5 ofrece una resolución espacial de $0.25^\circ \times 0.25^\circ$ (aprox. 23 km $\times$ 27 km en Ladakh) y una resolución temporal mensual.
• Estimación del PWV: Se compararon tres métodos para estimar el PWV a partir de los datos de ERA5:
  1. M1 (TCWV nativo): Uso directo del valor de vapor de agua total de la columna proporcionado por ERA5.
  2. M2 (Específico del sitio): Interpolación de la altura geopotencial para determinar la presión superficial exacta del sitio ($p_s$) y luego integrar el perfil de humedad específica desde la cima de la atmósfera hasta $p_s$. Este método corrige los errores de representación topográfica.
  3. M3 (Rejilla nativa): Uso de la presión superficial promedio de la celda de la rejilla ERA5 como límite inferior de integración.
  • Validación: Los métodos se validaron comparándolos con mediciones de radiómetros de inclinación (220 GHz) en Hanle (India), Mauna Kea (Hawái) y Chajnantor (Chile). Se determinó que el Método 2 (M2) es el más preciso para sitios específicos, mientras que el Método 3 (M3) es adecuado para el escaneo regional.
• Modelado de Transferencia Radiativa: Se utilizaron los valores de PWV derivados para calcular la transmitancia atmosférica y la temperatura de brillo del cielo en el rango de 10–1000 GHz utilizando el código de transferencia radiativa am (desarrollado por el Observatorio Astrofísico Smithsonian). Se evaluaron estados atmosféricos medianos (todos los meses) y óptimos (meses de invierno: noviembre-febrero).

3. Contribuciones Clave

• Validación de Metodología: Demostración de que el uso de la presión superficial específica del sitio (M2) es crucial para evitar sobreestimaciones masivas del PWV en terrenos complejos y de alta altitud, donde la presión de la celda de la rejilla nativa no representa la elevación real del observatorio.
• Mapeo Regional de Ladakh: Primera evaluación sistemática de la variabilidad espacial del PWV en toda la meseta de Ladakh, identificando regiones que superan a los sitios de observación existentes (Hanle y Merak).
• Identificación de Nuevos Sitios Candidatos: Detección de dos áreas específicas (Sitio A y Sitio B) que ofrecen condiciones atmosféricas significativamente mejores que los sitios actuales en la región.

4. Resultados Principales

• Comparación Global: El sitio Dome A (Antártida) mantiene el mejor PWV (< 1 mm todo el año). En comparación, Hanle y Merak tienen PWV medianos más altos, con Hanle mostrando un mínimo mensual de ~1.05 mm.
• Identificación de Sitios Prometedores en Ladakh:
  • Sitio A ($\approx 34.25^\circ$N, $78.75^\circ$E) y Sitio B ($\approx 32.50^\circ$N, $79.00^\circ$E) emergen como los candidatos óptimos.
  • Frecuencia de PWV $\le$ 1 mm:
    • Sitio A: Cumple el criterio en el 23% de los meses del estudio.
    • Sitio B: Cumple el criterio en el 19% de los meses.
    • Hanle (celda de rejilla): Solo el 5%.
    • Merak (celda de rejilla): Solo el 8%.
  • Estacionalidad: Durante los meses de invierno (NDJF), el rendimiento mejora drásticamente. El Sitio A mantiene PWV $\le$ 1 mm aproximadamente el 70% del tiempo en invierno, frente al ~15-20% de Hanle y Merak.
• Rendimiento en Frecuencias Submilimétricas:
  • Se calculó la frecuencia límite ($\nu_{0.8}$) donde la transmitancia es $\ge$ 0.8.
  • En condiciones de invierno (mediana), el Sitio A alcanza una frecuencia límite de ~351 GHz, permitiendo el uso viable de la ventana de 325–345 GHz.
  • En comparación, Hanle y Merak se limitan a ~313 GHz en invierno.
  • En la ventana de 850 $\mu$m, el Sitio A y B alcanzan transmitancias de 0.80 y ~0.79 respectivamente, superando significativamente a Hanle (0.71) y Merak (~0.67).
• Análisis Topográfico: El análisis de elevación (DEM) muestra que dentro de las celdas de la rejilla ERA5, la variación de elevación puede ser de 300-600 m. Esto implica que existen micro-sitios dentro de estas celdas que podrían ser aún más secos que el promedio de la celda.

5. Significado y Conclusiones

El estudio concluye que la región de Ladakh es favorable para la astronomía submilimétrica, ofreciendo condiciones que superan a los sitios de observación actuales en la India (Hanle y Merak) y son competitivos con otros sitios del hemisferio norte.

• Potencial Científico: La identificación del Sitio A y Sitio B sugiere que es posible construir un telescopio submilimétrico de clase científica en el hemisferio norte que pueda operar en ventanas de frecuencia más altas (hasta ~350 GHz en invierno) con menor ruido atmosférico.
• Siguientes Pasos: Los autores enfatizan que este análisis basado en reanálisis es necesario pero no suficiente. Se recomienda:
  1. Realizar mediciones in situ con radiómetros y GNSS en los sitios seleccionados durante al menos un ciclo anual completo para validar y corregir el sesgo positivo de humedad de ERA5.
  2. Evaluar factores logísticos, de infraestructura y de turbulencia atmosférica (seeing).
  3. Desarrollar una configuración de telescopio "esbozo" alineada con los ciclos de duty-cycle atmosférico inferidos.

En resumen, este trabajo proporciona la base científica y técnica para justificar la inversión en una campaña de medición de campo en Ladakh, con el objetivo final de establecer un nuevo observatorio submilimétrico que complemente la infraestructura global existente.