INTRODUCCIÓN
Como es bien conocido, las corrientes ionosféricas de la capa E son las fuentes de las variaciones diarias de los días calmos Sq, estas variaciones toman valores del orden de la decena de nT [1]. Además de ser estudiadas para intentar comprender el comportamiento ionosférico, las variaciones geomagnéticas de los días calmos son clave a la hora de realizar investigaciones vinculadas con disciplinas afines, por ejemplo, con la actividad de la radiación solar [2]. El sistema de corrientes ionosférica no se puede determinar a partir del análisis de datos de estaciones geomagnéticas solamente. Debido a esto se usa definir un sistema de corrientes bidimensional equivalente que fluye en una capa esférica delgada situada supuestamente a 110 km de altura. Una forma de evaluar las corrientes equivalentes es suponer una hoja de corriente horizontal en la ionósfera [3]. La ley de Biot- Savart permite calcular el campo ∆B [T] que produce una hoja de corriente de intensidad I [A/m]. Este toma la forma ∆B =µ 0I/4π. Cuando el campo magnético se da en nT, la expresión anterior puede reescribirse tomando el siguiente aspecto: ∆B = 2πK/10, donde K es la densidad de corriente en mA/m [1]. Para estudiar los sistemas de corrientes equivalentes superficiales ionosféricas hay que tener presente que en los registros de los días calmos se tendrán las contribuciones magnetosféricas, ionosféricas e inducidas en la corteza por las variaciones anteriormente mencionadas. Distintos autores proponen diferentes tipos de correcciones a los datos registrados cuando analizan las variaciones Sq con el fin de eliminar las contribuciones debidas a efectos magnetosféricos, como son las propuestas de Takeda [4] y de Choudhary [5]. Sin embargo, como se mostrará más adelante en este trabajo, ambos métodos presentan desventajas por no considerar la variabilidad del índice Dst.
Por otra parte, los índices para la corriente anillo SMR [6-8] obtenidos desde https://supermag.jhuapl.edu/, son elaborados a partir de los registros de 98 magnetómetros ubicados en latitudes medias y bajas.
Además, al momento de evaluar las corrientes superficiales equivalentes, la contribución de las corrientes corticales inducidas se consideran de un tercio de la variación Sq horizontal medida. De esta manera, las densidades superficiales de corriente equivalente norte y este, KN y KE sobre una estación [3], son:
donde ∆X, ∆Y se calculan restando al dato registrado y corregido, el valor de referencia nocturno. Si bien este modelo tiene sus limitaciones espaciales, algunos autores lo consideran adecuado cuando se analizan los datos en una única estación [3].
En el presente estudio se emplean los registros de la red de estaciones magnéticas permanentes de la República Argentina (ROMP). Esta red está conformada por las estaciones de PIL (Pilar, -31.4◦, 296.12◦), LAS (Las Acacias, -35.◦006, 302.◦31), TRW (Trelew, -43.◦267, 294.◦617) y ORC (Orcadas, -60.◦737, 315.◦26). La ubicación de las mismas se reprenta esquemáticamente en el mapa (Fig. 1). Debido a que en Las Acacias se registra solamente el módulo del campo, esta estación ha quedado excluida del presente análisis. Los datos relativos a las estaciones geomagnéticas fueron tomados de la red Intermagnet (www.intermagnet.org) en el intervalo de tiempo que abarca desde mayo de 2019 hasta mayo de 2021.
En este trabajo se procede por primera vez a producir la corrección magnetosférica mediante el empleo de los índices para la corriente anillo SMR. En la sección 1 se discute la selección de los días de estudio en el intervalo 2019-2021, en la Sec. 2 se realizan las correcciones magnetosféricas por los métodos de Takeda [4] y de Choudhary [5], posteriormente se discuten las desventajas que presentan ambos métodos. Se presenta la metodología de aplicación con los índices SMR y resultando en la remoción eficiente de los efectos magnetosféricos por la corriente anillo. En la Sec. 3, finalmente, se calculan las corrientes ionosféricas superficiales equivalentes por el método de Stening [3] y se describen las características salientes de las componentes KE, KN, corrientes superficiales este y norte, respectivamente.
MÉTODOS
Elección de los días de estudio
Para este estudio se dividieron los registros en tres grupos, siguiendo a Yamazaki y Maute [1]: “verano” (noviembre- febrero), “equinoccios” (marzo-abril-septiembre-octubre) e “invierno” (mayoagosto). De entre los cinco días internacionalmente calmos, se seleccionaron los de índice Kp menores o iguales a 2. De este modo se tomaron los días 19.10.19, 19.04.20 y 09.09.20 para el grupo “equinoccios”, 16.11.20 y 04.12.20 para el grupo “verano”, 9.05.20 y 5.05.21 para el grupo “invierno”. Este último grupo se descartó del estudio comparativo por encontrarse incompletos los registros en PIL y TRW. Tanto los días internacionalmente calmos (QDays por su denominación en inglés) como los valores del índice Kp fueron obtenidos de la página del World Data Center C de Japón.
Correcciones magnetosféricas
Al discutir las contribuciones de las corrientes magnetosféricas se suele tener en consideración principalmente las correspondientes a la corriente anillo y al sistema de la cola de la magnetósfera. La contribución perteneciente a la corriente de Chapman-Ferraro suele despreciarse pensando en que se encuentra muy alejada del punto de observación.
Por otra parte, este análisis es válido en estaciones geomagnéticas que se encuentren fuera de las regiones aurorales y del ecuador magnético, debido a que, como es bien conocido, el sistema de corrientes equivalentes en las regiones aurorales tiene un comportamiento peculiar, incluyendo los aspectos transitorios de reconexión magnética con el sistema de corrientes de la cola de la magnetósfera. En contraste, las estaciones del ecuador magnético se ven afectadas por el sistema de corrientes del electrojet ecuatorial cuya dependencia del sistema de Sq se encuentra en discusión actualmente [1].
En el presente trabajo se adoptó el criterio de realizar esta corrección solamente cuando la correlación entre el índice Dst y los promedios horarios de los datos registrados fuera alta. El intervalo que se usó para estimar el coeficiente de correlación por el método de Pearson es de 24 horas para cada día seleccionado. En la Tabla 1 se muestran los valores del coeficiente de correlación anteriormente mencionado para las fechas seleccionadas y las estaciones. Se encuentran escritas en rojo los coeficientes de las componentes que se deben corregir.
Método de Takeda
Distintos autores proponen distintos tipos de correcciones para la corriente anillo. Una corrección propuesta por Takeda [4] consiste en simular los efectos de la corriente anillo sobre la componente horizontal H mediante el índice Dst, resultando un ∆H por efectos magnetosféricos igual al valor del índice Dst multiplicado por el coseno de la latitud geomagnética de la estación (L), y estudiar el valor corregido utilizando
Se procedió al análisis de la corrección propuesta por Takeda. Para la misma se obtuvieron los valores del índice Dst desde la página Web del World Data Center C http://wdc.kugi.kyotou.ac.jp/dst_realtime/202001/index.html. Se debe recordar que el índice Dst evalúa el comportamiento global de la corriente anillo y está calculado con los datos de cuatro estaciones geomagnéticas.
El valor del coeficiente de correlación entre los datos de la componente horizontal en PIL y el índice Dst toma el valor -0,73. Pudo observarse que esta corrección no resulta eficiente en remover los efectos magnetosféricos, esta afirmación se fundamenta en el hecho de que el coeficiente de correlación entre los valores del índice Dst y los valores corregidos permanece alto, con valor -0.78. Los valores corregidos de la componente H por el método de Takeda se indican como Hmodelo1 (ver Fig. 2).
Método de Choudhary
Otro método, propuesto por Choudhary et al. [5], consiste en estimar la corrección mediante el ajuste lineal entre los promedios horarios de H y el índice Dst. En el caso de la componente horizontal en PIL la relación lineal entre ésta y Dst viene dada por
La corrección realizada con este método remueve eficientemente la contribución de la corriente anillo, el coeficiente de correlación entre los valores del índice Dst y los valores corregidos es cercana a cero (-0,009). Los valores corregidos de la componente H por el método de Choudhary se indican como Hmodelo2 (ver Fig. 2).
Método utilizando el índice SMR
Con la intención de evitar la introducción de variaciones ficticias, se recurrió al análisis de los índices de la corriente anillo SMR [7, 8]. Los valores del índice SMR fueron obtenidos desde la red SuperMAG.
En la Fig. 3 se muestra las series de los índices SMR y Dst correspondientes al día 4.12.20. El valor del coeficiente de correlación entre la componente horizontal y SMR toma el valor -0,55. Se realiza la corrección mediante el ajuste lineal
El valor del coeficiente de correlación entre los datos corregidos y el índice SMR toma el valor -0,003 en consecuencia, se produce la remoción de los efectos de la corriente anillo en forma eficiente.
A forma de comparación, se calculan los valores de la componente H corregidos mediante SMR para PIL y se indican como Hmodelo3 (ver Fig. 2).
Comparación de las correcciones
En la Fig. 2 se muestran los valores de la componente H en PIL corregida por los distintos métodos. También se grafica el modelo que contiene la corrección por el método de Takeda, teniendo en cuenta las observaciones comentadas en la subsección correspondiente. Para comparar los modelos del ajuste lineal con Dst (modelo2) y del ajuste lineal con SMR (modelo3), se agregó una corrección a los mismos por una pequeña variación no cíclica (Hmodelo2∗ y Hmodelo3∗).
Tanto el índice Dst como el índice SMR son susceptibles de contener contribuciones debidas al sistema de corrientes de la cola de la magnetósfera [1]. En el caso del índice Dst, tres de las estaciones que contribuyen con sus datos podrían percibir la influencia de la corriente de la cola de la magnetósfera en el intervalo entre las 3 y las 9 hs TU. Si bien el análisis de esta influencia está fuera del objetivo del presente trabajo, un estudio somero arroja que el valor del coeficiente de correlación para este intervalo entre SMR y PCS (índice de actividad polar sur, calculado a partir de los registros de la estación Vostok, situada en el continente antártico) es de 0,6 mientras que entre Dst y AE (índice auroral calculado a partir de los registros de estaciones del hemisferio norte) es 0,9. Por este motivo se estima que existe la posibilidad de observar una contribución no ionosférica y que la misma resultaría de mayor intensidad cuando se corrige mediante el índice Dst que cuando se emplea SMR. Los índices que se han utilizado para contrastar son diferentes (PCS, AE) debido a su discretización. Los índices PCS y AE pueden descargarse desde la página web del Servicio Internacional de índices http://isgi.unistra.fr/. En consecuencia, se observa que por el método de Choudhary se introduciría un aumento ficticio de las variaciones del orden de la decena de nT durante cuatro horas, entre las 5 y las 9 hs TU, con un máximo de 15 nT a las 8 hs TU. También se analizó la posibilidad de introducir variaciones ficticias mediante este procedimiento con SMR, estimando las mismas en el orden de 7 nT, entre las 5 y las 9 hs TU.
En lo que sigue de este trabajo todas las correcciones a las componentes que se indicaron en la Tabla 1 se han producido por el método del ajuste lineal con SMR.
EVALUACIÓN DEL SISTEMA DE CORRIENTES EQUIVALENTES CON LOS DATOS DE UNA ESTACIÓN
Posteriormente a la realización de las correcciones, se estimaron los valores de las componentes norte y este de las densidades de corriente superficiales ionosféricas sobre las estaciones en los días seleccionados.
Como se mencionó anteriormente, el sistema de corrientes ionosférica no se puede determinar a partir del análisis de datos de estaciones geomagnéticas solamente, se usa definir un sistema de corrientes bidimensional equivalente que fluye en una capa esférica delgada situada supuestamente a 110 km de altura, de acuerdo con Yamazaki y Maute. Una forma de evaluar las corrientes equivalentes es suponer una hoja de corriente horizontal en la ionósfera siguiendo a Stening, donde la contribución de las corrientes corticales inducidas se consideran de un tercio de la variación Sq horizontal medida.
De esta manera, las densidades superficiales de corriente equivalente norte y este (KN y KE) sobre una estación vienen dadas por las Ecs. 1 y 2.
Se realizaron las estimaciones de las densidades superficiales de las corrientes equivalentes norte y este (KN y KE) para los ejemplos de estudio sobre la estaciones geomagnéticas, de ellos surge el contenido de las Tablas 2-4.
Por otra parte, para el horario de las 14 hs tiempo local TL, (17 hs TU), los valores máximos y mínimos de KN y KE resultaron como se muestran en la Tabla 5, como puede verse en las Figs. 4-6.
Comparaciones entre las estaciones por fecha
19.10.19: KNMAX y KEMAX tienen lugar en TRW, mientras que KNMIN y KEMIN en ORC. KNMIN de ORC se adelanta al de TRW mientras que KEMIN de TRW se adelanta al de ORC (Fig. 7a y 7b).
09.09.20: Los valores de KNMAX coinciden en magnitud. KNMIN tiene lugar en PIL, KEMAX en PIL Y KEMIN en ORC (Fig. 8a y 8b).
16.11.20 y 04.12.20: KNMAX y KNMIN tienen lugar en ORC. KEMAX tiene lugar en PIL y KEMIN en ORC. KNMAX y KNMIN de ORC adelanta al de PIL. Ver las Figs. 9a y 9b para el 16.11.20 y Figs. 10a y 10b para el 04.12.20.
En un caso especial en el que pueden compararse los valores calculados para las tres estaciones, 19.04.20: KNMAX y KNMIN tienen lugar en TRW. KEMAX tiene lugar en PIL y KEMIN en ORC. Además, KNMAX ORC adelanta a KNMAX de PIL y de TRW (Figs. 11a y 11b).
Evolución del sistema de corrientes superficiales equivalentes
Puede apreciarse mayor intensidad en el intervalo de tiempo entre las 12 y las 24 hs TU. Para PIL la evolución es horaria, para ORC es antihoraria. Para TRW puede verse que de 12 - 15 hs TU sigue a ORC (evolución antihoraria) y en adelante mayoritariamente sigue a PIL (evolución horaria), ver Fig. 12. Este comportamiento podría deberse al desplazamiento relativo de TRW respecto del foco del sistema de corrientes superficiales equivalentes del hemisferio sur.
Esta suposición puede confirmarse calculando la posición del foco del sistema de corrientes equivalentes, siguiendo los pasos aplicados por Torta [9]. Primero se determina el tiempo t 0 , cuando la variación de D (o Y) pasa por cero en cada observatorio que en el caso del hemisferio sur corresponde a cuando pasa de negativo a positivo (ver Stening [10]). Posteriormente se grafica el valor de la variación de H (o de X) evaluada en el tiempo t 0 en función de la latitud. La posición del foco se obtiene a partir de un ajuste por mínimos cuadrados, siendo el valor de la latitud correspondiente donde la línea cruza el cero. En este caso, aunque se cuenta con poca cantidad de estaciones [9], la latitud del foco resulta -41◦55’ 26.4”. De esta manera puede interpretarse que mientras desde el punto de vista de TRW para el mes de abril, mientras el foco se encuentra al noreste, TRW presenta un comportamiento similar a ORC. Luego, cuando el foco se encuentra al suroeste de TRW, presenta un comportamiento similar al de PIL.
CONCLUSIONES
En este trabajo se estudia el sistema de corrientes superficiales equivalentes sobre las estaciones de la red ROMP (PIL, TRW y ORC) entre 2019 y 2021. Se remueve la contribución correspondiente a la corriente anillo mediante el ajuste lineal entre el índice SMR y los datos registrados en la estación. Se estiman los valores de las densidades superficiales de corrientes equivalentes componentes este y norte sobre cada estación siguiendo el método de Stening. Se encuentra que para las 14 hs tiempo local los valores de KNMAX y KEMAX tienen lugar en PIL para el 04.12.20, mientras que los de KNMIN y KEMIN tienen lugar en ORC, el 19.10.19 y el 16.11.20 respectivamente.
Por otra parte, los máximos y mínimos absolutos se producen así: KNMAX y KNMIN en TRW el 19.04.20, mientras que KEMAX tiene lugar en PIL y KEMIN en ORC, el 04.12.20 y 16.11.20 respectivamente.
En relación con la evolución, queda en evidencia el desplazamiento relativo de TRW y el foco del sistema de corrientes superficiales equivalentes del hemisferio sur.