首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
于子棚  刘海龙  林鹏飞 《大气科学》2017,41(5):1087-1100
海洋中的潮汐混合对大西洋经圈翻转环流AMOC(Atlantic Meridional Overturning Circulation)模拟的影响是海洋环流模式研究的热点问题之一。本文采用IAP/LASG发展的气候系统海洋模式LICOM(LASG/IAP Climate system Ocean Model)及与海冰耦合模式进行了有无潮汐混合方案的试验,重点探讨了潮汐混合对AMOC强度模拟的影响。结果显示,引入潮汐混合后模拟的AMOC强度极大值比对照试验增加约1倍,更接近RAPID(Rapid Climate Change Programme)观测。而且,潮汐混合试验中模拟的AMOC上层环流深度(3200 m)比对照试验加深1000 m左右,同样更接近RAPID观测。海洋底部的垂直混合增强,使海洋层结变得更加不稳定,加强了北大西洋高纬地区,特别是拉布拉多海等地区的深对流,这是AMOC加强的直接原因。同时,潮汐混合试验中上层海洋环流也加强,增加了中低纬副热带高盐海水向高纬输送,使表层增密,海洋层结更加不稳定,也可以进一步增强AMOC。  相似文献   

2.
为了定量评估北京气候中心(BCC)发展的BCC_CSM对当代全球海表温度和混合层深度的模拟能力,以WOA09(World Ocean Atlas 2009)观测资料作为检验模式的气候态实况场,提取包括BCC_CSM在内的CMIP5中的17个海气耦合模式的模拟结果,评估BCC_CSM模拟的全球海表温度和混合层深度的气候平均态并分析造成偏差的可能原因。结果表明:BCC_CSM模拟的海表温度在北半球中高纬的误差较大,而在其余纬度的模拟性能较佳。偏差的产生主要归因于海洋环流偏差。BCC_CSM模拟的最深混合层在北半球中高纬和南半球高纬地区的误差较大,同时这些区域也是多模式模拟差异最大的区域;其模拟的最浅混合层在南半球中高纬的偏差较大。冬季大西洋经向翻转环流的模拟在北大西洋下沉的位置偏南导致北半球高纬地区海表温度偏冷。由此认为包括BCC_CSM在内的许多海气耦合模式需重点改进对南、北半球深对流海域物理过程的描述,以提高气候预测的可信度。  相似文献   

3.
全球海洋模式对不同强迫场的响应   总被引:1,自引:0,他引:1  
使用中国科学院大气物理研究所研制的全球海洋环流模式(LASG/IAP Climate system Ocean Model,LICOM),通过设计三个试验,即以德国马克斯—普朗克气象研究所整理的海洋模式比较计划(OMIP)资料和美国国家海洋资料中心(NODC)发布的《世界海洋图集2009》(WOA09)资料为强迫场的试验W,用美国环境预报中心(NCEP)和国家大气研究中心(NCAR)联合推出的NCEP/NCAR再分析资料(简称NCEP资料)中的风应力资料代替试验W中的风应力资料的试验M,以及用NCEP资料中的热力强迫代替试验M中的热力强迫资料的试验N,来研究不同的热力和动力强迫场对模式的影响。三个试验的模拟结果均模拟出了水团和流场的分布型和极值区。从三个试验的结果对比可以看出,NCEP资料较弱的风应力使得试验M环流场明显偏弱,减弱了大洋内部的温盐输送,加大了深海温盐模拟结果与观测资料的偏差,但对原模式过强的南极中层水的输送有所改善。NCEP的短波辐射通量和非短波热通量弱于OMIP,且在两极区域NCEP资料的海表温度比WOA09资料最多低 4℃以上。试验N的模拟结果改善了南大洋60°S以南海区试验W 模拟的海表温度偏高问题,减小了北冰洋部分海域以及副热带大洋东部海表温度的偏差。此外,试验N高纬度较低的海表温度增强了北大西洋深水以及南极底层水的输送,因而改善了深海的温盐模拟结果。三个试验在一些关键海区得到的经向热输送在观测估计及前人模拟结果的范围中,总体上试验M的输送最弱。综合三个试验的模拟结果,可以认为OMIP风应力资料和NCEP海表温度资料更适合作为LICOM模式的强迫场。  相似文献   

