首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 625 毫秒
1.
To refine the previous hypothesis of turbulent energy spectrum (Nan-niti, 1970), further theoretical consideration and a simple experiment using electric fans was carried out. From the hypothesis ofq=δε/δk andεq=const., the following relations are derived for the spectrumF, and the autocorrelation functionR:F(k)k −(5−2α)/3 andR(t)∼1−(t/τ)2(1−α)/3, whereq is the local energy supply,ε the energy dissipation rate, andk the wavenumber,t the time,τ the lifetime of the largest eddy. The following further relations are also obtained:εl −α,vl β,vτ γ andγ=(1−α)/(2+α)=β/(1−β), wherel is the length scale of the phenomenon,v turbulent velocity,α=1/2,β=1/6, andγ=1/5. The experimental results support these relations.  相似文献   

2.
The atmosphere-ocean exchange of climatically important gases is determined by the transfer velocity (k) and concentration gradient across the interface. Based on observations in the northwestern subarctic Pacific and Sagami Bay, we report here the results of the first ever application of the natural abundance of triple isotopes of dissolved oxygen (16O, 17O and 18O) for direct estimation of k and propose a new relationship with wind speed. The k values estimated from nighttime variations in oxygen isotopes are found to be higher than the direct estimations at low wind speed (<5 m s−1) and lower at high wind speeds (>13 m s−1) and showed significant spatial variability. The method presented here can be used to derive seasonal and spatial variations in k and the influence of surface conditions on the value, leading to improved estimates of biogenic/anthropogenic gas exchange at the air-sea interface.  相似文献   

3.
In the present study, we have investigated the conditions influencing encystment and excystment in the dinoflagellate Gyrodinium instriatum under laboratory conditions. We incubated G. instriatum in modified whole SWM-3 culture medium and in versions of modified SWM-3 from which NO3 , PO4 3−, NO3 + PO4 3−, or Si had been omitted and observed encystment. Percentage encystment was high in the media without N and without P, while the percentage encystment in the medium lacking both N and P was highest. Moreover, to investigate N or P concentration which induced the encystment, Gyrodinium instriatum was also incubated in media with different concentrations of inorganic N and P; the concentrations of NO2 + NO3 and PO4 3− were measured over time. The precursors of cysts appeared within 2 or 3 days of a decrease in NO2 + NO3 or PO4 3− concentration to values lower than 1 μM or 0.2 μM, respectively. When cysts produced in the laboratory were incubated, we observed excystment after 8–37 days, without a mandatory period of darkness or low temperature. We incubated cysts collected from nature at different temperatures or in the dark or light and observed excystments. Natural cysts excysted at temperatures from 10 to 30°C, in both light and dark, but excystment was delayed at low temperatures. These studies indicate that G. instriatum encysts in low N or P concentration and excysts over a wide temperature range, regardless of light conditions, after short dormancy periods.  相似文献   

4.
Uncertainties in global mapping of Argo drift data at the parking level   总被引:1,自引:0,他引:1  
We used Argo float drift data to estimate average ocean currents at 1000 dbar depth from early 2000 to early 2010. Our estimates cover the global oceans, except for marginal seas and ice-covered regions, at a resolution of 1 degree in latitude and longitude. The estimated flow field satisfies the horizontal boundary condition of no flow through the topography, and is in geostrophic balance. We also estimated the uncertainty in the average flow field, which had a typical magnitude of 0.03 ms−1. The uncertainty is relatively large (>0.03 ms−1 in both the zonal and meridional directions) near the Equator and in the Southern Ocean. The array bias, which is the bias due to the horizontal gradient in the spatial density of the float data, is generally negligible, with an average magnitude outside the equatorial region of 0.007 ms−1, becoming relatively large (>0.01 ms−1) only near the coastal regions. The measurement uncertainty is assumed to be spatially uniform and includes errors due to the Argos positioning system, internal clock drift, unknown surface drift before submerging or after surfacing, and unknown drifts during ascent and descent between the surface and the parking depth. We found that the overall uncertainty was not sensitive to the assumed value of the measurement uncertainty (ɛ m )1/2 when (ɛ m )1/2 < 0.01 ms−1 but it increased with (ɛ m )1/2 for (ɛ m )1/2 > 0.01 ms−1.  相似文献   

