首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 250 毫秒
1.
The Community Atmosphere Model version 3 (CAM3) temperature simulation bias is examined in this paper. We compare CAM3 output with European Centre for Medium-Range Weather Forecasts (ECMWF) 40 year reanalysis (ERA-40) data. We formulate a time mean temperature bias equation then evaluate each term in the equation. Our focus is on the Northern Hemisphere winter time. We group the temperature equation terms into these categories: linear advection terms, nonlinear advection terms, transient eddy terms and diabatic heating, and find that linear advection and diabatic bias are the largest. The nonlinear terms (velocity bias advection of temperature bias) are much smaller than each of the other groups of terms at all levels except near the surface. Linear advection terms have dipolar pattern in the Atlantic (negative NW of positive) which reflects the shift of the CAM3 model North Atlantic storm track (NAST) into Europe, especially in the upper troposphere; opposite sign dipolar structure occurs over Alaska (positive) and the north Pacific storm track (negative). The transient advection terms in middle latitudes are larger in the upper troposphere and generally positive along the Atlantic storm track. Along the north Pacific storm track (NPST), the transient terms are negative in the mid and lower troposphere over much of the NPST (positive in upper troposphere). The diabatic heating bias has large values in the tropics along the Intertropical Convergence Zone (ICZ) and along the midlatitude storm tracks. During this time of year the ICZ is mainly in the Southern Hemisphere, but CAM3 emphasizes an ICZ-like heating in the northern hemisphere of the Atlantic and Pacific Oceans. CAM3 tends to have a weaker ICZ, especially in the Atlantic. In midlatitudes, we find large bias in heating by precipitation and vertically averaged net radiation over the NAST, Europe, and the Middle East.  相似文献   

2.
This study builds upon two prior papers, which examine Arctic region bias of CAM3 (NCAR Community Atmosphere Model version 3) simulations during winter. CAM3 output is compared with ECMWF (European Centre for Medium-Range Weather Forecasts) 40?year reanalysis (ERA-40) data. Our prior papers considered the temperature and the vorticity equation terms and demonstrated that diabatic, transient, and linear terms dominate nonlinear bias terms over most areas of interest. Accordingly, this paper uses a linearized form of the model??s dynamical core equations to study aspects of the forcing that lead to the CAM3 biases. We treat the model??s long term winter bias as a solution to a linear stationary wave model (LSWM). Key features of the bias in the vorticity, temperature, and ln of surface pressure (=q) fields are shown at medium resolution. The important features found at medium resolution are captured at the much lower LSWM resolution. The Arctic q bias has two key features: excess q over the Barents Sea and a missing Beaufort High (negative maximum q bias) to the north of Alaska and eastern Siberia. The forcing fields are calculated by the LSWM. Horizontal advection tends to create multi-polar combinations of negative and positive extrema in the forcing. The positive and negative areas of forcing approximately match corresponding areas in the bias. There is a broad relation between cold bias with elevated q bias, as expected from classical theory. Forcing in related quantities: near surface vorticity and surface pressure combine to produce the sea level pressure bias.  相似文献   

3.
Streamwise vorticity equation   总被引:5,自引:0,他引:5  
l.IntroductionTheverticalcomPOnentofthevorticityisthedominantpartinthelarge-scalesystems,therefore,itisgoodenoughtoonlyconsidertheverticalcomPOnentofthevorticityinthevorticityequationforstudyinglarge-scalesystems.Withtheapproximationthatstreamlinesaremainlyquasi-horizontalincontrasttotheirverticalcomponents,hence,thecomPOnentofabsolutevorticityalongthestreamline,so-calledsecondaryvorticity,isverysmallandnearlyomitted.ButinthemesoscaleandsmaIl-scalesystems,especiallyinsmall-scalesys-tems,ver…  相似文献   

4.
利用1979~2013年NCEP再分析数据,通过经验正交分解对比了前冬时期北大西洋风暴轴的高低空分布,并用涡动动能(Eddy Kinetic Energy,EKE)方程对风暴轴高低空分布型差异进行了诊断。研究结果表明:上层和下层第一空间分布型差异巨大,对流层下层风暴轴中心偏北,靠近极地,而上层风暴轴中心偏西南,靠近北美沿岸。EKE方程诊断结果表明:正压转换项在高低空符号相反,导致了EKE在上、下层分布出现显著差异,即上层正压转换项为负,在扰动发展中起能量耗散作用,而下层正压转换项为正,且极大值区域对应下层EKE极大值区域,为风暴轴下层向极区域增强的主要原因。而斜压转换和非地转位势通量散度在上层均为正,且远大于下层,为风暴轴上层涡动能量维持的原因,也从涡动能量收支上解释了风暴轴的主体出现在上层。  相似文献   

