首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 57 毫秒
1.
ABSTRACT: About 50 to 80 percent of precipitation in the southeastern United States returns to the atmosphere by evapotranspiration. As evapotranspiration is a major component in the forest water balances, accurately quantifying it is critical to predicting the effects of forest management and global change on water, sediment, and nutrient yield from forested watersheds. However, direct measurement of forest evapotranspiration on a large basin or a regional scale is not possible. The objectives of this study were to develop an empirical model to estimate long‐term annual actual evapotranspiration (ART) for forested watersheds and to quantify spatial AET patterns across the southeast. A geographic information system (GIS) database including land cover, daily streamflow, and climate was developed using long term experimental and monitoring data from 39 forested watersheds across the region. Using the stepwise selection method implemented in a statistical modeling package, a long term annual AET model was constructed. The final multivariate linear model includes four independent variables—annual precipitation, watershed latitude, watershed elevation, and percentage of forest coverage. The model has an adjusted R2 of 0.794 and is sufficient to predict long term annual ART for forested watersheds across the southeastern United States. The model developed by this study may be used to examine the spatial variability of water availability, estimate annual water loss from mesoscale watersheds, and project potential water yield change due to forest cover change.  相似文献   

2.
Wetland ecosystems are profoundly affected by altered nutrient and sediment loads received from anthropogenic activity in their surrounding watersheds. Our objective was to compare a gradient of agricultural and urban land cover history during the period from 1949 to 1997, with plant and soil nutrient concentrations in, and sediment deposition to, riparian wetlands in a rapidly urbanizing landscape. We observed that recent agricultural land cover was associated with increases in Nitrogen (N) and Phosphorus (P) concentrations in a native wetland plant species. Conversely, recent urban land cover appeared to alter receiving wetland environmental conditions by increasing the relative availability of P versus N, as reflected in an invasive, but not a native, plant species. In addition, increases in surface soil Fe content suggests recent inputs of terrestrial sediments associated specifically with increasing urban land cover. The observed correlation between urban land cover and riparian wetland plant tissue and surface soil nutrient concentrations and sediment deposition, suggest that urbanization specifically enhances the suitability of riparian wetland habitats for the invasive species Japanese stiltgrass [Microstegium vimenium (Trinius) A. Camus].  相似文献   

3.
ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.  相似文献   

4.
ABSTRACT: Simulated water quality resulting from three alternative future land‐use scenarios for two agricultural watersheds in central Iowa was compared to water quality under current and historic land use/land cover to explore both the potential water quality impact of perpetuating current trends and potential benefits of major changes in agricultural practices in the U.S. Corn Belt. The Soil Water Assessment Tool (SWAT) was applied to evaluate the effect of management practices on surface water discharge and annual loads of sediment and nitrate in these watersheds. The agricultural practices comprising Scenario 1, which assumes perpetuation of current trends (conversion to conservation tillage, increase in farm size and land in production, use of currently‐employed Best Management Practices (BMPs)) result in simulated increased export of nitrate and decreased export of sediment relative to the present. However, simulations indicate that the substantial changes in agricultural practices envisioned in Scenarios 2 and 3 (conversion to conservation tillage, strip intercropping, rotational grazing, conservation set‐asides and greatly extended use of best management practices (BMPs) such as riparian buffers, engineered wetlands, grassed waterways, filter strips and field borders) could potentially reduce current loadings of sediment by 37 to 67 percent and nutrients by 54 to 75 percent. Results from the study indicate that major improvements in water quality in these agricultural watersheds could be achieved if such environmentally‐targeted agricultural practices were employed. Traditional approaches to water quality improvement through application of traditional BMPs will result in little or no change in nutrient export and minor decreases in sediment export from Corn Belt watersheds.  相似文献   

