DC Langmuir Probe for Measurement of Space Plasma: A Brief Review
DC Langmuir Probe for Measurement of Space Plasma: A Brief Review
Journal of Astronomy and Space Sciences. 2015. Sep, 32(3): 167-180
Copyright © 2015, The Korean Space Science Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : April 21, 2015
  • Accepted : June 27, 2015
  • Published : September 15, 2015
Export by style
Cited by
About the Authors
Koichiro, Oyama
Institute of Plasma and Space Science, National Cheng Kung University, Tainan 70101, Taiwan

Herein, we discuss the in situ measurement of the electron temperature in the ionosphere/plasmasphere by means of DC Langmuir probes. Major instruments which have been reported are a conventional DC Langmuir probe, whose probe voltage is swept; a pulsed probe, which uses pulsed bias voltage; a rectification probe, which uses sinusoidal signal; and a resonance cone probe, which uses radio wave propagation. The content reviews past observations made with the instruments above. We also discuss technical factors that should be taken into account for reliable measurement, such as problems related to the contamination of electrodes and the satellite surface. Finally, we discuss research topics to be studied in the near future.
A DC Langmuir probe measures two basic parameters: electron density ( Ne ) and electron temperature ( Te ) in plasma ( Langmuir & Mott-Smith 1924 ). Since a DC Langmuir probe was first flown in a V-2 rocket from White Sands in New Mexico, USA ( Reifman & Dow 1949 ) to measure the lower E-region ionosphere, the probe has been flown in many sounding rockets and satellites there ( Boggess et al. 1959 ; Spencer et al. 1965 ; Brace et al. 1973 ; Brace 1998 ) and in Russia ( Andreyeva et al. 1971 ). Even now, DC Langmuir probes are generally employed in low earth orbit satellites, such as Astrid ( Holback et al. 2001 ), Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions (DEMETER) ( Lebreton et al. 2006 ), Challenging Minisatellite Payload (CHAMP) ( Cooke et al. 2003 ), SWARM ( Buchert et al. 2015 ), and Communications/Navigation Outage Forecasting System (C/NOFS) ( Costa et al. 2014 ). This is because Ne and Te are two basic parameters that are fundamental for study of the ionosphere/plasmasphere. It goes without saying that without Langmuir probe measurements, the detailed understanding of the ionosphere that we now have would not have been possible.
Fig. 1 shows the principle of the measurement. Varying voltage is applied to an electrode immersed in ionospheric plasma, with respect to the satellite (rocket) body, the electron and ion currents, I, that flow into and/or out from the electrode are measured. The I-V curve thus measured is fitted to a simple equation I = Ieexp (- eV/kTe )– Ip for the electron retarding region. V is measured from the potential of ambient plasma. In the equation, the electron saturation current, Ie = Se · Ne · e ·( kTe/2πme ) 1/2 , and Ip = bV + c . Se , Ne , e , k , and me are the surface area of the electrode, electron density, electron charge, Boltzmann’s constant, and electron mass, respectively. b is the slope of ion current, and c is the ion saturation current. Fig. 2 shows two I-V curves. One is the curve that was measured after the electrode was bombarded by highly negative voltage, and the other was measured by an electrode that was only washed by a chemical liquid, such as acetone. Here, Te and Ne cannot be obtained if the electrode is not ion bombarded or is not heated under vacuum. Thus, the principle of the DC Langmuir probe is simple and measurement appears to be easy. However, despite this simplicity, it does not always provide reliable results, especially for Te measurement. Even now we encounter problems in using Te data that are erroneously measured by satellites ( Kakinami et al. 2013 ). Ne is automatically calculated incorrectly because Ne is a function of Te as is related to by the equation Ie = Se · Ne · e ( kTe /2πme ) 1/2 in Langmuir probe theory, in addition to the reduction of probe current caused by electrode contamination.
PPT Slide
Lager Image
Electric circuit for probe measurement in space: An electrode of any shape (sphere, cylinder, or plane) is injected into the ionosphere/plasmasphere. The current (electron and ion current) flows into the electrode, satellite frame (rocket body), and ambient plasma, and it returns again to the electrode. A DC Langmuir probe assumes an infinite conductive surface area of the counter electrode (satellite surface, or conductive part of the rocket). For tiny satellites, the surface area becomes insufficient for DC probe measurement.
PPT Slide
Lager Image
I-V curve (Langmuir probe characteristic) when the electrode is clean and unclean: The horizontal axis is the probe voltage, and the vertical axis is the current. Note that the I-V curve provided by a contaminated electrode is not the same when the probe voltage is increased from low voltage to higher voltage. An I-V curve provided by a clean (ion bombarded in plasma) electrode shows no serious effects due to surface contamination.
We must provide reliable and accurate data to users so that the effort, time, and money spent for rocket and satellite missions are not spent in vain or inefficiently. In this review, we mainly discuss Te measurement because calculation of Te from DC Langmuir probe data very often provides inaccurate values, even if the instrument and the sensor are carefully maintained. It is noted however, that the electron collection mechanism is not simple when an electrode is moving because the wake behind the electrode makes the theory complicated ( Séran et al. 2005 ). The purpose of this paper is to make scientists aware of the problems that occur when Langmuir probes are used and to provide the information needed to obtain reliable and accurate electron temperature measurements. We will start with measurement in the E region; then proceed to measurement in the F /topside region and the plasmasphere and discuss two major problems associated with measurements. We finally end with a discussion of the subjects that need to be studied in the future.
- 2.1 E region
In the lower E region, we still cannot fully confirm the value of Te in the lower E region ( Oyama et al. 1980 ; Oyama 2000 ) even at mid-latitude, and in spite of the long history of sounding-rocket experiments. Most of the measurements report daytime Te , which is higher than the model neutral temperature Tn ( Smith et al. 1968 ; Schutz & Smith 1976 ; Rohde et al. 1993 ). During the daytime, the ionospheric electron density is about 10 5 els/cm 3 at the altitude of 95-120 km, while the density of neutral particles is 10 12 particles/cm 3 . That is, neutral density is 10 7 times higher than Ne . Accordingly, Te is expected to be nearly equal to Tn . We know that most of the measurements provided in 1970 are not correct ( Oyama et al. 1980 ) mainly due to electrode contamination ( Oyama 1976 ; Oyama et al. 2012 ). However, even when an electrode with a clean surface is used ( Hirao & Oyama 1972 ; Oyama & Hirao 1976 ; Oyama et al. 2008a , b ) or another instrument is used whose measurement is not influenced by electrode contamination, the measurements still show Te to be higher than Tn during the daytime. The Te value measured by a resonance cone instrument agrees with this measurement ( Piel et al. 1992 ), as well as that obtained by incoherent scatter radar ( Smith et al. 1978 ).
Electron temperature inside and/or around the sporadic E layer ( Es ) is controversial. The Te value measured by Aubry was about 600K in the Es layer ( Aubry et al. 1966 ), which is higher than that of outside plasma. The value measured by Schutz & Smith (1976) was about 380K inside Es , while outside Es , the Te value at the lower edge and upper edge of the Es layer were found to be 500K and 550K, respectively. The measurement was done by using a pulsed probe. According to Yoshimura et al. (1999) , Te was about 1,000K in the Es which appeared at the height of 92.5 km, while outside Te was about 600K. The Ne value changed from 10 4 els/cm 3 to 4×10 4 els/cm 3 . The measurement was done by using a glass-sealed Langmuir probe. The measurement was repeated using a glass-sealed Langmuir probe, and the Te value inside was very close to the neutral temperature Tn , while outside Es , the Te value was several hundreds of degrees higher than the neutral temperature ( Oyama et al. 2008a , b ). The Te values reported by Barjatya (2007) was similar to these results. Theoretical discussion on Te inside Es was conducted by Gleeson & Axford (1967) : Te is about 15K higher than the model neutral temperature, Tn , while outside Es , Te is higher than Tn by about 80 degrees. As such, Te measurement in the sporadic E layer is scarce ( Szuszczewicz & Holmes 1977 ), and we still have not established the general behavior of Te in it.
If we can guess the energetics of the usual E layer from the measurement of Te inside Es which is very close to Tn , and that of Te outside of Es , which is higher than Tn , Te in the lower E region seems to be unequal to Tn . If so, then a question arises regarding high Te in the lower E region. Have we missed a heating mechanism in the theory? Theoretical computations by Gorse et al. (1985) and Paniccia et al. (1986) for laboratory plasma show that vibrationally excited N2 contributes to the heating of thermal electrons. The vibrational temperature of N2 (measured by a sounding rocket) was close to Te ( Kurihara & Oyama 2005 ). When N2 - O2 mixed-gas plasma was irradiated with an extreme ultraviolet (EUV) lamp, which mimics the solar EUV, an increased peak occurred around 0.4 eV ( Oyama et al. 1991 , 2011 ). This peak increases most when the density ratio of N2 to O2 ([ N ( O2 )]/ N [( N2 )] is ~20%. On the other hand, according to Ashihara & Takayanagi (1974) , who took into account the interaction between thermal electrons and vibrationally excited N2 , the heating of thermal electrons seems to be small. Pavlov (1994) draws the same conclusion regarding the heating of thermal electrons.
It was found that a sharp layer of high Te (thickness of ~20 km) appeared at 11 hours ±30 minutes in local time (LT) in winter (January and February) in the midst of the Sq -current focused ( Oyama & Hirao 1979a ) around the altitude of 105 km as Fig. 3 shows. Higher Te was found as the sounding rocket got closer to the Sq focus, and when there were more sunspots. The measurement of energy distribution indicated that a non-Maxwellian distribution of thermal electrons, or a high-energy tail, exists around the Sq focus ( Oyama & Hirao 1979b ), as Fig. 4 shows. The measurement principle of the energy distribution function is described in section 2.3. It was also observed strong plasma irregularities up to an altitude around 160 km ( Oyama & Hirao 1985 ). One mechanism to explain these features is the potential difference between the summer hemisphere and the winter hemisphere. The potential in the winter hemisphere is higher than that in the summer hemisphere ( Fukushima 1994 ), and the potential difference might be concentrated at altitudes < 150 km, where parallel conductivity becomes smaller than other conductivities ( Oyama 2000 ). This might be related to the field-aligned current, which flows from winter to summer hemispheres ( Takeda 1982 ). Efforts to find a heating mechanism are ongoing (Abe 2015, private communication). It is surely worth trying to observe plasma irregularities from the ground by HF radar when the Sq focus is near. It is noted that the location of the Sq focus changes depending on the day and season in mid-latitudes ( Tarpley 1973 ).
PPT Slide
Lager Image
Anomalous heating of thermal electrons around the Sq focus (Sq focus anomaly): (a) Height profiles of electron temperature measured by an electron temperature probe (ETP) show a sharp maximum increase in the midst of the Sq current system at an altitude of 110 km. It is noted that the high Te measured by ETP reflects the existence of a high-energy tail in the energy distribution of thermal electrons. (b) Possible mechanism of the Sq focus anomaly: The potential in the winter hemisphere is higher than that in the summer hemisphere. Electrons that exist below ~150 km are accelerated by an electric field which might exist along a magnetic field line, originally produced by the potential difference. It may be noted that electrons in the summer hemisphere can travel along the magnetic line of force to the winter hemisphere.
PPT Slide
Lager Image
The second-harmonic components of probe current. The secondharmonic component is in semi-logarithmic scale. (a) In the height range between 88.6 km and 124.6 km, the curve shows the existence of a high energy tail. The curve at the height of 98.4 km shows especially high Te as well as a high energy tail. (b) The linear part of the curve is reduced as the height increases from 128.8 km to 250.5 km, which suggests that the electron energy of the thermal range becomes Maxwellian as the height increases.
- 2.2Fregion Ionosphere and Topside
Since the 1980s, not many aeronomy satellites have been in orbit within the F and topside ionosphere ( Schunk & Nagy 1978 ; Krehbiel et al. 1981 ; Sharma et al. 2010 ). Fig. 5 shows an example of the global distribution of Te at ~600 km at longitudes of 210°–285° and 285°–360°, which were observed by the HINOTOPRI sun observation satellite (to download data, visit the NSSDC website). The measurement was done using an electron temperature probe (modification of a DC probe), which is not influenced by electrode contamination ( Oyama & Cheng 2013 ). The difference between the two figures is especially related to morning overshoot of electron temperature in three seasons, which is considered to be due to zonal wind. Zonal wind becomes effective when the magnetic meridional plane tilts like in the longitude band of 285°–360° ( Oyama et al. 1996 ; Watanabe & Oyama 1996 ).
PPT Slide
Lager Image
Example global map of electron temperature at 600 km for three different seasons in the longitude bands of 210°–285° and 285°–360°. The longitude band of 210°–285° is a region where the magnetic meridional plane tilts eastward about 10°, while in the longitude band of 285°–360°, the magnetic meridional plane tilts tremendously. Due to the tilting of the magnetic meridional plane in the longitude band of 285°–360°, zonal wind, as well as meridional wind, influences the plasma behavior.
Recently, a number of satellites (Astrid, CHAMP, DEMETER, Proba-2, and SWARM) have been launched. Although these satellites are not purely dedicated to ionosphere study, all of these satellites carry a DC Langmuir probe ( Holback et al. 2001 ; Lebreton et al. 2006 ; Buchert et al. 2014 ). The data measured by Astrid has shown unexpectedly high Te ( Holback et al. 2001 ), and the Te measured by CHAMP at night is sometimes overestimated. The Te values at night (period 20 to 03 LT) during the March equinox are sometimes about 500K higher than the HINOTORI satellite data indicates ( Su et al. 1995 , 1996 ; Balan et al. 1997 ). We presume that this high temperature observed by CHAMP could be due to the low sampling rate ( Cooke et al. 2003 ) in addition to insufficient amplifier gain. At night, there are only three or four measurements in the electron retarding region, where electron density is low. If we fit the theoretical curve to the small number of observed data points, where some of the data points are near space potential and are not on the fitting curve, the result is that Te is higher than the true value. The Te values measured by the DEMETER satellite are 2–3 times higher than what we think is the true Te ( Kakinami et al. 2013 ). The high temperature is due to contamination of the electrode. The contamination effect appears as hysteresis in the I-V curve of the DC Langmuir probe, which we discuss in more detail in Section 3.1.
The Te in the plasma bubble observed by the HINOTORI satellite ( Oyama et al. 1988 ) shows clear local time dependence. As Fig. 6 shows, Te is lower than, equal to, and higher than that of outside plasma, before midnight, at midnight, and early morning, respectively. In Fig. 6 , there are several values that show positive deviation ( Te is higher than that of ambient plasma) at night. These points are seen in the South Atlantic Geomagnetic Anomaly, which suggests the heating of thermal electrons by the precipitation of high energy particles. One more elevated Te region is over Hawaii. The reason, we assume, is that energetic particles precipitated by the interaction between Very Low Frequency (VLF) and particles in the lower edge of the radiation belt heat thermal electrons. It should be noted that a VLF transmitter facility is located in Hawaii.
PPT Slide
Lager Image
(a) Local time variation of Te in the plasma bubble. The letters A, B, C, D, E , and F denote the data obtained in the longitude bands of 0˚- 60˚, 60˚ - 120 ˚, 120˚ - 180˚, 180˚ - 240˚, 240˚- 300˚, and 300˚ - 360˚, respectively. In (b), the upper panel shows the total number of plasma bubble (grey shade), and that of the bubble inside which Te is higher than that of ambient plasma (black shade). The lower panel shows the locations where Te inside the plasma bubble is higher (triangular), and lower (circle). The size of the triangle and circle shows the degree of the Te difference from that of ambient plasma (after Oyama et al. 1988).
- 2.3 Plasmasphere
Above the height of 3,000 km, the measurement of Te and Ne using a conventional DC probe is only possible at night, and on rare occasions during the daytime in the shadow of a satellite because the number of secondary electrons from the electrode used for measurement is higher than that of ambient plasma ( Rich et al. 1979 ). Secondary electrons from the satellite are mixed with electrons of the ambient plasma, although the energy of secondary electrons is higher than that of thermal electrons in ambient plasma ( Grard 1973 ). This means that the current caused by secondary electrons is dominant on the I-V curve, on which current due to electrons in the ambient plasma is superposed. Past measurement was done at night, or when the electrode was in the shadow of the satellite during the daytime ( Rich et al. 1979 ). As a result of this situation, the amount of data accumulated was small, and no systematic measurement was done before 1989 ( Titheridge 1998 ). To overcome this problem, we applied a second-harmonic method to measure energy distribution function and finally calculated Te (with an instrument called TED, Abe et al. 1990 ).
To measure the energy distribution of thermal/nonthermal electrons, the Druyvesteyn (1930) principle is applied, which states that energy distribution is proportional to the second derivative of the I-V curve, d 2 I/dV 2 . To get the second derivative, we superpose the sinusoidal signal (frequency of 2 kHz) of small amplitude to the probe voltage and pick up the second-harmonic components (4 kHz) from the distorted DC current. The second harmonic component I is proportional to the second derivative of the I-V curve by a relation I ~ a 2 /4 d 2 I/dV 2 , where a is the amplitude of the sinusoidal signal. Reliable measurement of energy distribution can only be obtained using a glass-sealed Langmuir probe ( Oyama & Hirao 1976 ). The system makes it possible to measure Te even at high altitudes except mid-latitudes or high latitudes in the existence of secondary electrons from the electrode due to EUV. The secondary electrons from the satellite itself have much higher energy than thermal electrons, and they do not affect the Te value. The instrument was deployed in the Japanese auroral study satellite “AKEBONO”, and it produced the first systematic Te measurement up to 8,000 km ( Balan et al. 1996a , b ; Titheridge 1998 ; Kutiev et al. 2002). Figs. 7 and 8 show some of the results. As the satellite followed an eccentric orbit, data was accumulated for a long time. Fig. 7 shows the height profile during the day and night, respectively, over the geomagnetic equator. Fig. 8 shows the average local time variation for heights from 1,000 to 8,000 km in two geomagnetic latitude zones, namely, the 15°-30° (left-side) and 0°-15° (right-side) meridian planes.
PPT Slide
Lager Image
Height profile of the electron temperature up to 8,000 km, which was calculated from the second-harmonic components measured by the AKEBONO satellite (Balan et al. 1996a). The second harmonic component, which is picked up from the probe current when a small sinusoidal signal is superimposed on an electrode, provides the second derivative of the I-V curve. The second derivative of the I-V curve still indicates the exponential function from which Te is calculated, if the thermal electrons of the ambient plasma are Maxwellian (after Balan et al. 1996a).
PPT Slide
Lager Image
Plot of altitude/local time of Te at two geomagnetic latitude zones (GLAT 0°–15°, and GLAT 15°–30°), observed by the Akebono satellite. Altitude is indicated by circles every 2,000 km. Note that Te is about 8,000K at 22 LT at 4,000 km for GLAT 15°–30°, while for GLAT 0°–15° at 22 LT, Te at 4,000 km is < 5,000K. This suggests that the magnetic tube at GLAT 15°–30° is larger than that at GLAT 0°–15°, causing slow Te increase in the early morning, and slow decrease of Te at night at GLAT 15°–30˚, where heat capacity is larger than at GLAT 0°-15˚.
As described above, Te is still not properly measured on some occasions despite the long history of Langmuir probe measurement in space. As space missions require a tremendous amount of manpower, energy, time, and money, we should make the maximum effort to obtain accurate and reliable electron temperatures. This attitude is strongly required because some scientists are trying to improve International Reference Ionosphere (IRI) or to make ionosphere models with the belief that the data measured is accurate and reliable ( Truhlik et al. 2009 , 2012 ). Erroneous results can be avoided only by paying special attention to electrodes as well as to the electronics. Next, we present the issues that must be taken into account to achieve reliable measurements.
- 3.1 Contamination of the Electrode Surface
Ne in the ionosphere ranging from about 60km to above, is on the order of 10 2 to 10 6 els/cm 3 , depending on the altitude, local time, and latitude. In the mid-latitude trough, where Ne drops below 10 3 els/cm 3 even at 2,000 km, Te appears to be in the range from 200K to about 3,000K. At higher altitudes above 3,000 km, Te can be higher, showing a gradual increase. Te and Ne are usually obtained by a DC probe. When sweep voltage is applied to an electrode in ionospheric plasma, current will flow into or out from the electrode, depending on electron or ion retardation. The current-voltage characteristic curve should take the same path when the direction of the sweep voltage is changed (increasing and decreasing directions). However, in almost all cases the current-voltage characteristic takes a different path, called hysteresis, as shown in Fig. 2 . This phenomenon has been well known to researchers for a long time ( Van Berkel 1938 ; Wehner & Medicus 1952 ; Winkler et al. 2000 ; Stamate & Ohe 2002 ), even when very few ionosphere scientists were aware of the problem of serious electrode contamination ( Sturges 1973 ; Amatucci et al. 2001 ; Hirt et al. 2001 ; Steigies et al. 2005 ; Barjatya et al. 2007 ). The hysteresis decreases as the frequency increases, and disappears at higher frequency ranges (frequency dependence); hysteresis also increases at higher densities (density dependence). The mechanism of hysteresis is thought to be related to contaminants on the electrode surface, whose main constituent is water ( Oyama 1976 ). The combination of the large dielectric constant of water (~80) and the thin contaminant layer in the range of 20-80 nm provides large capacitance of the order of 10-100 μF. The contamination layer is a parallel connection of the capacitor and a resister. The resistance is about 500 kΩ - 1 MΩ for the case of a cylindrical electrode 3mm in diameter and 20 cm in length.
To avoid undesirable effects of the contaminated layer on the measurement, Oyama & Hirao (1976) developed a way to keep the probe surface clean by sealing the cylindrical electrode within a glass tube until starting measurement. This device is called the glass-sealed cylindrical probe. Another way to remove the effect of electrode contamination is to sweep the DC probe voltage with a frequency high enough to remove the hysteresis. Usually a frequency of ~10 Hz is high enough. When the sweep voltage is changed quickly, a correct I-V (I in the vertical axis, and V in the horizontal axis) curve is obtained (the hysteresis disappears), although the horizontal axis shifts to a more positive direction. The Te values obtained by a fast sweep voltage agree with values obtained by a clean probe. A third way is to modulate a triangular wave by a short pulse, called a “pulse probe”. Its design was based on an equivalent circuit model. Because the pulse contains a high frequency component, the contamination effect disappears ( Szuszczewicz & Holmes 1975 ). However, this method cannot be used when we try to measure the energy distribution function using the second-harmonic method ( Boyd & Twiddy 1959 ). Instead, a glass-sealed Langmuir probe was used to measure energy distribution in the E region ( Oyama & Hirao 1979a , b ; Oyama et al. 1980 ; Oyama & Hirao 1985 ).
There is one more way, namely, to use an electron temperature probe (a modification of a DC Langmuir probe, see Hirao & Oyama 1970 ). In this approach, a low frequency sinusoidal signal is applied to the electrode, and the floating potential shift is picked up due to the nonlinearity of the I-V curve. The value can be measured from the shift. The advantages of using this probe are that the instrument is not influenced by electrode contamination ( Oyama & Cheng 2013 ), it can be manufactured cheaply, and it consumes very little power. The probe was developed in the 1970s ( Hirao & Oyama 1970 ) and was installed in five scientific satellites in Japan as well as more than 60 rockets, including rockets deployed in the Antarctic, US, India, Brazil, and West Germany. The probe was also carried in satellites, including those of Korea, Russia, and Brazil. Recently, the fixed-frequency oscillator was replaced with a varying frequency circuit, which changes the frequency from 500 kHz to 10 MHz. This allows both Te and Ne to be measured successively with one simple instrument ( Oyama et al. 2015 ).
- 3.2 Surface Area of Satellite
Operation of the DC Langmuir probe is based on the assumption that a counter electrode is infinite, compared with the surface area of the electrode. As satellites have become smaller, the conductive surface area of the satellite used as a counter electrode has become insufficient. As the voltage of the electrode increases with respect to the satellite frame, currents due to thermal electrons of ambient plasma increase. Accordingly, the satellite potential goes further negative to collect more ions from the ambient plasma, trying to cancel the electron current flowing out from the electrode. As a result, the I-V curve is not adequate to measure either Te or Ne . To overcome this problem, especially for small satellites, several possible solutions have been considered. For example, a satellite can be equipped with an electron transmitter, which will emit thermal electrons when the satellite potential becomes more negative, and as a result the satellite potential recovers. Although double-probe measurement ( Johnson & Malter 1950 ) is free from the problem of conductive surface, problems associated with electrode contamination still remain. Another approach is to provide an instrument for which measurement is not influenced by a small conductive surface ( Oyama et al. 2015 ). This instrument combines a small inlet electrode with an electron multiplier, which effectively increases the conductive surface area of the counter electrode ( Shimoyama et al. 2011 ).
- 3.3 Stray Capacitance Effects and Others
When we want to increase the spatial resolution of Te measurement, such as in a plasma bubble, we need to increase the frequency of the sweep voltage of the electrode. For 10 Hz and 100 Hz of the frequency, the horizontal space resolutions are ~ 700 and 70 m, respectively. However, as the frequency increases, the effect of stray capacitance around the first preamplifier becomes more serious. The sweep voltage is differentiated by a stray capacitor, and a square wave signal is obtained at the output voltage even when the electrode is not in a plasma environment. The amplitude of the square wave increases as the frequency of the sweep voltage increases, and it is much greater than that due to the probe current. This feature is overcome by taking several steps. To reduce stray capacitance, the same voltage should be applied to the shielding line of the coaxial cable that connects the electrode and the electronics box as that used as the sweep voltage. At the same time, a small capacitor should be added to cancel the imbalance in the frequency response of the preamplifiers. To get accurate and reliable values, laboratory experiments should be repeated so that information on the electronics and their plasma measurement is accumulated. It is stressed that every scientist should be aware of the fundamental problems involved in obtaining reliable data under close communication with the engineers who are in charge of measurement systems.
- 4.2 F region and Top Side Ionosphere
Measurements of Te in this region are still scarce, although since 1978 when Schunk & Nagy (1978) reviewed Te measurements in the F region, our knowledge has gradually increased. Measurements of Te and Ne , especially Te , might provide information on neutral wind systems with the help of computer simulation ( Su et al. 1995 , 1996 ).
Measurement of Te inside a plasma bubble has so far been difficult because high-resolution measurement is required, and an accordingly large down link is needed. So far, only three papers have reported electron temperature in a plasma bubble ( Oyama et al. 1988 ; Oyama et al. 1993 ; Park et al. 2008 ). The measurement of Te at the plasma edge might provide a further interesting research subject regarding the current flowing at the edge of a plasma bubble. Recently, Huba et al. (2011) , conducted a computer simulation using the SAM13/ESF simulation code. The result shows that cooling occurs just below a plasma bubble at the equator, and heating occurs at the foot of each hemisphere. Taking these new findings into account, further high space resolution observation by satellites and analysis of data which have been measured so far should be conducted.
The effect of sudden stratosphere warming (SSW) on the ionosphere has recently been studied, mainly using Total Electron Content (TEC) data ( Jin et al. 2012 ; Oyama et al. 2014 ). Although SSW is a high latitude phenomenon that occurs around 10 hPa height, the ionosphere disturbance extends globally in longitude and latitude, maybe due to changes in the wind system of the thermosphere as well as the dynamo region atmosphere, which then changes the distribution of the electric field. Accordingly, Ne and Te satellite data ( Lühr et al. 2012 ) are needed to study the SSW effect on a global scale, which TEC data does not provide. Data on Te in particular, might help us to understand the mechanism of ionosphere modification, because the neutral wind in the thermosphere and lower E region seems to play an important role in modifying the electric field in the equatorial region as well.
Recently, some scientists have started using Te and Ne satellite data to study atmospheric features preceding large earthquakes. Oyama et al. (2008c) reported the disappearance of afternoon overshoot in conjunction with three earthquakes. Such overshoot is usually found in the afternoon at low and middle latitudes. Ryu et al. (2014) reported that the equatorial ionosphere is influenced by large earthquakes (M > 6.5). The area disturbed by large earthquakes extends roughly 40° to the west and east from the epicenter (the total region is 80°), and 20° to the north and south (total 40°). Depending on the latitude, local time, and height, the features of ionospheric disturbance differ ( Oyama et al. 1996 ); therefore, it is difficult to cover the entire area affected using ground-based experiments. Although many reports on the ionosphere precursor have appeared, we have not found a way to predict a large earthquake in the short-term (about 5-10 days). We believe that a constellation of satellites is needed to cover the whole globe, to establish the morphology of ionosphere features as we describe in detail in section 5.
- 4.3 Plasmasphere
We still do not have sufficient data regarding the downward flow of plasma that fills the inside of flux tubes during the daytime and flows down along magnetic line of force at night. The efficiency by which thermal electrons are heated by photoelectrons traveling along magnetic lines of force from the ionosphere below has not yet been measured ( Balan et al. 1996b ). There have been many theoretical papers on the heating of thermal electrons by electrons escaping from ionosphere; however, observation of these is rare and insufficient to study these issues. To do so, we first need to know the electron temperature from the F region upward to the plasmasphere, but systematic measurements are still rarely conducted. Inside the trough above 3,000 km is still not known because Ne drops to 10 2 els/cm 3 , making measurement by DC probe difficult. For measurement in low density we may combine a small retarding potential analyzer and an electron multiplier that operates by current measurement mode, not by particle counting mode ( Shimoyama et al. 2011 ). The second-harmonic method is applied to the instrument, which removes the problem of the secondary electrons from the electrode and from the satellite wall.
There are still many issues that are not clear, and some are totally unknown. However to answer many questions, we need sounding rockets as well as satellites. We are now rapidly moving into an era of small satellites. Using a DC Langmuir probe or some modification of it on board tiny satellites (e.g., 10-50 kg total mass), we might be able to find answers to the questions discussed. It may even become possible, in the short term, to contribute to human safety by the study of atmospheric indicators of earthquakes. The tiny satellites (10-50 kg in total weight) should be in the regime of constellation missions, and the altitude of each satellite should be 250-300 km with a circular orbit and an inclination of about 40°. Although the launching cost would still be expensive, the cost could be shared among countries that would like to increase our knowledge of space, and to contribute to human society by study of the ionosphere/plasmasphere ( Oyama et al. 2010 ). DC probes might be carried in orbiters of other planets, as planetary exploration evolves ( Krehbiel et al. 1980 ). The experience we have acquired in earth orbiting satellites surely constitutes a foundation for successful planetary missions.
The manuscript was completed during the author’s time at the Korean Advanced Institute of Science and Technology (KAIST). He expresses his since thanks to Professor K. W. Min of KAIST for his invitation to do so. He also expresses his sincere thanks to the Mitsubishi Foundation, Japan for supporting a project of earthquake precursor study. Part of the information onTemeasurement inEswas provided by Dr. T. Abe, Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency.
Abe T , Oyama KI , Amemiya H , Watanabe S , Okuzawa T 1990 Measurements of temperature and velocity distribution of thermal electrons by the Akebono (EXOS-D) satellite J. Geomag. Geolectr 42 537 - 554    DOI : 10.5636/jgg.42.537
Amatucci WE , Schuck PW , Walker DN , Kintner PM , Powell S 2001 Contamination-free sounding rocket Langmuir probe Rev. Sci. Instrum Contamination-free sounding rocket Langmuir probe 72 2052 - 2057    DOI : 10.1063/1.1357234
Andreyeva LA , Burakov YB , Katasev LA , Komrakov GP , Nesterov VP 1971 Rocket investigations of the ionosphere at mid-latitudes Space Res 11 1043 - 1050
Ashihara O , Takayanagi K 1974 Velocity distribution of ionospheric low-energy electrons Planet. Space Sci 22 1201 - 1217    DOI : 10.1016/0032-0633(74)90005-1
Aubry M , Blanc M , Clauvel R , Taieb C , Bowen PJ 1966 Some rocket results on sporadic E Radio Sci 1 170 - 177
Balan N , Oyama KI , Bailey GJ , Abe T 1996 Plasmasphere electron temperature studies using satellite observations and a theoretical model J. Geophys. Res. 101 15323 - 15330    DOI : 10.1029/96JA00823
Balan N , Oyama KI , Bailey GJ , Abe T 1996 Plasmasphere electron temperature profiles and the effects of photoelectron trapping and an equatorial high-altitude heat source J. Geophys. Res. 101 21689 - 21696    DOI : 10.1029/96JA01798
Balan N , Oyama KI , Bailey GJ , Fukao S , Watanabe S 1997 A plasma temperature anomaly in the equatorial ionosphere J. Geophys. Res 102 7485 - 7492    DOI : 10.1029/97JA00039
Barjatya A 2007 Langmuir probe measurements in the ionosphere, PhD Dissertation Utah State University
Boggess RL , Brace LH , Spencer NW 1959 Langmuir probe measurements in the ionosphere J. Geophys. Res 64 1627 - 1630    DOI : 10.1029/JZ064i010p01627
Boyd RLF , Twiddy ND 1959 Electron energy distribution in plasmas. I Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences 250 53 - 69    DOI : 10.1098/rspa.1959.0050
Brace LH , Pfaff RF , Borovsky JE , Young DT 1998 Measurement Techniques in Space Plasmas: Particles, Geophysical Monograph American Geophysical Union;AGU Washington, D.C Langmuir probe measurements in the ionosphere 23 - 35
Brace LH , Theis RF , Dalgarno A 1973 The cylindrical electrostatic probes for Atmosphere Explorer -C, -D, and -E Radio Sci 8 341 - 348    DOI : 10.1029/RS008i004p00341
Buchert SC , Eriksson A , Gill R , Nilsson T , Åhlen L First electron density and temperature estimates from the Swarm Langmuir probes and a comparison with IS measurements 2014 EGU General Assembly meeting Vienna, Austria 27 Apr-02 May 2014
Buchert S , Zangerl F , Sust M , André M , Eriksson A 2015 SWARM observations of equatorial electron densities and topside GPS track losses Geophys. Res. Lett. 42 2088 - 2092    DOI : 10.1002/2015GL063121
Cooke DL , Roth C , Luehr H , Jakowski N In situ plasma observation at 400 km on the CHAMP satellite EGS-AGUEUG Joint Assembly meeting Nice, France 6-11 Apr 2003
Costa E , Roddy PA , Ballenthin JO 2014 Statistical analysis of C/NOFS planar Langmuir probe data Ann. Geophys 32 773 - 791    DOI : 10.5194/angeo-32-773-2014
Druyvesteyn MJ 1930 Der Niedervolt Bogen Z. Physik 64 781 - 798    DOI : 10.1007/BF01773007
Fukushima N 1994 Some topics and historical episodes in geomagnetism and aeronomy J. Geophys. Res 99 19113 - 19142    DOI : 10.1029/94JA00102
Gleeson LJ , Axford WI 1967 Electron and ion temperature variations in temperature zone sporadic-E layers Planet. Space Sci. 15 123 - 136    DOI : 10.1016/0032-0633(67)90071-2
Gorse C , Capitelli M , Ricard A 1985 On the coupling of electron and vibrational energy distributions in H2, N2, and CO post discharges J. Chem. Phys 82 1900 - 1906    DOI : 10.1063/1.448374
Grard RJL 1973 Properties of the satellite photoelectron sheath derived from photoemission laboratory measurements J. Geophys. Res. 78 2885 - 2906    DOI : 10.1029/JA078i016p02885
Hirao K , Oyama KI 1970 An improved type of electron temperature probe J. Geomagn. Geoelectr. 22 393 - 402    DOI : 10.5636/jgg.22.393
Hirao K , Oyama KI 1972 A critical study on the reliability of electron temperature measurements with a Langmuir probe J. Geomagn. Geoelectr 24 415 - 427    DOI : 10.5636/jgg.24.415
Hirt M , Steigies CT , Piel A 2001 Plasma diagnostics with Langmuir probes in the equatorial ionosphere: II. Evaluation of DEOS flight F06 J. Phys. D: Appl. Phys. 34 2650 - 2657    DOI : 10.1088/0022-3727/34/17/312
Holback B , Jacksén Å , Åhlén L , Jansson S-E , Eriksson AI 2001 LINDA - the Astrid-2 Langmuir probe instrument Ann. Geophys 19 601 - 610    DOI : 10.5194/angeo-19-601-2001
Huba JD , Joyce G , Krall J 2011 Plasma bubbles in the post-sunset ionosphere, NRL Review, Space research and satellite technology 224 - 225
Jin H , Miyoshi Y , Pancheva D , Mukhatarov P , Fujiwara H 2012 Response of migrating tides to the stratospheric sudden warming in 2009 and their effects on the ionosphere studied by a whole atmosphere-ionosphere model GAIA with COSMIC and TIMED/SABER observations J. Geophys. Res 117 A10323 -    DOI : 10.1029/2012JA017650
Johnson EO , Malter L 1950 A floating double probe method for measurements in Gas discharges Phys. Rev 80 58 - 68    DOI : 10.1103/PhysRev.80.58
Kakinami Y , Kamogawa M , Onishi T , Mochizuki K , Lebreton JP 2013 Validation of electron density and temperature observed by DEMETER Adv. Space Res. 52 1267 - 1273    DOI : 10.1016/j.asr.2013.07.003
Krehbiel JP , Brace LH , Theis RF , Cutler JR , Pinkus WH 1980 Pioneer Venus Orbiter Electron temperature probe IEEE Trans. Geosci. Remote Sens GE-18 49 - 54    DOI : 10.1109/TGRS.1980.350260
Krehbiel JP , Brace LH , Theis RF , Pinkus WH , Kaplan RB 1981 The Dynamics Explorer Langmuir probe instrument Space Sci. Instrum 5 493 - 502
Kurihara J , Oyama KI 2005 Rocket-borne instrument for measuring vibrational-rotational temperature and density in the lower thermosphere Rev. Sci. Instrum. 76 083101 -    DOI : 10.1063/1.1988189
Kurihara J , Oyama KI , Iwagami N , Takahashi T 2006 Numerical simulation of 3-D flow around sounding rocket in the lower thermosphere Ann. Geophys. 24 89 - 95    DOI : 10.5194/angeo-24-89-2006
Kutiev I , Oyama KI , Abe T 2002 Analytical representation of the plasmasphere electron temperature distribution based on Akebono data J. Geophys. Res. 107 1459 -    DOI : 10.1029/2002JA009494
Langmuir I , Mott-Smith H 1924 Studies of electric discharges in gas at low pressures Part I-V Gen. Electr. Rev. 27 449 - 455, 538-548, 616-623, 762-771, 810-820
Lebreton J-P , Stverak S , Travnicek P , Maksimovic M , Klinge D 2006 The ISL Langmuir probe experiment processing onboard DEMETER: Scientific objectives, description and first results Planet. Space Sci. 54 472 - 486    DOI : 10.1016/j.pss.2005.10.017
Lühr H , Park J , Ritter P , Liu H 2012 In-situ CHAMP observations of the ionosphere-thermosphere coupling Space Sci. Rev. 168 237 - 260    DOI : 10.1007/s11214-011-9798-4
Oyama KI 1976 A systematic investigation of several phenomena associated with contaminated Langmuir probes Planet. Space Sci. 24 183 - 190    DOI : 10.1016/0032-0633(76)90104-5
Oyama KI 2000 In situ measurements of Te in the lower ionosphere- A review Adv. Space Res. 26 1231 - 1240    DOI : 10.1016/S0273-1177(99)01220-X
Oyama KI , Cheng CZ , Oyama Cheng 2013 An Introduction to Space Instrumentation Terra Pub Tokyo 91 - 105
Oyama KI , Hirao K 1976 Application of grass-sealed Langmuir probe to ionosphere study Rev. Sci. Instrum 47 101 - 107    DOI : 10.1063/1.1134476
Oyama KI , Hirao K 1979 Anomalous heating of the thermal electrons near the focus of the Sq current vortex J. Geomagn. Geoelectr 31 11 - 19    DOI : 10.5636/jgg.31.11
Oyama KI , Hirao K 1979 Distortion of the energy distribution of ionospheric thermal electrons near the focus of Sq current vortex Planet. Space Sci. 27 183 - 192    DOI : 10.1016/0032-0633(79)90049-7
Oyama KI , Hirao K 1985 Energy gain of supra-thermal electrons from the excited neutral gases J. Geomagn. Geoelectr 37 913 - 926    DOI : 10.5636/jgg.37.913
Oyama KI , Hirao K , Banks PM , Williamson PR 1980 Is Te equal to Tn at the height of 100 to 120 km? Planet. Space Sci. 28 207 - 211    DOI : 10.1016/0032-0633(80)90013-6
Oyama KI , Schlegel K , Watanabe S 1988 Temperature structure of plasma bubbles in the low latitude ionosphere around 600km altitude Planet. Space Sci. 36 553 - 567    DOI : 10.1016/0032-0633(88)90025-6
Oyama KI , Suzuki K , Kawashima N , Zalpuri KS , Teii S 1991 An extreme ultraviolet radiation source for the simulation of the ionosphere Rev. Sci. Instrum. 62 1721 - 1726    DOI : 10.1063/1.1142412
Oyama KI , Watanabe S , Oya H 1993 Energetics in the plasma bubble Adv. Space Res 13 293 - 297    DOI : 10.1016/0273-1177(93)90025-7
Oyama KI , Watanabe S , Su Y , Takahashi T , Hirao K 1996 Season, local time, and longitude variations of electron temperature at the height of -600 km in the low latitude region Adv. Space Res 18 269 - 278    DOI : 10.1016/0273-1177(95)00936-1
Oyama KI , Abe T , Mori H , Liu JY 2008 Electron temperature in nighttime sporadic E layer at mid-latitude Ann. Geophys. 26 533 - 541    DOI : 10.5194/angeo-26-533-2008
Oyama KI , Hibino K , Abe T , Pfaff R , Yokoyama T 2008 Energetics and structure of the lower E region associated with sporadic E layer Ann. Geophys. 26 2929 - 2936    DOI : 10.5194/angeo-26-2929-2008
Oyama KI , Kakinami Y , Liu JY , Kamogawa M , Kodama T 2008 Reduction of electron temperature in low-latitude ionosphere at 600 km before and after large earthquakes J. Geophys. Res. 113 A113 -    DOI : 10.1029/2008JA013367
Oyama KI , Kakinami Y , Liu JY , Chen CY , Kodama T 2010 Micro/mini satellites for earthquake studies - toward international collaboration Adv. Geosci 21 251 - 256
Oyama KI , Shimoyama M , Liu JY , Cheng CZ 2011 Possible interaction between thermal electrons and vibrationally excited N2 in the lower E-region Ann. Geophys. 29 583 - 590    DOI : 10.5194/angeo-29-583-2011
Oyama KI , Lee CH , Fang HK , Cheng CZ 2012 Means to remove electrode contamination effect of Langmuir probe measurement in space Rev. Sci. Instrum. 83 055113 -    DOI : 10.1063/1.4722167
Oyama KI , Jhou JT , Lin JT , Lin C , Liu H 2014 Ionospheric response to 2009 sudden stratospheric Warming in the Northern Hemisphere J. Geophys. Res. 119 10260 - 10275    DOI : 10.1002/2014JA020014
Oyama KI , Hsu YW , Jiang GS , Chen WH , Fang HK 2015 Electron temperature and density probe for small aeronomy satellites Rev. Sci. Instrum submitted
Paniccia F , Gorse C , Bretagne J , Capitelli M 1986 Electron energy distribution functions in molecular nitrogen: The role of superelastic electronic collisions in discharge and postdischarge conditions J. Appl. Phys 59 4004 - 4006    DOI : 10.1063/1.337020
Park J , Min KW , Kim VP , Kil H , Su S-Y 2008 Equatorial plasma bubbles with enhanced ion and electron temperatures J. Geophys. Res 113 A09318 -    DOI : 10.1029/2008JA013067
Pavlov AV 1994 The role of vibrationally excited oxygen and nitrogen in the D and E regions of the ionosphere Ann. Geophys 12 1085 - 1090    DOI : 10.1007/s00585-994-1085-1
Piel A , Rohde V , Thiemann H , Oyama KI 1992 Plasma drifts deduced from resonance cone asymmetries: II. Evaluation of COREX data Adv. Space Res 12 263 - 266    DOI : 10.1016/0273-1177(92)90069-A
Reifman A , Dow WG 1949 Dynamic probe measurements in the ionosphere Phys. Rev. 76 987 - 988    DOI : 10.1103/PhysRev.76.987
Rich FJ , Sagalyn RC , Wildman PJL 1979 Electron temperature profiles measured up to 8000 km by S3-3 in the late afternoon sector J. Geophys. Res. 84 1328 - 1332    DOI : 10.1029/JA084iA04p01328
Rohde V , Piel A , Thiemann H , Oyama KI 1993 In situ diagnostics of ionospheric plasma with the resonance cone technique J. Geophys. Res 98 19163 - 19172    DOI : 10.1029/93JA01920
Ryu K , Lee E , Chae JS , Parrot M , Pulinets S 2014 Seismo-ionospheric coupling appearing as equatorial electron density enhancements observed via DEMETER electron density measurements J. Geophys. Res. 119 8524 - 8542    DOI : 10.1002/2014JA020284
Sampath S , Sastry TSG , Oyama KI , Hirao K 1974 Joule heating due to the equatorial electrojet as observed by rocket borne probes Space. Res 14 253 - 258
Schunk RW , Nagy AF 1978 Electron temperatures in the F region of the ionosphere: theory and observations Rev. Geophys 16 355 - 399    DOI : 10.1029/RG016i003p00355
Schutz SR , Smith LG 1976 Electron temperature measurements in mid-latitude sporadic E layers J. Geophys. Res 81 3214 - 3220    DOI : 10.1029/JA081i019p03214
Seran E , Berthelier JJ , Saouri FZ , Lebreton J-P 2005 The spherical segmented Langmuir probe in a flowing thermal plasma: numerical model of the current collection Ann. Geophys. 23 1723 - 1733    DOI : 10.5194/angeo-23-1723-2005
Sharma PK , Pathak PP , Sharma DK , Rai J 2010 Variation of ionospheric electron and ion temperatures during periods of minimum to maximum solar activity by the SROSS-C2 satellite over Indian low and equatorial latitudes Adv. Space. Res. 45 294 - 302    DOI : 10.1016/j.asr.2009.08.001
Shimoyama M , Oyama KI , Abe T , Yau AW 2011 Suprathermal plasma analyzer for the measurement of low-energy electron distribution in the ionosphere Rev. Sci. Instrum. 82 074501 -    DOI : 10.1063/1.3606450
Smith LG , Weeks LH , McKinnon PJ 1968 Rocket observations of electron temperature in the E-region J. Atmos. Terr. Phys 30 1301 - 1312    DOI : 10.1016/S0021-9169(68)91141-0
Smith LG , Zimmerman RK , Hirao K , Oyama KI , Calderon C 1978 Electron temperature in the equatorial E region measured by two rocket experiments and by incoherent scatter Space Res 18 265 - 268
Spencer NW , Brace LH , Carignan GR , Taeusch DR , Niemann H 1965 Electron and molecular nitrogen temperature and density in the thermosphere J. Geophys. Res. 70 2665 - 2698    DOI : 10.1029/JZ070i011p02665
Stamate E , Ohe K 2002 Influence of surface condition in Langmuir probe measurements J. Vac. Sci. Technol. A 20 661 - 666    DOI : 10.1116/1.1460893
Steigies CT , Jr. Pfaff RF , Rowland DE Proceedings of the 17th ESA symposium on European Rocket and Balloon Programmes and Related Research Sandefjord, Norway 30 May - 2 Jun 2005
Sturges DJ 1973 An evaluation of ionospheric probe performance- I. Evidence of contamination and clean-up of probe surfaces Planet. Space. Sci. 21 1029 - 1047    DOI : 10.1016/0032-0633(73)90148-7
Su YZ , Oyama KI , Bailey GJ , Takahahsi T , Watanabe S 1995 Comparison of the satellite electron density and temperature measurements at low latitudes with a plasmasphere-ionosphere model J. Geophys. Res. 100 14591 - 14604    DOI : 10.1029/95JA00682
Su YZ , Oyama KI , Bailey GJ , Fukao S , Takahashi T 1996 Longitudinal variations of the topside ionosphere at low latitudes: Satellite measurements and mathematical modelings J. Geophys. Res 101 17191 - 17205    DOI : 10.1029/96JA01315
Szuszczewicz EP , Holmes JC 1975 Surface contamination of active electrodes in plasmas: Distortion of conventional Langmuir probe measurements J. Appl. Phys. 46 5134 - 5139    DOI : 10.1063/1.321572
Szuszczewicz EP , Holmes JC 1977 Observations of electron temperature gradients in mid-latitude Es layers J. Geophys. Res 82 5073 - 5080    DOI : 10.1029/JA082i032p05073
Takeda M 1982 Three dimensional ionospheric currents and field aligned currents generated by asymmetrical dynamo action in the ionosphere J. Atmos. Terr. Phys. 44 187 - 193    DOI : 10.1016/0021-9169(82)90122-2
Tarpley JD 1973 Seasonal movement of the Sq current foci and related effects in the equatorial electrojet J. Atmos. Terr. Phys 35 1063 - 1071    DOI : 10.1016/0021-9169(73)90005-6
Titheridge JE 1998 Temperatures in the upper ionosphere and plasmasphere J. Geophys. Res. 103 2261 - 2278    DOI : 10.1029/97JA03031
Truhlik V , Bilitza D , Triskova L 2009 Latitudinal variation of the topside electron temperature at different levels of solar activity Adv. Space Res 44 693 - 700    DOI : 10.1016/j.asr.2009.04.029
Truhlik V , Bilitza D , Triskova L 2012 A new global empirical model of the electron temperature with the inclusion of the solar activity variations for IRI Earth Planets Space 64 531 - 543    DOI : 10.5047/eps.2011.10.016
Van Berkel WPJ 1938 Einflusz von anderungen des sondenzustandes auf sondencharakteristiken nach Langmuir Physica 5 230 - 240    DOI : 10.1016/S0031-8914(38)80080-X
Watanabe S , Oyama KI 1996 Effects of neutral wind on the electron temperature at a height of 600 km in the low latitude region Ann. Geophys 14 290 - 296    DOI : 10.1007/s00585-996-0290-5
Wehner G , Medicus G 1952 Reliability of probe measurements in hot cathode gas diodes J. Appl. Phys 23 1035 - 1046    DOI : 10.1063/1.1702342
Winkler C , Strele D , Tscholl S , Schrittwieser R 2000 On the contamination of Langmuir probe surfaces in a potassium plasma Plasma Phys. Control. Fusion 42 217 - 223    DOI : 10.1088/0741-3335/42/2/311
Yamamoto M , Fukao S , Tsunoda RT , Pfaff R , Hayakawa H 2005 SEEK-2 (Sporadic-E experiment over Kyushu 2) - Project outline and significance Ann. Geophys. 23 2295 - 2305    DOI : 10.5194/angeo-23-2295-2005
Yoshimura R , Imamura T , Abe T 1999 Measurement of Ts in the Es layer using the DC probe onboard the rocket S-310-27, In Japanese, Uchu-kagaku-kenkyusho-houkoku 1 - 33