首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
The dynamics of the Lesina coastal lagoon (Italy) in terms of nutrients, phytoplankton and chemical–physical parameters were evaluated, together with their functional relationships with freshwater inputs, in order to identify ecosystem responses to changes in driving forces in a Mediterranean non-tidal lentic environment. Lesina Lagoon is a shallow coastal environment characterised by limited exchange with coastal waters, which favours enrichment of nutrients and organic matter and benthic fluxes within the system. Lagoon–sea exchanges are influenced by human management. There is a steep salinity gradient from East to West. High nitrogen and silica values were found close to freshwater inputs, indicating wastewater discharges and agricultural runoff, especially in winter. Dissolved oxygen was well below saturation (65%) near sewage and runoff inputs in the western part of the lagoon during summer. Classification in accordance with EEA (2001) guidelines suggests the system is of “poor” or “bad” quality in terms of nitrogen concentrations in the eastern zone during the winter rainy period. In terms of phosphate concentrations, the majority of the stations fall into the “good” category, with only two stations (close to the sewage and runoff inputs) classed as “bad”. In both cases, the raw nitrogen levels make the lagoon a P-limited system, especially in the eastern part. There was wide space–time variability in chlorophyll a concentrations, which ranged from 0.25 to 56 μg l−1. No relationships between chlorophyll a and nutrients were found, suggesting that autotrophic biomass may be controlled by a large number of internal and external forcing factors driving eutrophication processes. Water quality for this type of environment depends heavily on pressure from human activities but also on the management of sewage treatment plants, agricultural practices and the channels connecting the lagoon with the sea.  相似文献   

2.
In tropical and subtropical zones, coastal lagoons are surrounded by mangrove communities which are a source of high quantity organic matter that enters the aquatic system through litter fall. This organic matter decomposes, becoming a source of nutrients and other substances such as tannins, fulvic acids and humic acids that may affect the composition and productivity of phytoplankton communities. Sontecomapan is a coastal lagoon located in the southern Gulf of Mexico, which receives abundant litter fall from mangrove. To study the phytoplankton composition and its variation in this lagoon from October 2002 to October 2003, we evaluated the concentrations of dissolved folin phenol active substances (FPAS) as a measure of plant organic matter, salinity, temperature, pH, O2, N-NH4+, N-NO3, P-PO43−, Si-SiO2, and phytoplanktonic cell density in different mangrove influence zones including the three main rivers that feed the lagoon. Nutrients concentrations depended on freshwater from rivers, however these varied seasonally. Concentrations of P-PO43−, N-NH4+ and FPAS were the highest in the dry season, when maximum mangrove litter fall is reported. Variation of these nutrients seemed to depend on the internal biogeochemical processes of the lagoon. Blooms of diatoms (Skeletonema spp., Cyclotella spp. and Chaetoceros holsaticus) and dinoflagellates (Peridinium aff. quinquecorne, Prorocentrum cordatum) occurred seasonally and in the different mangrove influence zones. The high cell densities in these zones and the occurrence of certain species and its ordination along gradient of FPAS in a canonical correspondence analysis, suggest that plant organic matter (i.e. mangrove influence) may contribute to phytoplankton dynamics in Sontecomapan lagoon.  相似文献   

