首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We have studied the stability of the electrostatic electron cyclotron wave in a plasma composed of hydrogen, oxygen and electrons. To conform to satellite observations in the low latitude boundary layer we model both the ionic components as drifting perpendicular to the magnetic field. Expressions for the frequency and the growth rate of the wave have been derived. We find that the plasma can support electron cyclotron waves with a frequency slightly greater than the electron cyclotron frequency ω ce ; these waves can be driven unstable when the drift velocities of both the ions are greater than the phase velocity of the wave. We thus introduce another source of instability for these waves namely multiple ion beams drifting perpendicular to the magnetic field.  相似文献   

2.
The plasma wave instrument (PWI) on board the Polar spacecraft made numerous passages of the dayside magnetopause and several probable encounters with the magnetosheath during the years 1996 and 1997. During periods of relatively high density, the PWI antenna-receiver system is coupled to the plasma and oscillates. The oscillations have been shown (cf. Radio Sci. 36 (2001) 203) to be indicative of periods of higher plasma density and plasma flows, possibly associated with magnetic reconnection. We have studied the plasma waves observed on three distinct magnetopause passes distinguished by the presence of these oscillations of the PWI receivers, and we report on the data obtained near, but not during, the times of the oscillations and the possible role of these waves in magnetic reconnection. Sweep-frequency receiver and high-resolution waveform data for some of these times are presented. The plasma wave measurements on each of the passes are characterized by turbulence. The most stable waves are whistler mode emissions typically of several hundred hertz that are seen intermittently in these regions. The data indicate the presence of impulsive solitary-like wave structures with strong electric fields both parallel and perpendicular to the magnetic field near, but not always within, suspected reconnection sites. The solitary waves show the highest occurrence when observed with electrostatic electron cyclotron waves. These latter waves have been observed in the past in the cusp, polar magnetosphere, and auroral regions and therefore may represent excursions into the cusp, but also indicate the presence of low-energy electron beams. Turbulence near the lower hybrid frequency, low-frequency EM waves, and impulsive monopolar electrostatic pulses are seen throughout the magnetopause and particularly near regions of large decrease in the local magnetic field and enhanced field-aligned flows, the suspected reconnection sites. The absence of significant solitary wave signatures within suspected reconnection sites may require modifications to some reconnection models.  相似文献   

3.
A theoretical study is made of the whistler mode cyclotron instability both in linear and nonlinear regimes in conjunction with the generation of VLF emissions in the magnetosphere. For the nonlinear treatment, a well-established quasilinear method is used and some physical processes of the cyclotron instability viz. energy conservation, mechanism of instability and frequency change of the excited emissions are clarified. The results are applied to some types of the triggered VLF emissions; whistler triggered emissions and artificially stimulated emissions (ASE). It is found that whistler triggered emissions excited around the upper cutoff frequencies of whistlers may be explained by the whistler mode cyclotron instability by a model distribution function inferred from satellite data. In order to see a nonlinear evolution of the whistler mode cyclotron instability, computer simulations were carried out and it is shown that the change of frequency with time of whistler triggered emissions as well as characteristics of ASE are well explained by resonant nonlinear behaviour of whistler mode cyclotron instability considered in the present paper.  相似文献   

4.
Electromagnetic waves propagating transverse to the magnetic field, containing inhomogenous and loss cone plasma, may become unstable due to the excitation of resonant proton, resonant electron and drift cyclotron instabilities. Resonant proton instability gets excited in inhomogenous plasma, irrespective of the presence of temperature anisotropy, loss cone or temperature gradient. However, the growth rate of this instability is much smaller than the other two instabilities. The maximum growth rates of resonant electron instability are enhanced with the increase of loss cone index, gradients in transverse temperature and magnetic field, and with the decrease of temperature anisotropy and gradients in density and parallel temperature. The drift cyclotron instability exists in a bounded range of wave numbers and its growth rate increases with the increase of electron temperature, density and magnetic field gradient, and with the decrease of proton temperature and temperature anisotropy. In the region of ring current for beyond plasmapause the resonant proton and resonant electron instabilities have the characterstic frequencies around 0.1Ωp and growth rates ~10?6Ωp and 10?3Ωp, respectively. In the ring current region the drift cyclotron instability is not excited whereas in the plasma sheet region the frequency and growth rate of this instability are around Ωp and 10?2Ωp, respectively. These instabilities can accelerate the ring current particles along the magnetic field lines and dump them into the auroral region.  相似文献   

