首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The increasing concentration of greenhouse gas in the atmosphere and their resultant climatic and environmental changes have been drawing much attention of the governments of various countries in recent years. The sphere of global influence and the comp…  相似文献   

2.
We present an uncertainty analysis of ecological process parameters and CO2 flux components (R eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO2 fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO2 fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO2 flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO2 flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated Reco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m−2·a−1), 34.33% (79 g C·m−2·a−1) and 5.4% (92 g C·m−2·a−1). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m−2·a−1), 2.1% (5.7 g C·m−2·a−1), and 0.26% (4.3 g C·m−2·a−1), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO2 flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of Reco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)  相似文献   

3.
三峡水库澎溪河消落区土-气界面CO2和CH4通量初探   总被引:1,自引:0,他引:1  
李哲  张利萍  王琳  郭劲松  高旭  方芳  蒋滔 《湖泊科学》2013,25(5):674-680
水库近岸湿地(消落区)温室气体(CO2、CH4)产汇是水库温室气体效应问题的重要组成部分.本文以三峡水库支流澎溪河的白家溪、养鹿两处大面积消落区为研究对象,于2010年6 9月水库低水位运行期间,对近岸消落区土-气界面CO2、CH4通量进行监测.白家溪消落区土-气界面CO2通量均值为12.38±2.42 mmol/(m2·h);CH4通量均值为0.0112±0.0064 mmol/(m2·h).养鹿消落区CO2、CH4通量均值分别为10.54±5.17、0.14±0.16 mmol/(m2·h).总体上,6 9月土-气界面CO2通量呈增加趋势,而CH4通量水平呈现显著的递减趋势.消落区土地出露后植被恢复,在一定程度上促进了土壤有机质含量的增加,使得6 9月CO2释放通量的总体趋势有所增加.消落区退耕后,其甲烷氧化菌的活性得到恢复,加之在土地出露曝晒过程中土壤透气性增强,使得消落区土壤对大气中CH4吸收氧化潜势增强.尽管如此,仍需进一步的研究以明晰消落区土-气界面CO2、CH4产汇的主要影响因素.  相似文献   

4.
Soil H2 and CO2 surveys were carried out along seven active faults and around the aftershock region of the 2000 Tottori-ken Seibu earthquake in Japan. Diffuse CO2 effluxes were also measured along one fault and around the 2000 aftershock region. The results show highly variable H2 concentration in space and time and it seems that the maximum H2 concentration at each active fault correlates with fault activity as exemplified by the time of the latest big earthquakes. Even though observed H2 concentrations in four faults were markedly lower than those collected previously in the latter half of the 1970s, it is evident that the higher H2 concentrations in this study are due to the addition of the fault gases. Comparing the chemical composition of trapped gases (H2: 5–20% and CO2/H2: 0.5–12) in fractured rocks of drill cores bored at the Nojima fault, a soil gas sample with the highest H2 concentration showed large amounts of the trapped fault gas, diluted with atmospheric component. The profile experiment across a fracture zone at the Yamasaki fault showed higher H2 concentrations and lower CO2/H2 ratios as was observed in soil gas from the fracture zone. A few days after the 2000 Tottori-kei Seibu earthquake, no CO2 effluxes related to the occurrence of earthquakes were observed at the aftershock region. However, only above the epicenter zone, relatively high H2 concentrations in soil gases were observed.  相似文献   

5.
Use of eddy covariance (EC) techniques to map the spatial distribution of diffuse volcanic CO2 fluxes and quantify CO2 emission rate was tested at the Horseshoe Lake tree-kill area on Mammoth Mountain, California, USA. EC measurements of CO2 flux were made during September–October 2010 and ranged from 85 to 1,766 g m−2 day−1. Comparative maps of soil CO2 flux were simulated and CO2 emission rates estimated from three accumulation chamber (AC) CO2 flux surveys. Least-squares inversion of measured eddy covariance CO2 fluxes and corresponding modeled source weight functions recovered 58–77% of the CO2 emission rates estimated based on simulated AC soil CO2 fluxes. Spatial distributions of modeled surface CO2 fluxes based on EC and AC observations showed moderate to good correspondence (R 2 = 0.36 to 0.70). Results provide a framework for automated monitoring of volcanic CO2 emissions over relatively large areas.  相似文献   