3.
Benthic oxygen profiles were acquired using microsensors over two seasonal cycles (December 2001, April and August 2002, January and May 2003) at two stations differently affected by shellfish farming activity in the Thau lagoon (French Mediterranean coast). This study was part of the Microbent-PNEC Program on the study of biogeochemical processes at the sediment–water interface in an eutrophicated environment. We explored seasonal and spatial heterogeneity as well as the biogeochemical drivers of oxygen uptake, such as in situ temperature, bottom water oxygenation and organic matter deposition. O2 consumption rates were determined by using a transport-reaction model. Maximum rates were reached in August and May and minimum rates in December, April or January. The effect of oyster farming on oxygen fluxes was clearly identified with higher diffusive oxygen uptake in the station inside the oyster parks (C5; 36.8 ± 18.5–87.7 ± 40.8 mmol m−2 d−1), compared with the station lying outside the oyster parks (C4; 8.6 ± 2.1–30.7 ± 8.3 mmol m−2 d−1). At C5, the large spatial heterogeneity was statistically concealing temporal variation, whereas a clear statistical difference between cold and warm periods appeared at C4. In these lagoon sediments, the seasonal dynamics of diffusive oxygen demand and consumption rates were mainly driven by seasonal temperature variation at both stations, as well as by seasonal organic matter delivery to the sediment at the station located outside the oyster parks. In the station located below the oyster parks, seasonal variation of organic matter deposition was dampened by oyster filtering activity. Seasonal temperature variation thus appeared as the major driver of oxygen dynamics in this station. Measurements of total O2 uptake rates indicated a significant fraction of microbial recycling and diffusive transport in oxygen uptake at the station located close to the oyster parks. In the open water site, fauna-mediated O2 transport prevailed in April 2002 (cold conditions), whereas the microbial recycling seemed to dominate in May 2003 (warm conditions).  相似文献   

4.
Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.  相似文献   

5.
Brackish water ecosystems are often exposed to wide variations in environmental variables, including temperature and salinity, which may cause strong selective pressures on organisms modifying the genetic patterns of species. The aim of this work was to test whether there is a ‘divergence-with-gene flow’ in coastal lagoon populations of white seabream (Diplodus sargus) (Ria Formosa, S Portugal and Mar Menor, SE Spain) respect to four marine populations, by using partial sequences of cyt b mitochondrial gene and information from nine microsatellite loci. Genetic diversity was highest in both coastal lagoons (Mar Menor and Ria Formosa) considering mitochondrial and nuclear markers. Although some of FST population pairwise comparisons were not significant, analyses of molecular variance (AMOVAs) detected differences between groups (coastal lagoon and marine) close to significance. Also, only two haplotypes (Cytb-17 and Cytb-18) were detected in both coastal lagoon sampling sites and these localities (Mar Menor and Ria Formosa) showed the highest number of singletons, some of them with a high number of mutations, as has been already described for other Mar Menor populations (Pomatochistus marmoratus and Holothuria polii). Also, several tests detected significant positive and balancing selection considering mtDNA and microsatellite data. These data support the hypothesis of selection as one of the drivers of the genetic differences found between coastal lagoon and marine populations. The life strategy adopted by Diplodus sargus in coastal lagoons allows it to decrease its mortality rate and improve the heritability of its genes. Also, the increase time spent in coastal lagoons with different temperatures and salinities favours the fitness selection and the maintenance of exclusive haplotypes and genotypes in coastal lagoon inhabitants favouring the ‘divergence-with-gene-flow’.  相似文献   

6.
7.
Drana Lagoon, located at the NW site of Evros River Delta, was drained in 1987 and re-flooded in 2004 within the framework of an integrated wetland restoration project. This study presents the results of a monitoring program of the lagoon's oceanographic, water quality and fish fauna characteristics, during the pre- and post-restoration period. Results depict the presence of high salinity water (up to 41) due to seawater intrusion, strong evaporation in its interior and inadequate freshwater inflows. Overall, nutrient levels were low depicting local changes. Tidal variability at the mouth was approximately 0.2 m, producing high velocity tidal currents (up to 0.75 m/s). Eleven fish fauna species were collected; seven species were caught in both the inlet channel and the lagoon during the pre-restoration period and nine species in the post-restoration period. Atherina boyeri (37.6%) and Pomatoschistus marmoratus (31.7%) dominated the lagoon during the post-restoration period. Most of the A. boyeri specimens (88.5%) were caught inside the lagoon, while P. marmoratus had an almost equal distribution in the inlet channel and the lagoon (56.3% and 43.7% respectively). The presence of species of the Mugilidae family (5.2% total average catches after lagoon re-flooding) was mainly in the inlet channel (12.6% of the average catches) and not inside the lagoon (only 1.3% of the average catches). The small number of fish species inhabiting the lagoon might be the result of the recent restoration or it could be related with the increased water flow observed at the lagoon mouth during the flood and ebb tidal phases, and also in the presence of a smooth bank in the concrete waterspout that connects the entrance channel with the lagoon. The limited presence of the Mugilidae juveniles inside the lagoon could be related to the prevailing tidal inlet dynamics (i.e. strong ebb flow at lagoon inlet), thus preventing the species to enter the lagoon. In order to restore the lagoon environment, careful and gradual steps should be undertaken under the basis of continuous monitoring of hydrologic, environmental and fisheries system's status.  相似文献   

