首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 296 毫秒
1.
The biological integrity of stream ecosystems depends critically on human activities that affect land use/cover along stream margins and possibly throughout the catchment. We evaluated stream condition using an Index of Biotic Integrity (IBI) and a habitat index (HI), and compared these measures to landscape and riparian conditions assessed at different spatial scales in a largely agricultural Midwestern watershed. Our goal was to determine whether land use/cover was an effective predictor of stream integrity, and if so, at what spatial scale. Twenty-three sites in first-through third-order headwater streams were surveyed by electrofishing and site IBIs were calculated based on ten metrics of the fish collection. Habitat features were characterized through field observation, and site HIs calculated from nine instream and bank metrics. Field surveys, aerial photograph interpretation, and geographic information system (GIS) analyses provided assessments of forested land and other vegetation covers at the local, reach, and regional (catchment) scales. The range of conditions among the 23 sites varied from poor to very good based on IBI and HI scores, and habitat and fish assemblage measures were highly correlated. Stream biotic integrity and habitat quality were negatively correlated with the extent of agriculture and positively correlated with extent of wetlands and forest. Correlations were strongest at the catchment scale (IBI with % area as agriculture, r2=0.50, HI with agriculture, r2=0.76), and tended to become weak and non-significant at local scales. Local riparian vegetation was a weak secondary predictor of stream integrity. In this watershed, regional land use is the primary determinant of stream conditions, able to overwhelm the ability of local site vegetation to support high-quality habitat and biotic communities.  相似文献   

2.
We investigated the effects of urban environments on the chemical properties of forest soils in the metropolitan areas of Baltimore, New York, and Budapest. We hypothesized that soils in forest patches in each city will exhibit changes in chemistry corresponding to urbanization gradients, but more strongly with various urban metrics than distance to the urban core. Moreover, differences in parent material and development patterns would differentially affect the soil chemical response in each metropolitan area. Results showed that soil chemical properties varied with measures of urban land use in all three cities, including distance to the urban core, which was an unexpected result. Moreover, the results showed that the spatial extent and amount of change was greater in New York than in Baltimore and Budapest for those elements that showed a relationship to the urbanization gradient (Pb, Cu, and to a lesser extent Ca). The spatial relationship of the soil chemical properties to distance varied from city to city. In New York, concentrations of Pb, Cu, and Ca decreased to approximately background concentrations at 75 km from the urban core. By contrast, concentrations of these elements decreased closer to the urban core in Baltimore and Budapest. Moreover, a threshold was reached at about 75% urban land use above which concentrations of Pb and Cu increased by more than twofold relative to concentrations below this threshold. Results of this study suggest that forest soils are responding to urbanization gradients in all three cities, though characteristics of each city (spatial pattern of development, parent material, and pollution sources) influenced the soil chemical response.  相似文献   

3.
Since the mid eighties, agricultural development and increased population growth in Vietnam’s northern highlands have modified land use patterns and thus, increased the runoff process and soil degradation induced by water erosion. In the last decade, Vietnamese literature has focused on the computation of soil losses over large areas. Most of these spatial and quantitative soil erosion studies do not consider the impact of agricultural land use diversity (spatial heterogeneity), particularly at the watershed scale, and the annual variability of seasonal landscape factors on soil erosion vulnerability and hence, landscape dynamics. We present an integrated approach combining field measurements and observations, GIS and modeling to determine the spatial and temporal dynamics of soil erosion vulnerability according to watershed units and hence, the impact of physical environment components and agricultural land use patterns on landscape evolution. Tables and graphics showing the cropping systems, the periods within a year, and the watershed units that are most vulnerable are presented. The double cultivation cycles for paddy rice fields not only imply two periods of land preparation and establishment that expose the soil surface to raindrop impacts, but also increased soil management practices that decrease the soil’s resistance to detachment. Despite the low levels of soil management practices for the shifting cultivation system, the near absence of soil conservation practices clearly increases their vulnerability. Hence, rainfed cropping systems, mainly soya and cassava, cultivated on sloping lands (hills and mountains) where soil erosion vulnerability is the highest represent the watershed units which are the most prone to soil loss.  相似文献   