5.
This study quantified nonpoint source nitrogen (NPS‐N) sources and sinks across the 14,582 km2 Neuse River Basin (NRB) located in North Carolina, to provide tabular data summaries and graphic overlay products to support the development of management approaches to best achieve established N reduction goals. First, a remote sensor derived, land cover classification was performed to support modeling needs. Modeling efforts included the development of a mass balance model to quantify potential N sources and sinks, followed by a precipitation event driven hydrologic model to effectively transport excess N across the landscape to individual stream reaches to support subsequent labeling of transported N values corresponding to source origin. Results indicated that agricultural land contributed 55 percent of the total annual NPS‐N loadings, followed by forested land at 23 percent (background), and urban areas at 21 percent. Average annual N source contributions were quantified for agricultural (1.4 kg/ha), urban (1.2 kg/ha), and forested cover types (0.5 kg/ha). Nonpoint source‐N contributions were greatest during the winter (40 percent), followed by spring (32 percent), summer (28 percent), and fall (0.3 percent). Seasonal total N loadings shifted from urban dominated and forest dominated sources during the winter, to agricultural sources in the spring and summer. A quantitative assessment of the significant NRB land use activities indicated that high (greater than 70 percent impervious) and medium (greater than 35 percent impervious) density urban development were the greatest contributors of NPS‐N on a unit area basis (1.9 and 1.6 kg/ha/yr, respectively), followed by row crops and pasture/hay cover types (1.4 kg/ha/yr).  相似文献   

6.
Detecting Temporal Change in Watershed Nutrient Yields   总被引:2,自引:1,他引:1  
Meta-analyses reveal that nutrient yields tend to be higher for watersheds dominated by anthropogenic uses (e.g., urban, agriculture) and lower for watersheds dominated by natural vegetation. One implication of this pattern is that loss of natural vegetation will produce increases in watershed nutrient yields. Yet, the same meta-analyses also reveal that, absent land-cover change, watershed nutrient yields vary from one year to the next due to many exogenous factors. The interacting effects of land cover and exogenous factors suggest nutrient yields should be treated as distributions, and the effect of land-cover change should be examined by looking for significant changes in the distributions. We compiled nutrient yield distributions from published data. The published data included watersheds with homogeneous land cover that typically reported two or more years of annual nutrient yields for the same watershed. These data were used to construct statistical models, and the models were used to estimate changes in the nutrient yield distributions as a result of land-cover change. Land-cover changes were derived from the National Land Cover Database (NLCD). Total nitrogen (TN) yield distributions increased significantly for 35 of 1550 watersheds and decreased significantly for 51. Total phosphorus (TP) yield distributions increased significantly for 142 watersheds and decreased significantly for 17. The amount of land-cover change required to produce significant shifts in nutrient yield distributions was not constant. Small land-cover changes led to significant shifts in nutrient yield distributions when watersheds were dominated by natural vegetation, whereas much larger land-cover changes were needed to produce significant shifts when watersheds were dominated by urban or agriculture. We discuss our results in the context of the Clean Water Act.  相似文献   

7.
ABSTRACT: Many water bodies within the United States are contaminated by non‐point source (NPS) pollution, which is defined as those materials posing a threat to water quality arising from a number of individual sources and diffused through hydrologic processes. One such NPS pollutant that is of critical concern are pathogens derived from animal wastes, including humans. The potential presence of pathogens is identified by testing the water for fecal conform, a bacteria also associated with animal wastes. Water contaminated by animal wastes are most often associated with urban and agricultural areas, thus it is postulated that by utilizing land cover indicators, those water bodies that may be at risk of fecal coliform contamination may be identified. This study utilizes land cover information derived from the Multi‐Resolution Land Characterization (MRLC) project to analyze fecal coliform contamination in South Carolina. Also utilized are 14 digit hydro‐logic unit code (HUC) watersheds of the state, a digital elevation model, and test point data stating whether fecal coliform levels exceeded State Water Quality Standards. Proportions of the various land covers are identified within the individual watersheds and then analyzed using a logistic regression. The results reveal that watersheds with large proportions of urban land cover and agriculture on steep slopes had a very high probability of being impaired. (KEY TERMS: Geographic Information Systems; land use planning; nonpoint source pollution; statistical analysis; water quality; watershed management.)  相似文献   

8.
Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.  相似文献   

9.
Red alder (Alnus rubra), a nitrogen(N)‐fixing deciduous broadleaf tree, can strongly influence N concentrations in western Oregon and Washington. We compiled a database of stream N and GIS‐derived landscape characteristics in order to examine geographic variation in N across the Oregon Coast Range. Basal area of alder, expressed as a percent of watershed area, accounted for 37% and 38% of the variation in summer nitrate and total N (TN) concentrations, respectively. Relationships between alder and nitrate were strongest in winter when streamflow and landscape connections are highest. Distance to the coast and latitude, potential surrogates for sea salt inputs, and watershed area were also related to nitrate concentrations in an all‐subsets regression analysis, which accounted for 46% of the variation in summer nitrate concentrations. The model with the lowest Akaike's Information Criterion did not include developed or agricultural land cover, probably because few watersheds in our database had substantial levels of these land cover classes. Our results provide evidence, at a regional scale, that background sources and processes cause many Coast Range streams to exceed proposed nutrient criteria, and that the prevalence of a single tree species (N‐fixing red alder) exerts a dominant control over stream N concentrations across this region.  相似文献   

