首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
Vertical profiles of CCl4, CFCl3, and CF2Cl2 mixing ratios in the upper troposphere and lower stratosphere have been measured on four flights with chartered aircraft, type HS 125. The flights were carried out in November and December 1976 over Europe at latitudes between 50 and 60°N. At least eight air samples were taken during each ascent and descent of the aircraft at altitudes between 7 and 12.5 km. The samples were analysed in the laboratory using gaschromatographic procedures. The results indicate a decrease of the CCl4, CFCl3, and CF2Cl2 mixing ratios above the tropopause. The observed average gradients in the stratosphere are 14 pptv/km for CCl4, 12 pptv/km for CFCl3 and 27.8 pptv/km for CF2Cl2. With exception of CFCl3 these gradients are higher than those predicted by model calculations. Apparently, further sink mechanisms for CCl4 and CF2Cl2 exist in the lower stratosphere not yet included in the models.  相似文献   

2.
Available data on halogenated molecules in the stratosphere will be reviewed. Presently vertical profiles of CFCl3 and CF2Cl2 in the stratosphere exist to 50 km altitude. Only measurements in the lower stratosphere are reported for the other major halocarbons, CCl4 and CH3Cl. Profiles of the product species ClO, Cl, HCl, and HF exist to about 35 km.Comparison with theoretical profiles from 1-D models shows generally good agreement except for ClO where the earliest measured concentrations exceed the calculated ones considerably.  相似文献   

3.
It is proposed that a component of meteoric smoke, sodium bicarbonate (NaHCO3), provides particularly effective condensation nuclei for noctilucent clouds. This assertion is based on three conditions being met. The first is that NaHCO3 is present at sufficient concentration (104 cm–3) in the upper mesosphere between 80 and 90 km. It is demonstrated that there is strong evidence for this based on recent laboratory measurements coupled with atmospheric modelling. The second condition is that the thermodynamics of NaHCO3(H2O)n cluster formation allow spontaneous nucleation to occur under mesospheric conditions at temperatures below 140 K. The Gibbs free energy changes for forming clusters with n = 1 and 2 were computed from quantum calculations using hybrid density functional/Hartree-Fock (B3LYP) theory and a large basis set with added polarization and diffuse functions. The results were then extrapolated to higher n using an established dependence of the free energy on cluster size and the free energy for the sublimation of H2O to bulk ice. A 1-dimensional model of sodium chemistry was then employed to show that spontaneous nucleation to form ice particles (n > 100) should occur between 84 and 89 km in the high-latitude summer mesosphere. The third condition is that other metallic components of meteoric smoke are less effective condensation nuclei, so that the total number of potential nuclei is small relative to the amount of available H2O. Quantum calculations indicate that this is probably the case for major constituents such as Fe(OH)2, FeO3 and MgCO3.  相似文献   

4.
Weekly air samples were collected at Cape Kumakahi (0 km) and at nearby Mauna Loa Observatory (3.4 km) which is above the boundary layer. EC/GC and GC/FID techniques were used to measure CH3I, CHCl3, CO and CH4, which are largely natural in origin, and C2Cl4, CCl4, CH3, CCl3, CCl3F (F-11), CCl2F2 (F-12), CHClF2 (F-22) and C2F3Cl3 (F-113), which are due to anthropogenic activities. It was found that all these gases are significantly (0.05) more abundant in the boundary layer than above it.  相似文献   

5.
A systematic approach was used to characterize the biosorption of copper(II) onto two biosorbents, Posidonia oceanica and peat, focusing on the interaction mechanisms, the copper(II) sorption–desorption process and the thermal behavior of the biosorbents. Sorption isotherms at pH 4–6 were obtained and the experimental data were fitted to the Langmuir model with a maximum uptake (qmax) at pH 6 of 85.78 and 49.69 mg g?1, for P. oceanica and peat, respectively. A sequential desorption (SD) with water, Ca(NO3)2, and EDTA was applied to copper‐saturated biosorbents. Around 65–70% copper(II) were desorbed with EDTA, indicating that this heavy metal was strongly bound. The reversibility of copper(II) sorption was obtained by desorption with HCl and SD. Fourier transform IR spectroscopy (FTIR) analysis detected the presence of peaks associated with OH groups in aromatic and aliphatic structures, CH, CH2, and CH3 in aliphatic structures, COO? and COOH groups and unsaturated aromatic structures on the surface of both biosorbents, as well as peaks corresponding to Si? O groups on the surface of peat. The results of SEM‐EDX and FTIR analysis of copper‐saturated samples demonstrated that ion exchange was one of the mechanisms involved in copper(II) retention. Thermal analysis of biosorbent samples showed that copper(II) sorption–desorption processes affected the thermal stability of the biosorbents.  相似文献   