6.
The fundamental assumption of 210Pb sediment dating is the stable flux of 210Pbex, which was derived from atmosphere and then transferred into sediments via lake water. When the sedimentation rate is relatively constant, the 210Pbex activity in sediments will be exponentially reduced with sedimentation age. 210Pbex in lake water is incorporated into sediments mainly via organic particulates. If the sedimentation flux of organic matter in lake water is suddenly increased, 210Pbex will be significantly deposited and then transferred into sediments. On the one hand such sudden purification effect is obviously unfit for the fundamental assumption of 210Pb dating; on the other hand, the sudden enhancement of 210Pbex flux would be indicative of the conspicuous variation of primary productivity of lake water. This problem will be discussed in accordance with the variation trend of 210Pbex in the vertical profile of recent sediments of Lake Chenghai, Yunnan Province. The sediment core was collected from the deep-water area of Lake Chenghai in June 1997. The vertical profile of 137Cs activity is characterized by a tree-peak pattern. This profile gave reliable ages, and also showed the stability of sediment accumulation in the recent ten years. The vertical profile of 210Pbex activity displays a specific distribution of peaks, and is similar to the vertical profile of Corg. This phenomenon seems to be related to the mechanism of constraining the transfer of 210Pbex into lake sediments. The average atomic ratios of Horg/Corg and Corg/Norg in Lake Chenghai sediments are 5.51 and 7.04, respectively, indicating that the organic matter was predominantly derived from the remains of endogenic algae. In terms of the three-stage evolutionary characteristics of organic matter in sediments, i.e., “deposition-de-composition-accumulation”, the sedimentation fluxes (F(Corg)) of organic carbon (Corg) since 1970 were calculated by modeling. The sedimentation fluxes of 210Pbex (F(210Pbex)) in different years display good synchronous relations with the sedimentation flux of organic carbon (F(Corg)), especially in the years of 1972–1974 and 1986–1989. The variation of F(Corg) led to the variation of F(210Pbex); the variation of F(210Pbex) reflects, to some extent, the historical variation of lake productivity.  相似文献   

7.
This paper presents a robust H∞ output feedback control approach for structural systems with uncertainties in model parameters by using available acceleration measurements and proposes conditions for the existence of such a robust output feedback controller. The uncertainties of structural stiffness, damping and mass parameters are assumed to be norm-bounded. The proposed control approach is formulated within the framework of linear matrix inequalities, for which existing convex optimization techniques, such as the LMI toolbox in MATLAB, can be used effectively and conveniently. To illustrate the effectiveness of the proposed robust H∞ strategy, a six-story building was subjected both to the 1940 El Centro earthquake record and to a suddenly applied Kanai-Tajimi filtered white noise random excitation. The results show that the proposed robust H∞ controller provides satisfactory results with or without variation of the structural stiffness, damping and mass parameters.  相似文献   

8.
Soil CO2 concentration data were collected periodically from July 2001 to June 2005 from sampling site grids in two areas located on the lower flanks of Mt. Etna volcano (Paternò and Zafferana Etnea–Santa Venerina). Cluster analysis was performed on the acquired data in order to identify possible groups of sites where soil degassing could be fed by different sources. In both areas three clusters were recognised, whose average CO2 concentration values throughout the whole study period remained significantly different from one another. The clusters with the lowest CO2 concentrations showed time-averaged values ranging from 980 to 1,170 ppm vol, whereas those with intermediate CO2 concentrations showed time-averaged values ranging from 1,400 to 2,320 ppm vol, and those with the highest concentrations showed time-averaged values between 1,960 and 55,430 ppm vol. We attribute the lowest CO2 concentrations largely to a biogenic source of CO2. Conversely, the highest CO2 concentrations are attributed to a magmatic source, whereas the intermediate values are due to a variable mixing of the two sources described above. The spatial distribution of the CO2 values related to the magmatic source define a clear direction of anomalous degassing in the Zafferana Etnea–Santa Venerina area, which we attribute to the presence of a hidden fault, whereas in the Paternò area no such oriented anomalies were observed, probably because of the lower permeability of local soil. Time-series analysis shows that most of the variations observed in the soil CO2 data from both areas were related to changes in the volcanic activity of Mt. Etna. Seasonal influences were only observed in the time patterns of the clusters characterised by low CO2 concentrations, and no significant interdependence was found between soil CO2 concentrations and meteorological parameters. The largest observed temporal anomalies are interpreted as release of CO2 from magma batches that migrated from deeper to shallower portions of Etna’s feeder system. The pattern of occurrence of such episodes of anomalous gas release during the observation period was quite different between the two studied areas. This pattern highlighted an evident change in the mechanism of magma transport and storage within the volcano’s feeder system after June 2003, interpreted as magma accumulation into a shallow (<8 km depth) reservoir.  相似文献   

