首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
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.  相似文献   

2.
This study presents new data‐driven, annual estimates of the division of precipitation into the recharge, quick‐flow runoff, and evapotranspiration (ET) water budget components for 2000‐2013 for the contiguous United States (CONUS). The algorithms used to produce these maps ensure water budget consistency over this broad spatial scale, with contributions from precipitation influx attributed to each component at 800 m resolution. The quick‐flow runoff estimates for the contribution to the rapidly varying portion of the hydrograph are produced using data from 1,434 gaged watersheds, and depend on precipitation, soil saturated hydraulic conductivity, and surficial geology type. Evapotranspiration estimates are produced from a regression using water balance data from 679 gaged watersheds and depend on land cover, temperature, and precipitation. The quick‐flow and ET estimates are combined to calculate recharge as the remainder of precipitation. The ET and recharge estimates are checked against independent field data, and the results show good agreement. Comparisons of recharge estimates with groundwater extraction data show that in 15% of the country, groundwater is being extracted at rates higher than the local recharge. These maps of the internally consistent water budget components of recharge, quick‐flow runoff, and ET, being derived from and tested against data, are expected to provide reliable first‐order estimates of these quantities across the CONUS, even where field measurements are sparse.  相似文献   

3.
This study analyzed changes in hydrology between two recent decades (1980s and 2010s) with the Soil and Water Assessment Tool (SWAT) in three representative watersheds in South Dakota: Bad River, Skunk Creek, and Upper Big Sioux River watersheds. Two SWAT models were created over two discrete time periods (1981‐1990 and 2005‐2014) for each watershed. National Land Cover Datasets 1992 and 2011 were, respectively, ingested into 1981‐1990 and 2005‐2014 models, along with corresponding weather data, to enable comparison of annual and seasonal runoff, soil water content, evapotranspiration (ET), water yield, and percolation between these two decades. Simulation results based on the calibrated models showed that surface runoff, soil water content, water yield, and percolation increased in all three watersheds. Elevated ET was also apparent, except in Skunk Creek watershed. Differences in annual water balance components appeared to follow changes in land use more closely than variation in precipitation amounts, although seasonal variation in precipitation was reflected in seasonal surface runoff. Subbasin‐scale spatial analyses revealed noticeable increases in water balance components mostly in downstream parts of Bad River and Skunk Creek watersheds, and the western part of Upper Big Sioux River watershed. Results presented in this study provide some insight into recent changes in hydrological processes in South Dakota watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.  相似文献   

4.
Abstract: Natural forests in southern China have been severely logged due to high human demand for timber, food, and fuels during the past century, but are recovering in the past decade. The objective of this study was to investigate how vegetation cover changes in composition and structure affected the water budgets of a 9.6‐km2 Dakeng watershed located in a humid subtropical mountainous region in southern China. We analyzed 27 years (i.e., 1967‐1993) of streamflow and climate data and associated vegetation cover change in the watershed. Land use/land cover census and Normalized Difference of Vegetation Index (NDVI) data derived from remote sensing were used to construct historic land cover change patterns. We found that over the period of record, annual streamflow (Q) and runoff/precipitation ratio did not change significantly, nor did the climatic variables, including air temperature, Hamon’s potential evapotranspiration (ET), pan evaporation, sunshine hours, and radiation. However, annual ET estimated as the differences between P and Q showed a statistically significant increasing trend. Overall, the NDVI of the watershed had a significant increasing trend in the peak spring growing season. This study concluded that watershed ecosystem ET increased as the vegetation cover shifted from low stock forests to shrub and grasslands that had higher ET rates. A conceptual model was developed for the study watershed to describe the vegetation cover‐streamflow relationships during a 50‐year time frame. This paper highlighted the importance of eco‐physiologically based studies in understanding transitory, nonstationary effects of deforestation or forestation on watershed water balances.  相似文献   

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

