Acessibilidade / Reportar erro

Algoritmo para completar series de valores de la frecuencia crítica de la región F de la ionósfera sobre Concepción, Chile

Algorithm to fill in ionospheric F-region critical frequency series for Concepción, Chile

Resúmenes

Se usa la serie incompleta de valores horarios de la foF2 sobre Concepción (36.8º S; 73.0º O), observados en el intervalo 1958 - 1989, para determinar inicialmente una serie completa de valores horarios de foF2 que pueden considerarse representativos de condiciones geomagnéticas quietas, pero que toma en cuenta las dependencias de la foF2 con la estación del año y con el nivel de actividad solar. Luego se propone un algoritmo para calcular valores horarios faltantes en el intervalo, en que sí se toma en cuenta el nivel de actividad geomagnética observado en las horas correspondientes, usando un índice de actividad geomagnética ponderado. Como evidentemente no se puede determinar la bondad del algoritmo propuesto para las horas en que no ha sido posible medir la foF2, se compara los valores horarios calculados con los observados, para días en que no hay ningún valor faltante, y que se consideran representativos de diferentes estaciones del año y niveles de actividad solar y geomagnética. Se encuentra que las diferencias no son significativas a excepción de días de equinoccio y verano en condiciones geomagnéticas muy perturbadas.

Frecuencia crítica; Región F; Ionósfera


A complete series of hourly values of foF2 considered representative of geomagnetic quiet conditions is determined using the incomplete series of observed hourly values at Concepción (36.8º S; 73.0º W) in the 1 March 1958 - 1 July 1989 interval. This series takes into account known seasonal and solar activity level dependencies. Then a simple algorithm is proposed to calculate the missing hourly values of the original series, which takes into account the geomagnetic activity level for the corresponding hours, in terms of a weighted geomagnetic activity index. Since there is no way to assert the goodness of the proposed algorithm when there is no observed foF2 values, calculated and observed values are compared for a few days for which they do exist. Days considered representative of different seasons and levels of solar and geomagnetic activity are used. It is found that differences between calculated and observed values are not significant except for equinox and summer conditions at high geomagnetic activity level.

Critical frequency; F-region; Ionosphere


Algoritmo para completar series de valores de la frecuencia crítica de la región F de la ionósfera sobre Concepción, Chile

M.A. Arriagada1 & A. J. Foppiano2

1Departamento de Física, Universidad del Bío-Bío, Casilla 5-C Concepción,

Teléfono +56-41-261219

Fax +56-41-320068

marriga@ubiobio.cl

2Departamento de Física de la Atmósfera y del Océano, Universidad de Concepción,

Casilla 4009, Concepción, Chile,

Teléfono +56-41-203083

Fax +56-41-220104

foppiano@udec.cl

Se usa la serie incompleta de valores horarios de la foF2 sobre Concepción (36.8º S; 73.0º O), observados en el intervalo 1958 – 1989, para determinar inicialmente una serie completa de valores horarios de foF2 que pueden considerarse representativos de condiciones geomagnéticas quietas, pero que toma en cuenta las dependencias de la foF2 con la estación del año y con el nivel de actividad solar. Luego se propone un algoritmo para calcular valores horarios faltantes en el intervalo, en que sí se toma en cuenta el nivel de actividad geomagnética observado en las horas correspondientes, usando un índice de actividad geomagnética ponderado. Como evidentemente no se puede determinar la bondad del algoritmo propuesto para las horas en que no ha sido posible medir la foF2, se compara los valores horarios calculados con los observados, para días en que no hay ningún valor faltante, y que se consideran representativos de diferentes estaciones del año y niveles de actividad solar y geomagnética. Se encuentra que las diferencias no son significativas a excepción de días de equinoccio y verano en condiciones geomagnéticas muy perturbadas.

Palabras claves: Frecuencia crítica; Región F; Ionósfera

