首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 203 毫秒
1.
An algorithm is presented to retrieve the concentrations of chlorophyll a, suspended pariclulate matter and yellow substance from normalized water-leaving radiances of the Ocean Color and Temperature Sensor (OCTS) of the Advanced Earth Observing Satellite (ADEOS). It is based on a neural network (NN) algorithm, which is used for the rapid inversion of a radiative transfer procedure with the goal of retrieving not only the concentrations of chlorophyll a but also the two other components that determine the water-leaving radiance spectrum. The NN algorithm was tested using the NASA's SeaBAM (SeaWiFS Bio-Optical Mini-Workshop) test data set and applied to ADEOS/OCTS data of the Northwest Pacific in the region off Sanriku, Japan. The root-mean-square error between chlorophyll a concentrations derived from the SeaBAM reflectance data and the chlorophyll a measurements is 0.62. The retrieved chlorophyll a concentrations of the OCTS data were compared with the corresponding distribution obtained by the standard OCTS algorithm. The concentrations and distribution patterns from both algorithms match for open ocean areas. Since there are no standard OCTS products available for yellow substance and suspended matter and no in situ measurements available for validation, the result of the retrieval by the NN for these two variables could only be assessed by a general knowledge of their concentrations and distribution patterns. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

2.
Dimethylsulfide (DMS), chlorophyll a (Chl-a), accessory pigments (fucoxanthin, peridinin and 19-hexanoyloxyfucoxanthin), and bacterial production (BP) were measured in the surface layer (0–100 m) of the subarctic North Pacific, including the Bering Sea, during summer (14 July–5 September, 1997). In surface sewater, the concentrations of DMS and Chl-a varied widely from 1.3 to 13.2 nM (5.1 ± 3.0 nM, mean ± S.D., n = 48) and from 0.1 to 2.4 µg L–1 (0.6 ± 0.6 µg L–1, n = 24), respectively. In the subarctic North Pacific, DMS to Chl-a ratios (DMS/Chl-a) were higher on the eastern side than the western side (p < 0.0001). Below the euphotic zone, DMS/Chl-a ratios were law and the correlation between DMS and Chl-a was relatively strong (r 2 = 0.700, n = 27, p < 0.0001). In the euphotic zone, DMS/Chl-a ratios were higher and the correlation between DMS and Chl-a was weak (r 2 = 0.128, n = 50, p = 0.01). The wide variation in DMS/Chl-a ratios would be at least partially explained by the geographic variation in the taxonomic composition of phytoplankton, because of the negative correlation between DMS/Chl-a and fucoxanthin-to-Chl-a ratios (Fuc/Chl-a) (r 2 = 0.476, n = 26, p = 0.0001). Furthermore, there was a positive correlation between DMS and BP (r 2 = 0.380, n = 19, p = 0.005). This suggests that BP did not represent DMS and dimethylsulfoniopropionate (DMSP) removal by bacterial consumption but rather DMSP degradation to DMS by bacterial enzyme.  相似文献   

3.
Satellite-derived ocean color data of Coastal Zone Color Scanner (CZCS) on board the Nimbus-7 and Ocean Color and Temperature Scanner (OCTS) on board the Advanced Earth Observing Satellite (ADEOS) are jointly used with historical in situ data to examine seasonal and spatial distributions of chlorophyll a (Chl-a) and suspended particulate matter (SPM) concentrations in the East China Sea. Ocean color imagery showed that Chl-a concentrations on the continental shelf were higher than those of the Kuroshio area throughout the year. Satellite-derived Chl-a concentrations are generally in good accordance with historical in situ values during spring through autumn (although no shipboard in situ measurement was conducted at nearshore areas). In contrast, ocean color imagery in winter indicated high Chl-a concentrations (4–10 mg m–3) on the continental shelf where bottom depth was less than 50 m when surface water was turbid (2–72 g m–3 of SPM at surface), while historical in situ values were usually less than 1 mg m–3. This suggests that resuspended bottom sediment due to wind-driven mixing and winter cooling is responsible for the noticeable overestimation of satellite-derived Chl-a concentrations. The algorithm for ocean color needs to be improved urgently for turbid water.  相似文献   

