首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 437 毫秒
1.
Over the last 20 years there has been a surge of interest in paleolimnology and as a result a large accumulation of lake sedimentation records. This emerging archive has allowed us to develop empirical models to describe which variables explain significant variation in sedimentation rates over the past ∼150 years across large spatial scales. We hypothesized that latitude would be a significant explanatory variable of profundal zone lake sedimentation rates across a temperate to polar gradient. We further hypothesized that along a more longitudinally-constrained dataset (i.e. east coast of North America), latitude would explain a greater proportion of the variance. To test these hypotheses, we collated data from 125 natural, average-sized lakes (with surface area <500 km2) by recording authors’ estimates of sedimentation rates (measured as mm/year) or by digitizing recent sediment profiles and calculating sedimentation rates over the past ∼150 years. We found that, at both scales, latitude was the strongest predictor of lake sedimentation rates (full dataset: r 2 = 0.28, P = 0.001, n = 125; east coast dataset: r 2 = 0.58, P < 0.001, n = 43). By conducting a multiple linear regression analysis, we found that 70% of the variance in sedimentation rates from the east coast transect was explained by latitude and elevation alone. This latter model is of sufficient strength that it is a robust predictive tool. Given that climate and land-use strongly co-vary with latitude and that both of these factors have previously been shown to influence lake sedimentation rates, it appears that latitude is a surrogate measure for climate and land-use changes. We also show support for land-use as an important variable influencing sedimentation rates by demonstrating large increases in recent versus Holocene accumulation rates. These results indicate that it is possible to make generalizations about sedimentation rates across broad spatial scales with even limited geographic data.  相似文献   

2.
In the present study, spatio-temporal urban sprawl and land consumption patterns were analysed in seven capital cities located in the Himalayan region during 1972, 1991 and 2015 using multi-temporal satellite images. The study exhibits that capital Himalayan cities experienced rapid growth (830.92%) with high population increase (333.45%) during the observation period (1972–2015). The significant urban growth was observed in the cities of western and middle Himalayan region viz., Srinagar (9.36 km2–142.19 km2), Kathmandu (11.38 km2–92.58 km2) and Dehradun (4.1 km2–50.09 km2) in the higher altitudes due to remarkable increase in the population (0.5–1 million persons) during 1972–2015. On the contrary, Itanagar (7.19 km2), Gangtok (7.09 km2), Shimla (3.04 km2) and Thimphu (2.93 km2) observed less urban growth with moderate to low population growth (i.e., 0.05 to 0.15 million persons). The Shannon entropy based study exhibits that the cities viz., Kathmandu, Gangtok and Itanagar observed comparatively more dispersed urban growth during later period (1991–2015) as compared to the previous period (1972–1991) whereas, the remaining cities observed comparatively less dispersed urban growth during later period. The temporal land consumption pattern exhibits low density urban growth in Srinagar, Dehradun and Kathmandu, as observed with decrease in population density and increasing land consumption during 1972–2015 as compared to other cities, wherein urban densification was evident with increase in population density and decrease in land consumption. The cities in central and western Himalayan region observed high urban growth as compared to cities in eastern Himalayan region. The result shows that the capital cities contributes insignificant proportion (0.5%; 314 km2) of urban area in Himalayan region and accommodating large (ca. 4 million) population during 2015. The study indicates unplanned and haphazard growth in all capital Himalayan cities, leading towards urban densification as well as dispersion in the periphery with varied pattern and intensity. The specific trends and patterns of urban and population growth are governed by geographical as well as socio-economic-political factors at local to regional scale. The high population pressure induced higher risk to the urban residents as well as constrained urban growth over higher vulnerable zones. The study necessitates implementation of suitable urban planning methods considering socio-economic and physico-cultural characteristics of the region.  相似文献   