Algorithm to fill in ionospheric F-Region critical frequency series for Concepción, Chile - A complete series of hourly values of foF2 considered representative of geomagnetic quiet conditions is determined using the incomplete series of observed hourly values at Concepción (36.8º S; 73.0º W) in the 1 March 1958 – 1 July 1989 interval. This series takes into account known seasonal and solar activity level dependencies. Then a simple algorithm is proposed to calculate the missing hourly values of the original series, which takes into account the geomagnetic activity level for the corresponding hours, in terms of a weighted geomagnetic activity index. Since there is no way to assert the goodness of the proposed algorithm when there is no observed foF2 values, calculated and observed values are compared for a few days for which they do exist. Days considered representative of different seasons and levels of solar and geomagnetic activity are used. It is found that differences between calculated and observed values are not significant except for equinox and summer conditions at high geomagnetic activity level.

Key words: Critical frequency; F-region; Ionosphere.

INTRODUCCIÓN

En el análisis de la larga serie de valores de la frecuencia crítica de la región F, foF2, sobre Concepción (36.8º S; 73.0º O), invariablemente se presenta el problema de completar la serie por la falta de registros, debidos a fallas en los equipos, a errores humanos, o a condiciones ionosféricas que impiden la obtención de valores numéricos. Para completar la serie se puede hacer uso de diversos procedimientos de interpolación, pudiendo elegirse el más apropiado de entre ellos, dependiendo fundamentalmente del número de valores horarios faltantes en relación con las conocidas variaciones regulares, diurna, estacional y con el nivel de actividad solar. Estos procedimientos no toman normalmente en cuenta otras variaciones, de carácter episódico, como las que se observan en los registros ionosféricos, por ejemplo en asociación con la ocurrencia de llamaradas solares y de eyecciones de masa de la corona solar.

El objetivo del presente trabajo es proponer un algoritmo que permita determinar los valores horarios faltantes de la foF2, para cualquier hora del día en el intervalo 1958 – 1989, el que se fundamenta en la determinación empírica de las dependencias de la foF2 con los niveles de actividades solar y geomagnética. Se estima que este procedimiento da una mejor cuenta de las variaciones regulares y de las episódicas, además de permitir la completación de valores faltantes, incluso para varios meses consecutivos. Aunque el algoritmo es obviamente válido solo para Concepción durante el intervalo considerado, se podría sugerir su aplicación a otras localidades que cuenten con series de valores comparables.

METODOLOGÍA

La metodología usada sigue estrechamente la propuesta por Wrenn et al. (1987) y Rodger et al. (1989) para el estudio del efecto de tormentas geomagnéticas en la región F.

En primer lugar se determinó el índice ponderado de actividad geomagnética acumulada, ap(0.75), definido inicialmente por Wrenn (1987), para cada intervalo trihorario desde las 05:00 horas Tiempo Universal (TU) del 1 de marzo de 1958 (00:00 horas de Tiempo Local (TL) para el huso horario de 75° O) hasta las 04:00 horas TU del 1 de julio de 1989 (23:00 horas TL).

Luego se agruparon para cada hora los valores horarios de la foF2 de todos los días de un mes determinado, correspondientes al intervalo 1 de noviembre de 1968 a 30 de noviembre de 1979, para los cuales el nivel de actividad geomagnética es bajo (índice trihorario ap(0.75) < 6) y el nivel de actividad solar (flujo solar del día en 10.7 cm, F10.7) está en uno de cuatro rangos: 1 (50 d F10.7 < 100); 2 (100 d F10.7 < 150); 3 (150 d F10.7 < 200) y 4 (200 d F10.7 < 250). Este intervalo se eligió inicialmente para que los resultados del presente análisis sean directamente comparables con los resultados de Wreen et al., (1987). Seguidamente se calcularon valores medios mensuales para cada hora, usando los valores horarios de todos los días así agrupados para cada mes. En la Fig. 1 se muestran estos valores medios mensuales para el caso del mes de agosto, los que ilustran las dependencias empíricas de la foF2 con el nivel de actividad solar para el mes en cuestión. Se hace notar que, para algunas horas, la condición nivel de actividad geomágnética baja no se cumple en el rango 4 de nivel de actividad solar.

