首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
The Archean atmospheric oxygen concentration and sulfur cycle was long debated. The banded iron formation (BIF) is a special type of the sedimentary formation, which has truly recorded the atmospheric and oceanic conditions at that time. In this study, the composition of multiple sulfur isotope (δ 34S/δ 33S/δ 32S) for sulfides bedded in the Archean (~2.7 Ga) BIFs, in Anshan-Benxi area of Liaoning Province has been measured. The value of △33S varies from -0.89‰ to 1.21‰, which shows very obvious mass-independent fractionation (MIF) signatures. These non-zero △33S values indicate that the Archean sulfur cycles are different from what it is today, which have been deeply influenced by gas phase photochemical reactions. Algoma-type BIFs which are closely related to the volcanic activity have negative △33S value, however, Superior-type BIFs which are far away from the volcanic center have positive △33S value. The δ 34S varies in a large range from -22.0‰ to 11.8‰, which indicates that the bacteria reduction activity has already existed at that time, and that the oceanic sulfate concentration has at least reached 1 mmol/L in local areas. Combined with the contemporaneous existence of the hematite, magnetite and the occurrence and preservation of the sulfur MIF, it can be inferred that the Archean atmospheric oxygen level must be at 10-2―10-3 of the present atmospheric level (PAL).  相似文献   

2.
In-situ Hf isotope analyses and U–Pb dates were obtained by laser ablation-MC-ICP-MS for a zircon-bearing mantle eclogite xenolith from the diamondiferous Jericho kimberlite located within the Archean Slave Province (Nunavut), Canada. The U–Pb zircon results yield a wide range of ages (2.0 to 0.8 Ga) indicating a complex geological history. Of importance, one zircon yields a U–Pb upper intercept date of 1989 ± 67 Ma, providing a new minimum age constraint for zircon crystallization and eclogite formation. In contrast, Hf isotope systematics for the same zircons display an intriguing uniformity, and corresponding Hf depleted mantle model ages range between 2.1 ± 0.1 and 2.3 ± 0.1 Ga; the youngest Hf model age is within error to the oldest U–Pb date.

The Jericho eclogites have previously been interpreted as representing remnants of metamorphosed oceanic crust, and their formation related to Paleoproterozoic subduction regimes along the western margin of the Archean Slave craton during the Wopmay orogeny. Hf isotope compositions and U–Pb results for the Jericho zircons reported here are in good agreement with a Paleoproterozoic subduction model, suggesting that generation of oceanic crust and eclogite formation occurred between 2.0 and 2.1 Ga. The slightly older Hf depleted mantle model ages (2.1 to 2.3 Ga) may be reconciled with this model by invoking mixing between ‘crustal’-derived Hf from sediments and more radiogenic Hf associated with the oceanic crust during the 2 Ga subduction event. This results in intermediate Hf isotope compositions for the Jericho zircons that yield ‘fictitiously’ older Hf model ages.  相似文献   


3.
4.
The sand-hill section in Hongguang, Pengze County, Jiangxi Province is composed of eolian sand beds interbedded with 12 weakly-pedogenetic paleosols, and this sedimentary sequence was formed from about 25 to 13 kaBP. A detailed rock-magnetic study demonstrates that the magnetic minerals in the sequence are dominated by multi-domain magnetite grains in addition to small amounts of maghemite and hematite. Compared with the eolian sand beds, the paleosols are characterized by smaller magnetic particles, higher proportion of low-coercivity minerals and markedly increase in the magnetic mineral content. In the section, various magnetic parameters (K, ARM, SIRM and S ratios) display a roughly gradual increase from eolian sand beds to the overlying paleosols but an abrupt decrease from paleosols to the overlying eolian sand beds. The content of low-coercivity magnetic minerals shows the lowest values in the second eolian sand bed from the bottom of the section, and from the bed upwards it gradually increases on the whole in the eolian sand beds. These magnetic properties reflect regional climatic changes during the last glacial maximum.  相似文献   

