首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
海浪破碎对海洋上混合层中湍能量收支的影响   总被引:2,自引:1,他引:2  
海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。  相似文献   

2.
刘子龙  史剑  蒋国荣 《海洋科学》2017,41(3):122-129
基于海浪模式WAVEWATCH Ⅲ模拟北太平洋海浪要素,结合NDBC浮标资料进行验证,发现模拟出的有效波高与浮标测量值具有很好的一致性。基于改进型白冠覆盖率耗散模型,利用海浪模式模拟出的有效波高、有效波周期和摩擦速度等海浪要素计算出单位面积水柱内因海浪破碎产生的湍动能通量。通过改变环流模式sbPOM湍动能方程的上边界条件,引入海浪破碎产生的湍动能通量,并探究海浪破碎对北太平洋海表面温度模拟的影响。研究表明,由于海浪破碎的引入,环流模式sbPOM对北太平洋海表面温度模拟的准确程度得到提升,这为大气模式提供一个准确的北太平洋下边界条件具有重要意义。  相似文献   

3.
A new three-dimensional numerical model is derived through a wave average on the primitive N-S equations, in which both the"Coriolis-Stokes forcing" and the"Stokes-Vortex force" are considered. Three ideal experiments are run using the new model applied to the Princeton ocean model (POM). Numerical results show that surface waves play an important role on the mixing of the upper ocean. The mixed layer is enhanced when wave effect is considered in conjunction with small Langmuir numbers. Both surface wave breaking and Stokes production can strengthen the turbulent mixing near the surface. However, the influence of wave breaking is limited to a thin layer, but Stokes drift can affect the whole mixed layer. Furthermore, the vertical mixing coefficients clearly rise in the mixed layer, and the upper ocean mixed layer is deepened especially in the Antarctic Circumpolar Current when the model is applied to global simulations. It indicates that the surface gravity waves are indispensable in enhancing the mixing in the upper ocean, and should be accounted for in ocean general circulation models.  相似文献   

4.
Effect of Stokes drift on upper ocean mixing   总被引:1,自引:0,他引:1  
Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2. 5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, ilffluence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.  相似文献   

5.
Effect of Langmuir circulation on upper ocean mixing in the South China Sea   总被引:2,自引:0,他引:2  
Effect of Langmuir circulation (LC) on upper ocean mixing is investigated by a two-way wave-current coupled model. Themodel is coupled of the ocean circulationmodel ROMS (regional ocean modeling system) to the surface wave model SWAN (simulating waves nearshore) via the model-coupling toolkit. The LC already certified its importance by many one-dimensional (1D) research andmechanismanalysis work. This work focuses on inducing LC’s effect in a three-dimensional (3-D) model and applying it to real field modeling. In ROMS, theMellor-Yamada turbulence closuremixing scheme is modified by including LC’s effect. The SWAN imports bathymetry, free surface and current information fromthe ROMS while exports significant wave parameters to the ROMS for Stokes wave computing every 6 s. This coupled model is applied to the South China Sea (SCS) during September 2008 cruise. The results show that LC increasing turbulence and deepening mixed layer depth (MLD) at order of O (10 m) in most of the areas, especially in the north part of SCS where most of our measurements operated. The coupled model further includes wave breaking which will bringsmore energy into water. When LC works together with wave breaking,more energy is transferred into deep layer and accelerates the MLD deepening. In the north part of the SCS, their effects aremore obvious. This is consistent with big wind event in the area of the Zhujiang River Delta. The shallow water depth as another reasonmakes themeasy to influence the oceanmixing as well.  相似文献   

6.
A new model developed from the full-spectrum model of Lyzenga and Bennett ( 1988 ) is built up by us preliminarily through considering the impact of the ocean surface mixed layer turbulence on SAR remote sensing of internal solitary waves. In the partial differential equation of the action spectral density of the surface gravity-capillary waves the source function representing the contribution of the turbulence is added besides the usual source function representing the contribution of the wind. The source function is determined by applying the κ - ε model and adopting the Nasmyth spectrum of oceanic turbulence ( Nasmyth, 1970; Oakey, 1982; Fan, 2002) on the basis of the previous simulation results of internal tide transformation obtained in the companion paper (Fan et al. ,2008). Either under relatively high wind speed, or under low wind speed, our model predicts significant large modulations of radar backscatter at all three bands ( L, C and X bands) for both VV and HH polarization. These results prove that considering the impact of ocean surface mixed layer turbulence on SAR remote sensing of internal solitary waves is reasonable and appropriate for settlement of the well-known problem of contemporary radar imaging models.  相似文献   

