首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A radial velocity study is presented of the cataclysmic variable V378 Pegasi (PG 2337 + 300). It is found to have an orbital period of 0.13858 ± 0.00004 d (3.32592 ± 0.00096 h). Its spectrum and long-term light curve suggest that V378 Peg is a nova-like variable, with no outbursts. We use the approximate distance and position in the Galaxy of V378 Peg to estimate E(B − V) = 0.095, and use near-infrared magnitudes to calculate a distance of 680 ± 90 pc and MV = 4.68 ± 0.70, consistent with V378 Peg being a nova-like. Time-resolved photometry taken between 2001 and 2009 reveals a period of 0.1346 ± 0.0004 d (3.23 ± 0.01 h). We identify this photometric variability to be negative superhumps, from a precessing, tilted accretion disk. Our repeated measurements of the photometric period of V378 Peg are consistent with this period having been stable between 2001 and 2009, with its negative superhumps showing coherence over as many as hundreds or even thousands of cycles.  相似文献   

2.
We present a time-resolved differential photometric study and time series analysis of the nova-like cataclysmic variable star LQ Peg. We discover three periodicities in the photometry, one with a period of 3.42 ± 0.03 h, and another with a period of 56.8 ± 0.01 h. We interpret these to be the apsidal superhump and precessional periods of the accretion disk, respectively, and predict that the orbital period of LQ Peg is 3.22 ± 0.03 h. The third periodicity, with a period of 41.3 ± 0.01 h, we interpret to be the nodal precessional period of the accretion disk. We also report a flare that lasted four minutes and had an energy in visible light of (1.2 ± 0.3) × 1036 ergs, or 104-5 times more energetic than the largest solar flares, comparable to the most energetic visible-light stellar flares known. We calculate the absolute magnitude of LQ Peg to be MJ = 4.78 ± 0.54, and its distance to be 800 ± 200 pc.  相似文献   

3.
Here we show results from thermal-infrared observations of km-sized binary near-Earth asteroids (NEAs). We combine previously published thermal properties for NEAs with newly derived values for three binary NEAs. The η value derived from the near-Earth asteroid thermal model (NEATM) for each object is then used to estimate an average thermal inertia for the population of binary NEAs and compared against similar estimates for the population of non-binaries. We find that these objects have, in general, surface temperatures cooler than the average values for non-binary NEAs as suggested by elevated η values. We discuss how this may be evidence of higher-than-average surface thermal inertia. This latter physical parameter is a sensitive indicator of the presence or absence of regolith: bodies covered with fine regolith, such as the Earth’s moon, have low thermal inertia, whereas a surface with little or no regolith displays high thermal inertia. Our results are suggestive of a binary formation mechanism capable of altering surface properties, possibly removing regolith: an obvious candidate is the YORP effect.We present also newly determined sizes and geometric visible albedos derived from thermal-infrared observations of three binary NEAs: (5381) Sekhmet, (153591) 2001 SN263, and (164121) 2003 YT1. The diameters of these asteroids are 1.41 ± 0.21 km, 1.56 ± 0.31 km, and 2.63 ± 0.40 km, respectively. Their albedos are 0.23 ± 0.13, 0.24 ± 0.16, and 0.048 ± 0.015, respectively.  相似文献   

4.
《Astroparticle Physics》2012,35(7):435-448
MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50 − 60 GeV. Advanced stereo analysis techniques allow MAGIC to achieve a sensitivity as good as (0.76 ± 0.03)% of the Crab Nebula flux in 50 h of observations above 290 GeV. The angular resolution at those energies is better than ∼0.07°. We also perform a detailed study of possible systematic effects which may influence the analysis of the data taken with the MAGIC telescopes.  相似文献   

