Variations of scale height at F-region peak based on ionosonde measurements during solar maximum over the crest of equatorial ionization anomaly region.
Investigation of the ionospheric electron density profile Ne(h) is an important topic for many practical problems and contributes to our understanding of the structure of the ionosphere, ionospheric dynamics, and radio wave for telecommunications. The most important ionospheric region to investigate is the ionospheric F region, which is divided into three subregions: F1 and F2, which are considered the bottomside of the ionosphere, and the upper F region also known as the topside ionosphere. F1 is the lowest region and is dominated by photochemical processes, whereas the upper F region is dominated by diffusion. F2 is in between these two regions, where there is a transition from chemical to diffusion dominance. Examining variations in the ionosphere's shape will help us to study the physical processes in the ionosphere.
The ionospheric bottomside electron density profile was derived from ground-based observations using an ionosonde, digisonde, and incoherent scatter radar in the form of ionograms. The ionospheric topside electron density profile can be obtained from topside sounder measurements or by theorematic models. Recently, many analytical functions and empirical modeling techniques have been reported to study the topside ionosphere, such as Chapman's function and exponential, parabolic, and Epstein functions [1-8]. Reinisch and Huang  used scale heights derived from ionograms recorded at Millstone Hill (a middle latitude site) and Jicamarca (located at the geomagnetic equator) stations for modeling the topside ionospheric profile. They proposed a convenient method to derive the topside ionospheric profile based on the Chapman scale height, [H.sub.m], around the F2-layer peak height. They have shown that variations in scale height derived using the [alpha]-Chapman function above the F2-layer peak are very small and have assumed a constant scale height for the interpolation of the topside ionosphere.
Over the past few years, many studies have been conducted on this subject, like that of [9-15]. Tulasi Ram et al.  studied [H.sub.m] around the F2-layer at all 13 stations and have shown maximum values in summer during daytime. Furthermore, Zhang et al.  also studied the [H.sub.m] at the Fregion peak height over Hainan (19.4[degrees]N, 109[degrees]E) and showed two conspicuous peaks occurring at local noon and presunrise. Their results show a high correlation between [H.sub.m] and B0 at different latitudes and equivalent slab thicknesses, r, at low latitudes. In addition, Lee and Reinisch  presented a postsunset peak in [H.sub.m] during the equinox and summer at the equatorial region, Jicamarca (12.0[degrees]S, 76.9[degrees]W), during high solar activity. Although numerous studies of scale height have been carried out in many locations, it has not yet been examined at the northern crest of the equatorial ionization anomaly (EIA) area.
Since [H.sub.m] plays an important role in the ionospheric vertical profile, especially for deriving and studying the topside ionosphere, we have analyzed the diurnal and seasonal variations of [H.sub.m] in order to study the ionospheric dynamics in the crest of the EIA area. In this work, we examine the ionograms recorded in 1999 at Chung-Li (24.9[degrees]N, 121.1[degrees]E) to derive [H.sub.m] during high solar activity. Furthermore, this investigation also studies the correlation between [H.sub.m], F2-layer peak height (hmF2), B0, and [tau] to investigate ionospheric physics processes over the crest of the EIA in East Asia.
2. Material and Methods
In our study, we collected data on the F2-layer critical frequency (foF2), hmF2, B0, [H.sub.m], and the total electron content (TEC) from the Chung-Li ionosonde station (24.9[degrees]N, 121.1[degrees]E) and YMSM GPS receiver (25.2[degrees]N, 121.6[degrees]E). We processed the ionogram traces manually and transferred them to computers, capturing the ionospheric true height profiles for more than 34,000 ionograms taken in 1999 during high solar activity. These data are classified into three periods: equinox (March, April, September, and October), summer (May-August), and winter (November-February).
The true height plasma frequency profile is calculated using the standard true-height inversion POLynomial ANalysis program (POLAN) . POLAN is widely accepted for inversion of ionograms obtained by classical ionosonde  and allowed our measurements of foF2, hmF2, B0, and [H.sub.m]. This study processes a given set of virtual-height data from ionograms with a manual work. Results obtained by POLAN are therefore normally stored by arrays giving the scaled frequencies and corresponding real heights. The ionograms were recorded every 15 minutes at Chung-Li station. GPS TEC data are used to calculate r, which is defined as the ratio of TEC/NmF2.
The goal of this paper is focused on the day-to-day and seasonal variations in [H.sub.m] and the correlation between [H.sub.m] and other ionospheric parameters such as hmF2, B0, and [tau]. In addition, we tried to get a clearer picture of ionospheric dynamics in the northern crest of EIA area.
3. Results and Discussion
In order to obtain average diurnal variations in foF2, hmF2, B0, [H.sub.m], and [tau] values every month during the whole year, we calculate median values for all the days in each month in 1999 during periods of high solar activity and smooth out the data. These results are illustrated in Figure 1 and represent seasonal variations in [H.sub.m], with a maximum appearing during daytime in June and August; lower values appear during nighttime in the winter months (Figure 1(a)). In addition, the daytime peak of [H.sub.m] also shows a seasonal variation; during the equinox and summer, the peak appears during 0800-1030 LT and, in winter, the peak appears later around 1030-1500 LT. Figure 1(b) shows the contour plot in [H.sub.m]F2, which has an apparent higher value during 1000-1300 LT in summer (May-August). Moreover, another conspicuous peak was observed around 2200 LT in summer. The summer daytime peak in [H.sub.m] occurs at an earlier time than hmF2 since the increase in [H.sub.m] could be associated with the sunrise. That is to say, electron density begins to increase rapidly owing to photoionization at sunrise and then rises slowly during the day . It is known that the [[H.sup.+]]/[[O.sup.+]] ratio is very important at the F2 peak. After sunrise, the meridional neutral wind blows toward the pole and produces a diffusion effect from the topside ionosphere to lower altitudes. Furthermore, the bottomside ionospheric electron density still increases from photoionization during daytime, which is associated with the increase in [[O.sup.+]] and uplift of the F2 layer to a higher altitudes, and it reduces the loss rate and thereby still increases the electron content around the F2 peak.
Figures 2(a)-2(c) show the monthly median values of [H.sub.m], hmF2, and foF2 in the F2 region for equinoctial months: March (solid line), April (gray rectangle line), September (dotted line), and October (plus line) in 1999. Figure 2(a) represents the diurnal variation in [H.sub.m]; daytime values were greater than nighttime values. In addition, two apparent peak values occur at about 0500 and 0900-1100 LT. The first [H.sub.m] peak appears around 0400-0500 LT and increases to about 56 km and then falls rapidly to a trough value of approximately 38 km. The second [H.sub.m] peak appears during 0900-1100 LT at 60-65 km and then decreases gradually.
Figure 2(b) shows the variation in hmF2 during equinoctial months. In general, there are three peaks to be observed: presunrise, local noon, and postsunset during the equinox, but the postdusk peak is not obvious in March and October. Figure 2(c) shows a diurnal variation in foF2 and a peak and trough occurring during 1300-1600 LT and around 0500 LT, respectively, for all seasons. Comprehensively, the peak [H.sub.m] during the equinox, which appears at 0400-0500 LT, contrasts with the increase in hmF2 during presunrise period while foF2 decreases to its minimum during the same period. This result indicates that the increase in [H.sub.m] is associated with the equatorward neutral wind and leads to an increase in hmF2 and in the [N.sub.2]/O ratio and also causes a decrease in foF2 . In addition, the shape of ionosphere is very sensitive to variations in the ratio of [[H.sup.+]]/[[O.sup.+]] at the F2 peak . Moreover, the [O.sup.+] density decreases with altitude . Hence, the enhanced hmF2 decreases [O.sup.+] which leads to an increase in the [[H.sup.+]]/[O.sup.+] ratio and increases [H.sub.m] at the F2 peak. The main peak in [H.sub.m] starts increasing at 0600 LT and reaches its maximum value (~65 km) at about 1000 LT. During the same time period, hmF2 also increases from the trough (~250 km) to its maximum (~380 km) and then decreases gradually.
Fejer et al.  studied equatorial ionospheric vertical plasma drifts and showed a vertical upward plasma drift from 0600 to 1000 LT, and then the drift turned downward. The results seem to correspond to the variations in [H.sub.m] at the EIA, which indicates the effect of dynamo electric fields during daytime. The slight increase in [H.sub.m] during 1400-1600 LT is associated with the increase in plasma and decrease in hmF2 (Figures 2(b) and 2(c)).
Figures 3 and 4 show the same data as Figure 2, but for the summer and winter months. During summer, [H.sub.m] shows a diurnal variation with the main peak appearing during 09001000 LT (Figure 3(a)). The main peak in [H.sub.m] is higher in June (-100 km) and lower in August (-70 km). Figure 3(b) shows the variation in hmF2 in summer and three peak values to be observed at presunrise, local noontime, and postsunset periods. Obviously, the peak value in [H.sub.m] does not correspond to hmF2 during the daytime period. The variation in [H.sub.m] is also attributed to the vertical plasma drifts induced by dynamoelectric fields . In addition, a slight presunrise peak in [H.sub.m] is observed in summer. This increase in [H.sub.m] contrasted with a slight uplift in hmF2 during the presunrise period. A slight increase in [H.sub.m] during 1500-1600 LT is concerned with the increase in foF2 (Figure 3(c)). Furthermore, a slight peak in [H.sub.m] is observed during 1800-2100 LT and is also associated with an increase of hmF2, which is attributed to the equatorial prereversal enhancement caused by the occurrence of a strong postsunset eastward electric field .
In winter months, [H.sub.m] shows two apparent peak values during presunrise and local noontime periods, with a maximum value of about 60 km at noon (Figure 4(a)). The main peak in [H.sub.m] starts to increase from 40 km at 0600 LT and reaches its maximum value (-60 km) at 1200 LT, and foF2 also increases simultaneously and lasts until 1600 LT (Figure 4(c)). The variation in [H.sub.m] also shows a similar trend that compares with the results of Fejer et al. , who studied the global vertical plasma drift during daytime at the equator. The presunrise peak in [H.sub.m] is attributed to an uplift hmF2 (Figure 4(b)). If [H.sub.m] is correlated to classic scale height, H, the temperature should be tightly correlated with [H.sub.m].
Figure 5 illustrates the variations in temperature of neutral gas, ions, and electrons ([T.sub.n], [T.sub.i], and [T.sub.e], resp.) and is measured using an IRI-2007 model over the Chung-Li location during 1999. Figure 5(a) illustrates the variations in [T.sub.n] during the equinox, summer, and winter. It shows a trough at 0500 LT and then increases to a maximum at around 1600 LT for all three seasons. It corresponds highly with the variations in foF2 and minimally with [H.sub.m] during the daytime period. For the [T.sub.i] (Figure 5(b)), the results represent a tendency similar to [T.sub.n] for all three seasons. Figure 5(c) shows [T.sub.e] variations during different seasons. An abrupt increase can be seen from about 800[degrees]K to 2600[degrees]K during 0500-0700 LT and then a rapid decrease to 1200[degrees]K during the equinox. These results indicate that the temperature of neutral gas had stronger correlation with foF2 and weaker correlation with [H.sub.m].
In addition, we also examined the correlation between [H.sub.m] and ionospheric profile parameters such as ionospheric bottomside shape, B0, in the IRI model and the ionospheric equivalent slab thickness, [tau]. Figure 6 shows the linear regression analysis between [H.sub.m] and hmF2, B0, and [tau]. The top panels show the linear correlation between [H.sub.m] and hmF2. The result shows a weak correlation (r varies in the range of -0.096 to 0.136) between [H.sub.m] and hmF2, except for a middle correlation (r = 0.631) in winter. Furthermore, the resulting [H.sub.m] compares with the result of vertical plasma drift during daytime at the equator reported by Fejer et al.  and finds that trends were similar during daytime. Our results indicate that [H.sub.m] in the EIA area was strongly dominated by the E x B drift at the equator. The medium panel shows the correlation analysis between [H.sub.m] and B0 during three seasons and represents a high correlation coefficient (r varies in the range 0.962-0.991) during the high solar activity period. This result is in agreement with Liu et al.  who indicated that [H.sub.m] and B0 have a strong correlation at all times over Wuhan and over 13 other stations.
Further, Lee and Reinisch  studied [H.sub.m] variations at an equatorial station, Jicamarca, and also indicated that [H.sub.m] and B0 are highly correlated. Therefore, these results indicate that [H.sub.m] values can be estimated using B0 values, which also corresponds to previous results at mid- and equatorial latitudes [9,10,12,13]. Meanwhile, a study  indicated a strong correlation between [H.sub.m] and [tau] at Millstone Hill. Zhang et al.  studied the variation of [H.sub.m] at Hainan (19.4[degrees]N, 109.0[degrees]E) and also proposed a high correlation between [H.sub.m] and r during 2002-2004. Therefore, we also examined the comparison between [H.sub.m] and r. The seasonal values of [H.sub.m] show a similar trend to r over Chung-Li during high solar activity periods. Figure 6 (bottom panel) shows the correlation analysis between [H.sub.m] and r during all seasons and represents a moderate correlation (r = 0.584) in summer and a weak or negative correlation during equinox (r = 0.148) and winter (r = -0.210). This result differs from the studies at other latitudes, mentioned above, since the EIA is more complex in its contribution to the electron density in the topside ionosphere, which is affected by electric fields and thermospheric winds. Therefore, the correlation was weak between [H.sub.m] and [tau].
This paper presents the diurnal and seasonal variations in Chapman scale height [H.sub.m] around the F2-layer peak height at the northern crest of the EIA region, over the Chung-Li area station in Taiwan, during periods of high solar activity in 1999. [H.sub.m], hmF2, foF2, and B0 values were derived from electron density profiles, which were calculated using the POLAN program. The key results obtained are as follows.
(1) During the periods of high solar activity, the diurnal variations in [H.sub.m] are generally characterized by higher daytime values as compared to the nighttime values in all the three seasons. Furthermore, monthly median daytime values of [H.sub.m] are highest in summer and lowest in winter.
(2) A presunrise peak in [H.sub.m] is observed in all three seasons during high solar activity, particularly in the winter. This phenomenon is also observed in hmF2, whose increase is also more prominent during winter. These results are also observed at midlatitudes [10,13] and equatorial latitudes . Meanwhile, the presunrise peak in [H.sub.m] is attributed to the uplift hmF2, which was caused by the zonal neutral wind. Murthy et al.  have mentioned that the meridian neutral wind turns southward at 0300 LT and reverse at 0500 LT for low latitudes. Hence, the F2 layer rises to higher altitude and causes a decrease in foF2 such that there is an increase in the ratio of [N.sub.2]/O.
(3) During the day, [H.sub.m] is highest in summer and strongly dominated by the vertical plasma drift produced at the equator. Further, our results show a poor correlation between [H.sub.m] and hmF2, except for a middle correlation in winter. This vertical plasma drift due to a strong eastward electric field at the equator leads to plasma uplift and then diffuses down the magnetic field lines to the EIA region and causes an increase in the thickness (r) of the topside ionosphere and [H.sub.m].
(4) During the postsunset period, a slight increase in [H.sub.m] is observed in summer for high solar activity periods and is associated with an increase in hmF2 caused by prereversal enhancement at the equator.
(5) A strong correlation is found between [H.sub.m] and B0 values for all three seasons (r varies in the range 0.9620.991). This result is also observed at equatorial latitudes and midlatitudes and suggests that [H.sub.m] could be estimated by B0 [9, 12, 13]. In addition, [H.sub.m] shows a poor correlation with r, which differs from the results of Zhang et al. , who showed a strong correlation in a low-latitude station in Hai-nan.
(6) This study provides more information about the ionospheric shape parameters during the solar maximum period over the crest of EIA region. The [H.sub.m] parameter could be provided to compare with the ionospheric topside model and the radio occultation measurements on satellites. Besides, the bottomside ionospheric slab thickness, B0, could be estimated by the [H.sub.m] to be easy in measuring B0 values. Furthermore, the [H.sub.m] provides the information of ionospheric photoionization reactions and the ionospheric dynamics over the EIA region.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
The author greatly thanks the editor and the reviewers for their valuable comments and suggestions. The author would also like to thank the Chung-Li Ionosphere Station National Communications Commission and Central Weather Bureau of Taiwan for providing the ionograms and GPS data online. This work was supported by the Grants of the Ministry of Science and Technology NSC102-2111-M-275-001, NSC1022119-M-275-001, and MUST103-2m-M-275-001.
 H. G. Booker, "Fitting of multi-region ionospheric profiles of electron density by a single analytic function of height," Journal of Atmospheric and Terrestrial Physics, vol. 39, no. 5, pp. 619-623, 1977.
 K. Rawer, D. Bilitza, and T L. Gulyaeva, "New formulas for the IRI electron density profile in the topside and middle ionosphere," Advances in Space Research, vol. 5, no. 7, pp. 3-12, 1985.
 K. Rawer, "Synthesis of ionospheric electron density profiles with Epstein functions," Advances in Space Research, vol. 8, no. 4, pp. 191-199, 1988.
 G. di Giovanni and S. M. Radicella, "An analytical model of the electron density profile in the ionosphere," Advances in Space Research, vol. 10, no. 11, pp. 27-30, 1990.
 B. W. Reinisch and X. Huang, "Deducing topside profiles and total electron content from bottomside ionograms," Advances in Space Research, vol. 27, no. 1, pp. 23-30, 2001.
 S. M. Stankov, N. Jakowski, S. Heise, P. Muhtarov, I. Kutiev, and R. Warnant, "A new method for reconstruction of the vertical electron density distribution in the upper ionosphere and plasmasphere," Journal of Geophysical Research: Space Physics, vol. 108, no. A5, 2003.
 B. W. Reinisch, X. Q. Huang, A. Belehaki, J. H. Shi, M. L. Zhang, and R. Ilma, "Modeling the IRI topside profile using scale heights from ground-based ionosonde measurements," Advances in Space Research, vol. 34, no. 9, pp. 2026-2031, 2004.
 B. W. Reinisch, P Nsumei, X. Huang, and D. K. Bilitza, "Modeling the F2 topside and plasmasphere for IRI using IMAGE/RPI and ISIS data," Advances in Space Research, vol. 39, no. 5, pp. 731-738, 2007.
 L. Liu, W. Wan, and B. Ning, "A study of the ionogram derived effective scale height around the ionospheric hmF2," Annales Geophysicae, vol. 24, no. 3, pp. 851-860, 2006.
 M. Zhang, B. W. Reinisch, J. Shi, S. Wu, and X. Wang, "Diurnal and seasonal variation of the ionogram-derived scale height at the F2 peak," Advances in Space Research, vol. 37, no. 5, pp. 967971, 2006.
 M. Mosert, R. Ezquer, B. de la Morena, D. Altadill, G. Mansilla, and G. Miro Amarante, "Behavior of the scale height at the F2 layer peak derived from Digisonde measurements at two European stations," Advances in Space Research, vol. 39, no. 5, pp. 755-758, 2007.
 C.-C. Lee and B. W. Reinisch, "Quiet-condition variations in the scale height at F2-layer peak at Jicamarca during solar minimum and maximum," Annales Geophysicae, vol. 25, no. 12, pp. 2541-2550, 2007.
 F. Nambala, L. McKinnell, and E. Oyeyemi, "Variations in the ionospheric scale height parameter at the F2 peak over Grahamstown, South Africa," Advances in Space Research, vol. 42, no. 4, pp. 707-711, 2008.
 S. Tulasi Ram, S.-Y. Su, C. H. Liu, B. W. Reinisch, and L. Lee McKinnell, "Topside ionospheric effective scale heights (HT) derived with RUCSAT-1 and ground-based ionosonde observations at equatorial and midlatitude stations," Journal of Geophysical Research: Space Physics, vol. 114, no. 10, 2009.
 J. E. Titheridge, "Ionogram analysis with the generalized program POLAN," Report UAG-93, World Data Center A for SolarTerrestrial Physics, 1985.
 L. F. McNamara, D. L. Cooke, C. E. Valladares, and B. W. Reinisch, "Comparison of CHAMP and digisonde plasma frequencies at jicamarca, peru," Radio Science, vol. 42, no. 2, Article ID RS2005, 2007.
 R. W. Schunk and A. F. Nagy, in Ionospheres: Physics, Plasma Physics, and Chemistry, J. A. Dessler, J. T. Houghton, and M. J. Rycroft, Eds., Cambridge Atmospheric and Space Science Series, pp. 335-375, Cambridge University Press, Cambridge, UK, 2000.
 G. W. Prolss, "Storm-induced changes in the thermospheric composition at middle latitudes," Planetary and Space Science, vol. 35, no. 6, pp. 807-811, 1987
 K. Davies and X. M. Liu, "Ionospheric slab thickness in middle and low latitudes," Radio Science, vol. 26, no. 4, pp. 997-1005, 1991.
 B. G. Fejer, E. R. de Paula, R. A. Heelis, and W. B. Hanson, "Global equatorial ionospheric vertical plasma drifts measured by the AE-E satellite," Journal of Geophysical Research, vol. 100, pp. 5769-5776, 1995.
 D. T Farley, E. Bonelli, and B. G. Fejer, "The pre-reversal enhancement of the zonal electric field in the equatorial ionosphere," Journal of Geophysical Research, vol. 91, pp. 1372313728, 1986.
 X. Huang and B. W. Reinisch, "Vertical electron content from ionograms in real time," Radio Science, vol. 36, no. 2, pp. 335342, 2001.
 B. K. K. Murthy, S. S. Hari, and V V. Somayajulu, "Nighttime equatorial thermospheric meridional winds from ionospheric h'F data," Journal of Geophysical Research: Space Physics, vol. 95, no. A4, pp. 4307-4310, 1990.
Department of Information Technology, Ling Tung University, Taichung 408, Taiwan
Correspondence should be addressed to Yu-Jung Chuo; email@example.com
Received 1 June 2014; Accepted 3 July 2014; Published 4 August 2014
Academic Editor: Jianshe Lei
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2014|
|Previous Article:||Coordinated control of slip ratio for wheeled mobile robots climbing loose sloped terrain.|
|Next Article:||Essential oil content of the rhizome of Curcuma purpurascens Bl.|