4.
Influences of upland and riparian land use patterns on stream biotic integrity   总被引:11,自引:2,他引:11  
Snyder  C.D.  Young  J.A.  Villella  R.  Lemarié  D.P. 《Landscape Ecology》2003,18(7):647-664
We explored land use, fish assemblage structure, and stream habitat associations in 20 catchments in Opequon Creek watershed, West Virginia. The purpose was to determine the relative importance of urban and agriculture land use on stream biotic integrity, and to evaluate the spatial scale (i.e., whole-catchment vs riparian buffer) at which land use effects were most pronounced. We found that index of biological integrity (IBI) scores were strongly associated with extent of urban land use in individual catchments. Sites that received ratings of poor or very poor based on IBI scores had > 7% of urban land use in their respective catchments. Habitat correlations suggested that urban land use disrupted flow regime, reduced water quality, and altered stream channels. In contrast, we found no meaningful relationship between agricultural land use and IBI at either whole-catchment or riparian scales despite strong correlations between percent agriculture and several important stream habitat measures, including nitrate concentrations, proportion of fine sediments in riffles, and the abundance of fish cover. We also found that variation in gradient (channel slope) influenced responses of fish assemblages to land use. Urban land use was more disruptive to biological integrity in catchments with steeper channel slopes. Based on comparisons of our results in the topographically diverse Opequon Creek watershed with results from watersheds in flatter terrains, we hypothesize that the potential for riparian forests to mitigate effects of deleterious land uses in upland portions of the watershed is inversely related to gradient.This revised version was published online in May 2005 with corrections to the Cover Date.  相似文献   

5.
This study quantified nitrogen mineralization and nitrification potentials in soils of hardwood forests of southern Ohio at three spatial scales: (1) the regional scale, represented by four study areas of 90–120 ha separated by 3–65 km, (2) the local scale, represented by three contiguous watersheds within each study area, and (3) the topographic scale, represented by xeric, intermediate, and mesic sites within each watershed, as defined by a GIS-generated Integrated Moisture Index (IMI). Organic C, NO3 pool size, net N mineralization, proportional nitrification, and net nitrification potentials all varied among study sites (i.e. at the regional scale). Using path analysis, we were able to construct scale-independent causal models explaining 30–35% of the variance in organic C and potential net N mineralization and 70% of the variance in potential net NO3 production. Site- and scale-specific differences in geology and/or land use history among study sites were likely responsible for the variation not explained by the path analysis. At the local scale, there were significant variations in organic C and inorganic N pool sizes among watersheds within a study site in two of the four study sites. In addition, most parameters we measured varied significantly along the topographic gradient (i.e. with long-term soil moisture availability/IMI). Based on our results, scaling up models of nitrification from plot scale to the regional scale should be straightforward, whereas scaling up organic C storage and N mineralization will require incorporation of independent scaling paradigms at three (or more) spatial scales.  相似文献   

6.
Forest structure and composition influence patterns of insect outbreaks and can be explained on the Walker Branch watershed by past land use (timber harvest and agriculture), soils, aspect, and slope. In particular, pine bark beetles caused large losses of pine on sites that had been used for agriculture, on Fullerton silt loam soils, and on north-to-northeast and east-to-southeast exposures. Hickory bark beetles had a high impact on hickory biomass on Bodine soil areas that were forested in 1935 and sloped greater than 11%. Thus, prior land use can have an indirect effect on future disturbances.Because forest disturbances can affect nutrient distribution, land use can also indirectly affect nutrient availability. For example, locations of hickory bark beetle outbreaks experience a large flux of calcium from dead wood to soil because hickory accumulates large amounts of calcium in woody tissue. The research demonstrates a link between past land use, insect outbreaks, and calcium cycling.  相似文献   

