Solar activity variations of equatorial spread F occurrence and sustenance during different seasons over Indian longitudes: Empirical model and causative mechanisms

A comprehensive analysis using nearly two decades of ionosonde data is carried out on the seasonal and solar cycle variations of Equatorial Spread F (ESF) irregularities over magnetic equatorial location Trivandrum (8.5°N, 77°E). The corresponding Rayleigh Taylor (RT) instability growth rates (γ) are also estimated. A seasonal pattern of ESF occurrence and the corresponding γ is established for low solar (LSA), medium solar (MSA) and high solar (HSA) activity periods. For LSA, it is seen that the γ maximizes during post sunset time with comparable magnitudes for autumnal equinox (AE), vernal equinox (VE) and winter solstice (WS), while for summer solstice (SS) it maximizes in the post-midnight period. Concurrent responses are seen in the ESF occurrence pattern. For MSA, γ maximizes during post-sunset for VE followed by WS and AE while SS maximises during post-midnight period. The ESF occurrence for MSA is highest for VE (80%), followed by AE (70%), WS (60%) and SS (50%). In case of HSA, maximum γ occurs for VE followed by AE, WS and SS. The concurrent ESF occurrence maximizes for VE and AE (90%), WS and SS at 70%.The solar cycle variation of γ is examined. γ shows a linear variation with F10.7?cm flux. Further, ESF percentage occurrence and duration show an exponential and linear variation respectively with γ. An empirical model on the solar activity dependence of ESF occurrence and sustenance time over Indian longitudes is arrived at using the database spanning two solar cycles for the first time.     

F3 layer feature under low and medium solar activity observed at a Chinese low latitude station Fuke

Ionogram observations from the ionosonde at Fuke (9.5°N geomagnetic latitude), a Chinese low latitude station, in 2010–2012 are analyzed to present the features of F3 layer under low and moderate solar activity conditions. Structure of the ionogram, displaying the F3 layer, was more distinct and clear during MSA than LSA periods especially during spring to summer. Start time of occurrence of the F3 layer is about at 0830–0900 LT and is approximately the same for LSA and MSA conditions. The average duration time of the F3 layer occurrence was 181 min per day under F_10.7 = 75 condition, 263 min in F_10.7 = 99 and 358 min in F_10.7 = 125, respectively. The differences of h′F2 and h′F3 exhibited obvious semiannual variation observed at Fuke from March 2010 to June 2012 and increased with increasing solar activity. The difference of foF2 and foF3 in the months February, March, September, October and November is less evident in the middle solar activity period 2011–2012 than the low solar activity 2010 and in the other period it shows a slight increase (0.5 MHz) or keeps constant. The results show that the solar activity dependence of the F3 layer occurrence at low latitude away from the magnetic equator is different from that at near the magnetic equator.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

Dependence of the F-region peak electron density (foF2) on solar activity observed in the equatorial ionospheric anomaly region in the Brazilian sector

The long-term (solar cycle) changes in the Sun and how it affects the ionospheric F-region observed at São José dos Campos (23.2° S, 45.9° W), Brazil, a location under the southern crest of the equatorial ionospheric anomaly, have been investigated in this paper. The dependence of the F-region peak electron density (foF2) on solar activity during the descending phase of the 23rd solar cycle for the periods of high, medium, and low solar activity has been studied. The ionospheric F-region peak electron densities observed during high and medium solar activity show seasonal variations with maxima close to the equinox periods, whereas during the low solar activity the maxima during the equinox periods is absent. However, during the low solar activity only change observed is a large decrease from summer to winter months. We have further investigated changes in the different ionospheric F-region parameters (minimum virtual height of the F-region (h′F), virtual height at 0.834foF2 (hpF2), and foF2) during summer to winter months in low solar activity periods, 2006–2007 and 2007–2008. Large changes in the two ionospheric parameters (hpF2 and foF2) are observed during summer to winter months in the two low solar activity periods investigated.     

Statistical characterisation of spread F over South Africa

