首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 687 毫秒
1.
A western North Pacific tropical cyclone (TC) intensity prediction scheme (WIPS) is developed based on TC samples from 1996 to 2002 using the stepwise regression technique, with the western North Pacific divided into three sub-regions: the region near the coast of East China (ECR), the South China Sea region (SCR), and the far oceanic region (FOR). Only the TCs with maximum sustained surface wind speed greater than 17.2 m s-1 are used in the scheme. Potential predictors include the climatology and persistence factors, synoptic environmental conditions, potential intensity of a TC and proximity of a TC to land. Variances explained by the selected predictors suggest that the potential intensity of a TC and the proximity of a TC to land are significant in almost all the forecast equations. Other important predictors include vertical wind shear in ECR, 500-hPa geopotential height anomaly at the TC center, zonal component of TC translation speed in SCR, intensity change of TC 12 or 24 h prior to initial time, and the longitude of TC center in FOR. Independent tests are carried out for TCs in 4 yr (2004-2007), with mean absolute errors of the maximum surface wind being 3.0, 5.0, 6.5, 7.3, 7.6, and 7.9 m s-1 for 12- to 72-h predictions at 12-h intervals, respectively. Positive skills are obtained at all leading time levels as compared to the climatology and persistence prediction scheme, and the large skill scores (near or over 20%) after 36 h imply that WIPS performs especially better at longer leading times. Furthermore, it is found that the amendment in TC track prediction and real-time model analysis can significantly improve the performance of WIPS in the SCR and ECR. Future improvements will focus on applying the scheme for weakening TCs and those near the coastal regions.  相似文献   

2.
Comparison of Three Tropical Cyclone Intensity Datasets   总被引:3,自引:0,他引:3  
Analyzed in this paper are the 16-yr (1988-2003) tropical cyclone (TC) intensity data from three major forecast centers of the western North Pacific, i.e., China Meteorological Administration (CMA), Regional Specialized Meteorological Center Tokyo (RSMC Tokyo), and Joint Typhoon Warning Center (JTWC) of the United States. Results show that there are significant discrepancies (at 1% significance level) in the intensity of TCs among the three centers, with a maximum difference for the same TC over 30 m s-1. The flight reconnaissance over TC can minish the discrepancy to some extent. A climatic and persistent prediction model is set up to study the impact of initial data from different forecast centers on the prediction of TC intensity. It is obtained that the root mean square error (RMSE) of a 4-yr independent test is the largest using data from JTWC, while the smallest using data from RSMC Tokyo. Average absolute deviation in 24-h intensity prediction is 2.5 m s-1 between CMA and RSMC Tokyo data, and 4.0 m s-1 between CMA and JTWC data, with a maximum deviation reaching 21 m s-1. Such a problem in the initial value increases the difficulty in intensity prediction of TCs over the western North Pacific.  相似文献   

3.
A western North Pacific tropical cyclone (TC) intensity prediction scheme has been developed based on climatology and persistence (CLIPER) factors as potential predictors and using genetic neural network (GNN) model. TC samples during June–October spanning 2001–2010 are used for model development. The GNN model input is constructed from potential predictors by employing both a stepwise regression method and an Isometric Mapping (Isomap) algorithm. The Isomap algorithm is capable of finding meaningful low-dimensional architectures hidden in their nonlinear high-dimensional data space and separating the underlying factors. In this scheme, the new developed model, which is termed the GNN-Isomap model, is used for monthly TC intensity prediction at 24- and 48-h lead times. Using identical modeling samples and independent samples, predictions of the GNN-Isomap model are compared with the widely used CLIPER method. By adopting different numbers of nearest neighbors, results of sensitivity experiments show that the mean absolute prediction errors of the independent samples using GNN-Isomap model at 24- and 48-h forecasts are smaller than those using CLIPER method. Positive skills are obtained as compared to the CLIPER method with being above 12 % at 24 h and above 14 % at 48 h. Analyses of the new scheme suggest that the useful linear and nonlinear prediction information of the full pool of potential predictors is excavated in terms of the stepwise regression method and the Isomap algorithm. Moreover, the GNN is built by integrating multiple individual neural networks with the same expected output and network architecture is optimized by an evolutionary genetic algorithm, so the generalization capacity of the GNN-Isomap model is significantly enhanced, indicating a potentially better operational weather prediction.  相似文献   

