首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
The long-term mean (31-year mean) surface heat fluxes over the Japan Sea are estimated by the bulk method using the most of the available vessel data with the resolution of 1o×1o. The long-term annual mean net heat flux is about –53 W m–2 (negative sign means upward heat flux) with the annual range from 133 W m–2 in May to –296 W m–2 in December. The small gain of heat in the area near Vladivostok seems to indicate the existence of cold water flowing from the north. In that area in winter, the mean loss of heat attains about 200 W m–2, and the Bowen's ratio is over the unity. The largest insolation occurs in May in the Japan Sea, and the upward latent heat flux becomes the largest in November in this area. The heat flux of Haney type is also calculated, and the result, shows that the constantQ 1 has the remarkable seasonal and spatial variation, while the coefficientQ 2 has relatively small variation throughout all seasons. Under the assumption of constant volume transport of 1.35×106 m3s–1 through the Tsugaru Strait, the long-term averages of the volume transport through the Tsushima and Soya Straits are estimated to be about 2.20 and 0.85×106 m3s–1 from the result of the mean surface heat flux, respectively.  相似文献   

2.
Dissolved gaseous mercury (DGM) was measured continuously using two newly developed techniques and a manual technique. The continuous techniques were based on the equilibrium between the aqueous and gaseous phase (DGM = Hgextr / H', Hgextr is the measured mercury concentration in the gas phase, H' is the Henry's Law coefficient at the desired temperature). In order to calculate the annual mercury evasion from the Mediterranean Sea, diurnal and seasonal measurements of DGM, total gaseous mercury in air (TGM), water temperature and wind speed were performed. During August 2003, March–April 2004 and October–November 2004 measurements of these parameters were conducted on board the RV Urania. The continuous measurements of DGM showed a diurnal variation in concentration, at both coastal and off shore sites, with higher concentrations during daytime than nighttime. The concentration difference could be as large as 130 fM between day and night. The degree of saturation was calculated directly from the measurements, S = Hgextr / TGM and was found to vary between the different seasons. The highest average degree of saturation (850%) and the largest variation in saturation (600–1150%) was observed during the summer. The spring showed the lowest variation (260–360%) and the lowest average degree of saturation (320%). The autumn also showed a large variation in saturation (500–1070%) but a lower average (740%) compared to the summer cruise. This might be explained by the temperature difference between the different seasons, since that parameter varied the most. The flux from the sea surface was calculated using the gas exchange model developed by Nightingale et al. [Nightingale, P.D., Malin, G., Law, C.S., Watson, A.J., Liss, P.S., Liddicoat, M.I., Boutin, J., Upstill-Goddard, R. C., 2000. In situ evaluation of air–sea gas exchange parameterization using novel conservative and volatile tracers. Global Biogeochemical Cycles, 14(1):373–387]. The evasion varied between the different seasons with the highest evasion during the autumn, 24.6 pmol m− 2 h− 1. The summer value was estimated to 22.3 pmol m− 2 h− 1 and the spring to 7.6 pmol m− 2 h− 1. Using this data the yearly evasion from the Mediterranean Sea surface was estimated to 77 tons.  相似文献   

3.
In January–February 1987, an urgent cruise JENEX-87 was carried out in the central equatorial Pacific during the occurrence of the 1986–87 El Niño. This cruise, supported by the Japan Science and Technology Agency, supplied heat flux data through the sea surface, on the basis of direct measurements of short- and long-wave radiation fluxes.In the time average, the heat gain due to the radiation flux (153 W m–2) was almost compensated by the heat loss due to latent heat flux (130 W m–2), and thus the net heat gain was small in magnitude (20 W m–2). On the other hand, day-to-day changes of the net heat flux ranged within ±130 W m–2, mainly reflecting the downward short-wave radiation variations.The heat balance in the surface oceanic mixed layer was investigated in two quadrangle areas (160°E-180° and 180°-160°W between 2°N and 2°S), using the surface heat flux and estimating the advective heat fluxes due to the geostrophic and Ekman currents. In these two quadrangles, we respectively derived –187±88 W m–2 and +27±95 W m–2. The former value, which is equivalent to about 1°C month–1 drop of the mixed layer temperature, is evidence of the abnormal oceanic condition in the occurrence of the 1986–87 El Niño event.  相似文献   

4.
A new algorithm using a multivariate regression technique for retrieving sea surface specific humidity (Q) from remote sensing data from the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) is proposed. Daily and monthly specific humidity data from the National Center for Environmental Prediction (NCEP) reanalysis dataset and data of sea surface temperature, atmospheric total water vapor, and wind speed from AMSR-E oceanographic products were used to derive the regression coefficients of the algorithm, and all the data for derivation are from the year 2003. An F-test was applied to the regression, and small P-values indicate that the regressions are significant to a high level of confidence. The derived coefficients have been validated using similar data from the year 2004. The root mean square (rms) error of the algorithm for daily retrieved Q over the global oceans is 1.05 g kg−1, and the rms error for monthly retrieved Q is 0.61 g kg−1.  相似文献   

