首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
—In order to understand the physical mechanism for maintenance of low frequency waves, computations of nonlinear kinetic energy (KE) exchanges into individual triad interactions in frequency domain by use of cross-spectral technique over global tropics (20°S–30°N) at 850 hPa and 200 hPa for June, July and August of 1988 are carried out. The KE exchanges among different frequencies and the interactions between the transient eddies and the time mean flow are found to be an order of magnitude smaller in the lower troposphere than those in the upper troposphere. The results show two distinct spectral peaks of periods 45 and 25 days in nonlinear energy transfer in the lower troposphere. The former is more pronounced than the latter. The role of time mean flow on the low frequency transients is found to be secondary compared to the effect of the leading term due to nonlinear interactions in the lower troposphere. Low frequency waves suffer a net loss of energy in the upper troposphere. In the lower troposphere, north of 20°N low frequency waves lose energy through nonlinear triad interactions, unlike the upper troposphere where gain of energy is noticed. Longitude-frequency distributions suggest that wave-CISK process and strong gradient of SST are the possible mechanisms for the strong energy interactions associated with low frequency waves in the lower troposphere over the west Pacific and east coast of Africa, respectively. The study may aid investigation of the rapid loss of predictability of low frequency modes over the tropics.  相似文献   

2.
— The work deals with the computation and analysis of spectral energetics in the frequency domain at 850?hPa and 200?hPa over the tropics (20°S–20°N) and extratropics (20°N–60°N). The data for the winter months, i.e., November, December and January of 1995, 1996 and 1997 are selected for this purpose. The results suggest that much of the low frequency variability of the Northern Hemisphere wintertime general circulation is associated with disturbances which derive their energy from the time-mean flow through barotropic instability. Low frequency fluctuations tend to be larger in horizontal scale and their kinetic energy is largely confined to the upper troposphere. At 850?hPa, strong energy interaction south of 5°N is noticed due to a southward shift of major inflow channel, originating from the Bay of Bengal and entering the ITCZ from the western Arabian Sea. The energy balances in the tropics and the extratropics during winter have different characteristics from those during summer. In contrast to the summer circulation, instead of a downscale decascade as in the case of the extratropics, kinetic energy is transferred in an opposite sense, namely from transients of shorter to those of longer time scales in the tropics during winter. The strong nonlinear energy interactions associated with low frequency waves over the Indian Ocean (5°N–5°S) during winter is the manifestation of the deep convection due to warm water coupled with the crossequatorial low level flow along the ITCZ over this region. Forcing from this region readily excites a large response in terms of nonlinear energy interaction over the extratropical northeast Pacific.  相似文献   

3.
A mean climatology is studied to examine atmospheric circulation characteristics to assess the wintertime (December, January, February and March - DJFM) synoptic weather system affecting northern India. The main objective is to study the mean circulation and mean energetics distribution pertaining to the winter season, which are embedded with an eastward moving synoptic weather system in westerlies, called Western Disturbances (WDs). Forty years (1958–1997) of uninitialized daily re-analysis data of the National Center for Environmental Prediction - National Center for Atmospheric Research (NCEP- NCAR, henceafter NCEP), U.S. has been considered for this study. Winter circulations are considered over the domain 15°S–45°N and 30°E–120°E. This domain is considered particularly to illustrate the impact of wintertime synoptic weather system Western Disturbances (WDs), which travel towards the east over the western Himalayas during winter and yield an enormous amount of precipitation in the form of snow. Large-scale balances of kinetic energy, vorticity, angular momentum, heat and moisture budget terms are analyzed. The main findings of the study show that strong rising motion in the extratropical region brings a significant amount of precipitation over the region of study. Also, horizontal flux of kinetic energy converges in the tropical region and diverges over the extratropical region. It is seen that both the zonal and meridional component of kinetic energy contributes to the production of kinetic energy in the upper troposphere. Vorticity budget shows that wintertime circulation over the western Himalayas is characterized by a negative generation of vorticity. The relative and planetary vorticity advection contributes to the horizontal transport of vorticity. The moisture flux transported into the region shows that in the middle tropospheric levels moisture undergoes phase transformation due to turbulent exchange and hence releases latent heat.  相似文献   