7.
Lobo  Agustín  Moloney  Kirk  Chic  Oscar  Chiariello  Nona 《Landscape Ecology》1998,13(2):111-131
An important practical problem in the analysis of spatial pattern in ecological systems is that requires spatially-intensive data, with both fine resolution and large extent. Such information is often difficult to obtain from field-measured variables. Digital imagery can offer a valuable, alternative source of information in the analysis of ecological pattern. In the present paper, we use remotely-sensed imagery to provide a link between field-based information and spatially-explicit modeling of ecological processes. We analyzed one digitized color infrared aerial photograph of a serpentine grassland to develop a detailed digital map of land cover categories (31.24 m × 50.04 m of extent and 135 mm of resolution), and an image of vegetation index (proportional to the amount of green biomass cover in the field). We conducted a variogram analysis of the spatial pattern of both field-measured (microtopography, soil depth) and image-derived (land cover map, vegetation index, gopher disturbance) landscape variables, and used a statistical simulation method to produce random realizations of the image of vegetation index based upon our characterization of its spatial structure. The analysis revealed strong relationships in the spatial distribution of the ecological variables (e.g., gopher mounds and perennial grasses are found primarily on deeper soils) and a non-fractal nested spatial pattern in the distribution of green biomass as measured by the vegetation index. The spatial pattern of the vegetation index was composed of three basic components: an exponential trend from 0 m to 4 m, which is related to local ecological processes, a linear trend at broader scales, which is related to a general change in topography across the study site, and a superimposed periodic structure, which is related to the regular spacing of deeper soils within the study site. Simulations of the image of vegetation index confirmed our interpretation of the variograms. The simulations also illustrated the limits of statistical analysis and interpolations based solely on the semivariogram, because they cannot adequately characterize spatial discontinuities.  相似文献   

8.
The spatial distribution of soil carbon (C) is controlled by ecological processes that evolve and interact over a range of spatial scales across the landscape. The relationships between hydrologic and biotic processes and soil C patterns and spatial behavior are still poorly understood. Our objectives were to (i) identify the appropriate spatial scale to observe soil total C (TC) in a subtropical landscape with pronounced hydrologic and biotic variation, and (ii) investigate the spatial behavior and relationships between TC and ecological landscape variables which aggregate various hydrologic and biotic processes. The study was conducted in Florida, USA, characterized by extreme hydrologic (poorly to excessively drained soils), and vegetation/land use gradients ranging from natural uplands and wetlands to intensively managed forest, agricultural, and urban systems. We used semivariogram and landscape indices to compare the spatial dependence structures of TC and 19 ecological landscape variables, identifying similarities and establishing pattern–process relationships. Soil, hydrologic, and biotic ecological variables mirrored the spatial behavior of TC at fine (few kilometers), and coarse (hundreds of kilometers) spatial scales. Specifically, soil available water capacity resembled the spatial dependence structure of TC at escalating scales, supporting a multi-scale soil hydrology-soil C process–pattern relationship in Florida. Our findings suggest two appropriate scales to observe TC, one at a short range (autocorrelation range of 5.6 km), representing local soil-landscape variation, and another at a longer range (119 km), accounting for regional variation. Moreover, our results provide further guidance to measure ecological variables influencing C dynamics.  相似文献   

9.
Zhang  Minghua  Geng  Shu  Ustin  Susan L. 《Landscape Ecology》1998,13(1):37-54
Quantitative agricultural landscape indices are useful to describe functional relationships among climatic conditions, groundwater dynamics, soil properties and agricultural land use for mathematical models. We applied methods of regression statistics, variance component estimation and a Geographical Information System (GIS) to construct indices describing crops and soils and to establish functional relationships among these variables. This paper describes the development of indices and the partitioning of the spatial and temporal variation in groundwater models using the data from Tulare County, California, which was selected as the study area. Indices of ground surface elevation, total crop water demand, soil water infiltration rate, and soil production index explain 91% of the variation in average spring groundwater level. After relating spatial patterns of groundwater use to indices of crop and soil properties, we found that mean groundwater use is positively related to total crop water demand and soil water infiltration rate while the variation in groundwater use was negatively correlated with the crop water demand and soil water infiltration rate and positively related to soil water holding capacity. The spatial variation in groundwater use was largely influenced by crops and soil types while the temporal variation was not. We also found that groundwater use increased exponentially with decreasing annual precipitation for most townships. Based on these associations, groundwater use in each township can be forecast from relative precipitation under current methods of agricultural production. Although groundwater table depth is strongly affected by topography, the statistically significant indices observed in the model clearly show that agricultural land use influences groundwater table depth. These simple relationships can be used by agronomists to make water management decisions and to design alternative cropping systems to sustain agricultural production during periods of surface water shortages.  相似文献   