5.
The purpose of this study is to validate and improve satellite-derived downward surface shortwave radiation (DSSR) over the northwestern Pacific Ocean using abundant in situ data. The DSSR derivation model used here assumes that the reduction of solar radiation by clouds is proportional to the product of satellite-measured albedo and a cloud attenuation coefficient. DSSR is calculated from Geostationary Meteorological Satellite-5/Visible Infrared Spin-Scan Radiometer data in 0.05° × 0.05° grids. The authors first compare the satellite DSSR derived with a cloud attenuation coefficient table determined in past research with in situ values. Although the hourly satellite DSSR agrees well with land in situ values in Japan, it has a bias of +13∼+34 W/m2 over the ocean and the bias is especially large in the low latitudes. The authors then improve the coefficient table using the ocean in situ data. Usage of the new table successfully reduces the bias of the satellite DSSR over the ocean. The cloud attenuation coefficient for low-albedo cases over the ocean needs to be larger in the low latitudes than past research has indicated. Daily and hourly DSSR can be evaluated from the satellite data with RMS errors of 11–14% and 30–33%, respectively, over a wide region of the ocean by this model. It is also shown that the cloud attenuation coefficient over land needs to be smaller than over the ocean because the effect of the radiation reflected by the land surface cannot be ignored.  相似文献   

6.
The nitrogen inputs from atmospheric deposition and bottom water entrainment to the surface layer were modelled in the summer period (May–September) over a 11-year period (1989–1999) and compared to investigate the significance of these fluxes for generating blooms in the Kattegat. In the summer periods the average atmospheric deposition was 2.81 mg N m−2 d−1 compared to average entrainment fluxes of 5.42 mg N m−2 d−1, 1.21 mg N m−2 d−1 and 1.15 mg N m−2 d−1 for the northern, central and southern part of the Kattegat, respectively. Atmospheric nitrogen deposition alone could not sustain biomass increases associated with observed blooms and entrainment fluxes dominated the high nitrogen inputs to the surface layer. The potential for a bloom through growth was typically obtained after several days of high nitrogen inputs from entrainment in the frontal area of the northern Kattegat and to some extent from atmospheric deposition. The modelled nitrogen input in this area could account directly for 30% of the observed blooms in the Northern sub-basin, and through advective transport 24% and 19% of the observed blooms in the central and southern Kattegat. The direct nitrogen inputs through atmospheric deposition and entrainment to the central and southern sub-basins were small and could not be linked to any bloom observation.  相似文献   

7.
Surface solar radiation over the Pacific Ocean off the Sanriku coast has been estimated using Visible and Infrared Spin Scan Radiometer data supplied by the Geostationary Meteorological Satellite 5 for September, 1996 to June, 1997, when the Ocean Color and Temperature Scanner was functioning. The hourly and daily insolation is estimated with a spatial resolution of 0.01-degree grid. Thein situ surface short wave radiation obtained by the research vessel,Kofu-Maru belonging to the Japan Meteorological Agency is used for validation of the estimated insolation. It is shown that the estimated hourly and daily insolation has an rms (root mean square) error of 17.05% and 8.13%, respectively, which are the ratios between the rms error (W/m2) and the mean insolation (W/m2).  相似文献   