5.
The MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) instrument on the Mars Express spacecraft provides both local and remote measurements of electron densities and measurements of magnetic fields in the martian ionosphere. The density measurements show a persistent level of large fluctuations, sometimes as much as a factor of three or more at high altitudes. Large magnetic field fluctuations are also observed in the same region. The power spectrums of both the density and magnetic field fluctuations have slopes on a log-log plot that are consistent with the Kolmogorov spectrum for isotropic fluid turbulence. The fractional density fluctuation, Δne/ne, of the turbulence increases with altitude, and reaches saturation, Δne/ne ∼ 1, at an altitude of about 400 km, near the nominal boundary between the ionosphere and the magnetosheath. The fluctuations are usually so large that a well-defined ionopause-like boundary between the ionosphere and the solar wind is seldom observed. Of mechanisms that could be generating this turbulence, we believe that the most likely are (1) solar wind pressure perturbations, (2) an instability in the magnetosheath plasma, such as the mirror-mode instability, or (3) the Kelvin-Helmholtz instability driven by velocity shear between the rapidly flowing magnetosheath and the ionosphere.  相似文献   

6.
Comparative study of ion cyclotron waves at Mars, Venus and Earth   总被引:1,自引:0,他引:1  
Ion cyclotron waves are generated in the solar wind when it picks up freshly ionized planetary exospheric ions. These waves grow from the free energy of the highly anisotropic distribution of fresh pickup ions, and are observed in the spacecraft frame with left-handed polarization and a wave frequency near the ion’s gyrofrequency. At Mars and Venus and in the Earth’s polar cusp, the solar wind directly interacts with the planetary exospheres. Ion cyclotron waves with many similar properties are observed in these diverse plasma environments. The ion cyclotron waves at Mars indicate its hydrogen exosphere to be extensive and asymmetric in the direction of the interplanetary electric field. The production of fast neutrals plays an important role in forming an extended exosphere in the shape and size observed. At Venus, the region of exospheric proton cyclotron wave production may be restricted to the magnetosheath. The waves observed in the solar wind at Venus appear to be largely produced by the solar-wind-Venus interaction, with some waves at higher frequencies formed near the Sun and carried outward by the solar wind to Venus. These waves have some similarity to the expected properties of exospherically produced proton pickup waves but are characterized by magnetic connection to the bow shock or by a lack of correlation with local solar wind properties respectively. Any confusion of solar derived waves with exospherically derived ion pickup waves is not an issue at Mars because the solar-produced waves are generally at much higher frequencies than the local pickup waves and the solar waves should be mostly absorbed when convected to Mars distance as the proton cyclotron frequency in the plasma frame approaches the frequency of the solar-produced waves. In the Earth’s polar cusp, the wave properties of ion cyclotron waves are quite variable. Spatial gradients in the magnetic field may cause this variation as the background field changes between the regions in which the fast neutrals are produced and where they are re-ionized and picked up. While these waves were discovered early in the magnetospheric exploration, their generation was not understood until after we had observed similar waves in the exospheres of Mars and Venus.  相似文献   

7.
PIERRE  F.  SOLOMON  J.  CORNILLEAU-WEHRLIN  N.  CANU  P.  SCIME  E. E.  BALOGH  A.  FORSYTH  R. J. 《Solar physics》1997,172(1-2):327-334
We present a study of whistler-mode wave generation and wave particle interaction in the vicinity of interplanetary shocks as observed by the Ulysses spacecraft. Generally the whistler-mode waves (measured in the frequency range 0.22–448 Hz) are observed downstream of the shocks where they persist for some hours. From the electron distribution functions (EDF) in the energy range 1.6 to 862 eV measured by the spacecraft, we compute the wave growth rate of the electromagnetic electron cyclotron and Landau instabilities for the case of oblique propagation of the waves with respect to the interplanetary magnetic field (IMF) B. In general, in agreement with the wave measurements, the instability grows mostly downstream of the shock fronts. Following the wave activity, velocity space diffusion of the electrons results in a marginally stable state with some sporadic fluctuations.  相似文献   