5.
186Os enrichments in volcanic rocks and peridotite-derived iridosmine grains have been attributed to contributions from Earth’s outer core to the mantle, and apparently constrain the scale of mantle convection and an early timing for inner–outer core segregation more than 3.5 Gyr ago. Here, we highlight that marine ferromanganese crusts and nodules are characterised by high Pt/Os ratios and Pt–Os contents that develop much larger 186Os excesses over geological time (≥0.2%/Gyr) than those hypothesised for Earth’s outer core (<0.005–0.01%/Gyr). 187Os/188Os ratios in ferromanganese crusts are radiogenic due to sequestering of continental Os from seawater. Similarly, ancient ferromanganese materials may have had 186Os excesses (>0.1%) as a result of high Pt/Os ratios in continental crust, even prior to in-growth of 186Os after formation due to their high Pt/Os ratios. Past recycling of small amounts of these materials into the Earth’s mantle will produce coupled 187Os–186Os excesses and little change in Re and platinum-group-element concentrations, as observed in Hawaiian picrites, and in contrast to the predicted result of outer core addition to the mantle. 187Os and 186Os enrichments in the Hawaiian mantle source are potentially consistent with it comprising recycled oceanic lithosphere, pelagic sediments and ferromanganese materials, and questions the notion that Os isotopes can be used to uniquely identify core–mantle interactions and the depth at which mantle sources for volcanism originate.  相似文献   

6.
Summary Two different spectral methods have recently been used to model the flow driven by harmonic loads in a Newtonian mantle with laterally variable viscosity. The first method, by Zhang and Christensen (1993), transforms the problem with a general three-dimensional viscosity into a series of standard spherically symmetric problems. A different approach has been chosen by Martinec et al. (1993). Their method is based on integral formulation of the problem. The solution, which corresponds to a minimum of the dissipative energy, is found by means of the gradient search. We have tested the efficiency and numerical behaviour of both methods. The results of our tests favour the former method which is found more accurate and significantly faster than the gradient algorithm.  相似文献   

7.
Melt inclusions from four individual lava samples representing the HIMU (Mangaia Island), EMI (Pitcairn Island) and EMII (Tahaa Island) end member components, have heterogeneous Pb isotopic composition larger than that defined by the erupted lavas in each island. The broad linear trend in 207Pb/206Pb–208Pb/206Pb space produced by the melt inclusions from Mangaia, Tahaa and fPitcairn samples reproduces the entire trend defined by the Austral chain, the Society islands and the Pitcairn island and seamount groups. The inclusions preserve a record of melt composition of far greater isotopic diversity than that sampled in whole rock basalts. These results can be explained by mixing of a common depleted component with the HIMU, EMI and EMII lavas, respectively. We favor a model that considers the oceanic lithosphere to be that common component. We suggest that the Pb isotopic compositions of the melt inclusions reflect wall rock reaction of HIMU, EMI and EMII melts during their percolation through the oceanic lithosphere. Under these conditions, the localized rapid crystallization of olivine from primitive basalt near the reaction zone would allow the entrapment of melt inclusions with different isotopic composition.  相似文献   

8.
In a traditional analytical method, the convective features of Earth’s mantle have been inferred from surface signatures obtained by the geodynamic model only with depth-dependent viscosity structure. The moving and subducting plates, however, bring lateral viscosity variations in the mantle. To clarify the effects of lateral viscosity variations caused by the plate-tectonic mechanism, I have first studied systematically instantaneous dynamic flow calculations using new density-viscosity models only with vertical viscosity variations in a three-dimensional spherical shell. I find that the geoid high arises over subduction zones only when the vertical viscosity contrast between the upper mantle and the lower mantle is O(103) to O(104), which seems to be much larger than the viscosity contrast suggested by other studies. I next show that this discrepancy may be removed when I consider the lateral viscosity variation caused by the plate-tectonic mechanism using two-dimensional numerical models of mantle convection with self-consistently moving and subducting plates, and suggest that the observed geoid anomaly on the Earth’s surface is significantly affected by plate-tectonic mechanism as a first-order effect.  相似文献   

