ES2640934B2 - Procedimiento de obtención de la perturbación geomagnética local a latitudes medias - Google Patents

Abstract

La invención describe un dispositivo y un procedimiento para determinar la perturbación geomagnética local a latitudes medias mediante la eliminación de las diferentes componentes de campo geomagnético no perturbado, que resuelve los problemas existentes en el estado de la técnica conocido. Para ello en la presente invención se describe un modelo de día en calma geomagnética a latitudes medias para la componente horizontal (H) de campo geomagnético. El procedimiento implica la caracterización de día en calma o día perturbado mediante una valoración de la bondad del ajuste del modelo inventado a los datos medidos de la componente horizontal de campo magnético local, una vez sustraída una curva auxiliar de tendencia de los valores de noche. El procedimiento puede implementarse en el dispositivo para determinar la perturbación geomagnética local en tiempo real.

Description

PROCEDIMIENTO DE OBTENCiÓN DE LA PERTURBACiÓN GEOMAGNÉTICA LOCAL A LATITUDES MEDIAS

SECTOR DE LA TÉCNICA

La presente invención se encuadra en el campo técnico de los dispositivos y métodos de producción de índices geomagneticos con propósitos de meteorología espacial: más concretamente, la invención se refiere a los índices que ofrecen una medida de la perturbación geomagnética ocasionada en una zona localizada en la superiicie de la Tierra por fuentes externas al planeta. La técnica propuesta es capaz de ofrecer información de la perturbación geomagnética a latitudes medias (alrededor de 40°) de alta resolución y en tiempo real.

ESTADO DE LA TÉCNICA

Problema técnico a resolver

Es conocido que el Sol tiene una actividad continua que afecta a la vida en la Tierra. La manifestación más evidente del efecto que el Sol produce en el planeta es la energía que llega hasta su superficie en unos 8 minutos en forma de radiación electromagnética, distribuida en diferentes rangos del espectro electromagnético, aunque la mayor aportación se encuentra en el visible.

Ademas de radiación, la corona solar esta emitiendo de forma continua un flujo de particulas cargadas, conocido como viento solar, que se encuentra en estado de plasma y que viaja a través del espacio interplanetario ligado al propio campo magnético solar, interaccionando con todos los planetas, cometas y asteroides que se encuentra en el interior de la heliosfera. A la velocidad típica de 400 km/s este plasma magnetizado llega a la Tierra en unos 4 días e interacciona con el campo magnético dipolar terrestre distorsionandolo. La interacción resulta en un sistema complicado de movimientos de convección de plasma y corrientes eléctricas en el interior de la magnetosfera terrestre que dan lugar a la topología completa que ésta adquiere como situación de equilibrio o de "estado en calma~

Sin embargo, en el Sol, como estrella altamente dinámica con actividad ciclica de unos 11 años, también tienen lugar, de forma esporádica, fenómenos explosivos y violentos en las estructuras magnéticas y plasma que sostienen, dando lugar a distintos tipos de fenómenos como son: fulguraciones, prominencias, eyecciones de masa coronal (eME, acrónimo del inglés Coronal Mass Ejection). etc. En ellos se libera una cantidad ingente de energía, plasma y campo magnético que pueden afectar de forma importante el entorno terrestre cuando se dirigen hacia la Tierra. Asimismo, el viento solar rápido (HSS, del acrónimo inglés High Speed Stream) liberado de estructuras solares de mayor tiempo de vida que las CMEs, conocidas como agujeros coronales (CH. del acrónimo inglés Coronal Hale), pueden también producir efectos notables en el entorno terrestre.

Al ya complicado escenario presentado, hay que añadir que también pueden producirse partículas de alta energía (del orden de MeV) que viajan a velocidades cuasi-relativistas y que tardan unas pocas horas en alcanzar la Tierra. Es lo que se conoce como sucesos de partículas energéticas solares (SEP, del acrónimo inglés Solar Energetic Partides), que se producen en las fulguraciones solares o en las ondas de choque de las CMEs. A estos sucesos de particulas de origen solar muy energéticas se ana den partículas de origen extrasolar. los rayos cósmicos, cuyo rango energético se solapa con el de las partículas solares más energéticas si bien la tendencia de su fluencia está anti-correlacionada con el ciclo solar.

Bajo condiciones de perturbación solar se producen cambios importantes en todas las escalas espaciales del entorno terrestre: magnetosfera, ionosfera, atmósfera, etc., y los efectos son detectables también en la superficie. Uno de los procesos físicos que tiene lugar cuando una eME interplanetaria alcanza la magnetosfera terrestre es la entrada de partículas en su interior y la consiguiente modificación del campo geomagnético, lo que se conoce como tormenta geomagnética.

Los efectos que las particulas energéticas, la radiación electromagnética y la modificación del campo geomagnético producen en el entorno terrestre. son muy variados y con escalas de tiempo muy diferentes. Pueden afectar de forma directa por ejemplo a la salud, o indirectamente, a través de fallos en tecnología que la sociedad usa diariamente, por ejemplo, en las comunicaciones radio. sistemas de navegación terrestre y de

posicionamiento, e incluso en infraestructuras críticas como son eléctrica, el transporte ferroviario o las grandes conducciones combustible entre paises.

las redes de potencia portadoras de gas o

S

Antecedentes

10 15

