Lecture 17 – May 24, Ring Current Dynamics讲座1–5月24日环电流的动力学

1、ess 261 spring 2010advanced topicsr.l. mcpherronlecture 17 may 24, 2010ring current dynamicsorigin of concept of “magnetic storm”(d. stern - /education/wmagstrm.html) the term magnetic storm, meaning a world-wide magnetic disturbance, was coined by alexander von humboldt

2、(1769-1859). a naturalist who gained attention by exploring the jungles of venezuela, humboldt devoted much of his life to the promotion of science. he produced five volumes of kosmos (starting the modern usage of that term), an encyclopedic account covering the broad spectrum of the sciences. it wa

3、s kosmos which brought to the worlds attention the discovery of the sunspot cycle by heinrich schwabe. after journeying through siberia, humboldt convinced the czar to set up a network of magnetic observatories across the russian lands, and additional stations were established throughout the british

4、 empire, from toronto to tasmania. this network clearly showed that magnetic storms were essentially identical all over the world: a steep decrease of the field over 6-24 hours, followed by a gradual recovery which lasted 1-4 days. the change in the magnetic field was small, in modern units some 50-

5、300 nt (nanotesla) out of a total intensity of 30-60,000 nt, but its world-wide scale suggested that something quite big was happening out in space.interplanetary magnetic field, ae and dst indices during stormcoronal mass ejection produce intervals of strong southward bz at the earthmagnetic reconn

6、ection drives magnetospheric convectionconvection drives currents along field lines and through ionosphereground magnetometers record effects of ionospheric currents in h and other componentsh traces are used to construct the ae and dst index-1000100coupling for storm with min at 02:54 ut oct 11, 19

7、97 propagated wind for doy = 284sym-hasy-h-1000-5000500au/al (nt)0246pc index -4-20ey gsm(mv/m)0510pdyn (np) 00.020.040.060.08ptot (np) 1997/10/10 17.0 0.0 0.0 1997/10/11 17.0 0.0 0.0 1997/10/12 17.0 0.0 0.0 1997/10/13 17.0 0.0 0.0 1997/10/14 17.0 0.0 0.0 contributions to the variation in the h comp

8、onent00 ( ) ( )( )( )( )( )( )( ) ( ) secular variation of main field ( ) solar dynamo in ionosphere sqmpsrcprcrtscwsqh thththththththththtobserved magnetic field ( ) magnetopause (chapman-ferraro) ( ) symmetric ring current ( ) partial ring current ( ) tail current ( ) substorm current wedge mpsrcp

9、rctscwhththththt0 ( ) ( )( )( ) sqh th ththtdisturbance daily variationlongitudinal profile of bj from magnetospheric currentssymmetric ring should create nearly constant longitudinal profile in h componentlocal time average of h at equator approximates b at center of earthbut other magnetospheric c

10、urrents create local time dependent deviations from symmetryassume asymmetric component has zero mean when averaged over local timedefine the disturbance storm time index dst as local time average of observed h profile jnjjjstbhnd1cos1local time0012180612dstb ooooeffects of magnetopause on the dst i

11、ndexbalance magnetic pressure against dynamic pressure0510150-10-8-6-4-2246810x (re)z (re)solarwindneutralpointk vabbkavb ntpnpdyn 2202022235()()-6-4-20246-35-30-25-20-15-10-50-xgsm (re)bz (nt)normal tailinner edgetotaleartha sheet current model of effect of tail current on dsttail current modelmagn

12、etic effectsbzxxxxxxx xxx xxioszrrbblnriro30 ntr12 rer42 renninob 1000 ntr7 rer12 ressisob magnetic effects of a substorm current wedgetransverse currents in the magnetosphere are diverted along field lines to the ionosphereviewed from above north pole the projection of the current system has a wedg

13、e shapemidlatitude stations are primarily affected by field-aligned currents and the equatorial closure (an equivalent eastward current)the local time profile of h component is symmetric with respect to the central meridian of wedgethe d component is asymmetric with respect to center of wedgesolar w

14、ind pressure correctiondst is contaminated by the magnetopause currents and incomplete removal of secular variationdynamic pressure brings these currents closer to the earth and makes them strongerthe correction is usually presented as:other currents contribute to the dst index:partial ring currenti

