首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
The results of a statistical investigation of bright galaxies with a UV excess are given. The data obtained for them are compared with those for normal galaxies. A list of bright galaxies (mpg ≤14m.5 ) with a UV excess has been compiled, numbering 461. Their absolute stellar magnitudes and average geometrical linear diameters were determined. The relationships between Mpg and log D are given for galaxies with a UV excess and normal galaxies. The graphs of the relationships are expressed by second-degree equations. It is shown that the UV evolutionary stage of the galaxies has no influence on the relationship between Mpg and log D. Translated from Astrofizika, Vol. 43, No. 1, pp. 21-32, January–March, 2000.  相似文献   

2.
High resolution electric field and particle data, obtained by the S23L1 rocket crossing over a discrete prebreakup arc in January 1979, are studied in coordination with ground observations (Scandinavian Magnetometer Array—SMA, TV and all-sky cameras) in order to clarify the electrodynamics of the arc and its surroundings. Height-integrated conductivities have been calculated from the particle data, including the ionization effects of precipitating protons and assuming a steady state balance between ion production and recombination losses. High resolution optical information of arc location relative to the rocket permitted a check of the validity of this assumption for each flux tube passed by the rocket. Another check was provided by a comparison between calculated (equilibrium values) and observed electron densities along the rocket trajectory. A way to compensate for the finite precipitation time when calculating the electron densities is outlined. The height-integrated HalI-Pedersen conductivity ratio is typically 1.4 within the arc and about 1 at the arc edges, indicative of a relatively softer energy spectrum there. The height-integrated conductivities combined with the DC electric field measurements permitted calculation of the horizontal ionospheric current vectors (J⊥), Birkeland currents (from div J⊥) and energy dissipation through Joule heating (ΣpE2). An eastward current of typically 1 A m?1 was found to be concentrated mainly to the arc region and equatorward of it. A comparison has been made with the equivalent current system deduced from ground based magnetometer data (SMA) showing a generally good agreement with the rocket results. An intense Pedersen current peak (1.2 A m?1) was found at the southern arc edge. This edge constituted a division line between a very intense (> 10 μA m?1) and localized (~ 6 km) downward current sheet to the south, probably carried by upward flowing cold ionospheric electrons and a more extended upward current sheet (> 10 μA m?2) over the arc carried by measured precipitating electrons. Joule and particle heating across the arc were anticorrelated, consistent with the findings of Evans et al. (1977) with a total value of about 100mW m?2.  相似文献   

3.
An analytic model has been developed for toroidal quarter wave (QW) oscillations in the Earth’s magnetosphere using idealistic and highly asymmetric ionospheric boundary condition. The background magnetic field is dipolar and plasma density distribution is governed by a power law 1/r m where r is the geocentric distance of any point along the field line and m is the density index. The solution thus obtained has been compared with the numerical solutions. Earlier workers had developed the analytic model for trivial 1/r 6 (m=6) type of plasma density distribution along the field line for which the period of the fundamental is twice that of corresponding half wave (i.e. toroidal oscillation in the symmetric ionospheric boundary). The present analytic model does reproduce this feature. In addition, it is seen that this ratio decreases for lower values of m. Moreover, for a particular value of m, this ratio shows a decreasing trend with increased harmonic number. The spatial characteristics of QW obtained from present analytic model are in excellent agreement with those computed numerically, thereby validating the model. The departure of frequency computed analytically from that obtained numerically is significant only for the fundamental and this departure reduces sharply with the increased harmonic number. It should be noted that there is no such departure for 1/r 6 type of plasma density distribution and spatial structures as well as frequency computed from the present analytic model match perfectly with those computed numerically.  相似文献   

4.
The main objective of this paper is to provide researchers that investigate fossil phytolith assemblages and model/data comparisons a new tool for estimating C3/C4 grass composition over time. We tested the reliability of modern soil phytolith assemblages and phytolith indices for tracing the dominance of different grass subfamilies and tree cover density. We analyzed modern soil phytolith assemblages from sites over elevation gradients on Mount Kenya (Kenya), Mount Rungwe and around Lake Masoko (southern Tanzania). These data were compared with available botanical data. A phytolith index named Ic, proved to be an effective proxy of the proportions of Pooideae, Arundinoideae and Bambusoideae grasses (mainly C3 grasses) versus Panicoideae grasses (mainly C4 grasses), increasing with elevation in East-Africa. When tropical mountains are covered by open habitats (e.g. grasses and shrublands), Ic should be a reliable proxy of the C3/C4 grass composition. These results highlight the value of the phytolith index Ic, when interpreting paleo-environmental records from tropical mountains, to: 1) better understand past local and regional C3/C4 grass distributions and associated climatic changes and 2) increase the set of C3/C4 data available for model/data comparisons.  相似文献   

