首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 734 毫秒
1.
用11年的全球无线电掩星数据(COSMIC),无线电探空数据(IGRA)以及欧洲中心再分析资料(ERA-Interim)对全球大气边界层高度(PBLH)进行估算比较.结果表明:(1)在1200 UTC和0000 UTC,由ERA-Interim和IGRA数据估算得到的全球PBLH空间分布较为一致,相关性较好,在白天正午时候太阳辐射能力较强,对流活动频繁,估算得到的大气边界层高度较高.(2)由COSMIC掩星数据估算得到的边界层高度比探空数据和再分析数据估算结果整体偏大.(3) COSMIC掩星数据,IGRA探空数据以及ERA-Interim再分析资料估算结果都表明边界层高度在低纬度地区偏大,高纬度地区偏小.(4)分析不同数据估算边界层高度纬向季节性差异表明,IGRA探空数据和COSMIC数据间差异为-1700m至-500m,IGRA与ERA-Interim之间的差异为-500m至250m.此外,对于大多数纬度而言,三个数据集之间的差异在冬季较大,在夏季较小.  相似文献   

2.
介绍了如何利用Micaps数据服务器对T639全球模式06UTC和18UTC(世界时)进行本地化,有效解决了使用中国气象局数值预报中心在2012年4月1日下发的T639全球模式解码程序操作步骤处理产品数量不齐全的问题,给预报员提供数量完整的T639全球模式预报产品.  相似文献   

3.
IGRA全球站点探空数据集介绍及其资料质量状况分析   总被引:2,自引:0,他引:2       下载免费PDF全文
美国国家气候数据中心(NCDC)研制的全球站点无线电探空资料数据集(The Integrated Global Radiosonde Archive,简称IGRA)包含全球1538个探空站点从1938年至今的多层次、多要素逐日探空资料。该数据是目前全球范围内时空密度最高、资料最为完整的一套高质量探空站点资料数据集。对该数据集概况及探空资料的质量控制过程进行了概述,并对该资料存在的质量问题进行分析。结果表明:该数据集虽然进行了质量控制,但是该数据集中部分资料的质量仍然存在问题,如该套资料中1973年之前的中国区域探空温度资料存在较大的问题,与国家气象信息中心整理的中国区域探空资料相比,两套资料中各站在1973年以前的不一致率都达到了70 %以上。另外,该套资料中风速资料在1975-1978年以及1981年错误较多,这5 a中出现的错误资料占所有错误资料的50.74 %。而这5 a中错误资料的分布也比较集中,风速出错站点主要分布在美国南部和阿拉斯加、非洲中部和南部以及俄罗斯。全球大部分探空台站通过全球通讯系统进行传输资料集中在1980年以后,因此本数据集中1980年以后的探空资料质量和完整性相对较高。结论和分析为用户使用该数据集提供了很好的参考依据。  相似文献   

4.
基于南极18个站点探空气象观测数据对欧洲中期天气预报中心的再分析数据(ERA-Interim)和美国国家环境预报中心的再分析数据(NECP)在南极地区高层大气的适用性进行验证。结果表明:在南极上空,随着高度抬升,探空气象观测数据与两套再分析数据中四个气象要素的差值均逐渐变大,再分析数据数值愈加偏离实际观测数值。两套再分析数据的位势高度和温度与探空观测数据偏差较小;风向则和探空观测数据相差甚远;两套再分析数据的风速与探空观测数据在300 hPa偏差较大。在季节变化中,南极的春季,再分析数据中的位势高度和温度与探空观测数据相差较大,在其他季节相差相对较小。再分析数据中的风速与探空观测数据在南极的夏季相差较小。再分析数据中的风向与探空观测数据存在较大偏差,且差值没有明显的季节变化。尽管两套再分析数据都存在很大偏差,但ERA-Interim数据整体上优于NCEP数据。对比分析也表明,采用这些再分析资料作为初始条件和边界条件驱动南极区域大气模式将带来较大的误差。未来需要加强南极探空观测,改进再分析资料同化和数值模拟系统。  相似文献   

