首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
ABSTRACT: Three instruments commonly used to measure stream canopy cover were evaluated: the clinometer, a modified spherical convex densiometer Model A, and a hemispherical image system. The hemispherical image system was also used to model shade. At each of five locations throughout Oregon, canopy cover above a stream was measured from the center of the stream along a series of transects. Vegetation along the reaches sampled ranged from a densely vegetated coastal forest stand to a sparsely vegetated eastern Oregon meadow. When techniques were compared within each site, canopy cover measured with the clinometer was similar to that measured with hemispherical imagery. The densiometer measurements, however, were typically lower than, though still highly correlated with, those derived from hemispherical images and from the clinometer. Because of site‐specific factors, the differences between all methods were significant in the Willamette and John Day Provinces. Canopy cover and shade were not equivalent for the streams studied.  相似文献   

2.
We describe the development of a bird integrity index (BII) that uses bird assemblage information to assess human impacts on 13 stream reaches in the Willamette Valley, Oregon, USA. We used bird survey data to test 62 candidate metrics representing aspects of bird taxonomic richness, tolerance or intolerance to human disturbance, dietary preferences, foraging techniques, and nesting strategies that were affected positively or negatively by human activities. We evaluated the metric responsiveness by plotting each one against a measure of site disturbance that included aspects of land use/land cover, road density, riparian cover, and stream channel and substrate conditions. In addition, we eliminated imprecise and highly correlated (redundant) metrics, leaving 13 metrics for the final index. Individual metric scores ranged continuously from 0 to 10, and index scores were weighted to range from 0 to 100. Scores were calibrated using historical species information to set expectations for the number of species expected under minimally disturbed conditions. Site scores varied from 82 for the least disturbed stream reach to 8.5 for an urban site. We compared the bird integrity index site scores with the performance of other measures of biotic response developed during this study: a fish index of biointegrity (IBI) and two benthic macroinvertebrate metrics. The three assemblages agreed on the general level of disturbance; however, individual sites scored differently depending on specific indicator response to in-stream or riparian conditions. The bird integrity index appears to be a useful management and monitoring tool for assessing riparian integrity and communicating the results to the public. Used together with aquatic indicator response and watershed data, bird assemblage information contributes to a more complete picture of stream condition.  相似文献   

3.
The Bird Integrity Index (BII) presented here uses bird assemblage information to assess human impacts to 28 stream reaches in the Blue Mountains of eastern Oregon. Eighty-one candidate metrics were extracted from bird survey data for testing. The metrics represented aspects of bird taxonomic richness, tolerance or intolerance to human disturbance, dietary preferences, foraging techniques, and nesting strategies that were expected to be positively or negatively affected by human activities in the region. To evaluate the responsiveness of each metric, it was plotted against an index of reach and watershed disturbance that included attributes of land use/land cover, road density, riparian cover, mining impacts, and percent area in clearcut and partial-cut logging. Nine of the 81 candidate bird metrics remained after eliminating unresponsive and highly correlated metrics. Individual metric scores ranged from 0 to 10, and BII scores varied between 0 and 100. BII scores varied from 78.6 for a minimally disturbed, reference stream reach to 30.4 for the most highly disturbed stream reach. The BII responded clearly to varying riparian conditions and to the cumulative effects of disturbances, such as logging, grazing, and mining, which are common in the mountains of eastern Oregon. This BII for eastern Oregon was compared to an earlier BII developed for the agricultural and urban disturbance regime of the Willamette Valley in western Oregon. The BII presented here was sensitive enough to distinguish differences in condition among stream riparian zones with disturbances that were not as obvious or irreversible as those in the agricultural/urban conditions of western Oregon.  相似文献   

4.
Landscape characteristics of a watershed are important variables that influence surface water quality. Understanding the relationship between these variables and surface water quality is critical in predicting pollution potential and developing watershed management practices to eliminate or reduce pollution risk. To understand the impacts of landscape characteristics on water quality in mine waste-located watersheds, we conducted a case study in the Tri-State Mining District which is located in the conjunction of three states (Missouri, Kansas and Oklahoma). Severe heavy metal pollution exists in that area resulting from historical mining activities. We characterized land use/land cover over the last three decades by classifying historical multi-temporal Landsat imagery. Landscape metrics such as proportion, edge density and contagion were calculated based on the classified imagery. In-stream water quality data over three decades were collected, including lead, zinc, iron, cadmium, aluminum and conductivity which were used as key water quality indicators. Statistical analyses were performed to quantify the relationship between landscape metrics and surface water quality. Results showed that landscape characteristics in mine waste-located watersheds could account for as much as 77% of the variation of water quality indicators. A single landscape metric alone, such as proportion of mine waste area, could be used to predict surface water quality; but its predicting power is limited, usually accounting for less than 60% of the variance of water quality indicators.  相似文献   