8.
Tidal measurements and a depth-averaged 2D model are used to examine wave progression and circulation in a long, shallow, micro-tidal lagoon in Sri Lanka. Ranges and phase lags for different tidal constituents are used to calibrate the model. A single drag coefficient, Cd = 0.0032, gives almost perfect agreement with data. Current measurements are used for validation of the model. The lagoon tide consists of a combination of progressive and standing waves, where progressive waves dominate in the outer part and standing waves in the inner. A Lagrangian based particle-tracking method is developed to study tidally and wind induced residence times. If tides were the only factor affecting the residual circulation, the residence time inside the narrowest section would be approximately 100 days. Steady winds (of typical monsoon average) decrease the residence times to 60–90 days. Estuarine forcing due to net freshwater supply is not modelled (due to lack of reliable runoff data), but independent, long-term salinity observations and calculations based on volume and salt conservation during periods of negligible freshwater supply (the lagoon is seasonally hypersaline) indicate residence times ranging from 40 to 80 days. Model derived residence times based on tides alone represent a minimum exchange. Even weak forcing, through winds, excess evaporation or freshwater supply efficiently reduces residence times.  相似文献   

9.
Abstract. Three years of observations on a population of Notomastus latericeus Sars from Acquatina lagoon (Lecce, Italy) are reported. The population dynamics, life cycle, and reproduction were investigated in order to explain periodic density variations of this species. Some physical parameters of the water column were measured and sediment analyses were performed to characterize the environment. Secondary production was also estimated as part of a functional study of the benthic system in the Acquatina lagoon. The species is opportunistic; its density is apparently correlated to variation in salinity, but its presence in the lagoon is dependent on recruitment success and competition with other polychaete species.  相似文献   

10.
A one-dimensional model of temperature, salinity, nutrients, oxygen, carbon, and argon chemistry is used to hindcast the annual cycle of sea surface pCO2 at weathership Station Papa in the subarctic Pacific (50°N, 145°W), based on recent biological and chemical measurements made in conjunction with the SUPER program. Heat fluxes are calculated from the meteorological time series data from the Canadian weathership program. The Price, Weller and Pinkle (1986) model is used for predicting mixed layer dynamics. The rate of new production in the model is based on sediment trap data (Welschmeyer, personal communication) and a comparison of model oxygen and argon concentrations with the recent data of Emerson, Quay, Stump, Wilbur and Knox (1991). The balances of nutrients and oxygen in the permanent halocline require isopycnal ventilation on a time frame of approximately 10 years; this estimate is consistent with estimate of Van Scoy, Fine and Ostlund (1991) based in tritium data from Geosecs and Long Lines programs. The model results are compared with the 5 year time series data presented by Wong and Chang (1990). The model appears to capture the mean sea surface pCO2 and the magnitude and timing of the seasonal cycle. The data, howeber, contain much greater high frequency variation than the model results, which we believe is caused by patchiness in the horizontal distribution of NO3. The model pCO2 sensitivity to the chemistry of the upwelling water and the rate of biological new production are presented.Although the model simulation of pCO2 is satisfactory, this study illustrates the limitations of modelling the chemistry of the upper ocean in one dimension. The slow currents and horizontally homogeneous sea surface in the subarctic Pacific make Papa one of the best available candidates for modelling in 1-D. In spite of this, a 1-D formulation is inadequate to address the chemistry of the halocline (a crucial lower boundary condition to the mixed layer) and does not constrain the transport of the nutrients by wind-driven currents in the mixed layer. We conclude that further progress in modelling the upper ocean nutrient and carbon cycles will requires simulation in three dimensions.  相似文献   