Se sabe que la medida del campo magnético terrestre tiene varias contribuciones, como son: la generada por la propia Tierra por el efecto de dinamo interna, la generada por la actividad solar habitual (también denominada como ~Sol en calma') y la generada por los fenómenos solares mencionados anteriormente que suceden de manera esporádica. El objeto de la presente invención se centra en la medida de las perturbaciones magnéticas a latitudes medias, para cuya determinación y correcta interpretación existen principalmente dos problemas: uno relacionado con la componente de ~Sol en calma" y el otro relacionado con las perturbaciones causadas por la actividad solar que ocurren de forma irregular.

20

El primer problema esta relacionado con el campo magnético generado por la llamada componente de ~Sol en calma", producido por un sistema equivalente de corriente eléctrica de origen ionosférico que se produce en el lado diurno del planeta y que presenta vorticidad a escala global, produciéndose uno en cada hemisferio y cuyos focos están localizados a latitudes medias. Este sistema no permanece fijo sino que deriva en longitud y latitud a lo largo del día y presenta una variabilidad día a día Que provoca incertidumbre en la determinación correcta tanto de la intensidad como de la posición del foco.

25 30

El segundo problema viene causado porque las perturbaciones magnéticas que los fenomenos solares producen en la superficie terrestre a escala global en el rango de latitudes medias están bien representadas por indicadores globales disenados al respecto, pero no asi las perturbaciones magnéticas locales, cuyo desconocimiento en ocasiones da lugar a interpretaciones erróneas, incluso confundiéndolas en ocasiones con efectos globales. Por ejemplo, está profundamente asumido que a latitudes medias/bajas (rango de latitudes comprendido aproximadamente entre los 15(l y 55(l Norte y entre 15(l y 55° Sur) la única fuente capaz de producir perturbaciones geomagnéticas de intensidad importante a escala global es la debida a la corriente magnetosférica de anillo, producidas por la

intensificación de la población de partículas de baja energía que circulan alrededor de la Tierra bajo condiciones adversas de meteorología espacial (Gonzalez, W. D., et al. "What is a geomagnetic storm?" Journal of Geophysical Research: Space Physics 99, A4 (1994): 5771-5792). Sin embargo, publicaciones recientes (Cid, C., et al. "A Carrington-like geomagnetic storm observed in the 21st century." Journal of Space Weather and Space Climate 5 (2015): A16 y Saiz, E., et al. "Searching for Carrington-Jike events and their signatures and triggers.~ Journal of Space Weather and Space Climate 6 (2016): A6) demuestran Que a dichas latitudes medias pueden existir perturbaciones magnéticas locales de intensidad importante muy simHares, e incluso bastante mayores que las perturbaciones causadas por la corriente del anillo. Estos recientes resultados de investigación establecen con especial importancia el papel Que desempeñan los rndices geomagnéticos locales en el campo de la Meteorología Espacial y crean nuevos requisitos en la precisión y fiabilidad de los mismos .

La componente de "Sol en calma" viene determinada por la contribución al campo geomagnético del sistema equivalente de las corrientes ionosfericas generadas por un efecto de dinamo en la zona iluminada de la Tierra. Posee un patrón de circulación que presenta un vórtice en cada hemisferio terrestre que deriva a lo largo del día, y sus focos se encuentran localizados aproximadamente en el rango de latitudes medias cercanas a 40°, tanto en el hemisferio Norte como en el hemisferio Sur. Este sistema de corrientes se identifica normalmente de manera mensual mediante el promedio de los cinco días del mes que menos perturbacion geomagnética presentan, los cuales son internacionalmente establecidos (lOO, International Quiel Oays). Esta representación y determinación de la variación diaria del día en calma es la que se usa actualmente, pero presenta varias desventajas. Cuando el propósito es obtener la perturbación global, dichas desventajas pueden ser irrelevantes, sin embargo, cobran especial importancia cuando el objetivo final es obtener la perturbación geomagnetica localizada en una zona o región (índice geomagnético local) como el de la presente invención. Las principales desventajas vienen determinadas principalmente por tres motivos:

a) la perturbación magnética en una zona puede llegar a presentar diferencias

importantes respecto de otra y, por tanto, la definición internacional de tiempo

en calma, la cual esta basada en un índice global (índice Kp), puede conducir a error cuando la perturbación geomagnética afecta solamente a una zona muy localizada;

b) la determinación de días en calma usando como tiempo base el tiempo universal (Ur, Universal Time, por sus siglas en ingles). Aunque es la usada en los IQOs, no es, por definición, una forma practica de detallar un dia en calma para una zona con caracter local ya que, inevitablemente, el dra local (con tiempo base el tiempo local (ML T, Magnetic Local Time, por sus siglas en ingles)) estara desplazado respecto del tiempo universal;

el la variabilidad inherente que existe entre días de calma, que hace que la determinación mediante el promediado de varios días pueda resultar en una representación del día en calma local muy diferente de la real. Ademas, para latitudes medias que son las de localización de los focos de este sistema de corrientes ionosféricas este problema es aún mas importante.

El objeto de la presente invención es el desarrollo de un procedimiento y un dispositivo que, a partir de la medida del campo geomagnético local a latiludes medias, sea capaz de estimar la componente de perturbación debida a fenómenos asociados a la actividad solar

esporadica. Para ello se tienen en cuenta los problemas descritos anteriormente y en especial los que afectan a latitudes medias.

Se conocen antecedentes sobre dispositivos y procedimientos para predecir la variación del campo geomagnético local, como el que se describe en el documento WO 2015/013499 A1 . Sin embargo, como se expone en esta memoria, las técnicas de tratamiento de las medidas de campo geomagnético para la obtención de la perturbación geomagnética local existentes en el estado de la técnica presentan notables diferencias con respecto a las que se exponen en la presente invencion.

DESCRIPCiÓN DE LA INVENCiÓN

La presente invención describe un procedimiento y dispositivo capaces de obtener la componente de perturbación geomagnetica local a latitudes medias a partir de la medida del campo magnético en la superficie terrestre. Dicha medida debe de ofrecer al menos la componente horizontal del campo geomagnetico en la superficie de la Tierra y puede

obtenerse directamente desde sensores geomagnéticos debidamente calibrados y aislados de cualquier perturbación electromagnética artificial. Dicha medida también puede obtenerse de los datos que ofrecen en tiempo real los observatorios geomagnéticos a través de conexiones remotas.

En la medida del campo magnético en la superficie terrestre se pueden englobar principalmente tres componentes; una componente interna, componente 1), compuesta a su vez por la componente principal, cuyo origen está en el proceso de dinamo que tiene lugar en el interior de la Tierra y una componente debida a la corteza y manto terrestres; una componente de "Sol en calma" , componente 11), originada por el sistema equivalente de las corrientes ionosféricas de naturaleza diurna; y una componente de perturbación, componente 111), debida a fuentes externas al entorno terrestre, como es la que se produce de manera irregular y esporádica por causa de la actividad solar.

