首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We propose that the Taurid meteor shower may contain bodies able to survive and be recovered as meteorites. We review the expected properties of meteorite‐producing fireballs, and suggest that end heights below 35 km and terminal speeds below 10 km s?1 are necessary conditions for fireballs expected to produce meteorites. Applying the meteoroid strength index (PE criteria) of Ceplecha and McCrosky (1976) to a suite of 33 photographically recorded Taurid fireballs, we find a large spread in the apparent meteoroid strengths within the stream, including some very strong meteoroids. We also examine in detail the flight behavior of a Taurid fireball (SOMN 101031) and show that it has the potential to be a (small) meteorite‐producing event. Similarly, photographic observations of a bright, potential Taurid fireball recorded in November of 1995 in Spain show that it also had meteorite‐producing characteristics, despite a very high entry velocity (33 km s?1). Finally, we note that the recent Maribo meteorite fall may have had a very high entry velocity (28 km s?1), further suggesting that survival of meteorites at Taurid‐like velocities is possible. Application of a numerical entry model also shows plausible survival of meteorites at Taurid‐like velocities, provided the initial meteoroids are fairly strong and large, both of which are characteristics found in the Taurid stream.  相似文献   

2.
We calculated equilibrium chemical composition of a mixture of meteoritic vapor and air during fireball events, i.e. during penetration of large meteoroids into terrestrial atmosphere. Different types of fireballs were considered, and calculations were performed for wide ranges of temperatures and pressures. Chemical composition at the quenching point was estimated by comparison of hydrodynamic and chemical reaction time scales. For the typical fireball temperatures of 4000-5000 K, most elements are expected to be in the form of atoms and ions. Notable exceptions are Si and C, which are expected to be mainly in the form of SiO and CO. Other molecules abundant at these temperatures are N2 and NO. Metal monoxides are most abundant at 2000-2500 K and are formed during the cooling phase. Conditions for formation of other molecules such as , CN, C2 and OH were also considered. The composition of freshly ablated meteoroid material was studied using the MAGMA code.  相似文献   

3.
Oxygen isotope measurements of olivine in main group (MG) pallasites by traditional laser fluorination method are associated with some uncertainties including terrestrial weathering, incomplete olivine reaction, and sample state. We improved our laser fluorination approach by pretreating olivine grains with acid to remove terrestrial weathering products and by modifying the sample holder for an efficient and complete laser reaction. Our experiments on Brahin olivine demonstrate that acid-washing successfully removes the terrestrial weathering with <0.1‰ variation in δ18O value and, at the same time, improving the ∆17O value significantly. We also achieved a complete olivine fluorination by employing a custom-designed sample holder with “V”-shaped profile having rounded bottom because incomplete/partial reaction of olivine gives comparatively lighter δ18O values. Using these new techniques, we present precise triple oxygen isotope data (N = 72) of 25 olivine samples separated from main group pallasites. The data are, on average, ~0.5‰ heavier in δ18O relative to the values published in the literature for the same samples. Critically, the ∆17O values of MG pallasites and to some extent their Fo-contents suggest that there are at least two populations of olivine. Based on our improved data set, we propose that MG pallasites potentially have high-∆17O- and low-∆17O-bearing subgroups that are statistically distinct. The subgroups present average ∆17O values of −0.166 ± 0.003 (2SE;= 16) and −0.220 ± 0.003 (2SE; N = 9), respectively. Furthermore, the high-∆17O-bearing subgroup samples trend toward lower Fo-contents compared to the other subgroup. Taken together, our data provide evidence that argues against a single parent body origin for MG pallasites.  相似文献   

4.
A progress report on cosmic fireballs is presented. The main new results are: (a) the phenomenon should be almost universal, and most explosive -ray sources should show the characteristic fireball spectrum; (b) even if the radiation density is insufficient, pair production in electron-proton or electron-electron scattering might start the fireball; (c) some computed fireball spectra are shown. They all have in common a 1/E low-energy behaviour, a 100 keV flattening, and a 0.5 MeV cut-off.Paper presented at the Symposium on Cosmic Gamma-Ray Bursts, held at Toulouse, France, 26–29 November, 1979.  相似文献   