6.
The implicit time integration scheme of Stott and Harwood (1993) was proposed as an efficient scheme for use in three-dimensional chemical models of the atmosphere. The scheme was designed for chemistry schemes using chemical families, in which species with short lifetimes are grouped into longer-lived families. Further study with more complex chemistry, more species and reactions showed the scheme to be non-convergent and unstable under certain conditions; particularly for the perturbed chemical scenarios of polar stratospheric winters. In this work the scheme has been improved by revising the treatment of families and the convergence properties of the scheme. The new scheme has been named IMPACT (IMPlicit Algorithm for Chemical Time-stepping). It remains easy to implement and produces simulations that compare well with integrations using more accurate higher order schemes.  相似文献   

7.
We present a comparison of the observed behavior of the F-region ionosphere over Millstone Hill during the geomagnetically quiet and storm periods of 6–12 April 1990 with numerical model calculations from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere. The major enhancement to the IZMIRAN model developed in this study is the use of a new loss rate of O+(4S) ions as a result of new high-temperature flowing afterglow measurements of the rate coefficients K1 and K2 for the reactions of O+(4S) with N2 and O2. The deviations from the Boltzmann distribution for the first five vibrational levels of O2(v) were calculated, and the present study suggests that these deviations are not significant. It was found that the difference between the non-Boltzmann and Boltzmann distribution assumptions of O2(v) and the difference between ion and neutral temperature can lead to an increase of up to about 3% or a decrease of up to about 4% of the calculated NmF2 as a result of a respective increase or a decrease in K2. The IZMIRAN model reproduces major features of the data. We found that the inclusion of vibrationally excited N2(v > 0) and O2(v > 0) in the calculations improves the agreement between the calculated NmF2 and the data on 6, 9, and 10 April. However, both the daytime and nighttime densities are reproduced by the IZMIRAN model without the vibrationally excited nitrogen and oxygen on 8 and 11 April better than the IZMIRAN model with N2(v > 0) and O2(v > 0). This could be due to possible uncertainties in model neutral temperature and densities, EUV fluxes, rate coefficients, and the flow of ionization between the ionosphere and plasmasphere, and possible horizontal divergence of the flux of ionization above the station. Our calculations show that the increase in the O+ + N2 rate factor due to N2(v > 0) produces a 5–36% decrease in the calculated daytime peak density. The increase in the O+ + O2 loss rate due to vibrationally excited O2 produces 8–46% reductions in NmF2. The effects of vibrationally excited O2 and N2 on Ne and Te are most pronounced during the daytime.  相似文献   

8.
Sources and sinks of atmospheric methane   总被引：7，自引：0，他引：7  
In 1972 average mixing ratio of methane in the troposphere was 1.41 ppm and 1.3 ppmv for the northern and southern hemisphere, respectively, which corresponds to a total amount of 4×1015 g of CH4 present in the atmosphere. Most is of recent biologic origin.14C analyses show that no more than 20 percent is released by fossil sources. The various ecosystems producing CH4 are discussed and the total annual production is estimated to lie between 5.5×1014 g/yr and 11×1014 g/yr. The corresponding turnover times for atmospheric CH4 range from 4 to 7 yrs. The destruction of CH4 takes place mainly in the troposphere, most probably through the reaction of CH4 + OH CH3 + H2O. About 10 percent of the CH4 is destroyed in the stratosphere. The CH4 cycle contributes on the order of 1 percent to the atmospheric carbon cycle.  相似文献   

9.
We examine the oscillations of the meridional neutral wind in the F region as seen by the EISCAT radar. We propose an interpretation in term of eddies (tourbillons) of typical size of a few tens to a few hundreds of kilometers. The observed rotation velocity is a few hundreds of meters per second. We suggest that the tourbillons are a common feature of thermospheric movements. We propose an optical experiment to check the validity of this assumption.  相似文献   