La presente invención trata sobre la determinación de las primeras dos componentes, para obtener la tercera, la perturbación geomagnética local generada por la actividad solar, con objeto de poder predecir situaciones potencialmente peligrosas como las que se han enumerado en el epígrafe Problema tecnico a resolver.

La técnica usada en la presente invención determina en primera instancia la componente 11), que posee variación (de magnitud considerable, del orden de decenas de nanoteslas) de periodo de tiempo mas corto (un dial, frente a periodos del orden de meses, o mayores, de la componente 1). Este método contrasta con el modo tradicional de proceder, en el que se suele obtener en primera instancia la componente 1, por poseer variaciones mas lentas. Esto es posible gracias al método usado en la presente invención para la obtención de la componente 11), que mediante un modelo de comportamiento de día en calma es capaz de determinar de manera unívoca la curva de variación diaria para los días no afectados por perturbaciones externas.

Además, la correcta determinación de la componente 11) durante un conjunto de dias prolongado en el tiempo, permite obtener la componente 1). Aunque la componente 1) posee variaciones lentas, su comportamiento es difícil de pronosticar. Existen modelos geomagnéticos, como es el modelo IGRF (/nternational Geomagnetic Reference Field, Campo Geomagnetico de Referencia Internacional), que proporcionan estas componentes para valores pasados, presentes, y con proyección en algunos afias hacia el futuro (no mas de cinco), pero su precisi6n, tanto espacial como temporal, no suele ser suficiente para muchas aplicaciones.

Como consecuencia, en la presente invención, la componente 1) se obtiene de los datos medidos, pero no directamente, sino después de realizar una parte del procesamiento para la obtencion de la componente 11. Esta forma de proceder aporta a la invencion dos caracteristicas importantes, por un lado, independencia absoluta del procedimiento sobre datos externos, y por otro lado mayor precisión para la obtención de [as componentes para la localización concreta de la medida, pues no se depende de datos derivados con caracter global como ocurriria con el uso de Indices o modelos geomagnéticos.

Para determinar la componente 11), como se ha mencionado, se realiza un ajuste de un modelo de día en calma a los datos representativos de la variación diaria. Para ello, previamente se sustrae de manera auxiliar una curva suave de tendencia de los valores de noche, es decir, una curva que une los valores del campo geomagnético medido en las horas de noche, menos afectadas por la variación diaria. Analizando la bondad del ajuste del modelo a los datos se obtiene una valoración de la calma geomagnética de los mismos, es decir, es posible conocer cuando los datos poseen una curva de variación diaria representativa de día en calma. Esta característica hace posible la determinación de la componente 1) a partir de la acumulación de días calificados como de calma .

La componente 1) se puede determinar de varias formas, por ejemplo mediante una curva de interpolación de valores promedios mensuales o anuales de los días en calma geomagnética seleccionados de la forma descrita anteriormente. Una variante puede ser usar como alternativa sólo los días de calma oficiales (International Quie! Days) establecidos internacionalmente por GFZ-Potsdam. Otra variante puede ser usar los valores promedio mensuales o anuales de los observatorios geomagnéticos mas cercanos a la zona de medida.

La componente 111) de perturbación geomagnética (objetivo final de la presente invención) se obtiene eliminando las componentes 1) y 11) a los datos medidos. Para los dias de calma

geomagnética en los que el modelo ofrece un ajuste bueno a los datos, la componente 11) que debe de eliminarse es precisamente la curva de variación diaria obtenida del ajuste del modelo. Para los días en los que el ajuste del modelo no es bueno, se presentan varias alternativas. Se puede eliminar la curva de variación diaria del día en calma más próximo

5 en el tiempo, el promedio de varios días o incluso no eliminar ninguna curva de variación diaria.

En el estado de la técnica existen métodos de interpolación para datos de medida geomagnética, por ejemplo, en el documento WO 2015/013499 A1 , pero el objetivo suele 10 ser siempre la determinacion de la componente 1), de variación mas lenta y menor complejidad que la curva de variación diaria a latitudes medias.

El carácter innovador de la presente invención es posible gracias a la aplicación del modelo de dia en calma. Dicho modelo se define mediante una curva de interpolación suave que 15 pase por puntos que se encuentran fijados a determinada posición temporal del dia. El modelo representa la variacion geomagnética sufrida en la zona iluminada de la Tierra, zona de día, por tanto, se diferencian dos tipos de puntos, unos fijados a intensidad nula en las horas de noche, puntos (N), y otros libres en cuanto a intensidad. en la zona de día, puntos (P), y que ofrecen como parametros los grados de libertad necesarios para el ajuste

20 del modelo. Para que la curva de variación diaria pueda tomar cualquier forma posible típica de las latitudes magnéticas medias, y al mismo tiempo, impedir que el modelo ofrezca un ajuste bueno a otras variaciones geomagnéticas de perturbación externa que deben de quedar excluidas, los puntos (P) deben de ser cuatro. estar colocados en la horas de dia (entre las 6 h Y 18 h) Y con una distancia entre ellos entre una y tres horas.