9.
Abstract Recent geological investigations of the Isua Supracrustal Belt (3.8 Ga), southern West Greenland, have suggested that it is the oldest accretionary complex on earth, defined by an oceanic plate‐type stratigraphy and a duplex structure. Plate history from mid‐oceanic ridge through plume magmatism to subduction zone has been postulated from analysis of the reconstructed oceanic plate stratigraphy in the accretionary complex. Comparison between field occurrence of greenstones in modern and ancient accretionary complexes reveals that two types of tholeiitic basalt from different tectonic settings, mid‐oceanic ridge basalt (MORB) and oceanic island basalt (OIB), occur. This work presents major, trace and rare earth element (REE) compositions of greenstones derived from Isua MORB and OIB, and of extremely rare relict igneous clinopyroxene in Isua MORB. The Isua clinopyroxenes (Cpx) have compositional variations equivalent to those of Cpx in modern MORB; in particular, low TiO2 and Na2O contents. The Isua Cpx show slightly light (L)REE‐depleted REE patterns, and the calculated REE pattern of the host magma is in agreement with that of Isua MORB. Analyses of 49 least‐altered greenstones carefully selected from approximately 1200 samples indicate that Isua MORB are enriched in Al2O3, and depleted in TiO2, FeO*, Y and Zr at the given MgO content, compared with Isua OIB. In addition, Isua MORB show an LREE‐depleted pattern, whereas Isua OIB forms a flat REE pattern. Such differences suggest that the Early Archean mantle had already become heterogeneous, depending on the tectonic environment. Isua MORB are enriched in FeO compared with modern MORB. Comparison of Isua MORB with recent melting experiments shows that the source mantle had 85–87 in Mg? and was enriched in FeO. Potential mantle temperature is estimated to be approximately 1480°C, indicating that the Early Archean mantle was hotter by at most approximately 150°C than the modern mantle.  相似文献   

10.
The local meteoric water line (LMWL), the functional relationship between locally measured values of δ18O and δ2H in precipitation, represents the isotopic composition of water entering hydrologic systems. The degree to which the LMWL departs from the global meteoric water line (GMWL), moreover, can reveal important information about meteoric sources of water (e.g. oceanic or terrestrial) and atmospheric conditions during transport. Here we characterize the isotopic composition of precipitation within an experimental watershed in the Western US that is subject to large topographic and seasonal gradients in precipitation. Interpreting the hydrometeorologic and spatial controls on precipitation, we constructed a seasonally weighted LMWL for southwestern Idaho that is expressed by the equation δ2H = 7.40 × δ18O ? 2.17. A seasonally weighted LMWL that is based on weighting isotopic concentrations by climatic precipitation volumes is novel, and we argue better represents the significant seasonality of precipitation in the region. The developed LMWL is considerably influenced by the semiarid climate experienced in southwest Idaho, yielding a slope and y‐intercept lower than the GMWL (δ2H = 8 × δ18O + 10). Moderate to strong correlations exist between the isotopic composition of precipitation from individual events and surface meteorologic variables, specifically surface air temperature, relative humidity, and precipitation amount. A strong negative correlation exists between the annual average isotopic composition of precipitation and elevation at individual collection sites, with a lapse rate of ?0.22‰/100 m. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
Dissolved organic carbon (DOC) originating in peatlands can be mineralized to carbon dioxide (CO2) and methane (CH4), two potent greenhouse gases. Knowledge of the dynamics of DOC export via run‐off is needed for a more robust quantification of C cycling in peatland ecosystems, a prerequisite for realistic predictions of future climate change. We studied dispersion pathways of DOC in a mountain‐top peat bog in the Czech Republic (Central Europe), using a dual isotope approach. Although δ13CDOC values made it possible to link exported DOC with its within‐bog source, δ18OH2O values of precipitation and run‐off helped to understand run‐off generation. Our 2‐year DOC–H2O isotope monitoring was complemented by a laboratory peat incubation study generating an experimental time series of δ13CDOC values. DOC concentrations in run‐off during high‐flow periods were 20–30 mg L?1. The top 2 cm of the peat profile, composed of decaying green moss, contained isotopically lighter C than deeper peat, and this isotopically light C was present in run‐off in high‐flow periods. In contrast, baseflow contained only 2–10 mg DOC L?1, and its more variable C isotope composition intermittently fingerprinted deeper peat. DOC in run‐off occasionally contained isotopically extremely light C whose source in solid peat substrate was not identified. Pre‐event water made up on average 60% of the water run‐off flux, whereas direct precipitation contributed 40%. Run‐off response to precipitation was relatively fast. A highly leached horizon was identified in shallow catotelm. This peat layer was likely affected by a lateral influx of precipitation. Within 36 days of laboratory incubation, isotopically heavy DOC that had been initially released from the peat was replaced by isotopically lighter DOC, whose δ13C values converged to the solid substrate and natural run‐off. We suggest that δ13C systematics can be useful in identification of vertically stratified within‐bog DOC sources for peatland run‐off.  相似文献   