4.
The major sudden stratospheric warming (SSW) events of 2003–04 and 2005–06 are considered to investigate changes in equatorial convection due to circulation changes associated with the SSW events. It is observed that the SSW events are accompanied by a considerable decrease in Outgoing Longwave Radiation (OLR), a proxy for tropical convection, over equatorial latitudes (15°N–15°S) in the Indonesian sector (90°E–150°E). However, unlike noted by earlier observations, the zonal mean OLR does not show any notable relationship with the SSW events. It can be explained from the latitude–longitude map of potential vorticity (PV) at 100 hPa, which shows a tongue of high PV emanating from high latitudes towards equator and converges in the longitude band of 90°E–150°E on the day of peak warming at 1 hPa in the case of 2003–04 and 10 hPa in the case of 2005–06. The latitude-height map of Eliassen–Palm (EP) vector and its divergence show convergence of EP flux in the upper troposphere at latitudes even lower than 20°N on these days. Further, vertical winds computed from the convergence of momentum flux are upward indicating convective activity at low-latitudes and downward at mid-latitudes.  相似文献   

5.
This study is diagnostic in character and entails a systematic investigation of the kinematic parameters of the mean monthly 300 mb horizontal-motion field for the global tropical region during January, April, and October. The kinematic parameters analysed include the horizontal velocity divergence, relative vorticity, velocity potential, and stream function. The magnitude of relative voritcity has been observed to be nearly equal to that of the horizontal velocity divergence in the equatorial latitudes and an order of magnitude larger in the subtropics. The magnitudes of?u/?x, ?v/?y, and horizontal velocity divergence are comparable and vary very little with latitude, while that of?u/?y is an order of magnitude larger than that of?v/?x, confirming its dominance in the computation of relative vorticity. The latitudinal position of zero absolute vorticity isopleth fluctuates within 5 degrees of the equator and is confined more to the winter hemisphere than to the summer hemisphere. The rotational component of wind is large compared with the divergent wind component. The magnitude of the rotational vector wind is four times that of the divergent vector wind in the equatorial latitudes and an order of magnitude larger in the subtropical latitudes.  相似文献   

6.
非对称型强飓风中的准平衡流特征分析   总被引:5,自引:0,他引:5       下载免费PDF全文
本文在论述飓风发生发展生命史过程中平衡、准平衡和非平衡态动力学特征的基础上,应用PV-ω方法,对具有非轴对称和长时间强度维持特征的飓风Bonnie(1998)进行了反演诊断分析,结果表明:基于非线性平衡模式的平衡流能够描述飓风水平涡旋场的基本特征,而加入准平衡ω方程得到的准平衡流能反映边界层入流、高层出流、眼墙区的剧烈倾斜上升运动和眼心区域的下沉运动.准平衡流描述了具有较长生命史组织化过程的强对流系统,而与其相联系的辐散运动与涡度同量级,证明了飓风准平衡流场具有涡散运动共存的特征,但在边界层顶的入流急流区和高层出流区仍存在高度非平衡态的超梯度流.利用反演的准平衡流场分析发现,当由环境风场低层到高层存在顺切变时,飓风内中尺度对流带移动方向的左侧,有利于强对流单体的发展和新对流单体的形成,右侧则相反,同时强气旋式旋转流场的作用,使得对流单体形成后随基本气流传播至对流减弱区,造成了飓风非对称结构的形成和维持.  相似文献   