DESCRIPCiÓN DE LOS DIBUJOS

Para complementar esta descripción, y con objeto de ayudar a una mejor comprensión de las características de la invención, se acompaña, como parte integrante de dicha

30 descripción , un j uego de dibujos en donde. con carácter ilustrativo y no limitativo, se ha representado lo siguiente:

la figura 1 muestra en el panel al la evolución temporal de la componente horizontal del campo geomagnético local, componente H, en un periodo concreto de tiempo (del 9 al 13 de marzo de 2001 " en una localización concreta (observatorio de San

Pablo de los Montes, Toledo, Espai'la) a modo de ejemplo y además incluye las curvas implicadas en el procedimiento objeto de la invención y en el panel b) la figura muestra el resultado de aplicar el procedimiento objeto de la invención, es decir, la componente 111) de perturbación geomagnélica local: la figura 2 ilustra el diagrama de flujo del procedimiento para la determinación de la perturbación geomagnética local en la superlicie terrestre; la figura 3 es un detalle ampliado del ajuste de los datos a la curva de variación diaria en un intervalo temporal de un día de duración: la figura 4 representa un esquema del dispositivo asociado a la invención sobre una de las realizaciones preferentes de la misma.

A continuación se incluye una lista con los distintos elementos que se representan en las figuras:

•

11 (trazo continuo); componente horizontal del campo geomagnético medido:

•

12 (trazo de puntos): curva de tendencia de valores de noche calculada:

•

13 (trazo discontinuo de guion grueso): curva de variación diaria de la componente horizontal del campo geomagnético local calculada;

•

14 (trazo discontinuo de guion fino): línea base asociada a la componente horizontal del campo geomagnético local asociada al campo magnético generado en el interior. la corteza y el manto terrestres;

•

15: perturbación magnética de la componente horizontal del campo magnético local;

•

20: bloque del diagrama de flujo de pre-procesado de la componente horizontal de campo geomagnélico medido:

•

21: bloque del diagrama de flujo para la determinación de la curva de variación diaria:

•

22: bloque del diagrama de flujo para la determinación de la linea base;

•

23; bloque del diagrama de flujo para el calculo final que determina la perturbación geomagnetica local:

•

211: etapa de determinación y sustracción de la curva auxiliar de tendencia de los valores de noche:

•

212: etapa de ajuste del modelo de variacion diaria a los datos procedentes de la etapa anterior;

•

213: etapa de decision sobre la determinación o no de los datos actuales como pertenecientes a día geomagnetico en calma;

•

214: etapa de búsqueda del día geomagnético en calma mas cercano al actual en la base de datos 225, que se ejecuta tras una respuesta negativa en 213:

•

225: base de datos de los días geomagnéticos en calma guardados tras respuestas positivas en 213:

•

226: etapa de cálculo de la línea base;

•

31 (trazo continuo): curva de los datos a ajustar por el modelo de día en calma, calculados a partir de la componente horizontal de campo magnetico y una vez sustraida la curva auxiliar de tendencia de los valores de noche;

•

32 (trazo discontinuo de guion grueso): curva de variación diaria obtenida del ajuste del modelo de día en calma:

•

331 · 334: puntos de control del spline cúbico, de tipo (N), fijos en posición temporal (horizontal) y en intensidad (vertical), del modelo de día en calma ;

•

341·344: puntos de control del spline cúbico, de tipo (P), fijos en posición temporal pero libres en intensidad, del modelo de día en calma;

•

40: dispositivo de obtención de la perturbación geomagnética local;

•

41 : unidad de entrada de datos del d isposilivo 40:

•

42: CPU del dispositivo 40:

•

43: unidad de almacenamiento de datos del dispositivo 40;

•

44: unidad de salida de datos del dispositivo 40.

DESCRIPCiÓN DE UNA FORMA DE REALIZACiÓN DE LA INVENCiÓN

Como es sabido, la medida del campo magnético en la superficie terrestre posee diferentes contribuciones que, en determinadas situaciones, pueden ser identificadas como tres componentes distintas:

1) una componente interna, compuesta por la componente principal, cuyo origen está en el proceso de dinamo que tiene lugar en el interior de la Tierra y una componente debida a la corteza y manto terrestre:

, ,

11) una componente de MSol en calma", originada por el sistema equivalente de las corrientes ionosféricas de naturaleza diurna; 111) una componente de perturbación. debida a fuentes externas al entorno terrestre, como es la que se produce de manera irregular y esporádica por causa de la actividad solar.

La presente invención trata sobre la determinación de las primeras dos componentes 1), 11) para obtener la perturbación geomagnética local generada por la actividad solar, componente 111), con objeto de poder predecir situaciones potencialmente peligrosas como las que se han enumerado en el epígrafe Problema técnico a resolver.

Visión global del procedimiento

Para dar una visión global y una mejor comprensión del procedimiento se han realizado las figuras 1 y 2 que muestran las curvas implicadas en el procedimiento y un diagrama de flujo del mismo, respectivamente.

