首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 343 毫秒
1.
Khabarova  O.  Zastenker  G. 《Solar physics》2011,270(1):311-329
Analysis of the Interball-1 spacecraft data (1995 – 2000) has shown that the solar wind ion flux sometimes increases or decreases abruptly by more than 20% over a time period of several seconds or minutes. Typically, the amplitude of such sharp changes in the solar wind ion flux (SCIFs) is larger than 0.5×108 cm−2 s−1. These sudden changes of the ion flux were also observed by the Solar Wind Experiment (SWE), on board the Wind spacecraft, as the solar wind density increases and decreases with negligible changes in the solar wind velocity. SCIFs occur irregularly at 1 AU, when plasma flows with specific properties come to the Earth’s orbit. SCIFs are usually observed in slow, turbulent solar wind with increased density and interplanetary magnetic field strength. The number of times SCIFs occur during a day is simulated using the solar wind density, magnetic field, and their standard deviations as input parameters for a period of five years. A correlation coefficient of ∼0.7 is obtained between the modelled and the experimental data. It is found that SCIFs are not associated with coronal mass ejections (CMEs), corotating interaction regions (CIRs), or interplanetary shocks; however, 85% of the sector boundaries are surrounded by SCIFs. The properties of the solar wind plasma for days with five or more SCIF observations are the same as those of the solar wind plasma at the sector boundaries. One possible explanation for the occurrence of SCIFs (near sector boundaries) is magnetic reconnection at the heliospheric current sheet or local current sheets. Other probable causes of SCIFs (inside sectors) are turbulent processes in the slow solar wind and at the crossings of flux tubes.  相似文献   

2.
Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.  相似文献   

3.
At the Sun-Earth distance of one astronomical unit (1 AU), the solar wind is known to be strongly supersonic and super Alfvenic with Mach and Alfven numbers being on average 12 and 9 respectively. Also, solar wind densities (average ∼10cm-3) and velocities (average ∼450kms-1) at 1AU, are known to be inversely correlated with low velocities having higher than average densities andvice versa. However, on May 11 and 12 1999 the Earth was engulfed by an unusually low density (< 0.1cm-3) and low velocity (< 350km s-1) solar wind with an Alfven Mach number significantly less than 1. This was a unique low-velocity, low-density, sub-Alfvénic solar wind flow which spacecraft observations have shown lasted more than 24 hours. One consequence of this extremely tenuous solar wind was a spectacular expansion of the Earth’s magnetosphere and bow shock. The expanding bow shock was observed by several spacecraft and reached record upstream distances of nearly 60 Earth radii, the lunar orbit. The event was so dramatic that it has come to be known asthe solar wind disappearance event. Though extensive studies of this event were made by many authors in the past, it has only been recently shown that the unusual solar wind flows characterizing this event originated from a small coronal hole in the vicinity of a large active region on the Sun. These recent results have put to rest speculation that such events are associated with global phenomenon like the periodic solar polar field reversal that occurs at the maximum of each solar cycle. In this paper we revisit the 11 May 1999 event, look at other disappearance events that have ocurred in the past, examine the reasons why speculations about the association of such events with global phenomena like solar polar field reversals were made and also examine the role of transient coronal holes as a possible solar source for such events.  相似文献   