3.
We present a simple method to derive spatial precipitation (P) and evapotranspiration (ET) for the typical steppe of the Xilin river catchment at 1 km and 8-day resolution during the main vegetation period (23 April to 28 August) of 2006. The hydrological model BROOK90 was parameterised from eddy covariance measurements. The daily model input data, precipitation, minimum (Tamin) and maximum air temperature (Tamax), were derived by manipulating MODIS leaf area index (LAI) and surface temperature data. P was estimated based on a linear regression of P measured at several sites against the mean gain of the MODIS LAI of surrounding 3 × 3 pixels areas (R2 = 0.76). Tamin and Tamax were derived using a relationship between measured Tamin and Tamax and MODIS surface temperatures (R2 = 0.92 and R2 = 0.88, respectively). The mean precipitation was 145 mm; it varied between 52 mm in the north-western region and 239 mm in the eastern region. In spring, the modelled ET was low (<0.8 mm d−1); evaporation dominated over transpiration and spatial differences were small. At the end of June, the mean ET reached its maximum (2 mm d−1) and spatial differences were pronounced. From July on, transpiration dominated over declining evaporation, and spatial differences decreased in August.  相似文献   

4.
Mountaintop removal mining is a dominant driver of land use/land cover changes in the Appalachian Region of the eastern United States and is expected to increase in scale in the coming decades. While several studies quantify land use/land cover changes attributed to traditional surface mining and at regional scales, no studies we are aware of focus specifically on mountaintop removal/valley fill mining practices at the watershed scale. Further, despite its scale and extent, its impact on runoff, particularly at larger spatial scales (103 km2), is poorly understood due to the complex relationships between climate, land use, and hydrology. To explore the impacts of this practice at broader scales, we estimated land use/land cover changes using Landsat 5 TM imagery over five periods between 1994 and 2010; used a simple rainfall–runoff model to estimate hydrologic response time; and conducted non-parametric trend analyses on annual hydrologic metrics (streamflow, Q/P, response time) for the Big Coal River watershed located in the southern West Virginia coalfields. No statistically significant trends were detected in any of the timeseries. The lack of detectable trends and correlations between land use changes and hydrology at the basin scale are not entirely unexpected due to the history and mosaic of land cover changes that span timescales larger than our study period. Further interannual variation likely overwhelms our ability to detect potential changes using monotonic trend analysis at the annual time scale, particularly in light of strong streamflow seasonality. Future studies therefore should include different methods of change detection applied to different timescales to more appropriately account seasonal and interannual variations. Until the significance of this practice on water resources (quality and quality) are understood, efforts to reduce the environmental problems associated with mountaintop mining will be difficult to achieve.  相似文献   

5.
Floods are a frequently occurring calamity in deltaic Bangladesh. This paper aims to assess the temporal expansion of waterbodies during flooding using geospatial techniques. Several water indices were applied to classify the satellite images at various temporal scales. Among them, the Normalized Difference Water Index (NDWI) showed the highest correlation (r = 0.831; where p = 0.01) with rainfall data. Specifically, the NDWI results showed that perennial waterbodies measured 37 km2 and 60 km2 in Sunamganj District in 2017 and 2019, respectively. The area of waterbodies notably increased 52-fold from March to April (37 km2 to 1958 km2) during the pre-monsoon flash flood of 2017. During the July 2019 monsoon flood, waterbodies started to extend after May and flooded 2784 km2 in area. NDVI analysis showed that in 2019, floodwater submerged 361.7 km2 of vegetation cover. At the same time, the Surma River's flooding resulted in a 73.9 per cent inundation of the total area of the Sunamganj District. We hope that this study will provide better understanding of the varying nature of floods that occur in the low lying bowl shaped Haor region which will in turn assist the government with flood mitigation.  相似文献   