En la figura 1 se ha representado en el panel superior (la) la evolución de la componente horizontal del campo geomagnético local medido, curva (11), junto con las diferentes curvas implicadas en el procedimiento, (12)-(14). En el panel inferior (lb) se ha representado el resultado de aplicar el procedimiento descrito en la presente invención, curva (15). En el caso de la figura la, los datos medidos, representados por (11), han sido proporcionados por el observatorio geomagnetico de San Pablo de los Montes en Toledo, Espar'la, identificado internacionalmente con las siglas SPT (código dado por el organismo IAGA, International Association of Geomagnetism and Aeronomy). El número entre paréntesis que figura en el titulo del eje de ordenadas corresponde al valor (expresado en nanoteslas, nT) que hay que sumar a los valores ofrecidos en el eje de ordenadas para obtener el valor absoluto de la medida. En el eje de abscisas inferior se representa el tiempo en el que se hace la medición, dado en hora local magnética (ML n. El eje de abscisas superior representa tiempo expresado en tiempo universal (UT). La diferencia de 16 minutos que se aprecia entre ambas escalas de tiempo corresponde a la diferencia temporal eqUivalente a la diferencia entre la longitud geográfica de referencia del meridiano

de Greenwich y la longitud del observatorio SPT donde se realiza la medida. La ventana temporal de la figura 1 tiene una duración de cinco días.

En la figura 2 se muestra el diagrama de flujo del procedimiento completo. Como puede observarse existe un bloque inicial (20), dos bloques o sub-procedimientos principales orientados verticalmente (21) y (22), Y un bloque final (23). Los bloques principales llevan a cabo la determinación de la componente 11) en el bloque (21) y la determinación de la componente 1) en el bloque (22). El bloque final (23) no es mas que una sustracción de las componentes 1) y 11), resultados de (21) Y (22), a la componente H de los datos medidos.

El primer bloque (20) es un pre-procesado de los datos medidos donde se realiza, entre otras operaciones. una asociación de los valores a su referencia temporal en el sis1ema local magnetico ML T. Para ello es necesario que la medida del campo geomagnético se acompañe de alguna referencia temporal en la que se realiza la medida, como puede ser la hora universal (UT). Para realizar la conversión a tiempo ML T es preciso aportar la posición en longitud del lugar en la que la medida se ha realizado. Asimismo, en el preprocesado se incluye la transformación del sistema de coordenadas, si fuera necesario, desde las componentes de campo proporcionadas por el sistema de medida a un formato del que sea posible extraer la componente horizontal (H).

La evolución del campo geomagnético local medido, dado por la curva (11) de la figura 1, supone la superposición de las tres componentes: 1). 11) Y 111) previamente mencionadas. Para obtener lo que de ese valor corresponde a la componente de perturbación geomagnética (componente 111)). medida localmente en un punto concreto de la geografía , y debida a fuentes externas, es decir, como consecuencia de la actividad solar esporádica. es necesario eliminar las componentes 1) y 11) descritas anteriormente. En la figura 1a) también se representan otras tres curvas: (12) es la curva auxiliar de tendencia de valores de noche. usada en el bloque (211); (13) es la curva de variación diaria (componente 11), y (14) es la linea base representativa de la componente 1).

La técnica usada en la presente invención determina en primera instancia la componente 11) por ser la que posee variación (de magnitud considerable) de periodo de tiempo mas corto (un dial, frente a periodos del orden de meses o mayores de la componente 1).

Además, la correcta determinación de la componente 11) durante un conjunto de dias prolongado en el tiempo, permite obtener la componente 1). Aunque la componente 1) posee variaciones lentas, su comportamiento es difícil de pronosticar. Existen modelos geomagneticos, como es el modelo IGRF, que proporcionan estas componentes para valores pasados, presentes, y con proyección en algunos años hacia el futuro (no mas de cinco), pero su precisión, tanto espacial como temporal, no suele ser suficiente para muchas aplicaciones. Además, estos modelos se obtienen de procesos en los que previamente se ha eliminado la componente de perturbación externa (la componen1e 111) mencionada) mediante una interpretación global de la misma, por ejemplo, mediante el uso de índices indicadores de perturbación geomagnética no local sino global como son el índice Osi o el índice Kp. Por tanto, cuando el propósito es obtener la perturbación geomagnética local con precisión, como en la presente invención, no es adecuado hacer uso de dichos modelos geomagnéticos. Como consecuencia, en la presente invención, la componente 1) es obtenida a partir de los datos observados en la misma localización donde se realiza la medida. Esta característica se puede apreciar mejor observando el diag rama de flujo de la figura 2: la componente 1, obtenida en el bloque (22), se obtiene de los datos originales medidos (20), pero no directamente, sino después de realizar una parte del procesamiento para la obtención de la componente I1 en el bloque (21). Esta forma de proceder aporta a la invención dos características importantes, por un lado, independencia absoluta del procedimiento sobre datos externos, y por otro lado mayor precisión para la obtención de las componentes para la localización concreta de la medida, pues no se depende de datos derivados con carácter global como ocurriría con el uso de indices o modelos geomagnéticos.

Modelo de dfa en calma

La componente 11) de variación diaria se determina gracias a un modelo de día en calma (relativo a la componente horizontal de campo geomagnético) válido para una latitud magnética cercana a la de la localización del foco del sistema de corrientes ionosféricas (aproximadamente entre los 35° y 45° desde el ecuador en ambos hemisferios). de la presente invención. El modelo se caracteriza por estar definido mediante un spline cúbico, es decir, una curva diferenciable definida a tramos mediante pOlinomios cúbicos (polinomios de tercer grado). Los puntos de control son el lugar de unión de los tramos y el lugar donde se imponen ciertas condiciones que hacen que la transición ocurra de forma suave. Para un spline cubico, dichas condiciones son tres y se deben de cumplir en cada punto de control. Estas tres condiciones implican continuidad de las funciones en los puntos de control, y se traducen en que tanto la evaluación del polinomio en el punto de control, como el de su primera derivada y su segunda derivada deben ser 19 uales para los dos polinomios existentes a cada lado del punto de control.

