首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 750 毫秒
1.
Summary A new physical inversion-based algorithm for retrieving rain rate over the ocean with the Special Sensor Microwave Imager (SSM/I) is described. In a departure from other rain rate retrieval algorithms, the satellite observables inverted in the present algorithm are not the raw brightness temperatures but rather normalized polarizations for 19.35, 37.0, and 85.5 GHz, plus an 85.5 GHz scattering index which is sensitive primarily to ice particles aloft. The normalized polarizations are interpreted as footprint-averages of theoretically derived analytic functions of the column optical depth associated primarily with liquid water. The effective vertical depth of the rain layer is specified as a function of the SSM/I estimated column water vapor.The retrieval algorithm performs an iterative search for a high resolution (12.5 km) rain field which is simultaneously consistent with the 19.35 and 37.0 GHz normalized polarizations. The first-guess rain rate field is supplied by the 85.5 GHz scattering index. At gridpoints for which the rain column optical depth exceeds the dynamic range of the attenuation-based indices, the first-guess field is left essentially unmodified; elsewhere, the required consistency with the 19 and 37 GHz indices usually results in significant modification of the scattering-based rain rate estimates.The algorithm as described here is a prototype implementation which was developed with reference only to idealized theoretical models; empirical improvements to the numerical scheme and the model coefficients will be made in the future as results from the first Precipitation [algorithm] Intercomparison Project 1 (PIP-1) and the second phase of the GPCP (Global Precipitation Climatology Project) algorithm Intercomparison Project (AIP/2) are analyzed, as well as data from individual validation efforts. Although the algorithm is physically based and uses all SSM/I channels, it is computationally much less demanding than cloud/radiative transfer model-based inversion algorithms published else-where.With 9 Figures  相似文献   

2.
Using synthetic geometrical clouds and radiative microwave model, we examine the possibility to correct the estimations of liquid water path (LWP) or rain rate with cloud cover measurement. This information may be gotten by co-localized measurements of microwave and infrared/visible measurements on new satellites (TRMM, ADEOS 2, …). In a first step, the effects of fractional cloud cover on microwave brightness temperatures (TB) are investigated in three typical cases of nonprecipitating and precipitating (stratiform and convective) clouds. The beam-filling error (BFE) on brightness temperatures may be analyzed with the known spatial variability using 1D or 3D radiative transfer model. Relationships between BFE and subpixel cloud fraction (CF) are discussed according to the cloud type. We tested several parameters that characterize the horizontal cloud inhomogeneity within a radiometer field of view. BFE was found very sensitive to cloud type and inhomogeneity and is maximum for raining cloud with open spatial structure. In order to account for the uncertainty introduced by the spatial distribution, dependence of BFE on textural-based parameters is also discussed using homogeneity, entropy and an indicator of CF horizontal gradient.  相似文献   

3.
对流性降水云微波辐射特性   总被引:2,自引:1,他引:1       下载免费PDF全文
结合MM5模式和三维微波辐射传输模式, 对2003年7月9日宜昌地区一次典型的中尺度降水中心的对流性降水云微波辐射特性进行研究。结果表明:MM5模式模拟的降水量和落区与实况一致,模拟的水凝物廓线也与TMI反演值接近,85.5 GHz通道辐射亮温与TMI实测情况相近。85.5 GHz通道亮温与地面雨强相关性很弱, 受云中云冰和雪花的散射降温作用显著, 由于其他粒子的综合作用以及斜角观测造成的位置偏移, 霰粒子对该通道亮温散射作用不明显。19.35 GHz通道亮温随雨强增加先升温后降温; 与霰粒子含量表现出明显的负相关关系。37.0 GHz通道亮温随雨强的增加而降温, 雨强大于20 mm/h后达到饱和, 主要受雨水降温作用影响。倾斜观测比天顶垂直观测产生更低的亮温低值中心, 且频率越高, 低值中心的偏移越严重。  相似文献   

