首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO2 in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO2 to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO2 uptake to be 1.64 and 1.73 Pg C yr-1 for the 1990s, and that the global ocean contained 86.8 and 92.7 Pg C of anthropogenic CO2 at the end of 1994, respectively. Both the total inventory and uptake from our model are smaller than the data-based estimates. In this presentation, the vertical distributions of anthropogenic CO2 at three meridional sections are discussed and compared with the available data-based estimates. The inventory in the individual basins is also calculated. Use of large isopycnal diffusivity can generally improve the simulated results, including the exchange flux, the vertical distribution patterns, inventory, storage, etc. In terms of comparison of the vertical distributions and column inventory, we find that the total inventory in the Pacific Ocean obtained from our model is in good agreement with the data-based estimate, but a large difference exists in the Atlantic Ocean, particularly in the South Atlantic. The main reasons are weak vertical mixing and that our model generates small exchange fluxes of anthropogenic CO2 in the Southern Ocean. Improvement in the simulation of the vertical transport and sea ice in the Southern Ocean is important in future work.  相似文献   

2.
A stability analysis of the coupled ocean–atmosphere is presented which shows that the potential energy (PE) of the upper layer of the ocean is available to generate coupled growing planetary waves. An independent analysis suggests that the growth of these waves would be maintained in the presence of oceanic friction. The growing waves are a consequence of relaxing the rigid lid approximation on the ocean, thus allowing an upward transfer of energy across the sea surface. Using a two and a half layer model consisting of an atmospheric planetary boundary layer, coupled with a two layer ocean comprising an active upper layer and a lower layer in which the velocity perturbation is vanishingly small, it is shown that coupled unstable waves are generated, which extract PE from the main thermocline. The instability analysis is an extension of earlier work [Tellus 44A (1992) 67], which considered the coupled instability of an atmospheric planetary boundary layer coupled with an oceanic mixed layer, in which unstable waves were generated which extract PE from the seasonal thermocline. The unstable wave is an atmospheric divergent barotropic Rossby wave, which is steered by the zonal wind velocity, and has a wavelength of about 6000 km, and propagates eastward at the speed of the deep ocean current. It is argued that this instability, which has a multidecadal growth time constant, may be generated in the Southern Ocean, and that its properties are similar to observations of the Antarctic Circumpolar Wave (ACW).  相似文献   

3.
Sea ice is an important component in the Earth’s climate system. Coupled climate system models are indispensable tools for the study of sea ice, its internal processes, interaction with other components, and projection of future changes. This paper evaluates the simulation of sea ice by the Flexible Global Ocean-Atmosphere-Land System model Grid-point Version 2 (FGOALS-g2), in the fifth phase of the Coupled Model Inter-comparison Project (CMIP5), with a focus on historical experiments and late 20th century simulation. Through analysis, we find that FGOALS-g2 produces reasonable Arctic and Antarctic sea ice climatology and variability. Sea ice spatial distribution and seasonal change characteristics are well captured. The decrease of Arctic sea ice extent in the late 20th century is reproduced in simulations, although the decrease trend is lower compared with observations. Simulated Antarctic sea ice shows a reasonable distribution and seasonal cycle with high accordance to the amplitude of winter-summer changes. Large improvement is achieved as compared with FGOALS-g1.0 in CMIP3. Diagnosis of atmospheric and oceanic forcing on sea ice reveals several shortcomings and major aspects to improve upon in the future: (1) ocean model improvements to remove the artificial island at the North Pole; (2) higher resolution of the atmosphere model for better simulation of important features such as, among others, the Icelandic Low and westerly wind over the Southern Ocean; and (3) ocean model improvements to accurately receive freshwater input from land, and higher resolution for resolving major water channels in the Canadian Arctic Archipelago.  相似文献   

