首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
近地层能量闭合度对陆面过程模式影响   总被引:1,自引:0,他引:1       下载免费PDF全文
大量近地层观测试验表明,利用涡动相关法观测的湍流通量小于近地层可利用能量,即近地层能量是不闭合的,这种不闭合度一般为20%甚至更高.而陆面过程模式是基于地气间能量平衡建立,并且模式中的湍流边界层参数化方案通常根据实际观测的湍流通量来确定,因此能量不闭合必将对陆面过程模式造成一定的影响.本文利用2007年春季SACOL站的近地层观测资料,依据能量守恒将能量不闭合中的残余能量通过波文比分配到观测的湍流通量中,即修正涡动相关法观测的湍流通量使得近地层能量达到平衡;之后分别利用观测和修正的湍流通量,建立了能量不闭合和闭合情形下的湍流参数化方案,借助陆面过程模式SHAW,通过数值模拟和对比分析方法考察近地层能量闭合度对陆面过程模式的影响.研究结果表明近地层能量闭合对陆面过程模式有显著的影响:在陆面过程数值模拟中,当应用近地层能量不闭合的湍流通量形成的湍流参数化方案时,陆面过程模式会明显高估地表长波辐射及土壤温度;但当应用修正湍流通量使得近地层能量达到闭合形成的湍流参数化方案后,在不改变任何地表土壤物理生化属性的情况下,陆面过程模式能较好地模拟地表长波辐射和土壤温度.  相似文献   

2.
沙漠陆面过程参数化与模拟   总被引:4,自引:2,他引:2       下载免费PDF全文
郑辉  刘树华 《地球物理学报》2013,56(7):2207-2217
沙漠地区植被稀疏、干旱少雨,其陆面物理过程具有与全球其它地区显著不同的特点.本文利用巴丹吉林沙漠观测资料,分析和计算了地表反照率、比辐射率、粗糙度和土壤热容量、热传导系数等关键陆面过程参数,建立了适合于沙漠地区的陆面过程模式DLSM (Desert Land Surface Model),并与NOAH陆面过程模式的模拟结果和观测资料进行了比较.结果表明:巴丹吉林沙漠地表反照率为0.273,比辐射率为0.950,地表粗糙度为1.55×10-3 m,土壤热容量和热扩散系数分别为1.08×106 J·m-3·K-1和3.34×10-7 m2·s.辐射传输、感热输送和土壤热传导过程是影响沙漠地区地表能量平衡的主要物理过程.通过对这三种过程的准确模拟检验,DLSM能够较准确地模拟巴丹吉林沙漠地气能量交换特征;短波辐射、长波辐射和感热通量的模拟结果与观测值间的标准差分别为7.98,6.14,33.9 W·m-2,与NOAH陆面过程模式的7.98,7.72,46.6 W·m-2的结果接近.地表反照率是沙漠地区最重要的陆面过程参数,地表反照率增大5%,向上短波辐射通量随之增加5%,感热通量则减小2.8%.本文研究结果对丰富陆面过程参数化方案,改进全球陆面过程模式、气候模式具有参考意义.  相似文献   

3.
陆面过程模式中有关土壤水热传输、植被冠层和空气动力学等过程的关键参数的不确定性严重制约着地表-大气相互作用模拟能力的提高.本文利用架设在我国吉林通榆和甘肃榆中典型半干旱区陆-气相互作用野外试验站的观测资料,结合大气边界层理论,利用多种方法系统估算了上述试验站地表空气动力学粗糙度(z0m)以及热传输附加阻尼(?B?1)的量值.结果表明,z0m在半干旱区具有明显的季节变化和年际变化特征,且在植被低矮的下垫面与现行通用的陆面模式中的默认值相差较大;而?B?1的日变化和季节变化特征明显.将修正后的z0m和?B?1参数化方案引入陆面过程模式,发现能够明显改善模式对于半干旱区地表感热通量的模拟能力.这些结果说明有必要进一步结合半干旱区的野外观测试验对陆面模式在该地区的缺省参数设置进行更广泛的评估,而基于外场观测试验和大气边界层理论估算的关键地表参数对于改进陆气相互作用的模拟体现出较大的应用潜力.  相似文献   