11.
Factors controlling the CO2 system parameters, including the partial pressure of CO2 (PCO2) in coral reef waters, were investigated in three mid-oceanic reefs of the Indo-Pacific region. Surface water PCO2 in the lagoons of Majuro Atoll and Palau barrier reef in the Pacific were 25 µatm and 48 µatm higher than those of the offshore waters, respectively, while South Male Atoll lagoon of the Maldives in the Indian Ocean exhibited relatively small difference in PCO2 compared to the offshore water. Observations from Majuro Atoll and Palau barrier reef are consistent with the view that calcium carbonate production predominates in coral reefs. On the other hand, results from South Male Atoll can be attributed to the thorough flushing of the lagoon, which is connected to the open ocean by numerous deep channels. The offshore-lagoon PCO2 difference depends on system-level net organic-to-inorganic carbon production ratio while reef topography, especially residence time of the lagoon, has a secondary effect on the magnitude of the offshore-lagoon difference. A potential for releasing CO2 might be more evident in an enclosed atoll where the reef water has a longer residence time. Oceanic atoll and barrier reef lagoons, which are in the terminal stage of evolutionary history of oceanic volcanic islands, have the potential to release CO2 to the atmosphere.  相似文献   

12.
Using data taken from three long-term monitoring programs, we modelled the population dynamics of 13 common macrozoobenthic species of the Balgzand, an approximately 50 km2 intertidal area in the Dutch part of the Wadden Sea. In order to identify likely interactions among species, while accounting as much as possible for other factors influencing population dynamics, models included both ‘environmental’ (water temperature, and phytoplankton concentration) and ‘biological’ (biomass of potentially interacting species) variables. This approach appeared to be effective at identifying certain types of interactions acting over relatively short time-frames, such as predation and competition for food. Among the species we considered, the strongest effects of this kind appeared to be from what we presume to be predation by Nephtys hombergii, and from competition for food among bivalves. Several of the strongest results of our analysis corroborated expectations derived from previous small-scale surveys or experimental studies, but there were also a number of highly significant results indicating possible interactions which were not immediately explicable, including the apparent positive effect of Macoma balthica on Arenicola marina, and the apparent negative effect of Nephtys hombergii on Mya arenaria. Our results also suggest that year to year changes in Cerastoderma edule and Arenicola marina populations did not have substantial and widespread influences on the population dynamics of the majority of other common infauna, which is not to deny the likely critical importance of recruitment effects and biogenic transformation of the intertidal environment by these species over longer time frames. Our modelling approach could be applied to other long-term data sets, and will become more useful as longer time series become available.  相似文献   

13.
《Oceanologica Acta》1998,21(6):819-829
Measurements of oxygen fluxes at the water-sediment interface (transparent and opaque enclosures) were made on bare sediments inhabited by microphytes on the one hand, and on sediments diversely colonized by macrophytes (macroalgae and seagrasses) on the other hand. Five stations, typical of different biotopes of a Mediterranean shallow lagoon were studied from May to July 1993 in order to observe an anoxic event which usually occurs at that time of the year. Average diurnal respiration of benthic communities ranged from 40 mg (or 1.25 mmol) O2 m2 h−1 in bare sediments (31 % of the lagoon area) to 76 mg (or 2.37 mmol) O2 m−2 h−1 in sediments with a medium coverage of macroalgae (37 % of the total area) and, finally, to 100 mg (or 3.12 mmol) O2 m−2h−1 in the denser macrophytic area (32 % of the total area). The highest diurnal gross production was observed in the zone colonized by macroalgae and seagrasses, and especially in corridors between shellfish-cultivation tables (300 mg O2 m−2 h−1, or 9.37 mmol O2m−2 h−1, equivalent to 113 mg C m−2 h−1). Overall, during this period, net bottom oxygen production was close to nil in ca. 30 %, and positive in 70 % of the lagoon area. The average net oxygen production for the whole lagoon in summer is thus in the order of 100 mg O2 m−2h−1. In 1993, at the end of July, an anoxic event was avoided due to a period of strong wind.  相似文献   

