首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Longitudinal and local time variations in the structure of the equatorial anomaly under high solar activity in the equinox are considered according to the Intercosmos-19 topside sounding data. It is shown that the anomaly begins to form at 0800 LT, when the southern crest is formed. The development of the equatorial anomaly is associated with well-known variations in the equatorial ionosphere: a change in the direction of the electric field from the west to the east, which causes vertical plasma drift W (directed upward) and the fountain effect. At 1000 LT, both anomaly crests appear, but they become completely symmetrical only by 1400 LT. The average position of the crests increases from I = 20° at 1000 LT to I = 28° at 1400 LT. The position of the crests is quite strong, sometimes up to 15°, varies with longitude. The foF2 value above the equator and the equatorial anomaly intensity (EAI) at 1200–1400 LT vary with the longitude according to changes in the vertical plasma drift velocity W. At this time, four harmonics are observed in the longitudinal variations of W, foF2, and EAI. The equatorial anomaly intensity increases to the maximum 1.5–2 h after the evening burst in the vertical plasma drift velocity. Longitudinal variations of foF2 for 2000–2200 LT are also associated with corresponding variations in the vertical plasma drift velocity. The equatorial anomaly intensity decreases after the maximum at 2000 LT and the crests decrease in size and shift towards the equator, but the anomaly is well developed at midnight. On the contrary, after midnight, foF2 maxima in the region of the anomaly crests are farther from the equator, but this is obviously associated with the action of the neutral wind. At 0200 LT, in contrast to the morning hours, only the northern crest of the anomaly is clearly pronounced. Thus, in the case of high solar activity during the equinoxes, a well-defined equatorial anomaly is observed from 1000 to 2400 LT. It reaches the maximum at 2000 LT.  相似文献   

2.
This paper presents the results of simultaneous observations of narrow-band noise VLF emissions in the frequency range 4–10 kHz at Kannuslehto ground station in Northern Finland and by Van Allen Probes (previously RBSP) in the equatorial part of the magnetosphere. The event of December 25, 2015, is considered. During the event, narrow-band noise VLF emissions were detected on the Earth in two frequency ranges, f = 3.5–6 kHz and f = 8–10 kHz, between 1100 and 1300 UT. Narrow-band VLF emissions in the equatorial zone were also observed during that time by the RBSP-B satellite; their frequency was close to the electron equatorial half-gyrofrequency and gradually increased from 3 to 11 kHz during the satellite motion from L = 5.0 to L = 3.0. Analysis of the fine structure of the emissions on the ground showed that their spectral and temporal characteristics corresponded to emissions by the satellites in localized zones at different L-shells. The ground-based observations at lower frequencies correlated with the satellite observations at larger L-shells. In order to localize the regions of the generation of the VLF emissions observed at Kannuslehto auroral station at different frequencies, we calculated the ray trajectories of waves from the equator for the plasma density distributions detected by Van Allen Probes. The calculations of the trajectories showed that the VLF waves detected at Kannuslehto station could travel to the ground only if they propagated in the large-scale density ducts (700–900 km) observed by Van Allen Probes.  相似文献   

3.
The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.  相似文献   

4.
Global electron content (GEC) as a new ionospheric parameter was first proposed by Afraimovich et al. [2006]. GEC is equal to the total number of electrons in the near-Earth space. GEC better than local parameters reflects the global response to a change in solar activity. It has been indicated that, during solar cycle 23, the GEC dynamics followed similar variations in the solar UV irradiance and F 10.7 index, including the 11-year cycle and 27-day variations. The dynamics of the regional electron content (REC) has been considered for three belts: the equatorial belt and two midlatitude belts in the Northern and Southern hemispheres (±30° and 30°–65° geomagnetic latitudes, respectively). In contrast to GEC, the annual REC component is clearly defined for the northern and southern midlatitude belts; the REC amplitude is comparable with the amplitude of the seasonal variations in the Northern Hemisphere and exceeds this amplitude in the Southern Hemisphere by a factor of ～1.7. The dayside to nightside REC ratio, R(t), at the equator is a factor of 1.5 as low as such a GEC ratio, which indicates that the degree of nighttime ionization is higher, especially during the solar activity maximum. The pronounced annual cycle with the maximal R(t) value near 8.0 for the winter Southern Hemisphere and summer Northern Hemisphere is typical of midlatitudes.  相似文献   

