首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 12 毫秒
1.
Macrophytes are a critical component of lake ecosystems affecting nutrient and contaminant cycling, food web structure, and lake biodiversity. The long-term (decades to centuries) dynamics of macrophyte cover are, however, poorly understood and no quantitative estimates exist for pre-industrial (pre-1850) macrophyte cover in northeastern North America. Using a 215 lake dataset, we tested if surface sediment diatom assemblages significantly differed among lakes that have sparse (<10% cover; group 1), moderate (10–40% cover; group 2) or extensive (>40% cover; group 3) macrophyte cover. Analysis of similarity indicated that the diatom assemblages of these a priori groups of macrophyte cover were significantly different from one another (i.e., difference between: groups 1 and 3, R statistic = 0.31, P < 0.001; groups 1 and 2, R statistic = 0.049, P < 0.01; groups 3 and 2, R statistic = 0.112, P < 0.001). We then developed an inference model for macrophyte cover from lakes classified as sparse or extensive cover (145 lakes) based on the surface sediment diatom assemblages, and applied this model using the top-bottom paleolimnological approach (i.e., comparison of recent sediments to pre-disturbance sediments). We used the second axis of our correspondence analysis, which significantly divided sparse and extensive macrophyte cover sites, as the independent variable in a logistic regression to predict macrophyte cover as either sparse or extensive. Cross validation, using 48 randomly chosen sites that were excluded from model development, indicated that our model accurately predicts macrophyte cover 79% of the time (r 2 = 0.32, P < 0.001). When applied to the top and bottom sediment samples, our model predicted that 12.5% of natural lakes and 22.4% of reservoirs in the dataset have undergone a ≥30% change in macrophyte cover. For the sites with an inferred change in macrophyte cover, the majority of natural lakes (64.3%) increased in cover, while the majority of reservoirs (87.5%) decreased in macrophyte cover. This study demonstrates that surface sediment diatom assemblages from profundal zones differ in lakes based on their macrophyte cover and that diatoms are useful indicators for quantitatively reconstructing changes in macrophyte cover.  相似文献   

2.
Surficial sediments of three northern Egyptian lakes (Manzala, Burullus and Edku) show differences in diatom assemblages deposited in different sites of these lakes. A total of 172 species and varieties belonging to 58 genera were identified and counted from 62 samples. Of these, 163 diatom taxa were recorded from Manzala Lake sediments, 147 taxa were found in Burullus Lake sediments, and 117 taxa were identified in Edku Lake sediments. The considerable variation in the composition and distribution of the diatom assemblages among these lakes was mainly related to differences in the water quality, salinity, the concentration of nutrients and climatic changes. The planktonic diatom Cyclotella meneghiniana was dominant in the majority of the samples from Manzala Lake, but dominant in only a few samples from the middle parts of Burullus and Edku lakes. The non-planktonic epiphytic taxa Cocconeis placentula and Epithemia sorex were the subdominant species in the surface sediments, especially in shallow and marginal parts of the lakes. Multivariate statistical techniques (hierarchical ascending clustering and canonical correspondence analysis) were used to identify ecological groups of diatoms and to investigate which environmental variables were important in explaining the variation between these groups. Eight ecological groups containing distinctive diatom assemblages reflect current environmental conditions; especially saltwater intrusion in the north and nutrient-rich freshwater in the south.  相似文献   