4.
Synoptic ship and satellite observations were performed of the Kuroshio warm-core ring (KWCR) 93A and its adjacent waters, off Sanriku, northwestern North Pacific, between early April and late June 1997. The temporal and spatial distribution of chlorophylla (Chl-a) and sea surface temperature in the study area were analyzed using data from ADEOS Ocean Color and Temperature Scanner (OCTS) and NOAA Advanced Very High Resolution Radiometer (AVHRR). The objective of this study was to describe the temporal and spatial variability of the spring bloom and understand its relationship with the changes in the hydrographic structure of these waters in and around KWCR 93A. The maximum value of Chl-a concentration in the ring was less than 1 mg/m3 during April. The spring bloom in the ring occurred early in May and the relatively high maximum (>1.0 mg/m3) continued from early in May to mid-June. In late June, a ship-observed surface Chl-a concentration of less than 0.4 mg/m3 suggests that the spring bloom had already declined in and around KWCR 93A. Double spiral structures of warm and cold streamers appeared from late April to mid-May, which may have an influence on the occurrence of the spring bloom in and around the ring. In this episodic event, the warm streamer can maintain the available potential energy of the ring and the strength of upwelling around the ring. The cold streamer provided water with a high Chl-a concentration to the surface layer of the ring. In order to understand the temporal and spatial variability of Chl-a distribution in the ring, the behavior of the warm and cold streamers needs to be taken into consideration.  相似文献   

5.
To trace offshore surface low-salinity water (LSW) in the northeastern Gulf of Mexico, a proxy was developed using the surface water beam attenuation coefficient (c p), and salinity matched with synchronous Sea-viewing Wide Field-of-view Sensor (SeaWiFS) satellite data from three annual summer cruises (July 1998–August 2000) using a two-step empirical approach. First, a relationship between in-situ salinity and c p was obtained. Second, in-situ c p was matched with SeaWiFS radiance ratios of all available blue-to-green wavelengths. Finally, satellite-derived surface salinity was determined directly by combining the two empirical relationships, providing a robust estimate over a range of salinities (26–36). This significantly improves the limited spatial and temporal resolution of surface salinity distribution obtained by shipboard sampling. The resulting correlation is best explained as mixing between low-salinity plume waters and normal salinity Gulf waters. The empirical relationships were used to map satellite-derived salinity using the average of SeaWiFS images during each summer cruise. As expected for summer, spatial patterns of LSW plumes with high c p, particulate matter (PM), particulate organic carbon (POC), and chlorophyll-a (Chl-a) were connected to the mouth of the Mississippi River Delta and extended to the east-southeast. Normal salinity Gulf water with lower c p, PM, POC, and Chl-a was confined to the shelf and upper slope in the eastern part of the study area, outside the plumes. This proxy approach can be applied throughout the region of shipboard sampling for more detailed coverage and analysis.  相似文献   

6.
An algorithm was evaluated for estimating primary production in the Kuroshio front of the East China Sea. The algorithm involves three principal inputs: chlorophyll-a (Chl-a) vertical profile, water column photosynthetic parameter and underwater irradiance parameter. We construct empirical models for estimating those three inputs from sea surface variables accessible by means of satellite remote sensing. Types of Chl-a vertical profile were closely associated with water mass features, and sea surface temperature seemed to be a potential variable to distinguish profile types. The consequence of the deep Chl-a maximum (DCM) presence was investigated by comparing the vertically non-uniform Chl-a profile against the uniform Chl-a profile in terms of the integrated biomass and primary production estimates. Since the DCM in the frontal region occurred in a relatively shallow layer with relatively high light level, allowing intensive photosynthetic processes to take place in the DCM layer. In contrast, as the DCM in the Kuroshio region occurred in the low-light, deep layer, it contributed less to the integrated primary production. We inferred that the DCM formation is more important in the frontal region than in the Kuroshio region, as they potentially contributed within 30.9% (±9.1%) and 20.9% (±5.4%) to the integrated primary production in the frontal and the Kuroshio regions, respectively. Such different degrees in the DCM contribution was then revealed by the higher integrated primary production in the frontal region than in the Kuroshio region.  相似文献   