15、onospheric currentstail currentfield-aligned currentinduced earth currentsit is not known how well dst represents the current produced by particles on closed drift paths around the earthdstdstb pcsw*dessler-parker-sckopke derivationdrift velocity of an equatorial ion in dipole fieldthe equatorial ri

16、ng current of the ionthe magnetic effect of the equatorial ring where is the ion energy and is its charge and is the distance from the earths dipole ve rqmeqrmiqvrbiremdd 322342100the magnetic effect of ions gyrationthe total magnetic effect at earths centerexpress as fraction of surface field where

17、 is the total field energy outside earthbrembbbembbeuus203012044223 units mksin )(105 . 2 )(32)( note functioninjection current ring thecalled is )( e wher)()()(decay timecurrent ring theis and current ring theinput toenergy of rate theis )( e wher)()()(earth outside field dipolein energy total thei

18、s particlescurrent ring theofenergy total theis )( field surface equatorial average theis current ring theofeffect theis )( re whe3/ )(2/ )(140*0*0*tujnttuebtqtqtdtqdttdtutetudttdeetebtdetebtdmststmstmstthe dessler-parker-sckopke relation and the burton equationburton, r. k., r. l. mcpherron, and c.

19、 t. russell (1975), an empirical relationship between interplanetary conditions and dst, j. geophys. res., 80(31), 4204-4214.the burton equationburton, r., r. mcpherron, and c. russell (1975), an empirical relationship between interplanetary conditions and dst, journal of geophysical research, 80(31

20、), 4204-4214. if we assume the energy in the ring current is governed by injection and decay, the dynamic equation is: which becomes the burton equation: q is the injection term, is the decay timede tdtu te t( )( )( )ddstdtq tdstt*( )( )a nonlinear equation for the rate of change of dstthe burton et

21、 al. equation for ring current dynamics is derived from the dps relationassume that the ring current is driven by a solar wind coupling function depending on imf bzcorrect the measured dst for effects of magnetopause and incorrect base linesobtain a first-order differential equation for dst with two

22、 solar wind drivers that is nonlinear in the unknown parametersthe equation is non-linear in the model parametersctcpbdpdtdbdtddcpbdddtctctqdtqdtdddynstdynstdynststststst)(11rearrange and substitute with observedconvert input solar wind is )( where )()(let )(the model parameters include:b dynamic pr

23、essure constant exponential decay time constant of proportion to esc dst base line adjustmentparameterprobablevalueprobablevalueprobablevalueunitsburtonoptimumfor bsoptimumfor bnb1 (1/)7.710.7316.95hoursb2 (b)15.96.988.87nt/npab3 (c)207.195.8ntb4 ()5.43.23na(nt/h)/(mv/m)pef15.737.829.1%parameters in

24、 the burton equationassume the parameters in the burton equation are constantsuse the burton et al. values and calculate the prediction efficiencyuse non-linear inversion to obtain self consistent values for southward imfrepeat for northward imfthe burton parameters are far from optimum values!the s

25、elf consistent decay time is different for imf north or south, and are longer than burtonthe self consistent pressure constant is only half as large as burtonthe self consistent injection constant is lower than burton because the decay time is longer vbsppdtddstdtddstdyndyn4321211bbbbbbthe obrien-mc

26、pherron equation the b parameter represents the influence of magnetopause currents on dst the parameter modulates the reconnection driven convection electric field the v0 parameter scales the total convection electric field)()(*vbsdstvbsqdtddstppbdstdst*cccevbsevbsevbsq0vbsvvqe0 ceqpbddtpdbdtddydyns

27、tdynst11the obrien & mcpherron equationassume that all parameters depend on vbzbin the data into narrow bins of vbzthe injection term is a constant for each bin so combine with baseline termuse data to determine optimum b, , q(ey)-c/ in each binplot the dependence of parameters on vbzthe hourly

28、dst index predictor (1)1()stststb tc tdtdbppta vb 2 10exp /36.902.87exp 4.475 4.590.53 0.53() 0 0.53sssssvbbntmv mhrsvbcntntvbvbhra vbmvvbsvbsm ring current response to dynamic pressure use model for hourly changes in dst index at constant vbs for equinox and solstice fit an exponential function to

29、the pressure parameter as function of in interval -10vbs 3 mv/m the model is poorly determined and could be taken as an exponential decrease with a decay parameter eo = 9.27 mv/m and bo = 7.84 nt*(np)1/2 the pressure parameter decreases in strength with stronger vbsmv/m -10vbz3 formv/m npnt/ .28 , 1

30、1.6 , 7.7,7.5,ebvbs.39 97.87.8bbring current decay timeuse model for hourly changes in dst index at constant vbs for equinox and solsticefit a complex exponential function to the decay parameter as function of vbzin interval -10vbs 0)/( , )/( ,)/(nt/hr) 223123mmvmmvmmvevbzvbzqvbz0.752.556.58bbbbbbbi