3.
Studies addressing within-lake variability of fossil chironomid assemblages are very few, and all deal with hydrologically stable temperate lakes where the question of spatial integration mostly relates to the mixing of faunal assemblages associated with shallow, warm-water habitat and those associated with deeper, cold-water habitat. Here we study within-lake variability of surface-sediment chironomid assemblages in the fairly large (∼100–170 km2 since 1983) and shallow (Z max = 5–8 m) fluctuating tropical lake basin of Lake Naivasha, Kenya, and compare the patterns observed with those in two smaller adjacent basins, one similarly shallow (Lake Oloidien, 5.1–5.7 km2, 5–8 m), the other deep and stratified (Crescent Island Crater, 1.9 km2, 14–17 m). Chironomid assemblages were analysed in core-top samples and surface sediments along inshore to offshore transects, and how well individual samples represented the total (basin-wide mean) subfossil assemblage was considered both in terms of taxon richness and taxon percent composition. Within-lake variability of subfossil chironomid concentrations (with generally higher absolute values in nearshore samples) could be explained by effects of sediment winnowing and focusing, whereas between-lake variability reflected their relative susceptibility to wind-driven sediment disturbance or bottom anoxia. In all study lakes, but most significantly in lakes Naivasha and Oloidien, species distribution in the subfossil chironomid assemblages showed a strong nearshore to offshore gradient, which in these shallow lakes, reflects the dominant control of substrate and food quality on species distribution in the living community. Particularly in the larger basins, nearshore samples better represented the total lake assemblage than offshore samples, because the former always contained a component of mud-dwelling species whereas the latter often lacked a component of macrophyte-dwelling species. Our results show that although sedimentation dynamics in the shallow, wind-stressed Lake Naivasha is dominated by frequent resuspension and random sediment redistribution, the near- to offshore gradient in chironomid habitat remains imprinted on subfossil assemblages. We conclude that also in shallow fluctuating lakes, given sufficient size, incomplete pre-burial spatial integration of habitat-specific chironomid assemblages can be exploited for within-lake calibration of environmental gradients.  相似文献   

4.
Surficial sediments from 76 lakes from two western Quebec regions(Abitibi and Haute Mauricie) were sampled to identify the relationships betweendiatoms and environmental variables. Because the two regions containedradically different diatom communities, we then investigated which factors maybe responsible for the large community discrepancies in the two nearbygeographical areas. Standard lake chemistry variables showed little differencesbetween the regions, although epilimnetic light regimes were slightly lower inAbitibi. Nevertheless, lakes of the two regions with similar light regimes andchemistry still showed a clear separation in their diatoms, implying that otherimportant factors are influencing assemblages. We found that the calculatedconcentration of CO2 in the open water can explain some of thediscrepancy in diatom assemblages. A pCCA constrained to the concentration ofCO2 with alkalinity and pH as covariables explained 12.5% of speciesvariance and was significant. Given the lack of a relationship between DOC andCO2, and because the lakes are heavily supersaturated withCO2 in the calibration set, lake-to-lake variations inCO2 concentrations are likely due to groundwater inputs; thepossibility that this environmental variable may be influencing diatomcommunities might allow, in some cases, the reconstruction of historicalchanges in groundwater inputs to lakes. Finally, new calibration models werebuilt in Quebec by using weighted averaging partial least square(WA-PLS) techniques in order to infer pH, CO2, TP, TN, and, DOC fromdiatom assemblages preserved in the surface sediments.  相似文献   

5.
Diatom abundances in the surface sediment samples of 41 mountain lakes in the central Austrian Alps (Niedere Tauern) were related to environmental variables using multi-variate techniques. Canonical correspondence analysis (CCA) revealed that the pH, date of autumn mixing (A mix), mean August water temperature (T Aug), dissolved organic carbon (DOC), and relative water depth (Z rel) made significant contributions to explain the diatom assemblage variation in the lakes of the training set. A weighted averaging partial least square regression and calibration model was used to establish Di-pH (R 2 boot= 0.72, RMSEPboot= 0.131), and a thermistor measurements-based PLS model for A mix (R 2 boot= 0.71, RMSEPboot= 0.006 log10 Julian days). The latter showed a better prediction than T Aug, and was used in terms of climate change. These transfer functions, together with analyses of loss on ignition (LOI), the total carbon/nitrogen (C/N)-ratios, and selected pollen, were applied to an early to mid-Holocene (11.5–4 cal. ky BP) sediment core section from an Austrian Alpine treeline lake on crystalline bedrock. Additionally, passive sample scores in the CCA of the diatom training set were used to show trends in the variables DOC and Z rel. During the early Holocene, diatoms indicative of increased pH, extended warm summers, and low water levels dominated. Between 10.2 and 7.6 cal. ky BP it was followed by diatom assemblages that indicated an increase in lake water depth and an earlier A mix. The multi-proxy data suggest that the A mix decline is the result of a series of snow-rich summer cool and wet climate fluctuations, which were divided by climate warming at ∼9 cal. ky BP. Increased A mix, LOI and DOC, and the correspondent decline in the C/N-ratios, show subsequent climate warming between 7.3 and 6 cal. ky BP. The long-term trend in Di-pH indicates the impact of catchment-related processes during the early-Holocene, that were superimposed by climate.  相似文献   