7.
Ninety-four stations were sampled in the Atlantic subtropical gyres during 10 cruises carried out between 1995 and 2001, mainly in boreal spring and autumn. Chlorophyll a (Chl-a) and primary production were measured during all cruises, and phytoplankton biomass was estimated in part of them. Picoplankton (<2 μm) represented >60% of total Chl-a concentration measured at the surface, and their contribution to this variable increased with depth. Phytoplankton carbon concentrations were higher in the upper metres of the water column, whereas Chl-a showed a deep maximum (DCM). At each station, the water column was divided into the upper mixed layer (ML) and the DCM layer (DCML). The boundary between the two layers was calculated as the depth where Chl-a concentration was 50% of the maximum Chl-a concentration. On average DCML extends from 67 to 126 m depth. Carbon to Chl-a (C:Chl-a) ratios were used to estimate phytoplankton carbon content from Chl-a in order to obtain a large phytoplankton carbon dataset. Total C:Chl-a ratios averaged (±s.e.) 103±7 (n=22) in the ML and 24±4 (n=12) in the DCML and were higher in larger cells than in picoplankton. Using these ratios and primary production measurements, we derived mean specific growth rates of 0.17±0.01 d−1 (n=173) in the ML and 0.20±0.01 d−1 (n=165) in the DCML although the differences were not significant (t-test, p>0.05). Our results suggest a moderate contribution of the DCML (43%) to both phytoplankton biomass and primary production in the Atlantic subtropical gyres.  相似文献   

8.
A bio-optical dataset collected during the 1998?C2007 period in the Yellow and East China Seas (YECS) was used to provide alternative empirical ocean-color algorithms in the retrieval of chlorophyll-a (Chl-a), total suspended matter (TSM), and colored dissolved organic matter (CDOM) absorption coefficients at 440 nm (ag440). Assuming that remote-sensing reflectance (Rrs) could be retrieved accurately, empirical algorithms for TChl (regionally tuned Tassan??s Chl-a algorithm) in case-1 waters (TChl2i in case-2 waters), TTSM (regionally tuned Tassan??s TSM algorithm), and Tag440 or Cag440 (regionally tuned Tassan??s or Carder??s ag440 algorithm) were able to retrieve Chl-a, TSM, and ag440 with uncertainties as high as 35, 46, and 35%, respectively. Applying the standard SeaWiFS Rrs, TChl was not viable in the eastern part of the YECS, which was associated with an inaccurate SeaWiFS Rrs retrieval because of improper atmospheric correction. TChl behaved better than other algorithms in the turbid case-2 waters, although overestimation was still observed. To retrieve more reliable Chl-a estimates with standard SeaWiFS Rrs in turbid water (a proxy for case-2 waters), we modified TChl for data with SeaWiFS normalized water-leaving radiance at 555 nm (nLw555) > 2 mW cm?2 ??m?1 sr?1 (TChl2s). Finally, with standard SeaWiFS Rrs, we recommend switching algorithms from TChl2s (for case-2 waters) to MOCChl (SeaWiFS-modified NASA OC4v4 standard algorithm for case-1 waters) for retrieving Chl-a, which resulted in uncertainties as high as 49%. To retrieve TSM and ag440 using SeaWiFS Rrs, we recommend empirical algorithms for TTSM (pre-SeaWiFS-modified form) and MTag440 or MCag440 (SeaWiFS Rrs-modified forms of Tag440 or Cag440). These could retrieve with uncertainties as high as 82 and 52%, respectively.  相似文献   

9.
River plumes have important effects on marine ecosystems. Variation in the extent and dispersal of river plumes is often associated with river discharge, wind characteristics and ocean circulation. The objectives of this study were to identify the Tokachi River plume by satellite, determine its relationship with river discharge and clarify its temporal and spatial dynamics. SeaWiFS multispectral satellite data (normalized water-leaving radiance: nLw) with 1.1 km spatial resolution were used to determine the spatial and temporal variability of the plume during 1998–2002. Supervised maximum likelihood classification using six channels of nLw at 412, 443, 490, 510, 555 and 670 nm with each band's spectral signature statistic was used to define classes of surface water and to estimate the plume area. Supervised maximum likelihood classification separated three to four classes of coastal water based on optical characteristics as a result of wind stress events. The satellite-observed plume area was correlated with the amount of river discharge from April to October. The plume distribution patterns were influenced by wind direction and magnitude, the occurrences of a near-shore eddy field and surface currents. Empirical orthogonal function (EOF) was used to express the spatial and temporal variability of the plume using anomalies of nLw(555) monthly averaged images. The first mode (44% of variance) showed the turbid plume distribution resulting from re-suspension by strong wind mixing along the coast during winter. This mode also showed the plume was distributed along-shelf direction in spring to early autumn. The second mode (17% of variance) showed spring pattern across-shelf direction. EOF analysis also explained the interannual variability of the plume signature, which might have been affected by the flow of the Oyashio Current and the occurrence of a near-shore eddy field.  相似文献   