7.
The leading empirical orthogonal function (EOF) of the June-Sept. mean, rotational horizontal wind at 850 hPa and 200 hPa (over the region 12.5°S–42.5°N, 50°E–100°E) from 56 years (1948–2003) of reanalysis (from the National Centers for Environmental Prediction) shows strong anti-cyclonic circulation at upper levels, strong Indian Ocean cross-equatorial flow and on-shore flow over western India at lower levels . The associated principal component (PC) is correlated at the 0.75 level with the seasonal mean observed Indian Monsoon rainfall (IMR). Composite differences of vertically integrated divergence (surface to 800 hPa) and vorticity (surface to 500 hPa) between ``strong' years (PC-1 exceeds one standard deviation σ) and ``weak' years (PC-1 less than − σ) suggest increased rising motion and storminess over the Bay of Bengal and central India. Composite difference maps of station rainfall from the India Meteorological Department (IMD) between strong years and normal years (weak years and normal years) are statistically significant over central India, with strong (weak) years associated with increased (decreased) precipitation. In both cases the maps of rainfall anomalies are of one sign throughout India. The correlation of PC-1 with global seasonal mean SST is strong and negative over the eastern equatorial Pacific, but positive in a surrounding horse-shoe like region. Significant negative correlation occurs in the northwestern Indian Ocean. The lag/lead correlation between the NINO3 SST index and PC-1 is similar to but stronger than the NINO3/IMR correlation. Modest (but significant) negative correlation is seen when NINO3 leads PC-1 (or IMR) by one-two months. Strong negative correlation is seen when PC-1 (or IMR) leads NINO3. The projections of running five-day means of horizontal rotational winds at 850 and 200 hPa onto EOF-1 (after removing the seasonal mean for each year) were pooled for strong, normal and weak years. The strong and normal year probability distribution functions (pdfs) are nearly indistinguishable, but the weak year pdf has more weight for moderate negative values and in both extreme tails and shows some hint of bi-modality.  相似文献   

8.
A diagnostie method of cumulus parameterization is suggested in which vertical transport of horizontal momentum by cumulus-scale is derived by making use of large-scale vorticity as well as divergence budget equations. Data for composite monsoon depression over India available from our earlier studies used to test the method. As a first approximation, the results are obtained using only the vorticity budget equation.The results show that in the southwest sector of the monsoon depression, which is characterized by maximum cloudiness and precipitation, there is an excess of cyclonic vorticity in the lower troposphere and anticyclonic vorticity in the upper troposphere associated with the large-scale motion. The distribution of eddy vertical transport of horizontal momentum is such that anticyclonic vorticity is generated in the lower troposphere and cyclonic vorticity aloft. Cumulus-scale eddies thus work against the large-scale system and tend to off-set the large-scale imbalance in vorticity.  相似文献   

9.
The surface ocean circulation in the Caribbean Sea is characterized by the interaction between anticyclonic eddies and the Caribbean Upwelling System (CUS). These interactions lead to instabilities that modulate the transfer of kinetic energy up- or down-cascade. The interaction of North Brazil Current rings with the islands leads to the formation of submesoscale vorticity filaments leeward of the Lesser Antilles, thus transferring kinetic energy from large to small scales. Within the Caribbean, the upper ocean dynamic ranges from large-scale currents to coastal upwelling filaments and allow the vertical exchange of physical properties and supply KE to larger scales. In this study, we use a regional model with different spatial resolutions (6, 3, and 1 km), focusing on the Guajira Peninsula and the Lesser Antilles in the Caribbean Sea, in order to evaluate the impact of submesoscale processes on the regional KE energy cascade. Ageostrophic velocities emerge as the Rossby number becomes O(1). As model resolution is increased submesoscale motions are more energetic, as seen by the flatter KE spectra when compared to the lower resolution run. KE injection at the large scales is greater in the Guajira region than in the others regions, being more effectively transferred to smaller scales, thus showing that submesoscale dynamics is key in modulating eddy kinetic energy and the energy cascade within the Caribbean Sea.  相似文献   

10.
As early as in the 1980s, Chinese scientists hadfirst proposed that there exits two summer monsoonsystems in Asia, namely the East Asian summer mon-soon (EASM) and the Indian summer monsoon(ISM)[1-4]. The two monsoon systems are quite dif-ferent in characteristics. Since then, such issue andconclusion had been documented and approved by alot of studies in the past two decades, and was appliedin the guideline of the South China Sea summer mon-soon experiment (SCSMEX), which was undertak…  相似文献   

