首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Variations in the cosmic-ray vector anisotropy observed on Earth are closely connected with the state of the near-Earth interplanetary medium. Hourly characteristics of vector anisotropy for the period 1957–2013, which were obtained by the global survey method from the data of the worldwide network of neutron monitors, make it possible to study the relationship between the cosmic-ray anisotropy and solar wind parameters. In the present work, we have studied the connection between the equatorial component of anisotropy of cosmic rays with a rigidity of 10 GV and the following parameters: velocity and density of the solar wind; density of the interplanetary magnetic field; and cosmic-ray density variations, in which the spatial gradient of cosmic rays in the interplanetary medium is manifested. The characteristics of cosmic-ray anisotropy at various combinations of the interplanetary medium parameters are compared. The possibility of diagnosing the solar wind state from data on the cosmic-ray anisotropy is discussed.  相似文献   

2.
When the effect of a solar wind dynamic pressure pulse on the magnetospheric and ionospheric dynamics is studied, it is usually difficult to detect the effect of a sudden change in the density against the background of the other varying solar wind parameters, which often play a most pronounced role. Cases in which the solar wind plasma density gradient dominated in the dynamics of the different parameters of an interplanetary medium and its magnetic field are considered in this work. Variations in the Earth’s dayside magnetopause current caused by a change in the solar wind ion density are presented for two such cases (February 11 and January 11, 1997) based on the method developed by us previously. Variations in the dayside magnetopause current for collisions of the magnetosphere with corotating interacting flows in January 2004, studied in detail by us previously, are also presented for Saturn. The estimates are comparable with the current values in the transitional three-dimensional current systems of Saturn that were previously calculated by us.  相似文献   

3.
Coronal mass ejections (CMEs) and high-speed solar wind streams (HSS) are two solar phenomena that produce large-scale structures in the interplanetary (IP) medium. CMEs evolve into interplanetary CMEs (ICMEs) and the HSS result in corotating interaction regions (CIRs) when they interact with preceding slow solar wind. This paper summarizes the properties of these structures and describes their geoeffectiveness. The primary focus is on the intense storms of solar cycle 23 because this is the first solar cycle during which simultaneous, extensive, and uniform data on solar, IP, and geospace phenomena exist. After presenting illustrative examples of coronal holes and CMEs, I discuss the internal structure of ICMEs, in particular the magnetic clouds (MCs). I then discuss how the magnetic field and speed correlate in the sheath and cloud portions of ICMEs. CME speed measured near the Sun also has significant correlations with the speed and magnetic field strengths measured at 1 AU. The dependence of storm intensity on MC, sheath, and CME properties is discussed pointing to the close connection between solar and IP phenomena. I compare the delay time between MC arrival at 1 AU and the peak time of storms for the cloud and sheath portions and show that the internal structure of MCs leads to the variations in the observed delay times. Finally, we examine the variation of solar-source latitudes of IP structures as a function of the solar cycle and find that they have to be very close to the disk center.  相似文献   

4.
在一维球坐标系下模拟了1998年11月4日至5日3个连续日冕物质抛射(CME)在行星际空间的传播和相互作用并最终形成“复杂抛射”的日地传输过程.首先在磁流体力学(MHD)数值模拟中应用Harten总变差减小(TVD)格式,通过调节计算模型中的引力无量纲参数α、等离子体参数β和气体多方指数γ,构造出数值计算所需的初态背景,使之在拉格朗日点处L1的太阳风速度vr、质子数密度Np及质子热压力与磁压力的比值βp与ACE卫星的观测数据一致.接着仅采用速度脉冲的扰动形式,其输入的幅度和持续时间由Lasco/C2、GOES、LEAR的观测数据并结合Michalek等提出的CME“锥模型”来确定.数值计算结果得到的两个激波到达时间和ACE卫星观测值的时间误差分别是3h和4h.这表明该模型能估算续发CME在行星际空间演化后驱动激波的到达时间和大致强度,在空间天气的激波到达时间的预报方面有潜在的应用价值.  相似文献   