5.
Prior research has demonstrated the utility of metrics based on spatial velocity gradients to characterize and describe stream habitat, with higher gradients generally indicative of higher levels of physical heterogeneity and thus habitat quality. However, detailed velocity data needed to compute these metrics are difficult to obtain. Acoustic Doppler current profilers (ADCP) may be used to rapidly collect detailed representations of river velocity fields. Herein we demonstrate use of ADCP to obtain ecologically relevant data and compute associated metrics. Data were collected from four reaches of the Little Tallahatchie River in northern Mississippi. Sampled reaches were selected to observe velocity regimes associated with three distinctly different conditions: downstream from a major flow obstruction (a low weir), downstream from the apices of each of two bends, and within an extremely long, straight reach created by channelization. Three-dimensional velocity data sets from each site were used to compute metrics of habitat quality proposed by others. A habitat metric based on the presence of rotational flow in the vertical plane proved to be the best discriminator among conditions within the sampled reaches. Two of four habitat quality metrics computed from these measured velocities were greatest for the sharpest meander bend. ADCP hold great potential for study of riverine physical aquatic habitats, particularly at the reach scale. Additional work is needed to develop generally applicable field protocols and data reduction tools. Specifically, guidelines for ADCP settings and configuration appropriate for a range of riverine site conditions must be developed. Advances in instrumentation are needed to allow collection of information in closer proximity to the free surface and solid boundaries.  相似文献   

6.
The objective of this study was to evaluate which macroinvertebrate and deposited sediment metrics are best for determining effects of excessive sedimentation on stream integrity. Fifteen instream sediment metrics, with the strongest relationship to land cover, were compared to riffle macroinvertebrate metrics in streams ranging across a gradient of land disturbance. Six deposited sediment metrics were strongly related to the relative abundance of Ephemeroptera, Plecoptera and Trichoptera and six were strongly related to the modified family biotic index (MFBI). Few functional feeding groups and habit groups were significantly related to deposited sediment, and this may be related to the focus on riffle, rather than reach-wide macroinvertebrates, as reach-wide sediment metrics were more closely related to human land use. Our results suggest that the coarse-level deposited sediment metric, visual estimate of fines, and the coarse-level biological index, MFBI, may be useful in biomonitoring efforts aimed at determining the impact of anthropogenic sedimentation on stream biotic integrity.  相似文献   

7.
ABSTRACT: An assessment of physical conditions in urban streams of the Puget Sound region, coupled with spatially explicit watershed characterizations, demonstrates the importance of spatial scale, drainage network connectivity, and longitudinal downstream trends when considering the effects of urbanization on streams. A rapid stream assessment technique and a multimetric index were used to describe the physical conditions of multiple reaches in four watersheds. Watersheds were characterized using geographic information system (GIS) derived landscape metrics that represent the magnitude of urbanization at three spatial scales and the connectivity of urban land. Physical conditions, as measured by the physical stream conditions index (PSCI), were best explained for the watersheds by two landscape metrics: quantity of intense and grassy urban land in the subwatershed and quantity of intense and grassy urban land within 500 m of the site (R2= 0.52, p > 0.0005). A multiple regression of PSCI with these metrics and an additional connectivity metric (proximity of a road crossing) provided the best model for the three urban watersheds (R2= 0.41, p > 0.0005). Analyses of longitudinal trends in PSCI within the three urban watersheds showed that conditions improved when a stream flowed through an intact riparian buffer with forest or wetland vegetation and without road crossings. Results demonstrate that information on spatial scale and patterns of urbanization is essential to understanding and successfully managing urban streams.  相似文献   