4.
降水条件下的典型干旱区陆面特征模拟验证   总被引:6,自引:0,他引:6       下载免费PDF全文
利用“我国西北干旱区陆-气相互作用观测试验"在敦煌双墩子戈壁站取得的观测资料及最近的一些研究成果对陆面模式中反照率、相似性函数及地表粗糙度(动量粗糙度,标量粗糙度)3个方面的参数化方案进行了改进,然后对一次典型降水过程的陆面特征及近地层的风、温、湿进行了模拟.结果表明:改进的模式能对降水条件下的干旱区陆面特征进行较好的模拟,其中对辐射、地表温度的模拟相当好,而对能量的模拟虽然还有要改进的地方,但总体结果令人满意;另外改进后的模式对近地层的温度和湿度的模拟也有明显的改善.  相似文献   

5.
基于卫星遥感资料的中国区域土壤湿度EnKF数据同化   总被引:6,自引:0,他引:6       下载免费PDF全文
土壤湿度在陆气相互作用过程中扮演着重要的角色,是气候、水文、农业、林业等研究中重要的地球物理参数之一.土壤湿度影响地面蒸散,径流、地表反射率、地表发射率以及地表感热和潜热通量,从而对气候有重要影响,它对大气的影响在全球尺度上仅次于海面温度,在陆地尺度其影响甚至超过海面温度.本文介绍了基于EnKF及陆面过程模型的中国区域陆面土壤湿度同化系统(CLSMDAS,China Land Soil Moisture Data Assimilation System),以及该系统应用于中国区域陆面土壤湿度同化试验的结果.CLSMDAS包括以下几个部分:1)陆面模式采用美国国家大气研究中心NCAR的陆面过程模型Community Land Model Version3.0(简写为CLM3.0);2)大气驱动场数据中的降水和地面入射太阳辐射数据来自FY2静止气象卫星每小时产品;3)陆面数据同化方法采用EnKF(Ensemble Kalman Filter)同化方法;4)观测数据包括AMSR-E卫星反演土壤湿度产品以及地面土壤湿度观测资料.利用CLSMDAS对2006年6~9月的土壤湿度同化试验结果的分析表明:陆面模式模拟和同化结果都能比较合理地反映出土壤湿度时空分布,同化的土壤湿度分布与2006年8月重庆、四川发生建国以来最严重的夏伏旱有非常好的对应关系,与发生在9月的湖北东部、广西南部等地的干旱区也有非常好的对应关系.  相似文献   

6.
区域气候模式对末次盛冰期东亚季风气候的模拟研究   总被引:14,自引:2,他引:14  
末次盛冰期(LGM)是距今最近的一个与现代环境反差最大的气候时期. 利用包含较为详细陆面过程的区域气候模式, 通过分别加入现代植被和根据花粉化石资料转化的东亚地区LGM古植被, 模拟了LGM东亚季风气候并研究了植被变化对LGM东亚季风的影响. 由区域气候模拟得到的较为精细的气候演变图像表明: LGM东亚大陆全年降温是导致东亚冬季风强盛、夏季风萎缩的重要原因; 夏季西太平洋副热带高压的西伸、加强, 是造成中国东部LGM夏季降水减少的重要原因. LGM青藏高原及中亚地区的降水及有效降水均有所增加, 高原有效降水的增加主要由降水增加所造成, 地表蒸散对其贡献较小. LGM青藏高原的积雪也有所增加, 有利于高原地区的冰川、冻土发育, 使得该时期的多年冻土区可向南扩展到30°N以南. 在LGM模拟中加入恢复的古植被会放大由外强迫造成的气候影响, 对于模拟的降温、降水变化、地表热平衡量的变化、积雪及其他气候参量的变化都有进一步的强化作用, 使模拟结果与有关地质资料更为接近.  相似文献   