8.
A method of obtaining the operative sea surface temperature (SST)t using satellite scanner observations in the spectral ranges 3.53–3.94 m and 10.3–11.3 m is realized. The method represents a combination of McClainet al.'s formula (1983) and expressions suggested by the authors which describe the universal angular structures of the radiation temperature fields. The RMS error of reconstructingt at scanning angles of 0–55o is equal to 0.2–0.3°C for atmospheric states corresponding to the SST variation within the limit 6–28°C. An atlas of temperature maps on the grid 0.5×0.5o with temporal averaging from 5 days to 1 month is compiled using the data obtained on board the RVAkademik Vernadsky in the Atlantic Ocean in 1987–1989.Translated by Mikhail M. Trufanov.  相似文献   

9.
New satellite-derived latent and sensible heat fluxes are performed by using Wind Sat wind speed, Wind Sat sea surface temperature, the European Centre for Medium-range Weather Forecasting(ECMWF) air humidity, and ECMWF air temperature from 2004 to 2014. The 55 moored buoys are used to validate them by using the 30 min and 25 km collocation window. Furthermore, the objectively analyzed air-sea heat fluxes(OAFlux) products and the National Centers for Environmental Prediction-National Center for Atmospheric Research reanalysis 2(NCEP2) products are also used for global comparisons. The mean biases of sensible and latent heat fluxes between Wind Sat flux results and buoy flux data are –0.39 and –8.09 W/m~2, respectively. In addition, the rootmean-square(RMS) errors of the sensible and latent heat fluxes between them are 5.53 and 24.69 W/m~2,respectively. The RMS errors of sensible and latent heat fluxes are observed to gradually increase with an increasing buoy wind speed. The difference shows different characteristics with an increasing sea surface temperature, air humidity, and air temperature. The zonal average latent fluxes have some high regions which are mainly located in the trade wind zones where strong winds carry dry air in January, and the maximum value centers are found in the eastern waters of Japan and on the US east coast. Overall, the seasonal variability is pronounced in the Indian Ocean, the Pacific Ocean, and the Atlantic Ocean. The three sensible and latent heat fluxes have similar latitudinal dependencies; however, some differences are found in some local regions.  相似文献   

10.
We have observed the temporal variation of oxygen deficient water with short time scale (less than a few days) in the central area of Ohmura Bay, Kyushu, Japan, in summer, 1995 and 1996. The vertical profiles of temperature were similar to those of dissolved oxygen. We noticed a linear relation between temperature and dissolved oxygen in the bottom layer, and applied the T-DO relation to estimate the net oxygen consumption rate, rather than direct evaluation of the advection and diffusion. Oxygen consumption rate just above the bottom was estimated to be about 0.21 g O2 m–3day–1 in July 1995, and about 0.28 g O2 m–3day–1 in August 1996. The net oxygen consumption rate estimated for the bottom layer below the second thermocline was about 0.61 g O2 m–3day–1 with variability from 0.55 to 0.66 g O2 m–3day–1 during July 25 to 29, 1995. This is was about 0.64 g O2 m–3day–1 with variability from 0.18 to 1.4 g O2 m–3day–1 during August 22 to 30, 1996. The net oxygen consumption rates are about half of those measured with a closed system in the Seto Inland Sea.  相似文献   

11.
A three-dimensional hydrodynamic model has been developed to simulate water mass circulation in estuarine systems. This model is based on the primitive equation in Cartesian coordinates with a terrain-following structure, coupled with a Mellor–Yamada 2.5 turbulence scheme. A fractional-step method is applied and the subset of equations is solved with finite volume and finite element methods. A dry–wet process simulates the presence of the tidal flat at low water. River inputs are introduced using a point-source method. The model was applied to a partially mixed, macrotidal, temperate estuary: Southampton Water, UK. The model is validated by comparisons with sea surface elevation, ADCP measurements and salinity data collected in 2001. The mean spring range 2(M2 + S2) and the mean neap range 2(M2 − S2) are modelled with an error relative to observation of 12 and 16%, respectively. The unique tidal regime of the system with the presence of the ‘young flood stand’ corresponding to the slackening conditions occurring at mid flood and ‘double high water’ corresponding to an extension of the slackening conditions at high tide is accurately reproduced in the model. The dynamics of the modelled mean surface and bottom velocity closely match the ADCP measurements during neap tides (rms of the difference is 0.09 and 0.01 m s−1 at the bottom and at the surface, respectively), whereas at spring the difference is greater (rms of the difference is 0.25 and 0.20 m s−1 at bottom and surface, respectively). The spatial and temporal variation of the degree of stratification as indicated by salinity distributions compares well with observations.  相似文献   

