首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Monsoon depressions, that form during the Indian summer monsoon season (June to September) are known to be baroclinic disturbances (horizontal scale 2,000–3,000 km) and are driven by deep convection that carries a very large vertical slope towards cold air aloft in the upper troposphere. Deep convection is nearly always organized around the scale of these depressions. In the maintenance of the monsoon depression the generation of eddy kinetic energy on the scale of the monsoon depression is largely governed by the “in scale” covariance of heating and temperature and of vertical velocity and temperature over the region of the monsoon depression. There are normally about 6–8 monsoon depressions during a summer monsoon season. Recent years 2009, 2010 and 2011 saw very few (around 1, 0 and 1 per season respectively). The best numerical models such as those from ECMWF and US (GFS) carried many false alarms in their 3–5 day forecasts, more like 6–8 disturbances. Even in recent years with fewer observed monsoon depressions a much larger number of depressions is noted in ECMWF forecasts. These are fairly comprehensive models that carry vast data sets (surface and satellite based), detailed data assimilation, and are run at very high resolutions. The monsoon depression is well resolved by these respective horizontal resolutions in these models (at 15 and 35 km). These models carry complete and detailed physical parameterizations. The false alarms in their forecasts leads us to suggest that some additional important ingredient may be missing in these current best state of the art models. This paper addresses the effects of pollution for the enhancement of cloud condensation nuclei and the resulting disruption of the organization of convection in monsoon depressions. Our specific studies make use of a high resolution mesoscale model (WRF/CHEM) to explore the impacts of the first and second aerosol indirect effects proposed by Twomey and Albrecht. We have conducted preliminary studies including examination of the evolution of radar reflectivity (computed inversely from the model hydrometeors) for normal and enhanced CCN effects (arising from enhanced monsoon pollution). The time lapse histories show a major disruption in the organization of convection of the monsoon depressions on the time scale of a week to 10 days in these enhanced CCN scenarios.  相似文献   

2.
Performance of national centers for environmental prediction based global forecast system (GFS) T574/L64 and GFS T382/L64 over Indian region has been evaluated for the summer monsoon season of 2011. The real-time model outputs are generated daily at India Meteorological Department, New Delhi for the forecasts up to 7 days. Verification of rainfall forecasts has been carried out against observed rainfall analysis. Performance of the model is also examined in terms of lower tropospheric wind circulation, vertical structure of specific humidity and precipitable water content. Case study of a monsoon depression is also illustrated. Results obtained show that, in general, both the GFS T382 and T574 forecasts are skillful to capture climatologically heavy rainfall regions. However, the accuracy in prediction of location and magnitude of rainfall fluctuates considerably. The verification results, at the spatial scale of 50 km resolution, in a regional spatial scale and country as a whole, in terms of continuous skill score, time series and categorical statistics, have demonstrated superiority of GFS T574 against T382 over Indian region. Both the model shows bias of lower tropospheric drying and upper tropospheric moistening. A bias of anti-cyclonic circulation in the lower tropospheric level lay over the central India, where rainfall as well as precipitable water content shows negative bias. Considerable differences between GFS T574 and T382 are noticed in the structure of model bias in terms of lower tropospheric wind circulation, vertical structure of specific humidity and precipitable water contents. The magnitude of error for these parameters increases with forecast lead time in both GFS T574 and T382. The results documented are expected to be useful to the forecasters, monsoon researchers and modeling community.  相似文献   

3.
Summer and winter climates simulated with the ECMWF (cycle 33) model at spectral scales T21, T42, T63 and T106 are analyzed to determine the impact of changes in horizontal resolution on atmospheric water vapor, clouds, convection, and precipitation. Qualitative changes in many moist processes occur in the transition from T21 to T42, especially in the tropics; at higher resolutions mostly incremental variations from patterns established at T42 result. Large-scale tropical moist processes are simulated more realistically at T21 than at finer resolutions, possibly reflecting a mismatch between the finer-scale dynamics and the scales at which the underlying assumptions of the physical parameterizations apply. Global precipitation increases monotonically with resolution, as a consequence of increasing convection. Global cloud cover, however, decreases in the transition from T21 to T42 due to drying of the tropics, but then increases slightly at finer resolutions. These small global increases are an outcome of compensating changes in different regions: decreases in cloud cover due to drying of the atmosphere at low latitudes are offset by high-latitude increases resulting from enhanced relative humidity associated with an intensifying atmospheric cold bias at finer resolutions.  相似文献   