5.
Abstract— Among the three large camera networks carrying out fireball observations through the seventies and eighties, the “European Fireball Network” is the last one still in operation. The network today consists of more than 34 all-sky and fish-eye cameras deployed with ~100 km spacing and covering an area of ~106 km2, in the Czech and Slovak Republics, Germany, as well as parts of Belgium, Switzerland, and Austria. Network operation results in ~10 000 image exposures per year, which represent on average 1200 h of clear sky observations—as imaging periods are restricted due to daylight, moonlight, and clouds. The cameras detect currently large meteors at a rate of ~50 per year; this is in good agreement with the encounter rates determined in previous fireball studies. From sightings of “meteorite candidates” (fireballs that may have deposited meteorites) and meteorite recoveries in the network area, we estimate that 15% of the influx of meteoritic matter is currently observed by the cameras, whereas <1% is recovered on the ground. Issues to be addressed by future fireball observations include the study of very large meteoroids (>1000 kg) for which statistics are currently very poor and an examination of their relationship to NEOs (near-Earth objects) identified by current NEO search programs.  相似文献   

6.
Three bright fireballs belonging to the August θ‐Aquillid (ATA) meteor shower were photographed by the Tajikistan fireball network in 2009. Two of them are classified as the meteorite‐dropping fireballs according to the determined parameters of the atmospheric trajectories, velocities, masses, and densities. Detection of the more dense bodies among cometary meteoroids points to a heterogeneous composition of the parent comet, and supports the suggestion that some meteorites might originate in the outer solar system, in the given case from the Jupiter‐family comet reservoir. A search for the stream's parent was undertaken among the near‐Earth asteroids (NEAs); as a result, the asteroid 2004MB6 was identified as a possible progenitor of the ATA meteoroid stream. Investigation of the orbital evolution of the 2004MB6 and the fireball‐producing meteoroid TN170809A showed that both objects have similar secular variations in the orbital elements during 7 kyr. The comet‐like orbit of the 2004MB6 and its association with the ATA shower suppose a cometary origin of the asteroid.  相似文献   

7.
Minor bodies that, at the present stage, have orbital characteristics similar to those of the Pribram meteorite and that differ from it in the date of activity of the radiants by no more than a 1.2-month period in either direction were detected on the basis of catalog data on the orbits of asteroids, individual fireballs, and fireball and meteor streams obtained from photographic observations. The following objects comprising the Prbram family are among such minor bodies: three asteroids (1863, 4486, and J98S70J) of the Apollo group, three fireball streams, five meteor streams, and ten fireballs. This extensive system of bodies consists of three branches: a Northern (N), an Ecliptical (Q), and a Southern (S) Branch. The family of meteor streams associated with the periodic comet Pons–Winnecke appears to be related to this family. Thus, there emerges an intricate complex of small bodies that is similar to the well-known Taurid complex. In the distribution of various populations of minor bodies according to the quantity c of the Tisserand criterion, both of these complexes of minor bodies, like the group of cometoids (or cometlike asteroids), are situated in the region of unstable motion or, to be more exact, near the gap that arises in the L 2 and L 3 collinear points of libration.  相似文献   

8.
The Hamburg (H4) meteorite fell on 17 January 2018 at 01:08 UT approximately 10 km north of Ann Arbor, Michigan. More than two dozen fragments totaling under 1 kg were recovered, primarily from frozen lake surfaces. The fireball initial velocity was 15.83 ± 0.05 km s?1, based on four independent records showing the fireball above 50 km altitude. The radiant had a zenith angle of 66.14 ± 0.29° and an azimuth of 121.56 ± 1.2°. The resulting low inclination (<1°) Apollo‐type orbit has a large aphelion distance and Tisserand value relative to Jupiter (Tj) of ~3. Two major flares dominate the energy deposition profile, centered at 24.1 and 21.7 km altitude, respectively, under dynamic pressures of 5–7 MPa. The Geostationary Lightning Mapper on the Geostationary Operational Environmental Satellite‐16 also detected the two main flares and their relative timing and peak flux agree with the video‐derived brightness profile. Our preferred total energy for the Hamburg fireball is 2–7 T TNT (8.4–28 × 109 J), which corresponds to a likely initial mass in the range of 60–225 kg or diameter between 0.3 and 0.5 m. Based on the model of Granvik et al. (2018), the meteorite originated in an escape route from the mid to outer asteroid belt. Hamburg is the 14th known H chondrite with an instrumentally derived preatmospheric orbit, half of which have small (<5°) inclinations making connection with (6) Hebe problematic. A definitive parent body consistent with all 14 known H chondrite orbits remains elusive.  相似文献   