El modelo de día en calma para esta forma de realización de la presente invención esta compuesto por ocho puntos de control distribuidos en el tiempo a lo largo de un día, 24 horas, de manera fija en determinadas horas del día, tal y como se muestra en la figura 3, representados por rombos y asteriscos, diferenciando dos tipos, paramétricos (P) y neutros (N), respectivamente. La figura 3 muestra dichos puntos en un ejemplo de día en calma geomagnética; a los datos mostrados como elemento (31), se ha ajustado el modelo de variación diaria, obteniendo la curva (32). Los datos mostrados en (31) no son los datos originales medidos, resultado del bloque (20), sino que previamente se ha sustraído una curva auxiliar en la etapa (211) tal y como se describe en el siguiente apartado. Los cuatro puntos (N) corresponden a dos puntos en las primeras cuatro horas del día y a otros dos puntos en las cuatro últimas horas del día, es decir, ocupando las horas de noche local (de ahí su denominación con ta letra N). Estos puntos (N) estan caracterizados por estar situados a una intensidad de campo magnético fijada en un valor de cero respecto del valor de base del modelo de día en calma. Los restantes cuatro puntos de control, puntos (P), estan situados en la zona ceniral, correspondientes a las horas de dia del lugar y poseen la característica de tener una intensidad variable, ofreciendo cada uno de ellos un parámetro libre al modelo que permite el ajuste de los datos al mismo. Además, en los dos puntos (N) de los extremos (los correspondientes a O h Y 24 h) se impone la condición de que la primera derivada sea nula. El modelo definido de esta forma posee cuatro grados de libertad, debido a los cuatro parámetros libres (P) de la zona central en las horas de día.

Determinacion de la curva auxiliar de tendencia de los valores de noche

Puesto que se parte de datos observacionales medidos en la superficie terrestre (que, como se describió anteriormente, son el resultado de la superposición de cuatro componentes), para poder realizar un ajuste del modelo de día en calma, con el objetivo de obtener la componente 11), los datos a los que se ajusta deben de contener la mInima aportaci6n de las otras dos componentes restantes; 1) y 111). La técnica adoptada en la invención pasa por determinar primeramente una aproximación de las dos componentes 1) y 111) superpuestas, que se restan de los dalos originales de manera auxiliar, y que permite realizar un ajuste del modelo de día en calma s610 a los valores de variación diaria que mejor representen la componente 11). Esta curva de aproximación, a la que se ha denominado curva de tendencia de los valores de noche (elemento (12) en figura 1 b), es una curva de interpolación suave entre un s610 punlo diario obtenido a partir del valor medio de las horas de noche, descrita en detalle mas adelante.

Para entender la importancia y a su vez la validez de la determinación de la curva auxiliar hay Que tener en cuenta el posible comportamiento de las otras dos componentes. La componente 1), como ya se dijo, posee variaciones lentas comparadas con la variación diaria y, por tanto, una evaluación del tipo Que se plantea con un solo punto diario (representando el valor medio de la zona de noche), es mas que suficiente para asegurar su correcta descripción. La componente 111) de perturbación externa contempla variaciones geomagnéticas de muy diversa indole, Que para los propósitos de la presente invención podemos diferenciar en dos tipos: III-A) variaciones de intensidad importante que perturban el campo magnético de manera irregular y esporádica, dificilmente predecibles y Que impiden una correcta determinación de la componente 11) de variación diaria; y 111-8) variaciones suaves pero mantenidas, de mas de un dia de duración y que pueden estar presentes en dia de calma geomagnética, debida a perturbaciones en los días previos. Las variaciones de tipo III-A) poseen cambios rápidos e importantes que resultaran en un ajuste del modelo de día en calma no válido, y por tanto, aunque no sea posible una determinación de la misma, sí que se puede constatar su presencia mediante la medida de la bondad del ajuste, como se vera mas adelante. La variación III-B) se traduce en un aumento, o disminución suave, del nivel diario del campo geomagnético medido, y por tanto la evaluación medianle una curva obtenida por un punto diario también asegura su correcta delerminacion previa al ajuste.

En este punto hay que aclarar que el hecho de determinar, con carácter previo al ajuste del modelo de día en calma, la curva de tendencia (12) como aproximación de las componentes 1) y 111-8), no significa que dicha variacion sea eliminada de forma definitiva de los datos observacionales originales, pues como se vera mas adelante, ésta sólo se usa como curva auxiliar para determinar la componente de variaci6n diaria. Para aclarar en lo que viene la terminología usada, se diferenciara entre eliminación, término con el que se quiere sena lar que es la operación de sustracción definitiva de una contribución y con la que se obtiene un resultado intermedio o fina 1; yel término sustracción, con el Que se senala que se hace la operación aritmética de resta, pero sólo como proceso intermedio Que no implica la eliminación definitiva sobre los datos originales.

La curva auxiliar de tendencia de los valores de noche (12) se obtiene a partir del valor promedio alcanzado en un conjunto de horas de noche que se elige de antemano dentro de un rango simétrico (máximo 4 horas) de horas anteriores y posteriores a la medianoche local de tránsito entre un día y el siguiente. En el caso del ejemplo para cinco dias mostrado en la figura 1b), los valores medios para las seis horas de noche están representados por asteriscos, la zona de las horas de noche utilizadas, desde las -3h a las +3h respecto de la medianoche local, están marcadas con fondo de rayas oblicuas. La curva suave que pasa por dichos valores medios de noche se ha realizado mediante un spline cúbico, aunque podrían usarse otros métodos de interpolación mediante curva suave similares. Una vez determinada la curva de tendencia de los valores de noche, ésta se sustrae de los datos medidos y se procede al ajuste del modelo de día en calma sobre el resultado de la sustracción.

Calificación y almacenamiento de dias geomagnéticos en calma