10.
The effect of present-day and future NOx emissions from aircraft on the NOx and ozone concentrations in the atmosphere and the corresponding radiative forcing were studied using a three-dimensional chemistry transport model (CTM) and a radiative model. The effects of the aircraft emissions were compared with the effects of the three most important anthropogenic NOx surface sources: road traffic, electricity generation and industrial combustion. From the model results, NOx emissions from aircraft are seen to cause an increase in the NOx and ozone concentrations in the upper troposphere and lower stratosphere, and a positive radiative forcing. For the reference year 1990, the aircraft emissions result in an increase in the NOx concentration at 250 hPa of about 20 ppt in January and 50 ppt in July over the eastern USA, the North Atlantic Flight Corridor and Western Europe, corresponding to a relative increase of about 50%. The maximum increase in the ozone concentrations due to the aircraft emissions is about 3-4 ppb in July over the northern mid-latitudes, corresponding to a relative increase of about 3-4%. The aircraft-induced ozone changes cause a global average radiative forcing of 0.025 W/m2 in July. According to the ANCAT projection for the year 2015, the aircraft NOx emissions in that year will be 90% higher than in the year 1990. As a consequence of this, the calculated NOx perturbation by aircraft emissions increases by about 90% between 1990 and 2015, and the ozone perturbation by about 50-70%. The global average radiative forcing due to the aircraft-induced ozone changes increases by about 50% between 1990 and 2015. In the year 2015, the effects of the aircraft emissions on the ozone burden and radiative forcing are clearly larger than the individual effects of the NOx surface sources. Taking chemical conversion in the aircraft plume into account in the CTM explicitly, by means of modified aircraft NOx emissions, a significant reduction of the aircraft-induced NOx and ozone perturbations is realised. The NOx perturbation decreases by about 40% and the ozone perturbation by about 30% in July over Western Europe, the eastern USA and the North Atlantic Flight Corridor.  相似文献   

11.
Ionospheric electron content (IEC) observed at Delhi (geographic co-ordinates: 28.63°N, 77.22°E; geomagnetic co-ordinates: 19.08°N, 148.91E; dip Latitude 24.8°N), India, for the period 1975/80 and 1986/89 belonging to an ascending phase of solar activity during first halves of solar cycles 21 and 22 respectively have been used to study the diurnal, seasonal, solar and magnetic activity variations. The diurnal variation of seasonal mean of IEC on quiet days shows a secondary peak comparable to the daytime peak in equinox and winter in high solar activity. IECmax (daytime maximum value of IEC, one per day) shows winter anomaly only during high solar activity at Delhi. Further, IECmax shows positive correlation with F10.7 up to about 200 flux units at equinox and 240 units both in winter and summer; for greater F10.7 values, IECmax is substantially constant in all the seasons. IECmax and magnetic activity (Ap) are found to be positively correlated in summer in high solar activity. Winter IECmax shows positive correlation with Ap in low solar activity and negative correlation in high solar activity in both the solar cycles. In equinox IECmax is independent of Ap in both solar cycles in low solar activity. A study of day-to-day variations in IECmax shows single day and alternate day abnormalities, semi-annual and annual variations controlled by the equatorial electrojet strength, and 27-day periodicity attributable to the solar rotation.  相似文献   

12.
Methane (CH4) is a potent greenhouse gas that is emitted from paddy fields, and the large CH4 fluxes represent a worldwide issue for the rice production eco-compatibility. In this work a model is proposed to investigate the role of water flows on CH4 emissions from flooded paddy soils. The model is based on a system of partial differential mass balance equations of the chemical species affecting CH4 fate, and water flows are modeled by the Darcy equation. Moreover, in order to properly model the dynamics of CH4, a number of physico-chemical processes and features not included in currently available CH4 emission models are considered: paddy soil stratigraphy; nutrient adsorption and root water uptake; gas transport and respiration within root aerenchyma compartment. The proposed model allows to simulate the spatio-temporal dynamics of chemical compounds within paddy soil as well as to quantify the influence of different processes on nutrient input/output budgets. Simulations without water flow have shown a considerable overestimation of CH4 emissions due to a different spatio-temporal dynamics of dissolved organic matter (DOC – source of energy for CH4 production). In particular, when water fluxes have not been modeled the overestimation can reach 54%, 41% and 67% of daily minimum, daily maximum, and total over the whole growing season CH4 emission, respectively. Moreover, the model results suggest that roots influence CH4 dynamics principally due to their nutrient uptake, while root effect on advective flow plays a minor role. Finally, the analysis of CH4 transport fluxes has shown the limiting effect of upward dispersive transport fluxes on the downward CH4 percolation.  相似文献   