The occurrence of mid-latitude spread F (SF) over South Africa has not been extensively studied since the installation of the DPS-4 digisondes in 1996 and 2000 at Grahamstown (33.32 °S, 26.50 °E) and Madimbo (22.38 °S, 30.88 °E) respectively. This study is intended to quantify the probability of occurrence of F region disturbances associated with SF over South Africa. A study was conducted using data for 8 years (2001–2008) over Madimbo (with a time resolution of 30 min) and Grahamstown (with a variable time resolution of 15 and 30 min). In this study, SF has been classified into frequency SF (FSF), range SF (RSF) and mixed SF (MSF). The SF events were identified by manually identifying ionograms showing SF and tabulating them according to type for further statistical analysis. The results show that the diurnal pattern of SF peaks strongly between 01:00 and 02:00 local time, LT (LT = UT + 2 h), where UT is the universal time. This pattern is true for all seasons and types of SF at Madimbo and Grahamstown in 2001 and 2005, except for RSF which had peaks during autumn and spring in 2001 at Madimbo. The probability of both MSF and FSF tends to increase with decreasing solar activity, with a peak in 2005 (a moderate solar activity period). The seasonal peaks of MSF and FSF are more frequent during winter months at both Madimbo and Grahamstown. In this study, SF was evident in ∼0.03% and ∼0.06% of the available ionograms at Madimbo and Grahamstown respectively during the 8 years.     

Comparison of ionospheric F2 peak parameters foF2 and hmF2 with IRI2001 at Hainan

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with quarter-hourly time interval resolution for the diurnal variation, obtained with DPS4 digisonde at Hainan (19.5°N, 109.1°E; Geomagnetic coordinates: 178.95°E, 8.1°N) are used to investigate the low-latitude ionospheric variations and comparisons with the International Reference Ionosphere (IRI) model predictions. The data used for the present study covers the period from February 2002 to April 2007, which is characterized by a wide range of solar activity, ranging from high solar activity (2002) to low solar activity (2007). The results show that (1) Generally, IRI predictions follow well the diurnal and seasonal variation patterns of the experimental values of foF2, especially in the summer of 2002. However, there are systematic deviation between experimental values and IRI predictions with either CCIR or URSI coefficients. Generally IRI model greatly underestimate the values of foF2 from about noon to sunrise of next day, especially in the afternoon, and slightly overestimate them from sunrise to about noon. It seems that there are bigger deviations between IRI Model predictions and the experimental observations for the moderate solar activity. (2) Generally the IRI-predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the experimental results, but there is a relatively good agreement in summer at low solar activity. The deviation between the IRI-predicted hmF2 using CCIR M(3000)F2 and observed hmF2 is bigger from noon to sunset and around sunrise especially at high solar activity. The occurrence time of hmF2 peak (about 1200 LT) of the IRI model predictions is earlier than that of observations (around 1500 LT). The agreement between the IRI hmF2 obtained with the measured M(3000)F2 and the observed hmF2 is very good except that IRI overestimates slightly hmF2 in the daytime in summer at high solar activity and underestimates it in the nighttime with lower values near sunrise at low solar activity.     

Occurrence of equatorial plasma bubbles over Kolhapur

This paper reports the nightglow observations of OI 630.0 nm emissions, made by using all sky imager operating at low latitude station Kolhapur (16.8°N, 74.2°E and dip lat. 10.6°N) during high sunspot number years of 24th solar cycle. The images are analyzed to study the nocturnal, seasonal and solar activity dependence occurrence of plasma bubbles. We observed EPBs in images regularly during a limited period 19:30 to 02:30 LT and reach maximum probability of occurrence at 22:30 LT. The observation pattern of EPBs shows nearly no occurrence during the month of May and it maximizes during the period October–April. The equinox and solstice seasonal variations in the occurrence of plasma bubbles show nearly equal and large differences, respectively, between years of 2010–11 and 2011–12.     

Diurnal and seasonal variations of F2 layer characteristics over Irkutsk during the decrease in solar activity in 2003–2006: Observations and IRI-2001 model predictions

The diurnal and seasonal variations of F2 layer characteristics (critical frequency, peak height and bottomside thickness) over Irkutsk, Russia (52.3 N and 104.3 E) are studied by the method of running medians. The comparison with the IRI-2001 model during the decrease in solar activity in 2003–2006 revealed cases of both close agreement and systematic differences between predictions and observations. The systematic difference is not the only reason for disagreement between IRI and observations; there are also intrayear variations which are not associated with seasonal behavior. The period of observation was too short to make conclusions about solar activity dependence of the noon bottomside thickness and the modification of its diurnal behavior with decreasing solar activity.     

Variations in statistical parameters of the NmF2 equinoctial asymmetry with latitude and solar activity near noon