5.
The mechanism of the formation of horizontal temperature/density gradients above underwater coastal slopes of natural basins due to heating/cooling from the surface is considered. It is shown that the time required for formation of the gradients is rather small (tens of minutes for a thermocline depth of tens of meters), but the development of the corresponding flows may not be accomplished even in a day long cycle. The time dependence of the horizontal water exchange between the shallow and deep areas is analytically treated. The spatial scale of the problem is the main parameter that defines the resulting quasi-stationary value of the flow rate. The joint analysis of the published field, laboratory, and numerical data of many authors in the range of the above-slope depths of 10−2 m < d < 3 × 102 m (dD, where D is the thickness of the upper thermally active layer of a basin) indicates that the relation between the value of the horizontal quasi-stationary volumetric flow rate and the local depth looks like Q = 0.00l3 × d 1.37 (R 2 = 0.96). The horizontal convective water exchange is shown to be generally two-layered, ageostrophic, and exhibits its maximum flow rate at the end of the slope. The inferences agree well with the field data and conclusions of other authors.  相似文献   

6.
The seasonal abundance of the dominant dinoflagellate, Ceratium fusus, was investigated from January 2000 to December 2003 in a coastal region of Sagami Bay, Japan. The growth of this species was also examined under laboratory conditions. In Sagami Bay, C. fusus increased significantly from April to September, and decreased from November to February, though it was found at all times through out the observation period. C. fusus increased markedly in September 2001 and August 2003 after heavy rainfalls that produced pycnoclines. Rapid growth was observed over a salinity range of 24 to 30, with the highest specific rate of 0.59 d−1 measured under the following conditions: salinity 27, temperature 24°C, photon irradiance 600 μmol m−2s−1. The growth rate of C. fusus increased with increasing irradiance from 58 to 216 μmol m−2s−1, plateauing between 216 and 796 μmol m−2s−1 under all temperature and salinity treatments (except at a temperature of 12°C). Both field and laboratory experiments indicated that C. fusus has the ability to grow under wide ranges of water temperatures (14–28°C), salinities (20–34), and photon irradiance (50–800 μmol m−2s−1); it is also able to grow at low nutrient concentrations. This physiological flexibility ensures that populations persist when bloom conditions come to an end.  相似文献   

7.
A Geographic Information System (GIS)-aided flow-tracking technique was adopted to investigate nutrient exchange rates between specific benthic communities and overlying seawater in a fringing reef of Ishigaki Island, subtropical Northwestern Pacific. Net exchange rates of NO3 , NO2 , NH4 +, PO4 3−, Total-N and Total-P were estimated from concentration changes along the drogue trajectories, each of which was tracked by the Global Positioning System and plotted on a benthic map to determine the types of benthic habitat over which the drogue had passed. The observed nutrient exchange rates were compared between 5 typical benthic zones (branched-coral (B) and Heliopora communities (H), seaweed-reefrock zone (W), bare sand area (S), and seagrass meadow (G)). The dependence of nutrient exchange rates on nutrient concentrations, physical conditions and benthic characteristics was analyzed by multiple regression analysis with the aid of GIS. The spatial correlation between nutrient exchange rates and benthic characteristics was confirmed, especially for NO3 and PO4 3−, which were usually absorbed in hydrographically upstream zones B and W and regenerated in downstream zones H and G. NO3 uptake in zones B and W was concentration-dependent, and the uptake rate coefficient was estimated to be 0.58 and 0.67 m h−1, respectively. Both nutrient uptake in zone W and regeneration in zone H were enhanced in summer. The net regeneration ratio of NO3 /PO4 3− in zone H in summer ranged 5.2 to 34 (mean, 17.4), which was somewhat higher than previously measured NO3 /PO4 3− for sediment pore waters around this zone (1.1–8.5). Nutrient exchanges in zone S were relatively small, indicating semi-closed nutrient cycling at the sediment-water interface of this zone. NH4 + efflux from sediments was suggested in zone G. The data suggest that the spatial pattern of nutrient dynamics over the reef flat community was constrained by zonation of benthic biota, and that abiotic factors such as nutrient concentrations and flow rates, influenced nutrient exchange rates only in absorption-dominated communities such as zones B and W.  相似文献   