El ajuste del modelo de día en calma a los datos se puede realizar mediante la minimización de una función como la del error cuadrático medio y evaluando la bondad del ajuste por medio del chi cuadrado. El método usado es un proceso iterativo en el que se establecen unos valores iniciales para los cuatro parámetros del mismo y se evalúa la bondad del ajuste. En el proceso, se modifican los parámetros hacia el valor que minimice la función de error, de forma que estos deben de converger hacia un valor estable siempre que los datos posean una variación diaria típica de día geomagnético en calma . El método iterativo puede acabar mediante un número de iteraciones (si no converge) o mediante el establecimiento de un umbral en el incremento relativo del valor de la bondad del ajuste entre una iteración y la siguiente, de manera que para incrementos inferiores al umbral, asegure que se ha finalizado el proceso de convergencia del mismo.

El vaJor final de la bondad de ajuste (valor del chi cuadrado) se usa en esta etapa (213), ver figura 2, para discriminar entre dia en calma del que no Jo es. Si el valor es suficientemente elevado indica que el ajuste no es bueno, seguramente porque los datos poseen perturbación geomagnética y, por tanto, no puede ser considerado como dia en calma. La curva así determinada no debe de usarse como componente 11). Por otro lado, el resultado con un valor del ch; cuadrado pequeno confirma que la curva resultante del ajuste es una curva valida de día en calma y, por tanto, no sólo se debe de eliminar como componente 11), sino que debe de ser almacenada en una base de dalas (etapa (225) en figura 2) que permitirá determinar la componente 1).

Para realizar esta discriminación de forma adecuada, se establece un valor umbral para la bondad del ajuste, que puede ser determinado de varias formas dependiendo de la cantidad de datos previos conocidos. Si se tienen datos históricos de manera continuada durante varios af"los. es posible realizar un estudio de la distribución de probabilidad de los valores del chi cuadrado obtenidos para los ajustes del modelo a cada día y establecer el umbral en un valor suficientemente alto como para asegurar una probabilidad cercana a cero de encontrar un día en calma para valores por encima de ese umbral. Si no se poseen datos históricos de forma continuada pero si un periodo de varios días, para establecer el umbral se puede realizar un estudio comparativo del valor de la bondad de ajuste entre días calmados y perturbados.

En el caso mas desfavorable que no se posea ningún dato previo, se sabe que el modelo de día en calma descrito en la presente invención posee cualro grados de libertad, por tanto, la distribución del chi cuadrado teórica suponiendo sólo días en calma con una distribución de error normal, puede servir de referencia para establecer el umbral en un valor que indique una probabilidad de encontrar un día en calma superior al 99.9%.

En el caso particular de la presente forma de realización de la invención, el umbral ha sido obtenido a partir de la función de distribución del eh; cuadrado que resulta del ajuste de los datos desde 1997 hasta 2012, estableciendo el umbral en 40.

Cuando un valor de la bondad del ajuste supera el umbral establecido, como se ha comentado anteriormente, la curva oblenida en el ajuste no debe usarse como representativa de día en calma pues indica la presencia de perturbaci6n geomagnética. Esto no implica la no existencia de una componente de variación diaria, pues ésta esta siempre presente y, por tanto, debe de usarse otra curva alternativa (etapa (214) del procedimiento de la figura 2). Uno de los métodos mas extendidos en uso es el de promediar varios días en calma representativos del mes (los IQDs u otros equivalentes elegidos a partir de los datos del observatorio), pero como ya se describió anteriormente, esta no es la mejor opción a adoptar debido a la variabilidad día a dia que pueden presentar las curvas de variación diaria, especialmente a las latitudes para las que se usa la presente invención. La curva de variación diaria a usar de forma alternativa que implica una menor probabilidad de error es la curva de día en calma válida mas cercana en el tiempo. Como puede apreciarse a simple vista en la figura 1b), el día 12 de marzo posee una componente de perturbación mayor que el resto de días, el resultado del ajuste del modelo de día en calma para ese día resulta en un valor del chi cuadrado superior al umbral y, por tanto, se ha considerado como curva de variaci6n diaria válida la curva del día anterior, que fue catalogado como día en calma por la bondad de su ajuste, como puede apreciarse a simple vista al comparar las curvas (11) Y (13) en la figura 1 b).

Determinación de fa componente 1)

Como ya se ha mencionado anteriormente, la componente 1) posee variaciones lentas, del orden de unos pocos nanoteslas (nT) en un periodo superior al mes y, por lanto, pueden eliminarse de forma conjunta por tener periodicidades muy superiores a la de interés, que es de un día. El método usado para ello es mediante la generación de una curva suave de interpolación. línea base, que pasa por puntos separados una distancia temporal. Con objeto de tener en consideración, si se desea, variaciones magneticas con periodicidades grandes (estacionales, anuales, de ciclo, etc.) se considera que la distancia temporal mínima entre puntos debe ser al menos de un mes, aunque otras distancias superiores pueden ser elegidas.

Los valores de los puntos por los que pasa dicha linea base pueden obtenerse de los valores promedios del intervalo temporal considerado. Estos valores promedio pueden obtenerse de diferentes formas. En su forma prioritaria de funcionamiento éstos son calculados a partir de los días geomagnéticos en calma catalogados como validos y almacenados en la base de datos (225) como producto de la etapa (213). En el caso de la figura lb), la lInea base asT obtenida con valores anuales se muestra como elemento (14). En caso de que no se disponga de suficientes elementos en la base de datos (225) debido, por ejemplo, a la puesta en funcionamiento por primera vez del dispositivo, se pueden usar los valores promedio oficiales que ofrecen los observatorios geomagnéticos o, en su defecto, los valores promedio calculados a partir de los ofrecidos por modelos geomagnéticos como eI/GRF.