8.
利用北京天文台高时间和高频率分辨率的射电频谱仪对射电尖峰的测量,可以对背景等离子体参数进行的自洽诊断( 磁场,密度,温度,波矢,及非热电子的性质) 。该诊断基于电子回旋脉塞不稳定性和回旋共振吸收。最后从诊断结果和太阳日冕典型参数的比较以确定尖峰辐射的谐波数。  相似文献   

9.
A complete dispersion relation for a whistler mode wave propagation in an anisotropic warm ion-electron magnetoplasma in the presence of parallel electric field using the dispersion relation for a circularly polarized wave has been derived. The dispersion relation includes the effect of anisotropy for the ion and electron velocity distribution functions. The growth rate of electron-ion cyclotron waves for different plasma parameters observed atL = 6.6R E has been computed and the results have been discussed in detail in the light of the observed features of VLF emissions and whistlers. The role of the combination of ion-cyclotron and whistler mode electromagnetic wave propagation along the magnetic field in an anisotropic Maxwellian weakly-ionized magnetoplasma has been studied.  相似文献   

10.
A series of quasi-periodic magnetopause crossings were recorded by the MESSENGER spacecraft during its third flyby of Mercury on 29 September 2009, likely caused by a train of propagating Kelvin–Helmholtz (KH) vortices. We here revisit the observations to study the internal structure of the waves. Exploiting MESSENGER's rapid traversal of the magnetopause, we show that the observations permit a reconstruction of the structure of a rolled-up KH vortex directly from the spacecraft's magnetic field measurements. The derived geometry is consistent with all large-scale fluctuations in the magnetic field data, establishes the non-linear nature of the waves, and shows their vortex-like structure. In several of the wave passages, a reduction in magnetic field strength is observed in the middle of the wave, which is characteristic of rolled-up vortices and is related to the increase in magnetic pressure required to balance the centrifugal force on the plasma in the outer regions of a vortex, previously reported in computer simulations. As the KH wave starts to roll up, the reconstructed geometry suggests that the vortices develop two gradual transition regions in the magnetic field, possibly related to the mixing of magnetosheath and magnetospheric plasma, situated at the leading edges from the perspectives of both the magnetosphere and the magnetosheath.  相似文献   

11.
In an earlier paper, Bowers (1973), ion plasma oscillations were found to be unstable in the steady state developed by Cowley (1972) for the neutral sheet in the Earth's geomagnetic tail. In this paper a similar stability analysis is carried out but for a different steady state, suggested by Dungey, with the result that unstable waves with frequencies near the electron plasma frequency are found. In the Dungey steady state the current necessary for magnetic field reversal is carried by plasma originating from both the magnetosheath and the lobes of the tail. This modifies the steady state proposed by Alfvén and subsequently developed by Cowley in which all the current is carried by plasma from the lobes of the tail thereby fixing the cross-tail potential Φ. With magnetosheath plasma present the value of Φ is no longer fixed solely by parameters in the lobes of the tail but the cross-tail electric field is still assumed localised in the dusk region of the sheet as in the Cowley model due to the balance of charge required in the neutral sheet. The value of Φ can be expected to increase as magnetic flux is transported to the tail which inflates and causes flux annihilation because the magneto-sheath plasma in the neutral sheet has insufficient pressure to keep the two lobes of the tail apart. The Vlasov-Maxwell set of equations is perturbed and linearised enabling a critical condition for instability to be found for modes propagating across the tail. Typically, this condition requireseΦ≳KT m whereT m is the temperature of magnetosheath electrons. The instability occurs in the presence of cold plasma which hasE×B drifted into the neutral sheet from the lobes of the tail. This contrasts with the usual two stream instability which is stabilised by the cold plasma. Once precipitated the instability may be explosive provided current disruption occurs, for then a further increase in Φ will result which drives a greater range of wave numbers unstable thereby causing even more turbulence and an even larger cross-tail electric field. Because of this behaviour the instability may be a trigger for a substorm.  相似文献   