8.
The Macroinvertebrate Biotic Integrity Index (MBII) was developed from data collected at 574 wadeable stream reaches in the Mid-Atlantic Highlands region (MAHR) by the U.S. Environmental Protection Agency's (USEPA) Environmental Monitoring and Assessment Program (EMAP). Over 100 candidate metrics were evaluated for range, precision, responsiveness to various disturbances, relationship to catchment area, and redundancy. Seven metrics were selected, representing taxa richness (Ephemeroptera richness, Plecoptera richness, Trichoptera richness), assemblage composition (percent non-insect individuals, percent 5 dominant taxa), pollution tolerance [Macroinvertebrate Tolerance Index (MTI)], and one functional feeding group (collector-filterer richness). We scored metrics and summed them, then ranked the resulting index through use of independently evaluated reference stream reaches. Although sites were classified into lowland and upland ecoregional groups, we did not need to develop separate scoring criteria for each ecoregional group. We were able to use the same metrics for pool and riffle composite samples, but we had to score them differently. Using the EMAP probability design, we inferred the results, with known confidence bounds, to the 167,797 kilometers of wadeable streams in the Mid-Atlantic Highlands. We classified 17% of the target stream length in the MAHR as good, 57% as fair, and 26% as poor. Pool-dominated reaches were relatively rare in the MAHR, and the usefulness of the MBII was more difficult to assess in these reaches. The process used for developing the MBII is widely applicable and resulted in an index effective in evaluating region-wide conditions and distinguishing good and impaired reaches among both upland and lowland streams dominated by riffle habitat.  相似文献   

9.
Abstract: We compared benthic macroinvertebrate assemblage structure at snag and riffle habitats in 43 Wisconsin streams across a range of watershed urbanization using a variety of stream quality metrics. Discriminant analysis indicated that dominant taxa at riffles and snags differed; Hydropsychid caddisflies (Hydropsyche betteni and Cheumatopsyche spp.) and elmid beetles (Optioservus spp. and Stenemlis spp.) typified riffles, whereas isopods (Asellus intermedius) and amphipods (Hyalella azteca and Gammarus pseudolimnaeus) predominated in snags. Analysis of covariance indicated that samples from snag and riffle habitats differed significantly in their response to the urbanization gradient for the Hilsenhoff biotic index (BI), Shannon’s diversity index, and percent of filterers, shredders, and pollution intolerant Ephemeroptera, Plecoptera, and Trichoptera (EPT) at each stream site (p ≤ 0.10). These differences suggest that although macroinvertebrate assemblages present in either habitat type are sensitive to detecting the effects of urbanization, metrics derived from different habitats should not be intermixed when assessing stream quality through biomonitoring. This can be a limitation to resource managers who wish to compare water quality among streams where the same habitat type is not available at all stream locations, or where a specific habitat type (i.e., a riffle) is required to determine a metric value (i.e., BI). To account for differences in stream quality at sites lacking riffle habitat, snag‐derived metric values can be adjusted based on those obtained from riffles that have been exposed to the same level of urbanization. Comparison of nonlinear regression equations that related stream quality metric values from the two habitat types to percent watershed urbanization indicated that snag habitats had on average 30.2 fewer percent EPT individuals, a lower diversity index value than riffles, and a BI value of 0.29 greater than riffles.  相似文献   

10.
ABSTRACT: Land cover and land use change have long been known to influence the chemical, physical, and biological characteristics of streams. This study makes use of land cover maps derived from fine resolution satellite imagery and an extensive stream quality dataset to determine the relationship between small watershed health rankings and land cover composition and configuration. Landscape metrics were derived from digital impervious surface area (ISA), tree cover (percent), and agricultural crop maps within Montgomery County, Maryland. Watershed rankings were developed by state and county collaborators (MD‐DNR and MCDEP) using extensive biological and chemical measurements. In stepwise logistic regression models the factors accounting for the most variation in stream health ranking were the percent ISA, followed by the percent of tree cover. Riparian buffer zone tree cover was also a significant predictor. Of the metrics that considered the spatial configuration of the landscape, a contagion index and the percent of ISA in the flow path from the ISA to the stream were also found to be significant predictors of stream health. Despite limited ability to characterize landscape configuration or narrow riparian buffer zone vegetation with coarser resolution imagery (from Landsat), model results were not significantly different from those based on the use of fine‐resolution ISA information, suggesting that broader area applications of the approach are possible. The results indicate that management practices designed to improve stream water quality should focus on the amount of ISA and tree cover in both the watershed and within the buffer zone.  相似文献   

