Radar studies of mid-latitude ionospheric plasma drifts L. Scherliess, B. G. Fejer, J. Holt, L. Goncharenko, C. Amory-Mazaudier, M. J. Buonsanto To cite this version: L. Scherliess, B. G. Fejer, J. Holt, L. Goncharenko, C. Amory-Mazaudier, et al.. Radar studies of mid-latitude ionospheric plasma drifts. Journal of Geophysical Research, American Geophysical Union (AGU), 2001, pp.1771-1783. <hal-00979238> HAL Id: hal-00979238 http://hal.upmc.fr/hal-00979238 Submitted on 15 Apr 2014 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L’archive ouverte pluridisciplinaire HAL, est destin´ee au d´epˆot et `a la diffusion de documents scientifiques de niveau recherche, publi´es ou non, ´emanant des ´etablissements d’enseignement et de recherche fran¸cais ou ´etrangers, des laboratoires publics ou priv´es. JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 106, NO. A2, PAGES 1771-1783, FEBRUARY Radar studies of midlatitude 1, 2001 ionospheric plasma drifts L. Scherliess •, B. G. Fejer•, J. Holt2, L. Goncharenko 2, C. Amory-Mazaudier a, and M. J. Buonsanto 2,4 Abstract. We use incoherent scatter radar measurements from Millstone Hill and Saint Santin to study the midlatitude F regionelectrodynamicplasma drifts during geomagneticallyquiet and active periods. We present initially a local time, season, and solar flux dependent analytical model of the quiet time zonal and meridional ExB drifts over these stations. We discuss,for the first time, the Saint Santin drift patterns during solar maximum. We have used these quiet time models to extract the geomagneticperturbation drifts which were modeled •s • function of the time history of the auroral electrojet indices. Our results illustrate the evolution of the disturbance drifts driven by the combinedeffectsof prompt penetration •nd longer lasting perturbation electric fields. The meridion•l electrodynamic disturbance drifts have largest amplitudes in the midnight-noon sector. The zonal drifts •re predominantly westward, with largest amplitudes in the dusk-midnight sector •nd, followinga decreasein the high-latitude convection,they dec•y more slowlyth•n the meridional drifts. The prompt penetration •nd steady state zonal disturbance drifts derived from radar measurements•re in good agreement with results obtained from both the ion drift meter data on board the DynamicsExplorer 2 (DE 2) satellite •nd f?om the Rice Convection 1. Introduction Model. The midlatitude ionospheric plasma drifts can be simultaneously driven by several electric field sources. The two most important ones are the ionosphericwind dynamo and the direct penetration of high-latitude convection electric fields into the plasmasphere. The ionospheric wind dynamo generateselectric fields through the dynamo action of the thermospheric wind circulation driven by solar heating and during geomagneticactive periods, also by thermospheric disturbance winds driven by enhanced Joule heating in the auroral re- Over the last two decades, incoherent scatter radar observationshave been used extensively to study the midlatitude ionosphericelectrodynamic(ExB) plasma drifts during geomagnetically quiet and disturbedperiods. Blanc and Amayenc [1979] have used incoherent scatter radar observationsduring geomagnetically quiet days over Saint Santin during the 1973- 1975 solar minimum period to derive an analytical seasonal model of their quiet time F region plasma drifts. This model included a steady drift componentand the first four diurnal harmonicoscillations. Wand [1981] has used Millstone Hill observations during May 1976 to November 1977 to analytically model the upper mid1987; Mazaudier and Venkateswaran, 1990; Scherliess latitude solar minimum electrodynamic plasma drifts gions [Blanc and Richmond, 1980; Mazaudier et al., and Fejer, 1997]. Prompt penetrationelectric fields for geomagneticallyquiet and disturbed periods. The can affect the midlatitude plasma drifts nearly instan- disturbance taneouslyduring transient periodsof insufficientshield- Kp index. Buonsantoet al. [1993]haveuseda significantly larger databaseof Millstone Hill plasma drift ing by plasmasheetcurrents[Jaggiand Wolf,19731and drifts were modeled as a function of the also during periodsof strongsteady auroral convection. observationscomprising73 experimentsfrom February 1984 to February 1992 and determined the averageseasonal quiet time plasma drift patterns for both solar •Center for Atmospheric and Space Sciences,Utah State University, Logan. 2HaystackObservatory,MassachusettsInstitute of Technology, Westford. aCentre d'Etude des Environnements Terrestre minimum and maximum conditions. The characteris- tics of equatorial, low-latitude, and midlatitude plasma drifts measured with incoherent scatter radar were re- viewedby Richmond[1995]. More recently,Buonsanto and Witasse[1999]presentedan updatedstudy of the local time, season,solar cycle, and geomagneticactivity dependent climatologiesof the Millstone Hill F region plasma drifts and thermosphericwinds. In this study, magneticactivity effectswere again accountedfor as a et Plane- taires, Saint-Maur-de-Fosses, France. 4DeceasedOctober 20, 1999. Copyright2001 by the AmericanGeophysicalUnion. Paper nulnbcr 2000JA000229. 0148-0227/01/2000JA000229509.00 function of the Kp index. 1771 1772 SCHERLIESS ET AL.' MIDLATITUDE A detailed study of the middle- and low-latitude F region zonal plasma drifts measured by the Ion Drift Meter (IDM) on boardthe polar orbitingDynamicsExplorer2 (DE 2) satelliteduringthe 1981-1983solarmax- PLASMA DRIFTS three directions and with the assumption of a uniform velocityfield [e.g., Buonsantoet al., 1993; Buonsanto and Witasse,1999]. The measurements in our database correspondto an altitude of 300 kin, have a typical er- imum periodwaspresentedby Heelisand Coley[1992]. ror of •08 m/s, and a cycletime of •-l hour. We will use They derived the average latitudinal plasma drift pat- here drift componentsperpendicular to the Earth magterns during geomagneticallyquiet and disturbed peri- netic field. In the F region over Millstone Hill, a plasma ods correspondingto Kp _(2.0 and Kp_)3.0, respectively. velocityof 22 m/s corresponds to an electricfield of 1 They could not derive the seasonaldependenceof these drifts since the DE 2 measurements cal time had season and lo- interlocked. Recently,Fejer and $cherliess[1995]haveintroduced The Saint Santin incoherentscatterradar site (45øN, 2øE, 40ø magnetic latitude) operated originally on a bistatic mode with a transmitter in Saint Santin and a new methodology which was able to separate,for the a receiver at a distance of 302 km in Nancay. In 1973, first time, the effects of different electric field sources two additional receiverswere installed in Monpazier and on the plasma drifts during geomagnetic active peri- Mende (•-100 km from the transmittersite), allowing ods. $cherliessand Fejer [1998]and Fejer and $cher- for the determination of three-dimensional ion drift vecliess[1998]have usedthis methodologyto analytically tors at a typical heightof •0300km [Baueret al., 1974]. errorwasof the orderof 5-10m/s and model the temporal and latitudinal variations of the The measurement low- and middle-latitude F region zonal perturbation the integration time was •030 min. The Saint Santin drifts measured by the DE 2 satellite during geomagnetically disturbed periods as a function of the time history of the hourly averaged AE indices. They determined the drift patterns owing to transient prompt penetration electric fields and of longer lasting distur- radar operated until 1986. Over Saint Santin, a plasma velocityof 25 m/s in the F regioncorresponds to an electricfield of 1 mV/m. We have used 15 years of Millstone Hill incoherent scatter radar observations from 1978 to 1992, comprisbances(time constantslongerthan a few hours)because ing more than 4500 hours of plasma drift observations, of the combined effects of ionosphericdisturbance dy- which have been recently reanalyzed. We have also namo electric fields and leakageof high-latitude electric used 14 years of plasma drift observations from the Saint Santin incoherent fields to lower latitudes. We have used extensive incoherent scatter radar ob- which were extracted scatter radar from 1973 to 1986 from the National Center for At- mospheric Research Coupling, Energetics and Dynam- servationsfrom Millstone Hill and Saint Santin to study midlatitude F region zonal and meridional plasma drifts during geomagnetic quiet and active conditions. We initially present a comprehensivesummary of the local time, seasonal, and solar cycle variations of the quiet icsof AtmosphericRegions(NCAR CEDAR) database. We discarded Saint Santin late afternoon and nighttime drift observations(1600-0700 LT) obtained be- time drifts over these two midlatitude expectedlylarge nighttime drifts (westwardvelocities stations. This in- tween September 1981 and March 1982 owing to un- quiet cludes,for the first time, the seasonaldependenceof the exceeding200 m/s, evenduring geomagnetically average solar maximum plasma drifts over Saint Santin. periods). The reasonfor theselarge drifts is unknown The main focus of this work is the study of the tempo- to us. The average values of the solar decimetric flux ral evolution of the midlatitude plasma drifts during indices for the Millstone Hill and Saint Santin data sets geomagneticactive periods using the methodologyin- are 140 and 125, respectively. The number and dates troducedby Fejer and $cherliess[1995].We haveused of the Millstone Hill and Saint Santin experiments used our quiet time drift model to determine the disturbance in this study can be obtained from the NCAR CEDAR drifts which were then modeled database. as a function of the AE index. In the following sections,we will first briefly review the radar measurement techniques and describe our data analysis. Then we present our quiet time model, describe the storm time dependent disturbance drift patterns and compare them with results from earlier studies. 2. Measurement Analysis Techniques and Data 2.1. Quiet Time Analysis We have used these data to determine empirical analytical modelsfor both the quiet time and disturbance plasma drifts over these midlatitude stations. For our quiet time plasma drift models, we have included only data with a current 3-hour magnetic index Kp<3.0 (Kp -•1.8). We will discusslater the implicationsof this relatively relaxed quiet time criterion. The electrodynamiccomponentsof the plasma drifts (perpendicular/eastward and perpendicular/northward The MillstoneHill radar (42.6øN,288.5øE,Apex magto B) were modeledas: netic latitude 54ø) measuresthe line-of-sightplasma drift velocity from the Doppler shift of the 440 MHz backscatter signal. The three-dimensionaldrift velocity over the radar site is obtained from measurements in 8 6 v(t,d,Sa)- Z • ai,kfkNi,4(t) i=1 k=l (1) SCHERLIESS ET AL.: MIDLATITUDE with PLASMA DRIFTS 1773 equatorward edge of the shielding layer is poleward of the radar fl = 1 May-August f2 = 1 fa = 1 November-February March-April, September-October f4 -- f• (Sa - 140) f5 = f2 ($a- 140) f6 = fs ($a- 140) Here t denoteslocal time, d is day of the year, and $• is the daily 10.7 cm solar flux index. Outside their above definedlimits, the functionsf• to f3 are set to zero. The site. We further minimized the inclusion of auroral drifts, and also of unusuallylarge drifts (associatedwith subauroralion drifts (SAIDs), for example) into our database, by limiting the magnitudes of the perturbation drifts to values smaller than 150 and 80 m/s for the zonal and perpendicular/northward drift component, respectively. Following the procedure describedby Scherliessand Fejer [1998], we have characterizedthe level of enhanced high latitude geomagneticactivity, using modified hourly averagedauroral electrojet (JEd) indices, with AEd = AE- 160 nT. For the time period prior to July 1988, we have used the standard AE index based coefficients ai,1 to ai,3 determinethe basicseasonal patternsof the electrodynamic driftsfor a solarflux of 140, on measurements from 12 auroral stations and, for the and the coefficients ai,4 to ai,6 representthe linearvari- later period, a provisional index based on 11 stations. ations of these drifts with the solar flux. The local time For the secondhalf of 1988 and also for 1989 (with dependenceis describedby eight normalized cubic-B the exception of March) AE indicesare not available, splinesof order four with equallyspacednodesat 0, 3, and, consequently, data from this period have not been 6, ..., 21 hours[e.g.,Fejer and$cherliess, 1997].These used in our perturbation drift analysis. The AE in- basicfunctionsare well suitedto describeplasmadrift variations up to terdiurnal oscillations. The model coef- dex has been empirically related to both the polar cap crosspotential and to the hemispheric high-latitude en- ficientshavebeendeterminedby weightedleastsquares ergy input [Jhn et al., 1983, 1992]. This allowsus to fits to the individual data sets, with each individual compare our empirical disturbance patterns with results observationweightedby the inversesquareof the mea- from global convection models. surementerror. However,in order not to overemphasize For each radar station the ionospheric perturbation measurementswith relatively small errorbars, which do drift components were expressedas not necessarilyrepresentthe averagequiet time pattern, 6 we have set errorbarssmallerthan one half of the dayto-day variability, obtained by the standard deviation v(t)- • [ai,iAAE(t - 30rain)+ i--1 of binned data presentedbelow, equal to this value. We have observed that variations of this criterion do ai,2AAE(t - 90 rain) + not changeour model resultsin any significantway and ai,3AEd(1 - 3 hour) + that evenan unweightedfit leadsto essentiallythe same (2) ai,aAEd(q- 9 hour) I Ni,•(t). drift patterns. Finally, a simple linear interpolation schemeover a The first two ternas under the summation resemble the range of 4- 15 days was employedfor the transition be- parametrizationusedby Scherliessand Fejev [1998]to tween our 4-month seasons.This interpolation provides accountfor prompt penetration drifts related to changes a reasonably realistic transition between seasons. 2.2. Storm Time Analysis As mentioned above, our quiet time model drifts correspond to an average Kp index of •-1.8. In the next step, we have determined the perturbation drift patterns during geomagneticactive periods and their storm time evolution. To extract the perturbation drifts out of our data set, we have subtracted the seasonand solar cycledependentaveragequiet time drifts, usingour new empirical quiet time models, from the measured in the auroral current systemswith averagetime delays of 30 and 90 min: respectively. To account for longer lasting drift perturbations, which have been parame- terized by $cherliessand Fejer [1998]by an averageof AE• indices over the past nine hours, we have instead used two separate parameters, i.e., the averageof AE• indices over the past 1 - 3 hours and the average over the past 4- 9 hours. This parametrization significantly improves the temporal evolution of the longer lasting perturbations and provides a better representationof the perturbation drifts during the initial hours of geodrifts (for a detailed descriptionseeFejer and Scherliess magneticactive periods. After severalhoursof auroral [1997]). The resultingperturbationsare due to storm activity, however,this parametrization givesessentially time generatedelectric fields and alsoto the day-to-day the same results as presented by $cherliess and Fejer [1998]. The limitations of the use of the AE• indices variability of the ionosphericdynamo electric fields. Sinceduring large nighttime storm periods, Millstone in our analysiswere discussedby Fejer and $cherliess Hill can be in the auroral zone, we have excludedplasma [1997,1998]. Our empiricalmodels,basedon the AE• drift observationswhen the maximum Kp index over the index, do not take into accountsomepotentially imporpreceding period of 6 hours was above 5.0. This crite- tant processesand averageout the effectsof a number rion restricts our Millstone Hill storm time analysisto of ionosphericand magnetosphericprocessesthat could moderately disturbed periods, i.e., to periods when the play important roles on the magnitude and phase of 1774 SCHERLIESS ET AL' MIDLATITUDE PLASMA DRIFTS MILLSTONE HILL MAY-AUG MAR-APR 60 Kp<_3 SEP-OCT NOV-FEB Sa=180 - .. Sa=190 40 - Sa=200 20 .• • 0 -20 _ " • , T T _ _ -40 __ 40 Sa=85 _ Sa=85 2O -20 _ 111] '1 '111 _ ., _ ,._. -40 I IIIIIIIIIIIII111 00 04 08 i IIIIII 12 16 20 24 00 04 08 12 LOCAL 16 20 24 00 04 08 12 16 24 20 TIME Figure 1. Quiet time averageperpendicular northwardplasmadrifts overMillstoneHill for low and highsolarflux conditions.The solidcurvesindicatethe velocitypatternsobtainedfrom our empirical model. Here $a denotesthe averagedecimetricsolarflux index. the electric field perturbations as well as on the shielding time constants. For example, results presentedby 3.1. Quiet Time Electrodynamic Drifts Fosteret al. [1986]suggest that IMF By mightcause Figure I showsthe average seasonalquiet time pat- Millstone Hill. A comprehensivestudy of these additional processesrequires a significantly larger database than currently available. In summary, we have developed individual empirical analytical models for the Millstone Hill and Saint Santin electrodynamicalplasma drifts for geomagnetically quiet and disturbed periods. These models include the seasonaland solar cycle variations of the quiet time drifts as well as short-term and longer lasting perturbation drifts during geomagneticallydisturbed periods. Combined, these models can describethe averagemidlatitude plasma drifts during various geomagneticcon- ponent over Millstone Hill for low and high solar flux plasmadrift comlarge changes in the perturbation electric fields over ternsfor the perpendicular/northward Results and Discussion In this section we will discussbriefly our solar minimum and solar maximum average quiet time drift patterns obtained models. Then from our Millstone we will describe The data have been divided into summer (May-August),equinox(March, April, September,October), and winter (November-February). The average decimetric solar flux indices are •85 for the low solar flux intervals and 180, 200, and 190, for summer, winter, and equinox high solar flux periods, respectively. The standard deviation of our averagedrifts is between •20 and 40 m/s. The averagenumberof samplesper bin is •20, resulting in standard errors of the mean between•5 and 10 m/s. Figure I indicatesthat the perpendicular/meridionaldrifts are northward in the early morning and prenoon sector and southward at the later times, with a return to the northward direction in the late afternoon during low solar flux summer and equinox conditions. The northward morning drifts ditions. 3. conditions. Hill and Saint Santin our disturbance drift patterns and compare them with other observationsand also with predictionsfrom global convectionmodels. reachmaximum valuesbetween•20 and 30 m/s, and the afternoon southward drifts have maximum values of •10 m/s. Although, large solar cyclevariationsare observedin the perpendicular/northwarddrift component (e.g., equinoxearly morning), no significantsystematic effects are evident. The results in Figure I are SCHERLIESS ET AL.' MIDLATITUDE PLASMA DRIFTS 1775 solarflux data. The averageKp index is again • 1.8 with ported by Buonsantoet al. [1993]and Buonsantoand slightlylarger valuesduring summerand equinoxhigh Witasse[1999],whodiscussed their main characteristics solar flux conditions. The Saint Santin drifts resemble and large day-to-day variability. The latter work also many of the featuresobservedover Millstone Hill, with in good agreement with the average drift patterns re- reported differences between the vernal and autumnal drift patterns, which were most pronouncedin the early morning to noon period. These effectsare not being accounted for in our current study. The thick solid lines in Figure i represent our quiet time model predictions, which have been calculated us- northward drifts in the prenoon sector and southward drifts in the afternoon. In general, the solar minimum resultsshownin Figure 2 exhibit strongsemidiurnaloscillations and closely resemblethe solar minimum pat- ing (1) for eachdata point and binnedand averagedthe terns presentedearlier by Blanc and Amayenc[1979] who used 3 years of drift observations.However,they reportedlarge downward/southward drifts during De- same way as the observations. As expected, our model results are in good agreement with the average drift spondingresultsindicateupward/northwarddrifts. patternswith typical fluctuationsof --•1-5m/s. Systematic discrepanciesbetween model and data are observed only during low solar flux winter afternoon conditions when the empirical model underestimates the equator- ward/downwarddrift velocities. These differencesare largely due to our use of a linear solar cycle variation cember solsticenighttime periods, whereasour correThe Saint Santin northward drifts showstrong semidiurnal patterns near solar minimum and diurnal type variations during solar maximum. Saint Santin solar maximum drifts have not been studied earlier, and therefore ours is the first model representation for these drifts. Figure 2 indicates that from -•0800 to 1400 LT, the equinoctial drifts do not changesignificantly Figure 2 showsthe averageperpendicular/northward with the phase of the solar cycle, whereasthe solsti- in our empirical model. quiet time plasma drifts over Saint Santin for solar minimum and maximum conditions and the corresponding empirical model results. The average solar flux values for summer, winter, and equinox are •85, 85, and 75 for the low solar flux data and 180,205, and 180 for the high tial data show opposite variations with an increasein the northward drifts in the summer data and a decrease in the winter. The relatively large southwarddrifts in the morningfrom -•0400 to 0800 LT and in the afternoon from -•1400 to 1800 LT drastically decreasewith Saint Santin - MAR-APR Kp<3 NOV-FEB MAY-AUG • 60-SEP-OCT 40: Sa=180 Sa=180 • _ Sa=205 _ 20 _ 0 T _ -. _ / .- __ -20 _ -40 40 _ , Sa=75 _ Sa:85 20 Sa=85 _ 7- -- -20 _ _ 40 - - _ IIit11111111111tll 00 04 08 12 16 20 24 00 04 08 12 LOCAL 16 20 24 00 04 TI ME Figure 2. Same as Figure 1, but for Saint Santin. 08 12 16 20 24 1776 SCHERLIESS ET AL.' MIDLATITUDE PLASMA MILLSTONE HILL 6O -- MAR-APR - SEP-OCT - ., ß 40 TM - Kp<3 MAY-AUG - - - - NOV- FEB - - - Sa=180 '"•Sa=190 - 20 TM DRIFTS Sa-200 - ß '- _ ,. •/•-i• -2O • -4O __ _ ._ ß i -- _ _ ,. .. .. .. _ __ _ . ß 40 TM Sa-85 - - Sa--85 " -- Sa=85 20 -l , _ ,, ., . -2O ..... ./ -4O .. , __ ., I II 00 I[ 04 111 08 II IIII 12 I II 16 II I111 20 I I I I I I I I I I I I I I I I I I"1 I.I I I 24 0 04 OR 12 LOCAL Figure 3. 1R 20 IIIII :1. 0 ( IIIIIIIIIIIIIIII 04 08 II 12 16 20 2, TI ME Quiet time averagezonal plasmadrifts over Millstone Hill for low and high solar flux conditions. increasing solar flux during equinox and summer and even reverse toward the northward direction in the June solsticemorning sector. It is not clear if the December solstice drift pattern is entirely realistic or biased by instrumental offsets. First of all, the December solsticedaytime and nighttime drifts are significantly different from the equinoctial values. In addition, as pointed out by Takami et flux values given for the northward component. The zonal drifts over Millstone Hill are westward at night with maximumvaluesof •60 m/s in winterand 40 m/s in summer. During daytime, and particularly around noon, the zonal drifts have eastward amplitudes with largestmagnitudes(up to 35 m/s) duringequinox.During summer, the average zonal drifts are always westward. The standard deviation is between •25 and 50 al. [1996], the nighttime Saint Santin winter drifts m/s, indicatinglarge quiet time variability. We will shouldbe strongly affectedby the correspondingconjugate summer ionosphere,which lies at a geographiclatitude of 31øS. Therefore we would expect a closesimilarity between the Saint Santin December solsticeand the show later that some of this variability was introduced by the choice of our quiet time binning criterion. The thick solid lines in Figure 3 represent again the predictions of our quiet time model. Generally, the zonal Shigaraki (35øN dip latitude) June solsticenighttime plasma drifts at Millstone Hill tend to be more westward drifts measuredby the MU radar [Takamiet al., 1996]. at night and more eastwardduring the day with increasHowever, the June solstice solar maximum nighttime ing solar flux, although, systematic variations are relaaverage drifts measured by this radar show southward tively small (of the orderof 10 m/s). Theseresultsare drift velocitiesof only •5 - 10 m/s, which would sug- also consistentwith thosefrom earlier studies[Buonsanto et al., 1993;Buonsantoand Witasse,1999]. gest a DC bias. Clearly, our solar maximum December solstice Saint Santin drifts should be considered with a Figure 4 showsthe average seasonalsolar minimum ahd maximum quiet time zonal drift patterns over Saint high degree of caution. The seasonallyaveraged Millstone Hill eastward drift Santin. The average solar flux indices are the same as pattern for geomagneticallyquiet conditionsare shown for Figure 2. The solar minimum drift patterns are in in Figure 3. The data are again shown for low and good agreementwith the resultsof Blanc and Amayenc high solar flux conditions with the same averagesolar [1979],with the exceptionof the summerearly night- SCHERLIESS ET AL.' MIDLATITUDE PLASMA Saint Santin - MAR-APR "- SEP-OCT 80 -l 1777 Kp_<3 MAY-AUG - - NOV_ FEB - ' _ 60 I -- 40 20 - DRIFTS _ Sa=180 .... :: ?',,,•,•. _ __ _ .... -20 -40 - ß .. -60 .. z .. o N 40 _ - Sa=75 _ " - _ Sa=85 _ - Sa=85 _ _ - _ - 20 _ _ _ _ .. -20 ß -40 -60 -I I I I I I I ] I I I I I I I I I I I I I I 1- lJlllllllllllllllltJlll 00 04 08 12 16 20 24 00 04 08 12 L¸CAL Figure 4. reversal 20 24 -I I I I I I I'1'1 I I I I I I I I I I I J I I00 04 08 12 16 20 24 TI ME Same as Figure 3, but for Saint Santin. time drifts for which their study indicated a westward drift 16 which is not seen in our results based on hemisphere,which is located at a significantlylarger geographiclatitude (45øS),largelydeterminesthe Arecibo a significantly larger data set. The Saint Santin zonal drifts vary significantlywith the phase of the solar cy- nighttime plasma drifts [Fcj½r, 1993; Ta]carai½t al., 1996]. The conjugatepoint of Saint Santin, however, cle. The nighttime drifts increasefrom •20-30 m/s to lies at a lowerlatitude (18øS)than this station. There•40-80 m/s betweenour low and high flux conditions. fore, as expected, the solar cycle variations of the Saint Figure 5 shows the combined winter and equinox average solar cycle variation in the early nighttime period. In this case, we have used solar flux bins of 30 units for indices between 60 and 210 and a single bin for larger values. The vertical bars indicate the standard deviation in each bin, and the solid line is our empirical model representation. These drifts increase linearly with increasingsolarflux at a rate of •40 m/s per 100 flux units. The solar cycle variations during summer are •50• kOSaintSantin September-Apri • 'G' 100 I---' Empirical Model . smaller. Zonal plasma drift observationsat Jicamarca (12øS aA, Arecibo(18øNgeographic, magneticlatitude30øN)[F½jet, 1993],and from the IDM onboardthe DE 2 satellite [Sch½rliess, 1997] indicate increasingsolar cyclevariations with decreasinglatitude, in good agreementwith the observed solar cycle dependence at Millstone Hill and Saint Santin. It is interestingto note that solarcycle variations at Arecibo are smallest during December solstice,when the ionosphericdynamo in the conjugate 50 %o"i'b'o' "i56'"i5&"'i&&" Solar 5i6'";>6 FIux Figure 5. Solar cycle variation of the equinox and winter early night Saint Santin zonal drifts. 1778 SCHERLIESS i 100- i i i i i i [ i [ i i i i i i ET AL.: MIDLATITUDE i i i [ i 1 PLASMA i -- _ - Millstone Hill 4020- - ..... - Kp(0-6hrs)=l.8 Kp(0-6hrs)=0.7 i , , , ito,t1[t,2 , i , , it3,t4it5, , i ,tl6 - •-500 - LU 3001ø 8060- DRIFTS 100•- _- -3 -- ' , 00 ' , I , , , ;, ' I" , 0'3 06 09 , , 1'2 " 15 Storm-Time (Hours) 0 Figure 7. -20- Idealized scenarioof the variation of the AE index used in this work. 40- _ _ 2øfTTT II!TTT quiet drifts (Kp = 0.7) correctedfor seasonaland solar cycle effects. Under the same conditions, the amplitudes of the Saint -40 .................. 02 06 10 Local 14 Time Santin zonal disturbance drifts about a factor of two smaller than over Millstone 18 are Hill. The effectsof the steady state leakageof high-latitude zonal electric fields, which drive perpendicular northward disturbancedrifts, are essentiallynegligiblefor Kp 22 Figure 6. (top) Yearly averageMillstone Hill perpendicular eastward drifts for two levels of magnetic activity. (bottom) Averagezonal drifts by subtracting the drifts in Figure 6; top. The solid curve indicates the pattern from our disturbance model. Santin zonal drifts are smallest in the summer, when they are determined mostly by the dynamo action in the local hemisphere. The quiet time model described above allowed us to determine the perturbation drifts which were used in the storm time study presented in section 3.2. As we will see, high-latitude electric fields affect the lowerlatitude plasma drifts with decreasingamplitudes toward the equator, even under quasi-stationary condi- < 3 even over Millstone 3.2. Midlatitude Hill. Disturbance Plasma Drift Patterns Figure 7 shows an idealized storm scenariowith an increaseof the AE index by 300 nT above our quiet time value of 130 nT during a period of 9 hours. Figure 8 presentsthe Millstone Hill disturbancedrift patterns for the storm times defined in Figure 7. The data and the scatterbarswereobtainedby binningand averaging the data for the conditionsin Figure 7 and smoothedby a three-point running average. The averagevaluesfor the parametersdefinedin (2) are givenin Table 1. The solid curvespresent the results from our empirical analytical modelobtainedfrom (2). It is important to note tions. Therefore,for a givenquiet time criterion(based that our data binning can only approximately reproon the Kp index, for example), we expect our quiet duce the idealized storm scenario and that it does not time drifts to have latitudinally increasingcontributions completely separate the responsesof the prompt penedue to electric fields of high-latitude origin. Figure 6 tration and disturbance dynamo processes.In spite of shows that the zonal disturbance drifts over Millstone these limitations, the binned data is generally in good Hill havesignificantamplitudesevenfor Kp _<3 (Kp = agreement with the results from our simultaneousmul1.8). In this caseFigure 6 (top) was obtainedby bin- tiparameter fitting procedureand, as we will seelater, ning the yearly averaged Millstone Hill zonal drifts for also with theoretical patterns obtained from the Rice geomagnetically quiet (average6-hourKp<_3.0)and ex- ConvectionModel (RCM). Figure 8 indicates that at time to, following an intremely quiet (average6-hour Kp_<l.3) conditions,and Figure 6 (bottom) showsthe differenceof the two curves creasein the AE index by 400 nT, the prompt penetraand the results obtained from our disturbance model. The scatterbars indicate the variability of the extremely tion zonal electric fields drive perpendicular northward drifts during the day and larger amplitude southward Table 1. Average Storm Time Parameters Storm Time to t• t2 t3 t4 t5 t6 AAE(t- 30 min) (nT) AAE(t- 90 min) (nT) AEd(1-3 hour) (nT) AEd(4-9 hour) (nT) 243 48 24 60 -33 209 149 79 -58 -27 293 89 -34 -29 287 275 -225 46 271 126 -41 -230 249 141 11 -7 30 165 SCHERLIESS ET AL' MIDLATITUDE Millstone PLASMA DRIFTS 1779 Hill 1978-92 _ 2O 10 0 -10 [ [ [ [ I [ ] [ ] I [ _ 3O z 15 0 -15 Z .,,•• o z _.. 10 0 -10 z - t1/••,.•, • 15 : ..... _ 0 -15 t 15 0 -15 10 0 -10 "--...__; t 10 0 -10 15 10 0 -10 15 0 -15 • ,. ß( t5 ...-15 0 -15 10 0 -10 ;:.. )( • , ,, 10 0 -10 15 0 -15 z _ ] 02 Figure 8. 06 10 14 18 Local Time 22 02 ] • 06 t I I [ ] [ 10 14 18 Local Time I 22 Averageperpendicular northwardandeastward disturbance driftsoverMillstone Hill for the conditionsand storm times shownin Figure 7. The solid lines indicate the patterns from our disturbancemodel,and the dashedlinesat stormtimes t2 and t3 wereobtainedusing the parametrization of $cherliess andFejer[1998].The scatterbarsdenotethe standarderror of the means. 1780 SCHERLIESS ET AL.' MIDLATITUDE PLASMA DRIFTS driftsat nightwith a maximumvalueof •-25 m/s at - ' - 0400 LT. The initial time responseof the meridional electric fields generatesa small eastward disturbance drift in the earlymorning-noon periodandlargerwestwarddrifts at otherlocaltimes. The largestinitial time westwarddrift perturbationis •40 m/s, and it occurs near dusk. As stormtime increases,the disturbance drifts initially decreasein amplitude and shift to later local times but do not changemuch after •-2 hoursof continuousmagneticactivity, as shownby the patterns at storm times t2 and t3. Stormtime t4 (Figure8) showsthe disturbance patternsresultingfrom the suddendecrease in the AE by 400 nT. In this case,the nighttime perturbationdrifts ' ' ' i , , , , i , , , , i , , , , i i , , _ Radar Model AAE=400nT 80-- ..... DE-2Model - _ _ 60-- _ Millstone Hill - 40-- _ _ 20: ..-----•ial Time Response -_ -20 / •. _ •) 20 - E 0 -- t0+90min •., _ • -20 __ - _ _ t'- o U - -- Saint Santin - - _ _ changefromthe downward/southward directionto up20Initial Time Response ward/northward, andtheprenoon downward/southward . .._...• ;::-,_.-._ driftsarefurtherincreased. Themostnoticeable changes in the zonaldriftsoccurbeforenoonwhenthey become 20to+90min westward. These perturbation drifts are again due to the combinedeffectsof largetransientpromptpenetration drifts (with oppositesignto thosefor storm time to in Figure 8) and longerlastingtime-delayeddistur'0'2' ' '0'6' ' '1•0 14• ' ' 1'8'• '2'2' bances.The agreementbetweenthe averageddata and Local Time our model predictions(solidline) is very good,but a phaseshiftof roughly1-2 hoursin the promptpenetra- Figure 9. Comparisonof the prompt penetration tion responsetowardlater local timesfor both compo- zonal disturbancedrift patterns derived from the Mill- _ _ _ - _ _ _ _ nents would make the agreementevenbetter. It is inter- estingto notethat Blanc[1983]hasnotedthat a similar phaseshiftwould also bring the Saint Santin drift data in closeragreementwith his model results. Onehour after the decrease of the AE index,at storm stone Hill and Saint Santin radar measurements and from DE 2 observations after an increase in the AE in- dex by 400 nT. The thin and thick dashed lines in the Saint Santin part indicate the DE 2 resultsobtained for A=35ø and A=45ø, respectively. time t5, the meridionalperturbationdrifts have significantlysmalleramplitudes. The zonaldrifts, how- It is important to note that the DE 2 results were deever, still exhibit large nighttime westwardperturba- rivedfrom longitudinallyaveragedmeasurements. Figtions, and small eastward perturbations around sunrise. ure 9 (bottom) comparesthe Saint Santin resultsand Storm time t6 (Figure 8) showsthat after 6 hoursof the correspondingDE 2 results. Fejer and $cherliess quieting, the perpendicular/northwarddrift perturba- [1998]pointedout that in the postmidnight-morning tions have largely returned to their quiet time level, sector,the DE 2 initial time disturbancedrifts change but the zonaldrifts still displaylargewestwardpertur- noticeablyfor invariantlatitudesfrom45ø to 35ø. Figbationsin the eveningsector,whichare not fully cap- ure 9 showsthat surprisinglythe Saint Santin initial tured by our current model. This suggeststhat other time disturbancepattern is in better agreementwith disturbancemechanismsor drift perturbationswith dif- the DE 2 resultsfor an averageinvariantlatitude of 35ø ferenttimescales than considered in ourmodelmightbe thanfor A = 45ø. FejerandScherliess [1998]suggested operating during the recovery phase. that the large changeof the DE 2 initial zonal distur- The smaller database of drift measurements from Saint Santin did not allow us to determine their storm- bancedrift patternsfrom 45ø to 35ø was due in part to the largeruncertainties resultingfromthe largede- time dependence in asmuchdetail asgivenin Figure8. creasein the numberof postmidnightmeasurements at However, as will be shown below, we were still able to lowerlatitudes.The SaintSantinresultssuggest that capture the basic featuresof their prompt penetration this argument is questionable. The t0+90 rain DE 2 and longer lasting disturbancedrift patterns. andSaintSantindisturbance drift patternsarein good $cherliessand Fejer [1998]and Fejer and$cherliess agreement,althoughthe radar data showlarger west[1998]have used extensiveion drift meter data from ward drifts in the duskto early night period. the DE 2 satellite to determine the middle- and low- latitude ionosphericzonal disturbancedrift patterns. Figure 9 (top) showsthe good agreementbetweenthe Millstone Hill prompt penetration zonal drifts at storm times t0+30 rain and t0+90 rain and the corresponding DE 2 drifts for an averageinvariant latitude of 55ø. Figure 10 showsthe Millstone Hill and Saint Santin a•d DE 2 disturbancepatterns after an increasein the AE index by 400 nT over a period of 9 hours. The radar drifts are again consistentwith the satellitedrifts but havesystematically largermagnitudes, in particular duringthe duskto early nighttimehours.However,we SCHERLIESS ET AL.- MIDLATITUDE [ [ ! ! ] ! [ , [ i , , ! ! ] 60- 40- - ! , Radar AEd(1-9hrs)-400 nT 20- ! Millstone , i ! PLASMA DRIFTS 1781 Figure11 showsthe comparison of the MillstoneHill and SaintSantinpromptpenetrationdrift patternsand , Models --- DE-2Models the RCM results for an increase in the polar cap poten- Hill _ tial of 33 kV. These RCM patterns, which correspond 0 to run C of Spiroet al. [1988],werediscussed in detail by Fejer et al. [1990]and comparedto DE 2 model resultsby Fejer and Scherliess[1998]. Figure 11 in- -20 -40 dicates that the initial time responseof the Millstone Saint Santi 2O 0 _ -20 -40 .... 02 ]06.... 110 .... Local 14 [ ....18 t ,2,2 , Time Hill driftsare in goodagreement with the RCM results. Notice that the drift amplitudesdecreaseand shift to later local times as the plasma sheet inner edge adjuststo the new crosspolar cappotential.The Saint Santinperpendicular meridionaldisturbance driftsare in generalagreementwith the modelresults,whereas the zonal drifts near dusk are noticeably larger. Fe- Figure 10. Sameas Figure 9, but for the long-lasting jet and $cherliess[1998]havepointedout that these RCM quasisteadystate drifts underestimate the DE 2 drifts. zonal disturbance zonal disturbance drifts for latitudes A _< 45 ø consistent with our results. The results above indicate that have to reiterate here that the DE 2 results represent the radar and DE 2 perturbation drift patterns are in longitudinallyaveragedpatternsand haveseasonand reasonablygoodagreement.They are alsogenerally local time lockedtogether (the early nighttime period, consistentwith the predictionsfrom the RCM, as well as from other global convectionmodels. for example,corresponds to equinoctialconditions). 4. Comparison With Theoretical Models 5. Summary and Conclusions We have used incoherent scatter radar observations from Millstone Hill and Saint Santin to study the char- The penetrationof high-latitudeelectricfieldsto middle, low, and equatoriallatitudeshas alsobeen investi- acteristicsof the averagemidlatitude quiet time and gatedusingseveralglobalconvection models[e.g.,Spiro disturbanceelectrodynamicplasma drifts. Our quiet e• al., 1981, 1988; Senior and Blanc, 1984; Zakharove• time Millstone Hill drift patterns are in good agreeal., 1989; Fejer et al., 1990; Peymirat and Fontaine, ment with earlier studies based on smaller databases. wereusedto 1994;Peymirat,1998; Tsunomura,1999]. Thesemod- The Saint Santinquiet time measurements variationsof els calculatethe global ionosphericelectricfields and determine,for the first time, the seasonal currentsby solvingthe continuityequationfor iono- these drifts for both low and high solar flux conditions. sphericcurrentsfor a givenhigh-latitudeelectrostatic The Saint Santin eveningzonal drifts increasewith sothe corresponding MillstoneHill drifts potential or the field-alignedcurrent distributionon a lar flux, whereas two-dimensional (thin shell)ionosphere with giveniono- are essentiallysolar flux independent. We have modeled the Millstone Hill and Saint Santin sphericconductances.The middle- and low-latitude perturbationelectricfieldsobtainedfromdifferentmod- disturbancedrifts obtainedby removingfrom eachmeaels are similar sincethey depend mostly on the high- surement the season and solar flux dependent quiet latitude potential penetratingto lower latitudes and time values. These stormtime dependent perturbation on the ionosphericconductances and not on the details drifts were usedto study and model,for the first time, of magnetospheric processes. The mostcomprehensivethe temporalevolutionof the promptpenetrationand theoreticalprompt penetrationresultshave been pro- longerlastingzonaland meridionalE xB disturbance videdby the RCM whichaccounts for the coupledelec- drifts. The radar disturbance patterns are in good trodynamicsof the innermagnetosphere andionosphere agreementwith resultsderivedfrom zonaldrift mea[e.g.,Wolfet al., 1986;Spiroet al., 1988]. The local surementsby the DE 2 satellite. This, and alsoearlier time and latitude zonal and meridional perturbation studies,alsoindicatesthat the RCM, as well as other electricfields predictedby the RCM were studiedby detailed convectionmodels, can reproduce the average Fejer et al. [1990].We cancomparetheir resultswith disturbancedrift patterns with reasonableaccuracy. The present empirical model providesa relatively the radar disturbance drift patterns presented above by simple description of the time dependentevolutionof usingthe empirical relationshipbetweenthe polar cap potential drop and the AE index derivedby Ahnet al. [1992],i.e., •(kV): the midlatitude disturbance drifts based on the history 36 + 0.082AE•2(nT),whereAE•2 of the AE index. A copy of the current model can be is the auroral electrojet index using 12 stations. In this obtained from the authors. The major challengefor the of the large case,a changein the AE index by 400 nT corresponds next few yearslies in the understanding to a changein the crosspolar cap potentialof •33 kV. variability of these disturbancepatterns. 1782 SCHERLIESS ET AL.' Radar ] [ ] [ I i ] [ MIDLATITUDE Model [ PLASMA DRIFTS ....... RCM I [ [ ' ] I ] [ [ [ I [ _- AAE=400nT _- AAE=400nT _- - _- _- Millstone Hill - MillstoneHill Init al Time Response 40 - 20 •- ,,' ,,' -,... InitialTimeResponse 60 : \ 30 0 -20 '" '• ' " -...... -30 to+60min- to+60m'n "",,..,,'"• 30 0 -30 - Saint Santin _- Saint Santin _ _ •. •_. 10-- - Initial Time Response;, F :/'-"',, oL -10 '. ,,'/ ', ; '. - - - ,' "-•', ',.',_ '.' "'--:•'.--/' "•• , ,: Initial Time Response - 15 0 ',• _•' / '"'" '-'-15 to+60min • ", '- 0 ..,•_ _.,"•-; "/'••- - -10 02 06 10 Local 14 Time 18 22 02 06 10 Local 14 Time 18 15 0 -15 22 Figure 11. Comparison of the disturbance drift patternsobtainedfromour MillstoneHill and SaintSantinmodels withresults fromtheRiceConvection Model(RCM)foran increase in the polarcappotentialdropby 33 kV. Herethe RCM initial time responses correspond to a time (to+ e)immediately afterthechange in • (short-dashed line)and10minlater(long-dashed line). SCHERLIESS Acknowledgments. ET AL.: MIDLATITUDE This work was supported by the Aeronomy Program, Division of Atmospheric Sciences,of the National ScienceFoundation through grants ATM-97146 77 and ATM-9731704 and the National Aeronautics and SpaceAdministration through grant NAG5-4469. The Saint Santin data were obtained through the CEDAR database, which is supported by the National ScienceFoundation. We also thank WDC-C2 Kyoto AE index service for providing us with the preliminary AE indices. Janet G. Luhmann thanks Oliver G. Witasse and W. Robin Coley for their assistencein evaluating this paper. References Ahn, B.-H., S.-I. Akasofu,and Y. Kamide, The Jouleheat productionrate and the particle energyinjectionrate as a functionof the geomagnetic indicesAE and AL, J. Geophys. Res., 88, 6275-6287, 1983. Ahn, B.-H., Y. Kamide,H. W. Kroehl,and D. J. Gorney, Crosspolarpotentialdifference, auroralelectrojetindices, and solar wind parameters,J. Geophys.Res., 97, 13451352, 1992. Bauer, P., P. Waldteufel, and C. Vialle, The French multistatic incoherentscatter facility, Radio Sci., 9, 77-83, 1974. Blanc, M., Magnetosphericconvectionelectric fields at midlatitudes, 3, Theoretical derivation of the disturbance con- vectionpattern in the plasmasphere,J. Geophys.Res., 33, 235-251, 1983. Blanc, M., and P. Amayenc, Seasonalvariations of the iono- sphericExB drifts aboveSaint Santin on quiet days, J. Geophys.Res., 8•, 2691-2704, 1979. Blanc, M., and A.D. Richmond, The ionosphericdisturbancedynamo, J. Geophys.Res., 85, 1669-1688,1980. Buonsanto, M. J., and O. G. Witasse, An updated clima- PLASMA DRIFTS 1783 latitudes observedby Dynamics Explorer 2, J. Geophys. Res., 97, 19,461-19,469, 1992. Jaggi, R. K., and R. A. Wolf, Self-consistentcalculation of the motion of a sheet of ions in the magnetosphere,J. Geophys. Res., 78, 2852-2866, 1973. Mazaudier, C. A., and S. V. Venkateswaran,Delayed ionospheric effects of the geomagnetic storms of March 22, 1979 studied by the sixth co-ordinated data analysiswork- shop (CDAW-6), Ann. Geophys.,8, 511-518, 1990. Mazaudier, C. A., A.D. Richmond, and D. Brinkman, On thermosphericwinds producedby auroral heating during magnetic storms and associated dynamo electric fields, Ann. Geophys.,5, 443-448, 1987. Peymirat, C., A magnetosphere-thermosphere-ionosphere electrodynamics general circulation model, J. Geophys. Res., 103, 17,467-17,477, 1998. Peymirat, C., and D. Fontaine, Numerical simulation of magnetospheric convection including the effects of field aligned currents and electron precipitation, J. Geophys. Res., 99, 11,155-11,176, 1994. Richmond, A.D., Ionosphericelectrodynamics,in Handbook of Atmospheric Electrodynamics,vol. 2, edited by H. Volland, pp. 249-290, CRC Press, 1995. Scherliess,L., Empirical studiesof ionosphericelectric fields, Ph.D. dissertation, Utah State Univ., Logan, 1997. Scherliess,L., and B. G. Fejer, Storm time dependenceof equatorial dynamo zonal electric fields, J. Geophys. Res., 102, 24,037-24,046, 1997. Scherliess, L., and B. G. Fejer, Satellite studies of mid- and low-latitude ionospheric disturbance zonal plasma drifts, Geophys.Res. Lett., 25, 1503-1506, 1998. Senior, C., and M. Blanc, On the control of magnetospheric convection by the spatial distribution of ionosphericcon- ductivities, J. Geophys. Res., 89, 261-284, 1984. Spiro, R. W., M. Harel, R. A. Wolf, and P. H. Reiff, Quantitative simulation of a magnetospheric substorm, 3, Plastology of thermosphericneutral wind and F regionion maspheric electric fields and evolution of the plasmapause, drifts oboveMillstoneHill, J. Geophys.Res., 10•, 24,675J. Geophys. Res., 86, 2261-2272, 1981. 24,687, 1999. Spiro, R. W., R. A. Wolf, and B. G. Fejer, Penetration of Buonsanto, M. J., M. E. Hagan, J. E. Salah, and B. G. high-latitude electric field effects to low latitudes during Fejer, Solarcyclevariationsin Fregion electrodynamics at SUNDIAL 1984, Ann. Geophys.,6, 39-50, 1988. MillstoneHill, Y. Geophys.Res., 98, 15,677-15,683,1993. Takami, T., W. L. Oliver, A.D. Richmond, and S. Fukao, Fejer, B. G., F region plasma drifts over Arecibo: Solar Ionospheric drift similarities at magnetic conjugate and cycle,seasonal,and magneticactivity effects,J. Geophys. nonconjugate locations, J. Geophys. Res., 101, 15,773Res., 98, 13,645-13,652, 1993. Fejer, B. G., The electrodynamicsof the low-latitude iono- 15,782, 1996. Tsunomura, S., Numerical analysisof global ionosphericcursphere: Recent resultsand future challenges,J. Atmos. rent system including the effect of equatorial enhanceSolar Terr. Phys., 59, 1465-1482, 1997. ment, Ann. Geophys., 17, 692-706, 1999. Fejer, B. G., and L. Scherliess,Time dependentresponse Wand, R. W., A model representation of the ionospheric of equatorialionospheric electricfieldsto magnetospheric disturbances,Geophys.Res. Lett., 22, 851-854, 1995. Fejer, B. G., and L. Scherliess,Empirical models of storm time equatorialzonalelectricfields,J. Geophys.Res., 102, 24,047-24,056, 1997. Fejer, B. G., and L. Scherliess, Mid- and low-latitudeprompt. penetrationionosphericzonalplasmadrifts, Geophys.Res. Left., 25, 3071-3074, 1998. Fejer, B. G., R. W. Spiro, R. A. Wolf, and J. C. Foster, Latitudinal variationof perturbationelectricfieldsduring magnetically disturbaed periods: 1986 SUNDIAL obser- vations and model results, Ann. Geophys.,8, 441-454, 1990. Fejer, B. G., E. R. dePaula, S. A. Gonzales,and R. F. Woodman, Average vertical and zonal F region plasma drifts electric field over Millstone Hill (A = 56ø), J. Geophys. Res., 86, 5801-5808, 1981. Wolf, R. A., R. W. Spiro, and G. A. Mantjoukis, Theoretical comments on the nature of the plasmapause, Adv. Space Res., 6, 177-186, 1986. Zakharov, V. E., M. A. Nikitin, and O. A. Stairnov, Re- sponseof low-latitude electric field to the action of magnetic source, Geomagn. Aeron., 29, 344-349, 1989. C. Amory-Mazaudier, Centre d'Etude des Environnements Terrestre et Planetaires, 94107 Saint-Maur-deFosses, France. M. J. Buonsanto,L. Goncharenko,and J. Holt, Haystack Observatory,Massachusetts Institute of Technology,West- overJicamarca, J. Geophys. Res.,96,13,901-13,906, 1991. ford, MA 01886. Foster,J. C., J. M. Holt, R. G. Musgrove, and D. S. B. G. Fejer and L. Scherliess,Center for Atmosphericand Evans,Solarwinddependencies of high-latitude convec- Space Sciences,Utah State University, Logan, UT 84322tionandprecipitation, in SolarWindMagnetosphere Cou- 4405. ([email protected], [email protected]) pling,editedbyY. KamideandJ. A. Slavin,pp. 477-494, Terra Sci., Tokyo, 1986. Heelis, R. A., and W. R. Coley, East-west ion drifts at mid- (ReceivedJune 9, 2000; revisedAugust 11, 2000; acceptedAugust 14, 2000.)
© Copyright 2025 ExpyDoc