6.
Abstract: Limited information exists on pesticide use for nonagricultural purposes, making it difficult to estimate pesticide loadings from nonagricultural sources to surface water and to conduct environmental risk assessments. A method was developed to estimate the amount of pesticide use on recreational turf grasses, specifically golf course turf grasses, for watersheds located throughout the conterminous United States (U.S.). The approach estimates pesticide use: (1) based on the area of recreational turf grasses (used as a surrogate for turf associated with golf courses) within the watershed, which was derived from maps of land cover, and (2) from data on the location and average treatable area of golf courses. The area of golf course turf grasses determined from these two methods was used to calculate the percentage of each watershed planted in golf course turf grass (percent crop area, or PCA). Turf‐grass PCAs derived from the two methods were used with recommended application rates provided on pesticide labels to estimate total pesticide use on recreational turf within 1,606 watersheds associated with surface‐water sources of drinking water. These pesticide use estimates made from label rates and PCAs were compared to use estimates from industry sales data on the amount of each pesticide sold for use within the watershed. The PCAs derived from the land‐cover data had an average value of 0.4% of a watershed with minimum of 0.01% and a maximum of 9.8%, whereas the PCA values that are based on the number of golf courses in a watershed had an average of 0.3% of a watershed with a minimum of <0.01% and a maximum of 14.2%. Both the land‐cover method and the number of golf courses method produced similar PCA distributions, suggesting that either technique may be used to provide a PCA estimate for recreational turf. The average and maximum PCAs generally correlated to watershed size, with the highest PCAs estimated for small watersheds. Using watershed specific PCAs, combined with label rates, resulted in greater than two orders of magnitude over‐estimation of the pesticide use compared to estimates from sales data.  相似文献   

7.
Abstract: Information on evapotranspiration (ET) can help us understand water balance, particularly in forested watersheds. Previous studies in China show that ET was relatively low (30‐40% of total precipitation) in the Minjiang Valley located in the upper reach of the Yangtze River Basin. However, this conclusion was derived from research on small‐scale watersheds (<100 km2). The objective of this paper was to present ET information on meso‐scale watersheds in the Minjiang Valley. Four meso‐scale watersheds (1,700‐5,600 km2) located in the Minjiang Valley were used to estimate ET using the water balance approach. We first generated forest vegetation variables (coniferous forest percentage, forest cover percentage, and derived forest vegetation index) using remote sensing data. Landsat 5 TM satellite images, acquired on June 26, 1994, were selected for the vegetation classification. Actual annual ET was calculated based on 11‐year estimated precipitation and measured streamflow data (1992‐2002). We also calculated potential ET (PET) using an improved Thornthwaite model for all four watersheds for the period of 1992‐1998. PET can provide additional information about potential capacity of water flux to atmosphere in the region. Seasonal (dry and rainy) PET and ET for all studied watersheds were also estimated for comparison purposes as the water balance approach, at shorter than annual scales, would likely provide inaccurate estimates of ET. The dominant vegetations in the Minjiang Valley were grasslands, conifer forests, and shrub‐lands. Our results confirmed that both ET and PET for three studied meso‐scale watersheds in the Minjiang Valley is relatively low (39.5‐43.8 and 28.2‐47.7% for ET and PET, respectively), with an exception of ET in the Yuzixi watershed being 71.1%. This result is generally consistent with previous research at small watershed scales. Furthermore, the low ET across various scales in the Minjiang Valley may be related to the unique deeply cut valley environment.  相似文献   

8.
This study quantitatively explores whether land cover changes have a substantive impact on simulated streamflow within the tropical island setting of Puerto Rico. The Precipitation Runoff Modeling System (PRMS) was used to compare streamflow simulations based on five static parameterizations of land cover with those based on dynamically varying parameters derived from four land cover scenes for the period 1953‐2012. The PRMS simulations based on static land cover illustrated consistent differences in simulated streamflow across the island. It was determined that the scale of the analysis makes a difference: large regions with localized areas that have undergone dramatic land cover change may show negligible difference in total streamflow, but streamflow simulations using dynamic land cover parameters for a highly altered subwatershed clearly demonstrate the effects of changing land cover on simulated streamflow. Incorporating dynamic parameterization in these highly altered watersheds can reduce the predictive uncertainty in simulations of streamflow using PRMS. Hydrologic models that do not consider the projected changes in land cover may be inadequate for water resource management planning for future conditions.  相似文献   

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

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