4.
A low pressure system that formed on 21 September 2006 over eastern India/Bay of Bengal intensified into a monsoon depression resulting in copious rainfall over north-eastern and central parts of India. Four numerical experiments are performed to examine the performance of assimilation schemes in simulating this monsoon depression using the Fifth Generation Mesoscale Model (MM5). Forecasts from a base simulation (with no data assimilation), a four-dimensional data assimilation (FDDA) system, a simple surface data assimilation (SDA) system coupled with FDDA, and a flux-adjusting SDA system (FASDAS) coupled with FDDA are compared with each other and with observations. The model is initialized with Global Forecast System (GFS) forecast fields starting from 19 September 2006, with assimilation being done for the first 24 hours using conventional observations, sounding and surface data of temperature and moisture from Advanced TIROS Operational Vertical Sounder satellite and surface wind data over the ocean from QuikSCAT. Forecasts are then made from these assimilated states. In general, results indicate that the FASDAS forecast provides more realistic prognostic fields as compared to the other three forecasts. When compared with other forecasts, results indicate that the FASDAS forecast yielded lower root-mean-square (r.m.s.) errors for the pressure field and improved simulations of surface/near-surface temperature, moisture, sensible and latent heat fluxes, and potential vorticity. Heat and moisture budget analyses to assess the simulation of convection revealed that the two forecasts with the surface data assimilation (SDA and FASDAS) are superior to the base and FDDA forecasts. An important conclusion is that, even though monsoon depressions are large synoptic systems, mesoscale features including rainfall are affected by surface processes. Enhanced representation of land-surface processes provides a significant improvement in the model performance even under active monsoon conditions where the synoptic forcings are expected to be dominant.  相似文献   

5.
Summary  This study explores the nowcasting and short-range forecasting (up to 3 days) skills of rainfall over the tropics using a high resolution global model. Since the model-predicted rainfall is very sensitive to model parameters, four key model parameters were first selected. They are the Asselin filter coefficient, the fourth order horizontal diffusion coefficient, the surface moisture flux coefficient, and the vertical diffusion coefficient. The optimal values were defined as those which contributed to the best one day rainfall forecasts in the present study. In order to demonstrate and improve the precipitation forecast skill, several numerical experiments were designed using the 14-level Florida State University Global Spectral Model (FSUGSM) at a resolution of T106. Comparisons were also made of the short-range forecasts obtained from a control experiment subjected to normal mode initialization (NMI) versus experiments based on physical initialization (PI). The latter experiments were integrated using the original FSUGSM and a modified version. This modified FSUGSM was developed here by applying a reverse cumulus parameterization alorithm to the regular forecast model, which restructures the vertical humidity distribution and constrains the large-scale model’s moisture error growth during the model integration. An improved short-range rainfall prediction skill was achieved from the modified FSUGSM in this study. The results showed a better agreement between model-based and observed rainfall intensity and pattern. Received January 18, 1999  相似文献   

6.
采用中尺度数值模式MM5V2对1998年6月下旬发生在长江流域持续的暴雨过程进行分析研究。通过尺度分离与数值模拟对比试验，着重分析了暴雨过程中天气尺度与次天气尺度各物理量场的结构特征，提出本次暴雨过程形成的物理机制:天气尺度流场与水汽场为降水提供持续的远距离水汽输送通道，次天气尺度流场形成稳定的经向强辐合，为水汽的抬升与凝结提供动力条件；在有利的高、低空急流的配置下，暴雨区落在高空急流轴以南、低空急流轴以北；次天气尺度温度场下暖上冷的热力不稳定层结促进了热力不稳定的发展，促使暴雨增幅；特大暴雨发生地区上空的次天气尺度湿度的高值中心，有利于湿空气在上升运动中释放潜热，形成暴雨的反馈机制。数值试验分离模式初始场不同尺度系统信息，揭示了不同尺度系统在暴雨发生过程中的动力作用，没有中尺度系统的配合，仅有天气尺度系统信息，或只有次天气尺度系统信息，没有大尺度系统的配合，暴雨的强度及范围都将有所消减。分析及数值试验结果表明大暴雨是在天气尺度和次天气尺度系统的共同作用下才得以产生和维持的。  相似文献   