4.
表层洋流对外强迫响应敏感度的数值研究   总被引:2,自引:1,他引:1  
利用数值模拟研究了海表流场对外强迫(风应力和海表热通量)的响应特征,探讨了其对该类外强迫异常响应的敏感性以及较敏感区域。在确认本文所用的海洋环流模式能够较好地模拟表层海洋流场的气候状态之后,通过几个敏感性试验与控制(对照)试验结果的比较,发现海洋表层环流对海表风应力异常响应的敏感区域主要在赤道附近及大洋西边界海区;相对于热带外地区,热带海域(20°S~20°N)的风应力异常对于大洋表层环流的变化有着更重要的显著作用,它不仅会导致热带海域表层流场有较大的变化,对中高纬海区的表层流场特别是西边界流也有明显影响;海洋表层环流对海表热通量异常的响应除了在赤道附近海域明显之外,在中高纬海区也十分显著;在外强迫有同等异常幅度(20%)的情况下,大洋西边界海域对热通量的响应明显要强于对风应力的响应。此外,热通量异常还对南太平洋东海岸的洋流和南极大陆的绕极环流有较为明显的影响。  相似文献   

5.
对比两个同化资料GODAS(Global Ocean Data Assimilation System)和SODA(Simple Ocean Data Assimilation),考察中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的气候系统海洋模式LICOM(LASG/IAP Climate system Ocean Model)模拟的北印度洋经向环流及热输送的气候态。LICOM能抓住北印度洋大尺度环流的季节变化特征,模拟的年平均越赤道热输送为-0.24 PW (1 PW=1015W),较之以往的数值模式结果更接近观测和同化资料。与同化资料的差异主要体现在季节变化强度,北半球夏季在赤道以南偏弱0.5 PW,这与模式夏季的纬向风应力偏弱,热输送中的大项Ekman热输送模拟偏弱,从而模拟的经圈翻转环流较浅有关。  相似文献   

6.
金向泽  张学洪 《大气科学》1994,18(Z1):769-779
本文是用简单海一气耦合模型模拟温盐环流在全球增暖事件中作用的研究工作的第一部分。为了建立一个简单海一气耦合模型,我们首先根据Wright和Stoker等人的设计复制出一个包括大西洋、太平洋和南大洋在内的二维温盐环流模式,从等温、等盐和无运动的初始状态出发,在给定的年平均海表强迫下将模式积分了4000年,模拟出了和原作相似的温盐环流。对模拟结果的分析表明,相对于北太平洋而言,北大西洋北部的高盐、低温特点(后者是由两大洋在地理上的差别决定的)是形成当代温盐环流的主要原因;从与温盐环流相联系的海表热通量来看,北大西洋北部是向大气提供热量的主要源地;模式温盐环流对于海表盐度通量的敏感性试验的结果表明,对于纬圈平均的二维模式而言,要想模拟出合理的温盐环流就必须人为地提高北大西洋北部的海表盐度,文章分析了这种作法的物理根据;模式中的对流过程对于温盐环流的维持是至关重要的,对比有无季节循环的试验结果可以看出,虽然温度场的明显的季节变化只出现在模式的最上面两层,但由于引进季节循环后冬季高纬海洋的对流活动加强,后者直接影响到温盐环流,使更多的深海热量上传并向大气释放。这是使海洋温跃层得以保持合理.厚度的一个重要原因。  相似文献   

7.
王万秋 《大气科学》1993,17(5):555-562
本文分析了用IAP两层大气环流模式模拟的海表动量和热量通量,并将其同Han等和Esbensen等的气候资料比较.模拟的热量通量与观测估计值有类似的水平分布和季节变化,但在中低纬地区有偏多的热量由海洋向大气输送,尤其1月北半球中纬大气需要从海洋获得过分多的热量;模拟的海表动量通量和气候估计值也类似,但模拟的北半球冬季中高纬西风动量通量中心位置偏东,赤道中西太平洋和大西洋的东风动量通量偏弱,南半球环绕南极的西风带模拟得过分弱(尤其在7月份).本文还检验了基本变量的日变化和日际变化对计算海表动量和湍流扩散热量通量的影响,结果表明用月平均,日平均和每小时的基本量计算的动量和湍流扩散热量通量依次增大,尤其在中高纬地区更明显.  相似文献   

8.
海洋环流模式模拟自然 和核辐射14C的分布   总被引:9,自引:1,他引:8  
金心  石广玉 《大气科学》2000,24(3):341-354
放射性14C在海洋环流研究和人为CO2问题的研究中都有重要地位。本文用海洋环流模式模拟了海洋中自然14C的分布及海洋对核辐射产生的放射性14C的吸收, 以期对海洋吸收人为CO2的能力做一初步的研究。模拟的海洋环流结果与观测相比符合得较好,成功地模拟出了北大西洋深水(NADW)、南极底水 (AABW)等基本特征。对自然14C的模拟揭示出了海洋通风的基本特征。模拟出的沿GEOSECS 路径的南、北垂直截面与观测结果符合得较好。对核辐射14C的模拟表明:模式模拟的沿GEOSECS 路径的南、北垂直截面与观测结果符合得较好;模拟出的海洋表面核辐射浓度与观测值一致,但核辐射14C在海洋中的柱存量和平均穿透深度都比观测结果要小。文中分析了造成这种差异的可能原因。  相似文献   