6.
The aim of this study was to evaluate four metrics to define the spatially variable (regionalised) hillslope sediment delivery ratio (HSDR). A catchment model that accounted for gully and streambank erosion and floodplain deposition was used to isolate the effects of hillslope gross erosion and hillslope delivery from other landscape processes. The analysis was carried out at the subcatchment (~ 40 km2) and the cell scale (400 m2) in the Avon-Richardson catchment (3300 km2), south-east Australia. The four landscape metrics selected for the study were based on sediment travel time, sediment transport capacity, flux connectivity, and residence time. Model configurations with spatially-constant or regionalised HSDR were calibrated against sediment yield measured at five gauging stations. The impact of using regionalised HSDR was evaluated in terms of improved model performance against measured sediment yields in a nested monitoring network, the complexity and data requirements of the metric, and the resulting spatial relationship between hillslope erosion and landscape factors in the catchment and along hillslope transects. The introduction of a regionalised HSDR generally improved model predictions of specific sediment yields at the subcatchment scale, increasing model efficiency from 0.48 to > 0.6 in the best cases. However, the introduction of regionalised HSDR metrics at the cell scale did not improve model performance. The flux connectivity was the most promising metric because it showed the largest improvement in predicting specific sediment yields, was easy to implement, was scale-independent and its formulation was consistent with sedimentological connectivity concepts. These properties make the flux connectivity metric preferable for applications to catchments where climatic conditions can be considered homogeneous, i.e. in small-medium sized basins (up to approximately 3000 km2 for Australian conditions, with the Avon-Richardson catchment being at the upper boundary). The residence time metric improved model assessment of sediment yields and enabled accounting for climatic variability on sediment delivery, but at the cost of greater complexity and data requirements; this metric might be more suitable for application in catchments with important climatic gradients, i.e. large basins and at the regional scale. The application of a regionalised HSDR metric did not increase data or computational requirements substantially, and is recommended to improve assessment of hillslope erosion in empirical, semi-lumped erosion modelling applications. However, more research is needed to assess the quality of spatial patterns of erosion depicted by the different landscape metrics.  相似文献   

7.
There is a growing requirement for techniques to assess land susceptibility to wind erosion, i.e. land erodibility, over large geographic areas (>104 km2). This requirement stems from a lack of wind erosion research between the field (101 km2) and regional scales, and a need to evaluate the performance of spatially explicit wind erosion models across these scales. This paper addresses this issue by presenting a methodology for monitoring land erodibility at the landscape scale (103 km2). First, we define criteria suitable for evaluating land erodibility based on empirical relationships between soil texture, vegetation cover, geomorphology, and wind erosion. The criteria were used to visually assess land erodibility over long distances (103 km) using vehicle-based transects run through the rangelands of western Queensland, Australia. Application of the data for testing the performance of a spatially explicit land erodibility model (AUSLEM) is then demonstrated by comparing the visual assessments of land erodibility with the model output. The model performed best in the west of the study area in the open rangelands. In regions with higher woody shrub and tree cover the model performance decreases. This highlights the need for research to better parameterise controls on erodibility in semi-arid landscapes consisting of forested and rangeland mosaics.  相似文献   

8.
Jose Luis Antinao  John Gosse   《Geomorphology》2009,104(3-4):117-133
The distribution and age of large (> 0.1 km2) Pliocene to recent rockslides in the Chilean Cordillera Principal (32–34.5 S), the Southern Central Andes, has been analyzed to determine the rockslide triggering mechanisms and impact on regional landscape evolution. Most of the rockslides appear in the western Cordillera Principal and cluster along major geological structures. Variographic analyses show spatial correlation between rockslides, geological structures and shallow seismicity. A relative chronosequence was calibrated with existing 14C and 40Ar/39Ar dates and new cosmogenic nuclide exposure ages for selected rockslides. Rockslide-induced sediment yield was estimated with empirical relations for rockslide area distributions. Throughout the Quaternary, rockslides have delivered sediment to streams at rates equivalent to denudation rates of 0.10 ±0.06 mm a− 1, while estimates using short term (20 a) seismicity records are 0.3− 0.2+ 0.6 mm a− 1. The estimates of sediment transfer and the spatial distribution of rockslides reflect a landscape in which tectonic and geological controls on denudation are more significant than climate.  相似文献   

9.
We studied the relationship between primary productivity and species richness of small mammals at both large and small spatial scales in the arid and semi-arid grasslands of north China. The productivity (x)–species richness (y) pattern at a large spatial scale can be described by a unimodal quadratic regression curve (y=7·41+0·1*x−0·0003*x2,p =0·008, r2=0·788). At a small spatial scale, however, neither linear nor quadratic regression fit the data for 1980 and 1994 (p>0·25). Primary productivity may not be an appropriate predictor of the species richness of small mammals at a small spatial scale. We conclude, therefore, that the primary productivity–species richness pattern of small mammals may be scale-dependent in the arid and semi-arid grasslands of north China. Landscape complexity should be considered in future studies of productivity–richness relationships.  相似文献   