13.
On August 21st 1998, a sharp southward turning of the IMF, following on from a 20 h period of northward directed magnetic field, resulted in an isolated substorm over northern Scandinavia and Svalbard. A combination of high time resolution and large spatial scale measurements from an array of coherent scatter and incoherent scatter ionospheric radars, ground magnetometers and the Polar UVI imager has allowed the electrodynamics of the impulsive substorm electrojet region during its first few minutes of evolution at the expansion phase onset to be studied in great detail. At the expansion phase onset the substorm onset region is characterised by a strong enhancement of the electron temperature and UV aurora. This poleward expanding auroral structure moves initially at 0.9 km s-1 poleward, finally reaching a latitude of 72.5°. The optical signature expands rapidly westwards at ~6 km s-1, whilst the eastward edge also expands eastward at ~0.6 km s-1. Typical flows of 600 m s-1 and conductances of 2 S were measured before the auroral activation, which rapidly changed to ~100 m s-1 and 10–20 S respectively at activation. The initial flow response to the substorm expansion phase onset is a flow suppression, observed up to some 300 km poleward of the initial region of auroral luminosity, imposed over a time scale of less than 10 s. The high conductivity region of the electrojet acts as an obstacle to the flow, resulting in a region of low-electric field, but also low conductivity poleward of the high-conductivity region. Rapid flows are observed at the edge of the high-conductivity region, and subsequently the high flow region develops, flowing around the expanding auroral feature in a direction determined by the flow pattern prevailing before the substorm intensification. The enhanced electron temperatures associated with the substorm-disturbed region extended some 2° further poleward than the UV auroral signature associated with it.  相似文献   

14.
A laboratory inorganic carbonate precipitation experiment at high pH of 8.96 to 9.34 was conducted, and the boron isotopic fractionations of the precipitated carbonate were measured. The data show that boron isotopic fractionation factors (αcarb-3) between carbonate and B(OH)3 in seawater range 0.937 and 0.965, with an average value of 0.953. Our results together with those reported by Sanyal and collaborators show that the αcarb-3 values between carbonate and B(OH)3 in solution are not constant but are negatively correlated with the pH of seawater. The measured boron isotopic compositions of carbonate precipitation (δ11Bcarb) do not exactly lie on the best-fit theoretical δ11B4-pH curves and neither do they exactly parallel any theoretical δ11B4-pH curves. Therefore, it is reasonable to argue that a changeable proportion of B(OH)3 with pH of seawater should also be incorporated into carbonate except for the dominant incorporation of B(OH)4 − in carbonate. Hence, in the reconstruction of the paleo-pH of seawater from boron isotopes in marine biogenic carbonates, the use of theoretical boron isotopic fractionation factor (α4−3) between B(OH)4 − and B(OH)3 is not suitable. Instead, an empirical equation should be established. Supported by National Natural Science Foundation of China (Nos. 40573013 and 40776071), State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences (Grant No SKLLQG0502) and State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences  相似文献   

15.
Warneck  P.  Junge  C. E.  Seiler  W. 《Pure and Applied Geophysics》1973,106(1):1417-1430
The consumption of both methane and carbon monoxide in the lower stratosphere is due predominantly to reaction with OH radicals. The possibility of deriving OH concentration from measurements of the decrease of CH4 and CO mixing ratios above the tropopause is explored. The observations and the basic chemistry are briefly summarized. Simple one-dimensional diffusion models are employed to derive expressions for the decrease of CH4 and CO mixing ratios with altitude above the tropopause, and the influence of important parameters is discussed. Vertical air velocities resulting from large-scale organized mean motion and from synoptic variations are shown to distort the concentration altitude profiles of methane and carbon monoxide, respectively. Suitable averaging of observational data is required to eliminate the effects due to vertical motion. Then a reliable value for the effective OH number density should be obtainable. At present an estimate of 4×106 molecules/cm3 is derived.  相似文献   