Determinacion de la perturbación geomagnetica

Una vez determinada la linea base (14), representativa de la componente 1), resultado del bloque (22) de la figura 2, y la curva de variación diaria (13), representativa de la componente 11), obtenida como resultado del bloque (21), éstas son eliminadas de los datos medidos pre-procesados para obtener el valor de la perturbación geomagnética local, componente 111), representado en la figura 2 como bloque (23).

Dispositivo para la implementacion del procedimiento

El procedimiento descrito en la presente invención puede implementarse por ejemplo, en un dispositivo como el mostrado en la figura 4 como elemento (40). Los datos de entrada

(41) se obtienen mediante una conexión remota a ficheros que pueden ser actualizados en tiempo real y que contienen la componente horizontal del campo geomagnético medido al menos por dos sensores magnéticos debidamente alineados en el plano horizontaL Los sensores pueden ser de tipo fluxgate y la medida debe ser estable en el tiempo y estar aislada de interferencias por tecnología externa.

Además, se debe registrar el tiempo de referencia local magnético, para el cual se necesita conocer, al menos, la hora de tiempo universal (UT) en la que se realiza cada medida, y la coordenada de la longitud geomagnética del punto de medida. La entrada de datos puede realizarse directamente a un fichero de texto con los datos actualizados con la medida más

reciente.

El dispositivo (40) posee una unidad de control y proceso (42), CPU, del procedimiento descrito en la presente invención. La unidad de cálculo debe ser capaz de realizar valores promedios, diferencia aritmética de vectores de datos (curvas), cálculo de spline cúbicos y calculo de error cuadrático medio y de chi cuadrado. La unidad de proceso debe ser capaz

5 de implementar el procedimiento mostrado en la figura 2 y, por tanto, debe ser capaz de ejecutar sentencias de control del tipo if para la etapa (213) y bucles far para el procesado iterativo de ajuste del modelo de día en calma de la etapa (212). Ademas, para la base de datos descrita como etapa (225) se usa un sistema de almacenamiento de datos (43) al que la unidad de control puede acceder en cualquier momento.

10 Por último, la salida (44) con los datos de perturbación geomagnética de la componente horizontal local puede realizarse a ficheros de texto e imágenes que mantenga la última actualización y, a su vez, a una pantalla que muestre una grafica con los mismos.

15 Los expertos en la materia comprobarán las ventajas que se derivan de la invención descrita. La protección de esta invención se realiza al amparo de la descripción realizada , incluyendo los dibujos, y queda reflejada por lo expresado en las siguientes reivindicaciones.

Claims (4)

1. REIVINDICACIONES

   1,-Procedimiento de obtención de la perturbación geomagnetica local (15) a latitudes medias, que comprende las etapas de"

   a) obtención de la medida de la componente horizontal de campo geomagnetico local (11 l, mediante el tratamiento de los datos originales suministrados por los sensores magnéticos el magnetómetro ubicado en el lugar;

   b) determinación de una curva de variación diaria (13), que representa la contribución de origen ionosférico de la componente horizontal de campo geomagnético medido en la superlicie terrestre;

   e) determinación de una curva de linea base (14), que representa la contribución de campo magnético asociada al núcleo, corteza y manto terrestres de la componente horizontal de campo geomagnético medido en la superficie terrestre;

   d) obtención de la perturbación geomagnética local (15) como resultado de eliminar de la componente horizontal del campo geomagnético local medido (11), los valores de la contribución ionosférica (13) Y los valores de la contribución del interior de la Tierra (14);

   caracterizado por que para la obtención de la curva de variación diana (13) se analiza,

   incluso en tiempo real, si se está ante un día en calma o un día perturbado, con el Siguiente

   procedimiento'

   1) se determina una curva auxiliar de tendenCia de los valores de noche (12): 2) se sustrae la curva de tendencia de los valores de noche (12) a los valores de la componente horizontal de campo magnético medido (11):

   3) se ajusta un modelo de día en calma (212) a los datos obtenidos en el paso 2), estando el modelo de dia en calma geomagnética definido mediante una curva de spline cübico con cuatro puntos de control de intensidad variable. que funcionan como parámetros libres y que se encuentran situados con una separación fija entre ellos abarcando todas las horas de dia local, y con puntos de control con intensidad

   de valor cero para las horas de noche local;

   4) se clasifica, a partir de los valores de la bondad del ajuste realizado en el paso 3) y

   con los umbrales establecidos según el caso , si los datos pueden catalogarse como

   los de un día en calma geomagnética o los de un día perturbado;

   5) se establecen los datos ajustados del modelo de día en calma, para los días

   clasificados como día en calma, como curva de variación diaria.
2. 2.-Procedimiento de obtención de la perturbación geomagnética local (15) a latitudes medias. segun la reivindicación primera, caracterizado por que la contribución del interior de la Tierra (14) se determina mediante un valor promedio de los datos existentes en una base de datos del procedimiento de días clasificados previamente como en calma.
3. 3.-Procedimiento de obtención de la perturbación geomagnética local a latitudes medias, según la reivindicación primera o segunda, caracterizado por que la curva de variación diaria (13), usada para un día clasificado como perturbado, es la curva de variación diaria del día clasificada cama en calma mas próximo en el tiempo.
4. 4.-Procedimiento de obtención de la perturbación geomagnética local (15) a latitudes medias, según cualquiera de las reivindicaciones primera a tercera caracterizado por que la clasificación de dia en calma geomagnética (213) se realiza mediante una valoración de la bondad del ajuste del modelo de día en calma (212) a los datos de la componente horizontal de campo magnético medido (11), una vez sustraídos los datos de la curva auxiliar de tendencia de los valares de noche (12).