31、njection function and coupling constant the injection rate is determined from a model parameter fit the rate with a function having zero value and slope at vbz=0 determine slope of data and of fit to dataimproved estimates of dst model parametersmake taylor series about bin center for parameter qfin

32、d average value of this over bin of width and assign to center of massuse the probability distribution n of samples in the bin to determine center of mass of x in the binrepeat for five bins of different width and obtain an over determined set of equations for q(x0) and q(x0) as function of distance

33、 from algebraic centerfit a line to these five paired valuesthe intercept of this line is the value of the parameter q(x0) vbz0algebraiccentercenterofmass 00000000000000001but centerrbin is and let xxdxxnxqxqxqxdxnxdxnxxxqxdxnxqxqxdxnxxxqxdxnxqxqxxxqxqxqxvbxbinbinbinbinbinbinzcomparison of models fo

34、r all subsets ceqpbddtpdbdtddydynstdynst11conclusions from empirical modelthe dst pressure constant is not the value used by most numerical modelsthe dst pressure constant decreases to zero as vbz becomes more negativethe decay time is shorter for large negative vbzthe injection rate is not a linear

35、 function of vbz using better estimates of q at more negative valuesthe coupling parameter increases continuously towards negative vbzb 8 nt/(pdyn)1/2b goes to 0 near -10 mv/m3 tau 16 hrthe injection rate q goes to zero nonlinearly as vbz 0there is no detectable seasonal difference in couplingthere

36、is no detectable difference in coupling as function of alfven mach number or plasma betasmall & big storms050100150-120-100-80-60-40-20020dst comparison for storm 1980-285dst (nt)0501001500123456ec = 0.49 mv/m vbs mv/mepoch hoursdst model (1hr step) model (multi-step)vbs 020406080100120140160180

37、-250-200-150-100-50050dst comparison for storm 1982-061dst (nt)020406080100120140160180051015vbs mv/mepoch hoursdst model (1hr step) model (multi-step)vbs ec = 0.49 mv/m dst parameters for high and low alfven mach #low mach # solar wind has a weak b field so vbz does not reach large negative valuest

38、here are many high mach # values at strong negative vbzthe parameters are not well determined for low mach #in the overlapping range there is no significant difference!dst parameters for high and low betahigh beta solar wind has a weak b field so vbz does not reach large negative valuesthere are man

39、y low beta values at strong negative vbzthe parameters are not well determined for high betain the overlapping range there is no significant difference!simulated dst indexebihara, y., and m. ejiri (2000), simulation study on fundamental properties of the storm-time ring current, j. geophys. res., 10

40、5(a7), 15843-15859. a storm in 1997 was simulated using polar cap potential and plasma sheet density as drivers the thick line in bottom panel using actual density slightly overestimates dst a simulation with a constant smaller density seriously underestimates dst both convection field and plasma sh

41、eet density are important driversthe ring current coupling constantthe coupling constant is defined as the slope of the injection functioncalculate this numerically from injection function by a running linear fit to seven points (blue)compare this to the derivative of a fit to the injection function with constraints listed abovecoupling decreases to zero at vbz = +0.75 mv/mfor large negative vbz it approaches 6.6 (nt/hr)/(mv/m)2311bbbvbzevbzddqdependence of coupling on dipole tilt separate data i