9.
785 traces of vertical components from shallow earthquakes recorded by 10 CDSN (Chinese Digital Seismographic Network) stations and 5 GSN (Global Seismographic Network) stations were collected to study the attenuation characteristics ofL g coda in the Chinese continent and its adjacent regions. The records were processed first using the stack spectral ratio method to obtain the average values ofQ 0 (Q at 1Hz) and η, the frequency dependence, ofL g coda in the ellipses corresponding to the paths. The back-projection technique was then employed to obtain the tomographic maps ofQ 0 and η values, and the distribution of their errors. Results indicate that in the studied areaQ 0 varies between 200 and 500. The lowest value ofQ 0 exists in the Yun-nan-Tibetan region, while the highest value ofQ 0 occurs in the southern edge of Siberian platform. η varies between 0.3 and 0.8. For most part of the studied area η varies inversely withQ 0.  相似文献   

10.
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  相似文献   

11.
On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO2 and SO2) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO2 flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO2 fluxes of the degassing plume and the CO2 fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO2 flux) and in soil degassing (CO2 flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.  相似文献   

12.
This paper provides a generic equation for the evaluation of the maximum earthquake magnitude mmax for a given seismogenic zone or entire region. The equation is capable of generating solutions in different forms, depending on the assumptions of the statistical distribution model and/or the available information regarding past seismicity. It includes the cases (i) when earthquake magnitudes are distributed according to the doubly-truncated Gutenberg-Richter relation, (ii) when the empirical magnitude distribution deviates moderately from the Gutenberg-Richter relation, and (iii) when no specific type of magnitude distribution is assumed. Both synthetic, Monte-Carlo simulated seismic event catalogues, and actual data from Southern California, are used to demonstrate the procedures given for the evaluation of mmax.The three estimates of mmax for Southern California, obtained by the three procedures mentioned above, are respectively: 8.32 ± 0.43, 8.31 ± 0.42 and 8.34 ± 0.45. All three estimates are nearly identical, although higher than the value 7.99 obtained by Field et al. (1999). In general, since the third procedure is non-parametric and does not require specification of the functional form of the magnitude distribution, its estimate of the maximum earthquake magnitude mmax is considered more reliable than the other two which are based on the Gutenberg-Richter relation.  相似文献   

13.
High-resolution sampling, measurements of organic carbon contents and 14C signatures of selected four soil profiles in the Haibei Station situated on the northeast Tibetan Plateau, and application of 14C tracing technology were conducted in an attempt to investigate the turnover times of soil organic carbon and the soil-CO2 flux in the alpine meadow ecosystem. The results show that the organic carbon stored in the soils varies from 22.12×104 kg C hm−2 to 30.75×104 kg C hm−2 in the alpine meadow ecosystems, with an average of 26.86×104 kg C hm−2. Turnover times of organic carbon pools increase with depth from 45 a to 73 a in the surface soil horizon to hundreds of years or millennia or even longer at the deep soil horizons in the alpine meadow ecosystems. The soil-CO2 flux ranges from 103.24 g C m−2 a−1 to 254.93 gC m−2 a−1, with an average of 191.23 g C m−2 a−1. The CO2 efflux produced from microbial decomposition of organic matter varies from 73.3 g C m−2 a−1 to 181 g C m−2 a−1. More than 30% of total soil organic carbon resides in the active carbon pool and 72.8%281.23% of total CO2 emitted from organic matter decomposition results from the topsoil horizon (from 0 cm to 10 cm) for the Kobresia meadow. Responding to global warming, the storage, volume of flow and fate of the soil organic carbon in the alpine meadow ecosystem of the Tibetan Plateau will be changed, which needs further research. Supported by the National Natural Science Foundation of China (Grant Nos. 40231015, 40471120 and 40473002) and the Guangdong Provincial Natural Science Foundation of China (Grant No. 06300102)  相似文献   

14.
以三峡水库香溪河库湾为研究对象,采用原位加密采样(2021年5月)和室内培养方法,结合沉积物特性与水环境因子分析,探讨了香溪河库湾沉积物甲烷(CH4)释放潜力、沉积物-水界面CH4产生和氧化通量空间分布规律及其影响因素。结果表明:三峡水库泄水期间,香溪河库湾沉积物CH4释放潜力的变化范围是6.35-2029.37 mg/(kg·d),沉积物-水界面CH4产生通量和氧化通量的变化范围分别为0.04~0.73、0.03~0.62 mmol/(m2·d);空间上,沉积物CH4释放潜力、沉积物-水界面CH4产生及氧化通量在香溪河库湾和各典型横切面(XX02、XX05和XX06)间表现出空间差异性,主要受水深、TOC和温度的影响。垂向上,CH4产生速率随沉积物深度的增加而减小,表层20 cm沉积物CH4释放潜力占整柱沉积物的70%,可以用于估算库湾沉积物CH4释放潜力。此外,沉积...  相似文献   

15.
Sliding-window cross-correlation method is firstly adopted to identify sPn phase, and to constrain focal depth from regional seismograms, by measuring the time separation between sPn and Pn phases. We present the focal depths of the 17 moderate-sized aftershocks (M S⩾5.0) of the Wenchuan M S8.0 earthquake, using the data recorded by the regional seismic broadband networks of Shaanxi, Qinghai, Gansu, Yunnan and Sichuan. Our results show focal depths of aftershocks range from 8 to 20 km, and tend to cluster at two average depths, separate at 32.5°N, i.e., 11 km to the south and 17 km to the north, indicating that these aftershocks are origin of upper-to-middle crust. Combined with other results, we suggest that the Longmenshan fault is not a through-going crustal fault and the Pingwu-Qingchuan fault may be not the northward extension of the Longmenshan thrust fault. Supported by the National Natural Science Foundation of China (Grant Nos. 40604009 and 40574040) and Special Project for the Fundamental R & D of Institute of Geophysics, China Earthquake Administration (Grant No.DQJB08B20)  相似文献   

16.
A detailed analysis of published data on the N2, Ar, and He content and Ar and He isotopic composition of fumarolic fluids from Vulcano crater (south Italy) supports a model with two endmembers comprising magmatic and hydrothermal fluids with correspondingly low and high H2O content. The magmatic component with the highest 3He/4He and highest absolute concentrations of N2, Ar, and He also has the lowest N2/Ar and N2/He ratios (∼300 and ∼500, respectively). In contrast, the hydrothermal endmember, with the lower 3He/4He and lower absolute N2, Ar, and He abundances, has high N2/Ar (∼1,000) and high N2/He (>3,000) ratios. The hydrothermal component is also characterized by the highest 40Ar/36Ar ratios (>1,000) and is proposed to be the main carrier of metamorphic gases from the arc crust.  相似文献   

17.
Expeditions to Muztagata (in the eastern Pamirs) during the summer seasons of 2002 and 2003 collected precipitation samples and measured their oxygen isotopes. The δ 18O in precipitation displays a wide range, varying from −17.40‰ to +1.33‰ in June-September 2002 and from −22.31‰ to +4.59‰ in May-August 2003. The δ 18O in precipitation correlates with the initial temperature of precipitation during the observing periods. The positive correlation between δ 18O and temperature suggests that δ 18O can be used as an indicator of temperature in this region. The δ 18O values in fresh-snow samples collected from two snow events at different elevations on the Muztagata Glacier show a strong “altitude effect”, with a ratio of nearly −0.40% per 100 m from 5500 m to 7450 m.  相似文献   

18.
N2O concentrations and denitrification-related factors (NO3, SO4, dissolved organic carbon (DOC) and CO2) were investigated in the surface groundwater of a catchment in northern Germany, the Fuhrberger Feld Aquifer (FFA). We sampled 79 plots that were selected according to the three criteria of land use, historical land use conversion (1954–1995) and groundwater level. We sampled three sites within each plot. The sampling depth was 0.5 m below the groundwater surface.We found no indication for the occurrence of autotrophic denitrification in the surface groundwater. Heterotrophic denitrification was identified as the main process for N2O accumulation. The variability of N2O concentrations on the plot-scale was extremely high and was poorly explained by the three sampling criteria. Other denitrification-related variables such as NO3, SO4 and DOC were less variable. The selection criteria land use and groundwater level clearly influenced the order of magnitude of N2O concentrations in the surface groundwater. Under arable land, high NO3 concentrations resulted in high N2O concentrations. The surface groundwater under forest and pasture was almost NO3-free and had also very small N2O concentrations. Plots where the distance from the soil surface to the groundwater surface was large (>1 m up to 3.4 m) showed higher N2O concentrations in the surface groundwater than plots where the distance was small (<1 m). A larger distance from the soil surface to the groundwater leads to a longer residence time and more decomposition of DOC in the soil. Consequently the less bioavailable DOC could inhibit the efficiency of the heterotrophic denitrification in the groundwater, yielding more N2O. Elevated organic carbon levels in plots with historic land use conversion (pasture to arable) were very stable and did not influence N2O concentrations. The high within plot variability showed that an upscaling of N2O from the plot-scale to the catchment-scale is possible as long as the groundwater level regime and the land use do not change.  相似文献   

19.
The oxygen minimum zones (OMZs) are recognized as intense sources of N2O greenhouse gas (GHG) and could also be potential sources of CO2, the most important GHG for the present climate change. This study evaluates, for one of the most intense and shallow OMZ, the Chilean East South Pacific OMZ, the simultaneous N2O and CO2 fluxes at the air–sea interface. Four cruises (2000–2002) and 1 year of monitoring (21°–30°–36°S) off Chile allowed the determination of the CO2 and N2O concentrations at the sea surface and the analysis of fluxes variations associated with different OMZ configurations. The Chilean OMZ area can be an intense GHG oceanic local source of both N2O and CO2. The mean N2O fluxes are 5–10 times higher than the maximal previous historical source in an OMZ open area as in the Pacific and Indian Oceans. For CO2, the mean fluxes are also positive and correspond to very high oceanic sources. Even if different coupling and decoupling between N2O and CO2 are observed along the Chilean OMZ, 65% of the situations represent high CO2 and/or N2O sources. The high GHG sources are associated with coastal upwelling transport of OMZ waters rich in N2O and probably also in CO2, located at a shallow depth. The integrated OMZ role on GHG should be better considered to improve our understanding of the past and future atmospheric CO2 and N2O evolutions.  相似文献   

20.
On the basis of Parry’s method (1986), an improved method was established to determine the molar volume (Vm) and compositions (X) of the NaCl-H2O-CO2 (NHC) system inclusion. To use this method, the determination of Vm-X only requires three microthermometric data of a NHC inclusion: partial homog-enization temperature (Th ,CO2), salinity (S) and total homogenization temperature (Th). Theoretically, four associated equations are needed containing four unknown parameters: X CO2, XNaCl, Vm and F (volume fraction of CO2 phase in total inclusion when occurring partial homogenization). When they are released, the Vm-X are determined. The former three equations, only correlated with Th ,CO2, S and F, have simplified expressions:XCO2=f1(Th,CO2,S,F),XNaCl=f2(Th,CO2,S,F),Vm=f3(Th,CO2,S,F). The last one is the thermodynamic relationship of X CO2, XNaCl, Vm and Th:f4(XCO2,XNaCl,Vm,Th)=0.Since the above four associated equations are complicated, it is necessary to adopt iterative technique to release them. The technique can be described by:(i) Freely input a F value (0≤F≤1),with Th ,CO2 and S, into the former three equations. As a result,X CO 2,XNaCl and the molar volume value recorded as Vm1 are derived. (ii) Input the X CO2 and XNaCl gotten in the step above into the last equation, and another molar volume value recorded as Vm2 is determined. (iii) If Vm1 is unequal to Vm2, the calculation will be restarted from “(i)”. The iteration is completed until Vm1 is equal to Vm2, which means that the four associated equations are released. Compared to Parry’s (1986) solution method, the improved method is more convenient to use, as well as more accurate to determine X CO 2. It is available for a NHC inlusion whose partial homogenization temperature is higher than clatherate melting temperature and there are no solid salt crystals in the inclusion at parital homogenization.  相似文献   

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

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