7.
The mean seasonal cycle of mixed layer depth (MLD) in the extratropical oceans has the potential to influence temperature, salinity and mixed layer depth anomalies from one winter to the next. Temperature and salinity anomalies that form at the surface and spread throughout the deep winter mixed layer are sequestered beneath the mixed layer when it shoals in spring, and are then re-entrained into the surface layer in the subsequent fall and winter. Here we document this ‘re-emergence mechanism’ in the North Pacific Ocean using observed SSTs, subsurface temperature fields from a data assimilation system, and coupled atmosphere–ocean model simulations. Observations indicate that the dominant large-scale SST anomaly pattern that forms in the North Pacific during winter recurs in the following winter. The model simulation with mixed layer ocean physics reproduced the winter-to-winter recurrence, while model simulations with observed SSTs specified in the tropical Pacific and a 50 m slab in the North Pacific did not. This difference between the model results indicates that the winter-to-winter SST correlations are the result of the re-emergence mechanism, and not of similar atmospheric forcing of the ocean in consecutive winters. The model experiments also indicate that SST anomalies in the tropical Pacific associated with El Niño are not essential for re-emergence to occur.The recurrence of observed SST and simulated SST and SSS anomalies are found in several regions in the central North Pacific, and are quite strong in the northern (>50°N) part of the basin. The winter-to-winter autocorrelation of SSS anomalies exceed those of SST, since only the latter are strongly damped by surface fluxes. The re-emergence mechanism also has a modest influence on MLD through changes in the vertical stratification in the seasonal thermocline.  相似文献   

8.
Dissipation rate statistics in the near-surface layer of the ocean were obtained during the month-long COARE Enhanced Monitoring cruise with a microstructure sensor system mounted on the bow of the research vessel. The vibration contamination was cancelled with the Wiener filter. The experimental technique provides an effective separation between surface waves and turbulence, using the difference in spatial scales of the energy-containing surface waves and small-scale turbulence. The data are interpreted in the coordinate system fixed to the ocean surface. Under moderate and high wind-speed conditions, we observed the average dissipation rate of the turbulent kinetic energy in the upper few meters of the ocean to be 3–20 times larger than the logarithmic layer prediction. The Craig and Banner (J. Phys. Oceanogr. 24 (1994) 2546) model of wave-enhanced turbulence with the surface roughness length from the water side z0 parameterized according to the Terray et al. (J. Phys. Oceanogr. 26 (1996) 792) formula z0=cHs provides a reasonable fit to the experimental dissipation profile, where z is the depth (defined here as the distance to the ocean surface), c≈0.6, and Hs is the significant wave height. In the wave-stirred layer, however, the average dissipation profile deviates from the model (supposedly because of extensive removing of the bubble-disturbed areas close to the ocean surface). Though the scatter of individual experimental dissipation rates (10-min averages) is significant, their statistics are consistent with the Kolmogorov's concept of intermittent turbulence and with previous studies of turbulence in the upper ocean mixed layer.  相似文献   

9.
本文通过理想化的外部强迫以及海洋站点实测数据驱动普林斯顿海洋模式来研究海洋热力学效应和斯托克斯漂流对上混合层数值模拟的影响。在Mellor-Yamada湍流闭合方案中,经常出现夏季海表面温度偏暖和混合层深度偏浅的模拟误差。实验表明,斯托克斯漂流在冬季和夏季均能增强湍流动能,加深混合层深度。这种效应可以改善夏季的模拟结果,但与观测数据相比,将增大冬季混合层深度的模拟误差。斯托克斯漂流可以通过增强湍动能来加深混合层深度。结果表明,将斯托克斯漂流与冷皮层和暖层对上部混合层的热效应相结合,可以正确地模拟混合层深度。在夏季,海洋冷皮层和暖层通过“阻挡结构”和双温跃层结构模拟出更真实的上混合层变化。在冬季,海洋热力学效应通过增强上层海洋层结平衡了斯托克斯漂流的影响,并且由斯托克斯漂流引起的过度混合被校正。  相似文献   

10.
Based on the theoretical spectral model of inertial internal wave breaking(fine structure) proposed previously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior below the surface mixed layer in the ocean general circulation model(OGCM) is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes(including inertial internal wave breaking product) on oceanic interior mixing is emphasized. We suggest that adding the inertial internal wave breaking mixing scheme(F-scheme for short) put forward in this paper to the turbulence mixing scheme of Canuto et al.( T-scheme for short) in the OGCM, except the region from 15°S to 15°N. The numerical results of F-scheme by using WOA09 data and an OGCM(LICOM, LASG/IAP climate system ocean model) over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation(AMOC) by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.  相似文献   