10.
Annual above-ground net primary production (ANPP), evapotranspiration (ET) and water use efficiency (WUE) of rangeland have the potential to provide an objective basis for establishing pricing for ecosystem services. To provide estimates of ANPP, we surveyed the biomass, estimated ET and prepared a water use efficiency for dwarf shrublands and arid savanna in the Riemvasmaak Rural Area, Northern Cape, South Africa. The annual production fraction was surveyed in 33 MODIS 1 km2 pixels and the results regressed against the MODIS fPAR product. This regression model was used to predict the standing green biomass (kg DM ha−1) for 2009 (dry year). Using an approach which combines potential evapotranspiration (ET0) and the MODIS fPAR product, we estimated actual evapotranspiration (ETa). These two models (greening standing biomass and ETa) were used to calculate the annual WUE for 2009. WUE was 1.6 kg DM mm−1 ha−1 yr−1. This value may be used to provide an estimate of ANPP in the absence of direct measurements of biomass and to provide a comparison of the water use efficiency of this rangeland with other rangeland types.  相似文献   

11.
We evaluated nest site selection at two spatial scales (microsite, territory) and reproductive success of Western Burrowing Owls (Athene cunicularia hypugaea) at three spatial scales (microsite, territory, landscape) in the eastern Mojave Desert. We used binary logistic regression within an information-theoretic approach to assess factors influencing nest site choice and nesting success. Microsite-scale variables favored by owls included burrows excavated by desert tortoise (Gopherus agassizii), burrows with a large mound of excavated soil at the entrance, and a greater number of satellite burrows within 5 m of the nest burrow. At the territory scale, owls preferred patches with greater cover of creosote bush (Larrea tridentata) within 50 m of the nest burrow. An interaction between the presence or absence of a calcic soil horizon layer over the top of the burrow (microsite) and the number of burrows within 50 m (territory) influenced nest site choice. Nesting success was influenced by a greater number of burrows within 5 m of the nest burrow. Total cool season precipitation was a predictor of nesting success at the landscape scale. Conservation strategies can rely on management of habitat for favored and productive nesting sites for this declining species.  相似文献   

12.
Mapping ecosystem services (ES) over large scales is important for environmental monitoring but is often prohibitively expensive and difficult. We test a hybrid, low-cost method of mapping ES indicators over large scales in Pará State, Brazil. Four ES indicators (vegetation carbon stocks, biodiversity index, soil chemical quality index and rates of water infiltration into soil) were measured in the field and then summarized spatially for regional land-cover classes derived from satellite imagery. The regionally mapped ES values correlated strongly with independent and local measures of ES. For example, regional estimates of the vegetation carbon stocks are strongly correlated with actual measures derived from field samples and validation data (significant anova test – p-value = 4.51e?9) and differed on average by only 20 Mg/ha from the field data. Our spatially-nested approach provides reliable and accurate maps of ES at both local and regional scales. Local maps account for the specificities of an area while regional maps provide an accurate generalization of an ES’ state. Such up-scaling methods infuse large-scale ES maps with localized data and enable the estimation of uncertainty of at regional scales. Our approach is first step towards the spatial characterization of ES at large and potentially global scales.  相似文献   

13.
Although tropical deforestation bears a close relationship with climate change, its exact contribution to climate warming and its threshold of exerting a noticeable influence remain unknown. This study attempts to bridge this knowledge gap by analyzing deforestation data of Heilongjiang Province, China in relation to climate data. It is found that forest cover was reduced from 238,335 km2 in 1958 to 216,009 km2 in 1980, and further to 207,629 km2 in 2000. During this period the provincial annual temperature rose by 1.68 °C, against the nation-wide warming of 0.99 °C during the same period. At the provincial level the observed deforestation caused a warming in the vicinity of 0.69 °C. This warming does not bear any definite relationship with latitude and elevation. At the local scale, deforestation is related inversely to the rise in decadal temperature in the form of ΔT = −0.013ΔF + 0.4114 (R2 = 0.30). There is a positive relationship between the accuracy (R2 value) of predicting climate warming from deforestation and its severity. The critical threshold for deforestation to exert a noticeable impact on climate warming (e.g., R2 = 50%) appears to be 5 km2. The amount of forest cover at the beginning of a period can inhibit temperature rise, but its exact effect on climate warming is difficult to quantify.  相似文献   