12.
The similarity solutions of the equations of kinematic flow on a cone-shaped surface, provide a means for checking the numerical models used for converging and diverging flows. These models can be further improved if the initial excess rainfall curve can be described as a power of time. Illustration of this improvement and assessment of numerical models are given for a constant excess rainfall and for excess rainfall proportional to the inverse of the square root of time.  相似文献   

13.
Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

14.
To improve understanding of Ca isotope transport during water-rock interaction on the continents, we measured dissolved δ44Ca values along a 236 km flow path in the Madison aquifer, South Dakota, where fluids have chemically evolved according to dolomite and anhydrite dissolution, calcite precipitation, and Ca-for-Na ion-exchange over a timescale spanning ~ 15 kyr. We used a reactive transport model employing rate data constrained from major ion mass-balances to evaluate the extent to which calcite precipitation and ion-exchange fractionate Ca isotopes. Elevated δ44Ca values during the initial and final stages of water transport possibly result from calcite precipitation under supersaturated conditions and Ca-for-Na ion-exchange, respectively. However, for the bulk of the flow path, δ44Ca values evolve by mixing between anhydrite and dolomite dissolution, with no fractionation during calcite precipitation under saturated conditions. We attribute the absence of Ca isotope fractionation to the long timescale of water-rock interaction and slow rate of calcite precipitation, which have enabled fluids to chemically and isotopically equilibrate with calcite. We therefore conclude that the equilibrium Ca isotope fractionation factor between calcite and water (Δcal–w) is very close to zero. To the extent that the Madison aquifer typifies other groundwater systems where calcite slowly precipitates from solutions at or near chemical equilibrium, this study suggests that groundwater contributions to δ44Ca variability on the continents can be modeled according to simple mixing theory without invoking isotope discrimination.  相似文献   

15.
Chemical heterogeneities in the Martian mantle are believed to result from the crystallization of a magma ocean in the first 100 million years of its history. Shergottite meteorites from Mars are thought to retain a compositional record of such early differentiation and the resulting mineralogy at different depths. The coupled 176Lu–176Hf and 147Sm–143Nd isotope systematics in 9 shergottites are used here to investigate these issues. Three compositional groups in the shergottites display distinct isotope systematics. One group, commonly termed as depleted, is characterized by positive 176Hfi from + 46.2 to + 50.4 and 143Ndi from + 36.2 to + 39.1. Another, termed as enriched, has negative 176Hfi = − 16.5 to − 13.2 and 143Ndi = − 7.0 to − 6.5. The third group is intermediate between the depleted and enriched groups with positive 176Hfi = + 30.0 to + 33.4 and 143Ndi = + 16.9. Together, they describe mixing curves between 176Hf/177Hf, 143Nd/144Nd, Lu/Hf, and Sm/Nd, implying that they sample two distinct sources in the Martian mantle. All shergottites are characterized by (Sm/Nd)source < (Sm/Nd)sample, but (Lu/Hf)source > (Lu/Hf)sample. This decoupling can be explained by two successive partial melting episodes in the depleted shergottite source and localized in the Martian upper mantle. The genesis of shergottites can be modeled using non-modal equilibrium partial melting in a source initially composed of 60% olivine, 21% clinopyroxene, 9% orthopyroxene, and 10% garnet, with degrees of partial melting of 8.8% and 3.9%, respectively, for the two successive events. The enriched end-member of the shergottite mixing curve is best modeled by late-stage quenched residual melt resulting from the crystallization of a magma ocean. The depleted shergottite source may be modeled as a mixture of cumulates and residual melt, as convection in the Martian magma ocean is expected to reduce the incompatible trace element heterogeneity in the final solidified layers. Consequently, equilibrium crystallization is preferred to model the crystallization of the Martian magma ocean. The models that best explain the shergottite data are those where the magma ocean is at a depth of at least 1350 km in Mars.  相似文献   

