Un sobrevuelo de seis segundos identifica fuentes radiactivas mezcladas

Un helicóptero que cruza una zona contaminada a unos 10 m/s pasa solo unos segundos cerca de cada fuente. Un equipo del Instituto Paul Scherrer (PSI), en Suiza, y la Universidad de Míchigan afirma que eso basta. En un artículo publicado el 26 de octubre de 2026 en Nature Communications, emplean 6 s de datos de una sola pasada, contrastados con ensayos independientes in situ y de laboratorio. Los autores informan de una identificación concluyente de mezclas de fuentes (>5σ), una localización a escala métrica y una recuperación de la actividad a escala porcentual. El resumen de la universidad concreta más las cifras: fuentes localizadas con un margen de dos metros y actividades estimadas a menos de dos puntos porcentuales de su valor real.

El avance es sobre todo computacional. El detector es una matriz comercial de NaI(Tl). Lo que ha cambiado es cómo construye el equipo el modelo físico con el que se comparan los escasos recuentos, y la rapidez con que puede evaluarlo.

01Por qué es tan difícil interpretar el espectro de un sobrevuelo

Cada radionúclido emite rayos gamma a energías características. En principio, el espectro medido identifica el isótopo y su cantidad. En la práctica, cada material tiene una huella formada por su espectro de energías gamma, y esa huella se degrada en vuelo. Los fotones se dispersan en el aire y en la estructura de la aeronave, y pierden energía antes de llegar a los cristales. Además, los recuentos son escasos. Según el PSI, de unos 10 000 fotones emitidos por una fuente en el suelo, solo uno llega al helicóptero, y hasta ese uno se dispersa en las superficies metálicas del interior.

Los autores describen la inversión resultante como mal planteada. La espectrometría gamma móvil es el patrón de referencia para identificar fuentes, pero los métodos actuales son imprecisos porque los datos son escasos y la inversión está mal planteada. El remedio habitual es el análisis de espectro completo. En lugar de integrar unas pocas ventanas de energía, ajusta todo el espectro medido a plantillas simuladas. El problema es el coste. Unas plantillas fieles exigen simulaciones de transporte Monte Carlo a través del aire y del fuselaje, y para cada geometría de vuelo, atmósfera y posición de la fuente eso puede llevar horas o días. Nadie puede hacerlo durante una emergencia.

02Dividir el problema del transporte en dos

El truco central es la factorización. El equipo separa el recorrido de un fotón en el transporte atmosférico y la respuesta del helicóptero y del detector. Según el resumen de Míchigan, los cálculos atmosféricos tienen en cuenta la temperatura y la humedad del aire, mientras que la parte de la aeronave usa funciones precalculadas específicas del modelo de helicóptero. Los métodos del artículo describen un «banco de flujos de rayos gamma» generado con el módulo de desintegración radiactiva del código FLUKA. También describen simulaciones independientes de la plataforma que registran la deposición de energía en cada cristal y la reducen a espectros esperados de altura de pulso.

El grueso del trabajo Monte Carlo se hace, por tanto, una sola vez y por adelantado. El preprint indica unas 10² horas-núcleo por fuente, con 2 × 10⁹ partículas primarias repartidas en 10² ejecuciones independientes. Después, cada geometría de vuelo solo requiere una consulta y un paso de combinación. Breitenmoser sostiene que esta reorganización, y no una mayor potencia de cálculo, es lo que redujo a microsegundos unas simulaciones que antes tardaban horas o días.

El modelo de la plataforma se apoya en trabajos anteriores del PSI. El sistema suizo tiene cuatro cristales prismáticos de NaI(Tl) con un volumen total de 16,8 l, montados en la bodega de carga de un Super Puma. En un estudio de validación previo, el equipo notificó un error espectral relativo mediano inferior al 10 % para la mayoría de los radionúclidos en el modelo de respuesta del detector.

03Inferencia bayesiana sobre miles de plantillas

Con plantillas tan baratas, la inferencia puede permitirse ser exhaustiva. El marco compara miles de espectros predichos con el medido y actualiza de forma continua su confianza estadística sobre lo que hay debajo. Las palabras clave del preprint señalan el método de Monte Carlo basado en cadenas de Markov (MCMC) como muestreador. Las incógnitas se estiman conjuntamente: qué nucleidos hay, dónde está la fuente y cuál es su actividad. También se ajusta el fondo natural del potasio, el torio y el uranio del suelo.