6.
The relationship between surface sediment diatom assemblages and measured limnological variables in thirty-three coastal Antarctic lakes from the Vestfold Hills was examined by constructing a diatom-water chemistry dataset. Previous analysis of this dataset by canonical correspondence analysis revealed that salinity accounted for a significant amount of the variation in the distribution of the diatom assemblages. Weighted-averaging regression and calibration of this diatom-salinity relationship was used to establish a transfer function for the reconstruction of past lakewater salinity from fossil diatom assemblages. Weighted-averaging regression and calibration with classical deshrinking provided the best model for salinity reconstructions and this was applied to the fossil diatom assemblages from one of the saline lakes in the Vestfold Hills in order to assess its potential for palaeosalinity and palaeoclimate reconstruction.  相似文献   

7.
8.
Diatom analyses were undertaken of sediment cores covering a range of water depths in a small eutrophic lake (Lough Augher, Co. Tyrone, N. Ireland). The significance of between-core variability in diatom relative frequency stratigraphy was assessed by Canonical Correspondence Analysis (CCA) where the ordination axes were constrained to external environmental variables (sediment depth, core location coordinates, water depth, effective fetch, distance-from-shore and distance-from-inflow). After the removal of the effect of sediment age by partialling it out, the resultant first two axes from the partial-CCA were significantly correlated with water depth and distance-from-shore, indicating non-uniform diatom stratigraphies across the lake. Despite this variability, all cores show the same succession of species and, therefore, record the eutrophication of the lake. Diatom-inferred total phosphorus (DI-TP) was inferred for six cores using weighted averaging regression and calibration. Apart from considerable differences of DI-TP in surficial sediment samples, there was good between-core repeatability of DI-TP profiles. These data support the use of DI-TP for establishing background nutrient concentrations for lakes, and associated implications for lake restoration schemes using single cores. Comparisons of DI-TP profiles and total diatom accumulation rate data for the individual cores indicate that diatom production peaked prior to the maximum TP concentrations in the lake.  相似文献   