Figura 1
- Frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O). (···) Valores medios mensuales calculados usando los valores horarios de todos los días de los meses de agosto del intervalo 1 de noviembre de 1968 al 30 de noviembre de 1979, para los cuales el nivel de actividad geomagnética es bajo (índice trihorario ap(0.75) < 6) y el nivel de actividad solar (flujo solar del día en 10.7 cm) está en uno de cuatro rangos: 1 (50 d F10.7 < 100); 2 (100 d F10.7 < 150); 3 (150 d F10.7 < 200) y 4 (200 d F10.7 < 250). Para algunas horas las dos condiciones no se cumplen simultáneamente en el rango 4. (ooo) Valores mensuales representativos, foF2mr, para agosto de 1969 (F10.7 medio mensual 146.1) calculados por interpolación lineal entre los correspondientes a los niveles de actividad solar 2 y 3 (118.26 a 121.13, y 164.26 a 167.61, respectivamente, dependiendo de la hora del día). Las barras de variabilidad (± una desviación estándar) indican las variabilidades máximas encontradas entre las 00 y 06, y las 06 y 18 horas, respectivamente. A otras horas las variabilidades son similares.

Figure 1 - Ionospheric F-region critical frequency observed over Concepción (36.8º S; 73.0º W). (···) Monthly mean values determined using hourly values for all days from the 1 November 1968 – 30 November 1979 interval, corresponding geomagnetic activity low level (three-hourly index ap(0.75) < 6) and four solar activity levels (daily 10.7 cm solar flux): 1 (50 d F10.7 < 100); 2 (100 d F10.7 < 150); 3 (150 d F10.7 < 200) y 4 (200 d F10.7 < 250). Both conditions are not simultaneously met for some hours at 4 level. (ooo) Monthly values, foF2mr, for August 1969 determined by linear interpolation between corresponding values for solar activity levels 2 and 3 (118.26 to 121.13, and 164.26 to167.61, respectively, depending on time-of-day). Variability bars (± one standard deviation) indicate maxima variabilities found between 00 and 06, and 06 to 18 hours, respectively. Variabilities for other times-of-day are similar.


A continuación se calculó el valor medio mensual del F10.7 para cada mes del intervalo considerado. Se determinó enseguida, para cada hora del día, valores de la foF2 representativos para cada mes del intervalo, foF2mr, por interpolación lineal entre los valores de la foF2 correspondientes a dos niveles de actividad solar sucesivos del mes considerado. En la Fig. 1 se muestran los valores representativos así calculados para agosto de l969. El valor medio del F10.7 para agosto de 1969 es 146.1, que está entre los valores del F10.7 correspondientes a agosto nivel 2 (118.26 a 121.13, dependiendo de la hora del día) y a agosto nivel 3 (164.26 a 167.61).

El siguiente paso fue asignar esos valores representativos de cada mes al día mediano del mes, por ejemplo al día 16 en los meses de 31 días, o al "día 15.5" para los meses de 30 días. En la Fig. 2 se muestran los valores foF2mr correspondientes a julio, agosto y septiembre de 1969. Luego se determinaron valores horarios, foF2dr, para todos los días entre los días medianos de dos meses sucesivos, por interpolación lineal entre los valores foF2mr asignados a esos meses. En la Fig. 2 se muestra estos valores foF2dr calculados para el día 1 de agosto de 1969. En la Fig. 3 se muestran los valores foF2dr para los días 29 de julio al 4 de agosto de 1969, los que ilustran los pequeños cambios debidos a las dependencias con la estación del año y con el nivel de actividad solar. De la misma manera se determinó valores foF2dr para cada hora de cada día del intervalo considerado, que se pueden considerar como típicos de condiciones geomagnéticas quietas, para el nivel de actividad solar que le corresponde, pues las dependencias empíricas con el nivel de actividad solar ya han sido incluidas.