14.
The aim of the present study was to investigate seasonal and spatial patterns of soil oxygen consumption, nitrification, denitrification and fluxes of dissolved inorganic nitrogen (DIN) in a tidal salt marsh of the Lagoon of Venice, Italy. In the salt marsh, intact soil cores including overlying water were collected monthly at high tide from April to October in salt marsh creeks and in areas covered by the dominant vegetation, Limonium serotinum. In May, cores were also collected in areas with vegetation dominated by Juncus maritimus and Halimione portulacoides. In laboratory incubations at in situ temperature in the dark, flux rates of oxygen and DIN were monitored in the overlying water of the intact cores. 15N-nitrate was added to the overlying water and nitrification and denitrification were measured using isotope-dilution and -pairing techniques. The results show that highest soil oxygen consumption coincided with the highest water temperature in June and July. The highest denitrification rates were recorded in spring and autumn coinciding with the highest nitrate concentrations. Soil oxygen consumption and nitrification rates differed between sampling sites, but denitrification rates were similar among the different vegetation types. The highest rates were recorded in areas covered with L. serotinum. Burrowing soil macrofauna enhanced oxygen consumption, nitrification and denitrification in April and May. The data presented in this study indicate high temporal as well as spatial variations in the flux of oxygen and DIN, and nitrogen transformations in the tidal salt marshes of the Venice lagoon during the growth season. The results identify the salt marshes of the Venice lagoon as being metabolically very active ecosystems with a high capacity to process nitrogen.  相似文献   

15.
Long-term macrobenthos data from Kiel Bight in the Western Baltic collected between 1968 and 2000 have been correlated with the winter NAO index (North Atlantic Oscillation Index) and other environmental data such as temperature, salinity and oxygen content in the bottom water in order to detect systematic patterns related to so far unexplained abiotic signals in the dynamics of zoobenthic species assemblages. The benthos data come from a cluster of five stations (Süderfahrt/ Millionenviertel) in Kiel Bay. Our investigations concentrated on the macrobenthic dynamics with a focus on the number of species m 2 (species richness). Using logarithms and the time series analysis approach of Box/Jenkins (ARIMA modelling, transfer function modelling) it was shown that species richness was strongly influenced by the winter NAO (adjusted for a linear time trend within the 1968-2000 period) and salinity (with a shift/lag of four years). Bootstrapping experiments (i.e. sampling from the error process) and analysis of prediction power (by means of the one- or more-years leaving-out method) showed that the parameter estimates behaved in a stable way, leading to a relatively robust model.  相似文献   

16.
Hydrodynamic modeling can be used to spatially characterize water renewal rates in coastal ecosystems. Using a hydrodynamic model implemented over the semi-enclosed Southwest coral lagoon of New Caledonia, a recent study computed the flushing lag as the minimum time required for a particle coming from outside the lagoon (open ocean) to reach a specific station [Jouon, A., Douillet, P., Ouillon, S., Fraunié, P., 2006. Calculations of hydrodynamic time parameters in a semi-opened coastal zone using a 3D hydrodynamic model. Continental Shelf Research 26, 1395–1415]. Local e-flushing time was calculated as the time requested to reach a local grid mesh concentration of 1/e from the precedent step. Here we present an attempt to connect physical forcing to biogeochemical functioning of this coastal ecosystem. An array of stations, located in the lagoonal channel as well as in several bays under anthropogenic influence, was sampled during three cruises. We then tested the statistical relationships between the distribution of flushing indices and those of biological and chemical variables. Among the variables tested, silicate, chlorophyll a and bacterial biomass production present the highest correlations with flushing indices. Correlations are higher with local e-flushing times than with flushing lags or the sum of these two indices. In the bays, these variables often deviate from the relationships determined in the main lagoon channel. In the three bays receiving significant riverine inputs, silicate is well above the regression line, whereas data from the bay receiving almost insignificant freshwater inputs generally fit the lagoon channel regressions. Moreover, in the three bays receiving important urban and industrial effluents, chlorophyll a and bacterial production of biomass generally display values exceeding the lagoon channel regression trends whereas in the bay under moderate anthropogenic influence values follow the regressions obtained in the lagoon channel. The South West lagoon of New Caledonia can hence be viewed as a coastal mesotrophic ecosystem that is flushed by oligotrophic oceanic waters which subsequently replace the lagoonal waters with water considerably impoverished in resources for microbial growth. This flushing was high enough during the periods of study to influence the distribution of phytoplankton biomass, bacterial production of biomass and silicate concentrations in the lagoon channel as well as in some of the bay areas.  相似文献   