5.
西伯利亚风暴轴的气候特征及其可能维持机制   总被引:2,自引:1,他引:1  
马阳  朱伟军  李天宇  王森  李欣 《气象科学》2017,37(5):587-597
基于1959—2014年NCEP/NCAR的逐日再分析资料,首先研究了西伯利亚风暴轴各季节的气候平均特征,然后以冬季为例,利用能量诊断方程,从能量学的角度对其的可能维持机制进行了探讨,并在上述分析过程中与北半球两大洋风暴轴的特征进行了对比。结果表明:(1)西伯利亚风暴轴一年四季都独立存在,虽强度要比两大洋风暴轴的强度弱很多且位置偏北,但可以定义为一个弱风暴轴。(2)比较来看,西伯利亚风暴轴强度的季节变化与北太平洋风暴轴的季节变化类似。与两大洋风暴轴位于急流东北侧不同,冬季西伯利亚风暴轴位于东亚温带急流的西侧。(3)进一步的能量分析结果表明,与两大洋风暴轴一样,斜压不稳定的能量转换(Ke4)也是西伯利亚风暴轴区域天气尺度扰动动能的主要来源;而扰动非地转位势通量散度项(Ke3)和时间平均气流对扰动动能的平流输送项(Ke1)也是风暴轴下游发展所需的扰动动能来源之一。  相似文献   

6.
With 40 years integration output of two atmospheric general circulation models (GAMIL/IAP and HadAM3/UKMO) forced with identical prescribed seasonally-varying sea surface temperature, this study examines the effect of the observed Indian-western Pacific Ocean (IWP) warming on the Northern Hemisphere storm tracks. Both models indicate that the observed IWP warming tends to cause both the North Pacific storm track (NPST) and the North Atlantic storm track (NAST) to move northward. Such a consistent effect on the two storm tracks is closely associated with the changes in the low-level atmospheric baroclinicity, high-level jet stream and upper-level geopotential height. The IWP warming can excite a wavelike circum-global teleconnection in the geopotential height that gives rise to an anticyclonic anomaly over the midlatitude North Pacific and a positive-phase NAO anomaly over the North Atlantic. These geopotential height anomalies tend to enhance upper-level zonal westerly winds north of the climatological jet axes and increase low-level baroclinicity and eddy growth rates, thus favoring transient eddy more active north of the climatological storm track axes, responsible for the northward shift of the both storm tracks. The IWP warming-induced northward shift of the NAST is quite similar to the observed, suggesting that the IWP warming can be one of the key factors to cause decadal northward shift of the NAST since the 1980s. However, the IWP warming-induced northward shift of the NPST is completely opposite to the observed, implying that the observed southward shift of the NPST since the 1980s would be primarily attributed to other reasons, although the IWP warming can have a cancelling effect against those reasons.  相似文献   

7.
The complete form of the vertical vorticity tendency equation (the complete-form vorticity equation) is derived from the Ertel potential vorticity equation to contain thermodynamic factors. In this study, a new complete-form vorticity equation, which has the same form as the original complete-form vorticity equation, is deduced from the absolute vorticity vector equation combined with the continuity equation and the expression of three-dimensional (3D) entropy gradient. By comparing the complete-form vorticity equation with the classical vertical vorticity equation, it is found that regardless of whether or not the isentropic surface is tilting, the two vorticity equations are in essence the same. The “baroclinic term” of the complete-form vorticity equation is exactly equal to the solenoidal term of the classical one, and there is a significant amount of cancellation between the two baroclinic items (the “slantwise term” and the horizontal vorticity change term) in the complete-form vorticity equation. In operational weather analysis, the tilt of the isentropic surface can be diagnosed according to the density of the isotherm on the upper-level isobaric map. For synoptic-scale motion, the vertical vorticity produced by the tilt of the isentropic surface is due to the contribution of atmospheric baroclinicity, which is measured by the solenoid. The 3D solenoid is parallel to the isentropic surface, so the more tilted the isentropic surface, the bigger the projection of the 3D solenoid in the vertical direction. The baroclinic contribution can be interpreted based on the PV thinking theory, but the relationship between the vorticity field and the potential vorticity field is not immediate.  相似文献   