5.
As part of a multiproxy investigation, phytoliths were extracted from sediments in a 197-m core in Hawkes Bay, New Zealand. They provide a continuous vegetation–climate record spanning the time period from at least the last interglacial (marine oxygen isotope stage 5) to the present. The phytolith record demonstrates that grass/cyperaceae grew during warmer periods and woody taxa dominated the site during colder periods.During the present interglacial, the Poukawa basin is occupied by a shallow lake surrounded by an extensive fen. During colder–drier periods, the floor of the basin dried out and woody taxa occupied the basin floor. This contrasts with the pollen record, which demonstrates a converse pattern. The apparent discrepancy reflects the purely local provenance of the phytolith assemblage.Significant changes in phytolith assemblages occur at the same depth as major tephras, indicating a sharp decline in trees and shrubs and a surge in grass and cyperaceae. A series of successional changes follow each major tephra fall. Initially, the woody taxa are killed off and replaced by grass and cyperaceae that rapidly colonise the fresh surface. Trees and or shrubs succeed the grass and cyperaceae after a significant lag.  相似文献   

6.
Abstract— Video observations of the Leonid shower aboard two aircraft in the 1998 Leonid multi-instrument aircraft campaign and from ground locations in China are presented. Observing at altitude proved particularly effective, with four times higher rates due to low extinction and low angular velocity at the horizon. The rates, derived from a total of 2500 Leonid meteors, trace at least two distinct dust components. One dominated the night of 1998 November 16/17. This two-day wide component was rich in bright meteors with r = N (m + 1)/N (m) ≈ 1.5 (s = 1.4) and peaked at an influx of 3.1 ± 0.4 × 10?12 m?2 s?1 (for particles of mass <7 × 10?5 g) at solar longitude Λ0 ≈ 234.52 (Eq. J2000). The other more narrow component peaked on 1998 November 17/18 at Λ0 = 235.31 ± 0.01. Rates were elevated above the broad component between Λ0 = 235.15 and 235.40, symmetric around the current node of the parent comet 55P/Tempel-Tuttle, peaking at 5.1 ± 0.2 × 10?12 m?2 s?1. The population index was higher, r = 1.8 ± 0.1 (s = 1.7), but not as high as in past Leonid storms (r = 3.0). The flux profile of this component has an unusual asymmetric shape, which implies a blend of contributions from at least two different but relatively recent epochs of ejection. The variation of r across the profile might be due to mass-dependent ejection velocities of the narrowest component. High rates of faint meteors occurred only in an isolated five-minute interval at Λ0 = 235.198, which is likely the result of a single meteoroid breakup in space.  相似文献   

7.
T. Takakura 《Solar physics》1984,91(2):311-324
In some gradual hard X-ray bursts with high intensity, hard X-ray source (15–40 keV) is steadily located in the corona along with softer X-ray source (5–10 keV).Two stationary models, high density and high temperature models, are proposed to solve the difficult problem of confinement of hot (or nonthermal) plasma in the direction of the magnetic field along the loops in the corona. In both models, an essential point is that the effective X-ray source is composed of fine dense filamentary loops imbeded in a larger rarefied coronal loop, and the electron number density in the filaments is so high as 1011–1012 cm-3. If the density is so high heat conduction can be as reasonably small as of the order of 1027 erg s -1 for the given emission measures of observed X-rays, since the required cross-sectional area is small and also classical conduction is valid. Collisional confinement of thermal tail, and nonthermal electrons if any, up to 50–60 keV in the filaments is also possible, so that the hard X-ray images can be loop like structure instead of double source (foot points).High density model is applicable to the coronal filamentary loops with temperature T m < 5 × 107 K at the loop summit. The heat flow from the summit downwards is lost almost completely by the radiation from the loop during the conduction to the foot points. A continuous energy release is assumed near the summit to maintain the stationary temperature T m, and pressure balance is maintained along the loop. In this model, the number density at the summit is given by n m - 106 T m 2 /sm, where s m is the length of the loop from the summit to the foot point, and the distribution of temperature and density along the loop are given by T = T m(s/sm)1/3 and n = n m(s/sm)-1/3, respectively.High temperature model is applicable to the filamentary loops with higher temperature up to about 108.5 K and comparatively lower number density as 1011 cm-3 for the requirement of magnetic confinement of the hot plasma in radial direction. The radiation from the loop is negligibly small in this model so that the heat flux is nearly conserved down to the foot points. In this case, temperature gradient is smaller than that of the high density model, depending on the tapering of the magnetic bottle.In both models, the differential emission measure is maximum at the highest temperature T m and the brightness distribution along the loop shows a maximum around the summit of the loop if some magnetic tapering is taken into account.  相似文献   