11.
Riparian areas link aquatic and terrestrial habitats, supporting species-rich bird communities, which integrate both terrestrial and aquatic processes. For this reason, inclusion of riparian birds in stream bioassessment could add to the information currently provided by existing programs that monitor aquatic organisms. To assess if bird community metrics could indicate stream conditions, we sampled breeding birds in the riparian zone of 37 reaches in 5 streams draining watersheds representing a gradient of agricultural intensity in central Italy. As a more direct indicator of water quality, stream macroinvertebrates were also sampled for computation of the Italian Extended Biotic Index (IBE). An anthropogenic index was calculated within 1 km of sampled reaches based on satellite-derived land-use classifications. Predictive models of macroinvertebrate integrity based on land-use and avian metrics were compared using an information-theoretic approach (AIC). We also determined if stream quality related to the detection of riverine species. Apparent bird species diversity and richness peaked at intermediate levels of land-use modification, but increased with IBE values. Water quality did not relate to the detection of riverine species as a guild, but two species, the dipper Cinclus cinclus and the grey wagtail Motacilla cinerea, were only observed in reaches with the highest IBE values. Small-bodied insectivorous birds and arboreal species were detected more often in reaches with better water quality and in less modified landscapes. In contrast, larger and granivorous species were more common in disturbed reaches. According to the information-theoretic approach, the best model for predicting water quality included the anthropogenic index, bird species diversity, and an index summarizing the trophic structure of the bird community. We conclude that, in combination with landscape-level information, the diversity and trophic structure of riparian bird communities could serve as a rapid indicator of stream-dwelling macroinvertebrates and, therefore, degradation of in-stream biotic integrity.  相似文献   

12.
Riparian zones in semi‐arid, mountainous regions provide a disproportionate amount of the available wildlife habitat and ecosystem services. Despite their importance, there is little guidance on the best way to map riparian zones for broad spatial extents (e.g., large watersheds) when detailed maps from field data or high‐resolution imagery and terrain data are not available. Using well‐established accuracy metrics (e.g., kappa, precision, computational complexity), we evaluated eight methods commonly used to map riparian zones. Focusing on a semi‐arid, mountainous watershed, we found that the most accurate and robust method for mapping riparian zones combines data on upstream drainage area and valley topography. That method performed best regardless of stream order, and was most effective when implemented with fine resolution topographic and stream line data. Other commonly used methods to model riparian zones, such as those based on fixed‐width buffers, yielded inaccurate results. We recommend that until very‐high resolution (<1 m) elevation data are available at broad extents, models of riparian zones for semi‐arid mountainous regions should incorporate drainage area, valley topography, and quantify uncertainty.  相似文献   

13.
Stream fish bioassessment methods assume that fish assemblages observed in sample sites reflect responses to local stressors, but fish assemblages are influenced by local factors as well as regional dispersal to and from connected streams. We hypothesized that fish movement to and from refugia and source populations in connected rivers (i.e., riverine dispersal) would weaken or decouple relations between fish community metrics and local environmental conditions. We compared fish-environment relations between streams that flow into large rivers (mainstem tributaries) and streams that lack riverine confluences (headwater tributaries) at multiple spatial grains using data from the USEPA's Environmental Monitoring and Assessment Program in the mid-Atlantic highlands, USA (n = 157 sites). Headwater and mainstem tributaries were not different in local environmental conditions, but showed important differences in fish metric responses to environmental quality gradients. Stream sites flowing into mainstem channels within 10 fluvial km showed consistently weaker relations to local environmental conditions than stream sites that lacked such mainstem connections. Moreover, these patterns diminished at longer distances from riverine confluences, consistent with the hypothesis of riverine dispersal. Our results suggest that (1) the precision of fish bioassessment metrics may be improved by calibrating scoring criteria based on the spatial position of sites within stream networks and (2) the spatial grain of fish bioassessment studies may be manipulated to suit objectives by including or excluding fishes exhibiting riverine dispersal.  相似文献   