5.
We present values from the Cassini Visual and Infrared Mapping Spectrometer (VIMS) of four fundamental disk-integrated spectrophotometric properties (bolometric Bond albedo, solar phase curve, phase integral, and geometric albedo at 7-15 different wavelengths in the λ = 0.35-5.1 μm range) for five mid-sized saturnian icy satellites: Rhea, Dione, Tethys, Mimas, and Enceladus. These values, which include data from the period 2004-2008 and add to past VIMS phase curves, include opposition surge effects at down to fractions of a degree in solar phase angle for several moons and extend to over double the solar phase angle coverage of the Voyager mission. We also present new rotational light curves for Rhea and Dione at 7 near-infrared bands not previously available in ground-based or spacecraft studies. The bolometric Bond albedos we derive are as follows: 0.48 ± 0.09 (Rhea), 0.52 ± 0.08 (Dione), 0.61 ± 0.09 (Tethys), 0.67 ± 0.10 (Mimas), and 0.85 ± 0.11 (Enceladus). We also provide breakdowns of the major photometric quantities in both leading and trailing hemispheres. These refined parameters can be used to construct future bolometric Bond albedo maps that will contribute to surface composition identification studies, as well as models of volatile transport and sublimation. Through such applications, these data will help to determine the physical properties of surface particles, how the E-ring affects the inner saturnian moons, what is responsible for the dark albedo patterns seen on Tethys, and if these moons (e.g., Dione) are geologically active.  相似文献   

6.
The atmospheres of Mars and Titan are loaded with aerosols that impact remote sensing observations of their surface. Here we present the algorithm and the first applications of a radiative transfer model in spherical geometry designed for planetary data analysis. We first describe a fast Monte-Carlo code that takes advantage of symmetries and geometric redundancies. We then apply this model to observations of the surface of Mars and Titan at the terminator as acquired by OMEGA/Mars Express and VIMS/Cassini. These observations are used to probe the vertical distribution of aerosols down to the surface. On Mars, we find the scale height of dust particles to vary between 6 km and 12 km depending on season. Temporal variations in the vertical size distribution of aerosols are also highlighted. On Titan, an aerosols scale height of 80 ± 10 km is inferred, and the total optical depth is found to decrease with wavelength as a power-law with an exponent of −2.0 ± 0.4 from a value of 2.3 ± 0.5 at 1.08 μm. Once the aerosols properties have been constrained, the model is used to retrieve surface reflectance properties at high solar zenith angles and just after sunset.  相似文献   

7.
We have carried out deep (V ∼ 21 mag) UBVRI photometric study of the star cluster Stock 18. These along with archival Infrared data have been used to derive the basic cluster parameters and also to study the star formation processes in and around the cluster region. The distance to the cluster is derived as 2.8 ± 0.2 kpc while its age is estimated as 6.0 ± 2.0 Myr. Present study indicates that interstellar reddening is normal in the direction of the cluster. The mass function slope is found to be −1.37 ± 0.27 for the mass range 1 < M/M < 11.9. There is no evidence found for the effect of mass segregation in main-sequence stars of the cluster. A young stellar population with age between 1-2 Myr have been found in and around the cluster region. The presence of IRAS and AKARI sources with MSX intensity map also show the youth of the Sh2-170 region.  相似文献   