9.
We report simultaneous multicolour observations in 5 bands (UBVRI) of the flickering variability of the cataclysmic variable AE Aqr. Our aim is to estimate the parameters (colours, temperature, size) of the fireballs that produce the optical flares. The observed rise times of the optical flares are in the interval 220‐440 s. We estimate the dereddened colours of the fireballs as (UB)0∼0.8‐1.4, (BV)0∼0.03‐0.24, and (VI)0∼0.26‐0.78. We find for the fireballs temperatures of 10000‐25000 K, masses of (7‐90)x1019 g, and sizes of (3‐7)x109 cm (using a distance of d = 86 pc). These values refer to the peak of the flares observed in the UBVRI bands. The data are available upon request from the authors (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

10.
Since 2006, systematic double-station photographic observations of fireballs using all-sky cameras equipped with Zeiss Distagon “fisheye” objectives (f/3.5, f = 30 mm) with a 180° field of view have been carried out at two observatories, Gissar (GisAO) and Sanglokh (IAOS), of the Institute of Astrophysics of the Tajik Academy of Sciences. In the method of astrometric reduction of fireball photographs, the empirical formulae for converting the measured coordinates to horizontal celestial coordinates are used. These formulae contain 12 unknown constants to be determined by the least-squares method and the iteration method. Such an approach enables the determination of the coordinates of an object at any point of the celestial hemisphere with a precision close to the theoretical limit whose value is quite comparable with the measurement errors. In the photometric reduction, the dependence of the measured width of the diurnal star trails on their magnitudes was used. As a result of astrometric and photometric reduction of the double-station photographs of five fireballs, the data on atmospheric trajectories, the coordinates of radiants, orbits in interplanetary space, light curves, and photometric masses of meteoroids which produced fireballs were obtained, and the belonging of fireballs to the known meteor showers was determined as well.  相似文献   

11.
In this paper we discuss two methods, one analytical and the other graphical, to determine the trajectory of a fireball using the arrival times of atmospheric shock waves recorded by a seismic network. In the analytical method the trajectory and the raypaths are assumed to be straight and we solve for the fireball velocity, the azimuth (?) and elevation angle (δ) of the trajectory, the coordinates of the intersection of the trajectory with the earth's surface, and the corresponding intersection time (t0). Because the problem is nonlinear, we solve it iteratively. The fireball velocity cannot be determined uniquely, and trades off with t0. The graphical method is based on the drawing of contours of arrival times, which should be elliptical for fireball shock waves. If the distribution of seismic stations is appropriate, the horizontal projection of the fireball is given by the axis of symmetry of the contours, which allows the estimation of ?, while δ can be estimated from the spacing between contours along the symmetry axis. Application of the two methods to data from four fireballs shows that the graphically derived parameters can be within a few degrees of the analytical parameters. In addition, a fireball recorded in the Czech Republic has reliable trajectory parameters derived from video recordings, which allows an independent assessment of the quality of the parameters determined analytically. In particular, ? and δ have errors of and , respectively, which are not particularly large considering that the station distribution was not favorable.  相似文献   

12.
Abstract— Precise atmospheric trajectories including dynamic and photometric data on thirteen of the brightest Leonid fireballs have been determined from the double‐station photographic observations of Leonid meteors during the ground‐based expedition to China in 1998 November. the expedition was organized as a collaboration between the dutch and chinese academy of sciences and was supported by the leonid multi‐instrument aircraft campaign (mac) program (jenniskens and butow, 1999). All data presented here were taken at Xinglong Observatory and at a remote station, Lin Ting Kou near Beijing, on the night of 1998 November 16/17. At the Xinglong station, photographic cameras were accompanied by an all‐sky television camera equipped with an image intensifier and 15 mm fish‐eye objective in order to obtain precise timings for all observed meteors up to magnitude +2. Whereas beginning heights of photographed meteors are all lower than 130 km, those observed by the all‐sky television system are at ~160 km, and for three brightest events, even > 180 km. Such high beginnings for meteors have never before been observed. We also obtained a precise dynamic single‐body solution for the Leonid meteor 98003, including the ablation coefficient, which is an important material and structural quantity (0.16 s2 km?2). From this and from known photometry, we derived a density of this meteoroid of 0.7 g/cm3. Also, all PE coefficients indicate that these Leonid meteors belonged to the fireball group IIIB, which is typical for the most fragile and weak interplanetary bodies. From a photometric study of the meteor lightcurves, we found two typical shapes of light curves for these Leonid meteors.  相似文献   

13.
During the 5-year-observational phase of Gaia several one-dimensional observations will be obtained from the same objects at different epochs and transit angles. Combining the information from all the transits to produce a single two-dimensional image is scientifically important as this will allow to detect much fainter stars in the vicinity of the main source than otherwise possible from single observations on-board the satellite or from individual transmitted patches. Most of the emerging objects will be physical companions, but some background/foreground objects will also be found. The combination procedure further increases the image resolution in the across scan-direction. The procedure consists of combining patches that are stretched back to 2-D windows, only possible since the centroid positions of the windows are known. The initial image window from which the patches originate is not square, but elongated along the scan-direction if the Astro AF11 field of view is used in the current calculations. Transit angles can take any value between 0 and 360. The window combining method is a simple mapping procedure that is fast and should handle all realistic complexities occurring during the Gaia data reduction. The method is not intended to be a perfect image restoration technique, but its purpose is to fulfill the scientific requirement of detecting faint neighboring sources existing in the close vicinity of the main targets. In this study we show examples of the combined windows and illustrate how faint companions around primary targets could be detected. We investigate the performances and limitations of the method and discuss possible complexities that might occur during the final data analysis.  相似文献   

14.
A dm-radio emission with fiber bursts observed on 11 July 2005 was analyzed using wavelet filtration and spectral methods. In filtered radio spectra we found structures with different characteristic period P and frequency drift FD: i) fiber substructures (composed of dot emissions) with P 1≈ 0.5 s, FD1=− 87 MHz s−1 on average, ii) fiber structures with P 2≈1.9 s, and iii) drifting structures with P 3≈81.4 s, FD2=− 8.7, + 98.5, and − 21.8 MHz s−1. In the wavelet spectra we recognized patterns having the form of tadpoles. They were detected with the same characteristic periods P as found for the filtered structures. The frequency drift of the tadpole heads is found to be equal to the frequency drift of some groups of fibers for the long-period wavelet tadpoles (P 3) and to the frequency drift of individual fibers for the short-period tadpoles (P 2). Considering these wavelet tadpoles as signatures of propagating magnetoacoustic wave trains, the results indicate the presence of several wave trains in the fibers’ source. While the long-period wave trains trigger or modulate a whole group of fibers, the short-period ones look like being connected with individual fiber bursts. This result supports the model of fibers based on magnetoacoustic waves. Using a density model of the solar atmosphere we derived the velocities of the magnetoacoustic waves, 107 and 562 km s−1, and setting them equal to the Alfvén ones we estimated the magnetic field in the source of fiber bursts as 10.7 and 47.8 G.  相似文献   