4.
Summary In this paper a retrieval technique for estimating rainfall rates is introduced. The novel feature of this technique is the combination of two satellite radiometers — the Special Sensor Microwave/Imager (SSM/I) and the Advanced Very High-Resolution Radiometer (AVHRR) — with mesoscale weather prediction model data. This offers an adjustment of the model atmospheres to reality which is necessary for calculating brightness temperatures that can be compared with microwave satellite measurements.In sensitivity studies it was found that the estimation of precipitation is determined to a high degree by the particle size distribution of rain and snow, especially by the size distribution of solid hydrometeors. These studies also reveal the influence of the knowledge of the correct cloud coverage inside a SSM/I pixel and the importance of using a realistic temperature profile instead of using standard atmospheres.The retrieval technique is based on radiative transfer calculations using the model of Kummerow et al. (1989). The algorithm consists of two parts: First Guess (FG) brightness temperatures for the SSM/I frequencies are generated as a function of the cloud top height and the cloud coverage, derived from AVHRR data and predictions from a meso-scale model. The rainfall rate of different types of clouds containing raindrops, ice particles and coexisting ice and water hydrometeors is then calculated as a function of the cloud top height. As an example, a strong convective rain event over the western part of Europe and over the Alps is taken to evaluate the performance of this technique. Good agreement with radar data from the German Weather Service was achieved. Compared to statistical rainfall algorithms, the current algorithm shows a better performance of detecting rainfall areas.With 12 Figures  相似文献   

5.
微波成像仪通道对降水云参数响应的数值模拟研究   总被引:3,自引:2,他引:1  
王雨  傅云飞 《气象学报》2010,68(3):315-324
由于降水云类型和结构等的差异,将导致微波信号出现不同的响应,因此基于星载微波仪器观测结果对降水云参数进行反演仍是目前国际上的难题之一。为了解层云和对流云等两类降水云的TMI微波信号特点,本文以相应的结构模型为基础,利用微波辐射传输模式MWRT分别模拟了不同下垫面和降水云参数条件下的TMI各通道微波亮温,并分析了各种组合通道信号,包括反映发射信息的极化差D(洋面)及反映散射信息的高频极化校正亮温PCT85(洋面)和高低频率垂直极化差VFD(陆面)的变化情况,以此从理论上探讨可能的降水反演方法。研究结果表明:在洋面,D和PCT85均随降水率的升高而减小,其中D_(10)和PCT85都能很好地反映降水率的变化,但分别受地表状况和云体性质(冻结层高度和霰粒大小)的影响较大;D_(19)和D_(37)虽然对下垫面状况和冻结层高度都不甚敏感,但随降水率的增加,存在极化丧失现象,因此不适合反演降水。在陆面,较高频取85 GHz的VFD与降水率有很好的对应关系,随降水率的增长而增大,其中VFD37-85受地表状况影响最小,但受云体性质影响较大;37 GHz对液态水变化敏感,其发射效应易混淆散射信号,故较高频取37 GHz的VFD10-37和VFD19-37不适合研究降水与微波信号之间的关系;层云与对流云存在有无霰粒的差异,会对微波散射性质造成影响,导致PCT85和VED对层云降水率的响应要高于对流云。  相似文献   

6.
RTTOV和CRTM对“罗莎”台风卫星微波观测的模拟研究与比较   总被引:3,自引:1,他引:2  
快速辐射传输模式中水成物辐射效应的考虑对卫星观测模拟具有重要的影响。利用区域中尺度数值模式WRF预报输出水成物含量和温湿廓线等参数,使用欧洲中期数值天气预报中心和美国卫星资料同化联合中心发展建立的快速辐射传输模式RTTOV和CRTM,结合雅可比模式计算得到的响应函数,分析了各种水成物辐射效应对AMSUA/B卫星微波观测各通道亮温模拟的影响,并比较了RTTOV和CRTM结果间的异同。结果表明,RTTOV和CRTM两个快速辐射传输模式计算得到的水成物响应函数和水成物辐射效应对模拟亮温影响特征基本一致,但量级上存在差异。雅可比表征的水成物响应函数和水成物本身的分布层次对应关系良好,云水和雨水的响应函数主要分布在中下层,冰水和雪则分布在中上层。相对于RTTOV,CRTM计算得到的水成物响应函数要大一些。对AMSUA,云水和雨水的辐射效应对通道模拟亮温主要以增温为主,尤其是在低频的窗区通道,冰晶和雪则体现散射效应为主的降温作用。对AMSUB,水成物辐射效应的考虑使得模拟亮温降低,尤其是通道2降温幅度最大。CRTM计算出的水成物辐射效应带来的模拟亮温变化幅度大于RTTOV。目前,CRTM除了考虑云、雨、冰和雪4种水成物外,还考虑霰和雹的辐射效应,对水成物辐射效应的考虑比RTTOV更完善一些。  相似文献   