Figura 2
- Frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O) para nivel de actividad geomagnética bajo (índice trihorario ap(0.75) < 6). (···) Valores mensuales representativos, foF2mr, para julio (J), agosto (A) y septiembre (S) de 1969, calculados como los que se muestran en la Fig. 1 para el caso de agosto de 1969, y asignados al día mediano del mes correspondiente. (ooo) Valores horarios representativos, foF2dr, para el 1 de agosto de 1969 calculados por interpolación lineal entre los correspondientes a los días medianos del julio y agosto.

Figure 2 - Ionospheric F-region critical frequency observed over Concepción (36.8º S; 73.0º W) for low level geomagnetic activity (three-hourly index ap(0.75) < 6). (···) Monthly values, foF2mr, for July (J), August (A) and September (S) 1969, determined as those shown in Fig. 1 for August 1969 and assigned to the median day of corresponding month. (ooo) Hourly values, foF2dr, for 1 August 1969 determined using linear interpolation between those corresponding to July and August median days.



Figura 3 - Frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O) para nivel de actividad geomagnética bajo (índice trihorario ap(0.75) < 6). Valores horarios representativos, foF2

dr, para los días 29 de julio al 4 de agosto de 1969, calculados como los que se muestran en la

Fig. 2 para el caso del 1 de agosto de 1969.

Figure 3 - Ionospheric F-region critical frequency observed over Concepción (36.8º S; 73.0º W) for low level geomagnetic activity (three-hourly index ap(0.75) < 6). Hourly values, foF2dr, for 29 July to 4 August 1969, determined as those shown in Fig. 2 for the 1 August 1969 case.

Seguidamente se determinó la dependencia de la foF2 con el nivel de actividad geomagnética. Para ello se calculó el cuociente, q, entre el valor de la foF2 observada y la foF2 representativa de condiciones geomagnéticas quietas, foF2dr, para cada hora de todo el intervalo inicialmente elegido (1 de noviembre de 1968 a 30 de noviembre de 1979), y se agrupó estos cuocientes de acuerdo con la hora del día, mes del año y cuatro rangos de actividad geomagnética seleccionados por el valor del ap(0.75) (muy quieto: 0 d ap(0.75) < 6; quieto: 6 d ap(0.75) < 18; perturbado: 18 d ap(0.75) < 30; y muy perturbado: 30 d ap(0.75)). En la Fig. 4a se muestran los valores de q así obtenidos para el mes de agosto, lo que ilustra las dependencias empíricas de la foF2 con el nivel de actividad geomagnética para el mes en cuestión. Se advierte que las dependencias para los distintos meses de una misma estación del año son muy similares, por lo que se calcula valores medios de q para equinoccio (marzo, abril, septiembre y octubre), invierno (mayo, junio, julio y agosto), y verano (noviembre, diciembre, enero y febrero). La Fig. 4b ilustra las dependencias para invierno. Estos resultados son consistentes con los publicados para un intervalo equivalente por Wrenn et al. (1987) y Rodger et al. (1989), correspondientes a otro lugar del hemisferio Sur, y por Field & Rishbeth (1997) para 53 lugares distribuidos en ambos hemisferios, incluyendo Concepción.

Figura 4
- Cuociente entre frecuencia crítica de la región F de la ionósfera observada sobre Concepción (36.8º S; 73.0º O) y la frecuencia crítica determinada para condiciones de nivel de actividad geomagnética bajo (índice trihorario ap(0.75) < 6), para el intervalo 1 de noviembre de 1968 al 30 de noviembre de 1979 (ver texto). ---- Muy quieto (0 d ap(0.75) < 6); LLL Quieto (6 d ap(0.75) < 18); ———— Perturbado (18 d ap(0.75) < 30); y - - - - Muy perturbado (30 d ap (0.75)). (a) Agosto. (b) Invierno.