4.
Using tropical cyclone (TC) best track and intensity of the western North Pacific data from the Joint TyphoonWarning Center (JTWC) of the United States and the NCEP/NCAR reanalysis data for the period of 1992-2002, the effects of vertical wind shear on TC intensity are examined. The samples were limited to the westward or northwestward moving TCs between 5°N and 20°N in order to minimize thermodynamic effects. It is found that the effect of vertical wind shear between 200 and 500 hPa on TC intensity change is larger than that of the shear between 500 and 850 hPa, while similar to that of the shear between 200 and 850 hPa. Vertical wind shear may have a threshold value, which tends to decrease as TC intensifies. As the intensifying rate of TC weakens, the average shear increases. The large shear has the obvious trend of inhibiting TC development. The average shear of TC which can develop into typhoon (tropical depression or tropical storm) is below 7 m s-1 (above 8 m s-1).  相似文献   

5.
2016年GRAPES_TYM改进及对台风预报影响   总被引:1,自引:0,他引:1       下载免费PDF全文
为了进一步提高国家气象中心区域模式台风数值预报系统(GRAPES_TYM)的预报能力,2016年对模式参考大气廓线以及涡旋初始化方案进行了改进:由模式初始场水平方向平均的一维参考大气代替原来的等温大气,涡旋初始化方案取消了原涡旋重定位并将涡旋强度调整半径由原来的12°减小到4°。对2014—2016年的生命史超过3 d的所有台风进行了回算,路径及近地面最大风速统计误差分析表明:参考大气的改进可以减小模式对台风预报路径预报的系统北偏和平均路径误差,尤其是140°E以东的转向台风。涡旋初始化方案中强度调整半径的减小会进一步减小模式预报路径的北偏趋势,从而进一步减小平均误差。同业务系统预报结果相比,改进后的GRAPES_TYM(包括参考大气和涡旋初始化)可以使平均路径误差分别减小10%(24 h),12%(48 h),16%(72 h),14%(96 h)以及15%(120 h)。同美国NCEP全球模式路径预报相比,GRAPES_TYM在西行、西北行登陆我国的台风路径预报有一定优势。  相似文献   

6.
The effects of vertical wind shear on tropical cyclone (TC) intensity change are examined based on the TC data from the China Meteorological Administration and the NCEP reanalysis daily data from 2001 to 2006. First, the influence of wind shear between different vertical levels and averages in different horizontal areas are compared. The results indicate that the effect of wind shear between 200 and 850 hPa averaged within a 200–800 km annulus on TC intensity change is larger than any other calculated vertical wind shear. High-latitude and intense TCs tend to be less sensitive to the effects of VWS than low-latitude and weak TCs. TCs experience time lags between the imposition of the shear and the weakening in TC intensity. A vertical shear of 8–9 m/s (9–10 m/s) would weaken TC intensity within 60 h (48 h). A vertical shear greater than 10 m/s would weaken TC intensity within 6 h. Finally, a statistical TC intensity prediction scheme is developed by using partial least squares regression, which produces skillful intensity forecasts when potential predictors include factors related to the vertical wind shear. Analysis of the standardized regression coefficients further confirms the obtained statistical results.  相似文献   

7.
热带气旋强度的卫星探测客观估计方法研究   总被引:3,自引:1,他引:3       下载免费PDF全文
运用2001~2002年6月10日到8月10日西北太平洋上15个热带气旋的日本静止气象卫星(GMS-5)高分辨率红外辐射亮温(TBB)资料,从代表热带气旋强度的热带气旋云系特征中选取多个TBB因子及热带气旋中心所在纬度等因子,结合《热带气旋年鉴》资料,采用逐步回归方法,经多次试验求得热带气旋强度的客观估计算式,复相关系数达到0.80以上。由此式估计的热带气旋中心气压,经24 h滑动平均后的结果与年鉴气压的复相关系数提高到0.89,且两者之间的气压差在±10 hPa之内的占整个样本数的83%以上,与美国新近展示的相关研究结果十分接近。用2000年两个热带气旋作检验,结果很好。该客观方法有望替代目前的Dvorak主观估计方法,成为一种新的业务方法。  相似文献   