14.
《Polar Science》2014,8(2):166-182
The larch forests on the permafrost in northeastern Mongolia are located at the southern limit of the Siberian taiga forest, which is one of the key regions for evaluating climate change effects and responses of the forest to climate change. We conducted long-term monitoring of seasonal and interannual variations in hydrometeorological elements, energy, and carbon exchange in a larch forest (48°15′24′′N, 106°51′3′′E, altitude: 1338 m) in northeastern Mongolia from 2010 to 2012. The annual air temperature and precipitation ranged from −0.13 °C to −1.2 °C and from 230 mm to 317 mm. The permafrost was found at a depth of 3 m. The dominant component of the energy budget was the sensible heat flux (H) from October to May (H/available energy [Ra] = 0.46; latent heat flux [LE]/Ra = 0.15), while it was the LE from June to September (H/Ra = 0.28, LE/Ra = 0.52). The annual net ecosystem exchange (NEE), gross primary production (GPP), and ecosystem respiration (RE) were −131 to −257 gC m−2 y−1, 681–703 gC m−2 y−1, and 423–571 gC m−2 y−1, respectively. There was a remarkable response of LE and NEE to both vapor pressure deficit and surface soil water content.  相似文献   

15.
This study aimed at characterizing land cover dynamics for four decades in Eastern Mau forest and Lake Nakuru basin, Kenya. The specific objectives were to: (i) identify and map the major land cover types in 1973, 1985, 2000 and 2011; (ii) detect and determine the magnitude, rates and nature of the land cover changes that had occurred between these dates, and; (iii) establish the spatial and temporal distribution of these changes. Land cover types were discriminated through partitioning, hybrid classification and spatial reclassification of multi-temporal Landsat imagery. The land cover products were then validated and overlaid in post-classification comparison to detect the changes between 1973 and 2011. The accuracies of the land cover maps for 1973, 1985, 2000 and 2011 were 88%, 95%, 80% and 89% respectively. Six land cover classes, namely forests-shrublands, grasslands, croplands, built-up lands, bare lands and water bodies, were mapped. Forests-shrublands dominated in 1973, 1985 and 2000 covering about 1067 km2, 893 km2 and 797 km2 respectively, but were surpassed by croplands (953 km2) in 2011. Bare lands occupied the least area that varied between 2 km2 and 7 km2 during this period. Overall, forests-shrublands and grasslands decreased by 428 km2 and 258 km2 at the annual average rates of 1% each, whereas croplands and built-up lands expanded by 660 km2 and 24 km2 at the annual rates of 6% and 16% respectively. The key hotspots of these changes were distributed in all directions of the study area, but at different times. Therefore, policies that integrate restoration and conservation of natural ecosystems with enhancement of agricultural productivity are strongly recommended. This will ensure environmental sustainability and socio-economic well-being in the area. Future research needs to assess the impacts of the land cover changes on ecosystem services and to project the future patterns of land cover changes.  相似文献   

16.
A large spatial variability in sediment yield was observed from small streams in the Ecuadorian Andes. The objective of this study was to analyze the environmental factors controlling these variations in sediment yield in the Paute basin, Ecuador. Sediment yield data were calculated based on sediment volumes accumulated behind checkdams for 37 small catchments. Mean annual specific sediment yield (SSY) shows a large spatial variability and ranges between 26 and 15,100 Mg km− 2 year− 1. Mean vegetation cover (C, fraction) in the catchment, i.e. the plant cover at or near the surface, exerts a first order control on sediment yield. The fractional vegetation cover alone explains 57% of the observed variance in ln(SSY). The negative exponential relation (SSY = a × e−b C) which was found between vegetation cover and sediment yield at the catchment scale (103–109 m2), is very similar to the equations derived from splash, interrill and rill erosion experiments at the plot scale (1–103 m2). This affirms the general character of an exponential decrease of sediment yield with increasing vegetation cover at a wide range of spatial scales, provided the distribution of cover can be considered to be essentially random. Lithology also significantly affects the sediment yield, and explains an additional 23% of the observed variance in ln(SSY). Based on these two catchment parameters, a multiple regression model was built. This empirical regression model already explains more than 75% of the total variance in the mean annual sediment yield. These results highlight the large potential of revegetation programs for controlling sediment yield. They show that a slight increase in the overall fractional vegetation cover of degraded land is likely to have a large effect on sediment production and delivery. Moreover, they point to the importance of detailed surface vegetation data for predicting and modeling sediment production rates.  相似文献   