10.
Abstract: The spatial scale and location of land whose development has the strongest influence on aquatic ecosystems must be known to support land use decisions that protect water resources in urbanizing watersheds. We explored impacts of urbanization on streams in the West River watershed, New Haven, Connecticut, to identify the spatial scale of watershed imperviousness that was most strongly related to water chemistry, macroinvertebrates, and physical habitat. A multiparameter water quality index was used to characterize regional urban nonpoint source pollution levels. We identified a critical level of 5% impervious cover, above which stream health declined. Conditions declined with increasing imperviousness and leveled off in a constant state of impairment at 10%. Instream variables were most correlated (0.77 ≤ |r| ≤ 0.92, p < 0.0125) to total impervious area (TIA) in the 100‐m buffer of local contributing areas (~5‐km2 drainage area immediately upstream of each study site). Water and habitat quality had a relatively consistent strong relationship with TIA across each of the spatial scales of investigation, whereas macroinvertebrate metrics produced noticeably weaker relationships at the larger scales. Our findings illustrate the need for multiscale watershed management of aquatic ecosystems in small streams flowing through the spatial hierarchies that comprise watersheds with forest‐urban land use gradients.  相似文献   

11.
Abstract: Impervious cover is a commonly used metric to help explain or predict anthropogenic impacts on aquatic resources; often it is used as a surrogate for intensity of human impacts when evaluating effects on aquatic resources. The most common way to estimate imperviousness is based on relationships with land use. Few studies have evaluated how the relationship between impervious surface and land use varies among geographies with different levels of development and between types of imagery used to assign land use type. In this study, we assess variability in estimates of imperviousness based on two locally available land use datasets: one based on aerial imagery (2‐m resolution) and another based on satellite imagery (30‐m resolution). The ranges and variability in imperviousness within land use categories were assessed at several spatial scales, including within counties, between counties, and between watersheds. Results indicate that there was considerable variability for all developed land use types. Estimated impervious cover often varied over a range of 20‐40% points within a land use category. Furthermore, there were clear spatial patterns both between and within counties, with impervious cover for a given land use type being higher near the urban centers and lower at the margins of development. Estimates of imperviousness for 12 study watersheds indicated that variability increased with increasing watershed development, making it difficult to confidently set management or regulatory targets based on impervious cover. This study suggests that locally derived, high resolution satellite or aerial imagery should be used to estimate imperviousness when a high level of accuracy and precision is required for regulatory or management decisions. Furthermore, the error associated with impervious land use relationships should be accounted for when using impervious cover in runoff or water quality models, or when making management decisions regarding stream health.  相似文献   

12.
ABSTRACT: We compared watershed land‐use and fish community data between the 1970s and 1990s in 47 small streams in southeastern Wisconsin. Our goal was to quantify effects of increasing urbanization on stream fishes in what had been a predominantly agricultural region. In the 43 test watersheds, mean surface coverage by agricultural lands decreased from 54 percent to 43 percent and urban lands increased from 24 percent to 31 percent between 1970 and 1990. Agriculture dominated the four reference watersheds, but neither agriculture (65–59 percent) nor urban (4.4–4.8 percent) land‐uses changed significantly in those watersheds during the study period. From the 1970s to the 1990s the mean number of fish species for the test stream sites decreased 15 percent, fish density decreased 41 percent, and the index of biotic integrity (IBI) score dropped 32 percent. Fish community attributes at the four reference sites did not change significantly during the same period, although density was substantially lower in the 1990s. For both the 1970s and 1990s test sites, numbers of fish species and IBI scores were positively correlated with watershed percent agricultural land coverage and negatively correlated with watershed urban land uses, as indexed by percent effective connected imperviousness. Numbers of fish species per site and IBI scores were highly variable below 10 percent imperviousness, but consistently low above 10 percent. Sites that had less than 10 percent imperviousness and fewer than 10 fish species in the 1970s suffered the greatest relative increase in imperviousness and decline in species number over the study period. Our findings are consistent with previous studies that have found strong negative effects of urban land uses on stream ecosystems and a threshold of environmental damage at about 10 percent imperviousness. We conclude that although agricultural land uses often degrade stream fish communities, agricultural land impacts are generally less severe than those from urbanization on a per‐unit‐area basis.  相似文献   