7.
地表覆盖是陆面及气候模式中的重要基础变量,其数据质量对气候模式性能有显著影响.本文基于一套中国自主研制的全球30m地表覆盖数据(GlobeLand30),利用北京气候中心气候系统模式(Beijing Climate CenterClimateSystemModel,BCC_CSM)中的陆面过程和大气环流分量模式,开展GlobeLand30数据对气候模式性能影响的研究.首先通过GlobeLand30和其他卫星遥感等辅助数据融合细化植被功能型(Plant FunctionType,PFT)类型数据满足BCC_CSM模式需求,然后利用基于面积比例的升尺度方法得到适合于全球模式的不同覆盖类型及其面积百分比数据.GlobeLand30与模式原有覆盖数据都能合理描述全球地表覆盖基本分布特征,但也存在一定差异,其中植被PFT类型的差别最明显.通过数值模拟试验全面评估GlobeLand30数据对模式结果的影响,结果表明在BCC_CSM的陆面模式和大气模式中采用GlobeLand30数据可以合理再现陆面与大气的基本气候特征;更新植被PFT数据的气候效应大于更新冰川和水体类型,综合更新所有覆盖类型数据对模拟结果的影响最大.引入GlobeLand30数据可降低BCC_CSM大气模式对北半球中高纬地区降水的正偏差和南美亚马逊等地区模拟降水的负偏差,以及南半球部分地区大气温度模拟的负偏差,因此GlobeLand30数据适用于BCC_CSM分量模式并对改善模式性能有正贡献.  相似文献   

8.
使用地球观测系统的中分辨率成像光谱仪(EOS.MODIS)提供的归一化植被指数(NDVI)产品估算植被覆盖度和航天飞机雷达地形测绘任务(SRTM)制成的数字高程模型(DEM)数据遥感产品替换WRF模式默认的植被覆盖度和地形高度,并且利用WRF模式及其先进的三维变分同化系统(WRF.3DVar)循环同化东北半干旱区自动气象站近地面气象要素,对东北半干旱区的温度场、湿度场、风场和能量场的结构及其日变化特征进行了较为细致的模拟研究.通过4组数值模拟试验分别探讨了同化气象要素与改变模式地表参数引起的不同下垫面潜热、感热的分配关系和降水、土壤湿度变化弓f起的地表能量通量模拟效果,并利用通榆站、奈曼站、锦州站、和密云站2009年6-8月的通量观测资料与模拟结果对比检验.结果表明,WRF模式能够较好地模拟出东北半干旱区夏季的近地面温度、风向、净辐射、感热和潜热等要素的变化特征及日变化规律.同化试验(Case2)模拟的近地面气温、相对湿度、风速相比控制性试验(Case1)有所改善;陆面参数试验(Case3)和集合试验(Case4)改善了感热和地表热通量的模拟.WRF模式能较好地模拟出下垫面土壤湿度随时间变化的规律,集合试验(Case4)土壤湿度模拟结果与4个通量站观测值相比无太大差别,但降水的模拟有待改善.本研究利用卫星遥感资料改善模式下垫面陆面参数,利用气象资料同化改善近地面大气要素模拟精度,这是将各种不同空间和时间尺度的多源数据与数值模拟融合的有益尝试.此研究生成的东北地区资料同化数据集可用于气候变化、干旱监测等方面,对深入了解半干旱区气候的形成和维持机理具有重要的意义.  相似文献   

9.
沙漠绿洲地区夏季地表能量收支的数值模拟   总被引:8,自引:4,他引:4       下载免费PDF全文
本文在MSPAS(Modified Soil-Plant-Atmosphere Scheme)的基础上,引入了一个有效的晴天大气辐射传输方案,建立了一个能在物理上真实地模拟陆气相互作用及其反馈机制的二维模式MLAIM(Modified Land Atmosphere Interaction Model).本文利用HEIFE实验的观测资料对MLAIM的模拟结果进行了检验,对其中不合理的部分进行了分析,指出了在干旱半干旱区陆面过程参数修正的必要性,对干旱半干旱区土壤水分传输以及大气近地面层湍流输送的参数化方案进行了改进.改进后的模式能够较好地模拟夏季连续晴天条件下沙漠的地表能量收支,因此,本文利用MLAIM研究了绿洲对其周围沙漠地表能量收支的影响,并对地表能量收支各分量之间的相互作用进行了分析.结果表明,绿洲向其下风向沙漠的水汽输送是导致其上下风向沙漠间地表能量收支差异的最重要的因子.  相似文献   