8.
Abstract— The grain size distribution and shapes of lechatelierite inclusions (silica glass inclusions) have been determined from 20 splash-form tektites from the Khorat Plateau, northeastern Thailand. The chemical and petrographic properties are reviewed, and the absence of any type of inclusion other than bubbles is confirmed. These data suggest that the parent material for the lechatelierite inclusions is not the conventional detrital quartz. One possible precursor is silica of plant origin in the form of biogenic opal-CT. According to this model, the lechatelierite inclusions are formed by shock melting of opal phytoliths in plants. These opal phytolith melts were included in the shock-melted soil and bedrock, jetted from the impact site. The expansion of the vapor plume ejected the melt droplets in ballistic trajectories. This model is extended to all tektite groups, because of the similarity between lechatelierite inclusions in them.  相似文献   

9.
David Andrew Fisher 《Icarus》2007,187(2):430-441
A time varying stable isotope model for the D/H history of Mars water cycle is developed with variable atmosphere, space loss rate, ground and ice cap flux rates. It considers coupled ground reservoirs and traces D/H in the air and reservoirs secularly and over obliquity cycles. The various flux rates are prescribed time variables that simulate surface flux, and solar driven space loss rates. Predicted bulk averages for the ice cap, ground ice reservoirs and atmosphere span the observed ranges reported by Mumma et al. [Mumma, M.J., Novak, R.E., DiSanti, M.A., Bonev, B., Dello Russo, N., Magee-Sauer, K., 2003. The Martian Atmosphere. Conference Reports of “Sixth International Conference on Mars Atmosphere,” No. 3186]. When the dominant obliquity cycle variations are scaled so that the model delivers present seasonal variations, the present long term bulk D/H average for the ice cap is ∼+2.7 (equivalent to +1700‰ in δ(D) wrt SMOW). The obliquity driven D/H cycle in the ice cap's layers varies between 3 and 6. The smaller more accessible reservoirs have larger bulk averages with the smallest being able to reach D/H values over 9 within ∼105 years. Small hypothetical solar activity driven variations in the escape rate to space and in the fractionation constant [Krasnopolsky, V.A., Feldman, P.D., 2001. Science 294, 1914-1917] for the escape process can produce a “solar wiggle” whose D/H amplitude can reach 0.1 (δ(D) amplitude of 100‰). Because of the temporal variability, a single modern measured atmospheric D/H ratio at a particular Ls cannot tell very much about the total water inventory of Mars. A bulk average for the Northern Ice Cap and better still a dated vertical profile of D/H from the ice cap would, however, go a long way towards illuminating the “modern” water history of Mars. The age and stability of the Northern Ice Cap and the D/H history locked in the layering is discussed. An ice cap that is very young and exchanges its mass through the atmosphere often will necessarily have a large D/H.  相似文献   

10.
In this paper, we present 50 surface water samples collected during the IMAGES III cruise (June–July 1997) along a transect from New Zealand to the China Sea (42°S–178°E, 21°N–120°E) covering a temperature range from 13.3 to 30.4 °C. A very worthwhile aspect of this study is a coupling of both biomarker (alkenone) and coccolithophorid counting. We show that the U37k′–temperature relationship is very similar to the Prahl et al. [Geochim. Cosmochim. Acta 52 (1988) 2203] culture calibration and to the global core top calibration of Müller et al. [Geochim. Cosmochim. Acta 62 (1998) 1757]. However, in the warmest surface waters of the Western Pacific ocean (>26.4 °C) where Gephyrocapsa oceanica is likely the most widespread species, the associated U37k′ has a constant value of 1.0. The consequence is that above this temperature threshold, U37k′ cannot be used as an accurate paleothermometer.  相似文献   