7.
Summary This study entails a reanalysis of the TOGA-COARE, CEPEX period covering the months December 1992 through March 1993. Four times daily data sets were reanalyzed for the period. The unique aspect of this reanalysis is that (i) it incorporates physical initialization i.e. assimilation of rainfall from raingauge and satellite (OLR and microwave radiometric data), (ii) this reanalysis is carried out at a very high resolution, T170, global spectral model (space resolution of roughly 70 km at equator) and (iii) ECMWF reanalysis is used as a first guess field for this study.This analysis is global and is carried out at 15 vertical levels using the FSU global spectral model. In this study we have analyzed a large number of meterological variables such as wind, temperature, humidity, surface pressure, vertical velocity, cloud, rainfall, surface fluxes and diabatic heating. In this part of the paper we present the monthly mean and samples of daily fields of the reanalysis. Some of the major results of the reanalysis include the relationships among the flow fields and the ITCZ convection; the shallow northeast monsoon current and its rapid turning towards westerlies with height; an overabundance of middle clouds over the region of the ITCZ; and the distribution of heat sources and sinks. A motivation for this reanalysis is the forthcoming Indian Ocean Experiment (INDOEX), which is an international field experiment to be conducted during the winter monsoon time frame in 1999.With 25 Figures  相似文献   

8.
１．ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｉｎｉｔｉａｔｉｏｎｏｆｔｈｅｃｒｏｓｅｑｕａｔｏｒｉａｌｆｌｏｗｏｆｔｈｅＳｏｍａｌｉａｃｏａｓｔｄｕｒｉｎｇＭａｙｉｎｒｅｓｐｏｎｓｅｔｏｔｈｅｈｅａｔｉｎｇｏｖｅｒｔｈｅＳｏｕｔｈＡｓｉａｎｃｏｎｔｉｎｅｎｔｍ...  相似文献   

9.
During the summer monsoon (1 June to 30 September) 2007, real-time district level rainfall forecasts in short-range time scale were generated for Indian region applying multimodel ensemble technique. The pre-assigned grid point weights on the basis of correlation coefficients (CC) between the observed values and forecast values are determined for each constituent model at the resolution of 0.5° × 0.5° utilizing two seasons datasets (1 June to 30 September, 2005 and 2006), and the multimodel ensemble forecasts (day 1 and day 2 forecasts) are generated at the same resolution on a real-time basis. The ensemble forecast fields are then used to prepare forecasts for each district taking the average value of all grid points falling in a particular district. In this paper we examined the performance skill of the multimodel ensemble-based real-time district level short-range forecast of rainfall. It has clearly emerged from the results that the multimodel ensemble technique reported in this study is superior to each ensemble member. District wise performance of the ensemble rainfall forecast reveals that the technique, in general, is capable of providing reasonably good forecast skill over most districts of the country, particularly over the districts where the monsoon systems are dominant. Though the procedure shows appreciable skill to predict occurrence or non-occurrence of rainfall at the district level, it always underestimates rainfall amount, particularly in heavy rainfall events. Possible reasons of this failure may be due to model bias and poor data assimilation procedure.  相似文献   

10.
To investigate the impact of increasing horizontal resolution on a simulated model climate, we conducted an experiment using the Japan Meteorological Agency (JMA) operational global atmosphere model (JMA-GSM0103). The models with four different horizontal resolutions ranging from T42 to T213 have been integrated over three years with prescribed climate sea surface temperature in the experiment. The distributions of 3-year averaged seasonal mean fields are basically similar among the models with different resolution, although there are some monotonic and systematic differences with increasing resolution. However, the climatology of synoptic scale phenomena is well represented in higher resolution models. The position and amount of precipitation in Baiu front (or Mei-yu) at higher resolution models agree well with observations. The start time of northward propagation of heavy precipitation over the Bay of Bengal, which is associated with Indian monsoon development, is also well simulated in higher resolution models. The number of tropical cyclones increases monotonically with resolutions. The simulated tropical cyclones become more realistic with increasing resolution.  相似文献   