5.
基于ERA5 PBLH数据集、PM2.5观测数据和气温、相对湿度等多种气象观测资料,采用统计方法对1980—2020年云南省大气边界层高度(PBLH)变化特征、变化趋势,PM2.5浓度变化与PBLH的相关性进行系统分析。结果表明,受太阳直接辐射在一天中的变化影响,云南PBLH夜间变化平稳,白天波动较强;地气温差年代际变化与PBLH具有较好的一致性,2000年代至2010年代云南PBLH呈稳步升高趋势。降雨年内分布、地理位置的特殊性和天气系统等影响差异使得云南高海拔区域年内PBLH变化差异较小,干季PBLH较中、低海拔区仍明显偏低。云南四季PBLH大值中心多出现在昆明至楚雄,变化百分率多表现为增大趋势,春季整体增大趋势最明显,冬季变化百分率空间分布差异性最强。近地面层低温、高湿、弱风造成傍晚至夜间(18:00—01:00)PBLH快速下降,致PM2.5堆积到近地面层,PBLH对PM2.5浓度影响随之增强。各州/市PBLH与PM2.5的相关性和相关程度随季节和时次不同存在一定差异,02:00和08:00 PBLH与PM2.5基本表现为负相关,相关程度整体高于14:00和20:00;人类活动...  相似文献   

6.
多套土壤温湿度资料在青藏高原的适用性   总被引:13,自引:0,他引:13  
刘川  余晔  解晋  周欣  李江林  葛骏 《高原气象》2015,(3):653-665
利用青藏高原中部和东部土壤温度和湿度观测资料,通过计算两套再分析资料(ERA-Interim和CFSR)和六套陆面模式资料(ERA/land、MERRA/land、GLDAS-NOAH、GLDAS-CLM、GLDAS-M OSAIC和GLDAS-VIC)分别与观测资料之间的平均偏差、偏差标准差、相关系数、标准差比等统计参数,结合Brunke排名法,综合评估了再分析资料和陆面模式资料中土壤温湿度数据在青藏高原的适用性。结果表明:对于土壤温度,CFSR与观测值最接近,其次是MERRA/land和GLDAS-CLM,而ERA-Interim和ERA/land与观测值相差较大;除GLDAS-CLM土壤温度比观测值偏高外,其他资料土壤温度在大部分站点比观测值偏低,其中ERA-Interim和ERA/land土壤温度比观测值偏低较多,部分站点平均偏差超过-20℃。对于非冻结期(5 10月)土壤湿度,GLDAS-CLM与观测值最接近,其次是GLDAS-NOAH或ERA-Interim;与观测值相比,CFSR、ERA-Interim和ERA/land的土壤湿度偏湿,平均偏差大部分在0.05~0.20 m3·m-3之间,而GLDAS-NOAH、GLDAS-CLM和GLDAS-M OSAIC的土壤湿度偏干。  相似文献   

7.
ERA-Interim地表温度资料在青藏高原多年冻土区的适用性   总被引:3,自引:0,他引:3  
地表温度综合反映了大气和地表植被、土壤等局地因素相互作用的能量交换状况,是许多冻土分布模型和寒区陆面过程模式的上边界条件,对多年冻土分布和活动层厚度估算具有重要意义。为了检验ERA-Interim再分析地表温度资料在青藏高原(下称高原)多年冻土区的适用性,综合比较了2011年1月1日至2012年12月31日期间高原不同类型多年冻土区3个综合观测场的观测地表温度和ERA-Interim再分析资料之间的偏差、均方差、相关系数、解释方差、均方根误差和平均绝对误差。结果表明,ERA-Interim再分析资料能较好地再现高原多年冻土区3个综合观测场地表温度的基本特征,并能较好地描述高原地表温度的季节变化。但ERA-Interim再分析年平均地表温度比观测值偏低,西大滩、五道梁和唐古拉站依次偏低1.7,1.0和0.9℃;地表温度的再分析值和观测值之间的相关系数和解释方差都较高,均方差也相近。ERA-Interim再分析地表温度资料对观测站点相对稀少且空间分布不均匀的高原多年冻土区具有较好的适用性,可以作为地表温度的有效代用资料。  相似文献   