4.
A global, flux-corrected climate model is employed to predict the surface wind stress and associated wind-driven oceanic circulation for climate states corresponding to a doubling and quadrupling of the atmospheric CO2 concentration in a simple 1% per year CO2 increase scenario. The model indicates that in response to CO2 increase, the position of zero wind stress curl in the mid-latitudes of the Southern Hemisphere shifts poleward. In addition, the wind stress intensifies significantly in the mid-latitudes of the Southern Hemisphere. As a result, the rate of water circulation in the subpolar meridional overturning cell in the Southern Ocean increases by about 6 Sv (1 Sv=106 m3 s−1) for doubled CO2 and by 12 Sv for quadrupled CO2, implying an increase of deep water upwelling south of the circumpolar flow and an increase of Ekman pumping north of it. In addition, the changes in the wind stress and wind stress curl translate into changes in the horizontal mass transport, leading to a poleward expansion of the subtropical gyres in both hemispheres, and to strengthening of the Antarctic Circumpolar Current. Finally, the intensified near-surface winds over the Southern Ocean result in a substantial increase of mechanical energy supply to the ocean general circulation.  相似文献   

5.
刘波  李阳春  徐永福  范广洲 《大气科学》2015,39(6):1149-1164
本文选用中国科学院大气物理研究所全球海洋模式(LICOM),对中尺度涡旋参数化方案(GM90方案)中等密度扩散系数和等密度面厚度扩散系数(统称为涡旋扩散系数Aρ)对物理场及CFC-11(一氟三氯甲烷)分布的影响进行了研究。本文做了两个试验,即涡旋扩散系数采用常系数方式(对照试验)和采用在非绝热层以下Aρ随海洋浮力频率垂直变化的参数化方案(浮力试验)。模拟结果表明,依浮力频率垂直变化的方案对模式物理场的模拟能力有一定程度的提升,如南极绕极流的输送强度比常系数方案增大了约20%~30%,与观测事实更接近;浮力试验对对照试验中过强的南极中层水有一定的削弱作用,使得模式对南大洋高纬次表层位密度的模拟有一定的改善。稍有不足的是,浮力试验对南极底层水也有一定的削弱,使得2000~3000 m深度位密度模拟较常系数方案偏低。通过对CFC-11分布、存储以及输送的研究发现,次网格参数取值的不同对南大洋CFC-11模拟情况有较大影响。浮力试验加大了南北高纬CFC-11海表的吸收通量,对南极大陆周边海域向南大洋主储藏区(34°S~60°S)的CFC-11输送能力有一定的增强,使得南大洋对CFC-11储藏量增大,大部分海区与观测资料更接近。通过CFC-11断面分析,浮力试验对南大洋上层海洋位密度模拟的改善使得CFC-11垂直结构与观测更接近。  相似文献   

6.
We evaluate the representation of dynamic sea surface height (SSH) fields of 33 global coupled models (GCMs) contributed to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). We use observations from satellite altimetry and basic performance metrics to quantify the ability of the GCMs to replicate observed SSH of the time-mean, seasonal cycle, and inter-annual variability patterns. The time-mean SSH representation has markedly improved from CMIP3 to CMIP5, both in terms of overall reduction in root-mean square differences, and in terms of reduced GCM ensemble spread. Biases of the time-mean SSH field in the Indian and Pacific Ocean equatorial regions are consistent with biases in the zonal surface wind stress fields identified with independent measurements. In the Southern Ocean, the latitude of the maximum meridional gradient of the zonal mean SSH CMIP5 models is shifted equatorward, consistent with the GCMs’ spatial biases in the maximum of the zonal mean westerly surface wind stress fields. However, while the Southern Ocean SSH gradients correlate well with the maximum Antarctic circumpolar current transports, there is no significant correlation with the maximum zonal mean wind stress amplitudes, consistent with recent findings that the eddy parameterisations in GCMs dominate over wind stress amplitudes in this region. There is considerable spread across the CMIP5 ensemble for the seasonal and interannual SSH variability patterns. Because of the short observational period, the interannual variability patterns depend on the time-period over which they are derived, while no such dependency is found for the time-mean patterns. The model performance metrics for SSH presented here provide insight into GCM shortcoming due to inadequate model physics or processes. While the diagnostics of CMIP5 GCM performance relative to observations reveal that some models are clearly better than others, model performance is sensitive to the spatio-temporal scales chosen.  相似文献   