8.
钱燕珍  张程明  孙军波  陈佩燕 《气象》2013,39(6):710-718
这是一个对GFS数值预报产品进行解释应用的方法.将支持向量机(SVM)回归方法应用于近海和登陆热带气旋(TC)的强度预报.从其本身强度,影响范围内气象因子情况,地形因子等三个方面,设计相关因子,建立预报模式,用来预报12、24、36、48、60和72 h的TC强度.总体上模式强度预报结果与中央气象台的预报结果相近,优于气候持续法的预报;趋势预报优势明显,可高出7~12个百分点.表明可以成为台风强度预报的另一个工具,投入业务应用.  相似文献   

9.
以三个西北太平洋热带气旋(TC)为例,利用WRF(Weather Research and Forecasting)模式进行了一系列海表温度(SST)敏感性数值试验,揭示了西北太平洋SST增暖对TC的强度、尺度及潜在破坏力的影响及其机理。结果表明,在距TC中心100 km以内区域的SST升高有利于TC强度增加,但会减小TC内核尺度;而在距TC中心100 km以外的SST升高并不会使TC强度明显增加甚至使TC强度减弱,但同时会增加TC内核尺度。伴随着低层向眼墙的入流,升高的外区SST会使TC区表层的大气温度和湿度升高,造成眼墙附近海气温差和湿度差及向内的气压梯度力减小,进而减少进入TC眼墙内的感热和潜热,不利于TC增强,但有利于眼墙向外移动,使TC内核尺度增加。内区SST升高与外区SST升高对TC强度及尺度变化的作用机理相反。因此,当TC移过冷或暖洋面时,TC的强度和尺度的变化不仅取决于局地洋面的冷或暖状况,还取决于TC内区和外区SST的差异。由于TC内区和外区SST对TC强度和内核尺度的不同作用,可能存在一个临界范围,当暖池范围在这个临界范围之内时TC潜在破坏力随暖池范围的扩大而增大,但当暖池范围超过这个临界范围时TC潜在破坏力便不会随着暖池范围的继续扩大而增大,甚至会有所减小。   相似文献   

10.
Summary Current understanding of tropical cyclone (TC) structure and intensity changes has been reviewed in this article. Recent studies in this area tend to focus on two issues: (1) what factors determine the maximum potential intensity (MPI) that a TC can achieve given the thermodynamic state of the atmosphere and the ocean? and (2) what factors prevent the TCs from reaching their MPIs? Although the MPI theories appear mature, recent studies of the so-called superintensity pose a potential challenge. It is notable that the maximum intensities reached by real TCs in all ocean basins are generally lower than those inferred from the theoretical MPI, indicating that internal dynamics and external forcing from environmental flow prohibit the TC intensification most and limit the TC intensity. It remains to be seen whether such factors can be included in improved MPI approaches.Among many limiting factors, the unfavorable environmental conditions, especially the vertical shear-induced asymmetry in the inner core region and the cooling of sea surface due to the oceanic upwelling under the eyewall region, have been postulated as the primary impediment to a TC reaching its MPI. However, recent studies show that the mesoscale processes, which create asymmetries in the TC core region, play key roles in TC structure and intensity changes. These include the inner and outer spiral rainbands, convectively coupled vortex Rossby waves, eyewall cycles, and embedded mesovortices in TC circulation. It is also through these inner core processes that the external environmental flow affects the TC structure and intensity changes. It is proposed that future research be focused on improving the understanding of how the eyewall processes respond to all external forcing and affect the TC structure and intensity changes. Rapid TC intensity changes (both strengthening and weakening) are believed to involve complex interactions between different scales and to be worthy of future research.The boundary-layer processes are crucial to TC formation, maintenance, and decaying. Significant progress has been made to deduce the drag coefficient on high wind conditions from the measurements of boundary layer winds in the vicinity of hurricane eyewalls by Global Positioning System (GPS) dropsondes. This breakthrough can lead to reduction of the uncertainties in the calculation of surface fluxes, thus improving TC intensity forecast by numerical weather prediction models.  相似文献   