5.
The properties of turbulent fluctuations of the solar wind plasma near the interplanetary shock observed at September 12, 2014 by the BMSW instrument are considered. The spectra of the density fluctuations in the solar wind and their statistical characteristics up-and downstream of the shock front are analyzed. They are compared with each other and with characteristics corresponding to different turbulence models. It is shown that the spectral and statistical characteristics of the density fluctuations in the solar wind conserve their basic properties after the arrival of an interplanetary shock. Intermittency is observed both before and after the front, but its level increases on average in the second case. In both regions, the scaling of the structure functions of the density fluctuations in the solar wind differ from the scaling of the classical Kolmogorov model and can be described by the log-Poisson turbulence model. Parameterization of the scaling of the structure functions revealed the presence of filamentary structures in the solar wind plasma, which provide the density intermittency in the studied space regions.  相似文献   

6.
王明  吕建永  李刚 《地球物理学报》2014,57(11):3804-3811
利用全球磁流体力学(MHD)的模拟结果,研究了太阳风压力系数与上游太阳风参数和日下点磁层顶张角的相关性.在识别出日下点附近磁层顶位置后,通过拟合得到日下点附近的磁层顶张角.在考虑上游太阳风中的磁压和热压以及磁层顶外侧的太阳风动压的情况下,计算了太阳风压力系数.通过分析行星际磁场不同方向时太阳风动压在日地连线上与磁压和热压的转化关系,详细研究了太阳风参数和日下点磁层顶张角对太阳风压力系数的影响,得到以下相关结论:(1) 在北向行星际磁场较大(Bz≥5 nT)时,磁层顶外侧磁压占主导,南向行星际磁场时磁层顶外侧热压占主导;(2) 太阳风压力系数随着行星际磁场的增大而增大,随着行星际磁场时钟角的增大而减小;并且在行星际磁场大小和其他太阳风条件相同时,北向行星际磁场时的太阳风压力系数要大于南向行星际磁场时的;北向行星际磁场时,太阳风压力系数随着太阳风动压的增大而减小,南向行星际磁场时,太阳风压力系数随着太阳风动压的增大而增大;以上结论是对观测结果的扩展;(3) 最后,我们还发现太阳风压力系数随着日下点磁层顶张角的增大而增大.  相似文献   

7.
Sharp changes of the solar wind parameters determining the dynamic pressure jump lead to strong magnetosphere-ionosphere disturbances. Here the effect on the Earth’s ionospheric high latitudes of the solar wind dynamic pressure pulse caused only by the increase of the interplanetary plasma density under southward constant IMF is considered. We investigate reaction of the cross-polar cap potential on the increase of AL index and/or jump of the solar wind density. It is found that for the case of 10 January 1997 the main contribution to the polar cap potential drop increase gave the growth of AL index relative to the input of the solar wind density jump. We also study the influence of the solar wind density increase on the crosspolar cap potential for the quiet magnetospheric conditions. It occurred that the polar cap potential difference decreases with the great increase of the interplanetary plasma density. For the disturbed magnetosphere the main role in the polar cap potential drop increase plays increase of AL. Thus, we found the change of the cross-polar cap potential due to the AL index variations and/or the solar wind density drop even in a case when the interplanetary electric field is constant.  相似文献   