8.
李智  刘宣飞  李传浩 《大气科学》2015,39(6):1081-1094
采用美国国家环境预测中心的CFSR(Climate Forecast System Reanalysis)再分析资料和QuickSCAT(Quick Scatterometer)、AVHRR(Advanced Very High Resolution Radiometer)、TRMM(Tropical Rainfall Measuring Mission)高分辨率卫星资料,研究了大气对春季东海黑潮锋响应的气压调整机制及其年际变化。结果表明,春季东海黑潮锋位于黑潮暖舌的西北侧,呈西南-东北走向,与大尺度气压背景场的等压线走向一致,锋区东南侧暖水与西北侧冷水之间产生的局地气压梯度与大尺度气压梯度形成同向叠加,使得锋区附近西北指向东南的气压梯度达到最大,造成该处的海表面10 m矢量风速也最大,在摩擦作用下形成东北偏北风(NNE)。锋区与其东南侧的NNE风之间沿锋区走向(跨锋区走向)的分量差,会在暖舌附近产生气旋性切变涡度(风速辐合),由此产生上升运动和强降水;而在锋区西北侧的冷水区情况正好相反,有反气旋性切变涡度(风速辐散),并伴有下沉运动和弱降水,从而形成跨锋区的次级环流圈。东海黑潮锋区偏强(弱)年,锋区东南侧暖水与西北侧冷水之间的局地气压梯度也偏强(弱),与大尺度气压梯度同向叠加后形成偏强(弱)的NNE风,造成锋区东南侧暖舌附近的气旋性切变涡度、风速辐合、上升运动和降水均偏强(弱),而锋区西北侧冷水区的反气旋性切变涡度、风速辐散和下沉运动均偏强(弱),跨锋区次级环流圈偏强(弱),这表明在年际时间尺度上气压调整机制仍起作用。  相似文献   

9.
利用NCEP/NCAR再分析资料,运用31点带通数字滤波、线性相关和合成分析方法,研究了1961/1962—2010/2011年冬季北太平洋风暴轴西部、东部区域强度指数的年际演变特征,划分了风暴轴的典型型态,并进一步探讨了与同期北半球500 hPa位势高度场和SLP的关系。结果表明:风暴轴气候态的极大值区域位于中纬度北太平洋中西部,最大值点的频数集中区域和均方差分布的异常中心都有两个。风暴轴西部和东部区域强度指数(WI和EI)的年际演变具有独立性,典型型态可分为单、双中心型两类。WI(EI)指数与北半球500 hPa位势高度场的相关分布类似于WP(PNA)遥相关型;单中心型风暴轴偏强时,极涡南扩,平均槽加深;呈双中心型时,极涡明显偏西。WI(EI)指数与SLP的相关分布类似于NPO(NAO)遥相关型;单中心型风暴轴偏强(弱)时,SLP距平场呈AO遥相关型的正(负)异常位相。  相似文献   

10.
A high resolution regional atmosphere model is used to investigate the sensitivity of the North Atlantic storm track to the spatial and temporal resolution of the sea surface temperature (SST) data used as a lower boundary condition. The model is run over an unusually large domain covering all of the North Atlantic and Europe, and is shown to produce a very good simulation of the observed storm track structure. The model is forced at the lateral boundaries with 15–20 years of data from the ERA-40 reanalysis, and at the lower boundary by SST data of differing resolution. The impacts of increasing spatial and temporal resolution are assessed separately, and in both cases increasing the resolution leads to subtle, but significant changes in the storm track. In some, but not all cases these changes act to reduce the small storm track biases seen in the model when it is forced with low-resolution SSTs. In addition there are several clear mesoscale responses to increased spatial SST resolution, with surface heat fluxes and convective precipitation increasing by 10–20% along the Gulf Stream SST gradient.  相似文献   