16.
Strontium, neodymium isotopic compositions and trace elements of the detrital sediments of Core NS90-103 from South China Sea were analyzed. The results show that the87Sr/86Sr ratios of the detritus during the last glacial range from 0.722 4 to 0.723 0. They are significantly higher than those during the Holocene and the maximum of the last interglacial, which range from 0.721 0 to 0.721 7. This indicates stronger continental weathering during the last glacial. On the other hand, the143Nd/144Nd ratios of these detritus are higher during the last glacial too, similar to the variation of the87Sr/86Sr ratios. The trace element geochemistry of these detritus indicates that more authigenic sediments, such as ferromanganese, during the last glacial may partly contribute to the increase of143Nd/144Nd ratios. Furthermore, much more detritus from continent of South China to the north of the South China Sea may probably contribute to143Nd/144Nd ratios increase during the last glacial, which was the result of the enhancement of northeast monsoon.  相似文献   

17.
利用地震面波频散重建川滇地区壳幔S波速度   总被引:3,自引:2,他引:3       下载免费PDF全文
张智  陈赟  李飞 《地球物理学报》2008,51(4):1114-1122
利用适配滤波频时分析技术分析覆盖川滇地区的长周期面波记录,计算了周期10~100 s内的面波群速度频散,对研究区进行划分尺度大小1.5°×1.5°分格后,采用射线追踪方法求取各分段射线的长度和时间,得到各个格子的纯路径频散.继而采用阻尼最小二乘法求解,反演得到该研究区壳幔S波速度分布.研究结果表明,川滇地区表现出地壳增厚和缩短,在地壳和上地幔顶部,川滇菱形块体内部与其外部相比,虽然存在局部速度负异常,总体上呈相对高速,其周边的走滑断裂带呈现深至上地幔顶部的负速度异常,这有助于地壳块体沿断裂的侧向挤出;此外,云南西部和四川西部壳内和上地幔高导层的存在被认为是与部分熔融的物质或与滑脱构造相关联;从纬向剖面和经向剖面可以得到四川盆地莫霍面平均深度大约为45 km,云南地区莫霍面深度南北方向不一致,云南地区最北端深度达到49 km,南端莫霍面深度大约为36 km,这说明不同构造块体在构造运动过程中受到影响的程度不同.  相似文献   