11.
The large-scale and small-scale errors could affect background error covariances for a regional numerical model with the specified grid resolution.Based on the different background error covariances influenced by different scale errors,this study tries to construct a so-called"optimal background error covariances"to consider the interactions among different scale errors.For this purpose,a linear combination of the forecast differences influenced by information of errors at different scales is used to construct the new forecast differences for estimating optimal background error covariances.By adjusting the relative weight of the forecast differences influenced by information of smaller-scale errors,the relative influence of different scale errors on optimal background error covariances can be changed.For a heavy rainfall case,the corresponding optimal background error covariances can be estimated through choosing proper weighting factor for forecast differences influenced by information of smaller-scale errors.The data assimilation and forecast with these optimal covariances show that,the corresponding analyses and forecasts can lead to superior quality,compared with those using covariances that just introduce influences of larger-or smallerscale errors.Due to the interactions among different scale errors included in optimal background error covariances,relevant analysis increments can properly describe weather systems(processes)at different scales,such as dynamic lifting,thermodynamic instability and advection of moisture at large scale,high-level and low-level jet at synoptic scale,and convective systems at mesoscale and small scale,as well as their interactions.As a result,the corresponding forecasts can be improved.  相似文献   

12.
We have evaluated the simulation of Indian summer monsoon and its intraseasonal oscillations in the National Centers for Environmental Prediction climate forecast system model version 2 (CFSv2). The dry bias over the Indian landmass in the mean monsoon rainfall is one of the major concerns. In spite of this dry bias, CFSv2 shows a reasonable northward propagation of convection at intraseasonal (30–60 day) time scale. In order to document and understand this dry bias over the Indian landmass in CFSv2 simulations, a two pronged investigation is carried out on the two major facets of Indian summer monsoon: one, the air–sea interactions and two, the large scale vertical heating structure in the model. Our analysis shows a possible bias in the co-evolution of convection and sea surface temperature in CFSv2 over the equatorial Indian Ocean. It is also found that the simulated large scale vertical heat source (Q1) and moisture sink (Q2) over the Indian region are biased relative to observational estimates. Finally, this study provides a possible explanation for the dry precipitation bias over the Indian landmass in the simulated mean monsoon on the basis of the biases associated with the simulated ocean–atmospheric processes and the vertical heating structure. This study also throws some light on the puzzle of CFSv2 exhibiting a reasonable northward propagation at the intraseasonal time scale (30–60 day) despite a drier monsoon over the Indian land mass.  相似文献   

13.
The performance of ECHAM5 atmospheric general circulation model (AGCM) is evaluated to simulate the seasonal mean and intraseasonal variability of Indian summer monsoon (ISM). The model is simulated at two different vertical resolutions, with 19 and 31 levels (L19 and L31, respectively), using observed monthly mean sea surface temperature and compared with the observation. The analyses examine the biases present in the internal dynamics of the model in simulating the mean monsoon and the evolution of the boreal summer intraseasonal oscillation (BSISO) and attempts to unveil the reason behind them. The model reasonably simulates the seasonal mean-state of the atmosphere during ISM. However, some notable discrepancies are found in the simulated summer mean moisture and rainfall distribution. Both the vertical resolutions, overestimate the seasonal mean precipitation over the oceanic regions, but underestimate the precipitation over the Indian landmass. The performance of the model improves with the increment of the vertical resolution. The AGCM reasonably simulates some salient features of BSISO, but fails to show the eastward propagation of the convection across the Maritime Continent in L19 simulation. The propagation across the Maritime Continent and tilted rainband structure improve as one moves from L19 to L31. The model unlikely shows prominent westward propagation that originates over the tropical western Pacific region. L31 also produces some of the observed characteristics of the northward propagating BSISOs. However, the northward propagating convection becomes stationary in phase 5–7. The simulation of shallow diabatic heating structure and the heavy rainfall activity over the Bay of Bengal indicate the abundance of the premature convection-generated precipitation events in the model. It is found that the moist physics is responsible for the poor simulation of the northward propagating convection anomalies.  相似文献   