8.
We performed a complete wavelet analysis of Saturn’s C ring on 62 stellar occultation profiles. These profiles were obtained by Cassini’s Ultraviolet Imaging Spectrograph High Speed Photometer. We used a WWZ wavelet power transform to analyze them. With a co-adding process, we found evidence of 40 wavelike structures, 18 of which are reported here for the first time. Seventeen of these appear to be propagating waves (wavelength changing systematically with distance from Saturn). The longest new wavetrain in the C ring is a 52-km-long wave in a plateau at 86,397 km. We produced a complete map of resonances with external satellites and possible structures rotating with Saturn’s rotation period up to the eighth order, allowing us to associate a previously observed wave with the Atlas 2:1 inner Lindblad resonance (ILR) and newly detected waves with the Mimas 6:2 ILR and the Pandora 4:2 ILR. We derived surface mass densities and mass extinction coefficients, finding σ = 0.22(±0.03) g cm−2 for the Atlas 2:1 ILR, σ = 1.31(±0.20) g cm−2 for the Mimas 6:2 ILR, and σ = 1.42(±0.21) g cm−2 for the Pandora 4:2 ILR. We determined a range of mass extinction coefficients (κ = τ/σ) for the waves associated with resonances with κ = 0.13 (±0.03) to 0.28(±0.06) cm2 g−1, where τ is the optical depth. These values are higher than the reported values for the A ring (0.01-0.02 cm2 g−1) and the Cassini Division (0.07-0.12 cm2 g−1 from Colwell et al. (Colwell, J.E., Cooney, J.H., Esposito, L.W., Srem?evi?, M. [2009]. Icarus 200, 574-580)). We also note that the mass extinction coefficient is probably not constant across the C ring (in contrast to the A ring and the Cassini Division): it is systematically higher in the plateaus than elsewhere, suggesting smaller particles in the plateaus. We present the results of our analysis of these waves in the C ring and estimate the mass of the C ring to be between3.7(±0.9) × 1016 kg and 7.9(±2.0) × 1016 kg (equivalent to an icy satellite of radius between 28.0(±2.3) km and 36.2(±3.0) km with a density of 400 kg m−3, close to that of Pan or Atlas). Using the ring viscosity derived from the wave damping length, we also estimate the vertical thickness of the C ring between 1.9(±0.4) m and 5.6(±1.4) m, comparable to the vertical thickness of the Cassini Division.  相似文献   

9.
Joseph Masiero 《Icarus》2010,207(2):795-799
We have investigated the effect of rotation on the polarization of scattered light for the near-Earth asteroid (1943) Anteros using the Dual Beam Imaging Polarimeter on the University of Hawaii’s 2.2 m telescope. Anteros is an L-type asteroid that has not been previously observed polarimetrically. We find weak but significant variations in the polarization of Anteros as a function of rotation, indicating albedo changes across the surface. Specifically, we find that Anteros has a background albedo of pv = 0.18 ± 0.02 with a dark spot of pv < 0.09 covering <2% of the surface.  相似文献   

10.
Jeremy Bailey  Linda Ahlsved 《Icarus》2011,213(1):218-232
We have obtained spatially resolved spectra of Titan in the near-infrared J, H and K bands at a resolving power of ∼5000 using the near-infrared integral field spectrometer (NIFS) on the Gemini North 8 m telescope. Using recent data from the Cassini/Huygens mission on the atmospheric composition and surface and aerosol properties, we develop a multiple-scattering radiative transfer model for the Titan atmosphere. The Titan spectrum at these wavelengths is dominated by absorption due to methane with a series of strong absorption band systems separated by window regions where the surface of Titan can be seen. We use a line-by-line approach to derive the methane absorption coefficients. The methane spectrum is only accurately represented in standard line lists down to ∼2.1 μm. However, by making use of recent laboratory data and modeling of the methane spectrum we are able to construct a new line list that can be used down to 1.3 μm. The new line list allows us to generate spectra that are a good match to the observations at all wavelengths longer than 1.3 μm and allow us to model regions, such as the 1.55 μm window that could not be studied usefully with previous line lists such as HITRAN 2008. We point out the importance of the far-wing line shape of strong methane lines in determining the shape of the methane windows. Line shapes with Lorentzian, and sub-Lorentzian regions are needed to match the shape of the windows, but different shape parameters are needed for the 1.55 μm and 2 μm windows. After the methane lines are modeled our observations are sensitive to additional absorptions, and we use the data in the 1.55 μm region to determine a D/H ratio of 1.77 ± 0.20 × 10−4, and a CO mixing ratio of 50 ± 11 ppmv. In the 2 μm window we detect absorption features that can be identified with the ν5 + 3ν6 and 2ν3 + 2ν6 bands of CH3D.  相似文献   