7.
云参数对微波亮温模拟计算的影响试验   总被引:1,自引:0,他引:1       下载免费PDF全文
利用CRTM (Community Radiative Transfer Model) 快速辐射传输模式对NOAA-K系列卫星的AMSU-A通道亮温进行正演模拟,重点研究云粒子类型、云高、云厚度等云参数对微波亮温模拟的影响。结果表明:改变云粒子类型时,云水和雨水对模拟亮温影响较大,模拟亮温值比晴空高1 K;霰、雪、冰、雹等固态粒子对模拟亮温的影响较小,模拟亮温值略低于晴空无云情况;云层光学厚度较大时,各通道亮温受云层影响的情况取决于权重函数峰值高度和云顶高度的配置;多个高度存在云时,若最上层云较厚 (2 km),光学厚度大,相应通道亮温取决于最上层云,较低层云对亮温不产生影响;云层变薄,光学厚度减小,高度低于云层或略高于云顶的通道亮温随云层厚度的变化明显,若通道高度远高于云顶,云层厚度的变化对于其亮温模拟的影响很小。  相似文献   

8.
Cloud radiative and microphysical effects on the relation between spatial mean rain rate, rain intensity and fractional rainfall coverage are investigated in this study by conducting and analyzing a series of two-dimensional cloud resolving model sensitivity experiments of pre-summer torrential rainfall in June 2008. The analysis of time-mean data shows that the exclusion of radiative effects of liquid clouds reduces domain mean rain rate by decreasing convective rain rate mainly through the reduced convective-rainfall area associated with the strengthened hydrometeor gain in the presence of radiative effects of ice clouds, whereas it increases domain mean rain rate by enhancing convective rain rate mainly via the intensified convective rain intensity associated with the enhanced net condensation in the absence of radiative effects of ice clouds. The removal of radiative effects of ice clouds decreases domain mean rain rate by reducing stratiform rain rate through the suppressed stratiform rain intensity related to the suppressed net condensation in the presence of radiative effects of liquid clouds, whereas it increases domain mean rain rate by strengthening convective rain rate mainly via the enhanced convective rain intensity in response to the enhanced net condensation in the absence of radiative effects of liquid clouds. The elimination of microphysical effects of ice clouds suppresses domain mean rain rate by reducing stratiform rain rate through the reduced stratiform-rainfall area associated with severely reduced hydrometeor loss.  相似文献   

9.
运用中尺度数值模式WRF3.5.1对2011年1月1日贵州境内的一次冻雨天气过程进行了数值模拟,研究了本次过程的大气层结、冻雨区云系的宏微观结构和云物理特征,初步分析了冻雨形成的云微物理过程和成因。结果表明,贵州境内的冻雨区(26°N~29°N)具有冷性和部分"冷—暖—冷"的温度层结,在高层没有显著的冰相粒子,冻雨区是相对较强的水汽辐合中心,丰富的水汽输送在冷性的环境条件下形成云滴,进而碰并产生雨滴,过冷雨水主要通过暖雨过程形成;雨滴继续下落至近地层并保持过冷雨水形式,最后接触到低于0°C的物体或地面,迅速冻结而产生地面冻雨。  相似文献   

10.
Cloud and precipitation parameterization schemes are evaluated, and their sensitivity to the method and/or parameters used to determine cloud physical processes is examined using a singlecolumn version of the Unified Model (SCUM). In the experiment for TWP-ICE, cloud fraction is overestimated (underestimated) in the upper (lower) troposphere due to the wet (dry) bias. The precipitation rate is well simulated during the active monsoon period, but overestimated during the suppressed monsoon and clear skies periods. In the moist convection scheme, trigger condition and entrainment process affect the lower tropospheric humidity through the impact on convective occurrence frequency and intensity, respectively. Strengthening the trigger condition and using the adaptive entrainment method alleviate the low-level dry bias. In the microphysics scheme, more large-scale precipitation is produced with prognostic rain, due to rain sedimentation considering vertical velocity of rain drop, than with diagnostic rain. Less ice/snow deposition with the prognostic two-ice category results in lower ice water content and upper-level cloud fraction than with the diagnostic splitting method for the twoice category. In the cloud macrophysics scheme, the prognostic cloud fraction and cloud/ice water content scheme produces a larger cloud fraction and more cloud/ice water content than the diagnostic scheme, mainly due to detrainment from moist convection (cloud source) that surpasses the effect of convective heating and drying (cloud sink). This affects temperature by influencing the radiative, convective, and microphysical processes. The experiment with combined modifications in cloud and precipitation schemes shows that interaction between modified moist convection and cloud macrophysics schemes results in more alleviation of the cold bias not only at the lower levels but also at the upper levels.  相似文献   