7.
 Recent observational and numerical studies of the maritime snow cover in the Antarctic suggest that snow on top of sea ice plays a major role in shaping the seasonal growth and decay of the ice pack in the Southern Ocean. Here, we make a quantitative assessment of the importance of snow accumulation in controlling the seasonal cycle of the ice cover with a coupled snow–sea-ice–upper-ocean model. The model takes into account snow and ice sublimation and snow deposition by condensation. A parametrisation of the formation of snow ice (ice resulting from the freezing of a mixture of snow and seawater produced by flooding of the ice floes) is also included. Experiments on the sensitivity of the snow–sea-ice system to variations in the sublimation/condensation rate, the precipitation rate, and the amount of snowfall transported by the wind into leads are discussed. Although we focus on the model response in the Southern Hemisphere, results for the Arctic are also discussed in some cases to highlight the relative importance of the processes under study in both hemispheres. It is found that the snow loss by sublimation can account for the removal of 0.45 m of snow per year in the Antarctic and that this loss significantly affects the total volume of snow ice. A precipitation decrease of 50% is conducive to large reductions in the Antarctic snow and snow-ice volumes, but it leads only to an 8% decrease in the annual mean ice volume. The Southern Ocean ice pack is more sensitive to increases in precipitation. For precipitation rates 1.5 times larger than the control ones, the annual mean snow, ice, and snow-ice volumes augment by 30, 20, and 180%, respectively. It is also found that the transfer to the ocean of as much as 50% of the precipitating snow as a result of wind transport has almost negligible effects on the total ice volume. All the experiments exhibit a marked geographical contrast in the ice-cover response, with a much larger sensitivity in the western sector of the Southern Ocean than in the eastern sector. Our results suggest that snow-related processes are of secondary importance for determining the sensitivity of the Arctic sea ice to environmental changes but that these processes could have an important part to play in the response of the Antarctic sea-ice cover to future, or current, climatic changes. Received: 30 June 1997/Accepted: 2 October 1998  相似文献   

8.
A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO2 in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO2 to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO2 uptake to be 1.64 and 1.73 ...  相似文献   

9.
The Kuroshio Extension region is characterized by energetic oceanic mesoscale and frontal variability that alters the air–sea fluxes that can influence large-scale climate variability in the North Pacific. We investigate this mesoscale air-sea coupling using a regional eddy-resolving coupled ocean–atmosphere (OA) model that downscales the observed large-scale climate variability from 2001 to 2007. The model simulates many aspects of the observed seasonal cycle of OA coupling strength for both momentum and turbulent heat fluxes. We introduce a new modeling approach to study the scale-dependence of two well-known mechanisms for the surface wind response to mesoscale sea surface temperatures (SSTs), namely, the ‘vertical mixing mechanism’ (VMM) and the ‘pressure adjustment mechanism’ (PAM). We compare the fully coupled model to the same model with an online, 2-D spatial smoother applied to remove the mesoscale SST field felt by the atmosphere. Both VMM and PAM are found to be active during the strong wintertime peak seen in the coupling strength in both the model and observations. For VMM, large-scale SST gradients surprisingly generate coupling between downwind SST gradient and wind stress divergence that is often stronger than the coupling on the mesoscale, indicating their joint importance in OA interaction in this region. In contrast, VMM coupling between crosswind SST gradient and wind stress curl occurs only on the mesoscale, and not over large-scale SST gradients, indicating the essential role of the ocean mesocale. For PAM, the model results indicate that coupling between the Laplacian of sea level pressure and surface wind convergence occurs for both mesoscale and large-scale processes, but inclusion of the mesoscale roughly doubles the coupling strength. Coupling between latent heat flux and SST is found to be significant throughout the entire seasonal cycle in both fully coupled mode and large-scale coupled mode, with peak coupling during winter months. The atmospheric response to the oceanic mesoscale SST is also studied by comparing the fully coupled run to an uncoupled atmospheric model forced with smoothed SST prescribed from the coupled run. Precipitation anomalies are found to be forced by surface wind convergence patterns that are driven by mesoscale SST gradients, indicating the importance of the ocean forcing the atmosphere at this scale.  相似文献   