11.
Restored annual streamflow (Qr) and measured daily streamflow of the Chaohe watershed located in northern China and associated long‐term climate and land use/cover data were used to explore the effects of land use/cover change and climate variability on the streamflow during 1961‐2009. There were no significant changes in annual precipitation (P) and potential evapotranspiration, whereas Qr decreased significantly by 0.81 mm/yr (< 0.001) over the study period with a change point in 1999. We used 1961‐1998 as the baseline period (BP) and 1999‐2009 the change period (CP). The mean Qr during the CP decreased by 39.4 mm compared with that in the BP. From 1979 to 2009, the grassland area declined by 69.6%, and the forest and shrublands increased by 105.4 and 73.1%, respectively. The land use/cover change and climate variability contributed for 58.4 and 41.6% reduction in mean annual Qr, respectively. Compared with the BP, median and high flows in the CP decreased by 38.8 and up to 75.5%, respectively. The study concludes that large‐scale ecological restoration and watershed management in northern China has greatly decreased water yield and reduced high flows due to the improved land cover by afforestation leading to higher water loss through evapotranspiration. At a large watershed scale, land use/cover change could play as much of an important role as climate variability on water resources.  相似文献   

12.
Schilling, Keith E., Thomas M. Isenhart, Jason A. Palmer, Calvin F. Wolter, and Jean Spooner, 2011. Impacts of Land‐Cover Change on Suspended Sediment Transport in Two Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):672‐686. DOI: 10.1111/j.1752‐1688.2011.00533.x Abstract: Suspended sediment is a major water quality problem, yet few monitoring studies have been of sufficient scale and duration to assess the effectiveness of land‐use change or conservation practice implementation at a watershed scale. Daily discharge and suspended sediment export from two 5,000‐ha watersheds in central Iowa were monitored over a 10‐year period (water years 1996‐2005). In Walnut Creek watershed, a large portion of land was converted from row crop to native prairie, whereas in Squaw Creek land use remained predominantly row crop agriculture. Suspended sediment loads were similar in both watersheds, exhibiting flashy behavior typical of incised channels. Modeling suggested that expected total soil erosion in Walnut Creek should have been reduced 46% relative to Squaw Creek due to changes in land use, yet measured suspended sediment loads showed no significant differences. Stream mapping indicated that Walnut Creek had three times more eroding streambank lengths than did Squaw Creek suggesting that streambank erosion dominated sediment sources in Walnut Creek and sheet and rill sources dominated sediment sources in Squaw Creek. Our results demonstrate that an accounting of all sources of sediment erosion and delivery is needed to characterize sediment reductions in watershed projects combined with long‐term, intensive monitoring and modeling to account for possible lag times in the manifestation of the benefits of conservation practices on water quality.  相似文献   