9.
Surface sediments from 68 small lakes in the Alps and 9 well-dated sediment core samples that cover a gradient of total phosphorus (TP) concentrations of 6 to 520 g TP l-1 were studied for diatom, chrysophyte cyst, cladocera, and chironomid assemblages. Inference models for mean circulation log10 TP were developed for diatoms, chironomids, and benthic cladocera using weighted-averaging partial least squares. After screening for outliers, the final transfer functions have coefficients of determination (r2, as assessed by cross-validation, of 0.79 (diatoms), 0.68 (chironomids), and 0.49 (benthic cladocera). Planktonic cladocera and chrysophytes show very weak relationships to TP and no TP inference models were developed for these biota. Diatoms showed the best relationship with TP, whereas the other biota all have large secondary gradients, suggesting that variables other than TP have a strong influence on their composition and abundance. Comparison with other diatom – TP inference models shows that our model has high predictive power and a low root mean squared error of prediction, as assessed by cross-validation.  相似文献   

10.
过去植被覆盖度重建在长尺度气候模拟和陆地生态环境演变机制研究等方面具有重要意义。基于现代过程的孢粉—植被覆盖度转换函数研究,可为利用孢粉地层数据重建过去植被覆盖度变化提供一条新的途径。以内蒙古高原等地39个湖泊和7个水库中心的表层沉积物孢粉组合与归一化植被指数(NDVI)为研究对象,利用加权平均偏最小二乘法(WA-PLS)、局部加权加权平均法(LWWA)和最佳类比法(MAT)分别建立了孢粉—NDVI转换函数;通过留一交叉检验法、自助法和空间自相关性等检验,筛选出最优模型并应用于古地层孢粉数据,实现定量重建过去植被覆盖度变化。结果表明:MAT和WA-PLS模型均受到空间自相关性的显著影响,而LWWA模型则是建立孢粉—NDVI转换函数的最优模型。与已有类似转换函数表现能力的对比表明,本研究建立的转换函数具有较强的可靠性和应用潜力。基于本研究转换函数重建的内蒙古辉腾锡勒(HTL)区域全新世植被覆盖变化与孢粉重建的植被变化具有很好的一致性,进一步印证了本研究转换函数的较强潜力,均表现出对降水变化的响应更为显著。  相似文献   

11.
We used statistical analyses to determine which subset of 36 environmental variables best explained variations in surface sediment δ13C and δ15N from 50 lakes in western Ireland that span a human-impact gradient. The factors controlling lake sediment δ13C and δ15N depended on whether organics in the lake sediment were mostly derived from the lake catchment (allochthonous) or from productivity within the lake (autochthonous). Lake sediments with a dominantly allochthonous organic source (high C:N ratio sediments) produced δ13C and δ15N measurements similar to values from catchment vegetation. δ13C and δ15N measurements from lake sediments with a dominantly autochthonous organic source (low C:N ratio sediments) were influenced by fractionation in the lake and catchment leading up to assimilation of carbon and nitrogen by lacustrine biota. δ13C values from lake sediment samples in agricultural catchments were more negative than δ13C values from lake sediment samples in non-impacted, bogland catchments. Hypolimnetic oxygen concentrations and methane production had a greater influence on δ13C values than fractionation due to algal productivity. δ15N from lake sediment samples in agricultural catchments were more positive than δ15N in non-impacted bogland catchments. Lower δ15N values from non-impacted lake catchments reflected δ15N values of catchment vegetation, while higher δ15N values in agricultural catchments reflected the high δ15N values of cattle manure and inorganic fertilisers. The influence of changing nitrogen sources and lake/catchment fractionation processes were more important than early diagenesis for lake sediment δ15N values in this dataset. The results from this study suggest a possible influence of bound inorganic nitrogen on the bulk sediment δ15N values. We recommend using a suitable method to control for bound inorganic nitrogen in lake sediments, especially when working with clay-rich sediments. This study confirms the usefulness of δ13C and δ15N from bulk lake sediments, as long as we are mindful of the multiple factors that can influence these values. This study also highlights how stable isotope datasets from lake surface sediments can complement site-specific isotope source/process studies and help identify key processes controlling lake sediment δ13C and δ15N in a study area.  相似文献   

12.
Diatom-based transfer functions for inferring epilimnetic total phosphorus (TP) have been developed from a data set of 33 southeastern Australian water storages. Regular institutional monitoring of these sites has allowed comparison of models developed from TP data covering different time periods. A model based on mean annual TP performs better than models derived from winter maximum TP, spring minimum TP or TP nearest the time of diatom sampling. A mean annual TP model (WA-PLS 2 component) has a jack-knifed diatom-inferred versus measured TP correlation coefficient (r 2 jack) of 0.69 and a root-mean-square-error of prediction (RMSEP) of 0.246 log10g TP l–1, while alternative models have RMSEP > 0.27. Deletion of two samples with uncharacteristic species composition and environmental conditions improved performance of the mean annual TP model (r 2 jack= 0.74; RMSEP = 0.233 log10g TP l–1). Comparison with other published diatom-TP calibration models indicates that this model performs relatively well, with possible contributing factors including the extensive characterisation of TP (with an average 15 determinations making up the annual mean) and the dominance of planktonic diatoms in most sites. Downcore application of the model will allow the reconstruction of reservoir nutrient histories since commissioning, and thus provide a basis for understanding and management of reservoirs.  相似文献   

13.
The fossil diatom records preserved in radiometrically dated sediment cores from four shallow lakes in the Norfolk Broads, UK (Barton Broad, Rollesby Broad, Wroxham Broad and Upton Broad) were analysed. A weighted-averaging partial least squares (WA-PLS) diatom-total phosphorus (TP) transfer function, based on a training set of 152 mostly shallow (maximum depth < 3 m) lakes in northwest Europe, was applied to the full diatom dataset for each core to reconstruct the past TP concentrations of the lakes. Owing to the dominance of non-planktonic Staurosira, Pseudostaurosira and Staurosirella spp. (formerly classified in the genus Fragilaria) throughout the diatom records, the quantitative diatom inferred TP (DI-TP) concentrations did not adequately reflect the changes that occurred in the lakes as indicated by shifts in the other diatom taxa, or as reported in the literature. This was most apparent at Barton Broad and Rollesby Broad, where there was a marked increase in the importance of planktonic taxa associated with highly nutrient-rich waters but no increase in DI-TP. The modern and fossil data were thus square-root transformed to downweight the dominant taxa and the new transfer function was applied to the cores. An improvement was seen only in the reconstruction for Barton Broad. Finally, the Staurosira, Pseudostaurosira and Staurosirella spp. were removed from the modern and fossil diatom data, and the transfer function was re-applied. The trends in DI-TP became less clear, particularly for Upton Broad and Barton Broad, owing to a paucity of data for calibration once these taxa were deleted from the counts data. The problems associated with reconstructing trophic status and determining TP targets for restoration from fossil diatom assemblages in these systems are discussed.  相似文献   

14.
Freshwater midges, consisting of Chironomidae, Chaoboridae and Ceratopogonidae, were assessed as a biological proxy for palaeoclimate in eastern Beringia. The northwest North American training set consists of midge assemblages and data for 17 environmental variables collected from 145 lakes in Alaska, British Columbia, Yukon, Northwest Territories, and the Canadian Arctic Islands. Canonical correspondence analyses (CCA) revealed that mean July air temperature, lake depth, arctic tundra vegetation, alpine tundra vegetation, pH, dissolved organic carbon, lichen woodland vegetation and surface area contributed significantly to explaining midge distribution. Weighted averaging partial least squares (WA-PLS) was used to develop midge inference models for mean July air temperature (r boot2 = 0.818, RMSEP = 1.46°C), and transformed depth (ln (x+1); r boot2 = 0.38, and RMSEP = 0.58).  相似文献   

15.
沉积物中的炭屑在探讨火(灾)的基本特征及其与植被和人类活动关系方面具有重要意义。采用孢粉学提取方法,对柴达木盆地北缘的湖泊表层沉积物进行了微体炭屑的分析工作,用以探讨其空间分布及其与区域植被和现代人类活动之间的关系。结果表明:不同粒径的炭屑浓度均呈现出自西向东升高的趋势,与区域植被覆盖度自西向东增加的趋势一致。其中,小粒径的炭屑以近圆形占多数,表明区域性火事件中燃烧的木本植被占主导地位;粒径较大的长条形炭屑在盆地东部区域占优势,表明地方性火事件向东增加,草本植物燃烧逐渐增多,与自西向东草本植物的增多相吻合,且燃烧由区域向近源转换。此外,该地区人口数量也有自西向东增加的趋势,推测人类生产生活用火对炭屑浓度有很大的贡献量,与大颗粒炭屑浓度向东升高有较好的对应关系。表层沉积物炭屑的空间变化特征不仅很好地反映了相应地区的火的信息,还对区域植被和人类活动等要素具有较好的指示意义。  相似文献   

16.
Diatom assemblages in surficial sediments, sediment cores, sediment traps, and inflowing streams of perennially ice-covered Lake Hoare, South Victorialand, Antarctica were examined to determine the distribution of diatom taxa, and to ascertain if diatom species composition has changed over time. Lake Hoare is a closed-basin lake with an area of 1.8 km2, maximum depth of 34 m, and mean depth of 14 m, although lake level has been rising at a rate of 0.09 m yr-1 in recent decades. The lake has an unusual regime of sediment deposition: coarse grained sediments accumulate on the ice surface and are deposited episodically on the lake bottom. Benthic microbial mats are covered in situ by the coarse episodic deposits, and the new surfaces are recolonized. Ice cover prevents wind-induced mixing, creating the unique depositional environment in which sediment cores record the history of a particular site, rather than a lake-wide integration. Shallow-water (<1 m) diatom assemblages (Stauroneis anceps, Navicula molesta, Diadesmis contenta var. parallela, Navicula peraustralis) were distinct from mid-depth (4–16 m) assemblages (Diadesmis contenta, Luticola muticopsis fo. reducta, Stauroneis anceps, Diadesmis contenta var. parallela, Luticola murrayi) and deep-water (26–31 m) assemblages (Luticola murrayi, Luticola muticopsis fo. reducta, Navicula molesta). Analysis of a sediment core (30 cm long, from 11 m water depth) from Lake Hoare revealed two abrupt changes in diatom assemblages. The upper section of the sediment core contained the greatest biomass of benthic microbial mat, as well as the greatest total abundance and diversity of diatoms. Relative abundances of diatoms in this section are similar to the surficial samples from mid-depths. An intermediate zone contained less organic material and lower densities of diatoms. The bottom section of core contained the least amount of microbial mat and organic material, and the lowest density of diatoms. The dominant process influencing species composition and abundance of diatom assemblages in the benthic microbial mats is episodic deposition of coarse sediment from the ice surface.  相似文献   

17.
This study investigated the distribution of subfossil diatom assemblages in surficial sediments of 100 lakes along steep ecological and climatic gradients in northernmost Sweden (Abisko region, 67.07° N to 68.48° N latitude, 17.67° E to 23.52° E longitude) to develop and cross-validate transfer functions for paleoenvironmental reconstruction. Of 19 environmental variables determined for each site, 15 were included in the statistical analysis. Lake-water pH (8.0%), sedimentary loss-on-ignition (LOI, 5.9% and estimated mean July air temperature (July T, 4.8%) explained the greatest amounts of variation in the distribution of diatom taxa among the 100 lakes. Temperature and pH optima and tolerances were calculated for abundant taxa. Transfer functions, based on WA-PLS (weighted averaging partial least squares), were developed for pH (r2 = 0.77, root-mean-square-error of prediction (RMSEP) = 0.19 pH units, maximum bias = 0.31, as assessed by leave-one-out cross-validation) based on 99 lakes and for July T (r2 = 0.75, RMSEP = 0.96 °C, max. bias = 1.37 °C) based on the full 100 lake set. We subsequently assessed the ability of the diatom transfer functions to estimate lake-water pH and July T using a form of independent cross-validation. To do this, the 100-lake set was divided in two subsets. An 85-lake training-set (based on single limnological measurements) was used to develop transfer functions with similar performance as those based on the full 100 lakes, and a 15-lake test-set (with 2 years of monthly limnological measurements throughout the ice-free seasons) was used to test the transfer functions developed from the 85-lake training-set. Results from the intra-set cross-validation exercise demonstrated that lake-specific prediction errors (RMSEP) for the 15-lake test-set corresponded closely with the median measured values (pH) and the estimations based on spatial interpolations of data from weather stations (July T). The prediction errors associated with diatom inferences were usually within the range of seasonal and interannual variability. Overall, our results confirm that diatoms can provide reliable and robust estimates of lake-water pH and July T, that WA-PLS is a robust calibration method and that long-term environmental data are needed for further improvement of paleolimnological transfer functions.  相似文献   

18.
A suite of surface sediment samples from three Icelandic lakes was analysed for subfossil chironomid head capsules, and a quantitative July air temperature inference model was applied to the data to investigate whether there was significant variability among samples taken from a lake. Ordination and simple regression methods were used to analyse the relationships between environmental and sedimentological variables and the chironomid assemblages and inferred temperature data. Substrate was the most important influence on the chironomid assemblages and inferred temperatures, while water depth at the sampling location had no relationship with the chironomid-inferred temperatures. Within-lake variability of the chironomid assemblages and their inferred temperatures, however, were not significant statistically, suggesting that in lakes of western and northwest Iceland within-lake sampling location has no effect on the data obtained, and therefore on training set samples.  相似文献   

19.
In the Ethiopian highlands, remarkable recovery of vegetation has been achieved using exclosures, protecting vegetation against livestock browsing and firewood harvesting. But these emerging forest resources require tools for sustainable use, implying knowledge on biomass stocks and growth. In this study we developed biomass functions estimating total, stem and branch biomass from diameter at stump height (DSH) and tree height (H) for an 11-year old exclosure in Tigray, Ethiopia. In a systematic grid of 55 plots, DSH and H of all trees and shrubs were recorded. 40 Acacia abyssinica trees were selected for destructive sampling. Allometric relationships using a natural log–log model were established between aboveground biomass, DSH and H. Models with only DSH were found best with R2 between 0.95 and 0.98. The functions were 10 fold cross-validated and R2_cv ranged from 0.94 to 0.97, indicating good model performance. The models were found well in range with those of other seasonal forests in East Africa. Total aboveground biomass was estimated 25.4 ton ha−1 with an annual production of 2.3 ton ha−1, allowing sustainable wood fuel use for 4 persons ha−1. The presented predictive functions help to harmonize between ecological and societal objectives and are as such a first step towards an integrated planning tool for exclosures.  相似文献   

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

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