12.
Concentrations of total carbonate, alkalinity and dissolved oxygen were obtained near the 1973 GEOSECS stations in the North Pacific subpolar region north of 40°N along 175°E between 1993 and 1994. A difference of excess CO2 content between the GEOSECS and our expeditions was estimated. The maximum difference in water column inventory of excess CO2 has increased by about 280 gC m–2 above 2000 m depth which apparently means an uptake of excess CO2 taken from air to sea during the last two decades. An averaged value of the annual flux of excess CO2 at 75–1000 m depth was 8.63±2.01 gC m–2yr–1 in the North Pacific subpolar region. By introducing the annual flux of excess CO2 into a two-box model for the North Pacific subpolar region, a penetration factor of excess CO2 from air to sea was obtained to be 1.08×10–2 gC m–3ppm–1 in the North Pacific subpolar region. Based on this factor, the surface concentration of excess CO2 in the North Pacific subpolar region was estimated to be 68 mole I–1, suggesting that the North Pacific subpolar region absorbed atmospheric excess CO2 more than the saturated concentration of excess CO2. Total amount of excess CO2 taken from the North Pacific subpolar region by 1993 was estimated to be 36.2×1015 gC, which was equal to about one tenth of that released by human activities after the preindustrial era.  相似文献   

13.
The shortwave radiative forcings of smoke aerosol in the cloudless atmosphere during the summer fires of 2010 in European Russia were quantitatively estimated for the land surface and the atmospheric upper boundary from measurement data obtained at the Zvenigorod Scientific Station of the Obukhov Institute of Atmospheric Physics (OIAP ZSS), Russian Academy of Sciences. Variations in the temperature of the surface air layer due to the smoke-induced attenuation of incoming solar radiation were estimated. The most intensive smoke generation in the atmosphere was observed on August 7–9, 2010, when the maximum aerosol optical thickness amounted to more than 4.0 at a wavelength of 550 nm. In this case, the albedo of single aerosol scattering amounted to ∼0.95–0.96 and the asymmetry factor amounted to ∼0.69–0.70. The maximum shortwave radiative forcing of aerosol amounted to about −360 W/m2 for the land surface and almost −150 W/m2 for the atmospheric upper boundary. During the period of intensive smoke generation, the cooling of the atmospheric surface layer over daylight hours (12 h) amounted, on average, to ∼6°C. The power character of the dependence of the shortwave radiative forcing of aerosol for the land surface on aerosol optical thickness up to its values exceeding 4.0, which was revealed earlier on the basis of data on aerosol optical thickness (up to 1.5) obtained at the OIAP ZSS during the summer forest and peatbog fires of 2002 in the region of Moscow, was supported.  相似文献   

14.
We selected surface flux datasets to investigate the heat fluxes during “hot events”; (HEs), defined as short-term, large-scale phenomena involving very high sea surface temperature (SST). Validation of the heat fluxes against in-situ ones, which are estimated from in-situ observation in HE sampling conditions, shows the accuracies (bias ± RMS error) of net shortwave radiation, net long wave radiation, latent heat and sensible heat fluxes are 20 ± 45.0 W m−2, −9 ± 12.3 W m−2, −2.3 ± 31.5 W m−2 and 1.5 ± 5.0 W m−2, respectively. Statistical analyses of HEs show that, during these events, net solar radiation remains high and then decreases from 246 to 220 W m−2, while latent heat is low and then increases from 100 W m−2 to 124 W m−2. Histogram peaks indicate net solar radiation of 270 W m−2 and latent heat flux of 90 W m−2 during HEs. Further, HEs are shown to evolve in three phases: formation, mature, and ending phases. Mean heat gain (HG) in the HE formation phase of 60 W m−2 is larger than the reasonably estimated annual mean HG range of 0–25 W m−2 in the Indo-Pacific Warm Pool. Such large daily HG in the HE formation phase can be expected to increase SSTs and produce large amplitudes of diurnal SST variations during HEs, which have been observed by both satellite and in-situ measurements in our previous studies.  相似文献   