10.
根据2011年6月茅尾海生态环境调查资料,对该海域海水和表层沉积物中叶绿素a的空间分布进行了分析。结果表明,海水叶绿素a变化范围1.384~4.060 mg/m3,平均值为2.143 mg/m3,表层沉积物叶绿素a范围为0.006~0.740 mg/kg(湿重),均值为0.124 mg/kg;海水与表层沉积物叶绿素a均呈现自河口向南逐渐降低的空间分布特征。单位面积表层沉积物叶绿素a平均含量为上方水柱叶绿素a平均含量的129.44%,沉积物对该海域初级生产力有显著的潜在贡献。相关分析表明,海水和沉积物叶绿素a均与无机氮及底栖动物栖息密度呈显著或极显著的正相关关系(P0.05或P0.01)。  相似文献   

11.
Both historic and currently operational chlorophyll algorithms of the satellite-borne ocean color sensors, such as SeaWiFS, were evaluated for in situ spectral radiation and chlorophyll data in some Case I waters, including the waters in the Indian Ocean sector of the Southern Ocean. Chlorophyll a concentration of the data set (n = 73) ranged from 0.04 to 1.01 mg m–3. The algorithms had higher accuracy for the low- and mid-latitude waters (RMSE: 0.163–0.253), specifically the most recently developed algorithms of OCTS and Sea WiFS showed 0.163 and 0.170 of Root Mean Square Errors, respectively. However, these algorithms had large errors (0.422–0.621) for the Southern Ocean data set and underestimated the surface chlorophyll by more than a factor of 2.6. The absorption coefficients in the blue spectral region retrieved from remote sensing reflectance varied in a nonlinear manner with chlorophyll a concentration, and the value in the Southern Ocean was significantly lower than that in the low- and mid-latitude waters for each chlorophyll a concentration. The underestimation of chlorophyll a concentration in the Southern Ocean with these algorithms was caused by the lower specific absorption coefficient in the region compared with the low- and mid-latitude waters under the same chlorophyll a concentration.  相似文献   

12.
Seasonal and interannual variations in physicochemical properties (i.e., temperature, salinity, dissolved oxygen and dissolved inorganic nutrients), chlorophyll a (Chl-a), particulate carbon and nitrogen (PC and PN, respectively), and primary production were investigated in the neritic area of Sagami Bay, Kanagawa, Japan, from January 2002 to December 2008. These abiotic/biotic variables, except for NH4 +–N, repeated similar seasonal variations for all 7 years. On the basis of the analysis of data obtained on 167 sampling dates, depth-integrated primary production in this water can be easily estimated from Chl-a at the surface using the regression equations obtained in the present study. Intermittently high values of dissolved inorganic nutrients, Chl-a, PC, PN and primary productivity at the surface during the summer stratified period were induced by high freshwater discharge from the rivers after rainfalls and by the expansion of nutrient-rich Tokyo Bay Water. Temperature, salinity and dissolved inorganic nutrients showed drastic variations within a scale of a few days and/or weeks, and these variations were related to sea levels that represent the intrusion of the Kuroshio Water, Intermediate Oyashio Water or deep water from the continental slope. However, there was no consistent trend in the variations in Chl-a, PC, PN and primary production due to the complex effects of these waters.  相似文献   