12.
A potentially promising way to gain knowledge about the internal dynamics of extrasolar planets is by remote measurement of an intrinsic magnetic field. Strong planetary magnetic fields, maintained by internal dynamo action in an electrically conducting fluid layer, are helpful for shielding the upper atmosphere from stellar wind induced mass loss and retaining water over long (Gyr) time scales. Here we present a whole planet dynamo model that consists of three main components: an internal structure model with composition and layers similar to the Earth, an optimal mantle convection model that is designed to maximize the heat flow available to drive convective dynamo action in the core, and a scaling law to estimate the magnetic field intensity at the surface of a terrestrial exoplanet. We find that the magnetic field intensity at the core surface can be up to twice the present-day geomagnetic field intensity, while the magnetic moment varies by a factor of 20 over the models considered. Assuming electron cyclotron emission is produced from the interaction between the stellar wind and the exoplanet magnetic field we estimate the cyclotron frequencies around the ionospheric cutoff at 10 MHz with emission fluxes in the range 10−4-10−7 Jy, below the current detection threshold of radio telescopes. However, we propose that anomalous boosts and modulations to the magnetic field intensity and cyclotron emission may allow for their detection in the future.  相似文献   

13.
The Io-controlled radio arcs are emissions in the decametric radio range which appear arc shaped in the time-frequency plane. Their occurrence is controlled by Io's position, so it has been for long inferred that they are powered by the Io-Jupiter electrodynamic interaction. Their frequency ranges correspond to the electron cyclotron frequencies along the Io Flux tube, so they are expected to be generated by cyclotron maser instability (CMI). The arc shape was proposed to be a consequence of the strong anisotropy of the decametric radio emissions beaming, combined with the topology of the magnetic field in the source and the observation geometry. Recent papers succeeded at reproducing the morphologies of a few typical radio arcs by modeling in three dimensions the observation geometry, using the best available magnetic field model and a beaming angle variation consistent with a loss-cone driven CMI. In the continuation of these studies, we present here the systematic modeling of a larger number of observations of the radio arcs emitted in Jupiter's southern hemisphere (including multiple arcs or arcs exhibiting abrupt changes of shape), which permits to obtain a statistical determination of the emitting field line localization (lead angle) relative to the instantaneous Io field line, and of the emitting particle velocities or energies. Variations of these parameters relative to Io's longitude are also measured and compared to the location of the UV footprints of the Io-Jupiter interaction. It is shown that the data are better organized in a reference frame attached to the UV spot resulting from the main Alfvén wing resulting from the Io-Jupiter interaction. It is proposed that the radio arcs are related to the first reflected Alfvén wing rather than to the main one.  相似文献   

14.
The magnetized solar wind carries a large amount of energy but only a small fraction of it enters the magnetosphere and powers its dynamics. Numerous observations show that the interplanetary magnetic field (IMF) is a key parameter regulating the solar wind-magnetosphere interaction. The main factor determining the amount of energy extracted from the solar wind flow by the magnetosphere is the plasma flow structure in the region adjacent to the sunward side of the magnetopause. While compared to the energy of the solar wind flow the IMF magnetic energy is relatively weak, it is considerably enhanced in a thin layer next to the dayside magnetopause variously called the plasma depletion layer or magnetic barrier. Important features of this barrier/layer are (i) a pile-up of the magnetic field with (ii) a concurrent decrease of density, (iii) enhancement of proton temperature anisotropy, (iv) asymmetry of plasma flow caused by magnetic field tension, and (v) characteristic wave emissions (ion cyclotron waves). Importantly, the magnetic barrier can be considered as an energy source for magnetic reconnection. While the steady-state magnetic barrier has been extensively examined, non-steady processes therein have only been addressed by a few authors. We discuss here two non-steady aspects related to variations of the magnetic barrier caused by (i) a north-to-south rotation of the IMF, and (ii) by pulses of magnetic field reconnection at the magnetopause. When the IMF rotates smoothly from north-to-south, a transition layer is shown to appear in the magnetosheath which evolves into a thin layer bounded by sharp gradients in the magnetic field and plasma quantities. For a given reconnection rate and calculated parameters of the magnetic barrier, we estimate the duration and length scale of a reconnection pulse as a function of the solar wind parameters. Considering a sudden decrease of the magnetic field near the magnetopause caused by the reconnection pulse, we study the relaxation process of the magnetic barrier. We find that the relaxation time is longer than the duration of the reconnection pulse for large Alfvén-Mach numbers.  相似文献   