11.
本文利用2010年1011号台风“凡亚比”登陆过程高分辨率数值模拟资料,诊断分析了“凡亚比”台风环流合力散度的垂直分布及其演变特征。结果指出,合力散度的显著区一直与台风系统相伴随,可以有效地示踪热带气旋(Tropical Cyclone,简称TC)的移动,并能较好地识别TC强度、结构的发展演变。台风中心偏东一侧流入层的合力散度异常信号首先出现并发展,反映出环流的非对称特征。随着台风趋于成熟,合力散度逐渐增强,高度扩展,对称性也逐渐增加;台风中心上空为合力辐合,外围为合力辐散,垂直方向上合力辐合与辐散相间的结构对应上升运动极值区及强降水,即对应台风眼墙位置。合力散度面积指数和强度指数的分析指出,垂直方向上辐合与辐散面积指数负相关;各层的合力辐合强度指数普遍大于辐散强度指数,垂直方向上两强度指数呈显著的正相关关系;结合面积指数与强度指数,可知垂直方向上合力辐合与辐散此消彼长。运用合力散度方程对该垂直分布特征的成因展开分析,发现风速u分量平流随经度变化项和风速v分量平流随纬度的变化项是TC眼区合力辐合部分的主要贡献项,垂直运动项决定了TC眼墙的合力辐合与辐散相间的垂直分布特征。  相似文献   

12.
沿海风工程设计风速中泊松-耿贝尔法的应用   总被引:2,自引:0,他引:2       下载免费PDF全文
每年西北太平洋热带气旋(TC)发生的次数、移动路径和强度都是随机的,我国东南沿海各地每年受TC影响的次数便构成了某种离散型分布,而TC影响下的最大风速则可以构成某种连续型分布。该文采用上海台风研究所提供的1961—2006年TC中心风速和TC影响期间各台站大风资料,利用泊松-耿贝尔联合极值风速计算方法,计算了沿海各气象站TC影响大风的多年一遇风工程设计最大风速。结果表明:当观测资料样本序列较短,特别是像TC这样随机性很强的天气事件,泊松-耿贝尔联合极值算法更具优势;我国沿海地区有53.9%的台站50年一遇最大风速在25 m/s以下,最大风速大于42.5 m/s以上的台站分布于浙江的大陈岛、嵊山、石浦,福建的北茭和台山,广东的遮浪、上川岛和海南的西沙岛,在这些地区进行风电开发风险较大,应引起足够重视。  相似文献   

13.
应用NOAA气候预测中心提供的热带大气季节内振荡(MJO)客观业务指数及中国气象局上海台风研究所提供的西北太平洋热带气旋(TC)最佳路径资料集,定量统计榆验了MJO对夏季西北太平洋TC活动的调制作用.结果表明:MJO对TC的生成、强度、路径和登陆活动都有显著的调节作用.当高空辐合中心位于120°E~160°E(MJO位相3~5)时,西北太平洋TC生成偏少,且生成位置偏北;而当高空辐合中心位于10°W~70°E(MJO位相8~10)时,西北太平洋TC生成偏多,且生成位置偏南;随着TC强度加强,能达到显著调节作用的MJO位相逐渐减少,当高空辐合辐散中心位于70°E(MJO位相10)时,对TC强度调制最显著.在路径调节方面,MJO位相1~4和10时,TC活跃于菲律宾以东的西北太平洋上,主要路径为西北偏北行,可能登陆华东、华北;而位相5~8时,TC主要活跃在菲律宾附近及以西到南海,以偏西行路径为主,可能登陆华南.MJO对登陆华南TC也有显著影响.该定量统计检验结果可为TC活动季节内预测提供依据.  相似文献   

14.
Atmospheric Infrared Sounder (AIRS) temperature profiles during the period 2003–2013 are used to examine the warm core structures and evolution characteristics associated with the formation and development of western North Pacific (WNP) tropical cyclones (TCs). The warm core with a steady 1.5-K warming in the layer of 500–300 hPa occurs 24 h prior to tropical storm formation. Apparent eye warming extends upward to upper troposphere and downward to near surface after tropical storm formation. TC intensity shows a robust positive correlation with the warm core strength and has a weaker but still significant positive correlation with the warm core height (the weaker correlation is primarily attributed to the scattered warm core heights of weak TCs). Future 24-h intensity change of TCs has little correlation with the warm core height while it has a significant negative correlation with the warm core strength. Weak to moderate warm core at 500–200 hPa may be a necessary but not sufficient initial condition for TC rapid intensification. AIRS-observed warm core structures, in combination with other environmental factors, have the potential to improve the prediction of tropical storm formation and rapid intensification of WNP TCs.  相似文献   