4.
We present a study of the plasma properties inside and dynamics of the low-latitude boundary layer (LLBL)/cusp during the ICME event on 7 November 2004 based on data from the four Cluster spacecraft. The interplanetary magnetic field (IMF) is predominantly strongly northward, up to 50 nT, with some short-duration rotations. The observed LLBL/cusp is very thick (∼6 – 7° invariant latitude (ILAT)) and migrates equatorward with rates of 0.55° and 0.04° ILAT per minute during quick southward IMF rotations and stable northward IMF, respectively. The LLBL/cusp observed by Cluster 1 and Cluster 4 is in a fast transition between different states and is populated by different types of plasma injection, presumably coming from multiple reconnection sites. During a period of extremely northward IMF, signatures of pulsed dual reconnection inside the LLBL/cusp are observed by Cluster 3, suggesting that at least part of the LLBL/cusp is on closed field lines. However, analysis of the ion data implies that the boundary layer is formed in the dawn sector of the magnetosphere and does not slowly convect from the dayside as has been suggested previously. A statistical study of the location of the LLBL/cusp equatorward boundary during the ICME events on 28 – 29 October 2003 and 7 – 10 November 2004 is performed. During extreme conditions the LLBL/cusp position is offset by −7° ILAT from the location under normal conditions, which might be explained by the influence of the high solar wind dynamic pressure. The LLBL/cusp moves equatorward with increasing southward and northward IMF. However, the LLBL/cusp position under strong southward IMF is more poleward than expected from previous studies, which could indicate some saturation in the dayside reconnection process or enhancement of the nightside reconnection rate. The LLBL/cusp position under strong northward IMF is extremely low and does not agree with the location predicted in previous studies. For the events with solar wind dynamic pressure >10 nPa, the LLBL/cusp position does not depend on the solar wind dynamic pressure. This might indicate some saturation in the mechanism of how the LLBL/cusp location depends on the solar wind dynamic pressure.  相似文献   

5.
An experimental simulation of planetary magnetospheres is being developed to investigate the formation of collisionless shocks and their effects. Two experimental situations are considered. In both, the solar wind is simulated by laser ablation plasmas. In one case, the “solar wind” flows across the magnetic field of a high-current discharge. In the other, a transverse magnetic field is embedded in the plasma flow, which interacts with a conductive obstacle. The ablation plasma is created using the “Tomcat” laser, currently emitting 5 J in a 6 ns pulse at 1 μm wavelength and irradiance above 1013 W/cm2. The “Zebra” z-pinch generator, with load current up to 1 MA and voltage up to 3.5 MV produces the magnetic fields. Hydrodynamic modeling is used to estimate the plasma parameters achievable at the front of the plasma flow and to optimize the experiment design. Particle-in-cell simulations reveal details of the interaction of the “solar wind” with an external magnetic field, including flow collimation and heating effects at the stopping point. Hybrid simulations show the formation of a bow shock at the interaction of a magnetized plasma flow with a conductor. The plasma density and the embedded field have characteristic spatial modulations in the shock region, with abrupt jumps and fine structure on the skin depth scale.  相似文献   

6.
The solar magnetic field maps every point in the corona to a corresponding place on the solar surface. Identifying the magnetic connection map is difficult at low latitudes near the heliospheric current sheet, but remarkably simple in coronal hole interiors. We present a simple analytic magnetic model (‘pseudocurrent extrapolation’) that reproduces the global structure of the corona, with significant physical advantages over other nearly analytic models such as source-surface potential field extrapolation. We use the model to demonstrate that local horizontal structure is preserved across altitude in the central portions of solar coronal holes, up to at least 30 Rs, in agreement with observations. We argue that the preserved horizontal structure may be used to track the magnetic footpoint associated with the location of a hypothetical spacecraft traveling through the solar corona, to relate in situ measurements of the young solar wind at ∼10–30 Rs to particular source regions at the solar surface. Further, we discuss the relationship between readily observable geometrical distortions and physical parameters of interest such as the field-aligned current density.  相似文献   

7.
Using a well-known method for calculating the propagation of waves in an inhomogeneous medium, we have managed to reduce the problem of wave propagation in pulsar magnetospheres to a system of two ordinary differential equations that allow the polarization characteristics of the radio emission to be quantitatively described for any magnetic field structure and an arbitrary density profile of the outflowing plasma. We confirm that for ordinary pulsars (period P ∼ 1 s, magnetic field B 0 ∼ 1012 G, particle production multiplicity parameter λ ∼ 104), the polarization is formed inside the light cylinder at a distance of the order of a thousand neutron star radii. For reasonable magnetic field strengths and plasma densities on the emission propagation path, the degree of circular polarization is found to be ∼5–20%, in good agreement with observations.  相似文献   