13.
ABSTRACT: Nonirrigated crop yields and forage production are limited by low and variable precipitation in the southern Great Plains. Precipitation variation involves production risks, which can be reduced by considering probability of precipitation, precipitation retention, and soil erosion under various production systems. The objective of this study was to probabilistically quantify the impact of precipitation variations, land use, cropping, and tillage systems on precipitation retention and soil erosion. Five 1.6 ha watersheds that had 3 to 4 percent slopes, and similar silt loam soils were selected. One was kept in native grass, and the others were planted into winter wheat (Triticum aestivum L.) under different cropping and tillage systems. Daily runoff and soil erosion were measured at the outlet of each watershed. Precipitation distributions exhibited great seasonal and interannual variations, and precipitation retention distributions resembled those of precipitation. Cropping and tillage systems affected precipitation retention but much less than did precipitation variations. Available soil water storage, which was largely controlled by ET, played an important role in retaining precipitation. This indicates that cropping systems should be adjusted to precipitation patterns, if predictable, for better soil water use. Land use and cropping and tillage systems had a much greater impact on soil erosion than on precipitation retention. Soil erosion risks, which were proportional to the levels of tillage disturbance, were mainly caused by a few large storms in summer, when surface cover was low. This study explored a novel approach for evaluating production risks associated with insufficient precipitation retention and excessive soil erosion for certain crops or cropping systems under assumed precipitation conditions.  相似文献   

14.
Assessment of water resources at a national scale is critical for understanding their vulnerability to future change in policy and climate. Representation of the spatiotemporal variability in snowmelt processes in continental‐scale hydrologic models is critical for assessment of water resource response to continued climate change. Continental‐extent hydrologic models such as the U.S. Geological Survey National Hydrologic Model (NHM) represent snowmelt processes through the application of snow depletion curves (SDCs). SDCs relate normalized snow water equivalent (SWE) to normalized snow covered area (SCA) over a snowmelt season for a given modeling unit. SDCs were derived using output from the operational Snow Data Assimilation System (SNODAS) snow model as daily 1‐km gridded SWE over the conterminous United States. Daily SNODAS output were aggregated to a predefined watershed‐scale geospatial fabric and used to also calculate SCA from October 1, 2004 to September 30, 2013. The spatiotemporal variability in SNODAS output at the watershed scale was evaluated through the spatial distribution of the median and standard deviation for the time period. Representative SDCs for each watershed‐scale modeling unit over the conterminous United States (n = 54,104) were selected using a consistent methodology and used to create categories of snowmelt based on SDC shape. The relation of SDC categories to the topographic and climatic variables allow for national‐scale categorization of snowmelt processes.  相似文献   

15.
Nitrogen (N) losses from agricultural lands in the Midwest United States are contributing to the expansion of the hypoxic zone in the Gulf of Mexico. This study evaluated the importance of inter‐annual variability in precipitation, land cover, and N fertilizer use on NO3 + NO2‐N loads in seven United States Midwestern Rivers using the backward stepwise regression analysis. At the annual scale, fluctuations in the current and previous years’ precipitations explained much of the variation in streamflow, baseflow, and N‐load. Previous years precipitation effects were associated with fillable soil porosity. In some years, higher residual soil N from previous dry years also contributed to an increase in N‐load. Area under soybean production (SOY), a surrogate for replacement of prairies and small grains was generally not a significant explanatory variable. Fertilizer use from 1987 to 2012 was also not a significant explanatory variable in the annual analysis. Precipitation in both the current and previous months and previous year were important in explaining variation in monthly streamflow, baseflow, and N‐load. SOY was significant in one or two months from June to August, but had a higher p‐value than precipitation. We conclude recent increases in river N‐loads are primarily due to wet climate and minimally due to the changes in land cover or N fertilizer use. Under current cropping systems and agronomic N application rates, tile water remediation will be necessary to reduce river N‐loads.  相似文献   

16.
ABSTRACT: Annual evapotranspiration from a watershed (ETws) is a function of annual precipitation (P) and fraction of the watershed covered by shrubs and trees (C). Other characteristics are not significant in explaining variance of ETws. A rational equation, ETws = (1-C) ETg + CETst, wherein ETg is the ET of herbaceous cover and ETst is the ET of shrubs and trees, is proposed. The equation has been calibrated for ET and P in inches for the watershed of Lake Cachuma on the Santa Ynez River. This equation, ETmax = 2.14(1-C)P0.647+ 4.53 C1.76 P0.68, is recommended for estimating maximum annual ET demand for conceptual models. Where C is not known, the upper limit of ET = f (P) may be approximated by use of 0.65 for C. The equation has been derived for large unmanaged watersheds. Applicability for evaluation of contemporary multiple purpose vegetation management should be determined by studies of the hydrology of small openings in shrub and tree cover.  相似文献   

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