10.
11.
Widespread degradation of wetlands has motivated the development of tools to evaluate wetland condition. The application of field-based tools over large regions can be prohibitively expensive; however, land cover data may provide a surrogate for intensive assessments, enabling rapid and cost-effective evaluation of wetlands throughout whole regions. Our goal was to determine if land cover data could be used to estimate the biotic integrity of wetlands in Alberta??s Beaverhills watershed. Biotic integrity was measured using both plant- and bird-based indices of biotic integrity (IBIs) in 45 wetlands. Land cover data were extracted from seven nested landscape extents (100?C3,000?m radii) and used to model IBI scores. Strong, significant predictions of IBI scores were achieved using land cover data from every spatial extent, even after factoring out the influence of location to address the spatial autocorrelation of land cover classes. Plant-based IBI scores were best predicted using data from 100?m buffers and bird-based IBI scores were best predicted using data extracted from 500?m buffers. Road cover or density and measures of the proportion of disturbed land were consistent predictors of IBI score, suggesting their universal importance to plant and bird communities. Simplified models using the proportion of undisturbed land were less accurate than more detailed models (reductions in r 2 of 0.31?C0.32). Regardless of the level of detail in land cover classification, our results emphasize the need to optimize landscape extent for the taxonomic group of interest: an issue that is typically poorly articulated in studies reporting on the development of GIS-based assessment methods. Our results also highlight the need to calibrate models in test areas before scaling up, to ensure predictive accuracy.  相似文献   

12.
Efforts to quantify the effects of riparian buffers on watershed nutrient discharges have been confounded by a commonly used analysis, which estimates buffer potential as the percentage of forest or wetland within a fixed distance of streams. Effective landscape metrics must instead be developed based on a clear conceptual model and quantified at the appropriate spatial scale. We develop new metrics for riparian buffers in two stages of increasing functional specificity to ask: (1) Which riparian metrics are more distinct from measures of whole watershed land cover? (2) Do functional riparian metrics provide different information than fixed-distance metrics? (3) How do these patterns vary within and among different physiographic settings? Using publicly available geographic data, we studied 503 watersheds in four different physiographic provinces of the Chesapeake Bay Drainage. In addition to traditional fixed-distance measures, we calculated mean buffer width, gap frequency, and measures of variation in buffer width using both “unconstrained” metrics and “flow-path” metrics constrained by surface topography. There were distinct patterns of relationship between watershed and near-stream land cover in each physiographic province and strong correlations with watershed land cover confounded fixed-distance metrics. Flow-path metrics were more independent of watershed land cover than either fixed-distance or unconstrained measures, but both functional metrics provided greater detail, interpretability, and flexibility than the fixed-distance approach. Potential applications of the new metrics include exploring the potential for land cover patterns to influence water quality, accounting for buffers in statistical nutrient models, quantifying spatial patterns for process-based modeling, and targeting management actions such as buffer restoration.  相似文献   

13.
The effects of time on the evolution of land use intensity and soil nutrients distribution were studied in a reclamation zone of the Yangtze Estuary. Land use types were grouped into five intensity levels according to the extent of human disturbance. We used the “space for time substitution” method to test the impact of time on changes in land use intensity after reclamation and found that land use levels increased quickly within the first 35 years, then slowed. Soil salinity was relatively higher in aquaculture ponds than that in areas with other types of land cover due to the use of saline water from underground and the sea. Soil organic matter, available phosphorous and nitrate nitrogen were relatively high in agricultural fields, while nitrate nitrogen was highly variable in agricultural fields. The variations of all four soil properties in the built-up zone were much higher than those in the other land use groups. The spatial distribution of different nutrients is the combined effect of time and land use post reclamation. The results will provide a sound basis for future land use planning of newly reclaimed land, and for further studies on ecological consequences of salt marsh reclamation.  相似文献   