Hourly values of the F2-layer peak density, NmF2, measured by 95 ionosondes near noon from 1957 to 2011 at low and middle geomagnetic latitudes of the northern and southern geographic hemispheres are used in a statistical study of the NmF2 equinoctial asymmetry. The ratios, R, of NmF2 measured during 61 days around the March equinox to NmF2 measured during 61 days around the September equinox at the same UT near noon during geomagnetically quiet daytime conditions for approximately the same solar activity conditions over the same ionosonde are analyzed. The conditional probability of the occurrence of R in an interval of R, the most probable value of R, and the mean expected value of R are calculated for the first time for the low, moderate, and high solar activity levels to study variations in these statistical parameters with latitude and solar activity. These statistical parameters are averaged over 5° geomagnetic latitude interval in the northern and southern geographic hemispheres to calculate and to study for the first time trends in latitude and solar activity of these averaged NmF2 equinoctial asymmetry statistical characteristics.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Effects observed in the ionospheric F-region in the South American sector during the intense geomagnetic storm of 14 December 2006

This investigation presents observations related to the generation of equatorial ionospheric irregularities (also known as equatorial spread F (ESF)) including ionospheric plasma bubbles and dynamic behavior of the ionospheric F-region in the South American sector during an intense geomagnetic storm in December 2006 (a period of low solar activity). In this work, ionospheric sounding observations and GPS data obtained between 13 and 16 December 2006 at several stations in the South American sector are presented. On the geomagnetically disturbed night of 14 and 15 December, ionospheric plasma bubbles were observed after an unusual uplifting of the F-region during pre-reversal enhancement (PRE) period. The unusual uplifting of the F-region during PRE was possibly associated with prompt penetration of electric field of magnetospheric origin. During the geomagnetic disturbance night of 14 and 15 December, strong oscillations due to the propagation of traveling ionospheric disturbances (TIDs) by the Joule heating in the auroral region were observed in the F-region at São José dos Campos (SJC, 23.2°S, 45.9°W; dip latitude 17.6°S), Brazil, and Port Stanley (PST, 51.6°S, 57.9°W; geom. latitude 41.6°S). The VTEC-GPS observations presented on the night of 14 and 15 December 2006 show both positive and negative storm phases in the South American sector, possibly due to changes in the large-scale wind circulation and changes in the O/N₂ ratio in the southern hemisphere, respectively.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

Features of the occurrence of the additional stratification on the bottom-side F region over the equatorial location of Trivandrum

The general features of occurrence of an additional layer on the bottom side of F region, referred to as F0.5 layer in the pre noon period, over the magnetic equatorial location of Trivandrum (8.5° N; 77° E; dip lat of 0.5° N) in India during the period from 2004 to 2007 are presented using ionosonde observations. The F0.5 layer has a June (northern summer) solsticial maximum probability of occurrence with secondary maxima during December (northern winter) solstice. The seasonal as well as the day-to-day variability in the occurrence of F0.5 layer as mentioned in this paper seems to be a result of the variations in the amplitude and phases of the tides and gravity waves, and inventory of the metallic ions of meteoric origin. This study brings out an important manifestation of morning time F layer base region dynamics.     

Radio occultation electron density profiles from the FORMOSAT-3/COSMIC satellites over the Brazilian region: A comparison with Digisonde data

This study aims to validate the electron density profiles from the FORMOSAT-3/COSMIC satellites with data from Digisondes in Brazil during the low solar activity period of the years 2006, 2007 and 2008. Data from three Brazilian Digisondes located in Cachoeira Paulista (22.7°S, 45°W), São Luís (2.5°S, 44.2°W) and Fortaleza (3.8°S, 38°W) were used in the comparisons. Only the profiles whose density peak have been obtained near the stations coordinates were chosen for the comparison. Although there is generally good agreement, some cases of discrepancies are observed. Some of these discrepancies cannot be explained simply by the differences in the position and local time of the measurements made by the satellite and the ground-based station. In such cases it is possible that local conditions, such as the presence of a trans-equatorial wind or electron density gradients, could contribute to the observed differences. Comparison of the F2 layer peak parameters, the NmF2 and hmF2, obtained from the two techniques showed that, in general, the agreement for NmF2 is pretty good and the NmF2 has a better correlation than hmF2. Cachoeira Paulista had the worst correlation for hmF2 possibly because this station is situated in the region under the influence of the equatorial ionization anomaly, a region where it is more difficult to apply the RO technique without violating the spherical symmetry condition.     

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Vertical TEC variations and model during low solar activity at a low latitude station,Xiamen