8.
利用中国西北地区2015年9月至2016年8月38个站点L波段探空观测、2016年7月加密探空观测和ERA-Interim边界层高度资料,对比分析了西北地区大气边界层高度变化特征。观测资料表明,在中国西北地区,08:00(北京时,下同)冬季边界层高度最高; 20:00春季边界层高度最高,边界层高度从西部到东部有显著降低的趋势。ERA-Interim资料基本能表现出边界层高度的区域分布,但相对于探空观测得到的边界层高度,除夏季20:00外,ERA-Interim再分析资料边界层高度均偏低。全年平均而言,08:00(20:00)偏低160 m(170 m),其中在08:00(20:00),冬季(春季)偏低最显著。08:00边界层高度与低层稳定度、近地层温度和风速相关更加显著; 20:00边界层高度与低层稳定度和相对湿度相关更加显著。2016年7月加密观测资料对比表明,ERA-Interim资料的对流(中性)边界层高度显著偏高;低层稳定度、相对湿度偏小,风速偏大可能是造成边界层高度偏高的原因; ERA-Interim资料的稳定边界层高度偏低,与低层稳定度和近地层温度偏低相关,但其影响因素相对更加复杂。  相似文献   

9.
研究利用实测数据对遥感产品(ECV CCI)、全球陆面同化系统产品(GLDAS/CLM、NOAH2.7、NOAH3.3、VIC、MOSAIC)、再分析资料(ERA-Interim、NCEP/DOE)等不同来源的8套长时间序列土壤湿度产品进行时空比较。结果表明:上述产品均可以模拟出中国不同区域土壤湿度的空间分布;从各产品平均态与实测数据的偏差、均方根误差、相关系数统计结果来看,NOAH模式和ERA-Interim表现最好,MOSAIC和NCEP/DOE与实测相差较大;从区域尺度上来看,所有产品在东北区域表现最好,在华南、西南南区和西北西区较差;ERA-Interim有效的模拟了土壤湿度的年际变化。  相似文献   

10.
利用1979-2009年NCEP/CFSR全球大气边界层高度(PBLH)、感热通量月平均资料,运用多种统计方法,探讨了东亚地区大气边界层高度和感热通量的变化特征,并研究了两个要素场之间的相互关系。结果表明:夏季大气边界层高度东部增高、西部降低,而冬季则与夏季相反。夏季青藏高原地区感热有减少的趋势,其余地区以增加为主;冬季东部地区及新疆西部感热减少,其余地区感热增加。夏季奇异值分析第一模态表明青藏高原地区的感热通量减少时,相应地区的PBLH降低;内蒙古东部和东北地区的感热增加时,相应地区的PBLH增加。该空间型在一定程度上反映了青藏高原、内蒙古东部及东北地区,PBLH变化主要受地表加热影响,且年代际变化显著。夏季奇异值分析第二模态表明当青藏高原南侧、华北的感热增加(减少)时,相应地区的PBLH升高(降低),该空间型年际变化较显著,一定程度上反映了大气边界层高度及感热在高原主体及其南侧分布的差异。冬季奇异值分析第一模态主要表现为在我国东西部反相变化的分布;冬季奇异值分析第二模态表现为华北和华东地区与其余地区反相相关的分布。冬季两个空间型均具有一定的年代际变化特征。在青藏高原及其南侧和西北干旱半干旱区部分地区,PBLH变化主要受地表加热的影响,而在东部季风区,感热变化仅是影响PBLH变化的因素之一。  相似文献   