5.
Ionospheric time delay (VΔt) variability using Global Positioning System (GPS) data over Akure (7.15°N, 5.12°E), Nigeria, has been studied. The observed variability of VΔt in comparison to older results of vertical total electron content (TEC) across similar regions has shown equivalent signatures. Higher monthly mean values of VΔt (MVΔt) were observed during daytime as compared to nighttime (pre- and post-midnight) hours in all months. The highest MVΔt observed in September during daytime hours range between ~6 and ~21 ns (~1.80 and ~6.30 m) and at post-midnight, they are in the range of ~1 to ~6 ns (~0.3 to ~1.80 m). The possible mechanisms responsible for this variability were discussed. Seasonal VΔt were investigated as well.  相似文献   

6.
The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N e and N i ) and temperatures (T e and T i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N e , N i , T e , and T i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.  相似文献   

7.
The study presents the results of the analysis of the F2-layer critical frequency variations obtained for the winter periods of 2008–2010, during which sudden stratospheric warmings were observed. The data were obtained at Kaliningrad ionospheric station (54.6° N, 20° E) with the Parus digital ionosonde in standard sounding mode. The mean daily foF2 values were used in the analysis. The results of spectral analysis based on continuous wavelet transform showed that, during all of the warmings that occurred in 2008–2010, the foF2 time variations demonstrated the presence of wave processes with periods of approximately 5?6 days, as well as more extended processes with periods of ~10?13 and 23?30 days. These periods coincide with the characteristic periods of planetary waves observed in the mesosphere during sudden stratospheric warmings, while the 13- and 30-day periods can be conditioned by the influence of the Sun.  相似文献   

8.
An alternative model for the nonlinear interaction term Snl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of Snl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of Snl, and to compare it to the existing models for Snl in WAVEWATCH; the DIA and WRT. Some differences between the different models for Snl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of Snl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.  相似文献   

9.
The modified scale M s(20R) is developed for the magnitude classification of the earthquakes of Russia’s Far East based on the surface wave amplitudes at regional distances. It extends the applicability of the classical Gutenberg scale M s(20) towards small epicentral distances (0.7°–20°). The magnitude is determined from the amplitude of the signal that is preliminarily bandpassed to extract the components with periods close to 20 s. The amplitude is measured either for the surface waves or, at fairly short distances of 0.7°–3°, for the inseparable wave group of the surface and shear waves. The main difference of the M s(20R) scale with the traditional M s(BB) Soloviev–Vanek scale is its firm spectral anchoring. This approach practically eliminated the problem of the significant (up to–0.5) regional and station anomalies characteristic of the M s(BB) scale in the conditions of the Far East. The absence of significant station and regional anomalies, as well as the strict spectral anchoring, make the M s(20R) scale advantageous when used for prompt decision making in tsunami warnings for the coasts of Russia’s Far East.  相似文献   

10.
The results of analysis of variations in the sporadic layer critical frequency (foEs) for winter periods of 2008–2010 in which sudden stratospheric warmings were observed are presented in the paper. The data were obtained at Kaliningrad ionospheric station (54.6° N, 20° E) by a Parus digital ionosonde under the usual sounding regime with an interval of 15 min. Daily mean values of foEs were used for the analysis. Solar and geomagnetic activity remained low during the periods under study, making it possible to relate the quasiwave time variations in foEs to the parameters of stratospheric warmings. The results of spectral analysis performed on the basis of continuous wavelet transform showed that, during all warmings occurring in 2008–2010, time variations in foEs show the presence of wave processes with a period of an order of 5 days and longer ones with a period of ~10—11 days. These periods coincide with characteristic periods of planetary waves observed in the atmosphere during sudden stratospheric warnings.  相似文献   

11.
We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h Am , which corresponds to the maximum increase, and layer peak height h m F, while h Am is always less than h m F. The difference between h Am and h m F linearly increases with increasing h m F. Whereas the difference is ~38 km for h m F = 280 km, it is ~54 km for h m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.  相似文献   

12.
Simultaneous observations of high-latitude long-period irregular pulsations at frequencies of 2.0–6.0 mHz (ipcl) and magnetic field disturbances in the solar wind plasma at low geomagnetic activity (Kp ~ 0) have been studied. The 1-s data on the magnetic field registration at Godhavn (GDH) high-latitude observatory and the 1-min data on the solar wind plasma and IMF parameters for 2011–2013 were used in an analysis. Ipcl (irregular pulsations continuous, long), which were observed against a background of the IMF Bz reorientation from northward to southward, have been analyzed. In this case other solar wind plasma and IMF parameters, such as velocity V, density n, solar wind dynamic pressure P = ρV2 (ρ is plasma density), and strength magnitude B, were relatively stable. The effect of the IMF Bz variation rate on the ipcl spectral composition and intensity has been studied. It was established that the ipcl spectral density reaches its maximum (~10–20 min) after IMF Bz sign reversal in a predominant number of cases. It was detected that the ipcl average frequency (f) is linearly related to the IMF Bz variation rate (ΔBz/Δt). It was shown that the dependence of f on ΔBz/Δt is controlled by the α = arctan(By/Bx) angle value responsible for the MHD discontinuity type at the front boundary of magnetosphere. The results made it possible to assume that the formation of the observed ipcl spectrum, which is related to the IMF Bz reorientation, is caused by solar wind plasma turbulence, which promotes the development of current sheet instability and surface wave amplification at the magnetopause.  相似文献   

13.
The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M0?=?1.12×1017 Nm, M w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω?2 source model) by attenuation parameters Q(f)?=?500f0.9, κ?=?0.04 s, and fmax?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.  相似文献   