10.
土壤湿度不仅是地表水循环的重要组成部分,而且对天气和气候也有重要影响,它的模拟误差严重阻碍了人们对水循环的认知.本文首先评估了1°×1°水平分辨率的全球陆面数据同化产品(Global Land Data Assimilation System,GLDAS)对青藏高原中部那曲地区和东部玛曲地区土壤湿度的模拟性能;鉴于GLDAS较粗的分辨率无法精细描述分析区域土壤湿度空间分布特征,于是我们基于通用陆面过程模式(Community Land Surface Model,Version 4.5),开展了高分辨率0.1°×0.1°的模拟,并对高分辨率模拟土壤湿度误差的原因进行了深入分析.结果表明:(1)GLDAS陆面数据同化产品和高分辨率陆面模式模拟结果都可以反映出土壤湿度的季节变化特征,但在非冻结期均存在不同程度的干偏差,尤其是在玛曲地区;(2)对比观测和模拟的土壤湿度发现,观测数据表现出强烈的空间异质性,而模拟结果呈现出的是空间均一性.按照模拟误差进行归类分组,对比模拟性能优劣的两组站点发现,模式物理过程不是模拟性能差异的主要因素,而两组站点间地表特征参数中的土壤质地和地形参数,以及驱动数据均没有体现出空间异质性,这可能是土壤湿度模拟结果没有表现出空间异质性的原因.  相似文献   

11.
城市地表潜热通量数值模拟方法研究   总被引:2,自引:0,他引:2  
近年来,中尺度气象模式中城市建筑物影响的数值模拟方法日趋完善,城市地表潜热通量数值模拟方法是目前城市地气耦合研究面临的瓶颈问题之一.采用北京325 m气象塔140 m高度处的能量平衡观测资料,用陆面-单层城市冠层耦合模式进行了一年连续的离线模拟.根据模式对城市下垫面潜热通量模拟的系统性误差较大的特点,基于观测资料研究了城市下垫面潜热通量的数值模拟方法.采用观测分析与城市陆面模式相结合的方法,求出了城市中绿地的绿洲效应系数、城市不透水下垫面水份有效率的时间变化公式以及人为潜热排放的日变化曲线和四季最大值,可直接应用于城市陆面模式,有效提高了模式的整体模拟性能.该方法对其他城市地区有一定的参考意义.  相似文献   

12.
东亚中全新世的气候模拟及其温度变化机制探讨   总被引:20,自引:0,他引:20  
陈星  于革  刘健 《中国科学D辑》2002,32(4):335-345
大量地质证据证实了东亚和北美中全新世(6 kaBP)全年和冬夏季气温高于现代.然而,国际上PMIP计划下18个模式的古气候模拟结果未能捕捉北半球大陆中低纬地区冬季升温的气候特征.这些古气候模拟的冬季降温与地质资料揭示的冬季升温存在着巨大差异,反映出仅仅在太阳辐射变化驱动下的古气候模拟存在重大缺陷. 使用含有陆面过程的全球9层大气环流谱模式(AGCM+SSiB),采用现代植被和中全新世植被预置的不同下垫面对6 kaBP气候及其植被影响进行了模拟试验.古植被强迫下的模拟结果表明,中全新世时东亚地区各季均出现升温.尤其是模拟的冬季增温与地质资料重建的气候特征接近,反映了除太阳辐射的变化外,植被的变化对东亚地区中全新世的增温有着重要作用.该模拟结果的意义在于:(1)用具有物理机制的数值化模型并采用实际下垫面边界条件,能够较好地模拟出中全新世气候特征;(2)从动力机制的角度揭示了中全新世气候与现代气候存在巨大差别的原因在于辐射变化和下垫面植被变化;(3)中全新世下垫面植被的改变引起的地表反射率变化使得东亚陆面与西太平洋表面的热力差异随季节发生变化,因而中全新世夏季风环流增强,而冬季风环流减弱,冷空气活动受到抑制,使得中国区域冬季温度增加, 形成暖冬气候特征.  相似文献   