11.
通过冰雹云模式模拟的一次冰雹云降水过程中降水粒子廓线和微波辐射传输模式结合,分析了冰雹云发展的不同阶段的微物理含量垂直结构变化及其对微波亮温的影响,得到以下几点结论:1)如果微波通道受到降水粒子散射和辐射的共同作用,如降水云早期的85 GHz亮温,成熟期的19 GHz亮温及消散期的37 GHz亮温,由于辐射和散射信息互相抵消,致使亮温随雨强的变化较复杂,这些通道亮温和雨强的相关性明显降低,不宜被用来反演地面雨强。2)根据19 GHz亮温随地面雨强或冰相粒子柱含量的改变,可以大致确定降雨云的不同阶段:在发展阶段,主要是降雨层以上的冰相粒子,尤其霰粒影响19 GHz亮温,致使其亮温与冰相粒子柱含量具有较好的负相关,而与地面雨强相关性较差;在成熟阶段,主要受雨水上层逐渐增加的辐射和冰相粒子散射共同作用,使得19GHz亮温与地面雨强和冰相粒子柱含量的相关性都不太好;在消散阶段,19 GHz亮温主要受较强的雨水辐射影响,与地面雨强和冰相粒子柱含量均有着较高的正相关。3)37 GHz是相对比较稳定的通道,其亮温与地面雨强有较好的线性关系,尤其与冰相粒子柱含量相关性更好,因此是反演地面雨强和冰相粒子柱含量的最佳通道。85 GHz亮温对降雨云体的中高层结构较为敏感,使得其亮温随地面雨强增加而降低的变化比较离散,不如37 GHz的集中。  相似文献   

12.
Accurate information on atmospheric temperature of tropical cyclones (TCs) is important for monitoring and prediction of their developments and evolution. For hurricanes, temperature anomaly in the upper troposphere can be derived from Advanced Microwave Sounding Unit (AMSU) and Advanced Technology Microwave Sounder (ATMS) through either regression-based or variational retrieval algorithms. This study investigates the dependency of TC warm core structure on emission and scattering processes in the forward operator used for radiance computations in temperature retrievals. In particular, the precipitation scattering at ATMS high-frequency channels can significantly change the retrieval outcomes. The simulation results in this study reveal that the brightness temperatures at 183 GHz could be depressed by 30–50 K under cloud ice water path of 1.5 mm, and thus, the temperature structure in hurricane atmosphere could be distorted if the ice cloud scattering was inaccurately characterized in the retrieval system. It is found that for Hurricanes Irma, Maria, and Harvey that occurred in 2017, their warm core anomalies retrieved from ATMS temperature sounding channels 4–15 were more reasonable and realistic, compared with the retrievals from all other channel combinations and earlier hurricane simulation results.  相似文献   

13.
我国新一代极轨业务气象卫星风云三号 (02) 批计划2012年发射。该文利用UWNMS模拟2005年Katrina飓风的结果作为基础数据集,借助VDISORT微波辐射传输模式对风云三号 (02) 批计划装载的微波探测仪器中50~60 GHz和新增的118.75 GHz频点的降水特性进行初步研究。首先通过晴空权重函数匹配,选择出50~60 GHz与118.75 GHz频点匹配关系较好的4对通道。敏感性分析表明:各通道对各种水凝物粒子均很敏感,可用于改进现有业务降水反演算法。分别选取50~60 GHz 4个通道、118.75 GHz 4个通道、50~60 GHz及118.75 GHz全部通道3种不同的通道组合进行反演试验。结果表明:将50~60 GHz及118.75 GHz通道联合起来进行降水反演可提高降水反演的精度,并可以更好地区分降水区与非降水区。  相似文献   