9.
参照Griffies et al.(2009)提出的海洋—海冰耦合模式参考试验(Coordinated Ocean-ice Reference Experiments,COREs),设计了一个800年积分的数值试验,对一个质量严格守恒的压力坐标海洋环流模式(Pressure Coordinate Ocean Model,PCOM1.0)的基本模拟性能进行了评估,并与观测资料和再分析资料进行了对比。结果表明,PCOM1.0模拟的温盐场和基本流场与COREs模式的模拟水平基本接近。其中,模拟的大西洋经向翻转流在45°N附近达到18 Sv(1 Sv=106 m3 s-1),与观测估计值接近;对海表面温度的模拟误差主要集中在北太平洋黑潮区和北大西洋湾流区等中高纬度急流区;模拟的热带太平洋温跃层过于深厚;模拟的经德雷克海峡的体积输送达130 Sv,比大部分COREs模式及再分析资料都更接近于观测估计值。  相似文献   

10.
评估了中国科学院大气物理研究所大气科学和地球流体动力学数值模拟国家重点实验室海洋环流模式L30T63和海气耦合模式FGCM-0模拟的热带太平洋年平均状态,资料取自L30T63由观测的大气强迫驱动的Control试验、由NCAR CCM3大气强迫驱动的Spinup试验、以及相应的海气耦合模式FGCM-0.主要的结论是:(1)在"准确"的海表强迫下,Control模拟的海面温度和温跃层与观测结果相当接近,模式的固有误差是赤道冷舌过分西伸和东南太平洋温跃层偏浅.(2)Spinup能模拟出合理的热带太平洋上层海洋环流,但存在两个问题,即:暖池区海面温度显著偏高、沿赤道的梯度过大;赤道温跃层偏浅、东西向坡度偏小,它们分别与CCM3提供的海表短波辐射通量和风应力的系统误差有关.这两个问题很可能是海气耦合模式FGCM-0运行初期误差迅速发展的重要原因.(3)FGCM-0模拟的赤道暖池区上层100 m的平均温度比观测低3℃.分析表明FGCM-0夸大了暖池区海洋动力过程的降温作用,使得模拟的"暖池"在一定程度上具有冷舌的属性.FGCM-0模拟的热带南太平洋温跃层比观测结果偏浅数十米到100 m,以致赤道两侧的上层海洋温度分布趋于对称,成为"double ITCZ"现象在上层海洋中的表现.风应力旋度的系统误差和垂直混合随深度衰减过快是温跃层偏浅的两个可能原因;FGCM-0中与北太平洋中高纬地区深厚冷偏差相关的经圈环流也有利于热带温跃层误差的维持.  相似文献   

11.
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.  相似文献   

12.
全球海洋CFC-11吸收对传输速度的敏感性   总被引:4,自引:2,他引:2  
赵琦  陈中笑  徐永福  李阳春 《大气科学》2012,36(6):1253-1268
本文在中国科学院大气物理研究所发展的全球海洋模式(LICOM)中使用五个不同的海气交换的气体 传输速度公式对CFC-11(一氟三氯甲烷)在海洋中的分布和吸收做了模拟。讨论了不同气体传输速度的差异, 发现差异最大的两个公式得到的全球年平均传输速度相差81%。对CFC-11的海表浓度分布、海气通量、水柱总量、海水累积吸收量以及在大洋断面的垂直浓度分布进行了分析。分析结果显示, 使用Liss and Merlivat (1986) 的传输速度公式的试验在海气通量和海洋吸收总量的模拟上均小于其他试验, Nightingale et al. (2000)、Ho et al. (2006) 和Sweeney et al. (2007) 等的公式虽然全球年平均值相近, 但在高风速地区Nightingale et al. (2000) 公式的传输速度要小于后两者, 导致了使用该公式的试验模拟结果在主吸收区和存储区的强度比后二者偏小。Wanninkhof (1992) 的公式在形式上与Ho et al. (2006) 以及Sweeney et al. (2007) 的公式一致, 但在系数上存在差别, 这使得模拟的水柱总量在南大洋的分布明显好于其他试验, 尽管其最大值仍比观测资料略小。在海洋累积吸收量的计算上, 使用Wanninkhof (1992) 传输速度公式得到的模拟结果比观测资料小8%左右。计算了Liss and Merlivat (1986) 和Wanninkhof (1992) 的传输速度公式模拟的单年吸收量相对差, 其总体上一直保持持续下降的趋势, 到2007年仅为2%。从该相对差变化趋势看, 在最初的前10年, 海气CFC-11交换通量对海气交换传输速度的敏感性更强, 而在更长时间的模拟上, 海洋对CFC-11的吸收则更依赖于物理模式的通风速率。通过对CFC-11垂直断面分布的分析可知, 不同的传输速度在主要吸收区的不同导致了一定的垂直分布差异。基于本文的结果可以认为Wanninkhof (1992) 的海气气体传输速度公式更适合本模式对CFC-11的模拟。  相似文献   