10.
邢如楠 《气象学报》1986,44(1):10-17
为研究大尺度海气相互作用、气候以及气候的变化,一个具有规则海岸线和无海底地形的三层斜压海洋模式已经建立起来。积分区域包括了赤道地区在内,南北方向从-2750 km到+4750 km,东西方向宽10000 km,总深度是4 km。模式方程为原始方程,用了静力近似和Bou-sinesq近似。模式数值积分分两个阶段,在第一个阶段,海洋环流从静止状态开始,以只随纬度变化的年平均风应力和海表热通量作为强迫边条件数值积分十年,这时海洋上层的环流和温度都趋于平稳,我们将第十年末的积分状态作为模式气候的准平衡态。计算结果模拟了热带太平洋上气候平均的温度分布以及主要海流的大尺度特征,如赤道地区的狭长冷水带和强上翻区、南赤道洋流、北赤道逆流以及表层下面沿赤道自西向东的潜流。在第二个阶段,以在第一个阶段末尾得到的准平衡态作为初值,分别用冬季和夏季的风应力作为强迫边条件再积分一年。比较这两种计算结果,我们看到海洋环流随信风系统的季节变化也呈现出明显的季节变化,特别是北赤道逆流,其强度在夏季强,在冬季弱甚至不出现。南赤道洋流无论哪个季节都穿过赤道,在赤道以北的南赤道洋流冬季比夏季强,而在赤道以南的南赤道洋流夏季比冬季强。  相似文献   

11.
Experiments with the coupled climate model CLIMBER-3α, which contains an oceanic general circulation model, show deep upwelling in the Southern Ocean to be proportional to the surface wind stress in the latitudinal band of Drake Passage. At the same time, the distribution of the Southern Ocean upwelling onto the oceanic basins is controlled by buoyancy distribution; the inflow into each basin being proportional to the respective meridional density difference. We observe approximately the same constant of proportionality for all basins, and demonstrate that it can be directly related to the flow geometry. For increased wind stress in the Southern Ocean, the overturning increases both in the Atlantic and the Indo-Pacific basin. For strongly reduced wind stress, the circulation enters a regime where Atlantic overturning is maintained through Pacific upwelling, in order to satisfy the transports set by the density differences. Previous results on surface buoyancy and wind stress forcing, obtained with different models, are reproduced within one model in order to distill a consistent picture. We propose that both Southern Ocean upwelling and meridional density differences set up a system of conditions that determine the global meridional overturning circulation.  相似文献   

12.
The ocean response to surface temperature transients is simulated with the use of the Hamburg large-scale geostrophic (LSG) ocean general circulation model (OGCM). The transition, from the present to a climate corresponding to a doubling of the atmospheric CO2 content, is compared with the reversed transition. For the Atlantic, the time scale for the deep ocean to adjust to the temperature changes was similar for both transitions. In the Pacific, the time scale is shorter for the present to warm transition than for the reverse case, a result of increased production of Antarctic bottom water (AABW) during the warm climate. While the transition from cold to warm climate shows no secular variability, the reversed transition generates considerable variability on time scales of 300–400 years. For the warm climate, oscillations with periods of 45 years are found in the Southern Ocean. Results of principal oscillation pattern (POP) analysis indicate that these oscillations are due to interaction between convection in the Southern Ocean and advected salinity anomalies in the Antarctic Circumpolar Current (ACC) and the Southern Pacific Ocean. Received: 19 September 1995 / Accepted: 15 March 1996  相似文献   