8.
Based on the experimental result, ε∝l −0.5, previously obtained, a hypothesis of turbulent energy spectrum in steady state is presented by introducing local energy supply,q. The foundamental relations areq=δε/δκ and εq=const., wherek is the wave number. Meteorological Research Institute  相似文献   

9.
The mixing processes in the Mixed Water Region (MWR) that lead to changes in the properties of North Pacific Intermediate Water (NPIW) have been studied using observational data sets obtained in May–June 1998. Neutral surfaces, the equation of water mass conversion rate on neutral surfaces and the equation of vertical velocity across neutral surfaces have been used to distinguish dominant processes by assuming the horizontal scale to be the streamer scale (under 100 km). The possibility of double diffusive convection is also discussed in relation to the density ratio. These results may be summarized as follows: (1) the difference between the potential density surface and the neutral surface may rise to −0.04 kg/m3 around the source water of NPIW; (2) horizontal diffusion causes strong modifications of the source water of NPIW; (3) the density range within which strong modification of the source water of NPIW occurs becomes dense from the northern part of MWR near the Oyashio Front to the southern part near the Kuroshio Front, and to the eastern part. Our modeling of these processes shows that cabbeling has effects on the density increment of the source water of NPIW in the northern and southern part of MWR. Double diffusive convection has effects on the density increment of the source water of NPIW, mainly in the northern part of MWR. The possible density increment due to cabbeling in these areas is estimated to be 0.01≈0.03 kg/m3. The possible density increment due to double diffusive convection is 0.01≈0.03 kg/m3. The total density increment due to cabbeling and double diffusive convection amounts to 0.06 kg/m3. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

10.
Two-dimensional temperature data observed by use of a 275 meter towed thermistor chain deployed from an oceanographic research vessel USS MARYSVILLE, which cruised with a speed of 6.2 knots in July 1966 across the Kuroshio Extension in the North Pacific, are investigated. Two-dimensional variations of the distribution of the isotherms along the ship's track are analyzed with special reference to their slope, wavelength and wave height. The results show that the slope and wave height of isotherms have a tendency to increase as the temperature decreases. Even if the contribution of wave heights smaller than 1.5 m is neglected, i.e., contribution of large scale slope with a horizontal scale of 5–30 km is subtracted, this tendency is still detected. In contrast to this, the wavelength evaluated by the crest to crest method has no dependency on the temperature. Power spectrum of the isotherm depth is proportional tok –1.87 for 13°C andk –2.13 for 27°C, wherek is the wave number. It is shown that the spectra of warmer isotherms are relatively well approximated by –2 power law (Garrett and Munk spectrum) for internal waves rather than the –5/3 power law (Kolmogorov spectrum) for three dimensional isotropic turbulence.  相似文献   

11.
We measured the ammonium excretion, phosphate excretion and respiration rates of the scyphomedusa Aurelia aurita from Ondo Strait, in the central part of the Inland Sea of Japan, at 28 and 20°C. The rates measured at 28°C were converted to those at 20°C using the Q10 values, i.e. 1.56, 1.57 and 2.80, for ammonium excretion, phosphate excretion and respiration rates, respectively. The composite relationships between metabolic rates and wet weight of a medusa (WW, g, range 11–1330 g) at 20°C were expressed by the following allometric equations. For ammonium excretion rate (N, μmoles N medusa−1d−1): N = 0.497WW 1.09, phosphate excretion rate (P, μmoles P medusa−1d−1): P = 0.453WW 0.84, and respiration rate (R, μmoles O2 medusa−1d−1): R = 96.9WW 1.06. Mean O:N ratios (i.e. atomic ratios of 16.9 and 11.0 at 28 and 20°C, respectively) indicated that the metabolism of A. aurita medusae was protein-dominated. These metabolic parameters enabled us to estimate the nitrogen and phosphorus regeneration rates of an A. aurita medusa population typical of early summer in the Ondo Strait (means of water temperature, medusa individual weight and population biomass: 20°C, 200 g WW and 50.8 g WW m−3, respectively). Regenerated nitrogen and phosphorus were equivalent to 10.0 and 21.6% of phytoplankton uptake rates, respectively, nearly twice that estimated for mesozooplankton, demonstrating that A. aurita medusae are key components of the plankton community, influencing the trophic and nutrient dynamics in the Ondo Strait during early summer.  相似文献   