11.
贾建颖  刘毅 《地球物理学报》1954,63(10):3626-3639
利用1979年1月至2016年12月ERA-interim月平均再分析资料和CAMP全球月降水资料,分析夏季(6—8月)南亚高压下方500 hPa到100 hPa暖中心的时空分布,从三维结构来揭示夏季南亚高压暖心特征.回归分析进一步探讨青藏高原上空暖中心对全球大气环流产生的可能影响.结果表明:南亚高压在150 hPa达到最强,这一层也是异常冷暖中心分界面,150 hPa以下有一强大异常暖中心,异常暖中心位于300 hPa附近,150 hPa以上为异常冷中心,中心位置位于70 hPa附近.异常暖中心从500 hPa向上逐渐向西向北倾斜,异常暖中心面积200 hPa达最大,150 hPa异常暖中心消失,100 hPa以上转变为异常冷中心.500~200 hPa异常暖中心表现出不断增暖的长期趋势(1979—2016),100 hPa异常冷中心则表现出不断变冷的长期趋势(1979—2016).去掉长期趋势的时间序列表现出明显的"准两年振荡"特征,异常暖中心位置在纬向上较稳定,在经向上表现出年际的"东西振荡".300 hPa异常暖中心是整个南亚高压的关键层.300 hPa异常暖中心对全球其他变量场进行回归分析.高度回归场表明,青藏高原上空异常暖中心在北半球中高纬度高度场上激发出3波的行星波,波特征在对流层中上层表现明显,波振幅随高度增高不断加强,在对流层中下层逐渐减弱并消失.纬向风回归场在对流层中上层表现出横跨南北半球的波列,这个波列在200 hPa振幅最大.经向风回归场在北半球中纬度(30°N—60°N)表现出7波型,说明南北能量交换频繁.降水回归场表明,东亚地区长江中下游至日本降水偏少,而其南北两侧降水偏多.  相似文献   

12.
The NCEP/NCAR R1 reanalysis data are employed to investigate the impact of forced and inertial instability in the lower troposphere over the Arabian Sea on the onset process of Indian summer monsoon(ISM),and to reveal the important role of zonal advection of zonal geostrophic momentum played in the forced unstable convection.Results show that during the ISM onset the zero absolute vorticity contour(??=0)shifts northward due to the strong cross-equatorial pressure gradient in the lower troposphere over southern Arabian Sea.Thus a region with negative absolute vorticity is generated near the equator in the northern hemisphere,manifesting the evident free inertial instability.When a southerly passes through this region,under the influence of friction a lower convergence that facilitates the convection flourishing at the lower latitudes appears to the north of zero absolute vorticity contour.However,owing to such a traditional inertial instability,the convection is confined near the equator which does not have direct influence on the ISM onset.On the contrary in the region to the north of the zero absolute vorticity contour and to the south of the low pressure center near the surface,although the atmosphere there is inertially stable,the lower westerly jet can develop and bring on the apparent zonal advection of zonal geostrophic momentum.Both theoretical study and diagnosing analysis present that such a zonal advection of geostrophic momentum is closely associated with the zonal asymmetric distribution of meridional land-sea thermal contrast,which induces a convergence center near and further north of the westerly jet in the lower troposphere over the southwestern coast of the Indian Peninsula,providing a favorable lower circulation for the ISM onset.It illustrates that the development of convection over the Arabian Sea in late spring and early summer is not only due to the frictional inertial instability but also strongly affected by the zonal asymmetric distribution of land-sea thermal contrast.Moreover,before the ISM onset due to the eastward development of the South Asian High(SAH)in the upper troposphere,high potential vorticity is transported to the region over the Arabian Sea.Then a local trumpet-shaped stream field is generated to cause the evident upper divergence-pumping effect which favors the ISM onset.When the upper divergence is vertically coupled with the lower convergence resulted from the aforementioned forced unstable convection development near the southwestern coast of Indian Peninsula,the atmospheric baroclinic unstable development is stimulated and the ISM onset is triggered.  相似文献   

13.
Recent observations of the tropical and subtropical atmosphere are interpreted in terms of scaling arguments and wave propagation theory advanced byCharney (1963, 1969).Charney’s idealizations describe the tropical atmosphere in terms of large regions of quasi-nondivergent flow containing small subdomains of heavy convection and divergence, and place emphasis upon the quasi-rotational regions. FGGE (First GARP Global Experiment) observations suggest that strongly divergent local tropical circulations are forced by latent heating and produce important direct modifications of the total wind field. We describe the extent to which the resulting field consists of divergent and rotational components in different analyses of the FGGE data, and present independent supporting documentation of the results in terms of heating estimates and rainfall observations. Local tropical heating rates on the order of 10°C/day are apparently due to latent heat release associated with precipitation rates as large as 6 cm/day during extended periods. The large contribution of the divergent wind is generally underestimated in models that do not retain such energetic local forcings, and this deficiency may be related to general underestimation of tropical-extratropical connections of many linear models. Such connections are commonly cited in relation to El Niño events, the Southern Hemisphere stationary-wave pattern, and in FGGE studies, but are not well simulated in most linear theories. It is not yet clear whether this is an inherent limitation of linear models, or whether the linear models have not yet explored all the potentially relevant ambient states. We explore the latter possibility by construction of a basic state that allows reasonable latitudinal evolution of the wave field. This basic state has zero absolute vorticity gradient throughout the tropics, and deviations linearized about this state are dynamically analogous to a “local” Hadley cell. To the extent that it is appropriate to regard the results in terms of wave propagation, our analysis suggests a prominent role for gravity-inertia waves in the tropics and for the extratropical connections. The relevance of gravity modes to observations and the theoretical explanation of the flat vorticity field remain to be established.  相似文献   