8.
Plasma and magnetic field data from the Helios 1/2 spacecraft have been used to investigate the structure of magnetic clouds (MCs) in the inner heliosphere. 46 MCs were identified in the Helios data for the period 1974–1981 between 0.3 and 1 AU. 85% of the MCs were associated with fast-forward interplanetary shock waves, supporting the close association between MCs and SMEs (solar mass ejections). Seven MCs were identified as direct consequences of Helios-directed SMEs, and the passage of MCs agreed with that of interplanetary plasma clouds (IPCs) identified as white-light brightness enhancements in the Helios photometer data. The total (plasma and magnetic field) pressure in MCs was higher and the plasma- lower than in the surrounding solar wind. Minimum variance analysis (MVA) showed that MCs can best be described as large-scale quasi-cylindrical magnetic flux tubes. The axes of the flux tubes usually had a small inclination to the ecliptic plane, with their azimuthal direction close to the east-west direction. The large-scale flux tube model for MCs was validated by the analysis of multi-spacecraft observations. MCs were observed over a range of up to 60° in solar longitude in the ecliptic having the same magnetic configuration. The Helios observations further showed that over-expansion is a common feature of MCs. From a combined study of Helios, Voyager and IMP data we found that the radial diameter of MCs increases between 0.3 and 4.2 AU proportional to the distance, R, from the Sun as R0.8 (R in AU). The density decrease inside MCs was found to be proportional to R–2.4, thus being stronger compared to the average solar wind. Four different magnetic configurations, as expected from the flux-tube concept, for MCs have been observed in situ by the Helios probes. MCs with left-and right-handed magnetic helicity occurred with about equal frequencies during 1974–1981, but surprisingly, the majority (74%) of the MCs had a south to north (SN) rotation of the magnetic field vector relative to the ecliptic. In contrast, an investigation of solar wind data obtained near Earths orbit during 1984–1991 showed a preference for NS-clouds. A direct correlation was found between MCs and large quiescent filament disappearances (disparition brusques, DBs). The magnetic configurations of the filaments, as inferred from the orientation of the prominence axis, the polarity of the overlying field lines and the hemispheric helicity pattern observed for filaments, agreed well with the in situ observed magnetic structure of the associated MCs. The results support the model of MCs as large-scale expanding quasi-cylindrical magnetic flux tubes in the solar wind, most likely caused by SMEs associated with eruptions of large quiescent filaments. We suggest that the hemispheric dependence of the magnetic helicity structure observed for solar filaments can explain the preferred orientation of MCs in interplanetary space as well as their solar cycle behavior. However, the white-light features of SMEs and the measured volumes of their interplanetary counterparts suggest that MCs may not simply be just H-prominences, but that SMEs likely convect large-scale coronal loops overlying the prominence axis out of the solar atmosphere.  相似文献   

9.
行星际背景太阳风的三维MHD数值模拟   总被引:2,自引:0,他引:2       下载免费PDF全文
杨子才  沈芳  杨易  冯学尚 《地球物理学报》2018,61(11):4337-4347
近地空间的太阳风参数预报具有重要的科学研究意义和实际应用价值,三维磁流体力学(MHD)数值模拟是太阳风参数预报的重要手段.本文建立了一套基于经验模型的三维MHD数值模型.模型的内边界设置在0.1天文单位(AU)处,在六片网格系统下利用TVD Lax-Friedrich格式求解理想MHD方程组,采用扩散法消除磁场的散度.模型以GONG的观测磁图作为输入数据,利用经验模型并结合卫星观测特征确定内边界条件.边界条件中保留了6个可调参数,以便适当调整参数使其方便适合模拟不同太阳活动期的太阳风.利用该模型分别模拟了2007年和2016年的背景太阳风,得到了太阳风速度、密度、温度和磁场强度,这些参数与ACE/WIND卫星观测符合较好.  相似文献   

10.
空间飞船Helios 1和Helios 2的磁场探测发现,太阳风高速流中磁场脉动在低频区的能谱的谱指数是随日心距离r而变化的,谱密度的空间梯度是随频率f而变化的,现有关于行星际磁场脉动的理论都不能解释上述现象。本文提出了一个计算行星际磁场能谱径向发展的理论模式。假设在不同频率的脉动之间有由低频向高频传输的能流存在,在这一基础上建立了谱方程,并得到了谱方程的数值解。数值解表明,由0.3AU至1AU,低频区的谱指数增加,而高频区的谱指数近似保持为常数(-1.6);低频区平均谱密度随着日心距离的变化为r-3.5,在高频区为r-4.1,所有这些都与观测相符合。串级的能量很可能最后传输到质子回旋频率范围,由于回旋共振而耗散,最后加热太阳风质子。这一模式有可能用来计算太阳风的加速问题。  相似文献   

11.
An analysis of the low frequency geomagnetic field fluctuations at an Antarctic (Terra Nova Bay) and a low latitude (L’Aquila, Italy) station during the Earth’s passage of a coronal ejecta on April 11, 1997 shows that major solar wind pressure variations were followed at both stations by a high fluctuation level. During northward interplanetary magnetic field conditions and when Terra Nova Bay is close to the local geomagnetic noon, coherent fluctuations, at the same frequency (3.6 mHz) and with polarization characteristics indicating an antisunward propagation, were observed simultaneously at the two stations. An analysis of simultaneous measurements from geosynchronous satellites shows evidence for pulsations at approximately the same frequencies also in the magnetospheric field. The observed waves might then be interpreted as oscillation modes, triggered by an external stimulation, extending to a major portion of the Earth’s magnetosphere.  相似文献   