14.
We present a detailed analysis of summer monsoon rainfall over the Indian peninsular using nonlinear spatial correlations. This analysis is carried out employing the tools of complex networks and a measure of nonlinear correlation for point processes such as rainfall, called event synchronization. This study provides valuable insights into the spatial organization, scales, and structure of the 90th and 94th percentile rainfall events during the Indian summer monsoon (June–September). We furthermore analyse the influence of different critical synoptic atmospheric systems and the impact of the steep Himalayan topography on rainfall patterns. The presented method not only helps us in visualising the structure of the extreme-event rainfall fields, but also identifies the water vapor pathways and decadal-scale moisture sinks over the region. Furthermore a simple scheme based on complex networks is presented to decipher the spatial intricacies and temporal evolution of monsoonal rainfall patterns over the last 6 decades.  相似文献   

15.
GRAPES-GFS模式暴雨预报天气学检验特征   总被引：5，自引：4，他引：1  
宫宇  代刊  徐珺  杨舒楠  唐健  张芳  胡宁  张夕迪  沈晓琳 《气象》2018,44(9):1148-1159
本文采用天气学检验方法,对2016年度国家气象中心GRAPES全球数值预报系统(GRAPES-GFS)业务预报暴雨过程及2013-2015年部分回算个例进行了检验,并结合对比欧洲中期天气预报中心确定性预报模式(EC模式)和国家气象中心全球谱模式T639L60(T639模式)降水预报,梳理总结业务GRAPES-GFS模式预报性能优势和系统性偏差特征。被检验暴雨过程共38次,其中南方暴雨过程20次,北方暴雨过程6次,热带扰动或台风降水过程12次。依靠预报员主观天气学检验分析,从降水预报效果检验出发,结合主要影响天气系统和示踪物理量检验,梳理总结模式预报系统性偏差,以期全面发掘该业务预报模式性能。结果表明对短期时效内的降水预报,GRAPES-GFS模式预报稳定性较好,整体明显优于T639模式。但还存在诸如对对流性降水预报较实况偏北或对主雨带南侧暖区降水预报不足的偏差特征;另对弱高空波动背景下的对流性降水预报偏弱;而在降水预报强度大致正确的情况下,对降水系统南侧偏南气流控制区域预报湿度偏大,对副热带地区的低涡系统预报偏强。  相似文献   

16.
采用(美国环保部的MODEL-3系统的)CMAQ源同化模型及4种不同空间分辨率的SO2、NO2实测资料,反演得到中国不同尺度的同化修正排放源,利用新一代中尺度气象模式WRF与多尺度空气质量模式CMAQ,模拟分析了中国不同观测信息密度对SO2、NO2源同化反演及其浓度预报的影响,重点分析了华北地区SO2、NO2浓度加密观测对改善SO2、NO2排放源和空气质量预报的重要影响。结果表明,采用不同分辨率的实测资料时,SO2、NO2的趋势预报效果改善程度有一定差异;采用较高分辨率的实测资料进行SO2、NO2源同化修正时,可明显减小SO2、NO2浓度的预报误差。华北地区较高分辨率的观测信息对于改进源同化修正效果及SO2、NO2浓度的趋势预报十分重要,尤其是对SO2浓度的预报尤为重要;采用经高分辨率的实测资料同化修正的排放源时,WRF-CMAQ模式对北京城市尺度SO2、NO2浓度的变化趋势、浓度水平和空间分布特征具有较好的预报效果。高分辨率的观测资料和区域源同化反演方法对于区域污染物浓度预报及排放源清单具有显著的改进作用。  相似文献   

17.
The multi-scale weather systems associated with a mei-yu front and the corresponding heavy precipitation during a particular heavy rainfall event that occurred on 4 5 July 2003 in east China were successfully simulated through rainfall assimilation using the PSU/NCAR non-hydrostatic, mesoscale, numerical model （MM5） and its four-dimensional, variational, data assimilation （4DVAR） system. For this case, the improvement of the process via the 4DVAR rainfall assimilation into the simulation of mesoscale precipitation systems is investigated. With the rainfall assimilation, the convection is triggered at the right location and time, and the evolution and spatial distribution of the mesoscale convective systems （MCSs） are also more correctly simulated. Through the interactions between MCSs and the weather systems at different scales, including the low-level jet and mei-yu front, the simulation of the entire mei-yu weather system is significantly improved, both during the data assimilation window and the subsequent 12-h period. The results suggest that the rainfall assimilation first provides positive impact at the convective scale and the influences are then propagated upscale to the meso- and sub-synoptic scales.
Through a set of sensitive experiments designed to evaluate the impact of different initial variables on the simulation of mei-yu heavy rainfall, it was found that the moisture field and meridional wind had the strongest effect during the convection initialization stage, however, after the convection was fully triggered, all of the variables at the initial condition seemed to have comparable importance.  相似文献   