Figure 4 - Ratio between ionospheric F-region critical frequency observed over Concepción (36.8º S; 73.0º O) and critical frequency determined for low geomagnetic activity level (three-hourly index ap(0.75) < 6) corresponding to 1 November 1968 – 30 November 1979 (see text). ---- very quiet (0 d ap(0.75) < 6); LLL quiet (6 d ap(0.75) < 18); ———— perturbed (18 d ap(0.75) < 30); and - - - - very perturbed (30 d ap (0.75). (a) August. (b) Winter.


Así, finalmente, el valor faltante de la foF2 para cualquier hora de un día determinado se puede calcular mediante la expresión:

foF2 = q foF2dr , (1)

en que el cuociente q se elige de la Tab. 1 que lo lista según hora del día, estación del año y el rango de ap(0.75) correspondiente, y foF2dr es el valor representativo de condiciones geomagnéticas quietas para la hora en cuestión (en la Fig. 4b sólo se muestran los valores de q para invierno). Cabe hacer notar que, aunque para una hora determinada las tablas listan el mismo q para los cuatro meses de una misma estación del año y un nivel de actividad geomagnética fijo, los valores de la foF2 para esos meses no son los mismos porque las correspondientes foF2dr, no son iguales (pequeñas diferencias por dependencias con la estación del año y con el nivel de actividad solar), y porque los q correspondientes también son en general diferentes (por la presencia o no de actividad magnética episódica a la respectiva hora).

Tabla 1
- Cuociente entre la frecuencia crítica de la región F de la ionósfera observada sobre Concepción (36.8º S; 73.0º O) y la frecuencia crítica determinada para condiciones de nivel de actividad geomagnética bajo (índice trihorario ap(0.75) < 6), correspondientes al intervalo 1 de noviembre de 1968 al 30 de noviembre de 1979, para distintas horas del día, estaciones del año y niveles de actividad geomagnética: Muy quieto (0 d ap(0.75) < 6); Quieto (6 d ap(0.75) < 18); Perturbado (18 d ap(0.75) < 30); y Muy perturbado (30 d ap (0.75).

Table 1 - Ratio between ionospheric F-region critical frequency observed over Concepción (36.8º S; 73.0º O) and critical frequency determined for low geomagnetic activity level (three-hourly index ap(0.75) < 6) corresponding to 1 November 1968 – 30 November 1979, for different times-of-day, seasons and geomagnetic activity levels: very quiet (0 d ap(0.75) < 6); quiet (6 d ap(0.75) < 18); perturbed (18 d ap(0.75) < 30); and very perturbed (30 d ap (0.75).


Como se mencionó anteriormente, en contadas ocasiones dentro del intervalo inicial considerado no existen valores de la foF2 que satisfagan simultáneamente las condiciones de actividad magnética quieta y actividad solar muy alta. En estos casos no es posible determinar valores de la foF2mr, por cuanto no es recomendable una extrapolación cuando sólo se tienen valores para los tres niveles de actividad solar mas bajos. En tal situación los valores foF2mr se determinan mediante un proceso de interpolación, considerando dos procedimientos, según la disponibilidad de datos. El primero consiste en calcular el valor foF2mr por interpolación mediante una función spline cúbica, usando los valores correspondientes a los dos meses consecutivamente anteriores y a los de los dos meses consecutivamente posteriores. El segundo se utiliza cuando también faltan valores foF2mr en los meses anteriores o en los posteriores, y consiste en una interpolación con una función spline cúbica con valores de la foF2mr, de meses del mismo nombre al del que falta, dentro del intervalo inicial, cuyos niveles de actividad solar sean similares al del mes faltante.

Todo lo anterior se hizo separadamente también para los intervalos aproximadamente undecenales del 1 de marzo de 1958 al 31 de octubre de 1968 y del 1 de diciembre de 1979 al 30 de junio de 1989. No se hizo el análisis para un intervalo único (1 marzo 1958 a 30 junio 1989) porque la dependencia de la foF2 con el nivel de actividad solar, pareciera ser diferente para cada uno de los tres intervalos considerados (Foppiano & Arriagada, 1993), como lo han sugerido otros autores para otros lugares e intervalos (e.g. Smith & King, 1981), aun cuando no hay consenso sobre su significancia (Bradley, 1994). La dependencia de la foF2 con el nivel de actividad geomagnética también pareciera diferir de un intervalo a otro.

RESULTADOS

Haciendo uso del algoritmo propuesto se determinaron los valores de la foF2 para todas las ocasiones en que no había registros, lográndose así una serie completa para el intervalo 1 de Marzo de 1958 al 30 de Junio de 1989. Cabe destacar que en una de esas ocasiones faltaban registros para tres meses completos (septiembre, octubre y noviembre de 1975) por falla de los equipos.

Como no es posible estimar la bondad del algoritmo para ocasiones en que no hay registros, se presentan ejemplos de días representativos de las distintas estaciones año y de los cuatro niveles de actividad geomagnética. Se eligieron aquellos días para los cuales existe el mayor número posible de registros.

En las Figs. 5 al 8 se muestran los valores de la foF2 calculadas con el algoritmo y los valores observados para todos los días seleccionados. En general, hay una buena predicción para los días de invierno en todos los niveles de actividad geomagnética. Lo mismo ocurre para los días de equinoccio y de verano, en condiciones geomagnéticas muy quietas y quietas, aunque en estos casos es posible observar diferencias no despreciables.


Figura 5

- Comparación entre los valores observados y los valores calculados con el algoritmo propuesto para la frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O). (···) Valores observados. (ooo) Valores calculados. Año 1974. Nivel de actividad geomagnética muy quieto (0 d ap(0.75) < 6). (a) Invierno (7 de junio). (b) Equinoccio (15 de abril). (c) Verano (30 de noviembre).

Figure 5 - Comparison between observed and calculated values using the proposed algorithm for the ionospheric F-region critical frequency over Concepción (36.8º S; 73.0º W). (···) Observed values. (ooo) Calculated values. 1974. Very quiet geomagnetic activity level (0 d ap(0.75) < 6). (a) Winter (7 July). (b) Equinox (15 April). (c) Summer (30 November).


Figura 6

- Comparación entre los valores observados y los valores calculados con el algoritmo propuesto para la frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O). (···) Valores observados. (ooo) Valores calculados. Año 1978. Nivel de actividad geomagnética quieto (6 d ap(0.75) < 18). (a) Invierno (1 de julio). (b) Equinoccio (1 de octubre). (c) Verano (16 de diciembre).

Figure 6 - Comparison between observed and calculated values using the proposed algorithm for the ionospheric F-region critical frequency over Concepción (36.8º S; 73.0º W). (···) Observed values. (ooo) Calculated values. 1978. Quiet geomagnetic activity level (6 d ap(0.75) < 18). (a) Winter (1 July). (b) Equinox (1 October). (c) Summer (16 December).


Figura 7

- Comparación entre los valores observados y los valores calculados con el algoritmo propuesto para la frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O). (···) Valores observados. (ooo) Valores calculados. Año 1974. Nivel de actividad geomagnética perturbado (18 d ap(0.75) < 30). (a) Invierno (26 de julio). (b) Equinoccio (1 de octubre). (c) Verano (20 de diciembre).