15.
The coronal magnetic field above a particular photospheric region will vanish at a certain number of points, called null points. These points can be found directly in a potential field extrapolation or their density can be estimated from the Fourier spectrum of the magnetogram. The spectral estimate, in which the extrapolated field is assumed to be random and homogeneous with Gaussian statistics, is found here to be relatively accurate for quiet Sun magnetograms from SOHO’s MDI. The majority of null points occur at low altitudes, and their distribution is dictated by high wavenumbers in the Fourier spectrum. This portion of the spectrum is affected by Poisson noise, and as many as five-sixths of null points identified from a direct extrapolation can be attributed to noise. The null distribution above 1500 km is found to depend on wavelengths that are reliably measured by MDI in either its low-resolution or high-resolution mode. After correcting the spectrum to remove white noise and compensate for the modulation transfer function we find that a potential field extrapolation contains, on average, one magnetic null point, with altitude greater than 1.5 Mm, above every 322 Mm2 patch of quiet Sun. Analysis of 562 quiet Sun magnetograms spanning the two latest solar minima shows that the null point density is relatively constant with roughly 10% day-to-day variation. At heights above 1.5 Mm, the null point density decreases approximately as the inverse cube of height. The photospheric field in the quiet Sun is well approximated as that from discrete elements with mean flux 〈|φ|〉=1.0×1019 Mx distributed randomly with density n=0.007 Mm−2.  相似文献   