14.
A vertically inhomogeneous mid-latitude mixed-phase altocumulus cloud was observed around 17:26 UTC on Oct. 14, 2001 during the 9th Cloud Layer Experiment (CLEX9). In this study the microphysical and optical properties of this cloud are investigated on the basis of in-situ observed vertical profiles of particle size and habit distributions. Two cloud models, assuming that the cloud properties were vertically homogeneous and inhomogeneous, are adopted to derive the bulk optical and radiative properties of this cloud. The observed microphysical properties are combined with the theoretical solutions to the scattering and absorption properties of individual cloud particles to determine the bulk optical properties at various heights within the cloud layer. The single-scattering properties of spherical liquid water droplets and nonspherical ice crystals are obtained from the Lorenz–Mie theory and an existing database, respectively. The bulk microphysical and optical properties associated with the inhomogeneous model depend strongly on the height above the cloud-base whereas the dependence is smoothed out in the case of the homogeneous model. Furthermore, the transfer of infrared radiation is simulated in conjunction with the two cloud models. It is shown that the brightness temperatures at the top of the atmosphere in the case of the homogeneous model can be 1.5% (3.8 K) higher than their counterpart associated with the inhomogeneous cloud model. This result demonstrates that the effect of the vertical inhomogeneity of a mixed-phase cloud on its radiative properties is not negligible.  相似文献   

15.
Direct assimilation of cloud-affected microwave brightness temperatures from AMSU-A into the GSI three-dimensional variational(3D-Var) assimilation system is preliminarily studied in this paper. A combination of cloud microphysics parameters retrieved by the 1D-Var algorithm(including vertical profiles of cloud liquid water content, ice water content, and rain water content) and atmospheric state parameters from objective analysis fields of an NWP model are used as background fields. Three cloud microphysics parameters(cloud liquid water content, ice water content, and rain water content) are applied to the control variable. Typhoon Halong(2014) is selected as an example. The results show that direct assimilation of cloud-affected AMSU-A observations can effectively adjust the structure of large-scale temperature, humidity and wind analysis fields due to the assimilation of more AMSU-A observations in typhoon cloudy areas, especially typhoon spiral cloud belts. These adjustments, with temperatures increasing and humidities decreasing in the movement direction of the typhoon,bring the forecasted typhoon moving direction closer to its real path. The assimilation of cloud-affected satellite microwave brightness temperatures can provide better analysis fields that are more similar to the actual situation. Furthermore, typhoon prediction accuracy is improved using these assimilation analysis fields as the initial forecast fields in NWP models.  相似文献   

16.
Presented is a review of the radiative properties of ice clouds from three perspectives: light scattering simulations, remote sensing applications, and broadband radiation parameterizations appropriate for numerical models. On the subject of light scattering simulations, several classical computational approaches are reviewed, including the conventional geometric-optics method and its improved forms, the finite-difference time domain technique, the pseudo-spectral time domain technique, the discrete dipole approximation method, and the T-matrix method, with specific applications to the computation of the singlescattering properties of individual ice crystals. The strengths and weaknesses associated with each approach are discussed.With reference to remote sensing, operational retrieval algorithms are reviewed for retrieving cloud optical depth and effective particle size based on solar or thermal infrared(IR) bands. To illustrate the performance of the current solar- and IR-based retrievals, two case studies are presented based on spaceborne observations. The need for a more realistic ice cloud optical model to obtain spectrally consistent retrievals is demonstrated. Furthermore, to complement ice cloud property studies based on passive radiometric measurements, the advantage of incorporating lidar and/or polarimetric measurements is discussed.The performance of ice cloud models based on the use of different ice habits to represent ice particles is illustrated by comparing model results with satellite observations. A summary is provided of a number of parameterization schemes for ice cloud radiative properties that were developed for application to broadband radiative transfer submodels within general circulation models(GCMs). The availability of the single-scattering properties of complex ice habits has led to more accurate radiation parameterizations. In conclusion, the importance of using nonspherical ice particle models in GCM simulations for climate studies is proven.  相似文献   