13.
The regional ocean circulation in the coastal areas or China (including a part of the western Pacific Ocean, the South China Sea and the Bay of Bengal el al.) is simulated by using the improved Princeton Uni-versity ocean circulation model (POM). Compared with the modeling results obtained by the large-scale ocean general circulation model (OGCM), the basic ocean circulation features simulated by the regional ocean circulation model are in good agreement with that simulated by OGCM and some detailed character-istics such as the regional ocean circulation, sea temperature, salinity and free sea surface height have also been obtained which are in good accord with the observations. These results indicate that the regional ocean circulation model has good capability to produce the regional ocean circulation characteristics and it can be used to develop coupled regional ocean-atmospheric model systems.  相似文献   

14.
An ocean biogeochemistry model was developed and incorporated into a global ocean general circulation model (LICOM) to form an ocean biogeochemistry general circulation model (OBGCM). The model was used to study the natural carbon cycle and the uptake and storage of anthropogenic CO2 in the ocean. A global export production of 12.5 Pg C yr-1 was obtained. The model estimated that in the pre-industrial era the global equatorial region within 15o of the equator released 0.97 Pg C yr-1 to the atmosphere, which was balanced by the gain of CO2 in other regions. The post-industrial air-sea CO2 flux indicated the oceanic uptake of CO2 emitted by human activities. An increase of 20-50 mol kg-1 for surface dissolved inorganic carbon (DIC) concentrations in the 1990s relative to pre-industrial times was obtained in the simulation, which was consistent with data-based estimates. The model generated a total anthropogenic carbon inventory of 105 Pg C as of 1994, which was within the range of estimates by other researchers. Various transports of both natural and anthropogenic DIC as well as labile dissolved organic carbon (LDOC) were estimated from the simulation. It was realized that the Southern Ocean and the high-latitude region of the North Pacific are important export regions where accumulative air-sea CO2 fluxes are larger than the DIC inventory, whereas the subtropical regions are acceptance regions. The interhemispheric transport of total natural carbon (DIC+LDOC) was found to be northward (0.11 Pg C yr-1), which was just balanced by the gain of carbon from the atmosphere in the Southern Hemisphere.  相似文献   

15.
基于ECMWF的ERA-40海表10m风场,对1958—2001年全球海表风速的变化趋势进行分析,主要分析了整体变化趋势、变化趋势的季节性差异、区域性差异、变化周期。结果表明:①近44年期间,全球海域海表风速整体上以0.0067m·s-1·a-1的速度显著性逐年线性递增。1958—1975年全球海域的海表风速变化较为平缓,1975—1983年递增趋势较为强劲,年平均海表风速的峰值出现在1999年,波谷出现在1975年。②全球海表风速的变化趋势表现出较大的区域性差异。递增趋势明显的区域主要分布于:南极、热带大西洋海域、北太平洋西风带海域、印度洋中低纬度海域、南半球60°S附近大面积带状海域;呈显著性逐年递减的区域主要分布于:赤道中东太平洋、胡安·费尔南德斯群岛附近海域、南大西洋西风带的中部海域,以及一些零星海域。③全球海表风速的变化趋势表现出较大的季节性差异。在各月均表现出显著的线性递增趋势,以1月的递增趋势最为强劲,达到0.0103m·s-1·a-1,7月的递增趋势弱于其余月份,约0.0033m·s-1·a-1。④全球海域海表风速存在明显的2.2~4.3年变化周期,以及6.5年以上长周期震荡。  相似文献   

16.
1.IntroductionInrecentyears,manystudiesonthemechanismoftheregionalclimateformationandvariationhavebeendonebyusingnestedhigh--resolutionregionalclimatemodels.TheregionalclimatemodelshavepresentedbetterPerformancesinsimulatingregionalclimatefeaturesthanlarge--scalegeneralcirculationmodels(GCM)becauseoftheaccuraterepresentationsofhigh--resolutiontopography,detailedunderlyingsurfacecharacteristics,landsurfaceprocessesandplanetaryboundarylayerparameterization.However,theoceanpartwithinthemodeldom…  相似文献   