8.
Theoretical model, explaining a phenomenon of formation of Intensive Magnetic Flux Tube (IMFT) in a converging flow of partially ionized solar photospheric plasma is considered. Special attention is paid to the fact of weak ionization (n/n n ∼ 10-4) of plasma in the photosphere. The cases of 2D magnetic slab and cylindric magnetic tube are considered. It was shown that in a converging flow of photospheric plasma thin magnetic tubes, or slabs with the characteristic scale L 0 ∼ (1 ÷ 5) ċ 107 cm and magnetic field 1000 ÷ 2000 G can be generated. By this 2D magnetic slabs could be unstable with respect to an exchange instability and appear as an intermediate step during IMFT formation on the boundary of two supergranulation cells. Formation of compact strong magnetic field structures, and their energy balance are discussed. Stationary Joule energy dissipation taking place on the photospheric levels in the models of magnetic slab or IMFT under consideration increases towards the periphery of these objects and can exceed radiation looses. This can cause the occurrence of magnetic tubes with hot external envelopes, and modification of plasma temperature and density distribution, with respect to ones in a quiet atmosphere. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

9.
Using minimum variance analysis of the circular mapping data from the Mars Global Surveyor (MGS) spacecraft during four selected weeks of observation, we identify 360 magnetic field structures in the Martian topside ionosphere with characteristic signatures of flux ropes. Physical parameters including size, peak field strength, helicity, orientation, and external conditions at the time of each observation are compiled for the events in each population. We observe that Martian flux ropes typically have a peak field amplitude of ∼15 nT and a diameter of ∼80–100 km assuming they are stationary. Flux ropes tend to be aligned approximately parallel to the planetary surface, and perpendicular to the direction from which the solar wind flows. They are more frequently observed during times of low solar wind pressure, but do not show a clear preference for a particular Interplanetary Magnetic Field (IMF) draping direction. Flux rope characteristics of peak field amplitude, diameter, and helicity vary with solar zenith angle. Amplitudes tend to be higher during periods of high solar wind pressure. The events are sorted into three populations based on the location at which they were observed, possibly corresponding to distinct formation mechanisms. Flux ropes observed in eclipse tend to have smaller peak amplitudes and are larger than those observed in sunlight, and are less likely to be oriented parallel to the planetary surface. Proximity to crustal fields does not appear to influence the characteristics of flux ropes observed at the 400 km spacecraft altitude. The frequent observation of flux rope structures near Mars in a variety of locations suggests that the low-altitude plasma environment is quite dynamic, with magnetic shear playing a prominent role in determining magnetic field structure near the planet.  相似文献   

10.
We introduce a method for constructing large-scale (0.25 AU) interplanetary magnetic field lines using only solar wind velocity from well-separated appropriately located spacecraft. The technique is based on labeling the field lines at each spacecraft with their coronal connection longitudes calculated in the EQRH (extrapolated quasi-radial hypervelocity) approximation (Nolte and Roelof, 1973). Even though the EQRH approximation is most applicable to quasi-steady solar wind, we propose that it should also be satisfactorily accurate for moderately evolving conditions. For strongly evolving conditions (e.g., flare-associated plasma) we propose a straightforward correction based on the inferred coronal longitudinal velocity profile. To illustrate the multispacecraft EQRH technique, we perform a calculation in which the interplanetary field lines in a model evolving solar wind disturbance are deduced from model observations at separated spacecraft. Since the expected agreement is found, we use data from Pioneers 8 and 9 and Vela to construct field lines for an unusually quiet period (April 26–30, 1969) and for a flare-associated disturbance accompanied by a Forbush decrease (March 23–25, 1969). The deduced field lines (even though strongly distorted by the disturbance), order the onsets of the Forbush decrease at the separated spacecraft, and the interplanetary plasma and field structures correspond to equatorial structures apparent in H synoptic charts of chromospheric magnetic features.  相似文献   