13.
This paper examines the mean annual cycle, interannual variability, and leading patterns of the tropical Atlantic Ocean simulated in a long-term integration of the climate forecast system (CFS), a state-of-the-art coupled general circulation model presently used for operational climate prediction at the National Centers for Environmental Prediction. By comparing the CFS simulation with corresponding observation-based analyses or reanalyses, it is shown that the CFS captures the seasonal mean climate, including the zonal gradients of sea surface temperature (SST) in the equatorial Atlantic Ocean, even though the CFS produces warm mean biases and underestimates the variability over the southeastern ocean. The seasonal transition from warm to cold phase along the equator is delayed 1 month in the CFS compared with the observations. This delay might be related to the failure of the model to simulate the cross-equatorial meridional wind associated with the African monsoon. The CFS also realistically simulates both the spatial structure and spectral distributions of the three major leading patterns of the SST anomalies in the tropical Atlantic Ocean: the south tropical Atlantic pattern (STA), the North tropical Atlantic pattern (NTA), and the southern subtropical Atlantic pattern (SSA). The CFS simulates the seasonal dependence of these patterns and partially reproduces their association with the El Niño-Southern Oscillation. The dynamical and thermodynamical processes associated with these patterns in the simulation and the observations are similar. The air-sea interaction processes associated with the STA pattern are well simulated in the CFS. The primary feature of the anomalous circulation in the Northern Hemisphere (NH) associated with the NTA pattern resembles that in the Southern Hemisphere (SH) linked with the SSA pattern, implying a similarity of the mechanisms in the evolution of these patterns and their connection with the tropical and extratropical anomalies in their respective hemispheres. The anomalies associated with both the SSA and NTA patterns are dominated by atmospheric fluctuations of equivalent-barotropic structure in the extratropics including zonally symmetric and asymmetric components. The zonally symmetric variability is associated with the annular modes, the Arctic Oscillation in the NH and the Antarctic Oscillation in the SH. The zonally asymmetric part of the anomalies in the Atlantic is teleconnected with the anomalies over the tropical Pacific. The misplaced teleconnection center over the southern subtropical ocean may be one of the reasons for the deformation of the SSA pattern in the CFS.  相似文献   

14.
重点分析了20世纪80年代以来国内外关于大气角动量平衡研究的一些代表性工作。总结了地气系统角动量交换、角动量输送及其与山脉和摩擦力矩异常、地转异常、ENSO等关系的研究结果。指出了一些应当深入研究的方向。  相似文献   

15.
A key question in studies of the potential for reducing uncertainty in climate change projections is how additional observations may be used to constrain models. We examine the case of ocean carbon cycle models. The reliability of ocean models in projecting oceanic CO2 uptake is fundamentally dependent on their skills in simulating ocean circulation and air–sea gas exchange. In this study we demonstrate how a model simulation of multiple tracers and utilization of a variety of observational data help us to obtain additional information about the parameterization of ocean circulation and air–sea gas exchange, relative to approaches that use only a single tracer. The benefit of using multiple tracers is based on the fact that individual tracer holds unique information with regard to ocean mixing, circulation, and air–sea gas exchange. In a previous modeling study, we have shown that the simulation of radiocarbon enables us to identify the importance of parameterizing sub-grid scale ocean mixing processes in terms of diffusive mixing along constant density surface (isopycnal mixing) and the inclusion of the effect of mesoscale eddies. In this study we show that the simulation of phosphate, a major macronutrient in the ocean, helps us to detect a weak isopycnal mixing in the upper ocean that does not show up in the radiocarbon simulation. We also show that the simulation of chlorofluorocarbons (CFCs) reveals excessive upwelling in the Southern Ocean, which is also not apparent in radiocarbon simulations. Furthermore, the updated ocean inventory data of man-made radiocarbon produced by nuclear tests (bomb 14C) enable us to recalibrate the rate of air–sea gas exchange. The progressive modifications made in the model based on the simulation of additional tracers and utilization of updated observational data overall improve the model’s ability to simulate ocean circulation and air–sea gas exchange, particularly in the Southern Ocean, and has great consequence for projected CO2 uptake. Simulated global ocean uptake of anthropogenic CO2 from pre-industrial time to the present day by both previous and updated models are within the range of observational-based estimates, but with substantial regional difference, especially in the Southern Ocean. By year 2100, the updated model estimated CO2 uptake are 531 and 133 PgC (1PgC?=?1015 gram carbon) for the global and Southern Ocean respectively, whereas the previous version model estimated values are 540 and 190 PgC.  相似文献   