15.
GRAPES-TCM业务试验结果分析   总被引:3,自引:2,他引:3  
朱振铎  端义宏  陈德辉 《气象》2007,33(7):44-54
应用数值模式GRAPES-TCM2.1(GT2.1)的升级版GRAPES-TCM2.6(GT2.6)对2005年西北太平洋及南海热带气旋路径进行了后报试验和性能检验。在检验预报路径时考察了平均距离误差、系统偏差、距离误差的地理分布和平均移速误差,并在总体检验的基础上,根据与热带气旋路径相关的重要因素(路径类型、强度、有无登陆过程以及有无移向或移速突变)进行了分类检验。结果表明,GT2.6表现出良好的预报性能,其对所有样本的24、48和72小时平均距离误差分别为135.8、230.7和336.0km,但前12小时距离误差较大(近100km)。用系统偏差订正GT2.6对转向类样本的预报路径可获得明显的改进。GT2.6对130°E以西区域特别是靠近我国东南沿海的热带气旋路径预报表现较好。GT2.6的总体平均预报移速在48小时内都比最佳路径移速偏快,初始的12小时内偏快最多(近1m·s^-1),在48小时后转为偏慢。GT2.6和GT2.1的前48小时路径预报性能的稳定性相当,总体误差特征相似,GT2.6对所有样本的48小时预报有显著改进,平均距离误差减小17km。初步分析表明,对GT2.6的改进应主要集中在数值预报初始场的改进方面,比如在初始化过程中采用较准确的涡旋分离方案、加入引导向量改进初始移动以及适当考虑加入非对称Bogus涡旋;扩大预报区域,并适当选取预报区域起始位置等。  相似文献   

16.
To improve understanding of essential aspects that influence forecasting of tropical cyclones (TCs), the National Key Research and Development Program, Ministry of Science and Technology of the People’s Republic of China conducted a five-year project titled “Key Dynamic and Thermodynamic Processes and Prediction for the Evolution of Typhoon Intensity and Structure” (KPPT). Through this project, new understandings of TC intensification, including outer rainband-driven secondary eyewall formation and the roles of boundary layer dynamics and vertical wind shear, and improvements to TC data assimilation with integrated algorithms and adaptive localizations are achieved. To promote a breakthrough in TC intensity and structure forecasting, a new paradigm for TC evolution dynamics (i.e., the correlations, interactions, and error propagation among the triangle of TC track, intensity, and structure) is proposed; and an era of dynamic-constrained, big-data driven, and strongly coupled data assimilation at the subkilometer scale and seamless prediction is expected.  相似文献   

17.
本文利用热带测雨卫星TRMM(Tropical Rainfall Measuring Mission)微波成像仪TMI(TRMM Microwave Imager)2A12 水汽凝结物(Hydrometeor)反演资料,对西北太平洋地区从1998~2009 年的236 个热带气旋个例的1776 个“快照”(snapshot)的水汽凝结物的结构特征进行了分析,并探讨了水汽凝结物的时空变化与热带气旋强度演变联系。研究结果表明:(1)TMI 2A12 水汽凝结物资料显示出了热带气旋内部的细致结构及变化特征,水汽凝结物的峰值集中于数十公里到一百多公里的热带气旋眼壁及云墙区;在热带气旋发展过程中,随着热带气旋强度的增强,水汽凝结物增多且往其中心靠拢,从发展阶段到成熟阶段,水汽凝结物的大值中心基本上集中在距离热带气旋中心约50 km 区域,而且强度越强的热带气旋,水汽凝结物的大值中心与热带气旋中心的距离越近;在热带气旋消亡的过程中,水汽凝结物不断减弱且往外围扩散,逐渐扩展到远离中心的区域;(2)热带气旋强度与水汽凝结物的分布关系密切,热带气旋强度变化与热带气旋中心附近200 km 范围内的水汽凝结物含量存在显著的正相关,而200 km 以外的外围水汽凝结物含量存在负相关;(3)热带气旋强度变化与水汽凝结物的变化存在时间差,水汽凝结物的变化超前于热带气旋强度的变化,在热带气旋迅速发展之前数小时,热带气旋中心0~50 km 环状区域的水汽凝结物含量就已经提前增加了,在热带气旋减弱前数小时到十数小时,即使热带气旋还处于它强度的鼎盛时期,其中心0~50 km 环状区域的水汽凝结物含量就已经提前显著减少了,这种水汽凝结物的变化超前于热带气旋强度的变化的现象,可能是热带气旋强度预报的潜在线索。  相似文献   