11.
We compared the variability of coronal hole (CH) areas (determined from daily GOES/SXI images) with solar wind (daily ACE data) and geomagnetic parameters for the time span 25 January 2005 until 11 September 2005 (late declining phase of solar cycle 23). Applying wavelet spectral analysis, a clear 9-day period is found in the CH time series. The GOES/SXI image sequence suggests that this periodic variation is caused by a mutual triangular distribution of CHs ∼120° apart in longitude. From solar wind parameters a 9-day periodicity was obtained as well, simultaneously with the 9-day period in the CH area time series. These findings provide strong evidence that the 9-day period in solar wind parameters, showing up as higher harmonic of the solar rotation frequency, is caused by the “periodic” longitudinal distribution of CHs on the Sun recurring for several solar rotations. The shape of the wavelet spectrum from the Dst index matches only weakly with that from the CH areas and is more similar to the wavelet spectrum of the solar wind magnetic field magnitude. The distinct 9-day period does not show up in sunspot group areas which gives further evidence that the solar wind modulation is strongly related to CH areas but not to active region complexes. The wavelet power spectra for the whole ACE data range (∼1998 – 2006) suggest that the 9-day period is not a singular phenomenon occurring only during a specific time range close to solar minimum but is occasionally also present during the maximum and decay phase of solar cycle 23. The main periods correspond to the solar rotation (27d) as well as to the second (13.5d) and third (9d) harmonic. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.  相似文献   

12.
The magnetic field of Mercury and the structure and dynamics of Mercury's magnetosphere, which will be studied by the spacecraft orbiting Mercury, are strongly influenced by the interaction of the solar wind with Mercury. In order to understand the internal magnetic field, it will be necessary to correct the observations of the external field for the distortions produced by the solar wind. Understanding of the solar wind interaction with Mercury is essential for understanding the structure and dynamics of the magnetosphere and phenomena such as magnetic storms. Helios 1 and 2 made a number of passes in the region traversed by the orbit of Mercury, and each pass provided a sample of the solar wind environment of Mercury. This paper reviews the plasma and magnetic field observations from Helios that provide a general basis for interpreting the observations of Mercury that will be made by orbiting spacecraft. The variables that govern the structure and dynamics of the magnetospheres of Mercury and Earth are approximately 5–10 times larger at Mercury than at Earth. Thus, the solar wind interaction with Mercury will be much stronger than the interaction with Earth. Moreover, the solar wind at Mercury is probably more variable than that at Earth. There is a clear need for measurements of the solar wind during the approach of spacecraft to Mercury and while they are in orbit around Mercury.  相似文献   

13.
In January 2004 the dust instrument on the Cassini spacecraft detected the first high-velocity grain expelled from Saturn - a so-called stream particle. Prior to Cassini’s arrival at Saturn in July 2004 the instrument registered 801 faint impacts, whose impact signals showed the characteristic features of a high-velocity impact by a tiny grain. The impact rates as well as the directionality of the stream particles clearly correlate with the sector structure of the interplanetary magnetic field (IMF). The Cosmic Dust Analyser (CDA) registered stream particles dominantly during periods when the IMF direction was tangential to the solar wind flow and in the prograde direction. This finding provides clear evidence for a continuous outflow of tiny dust grains with similar properties from the saturnian system. Within the compressed part of co-rotating interaction regions (CIRs) of the IMF, characterized by enhanced magnetic field strength and compressed solar wind plasma, CDA observed impact bursts of faster stream particles. We find that the bursts result from the stream particles being sped up inside the compressed CIR regions. Our analysis of the stream-particle dynamics inside rarefaction regions of the IMF implies that saturnian stream particles have sizes between 2 and 9 nm and exit the saturnian systems closely aligned with the planet’s ring plane with speeds in excess of 70 km s−1.  相似文献   