11.
分析 2 0 0 1年 7月 2 - 3日台风“榴莲”暴雨过程中湿位涡及其各分量的变化 ,发现对流层低层 85 0 h Pa湿位涡的负值中心、 70 0 h Pa湿位涡的正值区与强降水中心相对应 ;急流与层结稳定度的变化 ,影响着湿位涡的变化  相似文献   

12.
Spatiotemporal patterns of recent (1979–2008) air temperature trends are evaluated using three reanalysis datasets and radiosonde data. Our analysis demonstrates large discrepancies between the reanalysis datasets, possibly due to differences in the data assimilation procedures as well as sparseness and inhomogeneity of high-latitude observations. We test the robustness of arctic tropospheric warming based on the ERA-40 dataset. ERA-40 Arctic atmosphere temperatures tend to be closer to the observed ones in terms of root mean square error compared to other reanalysis products used in the article. However, changes in the ERA-40 data assimilation procedure produce unphysical jumps in atmospheric temperatures, which may be the likely reason for the elevated tropospheric warming trend in 1979–2002. NCEP/NCAR Reanalysis data show that the near-surface upward temperature trend over the same period is greater than the tropospheric trend, which is consistent with direct radiosonde observations and inconsistent with ERA-40 results. A change of sign in the winter temperature trend from negative to positive in the late 1980s is documented in the upper troposphere/lower stratosphere with a maximum over the Canadian Arctic, based on radiosonde data. This change from cooling to warming tendency is associated with weakening of the stratospheric polar vortex and shift of its center toward the Siberian coast and possibly can be explained by the changes in the dynamics of the Arctic Oscillation. This temporal pattern is consistent with multi-decadal variations of key arctic climate parameters like, for example, surface air temperature and oceanic freshwater content. Elucidating the mechanisms behind these changes will be critical to understanding the complex nature of high-latitude variability and its impact on global climate change.  相似文献   

13.
利用多普勒天气雷达资料及反演风场和常规观测资料,对2014年11月2日发生在北黄海(山东半岛北部海上)一次罕见海龙卷风暴的中尺度特征进行了分析。结果表明:冷空气、暖湿海面热力边界、山东半岛北部近海岸西北风与偏西风的辐合线是海龙卷风暴发生的天气背景。海龙卷风暴发生时雷达回波PPI最大分贝反射率因子为60 dBZ,高度为2.0 km,最高风暴顶为4.5 km,最大垂直累积液态水含量VIL为21 kg·m-2。利用雷达反演风场进行中尺度特征分析,结果表明:在海龙卷风暴发生发展过程中,低层风辐合对应4.0 km高度上是风辐散,海上有较强的偏南暖湿气流输送到雷暴区。中尺度动力特征:最大正涡度和散度辐合在1.0 km以下,低层正涡度和散度辐合、高层散度辐散是雷暴发生初期动力特征;低层没有正涡度和散度辐合、高层为正涡度和散度辐合是雷暴开始发展的动力特征;低层和高层为大的正涡度和散度辐合是雷暴成熟阶段的动力特征。高空冷空气叠加上低空强的偏南气流,造成局地涡度加大和低层辐合加强,使低层暖湿气流倾斜上升。海龙卷与辐合区的冷空气和暖湿气流有关。  相似文献   