17.
我国冬季冻雨和冰粒天气的形成机制及预报着眼点   总被引:14,自引:4,他引:10  
漆梁波 《气象》2012,38(7):769-778
利用探空和地面观测资料,通过对2001年冬季至2010年冬季我国不同区域(分为4个区域:北方、江南、华南、西南)的冻雨和冰粒天气形成的物理过程进行分析发现:(1)除北方区域外,我国其他区域的冻雨主要以暖雨机制为主。北方区域的融化类冻雨比例也仅为39%,但纬度越高,出现融化类冻雨的几率高于上述比例。暖层出现是冻雨天气的重要特征,但暖层作用主要是输送水汽和维持锋面系统,以保证降水的发生和持续,低层及地面气温普遍低于0℃可能是最重要的原因。(2)我国冰粒天气的形成机制主要以融化机制为主。冰粒天气的云顶高度普遍高于冻雨天气。冰粒天气的暖层厚度和强度均小于冻雨天气,这主要是由于弱暖层只是部分融化冰晶和雪花,使其重新冻结成为可能。冰粒天气的700hPa风速值普遍小于冻雨天气,这一方面说明冰粒天气对水汽输送条件要低一些,另一方面也反映了冰粒天气暖层较弱的特点。(3)云顶高度、暖层强度和厚度、低层冷层温度露点差、700hPa风速以及地面气温是甄别冻雨和冰粒天气的特征量,但不同区域,这些特征量的有效性不一样。西南区域冻雨和冰粒天气的主要差别在地面气温,其他特征量或差别不明显,或代表性不足,只可以作为辅助判断的因子。  相似文献   

18.
陈宝君  肖辉 《大气科学》2007,31(2):273-290
利用中国科学院大气物理研究所开发的三维全弹性冰雹云模式,对美国对流降水协作试验(CCOPE)期间观测的1981年8月1日雹云进行模拟,讨论在过冷雨水低含量条件下冰雹形成和增长机制及其碘化银催化效果。结果表明:(1) 自然云的模拟与观测事实一致,如最大上升气流速度、云顶高度、流场结构以及雹胚组成等。(2) 雹胚以霰为主,霰主要来自冰雪晶与过冷小水滴的碰冻,其次来自雪的积聚转化;霰、冻滴和冰雹在形成后主要靠碰并过冷云水增长。(3)人工催化试验表明,碘化银主要以凝华核(包括凝结-冻结)的作用产生大量的人工冰晶,加速了过冷水向冰晶的转化,过冷云水因而大量减少;催化后霰和冻滴的数浓度增大,对过冷云水的竞争增强,其平均尺度减小导致转化成雹的数量减少;冰雹碰冻过冷云水的增长在催化后也被削弱,导致冰雹总质量进一步减少。此外,催化后降雨量也显著减少。  相似文献   

19.
Aerosols emitted from volcanic activities and polluted mid-latitudes regions are efficiently transported over the Arctic during winter by the large-scale atmospheric circulation. These aerosols are highly acidic. The acid coating on ice nuclei, which are present among these aerosols, alters their ability to nucleate ice crystals. In this research, the effect of acid coating on deposition and contact ice nuclei on the Arctic cloud and radiation is evaluated for January 2007 using a regional climate model. Results show that the suppression of contact freezing by acid coating on ice nuclei leads to small changes of the cloud microstructure and has no significant effect on the cloud radiative forcing (CRF) at the top of the atmosphere when compared with the effect of the alteration of deposition ice nucleation by acid coating on deposition ice nuclei. There is a negative feedback by which the suppression of contact freezing leads to an increase of the ice crystal nucleation rate by deposition ice nucleation. As a result, the suppression of contact freezing leads to an increase of the cloud ice crystal concentration. Changes in the cloud liquid and ice water contents remain small and the CRF is not significantly modified. The alteration of deposition ice nucleation by acid coating on ice nuclei is dominant over the alteration of contact freezing.  相似文献   

20.
基于TMI产品资料对数值模式水凝物模拟能力的检验分析   总被引:4,自引:1,他引:3  
本文利用热带测雨卫星(Tropical Rainfall Measuring Mission,TRMM)搭载微波成像仪(TRMM Microwave Imager,TMI)的探测及反演结果,结合微波辐射传输模式,就2004年17号台风暹芭(Chaba)过程,对AREM(Advanced Regional Eta-coo...  相似文献   

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

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