12.
Ground-based geomagnetic Pc5 (2–7 mHz) pulsations, caused by the passage of dense transients (density disturbances) in the solar wind, were analyzed. It was shown that intensive bursts can appear in the density of the solar wind and its fluctuations, up to Np ~ 30–50 cm3, even during the most magnetically calm year in the past decades (2009). The analysis, performed using one of the latest methods of discrete mathematical analysis (DMA), is presented. The energy functional of a time-series fragment (called “anomaly rectification” in DMA terms) of two such events was calculated. It was established that fluctuations in the dynamic pressure (density) of the solar wind (SW) cause the global excitation of Pc5 geomagnetic pulsations in the daytime sector of the Earth’s magnetosphere, i.e., from polar to equatorial latitudes. Such pulsations started and ended suddenly and simultaneously at all latitudes. Fluctuations in the interplanetary magnetic field (IMF) have turned up to be less geoeffective in exciting geomagnetic pulsations than fluctuations in the SW density. The pulsation generation mechanisms in various structural regions of the magnetosphere were probably different. It was therefore concluded that the most probable source of ground-based pulsations are fluctuations of the corresponding periods in the SW density.  相似文献   

13.
A very strong magnetic storm of May 15, 2005, was caused by an interplanetary magnetic cloud that approached the Earths’ orbit. The sheath region of this cloud was characterized by a high solar wind density (~25–30 cm?3) and velocity (~850 km/s) and strong variations (to ~20 nT) in the interplanetary magnetic field (IMF). It has been indicated that an atypical bay-like geomagnetic disturbance was observed during the initial phase of this storm in a large longitudinal region at high latitudes: from the morning to evening sectors of the geomagnetic local time. Increasing in amplitude, the magnetic bay rapidly propagated to the polar cap latitudes up to the geomagnetic pole. An analysis of the global space-temporal dynamics of geomagnetic pulsations in the frequency band 1–6 mHz indicated that most intense oscillations were observed in the morning sector in the region of the equivalent ionospheric current at latitudes of about 72°–76°. The wavelet structure of magnetic pulsations in the polar cap and fluctuations in IMF was generally similar to the maximum at frequencies lower than 4 mHz. This can indicate that waves directly penetrated into the polar cap from the solar wind.  相似文献   

14.
The basic characteristics of the global distribution for the corona plasma and magnetic field near 2.5 Rs are analyzed with the statistical and numerical methods for 136 Carrington Rotations (CRs) covering four different phases of solar activity. By using the observational data and the velocity distribution model in the corona, the statistical average distribution of the magnetic field, density and the coronal mass outputs are analyzed for the four different phases. Then, a numerical study of the global distribution near 2.5 Rs has been made by solving a self-consistent MHD system. Finally, the solar wind speed at 1 AU is given by mapping the speed at 2.5 Rs to that near 1 AU, and the comparison of the numerical results with the observational measurements and the simulation result of the Wang–Sheeley–Arge (WSA) model are made during more than 5 years. The numerical results indicate that the global distributions on the source surface of 2.5 Rs at different phases of solar activity could be used to predict the change of the solar wind in interplanetary space.  相似文献   

15.
We investigate the inclinations of heliospheric current sheet at two sites in interplanetary space, which are generated from the same solar source. From the data of solar wind magnetic fields observed at Venus (0.72 AU) and Earth (1 AU) during December 1978–May 1982 including the solar maximum of 1981, 54 pairs of candidate sector boundary crossings are picked out, of which 16 pairs are identified as sector boundaries. Of the remainder, 12 pairs are transient structures both at Venus and Earth, and 14 pairs are sector boundaries at one site and have transient structures at the other site. It implies that transient structures were often ejected from the coronal streamer belt around the solar maximum. For the 16 pairs of selected sector boundaries, we determine their normals by using minimum variance analysis. It is found that most of the normal azimuthal angles are distributed between the radial direction and the direction perpendicular to the spiral direction both at Venus and Earth. The normal elevations tend to be smaller than ≈45° with respect to the solar equatorial plane, indicating high inclinations of the heliospheric current sheet, in particular at Earth. The larger scatter in the azimuth and elevation of normals at Venus than at Earth suggests stronger effects of the small-scale structures on the current sheet at 0.72 AU than at 1 AU. When the longitude difference between Venus and Earth is small (<40° longitudinally), similar or the same inclinations are generally observed, especially for the sector boundaries without small-scale structures. This implies that the heliospheric current sheet inclination tends to be maintained during propagation of the solar wind from 0.72 AU to 1 AU. Detailed case studies reveal that the dynamic nature of helmet streamers causes variations of the sector boundary structure.  相似文献   