16.
High-resolution (∼0.22 Å) spectra of the north jovian aurora were obtained in the 905-1180 Å window with the Far Ultraviolet Spectroscopic Explorer (FUSE) on October 28, 2000. The FUSE instrument resolves the rotational structure of the H2 spectra and the spectral range allows the study of self-absorption. Below 1100 Å, transitions connecting to the v?2 levels of the H2 ground state are partially or totally absorbed by the overlying H2 molecules. The FUSE spectra provide information on the overlying H2 column and on the vibrational distribution of H2. Transitions from high-energy H2 Rydberg states and treatment of self-absorption are considered in our synthetic spectral generator. We show comparisons between synthetic and observed spectra in the 920-970, 1030-1080, and 1090-1180 Å spectral windows. In a first approach (single-layer model ), the synthetic spectra are generated in a thin emitting layer and the emerging photons are absorbed by a layer located above the source. It is found that the parameters of the single-layer model best fitting the three spectral windows are 850, 800, and 800 K respectively for the H2 gas temperature and 1.3×1018, 1.5×1020, and 1.3×1020 cm−2 for the H2 self-absorbing vertical column respectively. Comparison between the H2 column and a 1-D atmospheric model indicates that the short-wavelength FUV auroral emission originates from just above the homopause. This is confirmed by the high H2 rovibrational temperatures, close to those deduced from spectral analyses of H+3 auroral emission. In a second approach, the synthetic spectral generator is coupled with a vertically distributed energy degradation model, where the only input is the energy distribution of incoming electrons (multi-layer model ). The model that best fits globally the three FUSE spectra is a sum of Maxwellian functions, with characteristic energies ranging from 1 to 100 keV, giving rise to an emission peak located at 5 μbar, that is ∼100 km below the methane homopause. This multi-layer model is also applied to a re-analysis of the Hopkins Ultraviolet Telescope (HUT) auroral spectrum and accounts for the H2 self-absorption as well as the methane absorption. It is found that no additional discrete soft electron precipitation is necessary to fit either the FUSE or the HUT observations.  相似文献   

17.
Abstract– On April 9, 2009, at 3:00 CEST, a very bright fireball appeared over Carinthia and the Karavanke Mountains. The meteoroid entered the atmosphere at a very steep angle and disintegrated into a large number of objects. Two main objects were seen as separate fireballs up to an altitude of approximately 5 km, and witnesses reported loud explosions. Three stones were found with a total weight of approximately 3.611 kg. The measured activity of short‐lived cosmogenic radionuclides clearly indicates that two specimens result from a very recent meteorite fall. All cosmogenic radionuclide concentrations suggest a rather small preatmospheric radius of <20 cm; a nominal cosmic‐ray exposure age based on 21Ne is approximately 4 Ma, but the noble gas and radionuclide results in combination indicate a complex irradiation. Jesenice is a highly recrystallized rock with only a few relic chondrules visible in hand specimen and thin section. The texture, the large grain size of plagioclase, and the homogeneous compositions of olivines and pyroxenes clearly indicate that Jesenice is a L6 chondrite. The bulk composition of Jesenice is very close to the published average element concentration for L ordinary chondrites. The chondrite is weakly shocked (S3) as indicated by the undulatory extinction in olivine and plagioclase and the presence of planar fractures in olivine. Being weakly shocked and with gas retention ages of >1.7 Ga (4He) and approximately 4.3 Ga (40Ar), Jesenice seems not to have been strongly affected by the catastrophic disruption of the L‐chondrite parent body approximately 500 Ma ago.  相似文献   