13.
Sanford, Ward E. and David L. Selnick, 2012. Estimation of Evapotranspiration Across the Conterminous United States Using a Regression with Climate and Land‐Cover Data. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12010 Abstract: Evapotranspiration (ET) is an important quantity for water resource managers to know because it often represents the largest sink for precipitation (P) arriving at the land surface. In order to estimate actual ET across the conterminous United States (U.S.) in this study, a water‐balance method was combined with a climate and land‐cover regression equation. Precipitation and streamflow records were compiled for 838 watersheds for 1971‐2000 across the U.S. to obtain long‐term estimates of actual ET. A regression equation was developed that related the ratio ET/P to climate and land‐cover variables within those watersheds. Precipitation and temperatures were used from the PRISM climate dataset, and land‐cover data were used from the USGS National Land Cover Dataset. Results indicate that ET can be predicted relatively well at a watershed or county scale with readily available climate variables alone, and that land‐cover data can also improve those predictions. Using the climate and land‐cover data at an 800‐m scale and then averaging to the county scale, maps were produced showing estimates of ET and ET/P for the entire conterminous U.S. Using the regression equation, such maps could also be made for more detailed state coverages, or for other areas of the world where climate and land‐cover data are plentiful.  相似文献   

14.
ABSTRACT: Multivariate analyses and correlations revealed strong relations between watershed and riparian‐corridor land cover, and reach‐scale habitat versus fish and macroinvertebrate assemblages in 38 warmwater streams in eastern Wisconsin. Watersheds were dominated by agricultural use, and ranged in size from 9 to 71 km2 Watershed land cover was summarized from satellite‐derived data for the area outside a 30‐m buffer. Riparian land cover was interpreted from digital orthophotos within 10‐, 10‐to 20‐, and 20‐to 30‐m buffers. Reach‐scale habitat, fish, and macroinvertebrates were collected in 1998 and biotic indices calculated. Correlations between land cover, habitat, and stream‐quality indicators revealed significant relations at the watershed, riparian‐corridor, and reach scales. At the watershed scale, fish diversity, intolerant fish and EPT species increased, and Hilsenhoff biotic index (HBI) decreased as percent forest increased. At the riparian‐corridor scale, EPT species decreased and HBI increased as riparian vegetation became more fragmented. For the reach, EPT species decreased with embeddedness. Multivariate analyses further indicated that riparian (percent agriculture, grassland, urban and forest, and fragmentation of vegetation), watershed (percent forest) and reach‐scale characteristics (embeddedness) were the most important variables influencing fish (IBI, density, diversity, number, and percent tolerant and insectivorous species) and macroinvertebrate (HBI and EPT) communities.  相似文献   

15.
ABSTRACT: The introduction of nutrients from chemical fertilizer, animal manure, wastewater, and atmospheric deposition to the eastern Iowa environment creates a large potential for nutrient transport in watersheds. Agriculture constitutes 93 percent of all land use in eastern Iowa. As part of the U.S. Geological Survey National Water Quality Assessment Program, water samples were collected (typically monthly) from six small and six large watersheds in eastern Iowa between March 1996 and September 1997. A Geographic Information System (GIS) was used to determine land use and quantify inputs of nitrogen and phosphorus within the study area. Streamliow from the watersheds is to the Mississippi River. Chemical fertilizer and animal manure account for 92 percent of the estimated total nitrogen and 99.9 percent of the estimated total phosphorus input in the study area. Total nitrogen and total phosphorus loads for 1996 were estimated for nine of the 12 rivers and creeks using a minimum variance unbiased estimator model. A seasonal pattern of concentrations and loads was observed. The greatest concentrations and loads occur in the late spring to early summer in conjunction with row‐crop fertilizer applications and spring nmoff and again in the late fall to early winter as vegetation goes into dormancy and additional fertilizer is applied to row‐crop fields. The three largest rivers in eastern Iowa transported an estimated total of 79,000 metric tons of total nitrogen and 6,800 metric tons of total phosphorus to the Mississippi River in 1996. The estimated mass of total nitrogen and total phosphorus transported to the Mississippi River represents about 19 percent of all estimated nitrogen and 9 percent of all estimated phosphorus input to the study area.  相似文献   