18.
A new Tropical Cyclone (TC) initialization method with the structure adjustable bogus vortex was applied to the forecasts of track, central pressure, and wind intensity for the 417 TCs observed in the Western North Pacific during the 3-year period of 2005–2007. In the simulations the Final Analyses (FNL) with 1° × 1° resolution of National Center for Environmental Prediction (NCEP) were incorporated as initial conditions. The present method was shown to produce improved forecasts over those without the TC initialization and those made by Regional Specialized Meteorological Center Tokyo. The average track (central pressure, wind intensity) errors were as small as 78.0 km (11.4 hPa, 4.9 m s?1) and 139.9 km (12.4 hPa, 5.5 m s?1) for 24-h and 48-h forecasts, respectively. It was found that the forecast errors are almost independent on the size and intensity of the observed TCs because the size and intensity of the bogus vortex can be adjusted to fit the best track data. The results of this study indicate that a bogus method is useful in predicting simultaneously the track, central pressure, and intensity with accuracy using a dynamical forecast model.  相似文献   

19.
热带气旋眼墙非对称结构的研究综述   总被引:2,自引:0,他引:2  
热带气旋的眼墙非对称结构与其发展过程密切相关。在热带气旋移动过程中,非对称风场伴随着边界层内非对称摩擦而引起的辐合,影响着热带气旋眼墙内的对流分布。此外,风垂直切变作为影响热带气旋强度的重要因子,将上层暖心吹离表层环流,引起眼墙垂直运动的非对称,导致云、降水在方位角方向的非均匀分布。当存在平均涡度的径向梯度时,罗斯贝类型的波动可以存在于涡旋内核区域,影响眼墙非对称结构。海洋为热带气旋提供潜热和感热形式的能量,是热带气旋发展的重要能量来源,关于海洋如何影响热带气旋眼墙非对称结构的相关研究较少。文中着重回顾了热带气旋与海洋相互作用的研究成果,并提出海洋影响热带气旋眼墙非对称结构的机制。海洋对热带气旋最显著的响应特征是冷尾效应,该效应通过降低海表温度,减少海洋向大气输送的潜热和感热,从而影响热带气旋眼墙非对称结构。此外,海浪改变海表粗糙度,通过边界层影响移动热带气旋的眼墙结构。  相似文献   

20.
Based on the data(including radius of maximum winds) from the JTWC(Joint Typhoon Warning Center),the tropical cyclones(TCs) radii of the outermost closed isobar, TCs best tracks from Shanghai Typhoon Institute and the Black Body Temperature(TBB) of the Japanese geostationary meteorological satellite M1 TR IR1, and combining13 tropical cyclones which landed in China again after visiting the island of Taiwan during the period from 2001 to2010, we analyzed the relationship between the number of convective cores within TC circulation and the intensity of TC with the method of convective-stratiform technique(CST) and statistical and composite analysis. The results are shown as follows:(1) The number of convective cores in the entire TC circulation is well corresponding with the outer spiral rainbands and the density of convective cores in the inner core area increases(decreases) generally with increasing(decreasing) TC intensity. At the same time, the number of convective cores within the outer spiral rainbands is more than that within the inner core and does not change much with the TC intensity. However, the density of convective cores within the outer spiral rainbands is lower than that within the inner core.(2) The relationship described above is sensitive to landing location to some extent but not sensitive to the structure of TC.(3) The average value of TBB in the inner core area increases(decreases) generally with increasing(decreasing) of TC intensity, which is also sensitive to landing situation to some extent. At the same time, the average value of TBB within the outer spiral rainbands is close to that within the entire TC circulation, and both of them are more than that within the inner core. However, they do not reflect TC intensity change significantly.(4) The results of statistical composite based on convective cores and TBB are complementary with each other, so a combination of both can reflect the relationship between TC rainbands and TC intensity much better.  相似文献   

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

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