El resultado es una distribución a posteriori, no un único mejor ajuste. Eso es importante desde el punto de vista operativo. El responsable de la intervención recibe una estimación de la actividad con su intervalo de credibilidad, y una preferencia cuantificada por «Cs-137 más Ba-133» frente a «solo Cs-137». El resultado de >5σ para la mezcla es una afirmación de este tipo.

04Qué se probó realmente

Los datos de campo proceden del ejercicio ARM22, realizado en el campo de entrenamiento militar de Thun el 16 de junio de 2022. Según el preprint, ambas configuraciones usaron las mismas dos fuentes selladas, de Cs-137 y Ba-133, con adquisición a 1 Hz. Las configuraciones fueron:

• Vuelo estacionario: unos 5 min a unos 90 m, usados para la comparación controlada de actividades y fondo.
• Pasada única: a unos 30 m de altitud y unos 10 m/s, en tres escenarios: solo Cs-137, solo Ba-133 y una fuente doble ubicada en el mismo punto, con una separación inferior a 20 cm.

Los valores de referencia fueron independientes. Las actividades de las fuentes procedían de los ensayos de laboratorio del fabricante, y el fondo natural se acotó con espectroscopía in situ en 23 posiciones del campo de entrenamiento. A 10 m/s, seis segundos equivalen a un recorrido de unos 60 m (cálculo propio). La inferencia se obtiene, por tanto, de un tramo corto de vuelo en torno a la fuente.

A modo de contexto, los trabajos bayesianos anteriores sobre prospecciones con vehículos por carretera no son directamente comparables. Un estudio que circulaba a 50 km/h con muestreo de 1 s halló que las desviaciones previstas de la actividad y la posición se mantenían en torno al 30 %, debido a la incertidumbre de los datos. Ese estudio usó otros detectores y otra geometría. Aun así, muestra la distancia entre un modelo aproximado de fuente puntual y un modelo de espectro completo específico de la plataforma.

05Por qué importa y qué queda sin demostrar

En un accidente de reactor o en una detonación nuclear, las aeronaves pueden ser el único medio de llegar a la zona afectada. El equipo sostiene que los datos obtenidos en vuelo pueden orientar decisiones sobre evacuación, confinamiento y distribución de yoduro de potasio. En seguridad nuclear, la propiedad útil es la resolución de mezclas. Distinguir una fuente médica de Cs-137 de un inventario mixto o industrial, y hacerlo en una sola pasada, cambia lo que hace a continuación un equipo de búsqueda. Los autores citan además como aplicaciones la vigilancia rutinaria de instalaciones, la cartografía geoquímica e incluso la abundancia de radionúclidos extraterrestres. Los materiales de prensa añaden la recuperación ambiental y la prospección minera.

Las limitaciones se aprecian en el propio diseño del artículo:

• Pocos nucleidos: Los ensayos usaron dos fuentes puntuales selladas. Una liberación real implica decenas de isótopos, penachos extensos y deposición en el suelo, por lo que tanto la biblioteca de plantillas como el problema de inferencia crecen.
• Una plataforma y un emplazamiento: Las funciones de la aeronave son específicas de un modelo de helicóptero. Otro fuselaje necesita su propio modelo Monte Carlo validado y su propio precálculo.
• Geometría favorable: Las trayectorias de vuelo pasaron justo por encima de las fuentes, cuya presencia se conocía. Buscar a ciegas en una zona desconocida es una prueba más difícil.
• Porcentajes: Las cifras de «dos puntos porcentuales» y «dos metros» proceden de los comunicados de prensa. El artículo habla de «escala porcentual» y «escala métrica». Conviene consultar las tablas por escenario del artículo antes de citar una cifra concreta.

Los próximos pasos no se detallan en las fuentes que he revisado. Los razonables serían más nucleidos, fuentes extensas y distribuidas, y pruebas en condiciones operativas. Los datos están archivados en Zenodo, de modo que grupos independientes pueden probar el método con los mismos espectros. El resultado principal es que una factorización consciente de la física hizo que un tratamiento bayesiano riguroso fuera lo bastante rápido como para aplicarlo a los datos de un vuelo de pocos segundos.