12.
To understand the processes transporting nitrate to the surface layer of the western and central equatorial Pacific, we measured the nitrogen isotopic ratio of nitrate (δ 15NO 3 ), which is a very useful tracer of the source of nitrate, above 200 m depth in this region in December 1999. δ 15NO 3 is higher (about 13.0‰) in the surface water than in the subsurface water (where it is about 6.5‰) due to isotopic fractionation during nitrate uptake by phytoplankton. The δ 15NO 3 value has a roughly linear relationship with the natural logarithm of nitrate concentration (ln[NO 3 ]). However, for values above 150 m depth, the intercept of this linear relationship varies with position from east to west. On the other hand, the data at 200 m depth at all observation stations are concentrated around a single point (ln[NO 3 ] = 2.5 and δ 15NO 3 = 6.5‰) and do not fit the linear relationships for the shallower values. To examine the meaning of the observed distributions of δ 15NO 3 and nitrate concentration we developed a box model including nitrogen and nitrogen isotopic cycles. By reproducing the observed relationship between δ 15NO 3 and nitrate concentration using this model we found that most nitrate is transported horizontally from the eastern equatorial Pacific. We also conducted case studies and investigated the effects of differences in pathways of nitrate transport on the distributions of δ 15NO 3 and nitrate concentration. From these studies we concluded that the observed linear relationships between δ 15NO 3 and ln[NO 3 ], having a common slope around 6‰ but different intercepts at each station, are evidence of the significant horizontal transport of nitrate to the surface water in this area.  相似文献   

13.
In order to determine quantitatively the reason for the high productivity in the Oyashio Region, which is the southwest part of the Pacific Subarctic Region, the annual-mean vertical circulation of nitrogen in the region was estimated from the vertical profiles of nitrate, dissolved oxygen and salinity, and sediment-trap data by adapting them to the balance equations. Estimates of the upwelling velocity (1.7×10−5cm sec−1) and the vertical diffusivity (2.1 cm2 sec−1) in the abyssal zone and the primary and secondary productivities (44 and 4 mgN m−2day−1, respectively) in the euphotic zone were close to those of previous works. The estimated vertical circulation of nitrogen strongly suggested that, since the divergence (5 mgN m−2day−1) is caused by the abyssal convergence (6 mgN m−2day−1) and the positive precipitation, the local new production (22 mgN m−2day−1) necessarily exceeds not only the sinking flux (10 mgN m−2day−1) itself but also the sum of the sinking flux and the downward diffusion of dissolved and particulate organic matter (7 mgN m−2day−1) produced probably in the euphotic zone. The important roles of the abyssal circulation, the winter convection, and the metabolic activity in the bathyal zone to support the high productivity in the euphotic zone were clarified quantitatively.  相似文献   

14.
Seasonal and interannual variations in physicochemical properties were investigated in the neritic area of Sagami Bay, Kanagawa, Japan, from December 2000 to December 2005. Physicochemical properties (i.e. temperature, salinity, density, dissolved oxygen and dissolved inorganic nutrient concentration) revealed clear seasonal variations, which were similar to each other during all 5 years. Temperature, salinity and dissolved inorganic nutrients showed rapid, drastic variations within a few days and/or weeks. These variations are related to sea levels, principally due to the shifting effects of the Kuroshio Current axis: they were strongly affected by the Kuroshio Water and other waters, when sea level difference was greater than ca. 35 cm and lower than ca. 15 cm, respectively. Temperature difference (DF T ) increased with sea level difference, and the difference of salinity and dissolved inorganic nutrients (NH4 +-N, NO3 +NO2 -N, NH4 ++NO3 +NO2 -N, PO4 3−-P and SiO2-Si) increased and decreased with DF T , respectively. All these correlations are significant. Total dissolved inorganic nitrogen (N), phosphate (P) and silicate (Si) revealed seasonal variations in the ranges of 0.57–16.08, 0.0070–0.91 and 0.22–46.38 μM, respectively. From the regression equations between these elements allowed the following relation to be obtained; Si:N:P = 14.8:13.4:1. Dissolved inorganic nutrients were characterized by Si and/or P deficiency, especially in the upper layer (0–20 m depth) during summer. Single and/or combined elements are discussed on the basis of potential and stoichiometric nutrient limitations, which could restrict phytoplankton (diatom) growth as a limiting factor.  相似文献   