13.
The variability of Chlorophyll-a (Chl-a) distribution derived from MODIS (on Aqua and Terra platforms) and MERIS sensors have been compared with SeaWiFS data in the Arabian Sea. MODIS Aqua has overestimated the SeaWiFS Chl-a within 25–32% in the coastal turbid (eutrophic) waters and underestimated in open ocean waters with error within 20%. However, there is no significant bias (?0.1 on log-scale) observed as the slope is well within 0.97-1.1 (log transformed). MODIS-Terra has underestimated the Chl-a concentration in open ocean waters by about 29–31%, which is higher than MODIS-Aqua. MODIS-Terra is observed to be more accurate than MODIS-Aqua in the coastal waters. MERIS is overestimating the SeaWiFS Chl-a with log RMS error of ~0.15 and log bias of ~0.13–0.2. The differences in the Chl-a estimates between each sensor are possibly due to differences in the sensor design, bio-optical algorithms and also due to the time differences between the satellites over passes. We have examined that the MERIS is performing similar to SeaWiFS and the MODIS-Aqua (Terra) data are reliable in open ocean (coastal) waters. However, Chl-a retrieval algorithms need to be improved especially for coastal turbid waters to continue with SeaWiFS data for long-term studies.  相似文献   

14.
The aim of this study was to identify the factors responsible for the differences in chlorophyll a concentration (Chl-a) observed between the California, Canary, Humboldt and Benguela upwelling areas. Monthly climatologic values of Chl-a obtained from satellite images, covering the years 1998–2004, revealed that this pigment was higher in the Benguela system than in the other areas. Upwelling intensity, as derived from offshore Ekman transport computations, was higher in the Benguela and Humboldt regions and, for the same upwelling intensity, Chl-a was higher in Benguela than in the other regions. Upwelling intensity appears to be able to drive Chl-a densities through nutrient supply, as nutrients are correlated to offshore Ekman transport. A linear regression model including the fraction of sea surface over the shelf in each 1° × 1° box, nitrate, silicate, turbulence and variability of offshore Ekman transport explained the 84.8% of the variance in Chl-a among the areas. Differences in offshore Ekman transport explained the lower Chl-a observed in Canary and California and the higher Chl-a observed in Benguela and Peru-Humboldt. A narrow continental shelf and low water column stability also contribute to reducing phytoplankton pigment biomass in the Canary and California areas. The higher Chl-a values observed in Benguela compared to Humboldt-Peru are due to a wider extension of the continental shelf in the Benguela region.  相似文献   

15.
As part of the environmental recovery plan for the Mar Piccolo (Ionian Sea), a basin suffering especially from urban pollution, a study was planned on the modifications induced by the human activities to distinguish them from those due to the alternation of the seasons and the exchanges with the neighbouring basin of Mar Grande. The present paper deals with the chemical-physical characteristics of sea water and its particulate matter content is investigated as a part of this study. The physical, chemical and biological variables (stochastic variables for the model) are: chlorophyll a, proteins, carbohydrates, organic matter, dry weight, temperature, salinity, dissolved oxygen, pH, ultraviolet absorbance, phosphate, nitrite, nitrate, ammonia, silicate, photosynthetically available radiation (PAR) and energy. These variables have been sampled at 15 different stations during one year of monthly sampling. A preliminary analysis was carried out on the variables for detecting possible interactions among all the investigated parameters. To this purpose, and for all the variables, an original computer program for non-parametric statistical-evaluation and data-correlation was used to find correlation-guided hierarchy of variables. The relations among: Prot and Chl-a, Carb and Chl-a and OM and DW have been studied. The slopes of linear regression of (Prot, Chl-a) and (Carb, Chl-a) allowed us to detect different relations at different space-time areas. A nonlinear relation between OM and DW has been found as well. This treatment is a first approach to the construction of a mathematical model of the basin. Different areas, for sea water characteristics, were observed to occur in the basin, despite its small extension. This may be due to the reduced water exchanges and the different residential development along the coast. Innermost areas, less involved in the water mass circulation, were more affected by eutrophication.  相似文献   