Figure 7 - Comparison between observed and calculated values using the proposed algorithm for the ionospheric F-region critical frequency over Concepción (36.8º S; 73.0º W). (···) Observed values. (ooo) Calculated values. 1974. Perturbed geomagnetic activity level (18 d ap(0.75) < 30). (a) Winter (26 July). (b) Equinox (1 October). (c) Summer (20 December).

Figura 8

- Comparación entre los valores observados y los valores calculados con el algoritmo propuesto para la frecuencia crítica de la región F de la ionósfera sobre Concepción (36.8º S; 73.0º O). (···) Valores observados. (ooo) Valores calculados. Año 1974. Nivel de actividad geomagnética muy perturbado (30 d ap(0.75)). (a) Invierno (27 de junio). (b) Equinoccio (16 de octubre). (c) Verano (12 de noviembre).

Figure 8 - Comparison between observed and calculated values using the proposed algorithm for the ionospheric F-region critical frequency over Concepción (36.8º S; 73.0º W). (···) Observed values. (ooo) Calculated values. 1974. Very perturbed geomagnetic activity level (30 d ap(0.75)). (a) Winter (27 June). (b) Equinox (16 October). (c) Summer (12 December).

La predicción desmejora considerablemente en los días de equinoccio y verano, en condiciones geomagnéticas perturbadas y muy perturbadas, donde para la mayoría de las horas del día se aprecian diferencias significativas.