11.
The global planetary boundary layer height (PBLH) estimated from 11 years (2007–17) of Integrated Global Radiosonde Archive (IGRA) data, Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) soundings, and European Center for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) data, are compared in this study. In general, the spatial distribution of global PBLH derived from ERA-Interim is consistent with the one from IGRA, both at 1200 UTC and 0000 UTC. High PBLH occurs at noon local time, because of strong radiation energy and convective activity. There are larger differences between the results of COSMIC and the other two datasets. PBLHs derived from COSMIC are much higher than those from radiosonde and reanalysis data. However, PBLHs derived from the three datasets all exhibit higher values in the low latitudes and lower ones in the high latitudes. The latitudinal difference between IGRA and COSMIC ranges from −1700 m to −500 m, while it ranges from −500 m to 250 m for IGRA and ERA-Interim. It is found that the differences among the three datasets are larger in winter and smaller in summer for most studied latitudes.摘要用11年的全球无线电掩星数据 (COSMIC) , 无线电探空数据 (IGRA) 以及欧洲中心再分析资料 (ERA-Interim) 对全球大气边界层高度 (PBLH) 进行估算比较. 结果表明: (1) 在1200 UTC和0000 UTC, 由ERA-Interim和IGRA数据估算得到的全球PBLH空间分布较为一致, 相关性较好, 在白天正午时候太阳辐射能力较强, 对流活动频繁, 估算得到的大气边界层高度较高. (2) 由COSMIC掩星数据估算得到的边界层高度比探空数据和再分析数据估算结果整体偏大. (3) COSMIC掩星数据, IGRA 探空数据以及 ERA-Interim 再分析资料估算结果都表明边界层高度在低纬度地区偏大, 高纬度地区偏小. (4) 分析不同数据估算边界层高度纬向季节性差异表明, IGRA探空数据和COSMIC数据间差异为-1700m至-500m, IGRA与ERA-Interim之间的差异为-500m至250m.此外, 对于大多数纬度而言, 三个数据集之间的差异在冬季较大, 在夏季较小.  相似文献   

12.
The CRA-Interim trial production of the global atmospheric reanalysis for 10 years from 2007 to 2016 is carried out by the China Meteorological Administration (CMA) in 2017. The structural characteristics of the horizontal shear line over the Tibetan Plateau (TPHSL) based on the CRA-Interim datasets are examined by objectively identifying the shear line, and compared with the analysis results of the European Centre for Medium-Range Weather Forecasts (ECMWF)the European Center (ECMWF) reanalysis data (ERA-Interim). The case occurred at 18UTC on July 5, 2016. The results show that both of the ERA-Interim and CRA-Interim datasets can well reveal the circulation background and the dynamic and thermal structure characteristics of TPHSL, and they have shown some similar features. The middle and high latitudes at 500 hPa are characterized by500hPa present the circulation situation of  相似文献   

13.
Accurate surface air temperature (T2m) data are key to investigating eco-hydrological responses to global warming. Because of sparse in-situ observations, T2m datasets from atmospheric reanalysis or multi-source observation-based land data assimilation system (LDAS) are widely used in research over alpine regions such as the Tibetan Plateau (TP). It has been found that the warming rate of T2m over the TP accelerates during the global warming slowdown period of 1998–2013, which raises the question of whether the reanalysis or LDAS datasets can capture the warming feature. By evaluating two global LDASs, five global atmospheric reanalysis datasets, and a high-resolution dynamical downscaling simulation driven by one of the global reanalysis, we demonstrate that the LDASs and reanalysis datasets underestimate the warming trend over the TP by 27%–86% during 1998–2013. This is mainly caused by the underestimations of the increasing trends of surface downward radiation and nighttime total cloud amount over the southern and northern TP, respectively. Although GLDAS2.0, ERA5, and MERRA2 reduce biases of T2m simulation from their previous versions by 12%-94%, they do not show significant improvements in capturing the warming trend. The WRF dynamical downscaling dataset driven by ERA-Interim shows a great improvement, as it corrects the cooling trend in ERA-Interim to an observation-like warming trend over the southern TP. Our results indicate that more efforts are needed to reasonably simulate the warming features over the TP during the global warming slowdown period, and the WRF dynamical downscaling dataset provides more accurate T2m estimations than its driven global reanalysis dataset ERA-Interim for producing LDAS products over the TP.  相似文献   

14.
大气边界层高度是天气、气候、大气环境研究中的一个重要参数,目前尚缺少基于激光雷达探测系统反演全天边界层高度的有效方法。文中利用北京朝阳站、大兴站的激光云高仪数据,首先评估了梯度法、标准偏差法、曲线拟合法和小波协方差法反演边界层高度的适用性和局限性,发现梯度法容易受环境噪声的影响,曲线拟合法稳定性较好,但在夜间弱湍流条件下会将残留层高度误判为夜间边界层高度。提出两步曲线拟合法,将一天中边界层结构划分为白天的对流边界层、夜间的残留层和稳定边界层,通过用不同的理想曲线对其进行两步拟合,获取全天边界层高度的变化。将两步曲线拟合法的反演结果与基于L波段探空雷达的位温梯度法的探测结果进行比较发现:两者相关系数为0.91,证明了两步曲线拟合法的可行性以及激光云高仪探测边界层高度的应用潜力。采用该方法反演2017年5—6月朝阳站与大兴站边界层高度,对比发现:城区特殊的下垫面性质使朝阳站日间对流边界层发展更早,边界层高度更高,全天朝阳站边界层高度的变化在308—1391 m,大兴站在197—1302 m。   相似文献   