14.
选择呼和浩特市城区公园绿地、居民小区绿地和道路绿化带为研究对象,采集城市中心和东南西北5个区域共45个样品。通过对不同绿地类型和不同区域土壤6种重金属含量的分析,采用单项污染指数法和内梅罗综合污染指数法对重金属的污染程度进行评价,采用主成分分析法分析重金属污染的来源,以期为呼和浩特市绿地土壤环境保护及其重金属污染防治提供科学依据。结果表明:呼和浩特市城区不同绿地类型土壤中4种重金属的平均含量均超过内蒙古的土壤背景值,重金属Cu和Cr含量分别达到背景值的2.20倍和1.68倍;单项污染指数评价结果显示,不同绿地类型土壤中重金属存在不同污染程度,重金属Cu在道路绿化带和居民小区绿地呈中度污染,在公园绿地呈轻度污染。不同土壤重金属元素污染程度依次为Cu>Cr>Zn>Pb>Ni>As。内梅罗综合污染指数法评价结果表明,道路绿化带呈中度污染,其次是居民小区绿地,污染水平为轻度污染,公园绿地污染程度最轻;城西呈中度污染,其余均为轻度污染。不同绿地类型和不同区域土壤重金属综合污染指数大小为道路绿化带(2.09)>居民小区绿地(1.83)>公园绿地(1.54);城西(2.04)>城南(1.78)>城东(1.77)>城中(1.76)>城北(1.73)。主成分分析结果显示,呼和浩特市绿地土壤不同重金属的来源存在差异,其中Cu、Cr和Ni主要来自交通污染源和人们日常生活废弃物的堆放;As主要来源于自然源。  相似文献   

15.
A method developed to evaluate the cumulative effect of wetland mosaics on water quality was applied to 33 lake watersheds in the seven-county region surrounding Minneapolis-St. Paul, Minnesota. A geographic information system (GIS) was used to record and measure landscape variables derived from aerial photos. Twenty-seven watershed land-use and land-cover variables were reduced to eight principal components which described 85% of the variance among watersheds. Relationships between lake water quality variables and the first six principal components plus an index of lake mixis were analyzed through stepwise multiple regression analysis. A combination of three landscape components (wetland/watershed area, agriculture/wetlands, and forest/soils components) explained 49% of the variance in a trophic state index, even though most of the lakes examined were already highly eutrophic, and thus were influenced by internal loading. The regression equations explained a range of 14 to 76% of the variation in individual water quality variables. Forested land-use was associated with lower lake trophic state, chloride, and lead. High lake trophic state was associated with agricultural land-use and with wetland distance from the lake of interest. The extent of wetlands was associated with low total lead and high color in lakes downstream. Wet meadows or herbaceous, seasonally-flooded wetlands contributed more to lake water color than did cattail marshes.  相似文献   

16.
为了解不同土地利用方式对冀西北坝上地区土壤化学性质的影响,对该地区4种土地利用方式(围封林地、围封采伐地、未围封林地和耕地)下的土壤化学性质进行了分析.结果 表明:耕地的pH在7.5~8.5之间,明显高于其他3种土地利用方式(在6.5~7.5之间)(P<0.05);耕地土壤的有机质、全氮、全磷、全钾及速效氮、磷、钾都明显高于3种林地(P<0.05);2种封育林地表层土壤(0~10 cm)各养分含量明显高于未封育林地(P<0.05),但深层土壤差异不明显.土地利用方式对土壤养分的垂直分布有显著影响,耕地土壤养分,除速效磷外,均表现为先下降后上升的趋势,养分含量变化曲线呈"V"型;3种林地土壤养分随土层深度的增加而降低,拐点均在10~20 cm土层,养分含量变化曲线呈"L"型.由以上结果得出结论,土地利用方式对土壤养分含量及其空间分布具有明显影响,耕地因为人工施肥的原因各养分含量明显高于林地,但土壤有明显的盐碱化趋势;封育可明显提高林地表层土壤养分含量,但对深层土壤影响较小.  相似文献   

17.
There is a pressing need to understand the consequences of human activities, such as land transformations, on watershed ecosystem services. This is a challenging task because different indicators of water quality and yield are expected to vary in their responsiveness to large versus local-scale heterogeneity in land use and land cover (LUC). Here we rely on water quality data collected between 1977 and 2000 from dozens of gauge stations in Puerto Rico together with precipitation data and land cover maps to (1) quantify impacts of spatial heterogeneity in LUC on several water quality indicators; (2) determine the spatial scale at which this heterogeneity influences water quality; and (3) examine how antecedent precipitation modulates these impacts. Our models explained 30–58% of observed variance in water quality metrics. Temporal variation in antecedent precipitation and changes in LUC between measurements periods rather than spatial variation in LUC accounted for the majority of variation in water quality. Urbanization and pasture development generally degraded water quality while agriculture and secondary forest re-growth had mixed impacts. The spatial scale over which LUC influenced water quality differed across indicators. Turbidity and dissolved oxygen (DO) responded to LUC in large-scale watersheds, in-stream nitrogen concentrations to LUC in riparian buffers of large watersheds, and fecal matter content and in-stream phosphorus concentration to LUC at the sub-watershed scale. Stream discharge modulated impacts of LUC on water quality for most of the metrics. Our findings highlight the importance of considering multiple spatial scales for understanding the impacts of human activities on watershed ecosystem services.  相似文献   