11.
For settlement of the well-known problem of contemporary radar imaging models,i.e.,the pmblem of a general underestimation of radar signatures of hydrodynamic features over oceanic internal waves and underwater bottom topography in tidal watels at at high radar frequency bands(X-band and C-band),the impact of the ocean surface mixed layer turbulence and the significance of strat-ified oceanic model on SAR remote sensing of internal solitary waves are proposed.In the north of the South China Sea by utilizing seme observed data of background field the nonlinearity coefficient,the dispersion coefficient,the horizontal variability coefficient and the phase speed in the generalized K-dV equation are determined approximately.Through simulations of internal tide transfor-mation the temporal evolution and spatial distribution of the vertical displacement and horizontal velocity of intereal wave field are obtained.The simulation results indicate that the maximum amplitudes of internal solitary waves occur at depth 35 m,but the maximum current speeds take place at depth 20 m in this area of the sea(about 20°30'N,114°E)in August.It was noticed that considering the effects of flood current and ebb current respectively is appropriate to investigate influence of the background shear flow on coefficients of the K-dV equation.The obtained results provide the possibility for the simulation of SAR signatures of inter-nal solitary waves under considering the impact of ocean surface mixed layer turbulence in the companion paper.  相似文献   

12.
Ocean temperature responses to Typhoon Mstsa in the East China Sea   总被引:1,自引:1,他引:0  
The MASNUM wave-tide-circulation coupled model, with 21 layers in the vertical and (1/8) °horizontal resolution, was employed to investigate the oceanic responses to Typhoon Mstsa which traversed the East China Sea (ECS) during the period of 4 - 6 August, 2005. Numerical experiment results are analyzed and compared with observation. The responses of the sea surface temperature (SST), in a focused area of (27° -29°N, 121° - 124°E), include heating and cooling stages. The heating is mainly due to warm Kuroshio water transportation and downwelling due to the water accumulation. In the cooling stage, the amplitude of the simulated cold wake ( -3℃ ), located on the right side of this typhoon track, is compared quite well with that of the satellite observed SST data. The wave-induced mixing(Bv) plays a key role for the SST cooling. Bv still plays a leading role, which accounts for 36%, for the ocean temperature drop in the upper ocean of 0 - 40 m, while the upwelling is responsible for 84% of the cooling for the lower layer of 40 - 70 m. The mixed layer depth (MLD) increased quickly from 28 to 50 m in the typhoon period. However, the simulated MLD without the wave-induced vertical mixing, evolution from 13 to 32 m, was seriously underestimated. The surface wave is too important to be ignored for the ocean responses to a typhoon.  相似文献   

13.
To explore the causes of the winter shallow mixed layer and high sea surface temperature (SST) along the strong Kuroshio jet from the East China Sea to the upstream Kuroshio extension (25.5°N–150°E) during 1988–1994 when the Japanese sardine stocks collapsed, high-resolution ocean general circulation model (OGCM) hindcast data are analyzed with a bulk mixed layer model which traces particles at the mixed layer base. The shallow mixed layer and high SST along the Kuroshio jet are mainly caused by the acceleration of the Kuroshio current velocity and the reduction of the surface cooling. Because the acceleration reduces the time during which the mixed layer is exposed to wintertime cooling, deepening and cooling of the winter mixed layer are restricted. The weaker surface cooling due to less severe meteorological forcing also causes the shallow mixed layer and the high SST. The impact of the strong heat transport along the Kuroshio extends to the southern recirculation gyre of the Kuroshio/Kuroshio extension regions; previous indications that the Japanese sardine recruitment is correlated with the winter SST and the mixed layer depth (MLD) in the Kuroshio extension recirculation region could be related to the velocity, SST, and MLD near the Kuroshio axis which also could affect the variability of North Pacific subtropical water.  相似文献   

14.
《Coastal Engineering》2001,42(1):53-86
A numerical model is used to simulate wave breaking, the large scale water motions and turbulence induced by the breaking process. The model consists of a free surface model using the surface markers method combined with a three-dimensional model that solves the flow equations. The turbulence is described by large eddy simulation where the larger turbulent features are simulated by solving the flow equations, and the small scale turbulence that is not resolved by the flow model is represented by a sub-grid model. A simple Smagorinsky sub-grid model has been used for the present simulations. The incoming waves are specified by a flux boundary condition. The waves are approaching in the shore-normal direction and are breaking on a plane, constant slope beach. The first few wave periods are simulated by a two-dimensional model in the vertical plane normal to the beach line. The model describes the steepening and the overturning of the wave. At a given instant, the model domain is extended to three dimensions, and the two-dimensional flow field develops spontaneously three-dimensional flow features with turbulent eddies. After a few wave periods, stationary (periodic) conditions are achieved. The surface is still specified to be uniform in the transverse (alongshore) direction, and it is only the flow field that is three-dimensional.The turbulent structures are investigated under different breaker types, spilling, weak plungers and strong plungers. The model is able to reproduce complicated flow phenomena such as obliquely descending eddies. The turbulent kinetic energy is found by averaging over the transverse direction. In spilling breakers, the turbulence is generated in a series of eddies in the shear layer under the surface roller. After the passage of the roller the turbulence spreads downwards. In the strong plunging breaker, the turbulence originates to a large degree from the topologically generated vorticity. The turbulence generated at the plunge point is almost immediately distributed over the entire water depth by large organised vortices. Away from the bed, the length scale of the turbulence (the characteristic size of the eddies resolved by the model) is similar in the horizontal and the vertical direction. It is found to be of the order one half of the water depth.  相似文献   

15.
仇颖  阳德华  李爽 《海洋科学》2019,43(11):103-110
Langmuir环流影响着海洋上层的能量输入,对海洋上混合层的形成和加深起着重要作用,对于海洋上混合层具有重要意义。近年来许多学者采用大涡模拟(LES)方法对Langmuir环流进行机制研究,并通过在雷诺平均模型中参数化Langmuir环流效应,将Langmuir环流过程引入到三维海洋环流或海洋耦合模式中,提出了一系列混合参数化方案。本文回顾了Langmuir环流在雷诺平均模式参数化中的研究进展,主要可分为以下几种方案:一种方法是用Langmuir数在KPP垂直混合参数化方案中引入湍流特征速度增强因子,并不断发展Langmuir数的定义;一种是在Mellor-Yamada2.5湍流闭合模型中增加斯托克斯漂流剪切效应项,此外还有通过修改模式中混合长方程来加入Langmuir效应等。通过在雷诺平均模式中应用的结果来看,现有的参数化方案在一定程度上改善了混合层深度和SST的模拟,肯定了Langmuir环流在加深混合等方面的作用,但仍存在一些问题需要在今后的研究中进一步改进。  相似文献   

16.
The significant underestimation of sea surface temperature (SST) and the temperature in the upper ocean is one of common problems in present climate models. The influence of the wave-induced mixing on SST and the temperature in the upper ocean was examined based on a global climate model. The results from the model coupled with wave-induced mixing showed a significant improvement in the simulation of SST and the temperature in the upper ocean compared with those of the original model without wave effects. Although there has still a cold bias, the new simulation is much closer to the climatology, especially in the northern ocean and tropical ocean. This study indicates that some important physical processes in the accurate simulation of the ocean may be ignored in present climate models, and the wave-induced mixing is one of those factors. Thus, the wave-induced mixing ( or the effect of surface waves) should be incorporated properly into climate models in order to simulate or forecast the ocean, then climate system, more accurately.  相似文献   

17.
The effects of biological heating on the upper-ocean temperature of the global ocean are investigated using two ocean-only experiments forced by prescribed atmospheric fields during 1990–2007, on with fixed constant chlorophyll concentration, and the other with seasonally varying chlorophyll concentration. Although the existence of high chlorophyll concentrations can trap solar radiation in the upper layer and warm the surface, cooling sea surface temperature (SST) can be seen in some regions and seasons. Seventeen regions are selected and classified according to their dynamic processes, and the cooling mechanisms are investigated through heat budget analysis. The chlorophyll-induced SST variation is dependent on the variation in chlorophyll concentration and net surface heat flux and on such dynamic ocean processes as mixing, upwelling and advection. The mixed layer depth is also an important factor determining the effect. The chlorophyll-induced SST warming appears in most regions during the local spring to autumn when the mixed layer is shallow, e.g., low latitudes without upwelling and the mid-latitudes. Chlorophyll-induced SST cooling appears in regions experiencing strong upwelling, e.g., the western Arabian Sea, west coast of North Africa, South Africa and South America, the eastern tropical Pacific Ocean and the Atlantic Ocean, and strong mixing (with deep mixed layer depth), e.g., the mid-latitudes in winter.  相似文献   

18.
《Ocean Modelling》2004,6(3-4):285-334
We have developed a general 1-D multi-component ecosystem model that incorporates a skillful upper ocean mixed layer model based on second moment closure of turbulence. The model is intended for eventual incorporation into coupled 3-D physical–biogeochemical ocean models with potential applications to modeling and studying primary productivity and carbon cycling in the global oceans as well as to promote the use of chlorophyll concentrations, in concert with satellite-sensed ocean color, as a diagnostic tool to delineate circulation features in numerical circulation models. The model is nitrogen-based and the design is deliberately general enough and modular to enable many of the existing ecosystem model formulations to be simulated and hence model-to-model comparisons rendered feasible. In its more general form (GEM10), the model solves for nitrate, ammonium, dissolved nitrogen, bacteria and two size categories of phytoplankton, zooplankton and detritus, in addition to solving for dissolved inorganic carbon and total alkalinity to enable estimation of the carbon dioxide flux at the air–sea interface. Dissolved oxygen is another prognostic variable enabling air–sea exchange of oxygen to be calculated. For potential applications to HNLC regions where productivity is constrained by the availability of a trace constituent such as iron, the model carries the trace constituent as an additional prognostic variable. Here we present 1-D model simulations for the Black Sea, Station PAPA and the BATS site. The Black Sea simulations assimilate seasonal monthly SST, SSS and surface chlorophyll, and the seasonal modulations compare favorably with earlier work. Station PAPA simulations for 1975–1977 with GEM5 assimilating observed SST and a plausible seasonal modulation of surface chlorophyll concentration also compare favorably with earlier work and with the limited observations on nitrate and pCO2 available. Finally, GEM5 simulations at BATS for 1985–1997 are consistent with the available time series. The simulations suggest that while it is generally desirable to employ a comprehensive ecosystem model with a large number of components when accurate depiction of the entire ecosystem is desirable, as is the prevailing practice, a simpler formulation such as GEM5 (N2PZD model) combined with assimilation of remotely sensed SST and chlorophyll concentrations may suffice for incorporation into 3-D prediction models of primary productivity, upper ocean optical clarity and carbon cycling.  相似文献   

19.
In order to investigate the formation mechanism of rapid decrease of maritime sea surface temperature (SST) observed by R/V Keifu Maru, the ocean response to Typhoon Rex is simulated using a mixed layer model. The rapid decrease of the maritime SST is successfully simulated with realistic atmospheric forcing and an entrainment scheme of which sources of turbulent kinetic energy (TKE) are production due to wind stress, generation during free convection, and production due to current shear. The rapid decrease at the observed station by R/V Keifu Maru is not produced by instant atmospheric forcing but is mainly produced by entrainment on the right side of the running typhoon as a part of cooling area during its passage, and remained during a few days. The sea surface cooling (SSC) is evident along the track and on the right side of the running typhoon, which is similar to the SSC of satellite observation by TRMM/TMI. The conspicuous SSC produced by both entrainment and upwelling is situated just under the track of typhoon when the typhoon moves slower. Intercomparison of entrainment schemes of the mixed layer model is implemented. Frictional velocity and buoyancy effects are effective for a gradual SSC covering the wide region. In contrast, the effect of current shear at the mixed layer base is related to the amount of SSC and the sharp horizontal gradient of SSC. The entrainment scheme including all three TKE sources has the best performance for SSC simulation.  相似文献   

20.
For settlement of the well-known problem of contemporary radar imaging models, i. e. , the problem of a general underestimation of radar signatures of hydrodynamic features over oceanic internal waves and underwater bottom topography in tidal waters at high radar frequency bands ( X-band and C-band), the impact of the ocean surface mixed layer turbulence and the significance of strat- ified oceanic model on SAR remote sensing of internal solitary waves are proposed. In the north of the South China Sea by utilizing some observed data of background field the nonlinearity coefficient, the dispersion coefficient, the horizontal variability coefficient and the phase speed in the generalized K-dV equation are determined approximately. Through simulations of internal tide transfor- mation the temporal evolution and spatial distribution of the vertical displacement and horizontal velocity of internal wave field are obtained. The simulation results indicate that the maximum amplitudes of internal solitary waves occur at depth 35 m, but the maximum current speeds take place at depth 20 m in this area of the sea (about 20°30'N, 114°E) in August. It was noticed that considering the effects of flood current and ebb current respectively is appropriate to investigate influence of the background shear flow on coefficients of the K-dV equation. The obtained results provide the possibility for the simulation of SAR signatures of internal solitary waves under considering the impact of ocean surface mixed layer turbulence in the companion paper.  相似文献   

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

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