11.
The first of a new class of objects now known as main belt comets (MBCs) or “activated asteroids” was identified in 1996. The seven known members of this class have orbital characteristics of main belt asteroids yet exhibit dust ejection like comets. In order to constrain their physical and orbital properties we searched the Thousand Asteroid Light Curve Survey (TALCS; Masiero, J.R., Jedicke, R., Durech, J., Gwyn, S., Denneau, L., Larsen, J. [2009]. Icarus 204, 145-171) for additional candidates using two diagnostics: tail and coma detection. This was the most sensitive MBC survey effort to date, extending the search from MBCs with H ∼ 18 (D ∼ 1 km) to MBCs as small as H ∼ 21 (D ∼ 150 m).We fit each of the 924 objects detected by TALCS to a PSF model incorporating both a coma and nuclear component to measure the fractional contribution of the coma to the total surface brightness. We determined the significance of the coma detection using the same algorithm on a sample of null detections of comparable magnitude and rate of motion. We did not identify any MBC candidates with this technique to a sensitivity limit on the order of cometary mass loss rate of about 0.1 kg/s.Our tail detection algorithm relied on identifying statistically significant flux in a segmented annulus around the candidate object. We show that the technique can detect tail activity throughout the asteroid belt to the level of the currently known MBCs. Although we did not identify any MBC candidates with this technique, we find a statistically significant detection of faint activity in the entire ensemble of TALCS asteroids. This suggests that many main belt asteroids are active at very low levels.Our null detection of MBCs allows us to set 90% upper confidence limits on the number distribution of MBCs as a function of absolute magnitude, semi-major axis, eccentricity, and inclination. There are ?400,000 MBCs in the main belt brighter than HV = 21 (∼150-m in diameter) and the MBC:MBA ratio is ?1:400.We further comment on the ability of observations to meaningfully constrain the snow line’s location. Under some reasonable and simple assumptions we claim 85% confidence that the contemporary snow line lies beyond 2.5 AU.  相似文献   

12.
We report on the detection of H2 as seen in our analysis of twilight observations of the lunar atmosphere observed by the LAMP instrument aboard NASA’s Lunar Reconnaissance Orbiter. Using a large amount of data collected on the lunar atmosphere between September 2009 and March 2013, we have detected and identified, the presence of H2 in the native lunar atmosphere, for the first time. We derive a surface density for H2 of 1.2 ± 0.4 × 103 cm−3 at 120 K. This is about 10 times smaller than originally predicted, and several times smaller than previous upper limits from the Apollo era data.  相似文献   

13.
We study the kinematics of the Galactic thin and thick disk populations using stars from the RAVE survey’s second data release together with distance estimates from Breddels et al. (2010). The velocity distribution exhibits the expected moving groups present in the solar neighborhood. We separate thick and thin disk stars by applying the X (stellar-population) criterion of Schuster et al. (1993), which takes into account both kinematic and metallicity information. For 1906 thin disk and 110 thick disk stars classified in this way, we find a vertical velocity dispersion, mean rotational velocity and mean orbital eccentricity of (σW, 〈VΦ〉, 〈e〉)thin = (18 ± 0.3 km s−1, 223 ± 0.4 km s−1, 0.07 ± 0.07) and (σW, 〈VΦ〉, 〈e〉)thick = (35 ± 2 km s−1, 163 ± 3 km s−1, 0.31 ± 0.16), respectively. From the radial Jeans equation, we derive a thick disk scale length in the range 1.5-2.2 kpc, whose greatest uncertainty lies in the adopted form of the underlying potential. The shape of the orbital eccentricity distribution indicates that the thick disk stars in our sample most likely formed in situ with minor gas-rich mergers and/or radial migration being the most likely cause for their orbits. We further obtain mean metal abundances of 〈[M/H]〉thin = +0.03 ± 0.17, and 〈[M/H]〉thick = −0.51 ± 0.23, in good agreement with previous estimates. We estimate a radial metallicity gradient in the thin disk of −0.07 dex kpc−1, which is larger than predicted by chemical evolution models where the disk grows inside-out from infalling gas. It is, however, consistent with models where significant migration of stars shapes the chemical signature of the disk, implying that radial migration might play at least part of a role in the thick disk’s formation.  相似文献   