17.
华莉娟  俞永强  尹宝树 《大气科学》2010,34(6):1046-1058
热带印度洋偶极子 (Indian Ocean Dipole) 是印度洋海域内海洋和大气环流年际变化的主要特征模态之一, 在热带海气耦合系统中起到非常重要的作用。同热带太平洋的ENSO现象类似, 热带印度洋偶极子也呈现出显著的不对称性。本文利用中国科学院大气物理研究所发展的全球海洋环流模式, 在观测风应力距平的强迫下, 评估了模式对热带印度洋季节变化、 热带印度洋偶极子 (IOD) 模态及其不对称性的模拟能力, 并且通过数值试验分析了IOD模态不对称性特征及其对气候平均态的影响。对照观测资料, 模式较好地再现了热带印度洋SST在季风驱动下的季节变化特征。在年际时间尺度上, 模式不仅能够再现IOD指数的变化趋势, 而且可以成功模拟出IOD模态的空间分布特征, 即表层和次表层海温在西印度洋表现为正异常, 在东印度洋表现为负异常。可见, 对于热带印度洋而言, IOD模态主要是对风应力异常的响应。热带印度洋海温与Niño3.4指数的相关性分析表明, 模式能够模拟出超前热带太平洋ENSO现象2~4个月时海温的偶极子型分布, 但是不能模拟出滞后ENSO现象2个月左右的全海盆增暖模态, 可能是因为模式试验中没有考虑热通量年际异常的强迫。同时, 模式模拟的IOD模态具有同观测结果相类似的不对称性, 进一步的敏感性试验表明风应力的不对称性对偶极子指数的不对称性贡献较小, 次表层及以下海温的不对称性可能主要受到海洋内部非线性动力过程的影响。通过数值试验, 本文还发现热带印度洋海温的不对称性对气候平均态会有影响, 而这种不对称性长期积累后, 会导致上层热带印度洋温度层结趋于稳定状态。  相似文献   

18.
Summary As a precursor to the embedding of a thermal stratification model of the ocean into an Ocean General Circulation Model (OGCM) we have compared the performance of an integral (or bulk) (Kraus-Turner) and a differential (Eddy Diffusion Dimension 1) thermal stratification model for the eastern Indian Ocean. Forced with the Comprehensive Ocean-Atmosphere Data Set (COADS) and Subjectively Estimated Cloudiness in the Australian Region (SECAR) data sets, both models successfully stimulated the annual cycle in sea surface temperature. On the other hand, their predictions of the mixed layer depth differed substantially. Whereas the differential model slightly underestimated the winter mixed layer depth, the winter mixed layer depths calculated by the bulk model were substantially larger than the observed values. An additional analysis of the simulated heat storage cycle of the eastern Indian Ocean with respect to the activity of tropical-extratropical cloudbands in the Australian region has reveaved a strong relationship between these two parameters. The results have shown that a stronger activity of the tropical-extratropical cloudbands coincides with periods of heat transfer oriented from the ocean to the atmosphere.With 11 Figures  相似文献   

19.
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  相似文献   

20.
1986—1987厄尔尼诺事件的数值模拟   总被引:1,自引:0,他引:1  
张荣华 《大气科学》1994,18(Z1):847-855
用高分辨率自由表面热带太平洋环流模式,在观测到的风应力和热量、水汽通量驱动下,对1986—1987厄尔尼诺(E1Nino)事件进行了数值模拟。各种变量场的时空结构及其演变表明,模式成功地模拟出1986—1987厄尔尼诺现象。始于1986年年中,赤道西太平洋的西风异常所推动的向东表层洋流不断向中、东太平洋输送暖水,至11月份,大量暖水在日界线附近堆积,造成海面上升(达32cm)和斜温层(用20℃等温线深度表示)加深。1986年年底的强西风异常激发出赤道Kelvin波,并向赤道东太平洋和南美沿岸传播,使那里的斜温层加深和海面上升,且具有双峰结构;Kelvin波所伴随的垂直冷平流的减弱造成赤道中、东太平洋海表温度上升;1987年春季在中、东太平洋和南美沿岸地区存在强的正海表温度异常,并伴随着整个赤道太平洋斜温层东西方向变平、赤道潜流弱而中心位置变浅。厄尔尼诺相伴随的热带太平洋环流异常首先于1987年年中从东太平洋开始消失,而中、西太平洋则一直维持到1988年初。  相似文献   

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

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