14.
Summary Two approaches are used to describe and distinguish between those regions of the atmosphere which most closely resemble heat engines and those which operate like refrigerators. Estimates of the conversion from eddy potential energy to mean zonal potential energy are presented as a function of height, latitude and season; the region of forced motion (or refrigerator) in the lower stratosphere appears to be closely related to the maximum westerly winds of the upper troposphere. A preliminary attempt is made to relate some values of mass transport derived from the observations of trace substances to the energy conversions.The research reported in this article was sponsored by the Atomic Energy Commission under Contract AT (30-1) 2241.  相似文献   

15.
利用美国航空航天局MERRA(Modern-Era Retrospective Analysis for Research and Applications)再分析资料和MODIS(Moderate-Resolution Imaging Spectroradiometer)卫星资料以及欧洲气象中心ECMWF-Interim(European Centre for Medium-Range Weather Forecasts)再分析资料,分析了发生于青藏高原北侧上空的一次地形重力波事件,并使用中尺度预报模式WRF-ARW.V3.0(Weather Research and Forecasting model,V3.0)对其进行了数值模拟.在此基础上,诊断分析了此次地形重力波在UTLS(Upper Troposphere and Lower Stratosphere)区域造成的物质和能量垂直传输特征.分析结果表明这一中尺度地形重力波信号的水平波长约为600km,与地形扰动水平尺度接近,重力波在对流层中传播的垂直波长约为3km,在垂直方向上随着高度的增加呈现出由东向西倾斜的结构特征.此次地形重力波上传进入平流层并在150hPa附近破碎,波破碎后动量通量在短时间内发生了强烈的衰减,重力波携带的能量在破碎高度附近释放.重力波破碎的同时垂直方向湍流混合变得异常强烈,湍流交换系数可在短时间内增加到背景值的8倍以上,剧烈湍流混合过程导致了对流层上层的空气进入平流层,使下平流层空气出现了位势涡度和臭氧的低值区,在浮力频率的垂直剖面中也可以看到由于地形重力波过程造成的平流层下层浮力频率异常低值区.  相似文献   

16.
The kinetic energy (KE) seasonality has been revealed by satellite altimeters in many oceanic regions. Question about the mechanisms that trigger this seasonality is still challenging. We address this question through the comparison of two numerical simulations. The first one, with a 1/10° horizontal grid spacing, 54 vertical levels, represents dynamics of physical scales larger than 50 km. The second one, with a 1/30° grid spacing, 100 vertical levels, takes into account the dynamics of physical scales down to 16 km. Comparison clearly emphasizes in the whole North Pacific Ocean, not only a significant KE increase by a factor up to three, but also the emergence of seasonal variability when the scale range 16–50 km (called submesoscales in this study) is taken into account. But the mechanisms explaining these KE changes display strong regional contrasts. In high KE regions, such the Kuroshio Extension and the western and eastern subtropics, frontal mixed-layer instabilities appear to be the main mechanism for the emergence of submesoscales in winter. Subsequent inverse kinetic energy cascade leads to the KE seasonality of larger scales. In other regions, in particular in subarctic regions, results suggest that the KE seasonality is principally produced by larger-scale instabilities with typical scales of 100 km and not so much by smaller-scale mixed-layer instabilities. Using arguments from geostrophic turbulence, the submesoscale impact in these regions is assumed to strengthen mesoscale eddies that become more coherent and not quickly dissipated, leading to a KE increase.  相似文献   