16.
An extended structure-function model is developed by including the new effect in the p-model of Meneveau and Sreenivasan which shows that the averaged energy cascade rate changes with scale, a situation which has been found to prevail in nonfullydeveloped turbulence in the inner solar wind. This model is useful for the small-scale fluctuations in the inner heliosphere, where the turbulence is not fully developed and cannot be explained quantitatively by any of the previous intermittency turbulence models. With two model parameters, the intrinsic index of the energy spectrum <alpha>, and the fragmentation fraction P 1, the model can fit, for the first time, all the observed scaling exponents of the structure functions, which are calculated for time lags ranging from 81 s to 0.7 h from the Helios solar wind data. From the cases we studied we cannot establish for P 1 either a clear radial evolution trend, or a solar-wind-speed or stream-structure dependence or a systematic anisotropy for both the flow velocity and magnetic field component fluctuations. Generally, P 1 has values between 0.7 and 0.8. However, in some cases in low-speed wind P 1 has somewhat higher values for the magnetic components, especially for the radial component. In high-speed wind, the inferred intrinsic spectral indices (<alpha>) of the velocity and magnetic field components are about equal, while the experimental spectral indices derived from the observed power spectra differ. The magnetic index is somewhat larger than the index of the velocity spectrum. For magnetic fluctuations in both high- and low-speed winds, the intrinsic exponent <alpha> has values which are near 1.5, while the observed spectral exponent has much higher values. In the solar wind with considerable density fluctuations near the interplanetary current sheet near 1 AU, it is found that P 1 has a comparatively high value of 0.89 for V x . The impact of these results on the understanding of the nature of solar wind fluctuations is discussed, and the limitations in using structure functions to study intermittency are also described.  相似文献   

17.
Recent observations have quantified the auroral wind O+ outflow in response to magnetospheric inputs to the ionosphere, notably Poynting energy flux and precipitating electron density. For moderate to high activity periods, ionospheric O+ is observed to become a significant or dominant component of plasma pressure in the inner plasma sheet and ring current regions. Using a global circulation model of magnetospheric fields and its imposed ionospheric boundary conditions, we evaluate the global ionospheric plasma response to local magnetospheric conditions imposed by the simulation and evaluate magnetospheric circulation of solar wind H+, polar wind H+, and auroral wind O+. We launch and track the motions of millions of test particles in the global fields, launched at randomly distributed positions and times. Each particle is launched with a flux weighting and perpendicular and parallel energies randomly selected from defined thermal ranges appropriate to the launch point. One sequence is driven by a two-hour period of southward interplanetary magnetic field for average solar wind intensity. A second is driven by a 2-h period of enhanced solar wind dynamic pressure for average interplanetary field. We find that the simulated ionospheric O+ becomes a significant plasma pressure component in the inner plasma sheet and outer ring current region, particularly when the solar wind is intense or its magnetic field is southward directed. We infer that the reported empirical scalings of auroral wind O+ outflows are consistent with a substantial pressure contribution to the inner plasma sheet and plasma source surrounding the ring current. This result violates the common assumption that the ionospheric load is entirely confined to the F layer, and shows that the ionosphere is often an important dynamic element throughout the magnetosphere during moderate to large solar wind disturbances.  相似文献   