14.
Properties of seismoelectric waves in relation to natural earthquakes have been investigated. The electromagnetic disturbances were analyzed to test the hypothesis that pulse-like electric variations are directly related to microcracks as source. Because variation is very difficult to detect, there have been few quantitative field investigations. We used selected events with clear S and P phases from the data catalog obtained before the Tohoku earthquake in 2011. The electric strength of the fast P wave (Pf), S wave (S), and electromagnetic wave (EM) associated with formation of cracks of tensile mode were estimated. The co-seismic electric signal accompanied by the S wave has the largest strength, well above the noise level, and the EM wave has the lowest strength. Analytical estimation of the ratio of the strengths of the Pf and EM phases to that of the S phase by use of Pride’s equations gave results partially in agreement with observation (the order was Apf > As > Aem). The strength of the observed electromagnetic mode is approximately two orders of magnitude larger than that estimated from the theory. We suggest this greater strength can be attributed to the converted modes at layer contracts or to the effect of the boundary between free atmosphere and crust. Overall agreement between observations and theoretical estimates suggests that electromagnetic anomalies, crustal deformation, and groundwater changes can be investigated on the basis of the unified equations for the coupled electromagnetics, acoustics, and hydrodynamics of porous media.  相似文献   

15.
Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s???1 (for Rayleigh waves) and 8 ×10???3 s???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.  相似文献   

16.
Lower Cretaceous C-isotope records show intermittent negative/positive spikes, and consistent patterns of coeval chemostratigraphic curves thus document shifts that signal simultaneous responses of temporal changes in the global carbon reservoir. The standard pattern registered by the δ 13Corg and δ 13Ccarb in Lower Aptian sediments includes distinct isotope segments C1 to C8 (Menegatti et al., 1998). In the El Pui section, Organyà Basin, Spain, C-isotope segment C2 is the longest interval preceding segments C3–C6 associated with oceanic anoxic event 1a (OAE 1a), and reveals a distinct negative shift of ~1.8‰ to ~2.23‰ defining the C-isotope pattern within that interval. Total inorganic carbon (TIC), total organic carbon (TOC), δ 13Corg, microfacies, n-alkanes show no difference before, during, or after the negative inflection. The biomarkers indicate that organic matter (OM) mainly originates from algal/microbial sources because short-chain length homologues (≤nC19) dominate. nC20 through nC25 indicate some contribution from aquatic vegetation, but little from higher plants (>nC25), as also suggested by the terrestrial/aquatic ratio of n-alkanes or (TAR) = [(nC27+nC29+nC31)/(nC15+nC17+nC19)] (averages 0.085). We suggest that conjoint pulses of contemporaneous LIPs (Ontong Java) and massive explosive volcanism in northeast Asia, the Songliao Basin (SB-V), best conform to plausible causes of the negative intra-C2 carbon isotopic excursion (CIE) at that time. Because of its apparent common occurrence the intra-C2 inflection could be a useful marker harbinger to the more pronounced CIE C3, the hallmark of OAE1a.  相似文献   