14.
This research investigates urban sprawl in the Greater Toronto Area (GTA) between 1985 and 2005 and the nature of the resulting landscape fragmentation, particularly with regard to the Oak Ridges Moraine (ORM), an ecologically important area for the region. Six scenes of Landsat TM imagery were acquired in summer of 1985, 1995, and 2005. These images and their texture measures were classified into eight land cover classes with very satisfactory final overall accuracies (93–95?%). Analysis of the classifications indicated that urban areas grew by 20?% between 1985 and 1995 and by 15?% between 1995 and 2005. Landscape fragmentation due to spatio-temporal land cover changes was evaluated using urban compactness indicators and landscape metrics, and results from the latter were used to draw conclusions about probable environmental impact. The indicator results showed that urban proportions increased in nearly all areas outside of the metropolitan center, including on portions of the ORM. The landscape metrics reveal that low density urban areas increased significantly in the GTA between 1985 and 2005, mainly at the expense of agricultural land. The metric results indicate increased vulnerability and exposure to adverse effects for natural and semi-natural land cover through greater contrast and lowered connectivity. The degree of urban perimeter increased around most environmentally significant areas in the region. Changes like these negatively impact species and the regional water supply in the GTA. Further investigation into specific environmental impacts of urban expansion in the region and which areas on the ORM are most at risk is recommended.  相似文献   

15.
ABSTRACT: Panchromatic black and white, color, and color infrared photographs and thermal infrared imagery are compared for a capability to show flood boundaries. In open agricultural and urban areas, these boundaries are easily delineated on all types of am. Boundaries are more difficult to see in wooded areas. In March, hardwood trees are dormant, but black and white photographs and color photographs show only the tops of these trees. Color infrared photographs in January and March have a distinctive color or tone in inundated woods; the limit of this tone is the flood boundary. Daytime thermal infrared imagery in March shows that inundated woods are cooler than dry land but warmer than open water. After about April 1, both color infrared photography and thermal infrared imagery show only the top of the tree canopy and do not reflect underlying flood water. Inundated areas can be delineated easily on ERTS satellite imagery from December through March. On imagery from May 4–5, 1973, however, most inundation boundaries had to be drawn as dashed lines; the tree canopy obscures flood waters in wooded areas. Despite this problem, the results of mapping flood boundaries on May imagery are believed to be reasonable for the scale of the imagery.  相似文献   

16.
To determine useful metrics for assessing stream water quality in the Southeastern Coastal Plain, we examined differences among two buffered and three unbuffered streams in an agricultural landscape in southwestern Georgia. Potential indicators included amphibian diversity and abundance, aquatic macroinvertebrate populations, riparian vegetative structure, water quality, and stream physical parameters. Variability among sites and treatments (buffered vs. unbuffered) existed, with sites in the same treatment as most similar, and disturbances from a nearby eroding gully strongly affecting one unbuffered site. Of the invertebrate metrics examined, percentages of clingers, Ephemeroptera-Plecoptera-Trichoptera (EPT), Elmidae (Coleoptera), Crustacea (Decapoda and Amphipoda), and dipterans were found to be possible indicators of stream health for perennial streams within this region. Overall, buffered sites showed higher percentages of sensitive invertebrate groups and showed lower and more stable concentrations of nitrate N, suspended solids, and fecal coliforms (FCs). Percent canopy cover was similar among sites; however, riparian vegetative coverage and percent leaf litter were greatest at buffered sites. No differences in amphibian abundance, presence, and absence within the riparian area were apparent between sites; however, instream larval salamanders were more abundant at buffered streams. In this study, stream buffers appeared to decrease nutrient and sediment loads to adjacent streams, enhancing overall water quality. Selected benthic macroinvertebrate metrics and amphibian abundance also appeared sensitive to agricultural influences. Amphibians show potential as indicator candidates, however further information is needed on their responses and tolerances to disturbances from the microhabitat to landscape levels.  相似文献   