17.
Mean tree biomass and soil carbon (C) densities for 39 map sheet grids (1° lat. × 1.5° long.) covering the Acacia woodland savannah region of Sudan (10–16° N; 21–36° E) are presented. Data from the National Forest Inventory of Sudan, Harmonized World Soil Database and FAO Local Climate Estimator were used to calculate C densities, mean annual precipitation (MAP) and mean annual temperature (MAT). Above-ground biomass C and soil organic carbon (SOC, 1 m) densities averaged 112 and 5453 g C m−2, respectively. Below-ground biomass C densities, estimated using root shoot ratios, averaged 33 g C m−2. Biomass C densities and MAP increased southwards across the region while SOC densities were lowest in the centre of the region and increased westwards and eastwards. Both above-ground biomass C and SOC densities were significantly (p < 0.05) correlated with MAP (rs = 0.84 and rs = 0.34, respectively) but showed non-significant correlations with MAT (rs = −0.22 and rs = 0.24, respectively). SOC densities were significantly correlated with biomass C densities (rs = 0.34). The results indicated substantial under stocking of trees and depletion of SOC, and potential for C sequestration. Up-to-date regional and integrated soil and forest inventories are required for planning improved land-use management and restoration.  相似文献   

18.
Ecological patterns and processes in dune ecosystems have been a research focus in recent years, however information on how dune stabilization influences vegetation and soil at different spatial scales...  相似文献   

19.
Land degradation and global warming are currently highly active research topics. Land degradation can both change land cover and surface climate and significantly influence atmospheric circulation. Researches have verified that carbon dioxide (CO2) and methane (CH4) are major greenhouse gases (GHG) in the atmosphere and are directly affected by human activity. However, to date, there is no research on the spatial distribution of GHG concentrations and also no research on how land degradations affect GHG concentrations in arid and semi-arid regions. In this study, we used GHG data from the ENVIronment SATellite (ENVISAT) and the Greenhouse gases Observing Satellite (GOSAT), the Normalized Difference Vegetation Index (NDVI) and Land Surface Temperature (LST) data from the MODerate resolution Imaging Spectroradiometer (MODIS) and precipitation data from ground stations to analyze the way land degradation affects GHG concentrations in northern China and Mongolia, which exhibit the most serious land degradation process in East Asia. Our research revealed that the CO2 and CH4 concentrations (XCO2 and XCH4) increased from 2003 to 2009 and then decreased into 2011. We used geostatistics to predict and simulate the spatial distribution of XCO2 and XCH4 and found that the distribution of XCO2 displays a seasonal trend and is primarily affected by plant photosynthesis, soil respiration and precipitation. As the distribution of XCH4 is mainly affected by the sources' distribution, microbial processes, LST and submarine hydrate, the CH4 concentration presents no obvious seasonal changes and the high XCH4 values are primarily found in northeast and southeast China. Land degradation increases the concentration of GHG: the correlation coefficient between NDVI and XCO2 is R2 = 0.76 (P < 0.01) and the value between NDVI and XCH4 is R2 = 0.75 (P < 0.01).  相似文献   

20.
Species-specific allometric models were developed to predict aboveground biomass (AGB) of eight woody species in the Borana rangelands, Ethiopia. The 23 equations developed (8 species; three biomass components: total aboveground, stem and branches) fit the data well to predict total AGB and by components for each of the species (r2 > 0.70; p < 0.001). The AGB of tree shaped species (e.g., Acacia bussei and Acacia etabaica) were significantly predicted from a single predictor (circumference of the stem at ankle height), with a high coefficient of determination (r2 > 0.95; p < 0.001). In contrast, the AGB of bushy shrubs (e.g., Acacia oerfota) was more effectively predicted by using the canopy volume (r2 = 0.84; p < 0.001). Shrubs with a tall stem and an umbrella-like canopy structure (e.g., Acacia mellifera) were most accurately predicted by a combination of both circumference of the stem at ankle height and canopy volume (r2 = 0.95; p < 0.001). Hence, our species-specific allometric models could accurately estimate their woody aboveground biomass in a semi-arid savanna ecosystem of southern Ethiopia. These equations will help in future carbon-trade discussions in times of climate change and CO2 emission concerns and mitigation strategies.  相似文献   

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

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