17.
Felice Arena  Diego Pavone   《Ocean Modelling》2009,26(3-4):217-225
This paper deals with the long-term modelling of high sea waves. The solution is given for the return period of sea storms during which an arbitrary chosen number of waves, with crest-to-trough heights exceeding a fixed threshold, occur. This return period is derived starting from the Equivalent Triangular Storm (ETS) model, which associates a triangle to each actual storm and thus represents a significant wave height time series at a fixed location by means of a sequence of triangular storms. The short-term statistics is then applied to investigate the occurrence of large crest-to-trough wave heights during a given storm. Finally, by combining the statistical distribution of significant wave heights, the ETS model and the short-term wave statistics, the solution is given for the return periods RN and RN of a sea storm in which N or at least N waves higher than a fixed threshold occur. The values of RN are then calculated, starting from data of two buoys moored in the Pacific Ocean and in the Mediterranean Sea.  相似文献   

18.
There are gaps in what is known about the patterns of gastropod shell use by hermit crabs in Brazilian lagoon areas, especially in Northeastern Brazil. However, this is important because the understanding of selection patterns provides information on life history and eco‐evolutionary conceptions of paguroids. The present study investigated the use of gastropod shells occupied by Pagurus criniticornis and Clibanarius sclopetarius in a coastal lagoon seasonally connected to an estuary and to the sea, correlating eco‐evolutionary aspects and hydrodynamic characters. The study was carried out between February 2013 and January 2014. Hermit crabs and their shells were identified and measured (hermit crabs’ cephalothorax shield length and width, major propodus length and height, weight, shell total length and width, shell aperture length and width). Partial least squares regression was used to analyse the morphometric data. Additionally, a multinomial proportions test was performed to infer patterns (inter‐/intra‐specific) of shells’ occupation. Hermit crab species occupied a total of 13 types of shells, predominantly those of Neritina virginea (67.83%). The shell weight was the most important determinant of shell occupancy in the morphometric model (variable importance in projection >1). The proportions of N. virginea shells used were similar in both species of hermit crabs, except between the sexes of P. criniticornis. Presumably, the high utilization rate of N. virginea shells is related to its abundance and dispersal, and to the shells’ suitability for hermit crabs, which is reflected by the morphometric model, as well as by the hydrodynamics of the lagoon‐estuarine environment. The relative occupation of different types of shells for each species of hermit crabs studied appears to be associated with regulation by inter‐/intra‐specific competition, which fosters the co‐occurrence of those populations. These results endorse coastal lagoons as a refuge and recruitment area for aquatic fauna. In addition, it highlights a challenge in the management and conservation of paguroid species whose population dynamics depend on the resources coming from the Gastropoda community.  相似文献   

19.
A modelling approach is proposed to evaluate the environmental dynamics of coastal lagoons. The water, heat and salt balances are addressed simultaneously, providing a better estimation of evaporation and water exchanges. Compared to traditional approaches, the model presented accounts for the effects of water salinity, heat storage and net energy advected in the water body. The model was applied daily to the Mar Menor coastal lagoon (SE Spain) from 2003 through 2006. Water exchanges with the Mediterranean Sea were estimated based on the monthly trend of the lagoon salinity and were correlated with monthly averages of wind speed. The mean daily water exchange with the sea was 1.77 hm3 d−1. This exchange accounted for only 1% of the heat losses in the lagoon heat balance, and it is the most important flow in the water balance. The mean annual evaporation flux amounted to 101.3 W m−2 (3.55 mm d−1), while the sensible heat flux amounted to 19.7 W m−2, leading to an annual Bowen ratio on the order of 0.19. To validate the model, daily water temperatures were predicted based on the daily heat balance of the water body and were compared with remote sensing data from water surface standard products.  相似文献   

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

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