17.
The large-scale terms in the vorticity equation are evaluated usingKrishnamurti's (1971a, b) summer mean winds at 200 mb for a global belt from 25°S to 45°N. The production of vorticity by the divergent wind field is found to be imbalanced over all of the tropical and subtropical belt. As a result there is a requirement for a sub-grid scale (space or time) mechanism which removes negative vorticity from the regions of strong divergence (Tibetan and Mexican highlands) and removes positive vorticity from the regions of strong convergence (mid-oceanic troughs) at 200 mb during northern summer at a rate of approximately 4×10–10 sec–2. As suggested byHolton andColton (1972), in regions of strong and persistent convection, such as the Tibetan Plateau, deep cumulus clouds can account for this transport. However, the mechanism for removing positive vorticity in the vicinity of the upper tropospheric mid-oceanic troughs is still an intriguing and open question.On leave-of-absence at the National Science Foundation, Climate Dynamics Research Section.  相似文献   

18.
We present validation studies of MLS V2.2 and V3.3 water vapor(WV) and ozone profiles over the Tibetan Plateau(Naqu and Lhasa) and its adjacent region(Tengchong) respectively by using the balloon-borne Cryogenic Frost point Hygrometer and Electrochemical Concentration Cell ozonesonde. Coincident in situ measurements were selected to compare the MLS V2.2 and V3.3 WV and ozone profiles for understanding the applicability of the two version MLS products over the region. MLS V2.2 and V3.3 WV profiles respectively show their differences within ?2.2±15.7%(n=74) and 0.3±14.9%(n=75) in the stratosphere at and above 82.5 h Pa. Accordingly, at 100 h Pa, the altitude approaching the tropopuase height, differences are within 9.8± 46.0%(n=18) and 23.0±45.8%(n=17), and they are within 21.5±90.6%(n=104) and 6.0±83.4%(n=99) in upper troposphere. The differences of MLS ozone are within ?11.7±16.3%(n=135, V2.2) and 15.6±24.2%(n=305, V3.3) at and above 82.5 h Pa. At 100 h Pa, they are within ?3.5±54.4%(n=27) and ?8.7±41.6%(n=38), and within 18.0±79.1%(n=47) and 34.2±76.6%(n=160) in the upper troposphere. The relative difference of MLS WV and ozone profile has significant oscillation and scatter at upper troposphere and lower stratosphere partly due to the stronger gradients of WV and ozone concentrations here as well the linear interpolation of sonde data for the intercomparison. At and below 70 h Pa, the relative differences of MLS ozone are significantly larger over Lhasa during the Tibetan Plateau "ozone valley" season, which is also the Asian Summer Monsoon period. The MLS ozone differences over the three sites are similar in their vertical distributions during that period. A simple linear correlation analysis between MLS and sonde profiles indicates that the sensitivity of MLS profile products is related to concentrations at each pressure level. The MLS V3.3 product sensitivity is slightly improved for WV at and above 82.5 h Pa, whereas it is not obvious for ozone. The possible factors contributing to the differences of the MLS profile products of WV and ozone are discussed.  相似文献   

19.
20.
Variation of atmospheric thermodynamical structure parameters between days of thunderstorm occurrence and non-occurrence is presented based on data sets obtained during Severe Thunderstorm-Observations and Regional Modeling (STORM) experiments conducted over Kharagpur (22.3°N, 87.2°E) in pre-monsoon season of 2009 and 2010. Potential instability (stable to neutral) is noticed in the lower layers and enhanced (suppressed) convection in the middle troposphere during thunderstorm (non-thunderstorm) days. Low-level jets are observed during all days of the experimental period but with higher intensity on thunderstorm days. Convective available potential energy (CAPE) builds up until thunderstorm occurrence and becomes dissipated soon after, whereas convective inhibition (CIN) is greatly decreased prior to the event on thunderstorm days. In contrast, higher CAPE and CIN are noticed on non-thunderstorm days. Analysis of thermodynamic indices showed that indices including moisture [humidity index (HI) and dew point temperature at 850 hPa (DPT850)] are useful in differentiating thunderstorm from non-thunderstorm days. The present study reveals that significant moisture availability in the lower troposphere in the presence of convective instability conditions results in thunderstorm occurrence at Kharagpur.  相似文献   

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

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