18.
ＰｅｃｕｌｉａｒＴｅｍｐｏｒａｌＳｔｒｕｃｔｕｒｅｏｆｔｈｅＳｏｕｔｈＣｈｉｎａＳｅａＳｕｍｍｅｒＭｏｎｓｏｏｎＢｉｎＷａｎｇ①ａｎｄＲｅｎｇｕａｎｇＷｕＤｅｐａｒｔｍｅｎｔｏｆＭｅｔｅｏｒｏｌｏｇｙ，ＵｎｉｖｅｒｓｉｔｙｏｆＨａｗａｉ，ＵＳＡＲｅｃ...  相似文献   

19.
Peculiar temporal structure of the south china sea summer monsoon   总被引：8，自引：0，他引：8  
Beijing located at the junction of four major components of the Asian-Australia monsoon system (the Indian, the western North Pacific, the East Asian subtropical, and the Indonesian-Australian monsoons), the monsoon cli-mate over the South China Sea (SCS) exhibits some unique features. Evidences are presented in this paper to reveal and document the following distinctive features in the temporal structure of the SCS summer monsoon:(1) pronounced monsoon singularities in the lower tropospheric monsoon flows which include the pre-onset and withdrawal easterly surges and the southwesterly monsoon bursts at Julian pentad 34-35 (June 15-24) and pentad 46-47 (August 14-23);(2) four prominent subseasonal cycles (alternative occurrences of climatological active and break monsoons);(3) considerably larger year-to-year variations in convective activity on intraseasonal time scale compared to those over the Bay of Bengal and the Philippine Sea;(4) the redness of the climatological mean spectrum of precipitation / deep convection on synoptic to intraseasona] time scales in the central SCS;(5) a remarkable asymmetry in the seasonal transitions between summer and winter monsoons and an extremely abrupt mid-May transition (the outburst of monsoon rain and the sudden switch in tie lower troposphere winds from an easterly to a westerly regime);(6) the bi-modal interannual variation of summer monsoon onset (normal and delayed modes).In addition, the monsoon rainfall displays enormous east-west gradient over the central SCS. Possible causes for these features are discussed. A number of specific science questions concerning some of the peculiar features are raised for the forthcoming SCS monsoon experiment to address  相似文献   

20.
This study examines the feasibility of using a variable resolution global general circulation model (GCM), with telescopic zooming and enhanced resolution (~35 km) over South Asia, to better understand regional aspects of the South Asian monsoon rainfall distribution and the interactions between monsoon circulation and precipitation. For this purpose, two sets of ten member realizations are produced with and without zooming using the LMDZ (Laboratoire Meteorologie Dynamique and Z stands for zoom) GCM. The simulations without zoom correspond to a uniform 1° × 1° grid with the same total number of grid points as in the zoom version. So the grid of the zoomed simulations is finer inside the region of interest but coarser outside. The use of these finer and coarser resolution ensemble members allows us to examine the impact of resolution on the overall quality of the simulated regional monsoon fields. It is found that the monsoon simulation with high-resolution zooming greatly improves the representation of the southwesterly monsoon flow and the heavy precipitation along the narrow orography of the Western Ghats, the northeastern mountain slopes and northern Bay of Bengal (BOB). A realistic Monsoon Trough (MT) is also noticed in the zoomed simulation, together with remarkable improvements in representing the associated precipitation and circulation features, as well as the large-scale organization of meso-scale convective systems over the MT region. Additionally, a more reasonable simulation of the monsoon synoptic disturbances (lows and disturbances) along the MT is noted in the high-resolution zoomed simulation. On the other hand, the no-zoom version has limitations in capturing the depressions and their movement, so that the MT zone is relatively dry in this case. Overall, the results from this work demonstrate the usefulness of the high-resolution variable resolution LMDZ model in realistically capturing the interactions among the monsoon large-scale dynamics, the synoptic systems and the meso-scale convective systems, which are essential elements of the South Asian monsoon system.  相似文献   

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

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