15.
Although there is no intrinsic magnetic field at Venus, the convected interplanetary magnetic field piles up to form a magnetic barrier in the dayside inner magnetosheath. In analogy to the Earth's magnetosphere, the magnetic barrier acts as an induced magnetosphere on the dayside and hence as the obstacle to the solar wind. It consists of regions near the planet and its wake for which the magnetic pressure dominates all other pressure contributions. The initial survey performed with the Venus Express magnetic field data indicates a well-defined boundary at the top of the magnetic barrier region. It is clearly identified by a sudden drop in magnetosheath wave activity, and an abrupt and pronounced field draping. It marks the outer boundary of the induced magnetosphere at Venus, and we adopt the name “magnetopause” to address it. The magnitude of the draped field in the inner magnetosheath gradually increases and the magnetopause appears to show no signature in the field strength. This is consistent with PVO observations at solar maximum. A preliminary survey of the 2006 magnetic field data confirms the early PVO radio occultation observations that the ionopause stands at ∼250 km altitude across the entire dayside at solar minimum. The altitude of the magnetopause is much lower than at solar maximum, due to the reduced altitude of the ionopause at large solar zenith angles and the magnetization of the ionosphere. The position of the magnetopause at solar minimum is coincident with the ionopause in the subsolar region. This indicates a sinking of the magnetic barrier into the ionosphere. Nevertheless, it appears that the thickness of the magnetic barrier remains the same at both solar minimum and maximum. We have found that the ionosphere is magnetized ∼95% of the time at solar minimum, compared with 15% at solar maximum. For the 5% when the ionosphere is un-magnetized at solar minimum, the ionopause occurs at a higher location typically only seen during solar maximum conditions. These have all occurred during extreme solar conditions.  相似文献   

16.
We consider the problem of long-time storage of high-energy protons, accelerated in the process of a flare, in coronal magnetic traps. From the viewpoint of the storage, one of the most important plasma instabilities is the kinetic cyclotron instability of the Alfvén waves. We carry out a detailed theoretical analysis of the instability for typical conditions of the solar corona. It is the refraction of the Alfvén waves in combination with a drastic decrease of the instability growth rate with an increase of the angle between the directions of the wave vector and the stationary magnetic field that leads to the possibility of the long-term storage of the flare protons. Sufficient conditions of the storage are determined.  相似文献   

17.
We investigated the angular direction and polarization of the solar radio millisecond spike emission in the model in which the spike emission is due to the second harmonic instability modes driven by electron cyclotron maser of loss cone distributed electrons during the propagation of a nonlinear plasma density wave near the magnetic mirror. We found that, when the angle θ between the wave vector and the magnetic field is > 60 °, the emission is in 100% X-mode polarization; when 40 ° < θ<60 °, the emission is in 100% O-mode polarization provided the amplitude of the density wave is below a certain limit; above that limit, the polarization will fall from 100% O-mode to even the X-mode. We also found that only 0.1% of the free energy of energy carrying electrons in the source region is converted into radiation wave energy.  相似文献   

18.
We describe analysis methods to estimate parameters of electromagnetic waves based on the multi-component measurements of the DEMETER spacecraft. Using the fact that the wave magnetic field is perpendicular to the wave vector, the wave normal direction can be estimated by different methods. We use these plane-wave estimates to interpret measurements of the observed wave emissions. For instance, we use the recently developed singular value decomposition (SVD) technique. The results of the plane-wave analysis have an advantage that they often allow a straightforward interpretation. These different methods have been successfully tested with the data of previous spacecraft. All these methods are also implemented in the analysis tools designed for the analysis of the DEMETER wave measurements.We show the first results of these analysis techniques for different types of wave emissions observed on board DEMETER. Obliquely propagating right-hand polarized electromagnetic waves at a few hundreds of Hz are usually connected with a multi-ion mode structure below the local proton cyclotron frequency and with a sharp lower cutoff of left-hand polarized waves, as well as with right-hand polarized waves tunelling below the multi-ion cross-over frequency. Electron and proton whistlers are also very frequently observed on DEMETER. An unusual narrow-band emission at 140 Hz (well below the local proton cyclotron frequency) serves us as another case for a detailed analysis. We find that these waves are right-hand polarized and obliquely propagating.Using this example case, we also present analysis methods to estimate continuous distribution of wave energy density as a function of wave vector directions. These techniques of wave distribution function (WDF) analysis need both wave and particle measurements. In the analyzed case, two different methods of WDF analysis give similar results consistent with the results of the plane-wave techniques. To identify the source region we use the backward ray-tracing method. The wave normal direction obtained by the analysis of multi-component data is used for a simulation of wave propagation from the point of measurement. By this procedure, we obtain an inverse trajectory of the wave ray. We can thus follow the ray path back to the anticipated source region which is in our case located a few degrees of latitude to the South from the spacecraft position.  相似文献   