16.
The forcing efficiency for the first and the second baroclinic modes by the wind stress in tropical oceans has been discussed by calculating equivalent forcing depth from annual mean, seasonal, and pentadal density profiles of the observational data. In the annual mean field, the first mode is forced preferentially in the western Pacific and the Indian Ocean, whereas the second mode is more strongly excited in the Atlantic and the eastern Pacific. This difference is mostly due to the pycnocline depth; the second mode is more dominantly forced where the pycnocline depth is shallower. We also revealed large seasonal variations of the second mode's equivalent forcing depth in the western Indian Ocean. The first mode is more dominantly forced during boreal spring and fall in the western Indian Ocean, while the second mode becomes more dominantly forced during boreal summer and winter. Those are due to seasonal variations of both the zonal wind and the pycnocline depth. Moreover, we show that the excitation of the second mode in the western Pacific increases after the late 1970s, which is associated with the decreasing trend of the zonal pycnocline gradient. Revealing the variation of the equivalent forcing depth will be useful for understanding the oceanic response to winds in tropical oceans and the improvement in the predictability of air-sea coupled climate variability in the tropics.  相似文献   

17.
Impacts on the coupled variability of the Indo-Pacific by including the effects of surface currents on surface stress are explored in four extended integrations of an experimental version of the Bureau of Meteorology??s coupled seasonal forecast model POAMA. The first pair of simulations differs only in their treatment of momentum coupling: one version includes the effects of surface currents on the surface stress computation and the other does not. The version that includes the effect of surface currents has less mean-state bias in the equatorial Pacific cold tongue but produces relatively weak coupled variability in the Tropics, especially that related to the Indian Ocean dipole (IOD) and El Ni?o/Southern Oscillation (ENSO). The version without the effects of surface currents has greater bias in the Pacific cold tongue but stronger IOD and ENSO variability. In order to diagnose the role of changes in local coupling from changes in remote forcing by ENSO for causing changes in IOD variability, a second set of simulations is conducted where effects of surface currents are included only in the Indian Ocean and only in the Pacific Ocean. IOD variability is found to be equally reduced by inclusion of the local effects of surface currents in the Indian Ocean and by the reduction of ENSO variability as a result of including effects of surface currents in the Pacific. Some implications of these results for predictability of the IOD and its dependence on ENSO, and for ocean subsurface data assimilation are discussed.  相似文献   

18.
J. Egger 《Climate Dynamics》2001,18(1-2):169-177
 A data set comprising the equatorial components of global atmospheric momentum and the related mountain and friction torques is used to demonstrate the usefulness of a novel analysis technique in climate research. The technique consists in deriving a Master equation in the phase plane of a pair of climate variables on the basis of time series of this pair. The Master equation predicts changes of the pair's probability distribution in this phase plane. A mean velocity and a diffusion coefficient can be derived from the coefficients of this equation. An analysis of the dynamics of the parameter pair is based on the Master equation. The circulation of the angular momenta in the phase plane is anticyclonic with an attraction towards the origin. The diffusion is strongest close to the center of the probability density distribution of the equatorial momenta. The anticyclonic motion reflects the westward propagation of those Rossby waves which contribute to the mass term. It is found that the mountain and orographic torques induce only a small fraction of the observed changes of the angular momenta. It is essentially the impact of both the Earth's deviation from sphericity and its rotation, which induces the variability of the angular momenta. The mass and wind terms are barely related to each other in contrast to the predictions of normal mode theory. Surprisingly, however, the mountain torque and the wind terms are closely linked. Received: 13 April 2000 / Accepted: 22 February 2001  相似文献   

19.
ResponseProcessofOceantoAtmosphericForcingandOptimalResponseFrequencyintheCZOceanModelNiYunqi(倪允琪),ZouLi(邹力)andWuAiming(吴爱明)(...  相似文献   

20.
王亚非  魏东  李琰 《高原气象》2011,30(5):1189-1194
利用NCEP/NCAR再分析资料计算了1998年5~6月区域角动量收支的主要分量,对比分析了范围较大的欧亚区域角动量矩和限制在欧亚大陆范围内的喜马拉雅山角动量矩的差异,发现两者有很好的相关,限制较小范围的计算表明,喜马拉雅山山脉力矩的作用显著。角动量收支项占较大量级的山脉力矩和摩擦力矩,并与梅雨锋区域的700hPa高度...  相似文献   

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

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