16.
Abstract:  It is critical that evapotranspiration (ET) be quantified accurately so that scientists can evaluate the effects of land management and global change on water availability, streamflow, nutrient and sediment loading, and ecosystem productivity in watersheds. The objective of this study was to derive a new semi‐empirical ET modeled using a dimension analysis method that could be used to estimate forest ET effectively at multiple temporal scales. The model developed describes ET as a function of water availability for evaporation and transpiration, potential ET demand, air humidity, and land surface characteristics. The model was tested with long‐term hydrometeorological data from five research sites with distinct forest hydrology in the United States and China. Averaged simulation error for daily ET was within 0.5 mm/day. The annual ET at each of the five study sites were within 7% of measured values. Results suggest that the model can accurately capture the temporal dynamics of ET in forest ecosystems at daily, monthly, and annual scales. The model is climate‐driven and is sensitive to topography and vegetation characteristics and thus has potential to be used to examine the compounding hydrologic responses to land cover and climate changes at multiple temporal scales.  相似文献   

17.
Changing climate and land cover are expected to impact flood hydrology in the Delaware River Basin over the 21st Century. HEC‐HMS models (U.S. Army Corps of Engineers Hydrologic Engineering Center‐Hydrologic Modeling System) were developed for five case study watersheds selected to represent a range of scale, soil types, climate, and land cover. Model results indicate that climate change alone could affect peak flood discharges by ?6% to +58% a wide range that reflects regional variation in projected rainfall and snowmelt and local watershed conditions. Land cover changes could increase peak flood discharges up to 10% in four of the five watersheds. In those watersheds, the combination of climate and land cover change increase modeled peak flood discharges by up to 66% and runoff volumes by up to 44%. Precipitation projections are a key source of uncertainty, but there is a high likelihood of greater precipitation falling on a more urbanized landscape that produces larger floods. The influence of climate and land cover changes on flood hydrology for the modeled watersheds varies according to future time period, climate scenario, watershed land cover and soil conditions, and flood frequency. The impacts of climate change alone are typically greater than land cover change but there is substantial geographic variation, with urbanization the greater influence on some small, developing watersheds.  相似文献   

18.
ABSTRACT: Existing land use data were used to estimate nonpoint source phosphorus loads to Lake Champlain (Vermont/New York/Quebec) in a loading function model that combined P concentration coefficients with regional hydrologic data. The estimates were verified against monitored loading data, then used to assess the relative magnitudes of contributions from major land uses and regions of the Lake Champlain Basin. The Basin is comprised of 62 percent forest, 28 percent agricultural land, 3 percent urban land, and 7 percent water. The best-fit model estimated an annual total P load of 457 mt/year, which did not differ significantly from the 458 metric tons/year measured for an average hydrologic year, and accurately predicted loads from major tributaries. Agriculture contributes 66 percent of the annual nonpoint source P load to Lake Champlain; urban and forest land contribute 18 percent and 16 percent, respectively. Because agricultural land contributes most nonpoint source P to Lake Champlain, load reduction effort must deal with agricultural sources. However, because the urban 3 percent of the basin contributes 18 percent of the estimated load, high load reduction efficiencies might be achieved by addressing urban sources. This assessment clearly demonstrated the relationship between land use and P loads in the Lake Champlain Basin, a prerequisite for policy-makers to endorse a P management strategy requiring changes in land use and management.  相似文献   

19.
ABSTRACT: A loading function methodology is presented for predicting runoff, sediment, and nutrient losses from complex watersheds. Separate models are defined for cropland, forest, urban and barnyard sources, and procedures for estimating baseflow nutrients are provided. The loading functions are designed for use as a preliminary screening tool to isolate the major contributors in a watershed. Input data sources are readily available and the functions do not require costly calibrations. Data requirements include watershed land use and soil information, daily precipitation and temperature records and rainfall erosivities. Comparison of predicted and measured water, sediment, and nutrient runoff fluxes for the West Branch Deleware River in New York, indicated that runoff was underpredicted by about 14 percent while dissolved nutrients were within 30 percent of observed values. Sediment and solid-phase nutrients were overpredicted by about 50 percent. An annual nutrient budget for the West Branch Delaware River showed that cornland was the major source of sediment, solid phase nutrients, and total phosphorus. Waste water treatment plants and ground water discharge contributed the most dissolved phosphorus and dissolved nitrogen, respectively.  相似文献   

20.
Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.  相似文献   

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

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