18.
Jansson  Åsa  Folke  Carl  Langaas  Sindre 《Landscape Ecology》1998,13(4):249-262
We estimate the nitrogen retention capacity of natural wetlands in the 1.7 million km2 Baltic Sea drainage basin, using a wetland GIS data base. There are approximately 138,000 km2 of wetlands (bogs and fens) in the Baltic Sea drainage basin, corresponding to 8% of the area. The input of nitrogen to natural wetlands from atmospheric deposition was estimated to 55,000–161,000 ton y1. A map of the deposition of both wet and dry nitrogen is presented. The input from the human population was estimated to 255,000 ton y1 in terms of excretory release in processed sewage water. There may also be leakage from forests and agricultural land into the wetlands. Due to lack of data on hydrology and topography, such potential nitrogen sources are not accounted for here. The capacity of the wetlands to retain the atmospheric deposition of nitrogen was estimated to 34,000–99,000 ton y1. The potential retention by wetlands was estimated to 57,000–145,000 ton y1 when the nitrogen input from the human population was added. If drained wetlands were to be restored and their area added to the present wetland area, the nitrogen retention capacity was estimated to increase to 196,000–261,000 ton y1. Our results indicate that existing natural wetlands in the Baltic Sea drainage basin annually can retain an amount of nitrogen which corresponds to about 5–13% of annual total (natural and anthropogenic) nitrogen emissions entering the Baltic Sea. The ecosystem retention service performed by wetlands accounts for a substantial nitrogen removal, thereby reducing the eutrophication of the Baltic Sea.  相似文献   

19.
Ants moderate the supply of critical resources such as water and nutrients in desert environments by changing the physical arrangement of soils during nest construction. We measured water infiltration and soil physical and chemical properties on and off the nests of two ant species (Pogonomyrmex rugosus, Aphaenogaster cockerelli) across five sites at differing landscape positions within a Chihuahuan Desert watershed. Our aim was to test whether the effects of these long-lived ant nests on water infiltration and soil chemistry varied between ant species or across sites within a watershed. Water flow was generally slowest at the site with the highest silt and clay contents, and fastest at the site with sandy soils. Flow was generally greater through ant nest soils than adjacent non-nest soils, and we attribute this to increases in macropores in the nests. However, the effects of both ant nests and species varied among sites. Despite wide variation in soil chemical properties across all sites, ant nests had a consistent effect on soil chemical properties, with higher levels of carbon, nitrogen, sulphur, phosphorus and electrical conductivity on nests compared with non-nest soils. Our results demonstrate that while we can generalise about the effects of ant nests on water flow and nutrient levels, differences in soil type, nest density and ant species across sites are likely to moderate these effects.  相似文献   

20.
There is increasing interest in developing criteria to evaluate the environmental implications of intensive agricultural land use. This implies discriminating between nature and man-made effects upon structural and functional attributes of agroecosystems. Adequate indicators of these combined effects should be cost efficient yet compatible with the core of ecological theory on biodiversity, spatial organization and ecosystem stability. We developed resistance-resilience metrics of plant growth to evaluate the intensity of agricultural use in a temperate irrigated basin in southern Argentina. The metrics are based on an analysis of the components of a temporal series of vegetation indices computed at a low resolution from available globally remote sensed reflectance imagery. We related the developed metrics to the properties of the soils and plant canopies observed at field scale and high-resolution imagery of the basin. Soil depth, soil erosion status and land fragmentation account for large fractions of the variance of the distribution of functional groups of the plant canopies and are also correlated with smaller scale attributes of land vegetation cover. Resistance-resilience indicators constitute a cost-efficient and adequate approach to evaluate the degree of intensification of land agricultural use.  相似文献   

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

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