17.
Hydroecological classification systems are typically based on an assemblage of streamflow metrics and seek to divide streams into ecologically relevant classes. Assignment of streams to classes is suggested as an initial step in the process of establishing ecological flow standards. We used two distinct hydroecological river classification systems available within North Carolina to evaluate the ability of a hydrologic model to assign the same classes as those determined by observed streamflows and to assess the transferability of such systems to ungaged streams. Class assignments were examined by rate of overall matches, rate of class matches, spatial variability in matches, and time period used in class assignment. The findings of this study indicate assignments of stream class: (1) are inconsistent among different classification systems; (2) differ between observed and modeled data; and (3) are sensitive to the period of record within observed data. One clear source of inconsistency/sensitivity in class assignments lies with the use of threshold values for metrics that distinguish stream classes, such that even small changes in metric values can result in different class assignments. Because these two hydroecological classification systems are representative of other classification systems that rely on quantitative decision thresholds, it can be surmised that the use of such systems based on stream flow metrics is not a reliable approach for guiding ecological flow determinations.  相似文献   

18.
We used invertebrate bioassessment, habitat analysis, geographic information system analysis of land use, and water chemistry monitoring to evaluate tributaries of a degraded northeast Nebraska, USA, reservoir. Bimonthly invertebrate collections and monthly water chemistry samples were collected for two years on six stream reaches to identify sources contributing to reservoir degradation and test suitability of standard rapid bioassessment methods in this region. A composite biotic index composed of seven commonly used metrics was effective for distinguishing between differentially impacted sites and responded to a variety of disturbances. Individual metrics varied greatly in precision and ability to discriminate between relatively impacted and unimpacted stream reaches. A modified Hilsenhoff index showed the highest precision (reference site CV = 0.08) but was least effective at discriminating among sites. Percent dominance and the EPT (number of Ephemeroptera, Plecoptera, and Trichoptera taxa) metrics were most effective at discriminating between sites and exhibited intermediate precision. A trend of higher biotic integrity during summer was evident, indicating seasonal corrections should differ from other regions. Poor correlations were evident between water chemistry variables and bioassessment results. However, land-use factors, particularly within 18-m riparian zones, were correlated with bioassessment scores. For example, there was a strong negative correlation between percentage of rangeland in 18-m riparian zones and percentage of dominance in streams (r 2 = 0.90, P < 0.01). Results demonstrate that standard rapid bioassessment methods, with some modifications, are effective for use in this agricultural region of the Great Plains and that riparian land use may be the best predictor of stream biotic integrity.  相似文献   

19.
Recovery of the benthic macroinvertebrate community in a small east Tennessee stream impacted by fly ash discharges from a power plant was investigated over a period of 6.5 years. The rate of recovery was greatest in the first 2 years after an initial 75% reduction in coal use led to a similar reduction in ash discharges and associated contaminants; further recovery followed after all fly ash discharges ceased. Recovery of the stream progressed through two phases. In the first phase, which lasted for approximately the first 2 years, most density and richness metrics increased considerably. In the second phase of recovery, the increases in metric values were followed by declines before fluctuating in and out of the lower reference ranges for the metrics. Detrended correspondence analyses and indicator species analyses showed that changes in species composition and community structure were ongoing throughout the second phase. Thus, the first phase was characterized by species additions, while the second phase involved species replacements and shifts in community dominants. Further recovery of the macroinvertebrate community will probably depend on additional flushing of fly ash deposits from the streambed and flood plain, because their continued presence reduces habitat quality in the stream and serves as a potential source of contaminants. Further recovery also may be limited by the availability of vagile species in nearby watersheds.  相似文献   

20.
Riparian forests attenuate solar radiation, thereby mediating an important component of the thermal budget of streams. Here, we investigate the relationship between riparian degradation, stream temperature, and channel width in the Chehalis River Basin, Washington State. We used lidar data to measure canopy opening angle, the angle formed between the channel center and trees on both banks; we assumed historical tree heights and calculated the change in canopy angle relative to historical conditions. We then developed an empirical relationship between canopy angle and water temperature using existing data, and simulated temperatures between 2002 and 2080 by combining a tree growth model with climate change scenarios from the NorWeST regional prediction. The greatest change between historical and current conditions (~7°C) occurred in developed portions of the river network, with the highest values of change predicted at channel widths less than ~40 m. Tree growth lessened climate change increases in maximum temperature and the length of river exceeding biologically critical thresholds by ~50%–60%. Moreover, the maximum temperature of channels with bankfull widths less than ~50 m remained similar to current conditions, despite climate change increases. Our findings are consistent with a possible role for the riparian landscape in explaining the low sensitivity of stream temperatures to air temperatures observed in some small mountain streams.  相似文献   

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

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