15.
Tropical cyclone(TC) genesis over the western North Pacific(WNP) is analyzed using 23 CMIP5(Coupled Model Intercomparison Project Phase 5) models and reanalysis datasets. The models are evaluated according to TC genesis potential index(GPI). The spatial and temporal variations of the GPI are first calculated using three atmospheric reanalysis datasets(ERA-Interim, NCEP/NCAR Reanalysis-1, and NCEP/DOE Reanalysis-2). Spatial distributions of July–October-mean TC frequency based on the GPI from ERA-interim are more consistent with observed ones derived from IBTr ACS global TC data. So, the ERA-interim reanalysis dataset is used to examine the CMIP5 models in terms of reproducing GPI during the period 1982–2005. Although most models possess deficiencies in reproducing the spatial distribution of the GPI, their multimodel ensemble(MME) mean shows a reasonable climatological GPI pattern characterized by a high GPI zone along 20?N in the WNP. There was an upward trend of TC genesis frequency during 1982 to 1998, followed by a downward trend. Both MME results and reanalysis data can represent a robust increasing trend during 1982–1998, but the models cannot simulate the downward trend after 2000. Analysis based on future projection experiments shows that the GPI exhibits no significant change in the first half of the 21 st century, and then starts to decrease at the end of the 21 st century under the representative concentration pathway(RCP) 2.6 scenario. Under the RCP8.5 scenario, the GPI shows an increasing trend in the vicinity of20?N, indicating more TCs could possibly be expected over the WNP under future global warming.  相似文献   

16.
基于1992~2010年全国778个农业气象站土壤湿度观测资料、ERA-Interim、JRA55、NCEP-DOE R2和20CR土壤湿度再分析资料,通过平均差值、相关系数、差值标准差、标准差比四个参数,利用Brunke排名方法和EOF(Empirical Orthogonal Function)分析,对四套土壤湿度再分析资料在中国西北东部—华北—江淮区域的适用性进行了分析。主要结论如下:不同季节的平均偏差空间分布上,JRA55资料同观测数据的平均偏差在±0.08m~3 m~(-3)之间,春、夏季西北东部JRA55土壤湿度偏小,ERA-Interim、NCEP-DOE R2、20CR资料较观测数据偏湿,华北南部、江淮地区平均偏差小于西北东部、华北北部。在年际变化上,各个季节ERA-Interim资料同观测资料最为接近,能稳定地再现西北东部、华北、江淮地区土壤湿度干湿变化趋势,反映出重要的旱涝年。整体而言,四套再分析资料中ERA-Interim资料同观测资料接近,JRA55、NCEP-DOE R2资料次之,20CR资料最差。  相似文献   

17.
Cloud radiative effects(CREs) at the top of the atmosphere(TOA) in three reanalysis datasets(the European Center for Medium-Range Weather Forecasts Interim Reanalysis(ERA-Interim), the Japanese 55-yr Reanalysis Project(JRA-55), and the Modern-Era Retrospective Analysis for Research and Applications(MERRA)) are evaluated using recent satellite-based observations. The reanalyses can basically capture the spatial pattern of the annual mean shortwave CRE, but the global mean longwave CRE in ERA-Interim and JRA55 is weaker than observed, leading to overestimations of the net CRE. Moreover, distinct CRE biases of the reanalyses occur in the Intertropical Convergence Zone(ITCZ), coastal Pacific and Atlantic regions, and East Asia. Quantitative examination further indicates that the spatial correlations of CREs and TOA upward radiation fluxes with corresponding observations in ERA-Interim are better than in the other two reanalyses. Although MERRA has certain abilities in producing the magnitudes of global mean CREs, its performance in terms of spatial correlations in winter and summer are worse than for the other two reanalyses. The ability of JRA55 in reflecting CREs lies between the other two datasets. Compared to the global mean results, the spatial correlations of shortwave CRE in East Asia decrease and the biases of regional mean CREs increase in the three reanalyses. This implies that, currently, it is still difficult to reproduce East Asian CREs based on these reanalyses. Relatively, ERA-Interim describes the seasonal variation of East Asian CREs well, albeit weaker than observed. The present study also suggests that in-depth exploration of the ability of reanalysis data to describe aspects relating to cloud properties and radiation is needed using more comprehensive observations.  相似文献   