15.
The biochemical effects of a cold-core eddy that was shed from the Kuroshio Current at the Luzon Strait bordering the South China Sea (SCS) were studied in late spring, a relatively unproductive season in the SCS. The extent of the eddy was determined by time-series images of SeaWiFS ocean color, AVHRR sea surface temperature, and TOPEX/Jason-1 sea surface height anomaly. Nutrient budgets, nitrate-based new production, primary production, and phytoplankton assemblages were compared between the eddy and its surrounding Kuroshio and SCS waters. The enhanced productivity in the eddy was comparable to wintertime productivity in the SCS basin, which is supported by upwelled subsurface nitrate under the prevailing Northeastern Monsoon. There were more Synechococcus, pico-eucaryotes, and diatoms, but less Trichodesmium in the surface water inside the eddy than outside. Prochlorococcus and Richelia intracellularis showed no spatial differences. Water column-integrated primary production (IPP) inside the eddy was 2–3 times that outside the eddy in the SCS (1.09 vs. 0.59 g C m−2d−1), as was nitrate-based new production (INP) (0.67 vs. 0.25 g C m−2d−1). INP in the eddy was 6 times that in the Kuroshio (0.12 g C m−2d−1). IPP and INP in the eddy were higher than the maximum production values ever measured in the SCS basin. Surface chlorophyll a concentration (0.40 mg m−3) in the eddy equaled the maximum concentration registered for the SCS basin and was higher than the wintertime average (0.29 ± 0.04 mg m−3). INP was 3.5 times as great and IPP was doubled in the eddy compared to the wintertime SCS basin. As cold core eddies form intermittently all year round as the Kuroshio invades the SCS, their effects on phytoplankton productivity and assemblages are likely to have important influences on the biogeochemical cycle of the region.  相似文献   

16.
Concentrations of particulate organic nitrogen (PN), dissolved inorganic nitrogen (DIN), and their nitrogen isotope ratios (δ 15N) in the Kiso-Sansen Rivers were determined from monthly observations over the course of a year to assess variations in the form and sources of riverine nitrogen discharged into Ise Bay. The δ 15N values of NO3 observed in the Kiso-Sansen Rivers showed a logarithmic decreasing trend from 8 to 0‰, which varied with the river discharge, indicating mixing between point sources with high δ 15N and non-point sources with low δ 15N. The influence of isotope fractionation of in situ biogeochemical processes (mainly DIN assimilation by phytoplankton) on δ 15N of NO3 was negligible, because sufficient concentrations of NH4 + for phytoplankton demand would inhibit the assimilation of NO3 . A simple relationship between river discharge and δ 15N of NO3 showed that the fraction of total NO3 flux arising from point sources increased from 4.0–6.3% (1.1–1.8 tN day−1) during higher discharge (>600 m3 s−1) to 30.2–48.3% (2.6–4.1 tN day−1) during lower discharge (<300 m3 s−1). Riverine NO3 discharge from the Kiso-Sansen Rivers can explain 75% of the variations in surface NO3 at the head of Ise Bay over the year.  相似文献   

17.
INTRODUCTIONTherehavebeenmanystudiesandcomputationsonVToftheKuroshiointheEastChinaThisprojectwassupportedbytheNationalNaturalScienceFoundationofChinaundercontractNo.49476278.Asanditsvacation.Forexample,(1)basedonhydrographicobservationsatactionG(PN)f...  相似文献   