15.
This study focused on the causes of the variation in microphytobenthic biomass and the effects of this variation on macrobenthic animals in the western Seto Inland Sea, Japan, where the importance of microphytobenthos as the primary food source for benthic animals has been recently reported. We investigated the microphytobenthic biomass together with light attenuation of seawater, phytoplanktonic biomass, macrobenthic density and biomass at eight stations (water depth = 5–15 m) during four cruises in 1999–2000. The increased light attenuation coefficient of the water column associated with increased concentration of the phytoplanktonic Chl-a caused a decrease in light flux that reached the seafloor. The biomass of the microphytobenthos within the upper 1 cm of the sediment, 1.9–46.5 mg Chl-a m−2, was inversely correlated with the phytoplanktonic biomass in the overlying water column, 10.9–65.0 mg Chl-a m−2. Thus, interception of light by phytoplankton is considered to be a main cause of the variation in the microphytobenthic biomass. The microphytobenthos biomass showed a significant positive correlation with the macrobenthic density (78–9369 ind. m−2) and biomass (0.4–78.8 gWW m−2). It appears that the increase in oxygen production by the microphytobenthos allowed macrobenthic animals to become more abundant, as a consequence of oxygenation of the organically enriched muddy sediments (14.5 ± 2.69 mg TOC g−1). This study suggests that the variation in the microphytobenthic biomass is influenced by the phytoplanktonic biomass due to shading effect, and the balance between these two functional groups might affect the variability in the macrobenthic density and biomass.  相似文献   

16.
Monthly seawater pH and alkalinity measurements were collected between January 1996 and December 2000 at 10°30′N, 64°40′W as part of the CARIACO (CArbon Retention In A Colored Ocean) oceanographic time series. One key objective of CARIACO is to study temporal variability in Total CO2 (TCO2) concentrations and CO2 fugacity (fCO2) at this tropical coastal wind-driven upwelling site. Between 1996 and 2000, the difference between atmospheric and surface ocean CO2 concentrations ranged from about − 64.3 to + 62.3 μatm. Physical and biochemical factors, specifically upwelling, temperature, primary production, and TCO2 concentrations interacted to control temporal variations in fCO2. Air–sea CO2 fluxes were typically depressed (0 to + 10 mmol C m 2 day 1) in the first few months of the year during upwelling. Fluxes were higher during June–November (+ 10 to 20 mmol C m 2 day 1). Fluxes were generally independent of the slight changes in salinity normally seen at the station, but low positive flux values were seen in the second half of 1999 during a period of anomalously heavy rains and land-derived runoff. During the 5 years of monthly data examined, only two episodes of negative air–sea CO2 flux were observed. These occurred during short but intense upwelling events in March 1997 (−10 mmol C m 2 day 1) and March 1998 (− 50 mmol C m 2 day 1). Therefore, the Cariaco Basin generally acted as a source of CO2 to the atmosphere in spite of primary productivity in excess of between 300 and 600 g C m 2 year 1.  相似文献   

17.
Measurement of dissolved oxygen using optodes in a FerryBox system   总被引:2,自引:0,他引:2  
Optode sensors can provide detailed information on concentrations of dissolved oxygen, which in turn may be used to quantify variations in net primary productivity. Throughout 2005 and 2006 the performance of commercially available oxygen optodes was examined, one in each year. The optode was part of an autonomous measurement system (FerryBox) on a ferry operating between Portsmouth (UK) and Bilbao (Spain). On crossings during which water samples were collected manually, the optode outputs were compared to measurements of dissolved oxygen made by Winkler titrations. The optodes maintained good stability with no evidence of instrumental drift during the course of a year. Over the observed concentration range (230–330 mM m−3) the optode data were approximately 2% low in both years. By fitting the optode data to the Winkler data the median difference between the optode and Winkler measurements is reduced to less than 1 mM m−3 (0.3%) in both years. The most appropriate calibration factor for 2005 was corrected O2 = Optode O2 × 1.018 and for 2006 the corresponding equation is corrected O2 = Optode O2 × 0.884 + 36.8. The standard deviation (95%) of the difference between the individual Winkler measurements was 5 mM m−3 and 3 mM m−3 in 2005 and 2006 respectively.Calculation of the oxygen saturation anomaly is required for calculation of the air sea exchange of oxygen and net biological production. This calculation requires the use of both salinity and temperature data. Salinity is measured to better than 0.1 so the corresponding error in anomaly is less than 0.2 mM m−3. Distortion of the temperature data is present due to warming of the water pumped to the optode. In winter this warming at the optode may be as great as 0.4 °C. The difference in the pumped water temperature can be corrected for by reference to other measurements of sea surface temperature reducing the error to less than 1 mM m−3.  相似文献   