18.
Temporal and spatial variations of stable oxygen (18O) and hydrogen (2H) isotope measurements in precipitation act as important proxies for changing hydro‐meteorological and regional and global climate patterns. Temporal trends in time series of the stable isotope composition in precipitation were rarely observed, and they are poorly understood. These might be a result of a lack of proper trend detection tools and effort for exploring trend processes. Here, we investigate temporal trends of δ18O in precipitation at 17 observation stations in Germany between 1978 and 2009. We test if significant trends in the isotope time series from different models can be observed. Mann–Kendall trend tests are applied on the isotope series, using general multiplicative seasonal autoregressive integrate moving average (ARIMA) models, which account for first and higher order serial correlations. Effects of temperature, precipitation, and geographic parameters on isotope trends are also investigated in the proposed models. To benchmark our proposed approach, the ARIMA results are compared with a trend‐free pre‐whitening procedure, the state of the art method for removing the first order autocorrelation in environmental trend studies. Moreover, we further explore whether higher order serial correlations in isotope series affects our trend results. Overall, three out of the 17 stations show significant changes when higher order autocorrelation are adjusted, and four show a significant trend when temperature and precipitation effects are considered. The significant trends in the isotope time series generally occur only at low elevation stations. Higher order autoregressive processes are shown to be important in the isotope time series analysis. Results suggest that the widely used trend analysis with only the first order autocorrelation adjustment may not adequately take account of the high order autocorrelated processes in the stable isotope series. The investigated time series analysis method including higher autocorrelation and external climate variable adjustments is shown to be a better alternative. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
Despite growing interest in soil erosion on agricultural land, relatively little attention has been paid to the influence of erosion processes on the pattern of contemporary landform evolution. This in part reflects the problems associated with up-scaling the results of short-term process studies to temporal and spatial scales relevant to the study of landform evolution. This paper presents a new approach to examining the influence of erosion processes on landform evolution on agricultural land which employs: caesium-137 (137Cs) measurements to provide medium-term (c. 40 years) estimates of rates of landform change; experimental data and a topographic-based model to simulate soil redistribution by tillage; a mass-balance model of 137Cs redistribution to separate the water erosion and tillage components of the 137Cs ‘signatures’; and field observations of water erosion for validation. This approach is used to examine the relative importance of water erosion and tillage processes for contemporary landform evolution at contrasting sites near Leuven, in Belgium, and near Yanan, in Shaanxi Province, China. This application of the approach provides good agreement between the derived water erosion rates and field observations, and hitherto unobtainable insights into medium-term patterns and rates of contemporary landform evolution. At Huldenberg in Belgium, despite rill incision of slope concavities and ephemeral gully incision of the valley floor, contemporary landform evolution is dominated by infilling of slope and valley concavities (rates >0.5 mm a−1) and gradual lowering of slope angles as a result of tillage. In contrast, at Ansai (near Yanan) the slope is characterized by increase in slope angle over most of the length, recession of the steepest section at a rate >5 mm a−1 and by increasing planform curvature. At this site, contemporary landform evolution is dominated by water erosion. The constraints on the approach are examined, with particular attention being given to limitations on extrapolation of the results and to the sensitivity of the models to parameter variation. © 1997 by John Wiley & Sons, Ltd.  相似文献   

20.
Picrites from the 61 million year old Vaigat Formation of the Nuussuaq Peninsula in West Greenland have 3He/4He ratios trapped in olivine phenocrysts which range up to 30 times the atmospheric ratio. These high values, measured during gas extraction by crushing in vacuum, are similar to the highest magmatic 3He/4He ratios found in young terrestrial volcanic rocks. By analogy with young basalts, in which crushing selectively extracts magmatic helium, any significant cosmogenic 3He appears to be absent in these picrites. Additional evidence for the absence of cosmogenic helium is provided by fusion results on the crushed olivine powders and by a single stepwise crushing experiment, in which only magmatic and radiogenic helium components are resolvable. The West Greenland picrites have Pb, Nd and Sr isotope compositions which overlap those found in picrites from Iceland and in basalts from Loihi Seamount, localities which today also have high 3He/4He ratios. Isotopic variations in He, Pb, Nd and Sr for the West Greenland picrites are interpreted to largely result from interaction of the early Iceland mantle plume with the upper mantle during plume ascent and dispersion beneath the continental lithosphere. The presence of high 3He/4He ratios in West Greenland, and the onset of magmatism across the North Atlantic Volcanic Province near 62 Ma, supports the hypothesis for very rapid dispersion (>1 m/year) of mantle plume head material during the earliest stages of plume impact, as predicted in recent numerical simulations of plume behavior during thermal mantle convection with non-Newtonian rheology.  相似文献   

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

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