19.
The electron–cyclotron maser is a process that generates coherent radiation from plasma. In the last two decades, it has gained increasing attention as a dominant mechanism of producing high-power radiation in natural high-temperature magnetized plasmas. Originally proposed as a somewhat exotic idea and subsequently applied to include non-relativistic plasmas, the electron–cyclotron maser was considered as an alternative to turbulent though coherent wave–wave interaction which results in radio emission. However, when it was recognized that weak relativistic corrections had to be taken into account in the radiation process, the importance of the electron–cyclotron maser rose to the recognition it deserves. Here we review the theory and application of the electron–cyclotron maser to the directly accessible plasmas in our immediate terrestrial and planetary environments. In situ access to the radiating plasmas has turned out to be crucial in identifying the conditions under which the electron–cyclotron maser mechanism is working. Under extreme astrophysical conditions, radiation from plasmas may provide a major energy loss; however, for generating the powerful radiation in which the electron–cyclotron maser mechanism is capable, the plasma must be in a state where release of susceptible amounts of energy in the form of radiation is favorable. Such conditions are realized when the plasma is unable to digest the available free energy that is imposed from outside and stored in its particle distribution. The lack of dissipative processes is a common property of collisionless plasmas. When, in addition, the plasma density becomes so low that the amount of free energy per particle is large, direct emission becomes favorable. This can be expressed as negative absorption of the plasma which, like in conventional masers, leads to coherent emission even though no quantum correlations are involved. The physical basis of this formal analogy between a quantum maser and the electron–cyclotron maser is that in the electron–cyclotron maser the free-space radiation modes can be amplified directly. Several models have been proposed for such a process. The most famous one is the so-called loss-cone maser. However, as argued in this review, the loss-cone maser is rather inefficient. Available in situ measurements indicate that the loss-cone maser plays only a minor role. Instead, the main source for any strong electron–cyclotron maser is found in the presence of a magnetic-field-aligned electric potential drop which has several effects: (1) it dilutes the local plasma to such an extent that the plasma enters the regime in which the electron–cyclotron maser becomes effective; (2) it generates energetic relativistic electron beams and field-aligned currents; (3) it deforms, together with the magnetic mirror force, the electron distribution function, thereby mimicking a high energy level sufficiently far above the Maxwellian ground state of an equilibrium plasma; (4) it favors emission in the free-space RX mode in a direction roughly perpendicular to the ambient magnetic field; (5) this emission is the most intense, since it implies the coherent resonant contribution of a maximum number of electrons in the distribution function to the radiation (i.e., to the generation of negative absorption); (6) it generates a large number of electron holes via the two-stream instability, and ion holes via the current-driven ion-acoustic instability which manifest themselves as subtle fine structures moving across the radiation spectrum and being typical for the electron–cyclotron maser emission process. These fine structures can thus be taken as the ultimate identifier of the electron–cyclotron maser. The auroral kilometric radiation of Earth is taken here as the paradigm for other manifestations of intense radio emissions such as the radiation from other planets in the solar system, from exoplanets, the Sun and other astrophysical objects.  相似文献   

20.
The suprathermal plasma analyser on the geostationary satellite Geos-2 can identify magnetospheric, boundary layer and magnetosheath electron distributions around the dayside equatorial magnetopause. As examples, data from two days when magnetopause crossings occurred, 28 August 1978 and 12 November 1978, are discussed. The boundary layer electrons are intermediate in temperature and density between those in the ring current and the magnetosheath but cannot be a simple admixture of the two populations. The transition from boundary layer to magnetosheath electrons is often sudden. We believe it to be coincident with the magnetopause where the magnetic field changes from terrestrial to interplanetary.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号