14.
Strong ultraviolet radiation from the Sun ionizes the upper atmosphere of Venus, creating a dense ionosphere on the dayside of the planet. In contrast to Earth, the ionosphere of Venus is not protected against the solar wind by a magnetic field. However, the interaction between charged ionospheric particles and the solar wind dynamic and magnetic pressure creates a pseudo-magnetosphere which deflects the solar wind flow around the planet (Schunk and Nagy, 1980). The combination of changing solar radiation and solar wind intensities leads to a highly variable structure and plasma composition of the ionosphere. The instrumentation of the Venus Express spacecraft allows to measure the magnetic field (MAG experiment) as well as the electron energy spectrum and the ion composition (ASPERA-4 experiment) of the upper ionosphere and ionopause. In contrast to the earlier Pioneer Venus Orbiter (PVO) measurements which were conducted during solar maximum, the solar activity was very low in the period 2006-2009. A comparison with PVO allows for an investigation of ionospheric properties under different solar wind and EUV radiation conditions. Observations of MAG and ASPERA have been analyzed to determine the positions of the photoelectron boundary (PEB) and the “magnetopause” and their dependence on the solar zenith angle (SZA). The PEB was determined using the ELS observations of ionospheric photoelectrons, which can be identified by their specific energy range. It is of particular interest to explore the different magnetic states of the ionosphere, since these influence the local plasma conductivity, currents and probably the escape of electrons and ions. The penetration of magnetic fields into the ionosphere depends on the external conditions as well as on the ionospheric properties. By analyzing a large number of orbits, using a combination of two different methods, we define criteria to distinguish between the so-called magnetized and unmagnetized ionospheric states. Furthermore, we confirm that the average magnetic field inside the ionosphere shows a linear dependence on the magnetic field in the region directly above the PEB.  相似文献   

15.
Energetic Particle Fluxes during the Bastille Day Solar Eruption   总被引:2,自引:0,他引:2  
Mäkelä  Pertti  Torsti  Jarmo 《Solar physics》2001,204(1-2):213-225
We report on our observations of solar energetic particle fluxes of p, He, C, O, Ne, Mg, Si, and Fe ions measured by the Energetic and Relativistic Nucleon and Electron (ERNE) experiment associated with the Bastille Day solar flare and coronal mass ejection (CME) on 14 July 2000. We observed two clear maxima of the Fe/O ratio at the energies 8.5–15 MeV nucl−1. The first Fe/O maximum occurred ∼ 3 hours after the beginning of the particle event, and the second maximum ∼ 22 hours after the first one at the arrival of the shock associated with the Bastille Day eruption. We also observed a change in the energy spectrum of oxygen concurrent with a change in the direction of the interplanetary magnetic field at the start of the second enhancement of the Fe/O ratio. We propose an interpretation of the particle event where observed interplanetary particle fluxes are associated with two different particle sources near the Sun and in interplanetary space. We suggest that heavy ions observed during the first period of the Fe/O enhancement were released when a coronal shock reached a magnetic foot point connected to 1 AU. The second maximum of Fe/O occurred when spacecraft encountered Fe-rich material stored in magnetic field flux tubes early in the event and was possibly reaccelerated by the interplanetary shock.  相似文献   

16.
The structure of the AGN object 1803+784 has been investigated at a wavelength of 7 mm with a limiting angular resolution reaching 20 μas. The ejector nozzle surrounded by a ring structure, an accretion disk, has been identified. The nozzle size is ∼0.1 pc, the diameter of the ring structure is ∼1.4 pc, and its width is ∼0.25 pc. The reaction of the plasma flow produces a multimode precession responsible for the conical helical structure of the jet with a variable step and a curved axis. The viewing angle of the flow ejection is ∼40°. The central part of the ejected flow moving along the axis accelerates to a relativistic velocity. The apparent velocity reaches 12 s at a distance of ∼1 mas or ∼6 pc from the ejector. The outer part of the flow moves along a helix around a high-velocity component whose step is a factor of 4 smaller, because the longitudinal velocity is relatively low. The plasma is ejected almost toward the observer, as confirmed by its high brightness temperature T b ≈ 8 × 1013 K and highly beamed emission. The polarized emission from the nozzle is axisymmetric. The orientation of the polarization of the flow along the whole length is aligned with the direction of its motion, suggesting the excitation of a ring magnetic field around it and self-focusing.  相似文献   