18.
We describe the fall of the Dingle Dell (L/LL 5) meteorite near Morawa in Western Australia on October 31, 2016. The fireball was observed by six observatories of the Desert Fireball Network (DFN), a continental-scale facility optimized to recover meteorites and calculate their pre-entry orbits. The 30 cm meteoroid entered at 15.44 km s−1, followed a moderately steep trajectory of 51° to the horizon from 81 km down to 19 km altitude, where the luminous flight ended at a speed of 3.2 km s−1. Deceleration data indicated one large fragment had made it to the ground. The four person search team recovered a 1.15 kg meteorite within 130 m of the predicted fall line, after 8 h of searching, 6 days after the fall. Dingle Dell is the fourth meteorite recovered by the DFN in Australia, but the first before any rain had contaminated the sample. By numerical integration over 1 Ma, we show that Dingle Dell was most likely ejected from the Main Belt by the 3:1 mean motion resonance with Jupiter, with only a marginal chance that it came from the ν6 resonance. This makes the connection of Dingle Dell to the Flora family (currently thought to be the origin of LL chondrites) unlikely.  相似文献   

19.
The rapid neutron capture process (r-process) is one of the major nucleosynthesis processes responsible for the synthesis of heavy nuclei beyond iron. Isotopes beyond Fe are most exclusively formed in neutron capture processes and more heavier ones are produced by the r-process. Approximately half of the heavy elements with mass number A>70 and all of the actinides in the solar system are believed to have been produced in the r-process. We have studied the r-process in supernovae for production of heavy elements beyond A=40 with the newest mass values available. The supernovae envelopes at a temperature >109 K and neutron density of 1024 cm−3 are considered to be one of the most potential sites for the r-process. We investigate the r-process in a site-independent, classical approach which assumes a chemical equilibrium between neutron captures and photodisintegrations followed by a β-flow equilibrium. We have studied the r-process path corresponding to temperatures ranging from 1.0×109 K to 3.0×109 K and neutron density ranging from 1020 cm−3 to 1030 cm−3. The primary goal of the r-process calculations is to fit the global abundance curve for solar system r-process isotopes by varying time dependent parameters such as temperature and neutron density. This method aims at comparing the calculated abundances of the stable isotopes with observation. The abundances obtained are compared with supernova explosion condition and found in good agreement. The elements obtained along the r-process path are compared with the observed data at all the above temperature and density range.  相似文献   

20.
Abstract— We present instrumental observations of the Tagish Lake fireball and interpret the observed characteristics in the context of two different models of ablation. From these models we estimate the pre‐atmospheric mass of the Tagish Lake meteoroid to be ?56 tonnes and its porosity to be between 37 and 58%, with the lowest part of this range most probable. These models further suggest that some 1300 kg of gram‐sized or larger Tagish Lake material survived ablation to reach the Earth's surface, representing an ablation loss of 97% for the fireball. Satellite recordings of the Tagish Lake fireball indicate that 1.1 times 1012 J of optical energy were emitted by the fireball during the last 4 s of its flight. The fraction of the total kinetic energy converted to light in the satellite pass band is found to be 16%. Infrasonic observations of the airwave associated with the fireball establish a total energy for the event of 1.66 ± 0.70 kT TNT equivalent energy. The fraction of this total energy converted to acoustic signal energy is found to be between 0.10 and 0.23%. Examination of the seismic recordings of the airwave from Tagish Lake have established that the acoustic energy near the sub‐terminal point is converted to seismic body waves in the upper‐most portion of the Earth's crust. The acoustic energy to seismic energy coupling efficiency is found to be near 10?6 for the Tagish Lake fireball. The resulting energy estimate is near 1.7 kT, corresponding to a meteoroid 4 m in diameter. The seismic record indicates extensive, nearly continuous fragmentation of the body over the height intervals from 50 to 32 km. Seismic and infrasound energy estimates are in close agreement with the pre‐atmospheric mass of 56 tonnes established from the modeling. The observed flight characteristics of the Tagish Lake fireball indicate that the bulk compressive strength of the pre‐atmospheric Tagish Lake meteoroid was near 0.25 MPa, while the material compressive strength (most appropriate to the recovered meteorites) was closer to 0.7 MPa. These are much lower than values found for fireballs of ordinary chondritic composition. The behavior of the Tagish Lake fireball suggests that it represents the lowest end of the strength spectrum of carbonaceous chondrites or the high end of cometary meteoroids. The bulk density and porosity results for the Tagish Lake meteoroid suggest that the low bulk densities measured for some small primitive bodies in the solar system may reflect physical structure dominated by microporosity rather than macroporosity and rubble‐pile assemblages.  相似文献   

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

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