18.
Direct observation of radiative flux in the southern yellow sea   总被引:2,自引:0,他引:2  
Direct measurements of four radiative components at air-sea boundary layer were conducted in the southern Yellow Sea during three cruises (seasons) in 2007. Simultaneous observations of meteorological (cloud cover, air temperature and humidity) and oceanographic (sea surface temperature) parameters were carried out. Observational results of radiative fluxes and meteorological and oceanographic parameters are presented. Mean diurnal cycles of four radiative components, net radiation, and sea surface albedo are calculated to achieve averages in different seasons. Net radiative fluxes in three seasons (winter, spring, autumn) are 8, 146, 60 W/m2, respectively. Comparisons between the observed radiative fluxes and those estimated with formulas are taken.  相似文献   

19.
Data are presented for the concentrations of Al, Fe, Mn, Ni, Co, Cr, V, Cu, Zn, Pb and Cd in aerosols collected over two contrasting regions of the Indian Ocean. These are: (1) the northern Arabian Sea (AS), from which samples were collected in the northeast monsoon, during which the region receives an input of crustal material from the surrounding arid land masses; and (2) the Tropical Southern Indian Ocean (TSIO), a remote region from which samples were collected from air masses for which there were no large-scale up-wind continental aerosol sources. The TSIO samples can be divided into two populations: Population I aerosols, collected from air masses which have probably impinged on Madagascar, and Population II aerosols, which have been confined to open-ocean regions to the south of the area.The data indicate that there are strong latitudinal variations in the chemical signatures of aerosols over the Indian Ocean. The input of crustal material to the Arabian Sea gives rise to an average Al concentration of about 1000 ng m−3 of air in the northeast monsoon regime. As a result, the concentrations of all trace metals are relatively high, and the values of crustal enrichment factors are less than 10 for most metals, in the AS aerosols. In contrast, TSIO Population II ‘open-ocean southern air’ sampled during the southwest monsoon season, has an average Al concentration of only about 10 ng m−3 of air. Trace metal concentrations in the TSIO ‘open-ocean southern air’ during the southwest monsoon season are representative of ‘clean’ remote marine air and are generally similar to those reported over the central North Pacific.Mineral dust concentrations over the Indian Ocean decrease in a north to south direction, from about 15–20 μg m−3 of air in the extreme north to about 0.01–0.25 μg m−3 of air in the far south. The deposition of mineral dust over the northern Arabian Sea can account for approximately 75% of the non-carbonate material incorporated into the underlying sediments.In the Arabian Sea the dissolved atmospheric inputs of all the trace metals, with the exception of Cu and Co, exceed those from fluvial run-off by factors which range from 9.6 for Pb to 1.6 for Cr.  相似文献   

20.
Production of the marine calanoid copepod Acartia steueri was measured from 2 October 1991 to 8 October 1992 at a station in Ilkwang Bay, on the southeastern coast of Korea. Phytoplankton standing stock ranged over 1.0 to 9.3 mg chl.a m−3, and annual primary productivity (by the C-14 method) at three stations was estimated at 200 gC m−2 yr−1. Acartia steueri (nauplii + copepodids + adults) were present in the plankton throughout the year, with seasonal variation in abundance. Biomass of A. steueri, excluding the NI stage, was 0.01–4.55 mgC m−3 (mean: 0.68 mgC m−3) with peaks in November, February, May and July-early August, and relatively low biomass in September– January. Instantaneous growth rates of the nauplius stages were higher than the copepodid stages. Annual production of A. steueri was 25.1 mgC m−3 yr−1 (or 166 mgC m−2 yr−1), showing peaks in November, May and July–August with a small peak in February, and low production in December–April and September–October. There were no significant relationships between the daily production rate of A. steueri and temperature or chlorophyll a concentration, indicating that unknown other factors might be related to the variation of the production rate.  相似文献   

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

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