16.
During the Arctic Tropospheric Ozone Chemistry (ARCTOC) campaigns at Ny-Ålesund, Spitsbergen, the role of halogens in the depletion of boundary layer ozone was investigated. In spring 1995 and 1996 up to 30 ppt bromine monoxide were found whenever ozone decreased from normal levels of about 40 ppb. Those main trace gases and others were specifically followed in the UV-VIS spectral region by differential optical absorption spectroscopy (DOAS) along light paths running between 20 and 475 m a.s.l. The daily variation of peroxy radicals closely followed the ozone photolysis rate J(O3(O1D)) in the absence of ozone depletion most of the time. However, during low ozone events this close correlation was no longer found because the measurement of radicals by chemical amplification (CA) turned out to be sensitive to peroxy radicals and ClOx. Large CA signals at night can sometimes definitely be assigned to ClOx and reached up to 2 ppt. Total bromine and iodine were both stripped quantitatively from air by active charcoal traps and measured after neutron activation of the samples. Total bromine increased from background levels of about 15 ppt to a maximum of 90 ppt during an event of complete ozone depletion. For the spring season a strong source of bromine is identified in the pack ice region according to back trajectories. Though biogenic emission sources cannot be completely ruled out, a primary activation of halogenides by various oxidants seems to initiate an efficient autocatalytic process, mainly driven by ozone and light, on ice and perhaps on aerosols. Halogenides residing on pack ice surfaces are continuously oxidised by hypohalogenous acids releasing bromine and chlorine into the air. During transport and especially above open water this air mixes with upper layer pristine air. As large quantities of bromine, often in the form of BrO, have been observed at polar sunrise also around Antarctica, its release seems to be a natural phenomenon. The source strength of bromine from halogen activation on the pack ice, as based on the measured inorganic bromine levels, averages about 1012 Br-atoms m−2 s−1 during sunlit periods in Arctic spring. The total source strength of inorganic bromine from sunlit polar regions may therefore amount to 30 kt y−1.  相似文献   

17.
The present study provides an electrocoagulation method, for the removal of NO3– from drinking water using magnesium as the anode and cathode. The experiments are carried out as a function of pH, temperature, and current density. The results show that the maximum removal efficiency of 95.8% was achieved at a current density of 0.25 A/dm2, at a pH of 7.0. The adsorption of NO3– preferably fitting the Langmuir adsorption isotherm suggests monolayer coverage of the adsorbed molecules. The adsorption process follows a second‐order kinetics model. Thermodynamic studies show that the adsorption was exothermic and spontaneous in nature.  相似文献   

18.
This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N+2, O+2, O+(4S), O+(2D), O−(2P), and O+(2P) ions, and O(1D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm−3 s−1 was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N2(v > 0) and O2(v > 0) in the calculations of the O+(4S) loss rate improves the agreement between the calculated Ne and the data on 4 August, 1972. The N2(v > 0) and O2(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O+ + N2 rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O+ + N2 loss rate due to vibrationally excited O2 produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O2(a1δg) and O2(b1σg+), and rotational excitation of O2, and found that the effect of these cooling rates on Te can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(3P) with thermal electrons with excitation of O(1D), and the electron cooling rates by vibrational excitation of O2 and N2 are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N2(B3πg+), N2(C3πu), and N2(A3σu+) at several vibrational levels, O(1S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm.  相似文献   

19.
Low frequency stochastic variations of the geomagnetic AE-index characterized by 1/fb-like power spectrum (where f is a frequency) are studied. Based on the analysis of experimental data we show that the Bz-component of IMF, velocity of solar wind plasma, and the coupling function of Akasofu are insufficient factors to explain these behaviors of the AE-index together with the 1/fb fluctuations of geomagnetic intensity. The effect of self-organized criticality (SOC) is proposed as an internal mechanism to generate 1/fb fluctuations in the magnetosphere. It is suggested that localized spatially current instabilities, developing in the magnetospheric tail at the initial substorm phase can be considered as SOC avalanches or dynamic clusters, superposition of which leads to the 1/fb fluctuations of macroscopic characteristics in the system. Using the sandpile model of SOC, we undertake numerical modeling of space-localized and global disturbances of magnetospheric current layer. Qualitative conformity between the disturbed dynamics of self-organized critical state of the model and the main phases of real magnetospheric substorm development is demonstrated. It is also shown that power spectrum of sandpile model fluctuations controlled by real solar wind parameters reproduces all distinctive spectral features of the AE fluctuations.  相似文献   

20.
Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10−3 moles/kg) the gas concentrations are: 3–4.5 × 10−5 cm3 STP/kg helium, 0.62–1.24 cm3 STP/kg CH4, 10.80–16.71 × 10−3 moles/kg CO2. The samples contain large quantities (95–126 cm3/kg) of H2 and some carbon monoxide (0.26–0.36 cm3/kg) which result from reaction with the titanium sampling bottles. δ13C in methane and CO2 (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ13CCO is compatible with formation by reduction of CO2 on Ti at 350°C close to the sampling temperature.3He/4He are very homogeneous at (7.5 ± 0.1)RA(3He/4He = 1.0 × 10−5) and very similar to already published data as well as CH4/3He ratios between 1.4 and 2.1 × 106.18O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct18O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The3He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and3He mantle flux.  相似文献   

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

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