18.
利用农业气象站观测资料对长江中下游地区1988-2010年遥感土壤湿度进行了验证,并与NCEP和ERA-Interim土壤湿度做了对比分析。研究表明,ECV遥感土壤湿度冬季平均土壤湿度最高,春季和秋季次之,夏季平均土壤湿度最低;这种季节性干湿变化与农业气象站观测资料一致。但是,NCEP和ERA-Interim土壤湿度再分析资料,则夏季平均土壤湿度高,春季和秋季次之,而冬季平均土壤湿度最低;这种季节性变化与ECV遥感土壤湿度和农业气象站观测资料呈反位相。就年际变化而言,ECV遥感土壤湿度与农业气象站观测资料和两套再分析资料均有较高的一致性,并在春季和秋季最高,尤其是在长江以北地区和长江以南洞庭湖、鄱阳湖两大湖区,相关系数达到0.7~0.9;而夏季一致性最低,相关系数仅为0.4左右。在研究时段,ECV土壤湿度在冬季明显增加,在夏季则有明显下降趋势。  相似文献   

19.
利用ERA-Interim、MERRA和NCEP/NCAR三套再分析资料,分析1979~2014年夏季青藏高原大气水汽含量的时空变化特征,同时对比了各套资料异同点,结果表明:(1) ERA-Interim和MERRA资料均显示出夏季青藏高原大气水汽含量呈现显著的上升趋势,在1994~1995年前后发生明显突变,大气水汽含量由偏低时期向偏高时期转变;而NCEP/NCAR资料并没有出现类似的显著上升趋势和突变年份;ERA-Interim资料与MERRA资料的夏季青藏高原湿池指数之间的相关性明显强于NCEP/NCAR资料与它们任何一个之间的相关性。(2)夏季青藏高原大气水汽含量呈现出自高原东南边缘地区向西北部递减的分布形式。其中,MERRA与ERA-Interim资料显示的水汽含量分布更为相似,而NCEP/NCAR资料反映的水汽含量在高原中部往北递减不明显,湿度中心较为分散。(3)从空间分布上,MERRA与ERA-Interim资料显示青藏高原大部分地区夏季水汽含量均呈显著的增加趋势,而NCEP/NCAR资料仅在高原东北部小部分区域存在显著的增加趋势。(4)从夏季青藏高原大气水汽含量的年际变化特征分析来看,ERA-Interim和MERRA资料相对于NCEP/NCAR资料也更为接近。   相似文献   

20.
《大气与海洋》2012,50(4):262-278
ABSTRACT

This study evaluates the 1981–2010 spatiotemporal differences in six available climate datasets (daily total precipitation and mean air temperature) over the Lower Nelson River Basin (LNRB) in ten of its sub-watersheds at seasonal and annual time scales. We find that the Australian National University spline interpolation (ANUSPLIN), and inverse distance weighted (IDW) interpolated observations from 14 Environment and Climate Change Canada (ECCC) meteorological stations show dry biases, whereas reanalysis products tend to overestimate precipitation across most of the basin. All datasets exhibit prominent disagreement in precipitation trends whereby the European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) and European Union Water and Global Change (WATCH) Forcing Data ERA-Interim (WFDEI) show exceptional wetting trends, while the IDW and ANUSPLIN data manifest drying trends. Mean air temperature trends generally agree across most of the datasets; however, the North American Regional Reanalysis (NARR) and IDW show stronger warming relative to other datasets. Overall, analyses of the different climate datasets and their ensemble reveal that the choice of input dataset plays a crucial role in the accurate estimation of historical climatic conditions, particularly when assessing trends, for the LNRB. Using the ensemble has the distinct advantage of preserving the unique strengths of all datasets and affords the opportunity to estimate the uncertainty for hydrologic modelling and climate change impact studies.  相似文献   

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

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