18.
We investigate the features of the planetary distribution of wave phenomena (geomagnetic pulsations) in the Earth’s magnetic shell (the magnetosphere) during a strong geomagnetic storm on December 14–15, 2006, which is untypical of the minimum phase of solar activity. The storm was caused by the approach of the interplanetary magnetic cloud towards the Earth’s magnetosphere. The study is based on the analysis of 1-min data of global digital geomagnetic observations at a few latitudinal profiles of the global network of ground-based magnetic stations. The analysis is focused on the Pc5 geomagnetic pulsations, whose frequencies fall in the band of 1.5–7 mHz (T ~ 2–10 min), on the fluctuations in the interplanetary magnetic field (IMF) and in the solar wind density in this frequency band. It is shown that during the initial phase of the storm with positive IMF Bz, most intense geomagnetic pulsations were recorded in the dayside polar regions. It was supposed that these pulsations could probably be caused by the injection of the fluctuating streams of solar wind into the Earth’s ionosphere in the dayside polar cusp region. The fluctuations arising in the ionospheric electric currents due to this process are recorded as the geomagnetic pulsations by the ground-based magnetometers. Under negative IMF Bz, substorms develop in the nightside magnetosphere, and the enhancement of geomagnetic pulsations was observed in this latitudinal region on the Earth’s surface. The generation of these pulsations is probably caused by the fluctuations in the field-aligned magnetospheric electric currents flowing along the geomagnetic field lines from the substorm source region. These geomagnetic pulsations are not related to the fluctuations in the interplanetary medium. During the main phase of the magnetic storm, when fluctuations in the interplanetary medium are almost absent, the most intense geomagnetic pulsations were observed in the dawn sector in the region corresponding to the closed magnetosphere. The generation of these pulsations is likely to be associated with the resonance of the geomagnetic field lines. Thus, it is shown that the Pc5 pulsations observed on the ground during the magnetic storm have a different origin and a different planetary distribution.  相似文献   

19.
This paper reviews quantitative analysis results of the energy transfer and dissipation processes in the GUMICS-4 global MHD simulation. Reconnection power dissipating magnetic energy, dynamo power transferring energy from plasma to the field, and energy flux transport across the magnetopause surface are all examined separately and shown to yield consistent results. This is used to argue that magnetic reconnection is the process controlling the energy transfer, even though it is not localized near the reconnection line. The most important factors controlling the reconnection efficiency are the interplanetary magnetic field (IMF) orientation and the solar wind speed, while the IMF magnitude and solar wind density play a lesser role. During northward IMF, the reconnection efficiency is larger for high speed and low IMF than for low speed and high IMF magnitude, even though the solar wind electric field in both cases is the same. Moreover, increasing pressure by increasing density has a different effect from equal increase of pressure by increasing the solar wind speed. Comparison with statistical observational results shows that the simulation results are in qualitative agreement with the observations, which significantly increases our confidence in interpreting the simulation results.  相似文献   

20.
The general features of the region of interaction of the solar wind with the ionosphere of Venus and Mars are compared using data obtained with the Mariner 5 and the Pioneer Venus Orbiter (PVO) spacecraft for Venus and with the Phobos II, the Mars Global Surveyor (MGS) and the Mars Express spacecraft for Mars. Despite the overall weak intrinsic global magnetic field that is present in both planets there are significant differences in the manner in which the interplanetary magnetic field accumulates and is organized around and within their ionosphere. Such differences are unrelated to the crustal magnetic field remnants inferred from the MGS measurements around Mars. In fact, while in Venus and Mars there is a region in which the magnetic field becomes enhanced as it piles up in their plasma environment it is shown that such a region exhibits different regimes with respect to changes in the ion composition measured outside and within the ionosphere. At Venus the region of enhanced magnetic field intensity occurs in general above the ionopause which represents the boundary across which there is a change in the ion composition with dominant solar wind protons above and planetary O+ ions below. At Mars the region of enhanced magnetic field is located below a magnetic pileup boundary across which there is also a comparable change in the ion composition (solar wind protons above and planetary O+ ions below). It is argued that this difference in the relative position of the region of enhanced magnetic field with respect to that of a plasma boundary that separates different ion populations results from the peculiar response of the ionosphere of each planet to the oncoming solar wind dynamic pressure. While at Venus the peak ionospheric thermal pressure is in general sufficient to withhold the incident solar wind kinetic pressure there is a different response in Mars where the peak ionospheric thermal pressure is in general not large enough to deviate the solar wind. In this latter case the ionosphere is unable to force the solar wind to move around the ionosphere and as a result the oncoming electron population can reach low altitudes where it is influenced by neutral atmospheric particles (the solar wind proton population is replaced at the magnetic pileup boundary which marks the upper extent of the region where the interplanetary magnetic field becomes enhanced). Peculiar conditions are expected near the magnetic polar regions and over the terminator plane where the solar wind is directed along the sides of the planet.  相似文献   

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

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