17.
The MESSENGER spacecraft flyby of Mercury on 14 January 2008 provided a new opportunity to study the intrinsic magnetic field of the innermost planet and its interaction with the solar wind. The model presented in this paper is based on the solution of the three-dimensional, bi-fluid equations for solar wind protons and electrons in the absence of mass loading. In this study we provide new estimates of Mercury’s intrinsic magnetic field and the solar wind conditions that prevailed at the time of the flyby. We show that the location of the boundary layers and the strength of the magnetic field along the spacecraft trajectory can be reproduced with a solar wind ram pressure Psw = 6.8 nPa and a planetary magnetic dipole having a magnitude of 210 RM3 − nT and an offset of 0.18 RM to the north of the equator, where RM is Mercury’s radius. Analysis of the plasma flow reveals the existence of a stable drift belt around the planet; such a belt can account for the locations of diamagnetic decreases observed by the MESSENGER Magnetometer. Moreover, we determine that the ion impact rate at the northern cusp was four times higher than at the southern cusp, a result that provides a possible explanation for the observed north-south asymmetry in exospheric sodium in the neutral tail.  相似文献   

18.
The MESSENGER mission to Mercury, to be launched in 2004, will provide an opportunity to characterize Mercury's internal magnetic field during an orbital phase lasting one Earth year. To test the ability to determine the planetary dipole and higher-order moments from measurements by the spacecraft's fluxgate magnetometer, we simulate the observations along the spacecraft trajectory and recover the internal field characteristics from the simulated observations. The magnetic field inside Mercury's magnetosphere is assumed to consist of an intrinsic multipole component and an external contribution due to magnetospheric current systems described by a modified Tsyganenko 96 model. Under the axis-centered-dipole approximation without correction for the external field the moment strength is overestimated by ∼4% for a simulated dipole moment of , and the error depends strongly on the magnitude of the simulated moment, rising as the moment decreases. Correcting for the external field contributions can reduce the error in the dipole term to a lower limit of ∼1-2% without a solar wind monitor. Dipole and quadrupole terms, although highly correlated, are then distinguishable at the level equivalent to an error in the position of an offset dipole of a few tens of kilometers. Knowledge of the external magnetic field is therefore the primary limiting factor in extracting reliable knowledge of the structure of Mercury's magnetic field from the MESSENGER observations.  相似文献   

19.
The behavior of solar energetic particles (SEPs) in a shock – magnetic cloud interacting complex structure observed by the Advanced Composition Explorer (ACE) spacecraft on 5 November 2001 is analyzed. A strong shock causing magnetic field strength and solar wind speed increases of about 41 nT and 300 km s−1, respectively, propagated within a preceding magnetic cloud (MC). It is found that an extraordinary SEP enhancement appeared at the high-energy (≥10 MeV) proton intensities and extended over and only over the entire period of the shock – MC structure passing through the spacecraft. Such SEP behavior is much different from the usual picture that the SEPs are depressed in MCs. The comparison of this event with other top SEP events of solar cycle 23 (2000 Bastille Day and 2003 Halloween events) shows that such an enhancement resulted from the effects of the shock – MC complex structure leading to the highest ≥10 MeV proton intensity of solar cycle 23. Our analysis suggests that the relatively isolated magnetic field configuration of MCs combined with an embedded strong shock could significantly enhance the SEP intensity; SEPs are accelerated by the shock and confined into the MC. Further, we find that the SEP enhancement at lower energies happened not only within the shock – MC structure but also after it, probably owing to the presence of a following MC-like structure. This is consistent with the picture that SEP fluxes could be enhanced in the magnetic topology between two MCs, which was proposed based on numerical simulations by Kallenrode and Cliver (Proc. 27th ICRC 8, 3318, 2001b).  相似文献   

20.
The POLAR Investigation of the Sun (POLARIS) mission uses a combination of a gravity assist and solar sail propulsion to place a spacecraft in a 0.48 AU circular orbit around the Sun with an inclination of 75° with respect to solar equator. This challenging orbit is made possible by the challenging development of solar sail propulsion. This first extended view of the high-latitude regions of the Sun will enable crucial observations not possible from the ecliptic viewpoint or from Solar Orbiter. While Solar Orbiter would give the first glimpse of the high latitude magnetic field and flows to probe the solar dynamo, it does not have sufficient viewing of the polar regions to achieve POLARIS’s primary objective: determining the relation between the magnetism and dynamics of the Sun’s polar regions and the solar cycle.
T. AppourchauxEmail:
  相似文献   

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

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