Probablemente una parte de las diferencias observadas se deba a que no se ha tomado en cuenta que la respuesta de la ionósfera a cambios del nivel de actividad geomagnética puede depender del nivel de actividad solar. Por definición, el cuociente q es independiente del nivel de actividad solar.

CONCLUSIONES

Se ha logrado obtener dos series completas de valores de la frecuencia crítica de la región F sobre Concepción (36.8º S; 73.0º O), para el intervalo 1958 – 1989. Una se puede considerar representativa de condiciones geomagnéticas quietas, y por consiguiente puede ser útil para determinar tendencias interdecadales que no incluyan los efectos de la actividad geomagnética. La otra puede ser de valor tanto para el estudio de efectos de carácter episódicos, como para construir una base para la preparación de modelos empíricos.

Se encuentra que, en general, las diferencias entre los valores de la frecuencia crítica calculados con el algoritmo propuesto y los valores observados, cuando estos existen, son mínimas, a excepción de los días de equinoccio y verano en condiciones geomagnéticas perturbadas y muy perturbadas.

AGRADECIMIENTOS

Se agradece los comentarios recibidos de tres revisores, los que condujeron a una significativa revisión del texto original. Este trabajo ha sido patrocinado por el Fondo Nacional de Desarrollo Científico y Tecnológico (Proyecto FONDECYT N° 1960628).

Manuscript submitted May 5, 1998

Revised version accepted July 20, 1999

Algorithm to fill in ionospheric F-Region critical frequency series for Concepción, Chile

A complete series of foF2 hourly values considered representative of geomagnetic quiet conditions is determined using the incomplete series of observed hourly values at Concepción (36.8º S; 73.0º W) in the 1 March 1958 - 1 July 1989 interval. This is carried out in four steps, closely following procedures published by Wrenn et al. (1987) and Rodger et al. (1989) as follows. foF2 hourly values, corresponding to low geomagnetic activity (as measured by a weighted geomagnetic activity index), for the 1 November 1968 - 30 November 1979 sub-interval are first selected. These are grouped according to time-of-day, month and four solar activity levels (as measured by solar flux in the 10.7 cm band) independent of year, and monthly mean values are determined (Fig. 1). Then representative monthly values, foF2mr, are determined for each month of each year by interpolation between monthly mean values corresponding to the two nearest solar activity levels using the monthly mean solar activity for the given month (Fig. 1). Assigning these foF2mr,values to the median day of the month, representative hourly values, foF2dr, for each day of the sub-interval are calculated by interpolation (Fig. 2). The complete foF2dr hourly series thus obtained for the whole sub-interval takes into account known seasonal and solar activity level dependencies (Fig. 3).

Computing the ratio, q, between observed foF2 and the corresponding foF2dr then permits the dependency of foF2 hourly values on geomagnetic activity to be determined. Values of q are grouped according to time-of-day, month and four levels of geomagnetic activity (again as measured by a weighted geomagnetic activity index) independent of year, and monthly mean values are determined (Fig. 4a). These are a measure of the foF2 dependency on geomagnetic activity. In practice, q monthly mean values for months corresponding to the same season are very similar, and thus only seasonally mean values are needed (Fig. 4b and Table 1).