13.
何文英  陈洪滨  李军 《地球物理学报》1954,63(10):3573-3584
复杂多变的陆地表微波比辐射率,造成陆面上星载微波观测反演大气参数较为困难,也使得许多卫星微波资料不易同化应用到数值模式,因此迫切需要提供准确可靠的陆面微波地表比辐射率信息.随着卫星观测技术的迅速发展,利用丰富的星载被动微波观测直接反演陆面微波比辐射率成为一种主要手段.国外针对星载微波成像仪和微波垂直探测器开展较为系统的陆面微波比辐射率研究,建立不同类型的地表比辐射率反演方法,开发地表比辐射率参数化方法并应用于辐射资料同化.对于卫星观测反演陆面微波比辐射率存在的问题,开展了评估分析和方法订正.国内利用卫星观测也开展了一些陆面微波比辐射率研究工作,尚需要系统、综合的提炼.对于地表特征复杂的中国地区,还需要评估认识不同陆面微波比辐射率反演方法在我国适用情况,需要增强陆面微波比辐射率数据质量的认识以及业务应用.  相似文献   

14.
青藏高原地表温度的比较分析   总被引:11,自引:0,他引:11       下载免费PDF全文
利用青藏高原改则、甘孜、拉萨和那曲气象站的地面水银温度计观测的地面0 cm温度和地面向上、大气向下长波辐射观测资料,计算和比较分析了地面水银表温度与地表辐射温度的差异.分析表明,由于地面水银温度计接触式测温方法的局限性,使得气象台站观测场地的地面水银表温度存在较大的误差:在裸露地表,夏季白天的正偏差可达4℃以上,夜间的负偏差在-2℃左右;在稠密植被地表,白天正偏差最高可达16℃以上,夜间的负偏差可达-3℃.由于偏差的不均衡性和离散性,即使对地面水银表温度作简单平均或积分平均处理,其代表性仍然很差.在对地表温度准确性要求较高的定量化研究与应用中,应避免直接采用地面水银表温度资料,即使由于资料缺乏而以其代之,也应对其采取适当的订正补偿等处理措施.随着科学研究对地表温度准确性的更高要求,推广使用地表辐射测温的观测方法是必要的.  相似文献   

15.
以行播冬小麦为例, 深入研究了陆面目标组分温度的反演方法. 结果表明:行播作物结构模式的选择能够显著地影响组分温度的反演精度, 适当选择行播作物的结构模式对提高反演稳定性有重要意义. 认为矩阵求逆反演与迭代反演相结合是一种较为稳健的、适于陆面目标的组分温度反演方法. 模拟试验和野外地面实验验证结果都显示, 利用综合反演方法可以显著提高土壤光照面和植被顶层温度的反演精度, 并明显增强组分温度反演的稳定性. 采用的多角度多波段热红外遥感大气效应纠正只需要两个基本参数:大气透过率和大气上行辐射亮度. 如果它们能与地表组分温度协同反演, 则可实现真正意义上的精确大气纠正. 通过对ATSRII数据的反演实验, 基本上证实这一方法的可行性.  相似文献   

16.
陆面过程模式是气候模式和天气模式的核心组成部分之一.在土壤—植被—大气耦合模式(Soil-PlantAtmosphere Model,SPAM)的基础上,发展了新一代北京大学陆面过程模式PKULM(Peking University Land Model).本文首先介绍了PKULM的辐射传输、湍流输送、光合作用、土壤水热输送等过程的参数化方案;采用隐式迭代计算框架,发展并应用了一个快速的线性方程组求解算法,提高了模式计算稳定性;提出并使用了二分搜索算法计算气孔阻抗,避免了CLM(Community Land Model)等使用的迭代方法在干旱区不稳定的情况,提高了模式的适用性;采用水势为基础的土壤水分扩散方程,使模式能够模拟土壤饱和区的水分输送过程,为进一步与水文过程模式耦合奠定了基础;还发展了一个地表积水与径流过程的机理模型,提高了模式对地表水分平衡过程的模拟能力;最后,使用"中国西北干旱区陆—气相互作用观测试验"平凉站的资料对模式进行了检验并与NOAH(National Center for Environmental Prediction,Oregon State University,Air Force,and Hydrology Lab model)陆面过程模式的模拟结果进行了比较,结果表明PKULM能够较好地模拟西北半干旱区农田下垫面地气交换过程.  相似文献   