18.
Compared to oxygen isotopes, the carbon isotope composition of biogenic carbonates is less commonly used as proxy for palaeoenvironmental reconstructions because shell δ13C is derived from both dissolved inorganic (seawater) and organic carbon sources (food), and interactions between these two pools make it difficult to unambiguously identify any independent effect of either. The main purpose of this study was to demonstrate any direct impact of variable food supply on bivalve shell δ13C signatures, using low/high rations of a 13C-light mixed algal diet fed to 14-month-old (adult) cultured Japanese Crassostrea gigas under otherwise essentially identical in vitro conditions during 3 summer months (May, June and July 2003, seawater temperature means at 16, 18 and 20 °C respectively) in experimental tanks at the Argenton laboratory along the Brittany Atlantic coast of France. At a daily ration of 12% (versus 4%) oyster dry weight, the newly grown part of the shells (hinge region) showed significantly lower δ13C values, by 3.5‰ (high ration: mean of −5.8  ± 1.1‰, n = 10; low ration: mean of −2.3  ± 0.7‰, n = 6; ANOVA Scheffe’s test, p < 0.0001). This can be explained by an enhanced metabolic activity at higher food supply, raising 13C-depleted respiratory CO2 in the extrapallial cavity. Based on these δ13C values and data extracted from the literature, and assuming no carbon isotope fractionation between food and shell, the proportion of shell metabolic carbon would be 26  ± 7 and 5  ± 5% for the high- and low-ration C. gigas shells respectively; with carbon isotope fractionation (arguably more realistic), the corresponding values would be 69  ± 14 and 24  ± 9%. Both groups of cultured shells exhibited lower δ13C values than did wild oysters from Marennes-Ol éron Bay in the study region, which is not inconsistent with an independent influence of diet type. Although there was no significant difference between the two food regimes in terms of δ18O shell values (means of 0.1  ± 0.3 and 0.4  ± 0.2‰ at high and low rations respectively, non-significant Scheffe’s test), a positive δ13C vs. δ18O relationship recorded at high rations supports the interpretation of a progressive temperature-mediated rise in metabolic activity fuelled by higher food supply (in this case reflecting increased energy investment in reproduction), in terms not only of δ13C (metabolic signal) but also of δ18O (seawater temperature signal). Overall, whole-shell δ18O trends faithfully recorded summer/winter variations in seawater temperature experienced by the 17-month-old cultured oysters.  相似文献   

19.
In Tokyo Bay the concentrations of dissolved gaseous mercury (DGM) in the surface seawater and total gaseous mercury (TGM) over the sea were measured during December 2003, October 2004 and January 2005. Based on these data, the evasional fluxes of mercury from the sea surface were estimated using a gas exchange model. In addition, an automatic wet and dry deposition sampler was used to measure the wet and dry depositional fluxes of mercury from December 2003 to November 2004 at three locations in and near Tokyo Bay. The results indicate that the average DGM and TGM levels of seven locations are 52 ± 26 ng m−3 and 1.9 ± 0.6 ng m−3, respectively, which shows that the surface seawater in Tokyo Bay is supersaturated with gaseous mercury, leading to an average mercury evasional flux of 140 ± 120 ng m−2d−1. On the other hand, the annual average wet and dry depositional fluxes of mercury at three locations were 19 ± 3 μg m−2yr−1 and 20 ± 9 μg m−2yr−1, respectively. These depositional fluxes correspond to the daily average total depositional flux of 110 ± 20 ng m−2d−1. Thus, it is suggested that in Tokyo Bay, the evasional fluxes of mercury are comparable to the depositional fluxes.  相似文献   

20.
The egg production rate (EPR) of Calanus sinicus was measured from March 2007 to April 2010 at three stations along the Korean coast of the Yellow Sea (in coastal waters off Saemangum, Yeongheungdo, and Asan Bay) to estimate in situ maximum egg production rate (MEPR) and to understand whether the females were limited in their growth or fecundity in the field. The mean EPR of C. sinicus at each sampling date ranged from 10.3–34.9 eggs female−1 d−1 (mean 23.4 eggs female−1 d−1), and the EPR of individual copepods ranged from 0–81 eggs female−1 d−1. The mean EPR was positively correlated with the body weight of female copepods. The MEPR at each sampling date ranged from 40–81 eggs female−1 d−1 (mean 50.4 eggs female−1 d−1). Over 84% of eggs spawned hatched successfully. The weight-specific growth rate (WSGR) ranged from 0.038–0.111 d−1 (mean 0.082 d−1), indicating that 3.8–11.1% of the carbon in an adult female was produced daily as female growth. The WSGR was negatively correlated with water temperature. The ratio of mean EPR to observed mean MEPR ranged from 20–70% (mean 46%), indicating that ∼54% of a female’s growth might be limited in the field. We suggest that the ratio of observed EPR to mean MEPR of copepod can be applied to understand how the copepod responds to environmental changes, as well as EPR and hatching success.  相似文献   

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

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