An anisotropic meandering model for atmospheric dispersion in low-wind-speed conditions

Resumen Técnico: Un Modelo de Meandro Anisotrópico para la Dispersión Atmosférica en Condiciones de Viento Débil

Planteamiento del Problema
La dispersión de contaminantes bajo condiciones de viento débil (típicamente $< 1.5\text{--}2 \text{ m s}^{-1}$) presenta un desafío crítico en el modelado de la calidad del aire. En estos regímenes, la dirección media del viento está mal definida debido a oscilaciones horizontales de gran escala y baja frecuencia conocidas como "meandro" (meandering). Estas oscilaciones dominan el transporte de contaminantes, extendiendo las plumas sobre amplios sectores angulares y, a menudo, provocando altas concentraciones a nivel del suelo debido a la débil advección. Los modelos operativos existentes suelen fallar en estas condiciones porque, o bien descuidan la difusión longitudinal (que se vuelve significativa en números de Péclet bajos), o bien absorben el meandro en una difusividad lateral efectiva incrementada, un procedimiento que distorsiona la distribución angular de la concentración. Además, la naturaleza mesoescala del meandro dificulta su parametrización utilizando únicamente variables locales.

Metodología
Los autores desarrollan un modelo analítico que trata el núcleo difusivo y el fenómeno de meandro como dos componentes distintos y analíticamente acoplados:

1. Núcleo Difusivo: El modelo utiliza la solución cerrada exacta de la ecuación de advección-difusión tridimensional con difusividades de remolino constantes pero distintas (anisotrópicas) ($K_x, K_y, K_z$). Crucialmente, se retiene el término de difusión longitudinal ($K_x \partial^2 C / \partial x^2$), ya que no puede ser despreciado en vientos débiles. La solución se deriva mediante isotropización de coordenadas y transformación en una ecuación de Helmholtz modificada, produciendo una función de Green que tiene en cuenta la reflexión en el suelo mediante el método de imágenes.
2. Representación del Meandro: En lugar de modificar la difusividad, el meandro se representa explícitamente mediante una convolución angular. La concentración promediada en el tiempo se calcula como la superposición de plumas instantáneas (el núcleo difusivo) rotadas por una dirección del viento $\theta$, ponderadas por una función de densidad de probabilidad gaussiana de las direcciones del viento. Este enfoque trata la dirección del viento como una variable aleatoria de variación lenta.
3. Parametrización:
   • Difusividades Horizontales ($K_x, K_y$): Derivadas de la teoría de la difusión estadística de Taylor aplicada al espectro de velocidad de viento bajo propuesto por Mortarini y Anfossi. Este espectro incluye explícitamente el "pico de meandro" (la transformada de Fourier de la función de autocorrelación de Frenkiel).
   • Difusividad Vertical ($K_z$): Parametrizada según el régimen de estabilidad de la capa límite planetaria (estable o convectiva). Para el modelo de coeficiente constante, $K_z$ se evalúa en la altura de la línea central de la pluma.
   • Ancho Angular ($\sigma_\theta$): Calculado a partir de la desviación estándar de la velocidad lateral ($\sigma_v$) y la velocidad media del viento ($U$), incorporando un límite inferior de viento bajo ($\sigma_v \ge 0.5 \text{ m s}^{-1}$) para asegurar un nivel mínimo de meandro.

Contribuciones Clave

• Separación de Mecanismos: El modelo separa analíticamente el transporte difusivo (regido por difusividades de remolino anisotrópicas y difusión longitudinal) del transporte por meandro (regido por la estadística de la dirección del viento), evitando la distorsión de los perfiles angulares de concentración comunes en los modelos que utilizan una única difusividad lateral incrementada.
• Solución Exacta en Forma Cerrada: El modelo proporciona una solución exacta a la ecuación de advección-difusión 3D con difusión longitudinal, evitando la necesidad de inversiones de Laplace o series infinitas, manteniéndose consistente con la pluma gaussiana clásica en límites de viento alto.
• Validación Frente a Diversos Regímenes: El modelo se valida contra datos experimentales del experimento INEL (capa límite estable) y el experimento IIT Delhi (capa límite convectiva), demostrando un desempeño satisfactorio en ambos regímenes de estabilidad.
• Comparación con Modelos Multicapa: Una contribución significativa es la comparación con el modelo GIADMT (Generalized Integral Advection–Diffusion Multilayer Technique), un modelo semi-analítico que retiene perfiles de difusividad vertical dependientes de la altura. El estudio aísla el efecto del tratamiento de la difusividad vertical alimentando a ambos modelos con las mismas difusividades horizontales y distribuciones de meandro.

Resultados

• Desempeño Estadístico: El modelo propuesto cumple con los criterios de aceptación estándar para índices estadísticos (Error Cuadrático Medio Normalizado, Sesgo Fraccional, Correlación, Factor de Dos) tanto para los conjuntos de datos de INEL como de IIT Delhi. Logra un Factor de Dos (FA2) de 0.87 para el caso estable de INEL y de 0.88 para el caso convectivo de IIT Delhi.
• Comparación con GIADMT: Al compararse con el modelo GIADMT (que utiliza un perfil de $K_z(z)$ dependiente de la altura), el modelo de coeficiente constante propuesto funcionó significativamente mejor:
  • Concentraciones Máximas: El GIADMT sobreestimó sistemáticamente las concentraciones máximas a nivel del suelo, fallando en los criterios de aceptación, mientras que el modelo propuesto se mantuvo bien calibrado.
  • Distribución Angular: Para la distribución angular completa a lo largo de los arcos de muestreo, el modelo propuesto redujo el Error Cuadrático Medio (RMSE) en un factor medio de aproximadamente 3.4 en comparación con GIADMT.
• Mecanismo de Mejora: El análisis revela que, en condiciones de viento bajo, la concentración a nivel del suelo está controlada por la difusividad vertical en la delgada capa entre la fuente y el receptor. El perfil dependiente de la altura utilizado en GIADMT desaparece en el suelo, proporcionando una difusividad vertical insuficiente en esta capa poco profunda, lo que suprime la dilución vertical e infla las concentraciones predichas. Evaluar $K_z$ en la línea central de la pluma (como hace el modelo propuesto) proporciona un valor mejor calibrado para esta capa específica.

Significancia y Reivindicaciones
El artículo sostiene que la complejidad matemática de resolver la ecuación de advección-difusión con coeficientes dependientes de la altura no necesariamente mejora la precisión predictiva en condiciones de viento bajo. El estudio demuestra que:

1. En vientos bajos, la dispersión lateral está dominada por el meandro, el cual es capturado igualmente bien por la convolución angular en ambos modelos.
2. La precisión de las predicciones de concentración a nivel del suelo está controlada primordialmente por la difusividad vertical en la capa fuente-receptor.
3. Una solución exacta simple de forma cerrada con una difusividad vertical constante (evaluada en la línea central de la pluma) describe la dispersión en condiciones de viento bajo tan bien como, o incluso mejor que, esquemas más sofisticados que retienen perfiles dependientes de la altura.

Los autores concluyen que, siempre que el meandro y la anisotropía de la difusión estén correctamente parametrizados, una solución exacta simple es suficiente para modelar la dispersión en vientos bajos, desafiando la necesidad de esquemas numéricos más complejos para este régimen específico. Se pretende aplicar el modelo a otros experimentos en trabajos futuros, abordar la variabilidad mesoescala en la varianza lateral y combinar la convolución angular con modelos de dispersión fraccional.