17.
我国黄土高原地区地处夏季风边缘,分布在气候和生态过渡带,气候环境的空间差异很大,对陆面能量的空间分布格局影响非常显著.然而,由于受该地区陆面过程观测站点较少的局限,对整个黄土高原区域陆面能量的空间分布规律及其影响机制的认识十分有限.在对CLM模式模拟的陆面能量平衡分量资料进行试验验证的基础上,利用CLM模式模拟的近30年黄土高原地区陆面能量平衡分量资料,分析了该地区近30年平均陆面能量平衡分量的空间变化特征以及与最干燥年和最湿润年的差异,研究了陆面能量平衡分量空间分布与经、纬度和海拔高度等地理因素及降水和气温等气候因子之间的关系.发现,黄土高原地区陆面能量平衡分量空间差异非常显著,地表净辐射和感热通量由南至北增加,潜热通量和土壤热通量从东南向西北减少;空间最干格点和最湿格点之间的地表感热、潜热和土壤热通量几乎相差1倍左右,地表能量分配由最干格点的感热通量主导转变为最湿格点的感热和潜热平分秋色;年际干湿波动对地表能量平衡分量的影响也相当显著,对感热和潜热通量的改变幅度最大接近30%.而且,经纬度和海拔高度等地理因素及温度和降水等气候要素均与陆面能量平衡分量空间分布有一定的相关性,但地表净辐射与海拔高度和纬度的关系更密切,感热通量与降水和纬度的关系更密切,而潜热和土壤热通量只与降水的关系比较密切.  相似文献   

18.
末次盛冰期东亚气候的数值模拟   总被引:16,自引:0,他引:16  
用美国国家大气研究中心(NCAR)的CCM3全球气候模式研究了末次盛冰期(LGM)和现代情景下的东亚季风和地面水分特征以及青藏高原冰川扩张. 结果表明: 在LGM时, 我国北方和西太平洋地区冬季风显著加强, 南方地区冬季风变化不大; 而对LGM时期的夏季风, 我国南方和南海地区显著减弱, 北方变化不显著; LGM时期季风的这种变化, 使我国东北、华北大部分地区、黄土高原和青藏高原东部年降水量比现代显著减少, 造成这些地区当时地面净失去更多水分, 使当地变干燥, 其中青藏高原东部、黄土高原西部地面变干燥最显著; 而在LGM时期青藏高原中部一些地区由于蒸发减少使地面变湿润, 有利于当时这些地区的湖面上升; 此外, LGM时期冬季青藏高原 绝大部分地区积雪明显比现代厚, 通过分析模拟资料计算的冰川平衡线高度发现: 尽管我们模拟出LGM时期较小的降温幅度, 但是通过模式中大气物理过程青藏高原降水和气温之间保持平衡, LGM时期当地冰川平衡线高度与现代相比降低了300~900 m, 即从现代的5400 m以上降为4600~5200 m, 指示着LGM时期青藏高原冰川的大规模扩张.  相似文献   

19.
利用中尺度模式WRFV2.2及NCEP再分析资料,以24h短期天气模拟的方式,模拟了我国2003年7月下旬江南及华南地区的高温天气,检验了模拟高温对不同陆面方案的敏感性.结果表明,模式总体上能模拟出≥35℃的气温分布形势,但模拟高温对不同陆面方案比较敏感,在模拟强度、TS评分及误差上存在较大差异,如从评分来看,SLAB方案模拟最好,RUC次之,NOAH再次之;不同陆面方案能模拟出方案之间的系统性差异,如三个方案中SLAB方案模拟平均绝对误差最小,模拟低层西太平洋副热带高压脊线最为偏东;不同的陆面方案主要通过模拟地表感热通量的变化影响地面气温变化,由不同方案引起的环流变化与地面气温变化之间的正反馈机制不重要,甚至在主要高温区存在"负"反馈.改进陆面方案对提高地面气温的数值预报能力具有重要意义.  相似文献   

20.
本文针对模式发展的需要, 在Farouki土壤热传导率参数化方案的基础上, 综合Johansen和Côté的参数化方案, 发展了一个用于青藏高原中部的土壤热传导率参数化方案, 用“全球协调加强观测计划之亚澳季风青藏高原试验计划(CAMP/Tibet)”中那曲布交(BJ)站实际资料对该方案进行了检验, 并将它用于公共陆面模式(CoLM)中, 对青藏高原那曲地区进行了单点数值模拟试验. 结果表明: 在未冻结及冻结土壤中, 新方案比Farouki方案计算的土壤热传导率小, 更接近实测值. 加入新方案的CoLM模式对土壤温度模拟的准确性比原模式有一定的提高.  相似文献   

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

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