14.
Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm−1 spectral resolution, adequate to resolve the lines at high pressure (pN2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν2+2ν3 band (7515-7620 cm−1), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm−1 sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH4 mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.  相似文献   

15.
Using a Markov chain model, we consider the regolith growth on a small body in orbit around Saturn, subject to meteoritic bombardment, and assuming all impact ejecta are re-collected. We calculate the growth of regolith and the fractional pollution, assuming an initial pure ice body and amorphous carbon as a pollutant. We extend the meteorite flux of Cuzzi and Estrada (Cuzzi, J., Estrada, P. [1998]. Icarus 132, 1-35) to larger sizes to consider the effect of disruption of the moonlet on other moonlets in the ensemble. This is a relatively small effect, completely negligible for moonlets of 1 m radius. For the given impact model, fractional pollution reaches 22% for 1 m bodies, but only 3% for 10 m bodies, 1.7% for 20 m bodies, and 1% for 30 m bodies after 4 byr. By considering an ensemble of moonlets, which have identical cross-sections for releasing and capturing ejecta, this analysis can be extended to a model of particles in Saturn’s rings, where the calculated spectra can be compared to observed ring spectra. The measured spectral reflectance of Saturn’s rings from Cassini observations therefore constrains the size and age of the ring particles. The comparison between 1 m, 10 m, 20 m, and 30 m particles confirms that for larger ring mass, the current rings would be less polluted; for the largest particles, we expect negligible changes in the UV spectrum after 4 byr of meteoritic bombardment. We consider two end members for mixing of the meteoritic material: areal and intimate. Given the uncertainties in the actual mixing of the meteoritic infall and in its composition (as a worst case, we assume the meteoritic material is 100% amorphous carbon, intimately mixed) initially pure ice 30 m ring particles would darken after 4 byr of exposure by 15%.  相似文献   

16.
In this paper, we present standard Johnson UBV photometry of the eclipsing binary BD+36 3317 which is known as a member of Delta Lyrae (Stephenson 1) cluster. We determined colors and brightness of the system, calculated E(B − V) color excess. We discovered that the system shows total eclipse in secondary minimum. Using this advantage, we found that the primary component of the system has B9 − A0 spectral type. Although there is no published orbital solution, we tried to estimate the physical properties of the system from simultaneous analysis of UBV light curves with 2003 version of Wilson-Devinney code. Then we considered photometric solution results together with evolutionary models and estimated the masses of the components as M1 = 2.5 M and M2 = 1.6 M. Those estimations gave the distance of the system as 353 pc. Considering the uncertainties in distance estimation, resulting distance is in agreement with the distance of Delta Lyrae cluster.  相似文献   

17.
Anderson and Schubert [2007. Saturn's Gravitational field, internal rotation, and interior structure. Science 317, 1384-1387 (paper I)] proposed that Saturn's rotation period can be ascertained by minimizing the dynamic heights of the 100 mbar isosurface with respect to the geoid; they derived a rotation period of 10 h 32 m 35 s. We investigate the same approach for Jupiter to see if the Jovian rotation period is predicted by minimizing the dynamical heights of its isobaric (1 bar pressure level) surface using zonal wind data. A rotation period of 9 h 54 m 29.7 s is found. Further, we investigate the minimization method by fitting Pioneer and Voyager occultation radii for both Jupiter and Saturn. Rotation periods of 9 h 55 m 30 s and 10 h 32 m 35 s are found to minimize the dynamical heights for Jupiter and Saturn, respectively. Though there is no dynamical principle requiring the minimization of the dynamical heights of an isobaric surface, the successful application of the method to Jupiter lends support to its relevance for Saturn.We derive Jupiter and Saturn rotation periods using equilibrium theory to explain the difference between equatorial and polar radii. Rotation periods of 9 h 55 m 20 s and 10 h 31 m 49 s are found for Jupiter and Saturn, respectively. We show that both Jupiter's and Saturn's shapes can be derived using solid-body rotation, suggesting that zonal winds have a minor effect on the planetary shape for both planets.The agreement in the values of Saturn's rotation period predicted by the different approaches supports the conclusion that the planet's period of rotation is about 10 h 32 m.  相似文献   

18.
Using a global numerical model, we have studied how the present Martian magnetosphere may have looked in the past when the planet had a global intrinsic magnetic field. A Mars version (HYB-Mars) of the self-consistent quasi-neutral hybrid model was used which treats the ions as particles and the electrons as a massless charge-neutralizing fluid. We compare four cases where an intrinsic dipole magnetic field was 0 nT (the present situation), 10, 30, and 60 nT at the surface of Mars along the magnetic equator. We find that the 10 nT dipolar magnetic field already results in a magnetosphere which in many respects is more Earth-like than, a non-magnetized, “induced” magnetosphere. However, the 10 nT dipole magnetosphere is still relatively strongly connected to the interplanetary magnetic field, while the 30 nT dipole case, and especially the 60 nT dipole case, results in a magnetosphere whose morphology is determined predominantly by the Martian intrinsic magnetic field. A change of the magnetosphere due to a decreasing dipole magnetic field strength from 60 to 0 nT could have happened during the history of Mars when a globally magnetized Mars turned into the present, globally non-magnetized, planet.  相似文献   

19.
Scattering and absorption of sunlight by aerosols are integral to understanding the radiative balance of any planetary atmosphere covered in a haze, such as Titan and possibly the early Earth. One key optical parameter of an aerosol is its refractive index. We have simulated both Titan and early Earth organic haze aerosols in the laboratory and measured the real and imaginary portion of their refractive index at λ = 532 nm using cavity ringdown aerosol extinction spectroscopy. This novel technique allows analysis on freely-floating particles minutes after formation. For our Titan analog particles, we find a real refractive index of n = 1.35 ± 0.01 and an imaginary refractive index k = 0.023 ± 0.007, and for the early Earth analog particles we find n = 1.81 ± 0.02 and k = 0.055 ± 0.020. The Titan analog refractive index has a smaller real and similar imaginary refractive index compared to most previous laboratory measurements of Titan analog films, including values from Khare et al. (Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callcott, T.A., Williams, M.W. [1984]. Icarus 60, 127-137). These newly measured Titan analog values have implications for spacecraft retrievals of aerosol properties on Titan. The early Earth analog has a significantly higher real and imaginary refractive index than Titan analogs reported in the literature. These differences suggest that, for a given amount of aerosol, the early Earth analog would act as a stronger anti-greenhouse agent than the Titan analog.  相似文献   

20.
Arecibo (2380 MHz, 12.6 cm) and Goldstone (8560 MHz, 3.5 cm) delay-Doppler radar images obtained in July and August of 2000 reveal that 4486 Mithra is an irregular, significantly bifurcated object, with a central valley ∼380 m deep and a long axis potentially exceeding 2 km. With its bimodal appearance, Mithra is a strong candidate for a contact binary asteroid. Sequences of Goldstone images spanning up to 3 h per day show very little rotation and establish that Mithra is an unusually slow rotator. We used Goldstone and Arecibo data to estimate Mithra’s 3D shape and spin state. We obtain prograde (λ = 337°, β = 19°) and retrograde (λ = 154°, β = −19°) models that give comparable fits, have very similar shapes roughly resembling an hourglass, and have a rotation period of 67.5 ± 6.0 h. The dimensions of these two models are very similar; for the prograde solution the maximum dimensions are X = 2.35 ± 0.15 km, Y = 1.65 ± 0.10 km, Z = 1.44 ± 0.10 km. Dynamical analysis of our models suggests that in the past, Mithra most likely went through a period of even slower rotation with its obliquity close to 90°. The spin rate is predicted to be increasing due to thermal torque (YORP), while the obliquity, which is currently +68° and +106° for the prograde and retrograde models, respectively, is predicted to move away from 90°.  相似文献   

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

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