17.
Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.  相似文献   

18.
The occurrence probabilities of the first and second anomalous nighttime local maximums in the diurnal variations in the electron density at a maximum of the ionospheric F 2 layer (NmF2) in the region where the crest (hump) of the equatorial anomaly originates in the northern geographic hemisphere have been studied using the data of the stations for vertical sounding of the ionosphere (Paramaribo, Dakar, Quagadougou, Ahmedabad, Delhi, Calcutta, Chongoing, Guangzhou, Taipei, Chung-Li, Okinawa, Yamagawa, Panama, and Bogota) from 1957 to 2004. It has been demonstrated that the anomalous nighttime NmF2 maximums are least frequently formed at ～53° geomagnetic longitude. The calculations have indicated that the studied probabilities are independent of solar activity. Geomagnetic activity weakly affects the rate of occurrence of the first nighttime NmF2 maximum at geomagnetic longitudes of approximately 140° to 358°. At geomagnetic longitudes of approximately 16° to 70° (i.e., in the longitudinal zone of a decreased occurrence frequency of anomalous nighttime maximums), the occurrence probability of the first anomalous nighttime NmF2 maximum under geomagnetically quiet conditions is pronouncedly lower than under geomagnetically disturbed conditions. The dependence of the occurrence probabilities of the first and second anomalous nighttime NmF2 maximums on the month number in a year has been studied.  相似文献   

19.
Simultaneous morning Pc5 pulsations (f ~ 3–5 mHz) in the geomagnetic field, aurora intensities (in the 557.7 and 630.0 nm oxygen emissions and the 471.0 nm nitrogen emission), and riometer absorption, were studied based on the CARISMA, CANMOS, and NORSTAR network data for the event of January 1, 2000. According to the GOES-8 satellite observations, these Pc5 geomagnetic pulsations are observed as incompressible Alfvén waves with toroidal polarization in the magnetosphere. Although the Pc5 pulsation frequencies in auroras, the geomagnetic field, and riometer absorption are close to one another, stable phase relationships are not observed between them. Far from all trains of geomagnetic Pc5 pulsations are accompanied by corresponding auroral pulsations; consequently, geomagnetic pulsations are primary with respect to auroral pulsations. Both geomagnetic and auroral pulsations propagate poleward, and the frequency decreases with increasing geomagnetic latitude. When auroral Pc5 pulsations appear, the ratio of the 557.7/630.0 nm emission intensity sharply increases, which indicates that auroral pulsations result from not simply modulated particle precipitation but also an additional periodic acceleration of auroral electrons by the wave field. A high correlation is not observed between Pc5 pulsations in auroras and the riometer absorption, which indicates that these pulsations have a common source but different generation mechanisms. Auroral luminosity modulation is supposedly related to the interaction between Alfvén waves and the region with the field-aligned potential drop above the auroral ionosphere, and riometer absorption modulation is caused by the scattering of energetic electrons by VLF noise pulsations.  相似文献   

20.
The characteristics and interplanetary excitation conditions of isolated bursts of Pi2 geomagnetic pulsations observed during the development of magnetospheric substorms (substorm Pi2) and in its absence (nonsubstorm Pi2) on the night side of the Earth are comparatively analyzed. It is shown that, regardless of the local time and season, the amplitude of isolated Pi2 substorm bursts is always higher than that of the nonsubstorm ones, and the periods and duration of the wave packets of substorm Pi2 bursts are less than those of nonsubstorms. Diurnal and seasonal variations in the characteristics of the two groups of Pi2 bursts differ in the form and position of maxima and minima. It is found that the start of excitation of isolated Pi2 bursts, during substorms and in its absence, is controlled by the preferred direction of the interplanetary magnetic field (IMF) vector perpendicular to the Sun–Earth line (angle θxB = arccos(Bx/B) → 90°). It is assumed that isolated Pi2 bursts of both groups are triggered by reorientation of the IMF vector in the ecliptic plane and the plane perpendicular to it ~15 min before their onset. The most likely source of midlatitude isolated Pi2 bursts during substorm development and in its absence are bursty bulk flows (BBFs) in the plasma sheet of the magnetospheric tail, the regularities of which coincide in many respects with the observed features of Pi2 bursts.  相似文献   

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

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