The proposed algorithm to compute the sub-interval foF2 missing values uses the simple expression

foF2 = q foF2dr, (1)

where values of q are selected from Table 1 and foF2dr are those referred to as the complete series of hourly values of foF2 considered representative of geomagnetic quiet conditions already obtained.

The same procedures indicated above are used independently for the other approximately undecenal sub-intervals 1 March 1958 - 31 October 1968 and 1 December 1979 - 30 June 1989 (the corresponding q tables are not reproduced here).

Since there is no way to assert the goodness of the proposed algorithm when there is no observed foF2 values, calculated and observed values are compared for a few days for which they do exist. Days considered representative of different seasons and levels of solar and geomagnetic activity are used (Figs. 5 to 8).

Summarizing: (a) Two complete series of foF2 hourly values have been obtained spanning 30 years. One can be considered representative of quiet geomagnetic conditions and thus can be useful to deriving long term trends that are independent of geomagnetic activity. The other can be of value both to study event effects and to serve as a base for empirical modelling; (b) Differences between calculated and observed values are not significant except for equinox and summer conditions at high geomagnetic activity level.

NOTAS SOBRE LOS AUTORES

A.J. Foppiano.

Profesor de Estado en Matemáticas y Física, Universidad de Chile y Doctor (Ph.D), Imperial College of Science and Technology, University of London. Se desempeña actualmente como Profesor Titular del Departamento de Física de la Atmósfera y del Océano de la Universidad de Concepción, Chile. Su área de interés es la Física de la Ionosfera, en particular absorción auroral de ondas radioeléctricas y morfología y dinámica de la región F.

M.A. Arriagada.

Profesor de Estado en Física, Universidad de Chile y Magister en Física, Universidad Austral de Chile. Se desempeña actualmente como Profesor Asociado del Departamento de Física de la Universidad del Bío-Bío, Chile. Su área de interés es la Física de la Ionosfera, en particular morfología y dinámica de la región F.

  • BRADLEY, P.A. -1994- Further study of foF2 and M(3000) F2 in different solar cycles. Annali di Geofisica, 37:201-208.
  • FIELD, R.P. & RISHBETH, H. -1997- The response of the ionospheric F2-layer to geomagnetic activity: an analysis of worldwide data. J. Atmos. Sol. Terr. Phys. 59:163-180.
  • FOPPIANO, A.J. & ARRIAGADA, M.A. -1993- Efecto del nivel de actividad geomagnética sobre la concentración electrónica máxima de la ionósfera sobre Concepción, para observaciones correspondientes a tres ciclos solares. Tercera Conferencia Latinoamericana de Geofísica Espacial, 1-5 de Noviembre, Habana, Cuba.
  • RODGER, A.S., WRENN, G.L. & RISHBETH, H. -1989- Geomagnetic storms in the Antarctic F-region. II Physical interpretation. J. Atmos. Terr. Phys. 51:851-866.
  • SMITH, P.A. & KING, J.W. -1981- Long-Term relationships between sunspots, solar faculae and the ionosphere. J. Atmos. Terr. Phys. 43:1057-1063.
  • WRENN, G.L. -1987- Time-weighted accumulations ap (t) and Kp(t). J. Geophys. Res. 92:10125-10129.
  • WRENN, G.L., RODGER, A.S. & RISHBETH, H. -1987- Geomagnetic storms in the Antarctic F-region I. Diurnal and seasonal patterns for main phase effects. J. Atmos. Terr. Phys. 49:901-913.

Fechas de Publicación

  • Publicación en esta colección
    18 Mayo 2000
  • Fecha del número
    Mar 1999

Histórico

  • Recibido
    05 Mayo 1998
  • Acepto
    20 Jul 1999
Sociedade Brasileira de Geofísica Av. Rio Branco, 156, sala 2510, 20043-900 Rio de Janeiro RJ - Brazil, Tel. / Fax: (55 21) 2533-0064 - São Paulo - SP - Brazil
E-mail: sbgf@sbgf.org.br