11.
The HRSC (image 0037) and MOC imagery and MOLA altimetry were used to determine the following parameters of the lava flows typical of the southern slope of the Martian volcano Olympus: the length (13–35 km), the width (0.2–4.8 km), and the angles of ground slopes along which these flows advanced (3.4°–6.9°). To measure the thickness of the flows, we applied a method which had never been used before for Mars. In this method, the apparent thickness obtained from the MOC images and the slope steepness obtained from the MOLA data are used to determine the true thickness. The average estimates of the thickness of lava flows vary from 4 to 11 m and from 4 to 26 m for the volcano flanks and caldera scarps, respectively. These values are close to those of terrestrial basalt flows and to the lower limit found for the Martian flows by other researchers. Based on the performed measurements, we estimated the lava yield strength (0.9 × 103?3.6 × 104 Pa), the supply rate (24–137 m3/s), and the viscosity (1.4 × 103?2.8 × 107 kg/m s). These values are close to the estimates found for the Martian lavas by other researchers and to the characteristic values of these parameters for terrestrial lava flows with basalt and basalt-andesite composition.  相似文献   

12.
A highly sensitive device has been developed for the detection of micrometeoroids with masses extending down to and below the expected classical solar radiation pressure limit. This limit occurs at a mass of about 10−16 kg for a particle density of 3 × 103 kg m−3. The detector operates on the principle that charge (Q) is released when a hypervelocity microparticle impacts on a solid surface. This process has previously been investigated at particle velocities (ν) greater than 1 km sec−1 and the empirical relationship Q = Kmανβ obtained where m is the particle mass and K, α and β are constants. In the present study the validity of this relationship has been demonstrated for iron particles with masses from 10−16 to 10−13 kg in the velocity range 0.05–1.4 km sec−1 impacting on a molybdenum target. A value of α ≈ 1 is indicated and a value of β of 3.2 ± 0.4 has been obtained.  相似文献   

13.
Air and ground temperatures measured in Eastern Siberia has been compiled and analyzed. The analysis of mean annual air temperatures measured at 52 meteorological stations within and near the East-Siberian transect during the period from 1956 through 1990 demonstrates a significant and statistically significant (at 0.05 level) positive trend ranging from 0.065 to 0.59 °C/10 yr. A statistically significant (at 0.05 level) positive trend was also observed in mean annual ground temperatures for the same period. The permafrost temperature reflects changes in air temperature on a decadal time scale much better than on an interannual time scale. Generally, positive trends in mean annual ground temperatures are slightly smaller in comparison with trends in mean annual air temperatures, except for several sites where the discordance between the air and ground temperatures can be explained by the winter snow dynamics. The average trend for the entire region was 0.26 °C/10 yr for ground temperatures at 1.6 m depth and 0.29 °C/10 yr for the air temperatures. The most significant trends in mean annual air and ground temperatures were in the southern part of the transect, between 55° and 65° N. Numerical modeling of ground temperatures has been performed for Yakutsk and Tiksi for the last 70 yr. Comparing the results of these calculations with a similar time series obtained for Fairbanks and Barrow in Alaska shows that similar variations of ground temperatures took place at the same time periods in Yakutsk and Fairbanks, and in Tiksi and Barrow. The decadal and longer time scale fluctuations in permafrost temperatures were pronounced in both regions. The magnitudes of these fluctuations were on the order of a few degrees centigrade. The fluctuations of mean annual ground temperatures were coordinated in Fairbanks and Yakutsk, and in Barrow and Tiksi. However, the magnitude and timing of these fluctuations were slightly different for each of the sites.  相似文献   