14.
1. Introduction There have been two di?erent approaches used fortracking extratropical cyclones. The traditional andmost common approach is to follow the minimum sur-face pressure of a cyclone (e.g., Petterssen, 1956; Car-nell and Senior, 1998; Serreze…  相似文献   

15.
利用实况资料和WRF中尺度数值模式对2008年6月12日18时—14日00时的华南双雨带暴雨过程进行了数值模拟与诊断分析。结果表明:随着锋面的南压,在锋面的西南方向(广西沿海)生成一低涡,该低涡作为位涡源在中高层表现稳定,分别为锋面雨带(北雨带)与暖区雨带(南雨带)提供正位涡。南雨带对北雨带的作用主要体现在中层(112~114°E附近),南雨带中有位涡的大值向北输送,其输送过程导致两条雨带在该处相连,而在115°E以东的南雨带则无明显的输送过程。同时,北部高空槽中也有大值位涡向北雨带输送,以维持北雨带。研究还发现,本次过程中暖区暴雨与锋面暴雨雨带的结构差异明显,锋面雨带的结构与传统雨带的结构比较一致;有利于暖区暴雨降水的形势主要表现在中高层。RIP轨迹模式的结果也表明,质点在运动过程中位涡的输送源是位于广西沿海的低涡,可见该位涡源对双雨带形成有重要的作用。  相似文献   

16.
北太平洋风暴轴“深冬抑制”现象的能量分析   总被引:2,自引:2,他引:0  
采用欧洲中期天气预报中心逐日再分析资料(ERA-40),从局地能量变化方程出发,通过分析北太平洋风暴轴区域对流层不同层次局地能量的季节演变过程,对风暴轴区域各能量项在“深冬抑制”现象中的作用进行了深入探讨。结果表明,天气尺度扰动动能的季节变化可以很好地反映北太平洋风暴轴的“深冬抑制”现象,并且该现象在对流层上层最为显著,其发生概率约为80%,其中20世纪70年代中后期到80年代前期抑制最强。从同期各能量项的变化来看,扰动动能的变化主要受斜压能量转换项、涡动非地转位势通量的散度项和正压能量转换项的影响。在深冬季节,由于消耗扰动动能的正压能量转换项虽有些微弱减少从而使得扰动动能有所增加,但为风暴轴提供扰动动能的斜压能量转换项和涡动非地转位势通量的散度项减少的幅度却更大,因而总的效果是扰动动能大为减小,这可能是造成北太平洋风暴轴“深冬抑制”现象的直接原因。  相似文献   

17.
北太平洋风暴轴的气候特征及其变化的初步研究   总被引:7,自引:7,他引:7  
丁叶风  任雪娟  韩博 《气象科学》2006,26(3):237-243
本文使用ECMWF再分析网格点资料(ERA-40),分析了不同季节和不同高度层上北太平洋风暴轴(天气尺度瞬变扰动活动)的气候特征及其时间演变规律。分析表明,气候平均而言,北太平洋风暴轴冬季强且偏西南,夏季弱且偏东北。在近45 a里,冬季及夏季风暴轴在整个对流层、尤其是对流层中上层具有整体一致的年际和年代际变化特征。在各个高度层上冬季风暴轴于1985年前后一致性地发生了由弱至强的年代际跃变;夏季风暴轴发生年代际由弱至强跃变的时间在各个高度层上略有不同,对流层中高层的跃变时间与北太平洋大气海洋系统1970年代的跃变时间一致。在冬季风暴轴活动偏强年里,中纬度天气尺度瞬变扰动增强的同时,位置也有所北抬,反之亦然。夏季风暴轴活动偏强(弱)年则主要表现为瞬变扰动在气候平均位置上的增强(减弱)。  相似文献   

18.
北太平洋东部风暴轴的时空演变特征   总被引:6,自引:2,他引:4  
本文利用欧洲中期天气预报中心(ECMWF)逐日再分析资料(ERA-40),以500 hPa位势高度滤波方差为代表,对1957年12月~2001年11月期间44年528个月北太平洋区域(30°N~60°N,120°E~120°W)月平均风暴轴的多中心数目和最强中心位置进行了客观统计,在此基础上,对北太平洋区域进行了分区,通过对比各区域风暴轴的时间演变和结构变化,重点揭示了北太平洋东部地区风暴轴的时空演变特征.主要结论如下:(1)逐月来看,北太平洋风暴轴“多中心”现象普遍存在,概率高达94.7%,最典型的分布呈2~3个中心分布;从季节上来看,春季是“多中心”现象最容易出现的季节,秋季和冬季相对较少,而4个及以上的“多中心”现象则更容易出现在夏季.(2)若把最强中心出现在160°W以东地区的北太平洋风暴轴定义为东部型风暴轴,那么从月份上来看,7月相对最容易出现东部型风暴轴,1月和2月最难;从季节上来看,夏季相对最容易出现东部型风暴轴,冬季最难;总的来看,出现东部型风暴轴的频数大约占总频数的三分之一.(3)从垂直结构上看,在北太平洋160°W以东地区,风暴轴的强度可以最强,但与斜压性密切联系的涡动向极和向上热量通量的最大值却并不是最强.(4)经验正交函数分解(EOF)分析的结果表明,在不同季节、不同区域以及是否单独考虑东部型风暴轴的情况下,风暴轴的变化虽然表现出了一定的差异,但都反映出在北太平洋东部区域风暴轴的变化特征有其独特特点,如在该区域风暴轴的主要变化模态并不一定时时与其他区域的主要变化模态一一对应.北太平洋东部区域风暴轴变化的原因和机制值得进一步深入探讨.  相似文献   

19.
运用NCEP、Had ISST再分析资料,北大西洋涛动(NAO)月指数序列,探讨了海表面温度(SST)锋的时空变化特征,揭示了北大西洋SST锋的主要气候变率及其与北大西洋风暴轴和大气大尺度环流异常的关系。研究表明,剔除季节循环后的SST锋显示其最主要变率为锋区的向南/北摆动,其对应的风暴轴发生相应的西南/东北移动,并同时在北大西洋上空对应一个跨海盆的位势高度负/正异常。这种环流异常可引起高纬度海平面气压(SLP)的反气旋/气旋式环流,这有利于增强海表面风对大洋副极地环流的负/正涡度异常输入,进一步减弱/加强了高纬度上层冷水向SST锋区的输送。北大西洋SST锋的另一主要模态为锋区在南北方向的分支和合并。当SST锋异常在40°N~45°N以单支形式加强时,对流层位势高度场和SLP南北梯度增大,对应NAO正位相,此时风暴轴也为单支型;同时SLP异常场促使冰岛附近具有气旋式风应力异常,亚速尔地区具有反气旋式风应力异常,导致副极地环流和副热带环流均加强,增加高纬度冷水和低纬度暖水在锋区的输入,从而进一步增强40°N~45°N附近的SST锋区。当SST锋异常在40°N~45°N纬带南北发生分支时,风暴轴也同时出现北强南弱的南北分支,此时对应了负位相NAO,来自北南的冷暖水输送减弱,SST锋也发生减弱分支。此外,位于大洋内区的SST锋东端也存在一个偶极子型的模态,尽管其解释方差相对较小,但仍与偏东北的NAO型具有显著相关。谱分析表明,北大西洋SST锋与风暴轴具有1~3年和年代际共振,与中高纬大尺度环流也存在周期1~3年的共变信号,其中准一年共变信号体现了SST锋和NAO之间的对应关系。进一步诊断分析表明,SST锋上空的近表层大气斜压性和经向温度梯度随着SST锋的增强而增强,经向热通量的向北输送导致涡动有效位能的增加;海洋的非绝热加热产生更强的垂直热量通量,这有利于涡动有效位能释放成为涡动动能,从而表现为该区域的风暴轴加强,并进一步影响风暴轴中的天气尺度扰动与下游大尺度环流异常的相互作用过程。  相似文献   

20.
利用欧洲中期天气预报中心ECMWF(European Center for Medium-range Weather Forecast)逐日再分析资料(ERA40),通过经验正交函数(empirical orthogonal function,EOF)分解发现,冬季北太平洋东部风暴轴有着显著的年际变化特征:第一变化模态为在气候平均位置南北相反的偶极子变化型,第二变化模态为在气候平均位置处一致增强或减弱的变化型,第三变化模态为三极子的变化型。进一步的回归分析发现:当东部风暴轴南压(北抬)时,同期冬季是一种厄尔尼诺(拉尼娜)年海温异常空间分布型,中纬度北太平洋海区以及赤道中、东海区,冬季冷(暖)异常的洋面上是异常低压(高压),对流层中层是太平洋—北美型(Pacific-North American Pattern,PNA)遥相关的正(负)位相;当东部风暴轴增强(减弱)时,同期冬季黑潮区海温偏暖(偏冷),对流层中层表现为西太平洋型(West Pacific Pattern,WP)遥相关的正(负)位相;当东部风暴轴呈现西北—东南+-+(-+-)相间三极子的分布时,同期冬季巴布亚新几内亚附近海温异常偏暖(冷),夏威夷附近海温异常偏冷(暖),冬季冷(暖)异常的洋面上是异常低(高)压,对流层中层表现出类似PNA正(负)位相。EOF分解各模态所对应时间系数与阿留申低压(Aleutian Low,AL)指数、PNA指数、Nino3指数、WP指数、黑潮海温(Kuroshio Current,KC)指数之间存在显著的相关,这些证明了东部风暴轴与同期大气环流及SST异常之间的联系。  相似文献   

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

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