GPS-derived vertical TEC recorded at Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), China, during year 2006 is analyzed for the first time and compared to that predicted by ionosphere model SPIM recommend by ISO. A manifest seasonal anomaly is found with the high value during equinoctial season and low value during summer and winter season. Relative standard deviation for VTEC shows high value at around midnight and before sunrise. The correlation analysis exhibits that the variation of VTEC has a very weak relation with geomagnetic and solar activities (Dst, AP, SSN and F10.7). Comparative results reveal that the SPIM overestimates the observed VTEC at most of the time.     

Comparison of diurnal,seasonal and latitudinal effect of MUF VR and NmF2 VR during some solar cycle epochs

The diurnal and seasonal changes of the variability (VR) of Maximum Useable Frequency (MUF) are compared with those of peak electron density (NmF2) at Ibadan (7.4°N, 3.9°E, 6°S dip) in the African sector. Also compared is the latitudinal effect on both characteristics by combining data from Singapore (1.3°N, 103.8°E, 17.6°S dip) in the East Asian sector and Slough (51.5°N, 359.4°E, 66.5°S dip) in the European sector. MUF VR is found to be about half of NmF2 VR at all the hours and seasons and during the solar cycle epochs considered for the three stations. While nighttime MUF VR is greater in June Solstice and September Equinox during both low and moderate solar activities and in September Equinox and December Solstice during high solar activity, nighttime NmF2 VR is greater in June Solstice and September Equinox during high solar activity and greater at the equinoxes during low and moderate solar activities. This signifies a shift in nighttime MUF peak VR from the middle six months during low and moderate solar activities to the last half of the year during high solar activity. Daytime VR of both characteristics are not observed to show any seasonal variation. MUF VR and that of NmF2 are found to increase and decrease alternately with the Zurich sunspot number (Rz) for Ibadan and Singapore. For Slough, the VR of both characteristics increases with Rz during the first half of the day. It then increases and decreases alternately with Rz during the remaining hours of the day. While nighttime MUF VR decreases with latitude, just like nighttime NmF2 VR, no latitudinal effect is found for daytime VR of both characteristics.     

Monthly mean foF2 model for an African low-latitude station and comparison with IRI

We have employed the hourly values of the ionospheric F-region critical frequency (foF2) obtained from Ouagadougou ionosonde, Burkina Faso (geographic coordinates 12° N, 1.8° W) during the interval of 1985–1995 (solar cycle 22) and solar radio flux of 10 cm wavelength (F10.7) to develop a local model (LM) for the African low-latitude station. The model was developed from regression analysis method, using the two-segmented regression analysis. We validated LM with foF2 data from Korhogo observatory, Cote d’Ivorie (geographical coordinates 9.3° N, 5.4° W). LM as well as the International Reference Ionosphere (IRI) agrees well with observations. LM gave some improvement on the IRI-predicted foF2 values at the sunrise (06 LT) at all solar flux levels and in all seasons except June solstice. The performance of the models at the representing the salient features of the equatorial foF2 was presented. Considering daytime and nighttime performances, LM and IRI are comparable in low solar activity (LSA), LM performed better than IRI in moderate solar activity (MSA), while IRI performed better than LM in high solar activity (HSA). CCIR has a root mean square error (r.m.s.e), which is only 0.10 MHz lower than that of LM while LM has r.m.s.e, which is about 0.05 MHz lower than that of URSI. In general, our result shows that performance of IRI, especially the CCIR option of the IRI, is quite comparable with the LM. The improved performance of IRI is a reflection of the numerous contributions of ionospheric physicists in the African region, larger volume of data for the IRI and the diversity of data sources, as well as the successes of the IRI task force activities.     

Predicting the probability of occurrence of spread-F over Brazil using neural networks

The probability of occurrence of spread-F can be modeled and predicted using neural networks (NNs). This paper presents a feasibility study into the development of a NN based model for the prediction of the probability of occurrence of spread-F over selected equatorial stations within the Brazilian sector. The input space included the day number (seasonal variation), hour (diurnal variation), sunspot number (measure of the solar activity), magnetic index (measure of the magnetic activity) and magnetic position. Twelve years of spread-F data from Brazil (covering the period 1978–1989) measured at the equatorial site Fortaleza (3.9°S, 38.45°W) and low latitude site Cachoeira Paulista (22.6°S, 45.0°W) are used in the development of an input space and NN architecture for the model. Spread-F data that is believed to be related to plasma bubble developments (range spread-F) was used in the development of the model. The model results show the probability of spread-F occurrence as a function of local time, season and latitude. Results from the Brazilian Sector NN (BSNN) based model are presented in this paper, as well as a comparative analysis with a Brazilian model developed for the same purpose.