14.
A study of the structure and kinematics of the Galaxy from Tautenburg Schmidt plates taken towards the Galactic centre (l = 17.0°, b = +0.8°) is presented. Proper motions and B, V magnitudes were determined for about 36 500 stars up to the limiting magnitude V = 16m.8 in a field of 8.95 square degrees. Proper motion accuracy of about 3 mas/year has been obtained for stars brighter than V = 15m. The rms errors of stellar magnitudes and (B–V) colours is about 0.1 mag. The majority of field stars in the survey are main sequence stars and red giants of the disk. They belong to the Local or Sagittarius-Carina arms, or they are located between these spiral arms. Comparing the modelled and observed distributions of magnitudes and colours, we specified the interstellar extinction determined in the preceding study of open clusters in this field. The luminosity function towards the Galactic centre was determined for stars with absolute magnitudes from -4m.35 to +9m. Kinematical and spatial distribution parameters up to 4 kpc from the Sun were obtained as a function of galactocentric distance.  相似文献   

15.
In 2003–2008, highly accurate photoelectric and CCD observations of the close binary system DI Her were performed in the V band. The light curves of three primary and three secondary eclipses were constructed. These observations, along with the highly accurate photoelectric observations of other authors obtained in different years from 1963 to 1986, have confirmed the difference between the observed (1 . o 3 ± 0 . o 1/100 yr) and theoretical (4 . o 3/100 yr) rates of apsidal motion. Our photometric data are indicative of a possible variability in the system with period P′ = 1.175 days and amplitude A′ = 0 . m 011, which is probably related to the pulsations of one of the components. There may be a third body in the system that produces in-phase variations in the times of primary and secondary minima with a period of 10.5 yr and an amplitude of 1 . m 5.  相似文献   

16.
A novel methodology for evaluating the field of anisotropically scattered radiation within a homogeneous slab atmosphere of arbitrary optical thickness is provided. It departs from the traditional radiative transfer approach in first considering that the atmosphere is illuminated by an isotropic light source. From the solution of this problem, it subsequently proceeds to that for the more conventional case of monodirectional illumination. The azimuthal dependence of the field is separated in the usual manner by an harmonic expansion, leaving a problem in four dimensions (=optical depth, 0=thickness, , =directions of incidence and scattering) which, as is well known, is numerically extremely inconvenient. Two auxiliary radiative transfer formulations of increasing dimensionality are considered: (i) a transfer equation for the newly introduced functionb m(,,0) with Sobolev's function m(,0) playing the role of a source-function. Because the incident direction does not intervene, m is simply expressed as a single integral term involvingb m. For bottom illumination, an analogous equation holds for the other new functionh m(,,0). However, simple reciprocity relations link the two functions so that it is only necessary to considerb m; (ii) a transfer equation for the other new functiona m(,,,0) with a source-function provided by Sobolev's functionD m(,,0). For bottom illumination, another functionf m(,,,0) is introduced; by a similar argument using reciprocity relations,f m is reduced toa m rendering necessary only the consideration ofa m. However, a fundamental decomposition formula is obtained which shows thata m is expressible algebraically in terms of functions of a single angular variable. The functions m andD m are shown to be the values in the horizontal plane ofb m anda m, respectively. The other auxiliary functionsX m andY m are also expressed algebraically in terms ofb m. These results enable one to proceed to the final step of evaluating the radiation field for monodirectional illumination. The above reductions toalgebraic relations involving only the functionb m appear to be more advantageous than Sobolev's (1972) recent approach; they also circumvent some basic numerical difficulties in it. We believe the present approach may likewise prove to be superior to most (if not all) other methods of solution known heretofore.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory under Contract No. NAS-7-100 sponsored by the National Aeronautics and Space Administration.  相似文献   

17.
A new prediction technique based on logarithmic values is proposed to predict the maximum amplitude (R m) of a solar cycle from the preceding minimum aa geomagnetic index (aa min). The correlation between lnR m and lnaa min (r=0.92) is slightly stronger than that between R m and aa min (r=0.90). From this method, cycle 24 is predicted to have a peak size of R m(24)=81.7(1±13.2%). If the suggested error in aa (3 nT) before 1957 is corrected, the correlation coefficient between R m and aa min (r=0.94) will be slightly higher, and the peak of cycle 24 is predicted much lower, R m(24)=52.5±13.1. Therefore, the prediction of R m based on the relationship between R m and aa min depends greatly on the accurate measurement of aa.  相似文献   