16.
We describe the oceanographic condition as observed by hydrographic data and phytoplankton spring bloom detected by OCTS images off Sanriku, northwestern Pacific, during the spring bloom period in 1997. The relationship between the two is discussed. OCTS images detected the bloom in early April in the coastal area around the Izu ridge north of the Kuroshio and the eastern coastal area of Hokkaido to the Oyashio front. The bloom areas were seen along the offshore Kuroshio Extension from the end of April, in the upstream region of the Oyashio south of the Kurile Islands, except for a part of coastal area from the end of May, and in the Kuroshio warm-core ring 93A (KWCR 93A) from early June. The temperature difference between the surface and subsurface layer is used as a stratification index. This was large in the upstream region of the Oyashio south of the Kurile Islands and KWCR 93A from early June. Previous research has pointed out that the spring bloom usually corresponds to the development of stratification in the water column due to seasonal warming. In addition to that, we suggest that the transportation of water containing a high chlorophylla concentration by advection due to strong currents, like the Kuroshio and the Oyashio, is important for the formation of an area of high chlorophylla concentration. These results indicate that the OCTS images are useful for a knowledge of the distribution and the change of chlorophylla concentration in the northwestern Pacific region.  相似文献   

17.
Determinations of Chlorophyll a (Chl a) in eutrophic coastal marine waters were compared using N, N-dimethylformamide (DMF) and 90% acetone techniques. Measured Chl a concentrations ranged from 0.89 to 10.65 μg l−1 for 90% acetone extracts and from 0.97 to 12.92 μg l−1 for DMF extracts, respectively, for 24 surface water samples from the Harima-nada, Seto Inland Sea, Japan. Chl a concentrations using DMF as a solvent were consistently higher than those found using 90% acetone (p < 0.001, n= 24). Chl a is thus likely to be underestimated (by 13%) when 90% acetone is used for eutrophic waters. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

18.
The Ocean Color and Temperature Scanner (OCTS) on board the Advanced Earth Observing Satellite (ADEOS) observed the chlorophyll-a concentration in the surface layer and sea surface temperature in global oceans from October 1996 to June 1997. The OCTS team was formed in the National Space Development Agency of Japan (NASDA) Earth Observation Research Center (EORC) to develop algorithm, calibrate and validate OCTS products and promote OCTS data usage. Intensive efforts to improve the quality of the OCTS products were made after the launch of ADEOS. Much sea-truth data was collected, and the algorithms to retrieve the ocean parameters have been revised several times. The OCTS data were distributed to the user community through the Internet while OCTS was functioning. An overview of the OCTS mission is presented in this paper. The OCTS Team is formed by NASDA personnel and supporting scientists listed in Appendix 1.  相似文献   

19.
Simultaneous acquisition of water samples, radiance and irradiance measurements were carried out from 40 stations in the Mandovi–Zuari estuaries during February to May 2002. From the samples collected, inherent and apparent optical properties (IOP and AOP) such as absorption coefficient (a), upwelling diffuse attenuation coefficient (ku) and subsurface reflectance (R) were derived. Using these optical properties, radiative transfer at each water column is examined. On the basis of the radiative transfer outcome, band-ratio algorithms are derived for three optically active substances (OAS), viz, chlorophyll-a, suspended sediment and coloured dissolved organic matter (CDOM). The respective algorithms are 670/555, 490/670 and 412/670 nm for chlorophyll-a, suspended sediment and CDOM. These algorithms are applied to Ocean Colour Monitor (OCM), onboard Indian Remote Sensing Satellite (IRS)-Polar Satellite Launch Vehicle (P4), scenes (digital data), to synoptically analyze these OAS. The synoptic analysis of OAS revealed different hydrodynamic characteristics of the estuaries during non-monsoon seasons.  相似文献   

20.
The distribution of transparent exopolymer particles (TEP) was investigated during a coccolithophorid bloom in the northern Bay of Biscay (North Atlantic Ocean) in early June 2006. MODIS chlorophyll-a (Chl-a) and reflectance images before and during the cruise were used to localize areas of important biological activity and high reflectance (HR). TEP profiles along the continental margin, determined using microscopic (TEPmicro) and colorimetric (TEPcolor) methods, showed abundant (6.1×106–4.4×107 L?1) and relatively small (0.5–20 μm) particles, leading to a low total volume fraction (0.05–2.2 ppm) of TEPmicro and similar vertical profiles of TEPcolor. Estimates of carbon content in TEP (TEP-C) derived from the microscopic approach yielded surface concentration of 1.50 μmol C L?1. The contribution of TEP-C to particulate organic carbon (POC) was estimated to be 12% (molar C ratio) during this survey. Our results suggest that TEP formation is a probable first step to rapid and efficient export of C during declining coccolithophorid blooms.  相似文献   

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

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