18.
The 2001 National Land Cover Database (NLCD) provides 30-m resolution estimates of percentage tree canopy and percentage impervious cover for the conterminous United States. Previous estimates that compared NLCD tree canopy and impervious cover estimates with photo-interpreted cover estimates within selected counties and places revealed that NLCD underestimates tree and impervious cover. Based on these previous results, a wall-to-wall comprehensive national analysis was conducted to determine if and how NLCD derived estimates of tree and impervious cover varies from photo-interpreted values across the conterminous United States. Results of this analysis reveal that NLCD significantly underestimates tree cover in 64 of the 65 zones used to create the NCLD cover maps, with a national average underestimation of 9.7% (standard error (SE) = 1.0%) and a maximum underestimation of 28.4% in mapping zone 3. Impervious cover was also underestimated in 44 zones with an average underestimation of 1.4% (SE = 0.4%) and a maximum underestimation of 5.7% in mapping zone 56. Understanding the degree of underestimation by mapping zone can lead to better estimates of tree and impervious cover and a better understanding of the potential limitations associated with NLCD cover estimates.  相似文献   

19.
The hydrologic response to statistically downscaled general circulation model simulations of daily surface climate and land cover through 2099 was assessed for the Apalachicola‐Chattahoochee‐Flint River Basin located in the southeastern United States. Projections of climate, urbanization, vegetation, and surface‐depression storage capacity were used as inputs to the Precipitation‐Runoff Modeling System to simulate projected impacts on hydrologic response. Surface runoff substantially increased when land cover change was applied. However, once the surface depression storage was added to mitigate the land cover change and increases of surface runoff (due to urbanization), the groundwater flow component then increased. For hydrologic studies that include projections of land cover change (urbanization in particular), any analysis of runoff beyond the change in total runoff should include effects of stormwater management practices as these features affect flow timing and magnitude and may be useful in mitigating land cover change impacts on streamflow. Potential changes in water availability and how biota may respond to changes in flow regime in response to climate and land cover change may prove challenging for managers attempting to balance the needs of future development and the environment. However, these models are still useful for assessing the relative impacts of climate and land cover change and for evaluating tradeoffs when managing to mitigate different stressors.  相似文献   

20.
Claggett, Peter R., Judy A. Okay, and Stephen V. Stehman, 2010. Monitoring Regional Riparian Forest Cover Change Using Stratified Sampling and Multiresolution Imagery. Journal of the American Water Resources Association (JAWRA) 46(2):334-343. DOI: 10.1111/j.1752-1688.2010.00424.x Abstract: The Chesapeake Bay watershed encompasses 165,760 km2 of land area with 464,098 km of rivers and streams. As part of the Chesapeake Bay restoration effort, state and federal partners have committed to restoring 26,000 miles (41,843 km) of riparian forest buffers. Monitoring trends in riparian forest buffers over large areas is necessary to evaluate the efficacy of these restoration efforts. A sampling approach for estimating change in riparian forest cover from 1993/1994 to 2005 was developed and implemented in Anne Arundel County, Maryland, to exemplify a method that could be applied throughout the Bay watershed. All stream reaches in the county were stratified using forest cover change derived from Landsat imagery. A stratified random sample of 219 reaches was selected and forest cover change within the riparian buffer of each sampled reach was interpreted from high-resolution aerial photography. The estimated footprint of gross change in riparian forest cover (i.e., the sum of gross gain and gross loss) for the county was 1.83% (SE = 0.22%). Stratified sampling taking advantage of a priori knowledge of locations of change proved to be a practical and efficient protocol for estimating riparian forest buffer change at the county scale and the protocol would readily extend to much broader scale monitoring.  相似文献   

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

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