18.
Du  Zhanle 《Solar physics》2011,270(1):407-416
The correlation coefficient (r) between the maximum amplitude (R m) of a sunspot cycle and the preceding minimum aa geomagnetic index (aa min), in terms of geomagnetic cycle, can be fitted by a sinusoidal function with a four-cycle periodicity superimposed on a declining trend. The prediction index (χ) of the prediction error relative to its estimated uncertainty based on a geomagnetic precursor method can be fitted by a sinusoidal function with a four-and-half-cycle periodicity. A revised prediction relationship is found between the two quantities: χ<1.2 if r varies in a rising trend, and χ>1.2 if r varies in a declining trend. The prediction accuracy of R m depends on the long-term variation in the correlation. These results indicate that the prediction for the next cycle inferred from this method, R m(24)=87±23 regarding the 75% level of confidence (1.2-σ), is likely to fail. When using another predictor of sunspot area instead of the geomagnetic index, similar results can be also obtained. Dynamo models will have better predictive powers when having considered the long-term periodicities.  相似文献   

19.
We revisit the appropriate energies to be used for computing heat production from 26Al decay. Due to the complexity of the decay scheme of this radioisotope, previous geophysical studies have used values ranging from 1.2 to 4 MeV per decay. The upper bound corresponds to the difference in mass energy between the 26Al and 26Mg ground states. This includes energy carried away by neutrinos, which does not contribute to heating planetary material. The lower bound does not account for the heating caused by the absorption of the γ rays from the excited 26Mg, or for the annihilation energy deposited in the material if the decay occurs inside even small planetesimals. Based on the calculations described by Schramm et al. [Schramm, D., Tera, F., Wasserburg, G.J., 1970. The isotopic abundance of 26Mg and limits on 26Al in the early Solar System. Earth Planet. Sci. Lett. 10, 44-59] updated with the most recent nuclear constants, we recommend using a heat production value of 3.12 MeV per decay, which is the total energy of disintegration minus the energy carried off by the neutrinos. This heat production value is higher than the value used in the modeling of Iapetus by Castillo-Rogez et al. [Castillo-Rogez, J., Matson, D.L., Sotin, C., Johnson, T.V., Lunine, J.I., Thomas, P.C., 2007. Iapetus’ geophysics: Rotation rate, shape, and equatorial ridge. Icarus 190, 179-202] by about a factor 2.5. The resulting estimate of the time of formation of Iapetus is shifted by about 1 Myr, to between ∼3.4 and 5.4 Myr after the production of the calcium-aluminum inclusions (CAIs).  相似文献   

20.
Systematic and uniform sets of photometric and polarimetric observations of comet 1P/Halley have been analyzed. The phase dependence of brightness for comet Halley was obtained at phase angles α ranging from 1.4° ≤ α ≤ 65°. The following parameters were determined: the amplitude of the opposition effect Δm = 0.75m ± 0.06m; the half-width at a half-maximum of intensity HWHM = 6.4° ± 1.6°; the linear phase coefficient β = 0.0045 ± 0.0001 mag/deg for α from 30° ≤ α ≤ 65°; and the phase angle at which a nonlinear increase in brightness starts, α opp ≈ 31°. For the first time, the phase-angle dependence was obtained for the color of the dust of comet Halley: the color index BC-RC systematically decreases with increasing phase angle. Such a phase behavior of the dust color can be caused by the decrease in the mean size of dust particles that occurs when the comet approaches the Sun. For comet Halley, the negative polarization branch is almost symmetric; the minimal value of polarization is P min = −1.54% at a phase angle αmin = 10.5°, and the inversion angle is αinv = 21.7°. A comparative analysis of the phase functions of brightness and polarization has been performed for the cometary dust and atmosphereless bodies. Among the latter are low-albedo asteroids of the P and C types (102 Miriam and 47 Aglaja, respectively), as well as Deimos; high-albedo objects, such as the E-type asteroid 64 Angelina and the icy satellite of Jupiter Europa; and the Moon with its intermediate albedo. The possibility of a weak depression in the negative polarization branch of comets Halley and 47P/Ashbrook-Jackson at phase angles smaller than 2° is discussed.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 353–363.Original Russian Text Copyright © 2005 by Rosenbush.  相似文献   

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

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