首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
A latent heating peak in the PBL was detected in a simulation by a global GCM that failed to reproduce Madden-Julian Oscillation(MJO).The latent heating peak in the PBL was generated by very shallow convection,which prevented moisture from being transported to the free troposphere.Large amount of moisture was therefore confined to the PBL,leading to a dry bias in the free atmosphere.Suffering from this dry bias,deep convection became lethargic,and MJO signals failed to occur.When the latent heating peak in the PBL was removed in another simulation,reasonable MJO signals,including the eastward propagation and the structure of its large-scale circulation,appeared.We therefore propose that the excessive latent heating peak in the PBL due to hyperactive shallow convection may be a reason for a lack of MJO signals in some simulations by other GCMs as well.  相似文献   

2.
 Experiments using a GCM with two different vertical resolutions show differences in the amount of variability in the tropical upper tropospheric zonal wind component associated with the Madden-Julian Oscillation (MJO). The GCM with lower vertical resolution shows very little variability in this quantity whereas when the vertical resolution is doubled in the free troposphere, the GCM produces variability which is of the same strength as observations. However, the eastward propagation of an enhanced convective region from the Indian Ocean into the west Pacific is not well represented in either simulation of this atmospheric GCM. A water-covered or “aqua-planet” version of the same GCM is used to investigate the behaviour of tropical convection when the vertical resolution is doubled. When the vertical resolution is increased, the spectrum of tropical cloud types changes from a bimodal distribution with peaks representing shallow cumulus and deep cumulonimbus clouds to a trimodal distribution with a third peak in mid-troposphere near the melting level. Associated with periods when these mid-level congestus clouds are dominant, the detrainment from these clouds significantly moistens the mid-troposphere. The appearance of these congestus clouds is shown to be partly due to improved resolution of the freezing level and the convective processes occurring at this level. However, due to the way in which convective detrainment is parametrized in this model, the vertical profile becomes rather noisy and this too contributes to the change in the nature of the convective clouds. The resulting cloud distribution more closely resembles observations, particularly during the suppressed phase of the MJO when cumulus congestus is the dominant cloud type. Received: 17 April 2000 / Accepted: 30 November 2000  相似文献   

3.
Vertical cumulus momentum transport is an important physical process in the tropical atmosphere and plays a key role in the evolution of the tropical atmospheric system. This paper focuses on the impact of the vertical cumulus momentum transport on Madden-Julian Oscillation (MJO) simulation in two global climate models (GCMs). The Tiedtke cumulus parameterization scheme is applied to both GCMs [CAM2 and Spectral Atmospheric general circulation Model of LASG/IAP (SAMIL)]. It is found that the MJO simulation ability might be influenced by the vertical cumulus momentum transport through the cumulus parameterization scheme. However, the use of vertical momentum transport in different models provides different results. In order to improve model's MJO simulation ability, we must introduce vertical cumulus momentum transport in a more reasonable way into models. Furthermore, the coherence of the parameterization and the underlying model also need to be considered.  相似文献   

4.
利用NCEP/NCAR再分析资料和中科院大气物理研究所PIAP3大气环流模式,分析了印度洋偶极子对夏季中国南海西南季风水汽输送的影响。结果表明,印度洋偶极子正位相期间夏季中国南海西南水汽输送较强,负位相期间则较弱。原因可归结为以下:正位相期间,MJO(Madden-Julian Oscillation)多活动于热带西印度洋,其向东传播受到阻碍,但经向传播明显,通常可传播至孟加拉湾地区,同时PIAP3显示印度洋季风槽位置偏北,且印尼以西过赤道气流较强,从而使得这一地区气旋性环流得到建立与加强。孟加拉湾地区对应着较强的对流活动以及深厚积云对流加热,从而通过对流加热的二级热力响应使西太平洋副热带高压位置向北推进,进而使得南海地区西南季风水汽输送得到建立与加强。在此期间孟加拉湾、中南半岛至南海地区对流活动较强,而苏门答腊沿岸对流活动受到抑制,由此增强了Reverse-Hadley环流,使低层经向风较强,进而增强了南海西南季风的水汽输送,PIAP3大气环流模式证实了Reverse-Hadley环流的增强。负位相期间,MJO多活动于热带东印度洋,在东传过程中受到Walker环流配置影响,在140°E赤道附近形成东西向非对称积云对流加热热源,其东侧Kelvin波响应加强了东风异常并配合副热带高压南缘东风压制了中国南海的西南季风水汽输送。在此期间,MJO在南海地区的经向传播较强,但经向传播常止步于南海地区15°N附近,虽携带大量水汽,但深厚积云对流强烈地消耗水汽使大气中水汽含量降低,PIAP3大气环流模式证实负位相期间深厚积云对流对水汽消耗加大,从而使得负位相期间南海地区水汽含量与正位相期间大体相近,但由于经向风不足使水汽向北输送较弱。  相似文献   

5.
Westerly Wind Events (WWEs) over the tropical Pacific are characterized by their spatiotemporal diversity and classified into six types based on their different locations and durations. Various types of WWEs exhibit quite different characteristics in terms of their amplitudes. The long events with a 10–30-day duration are stronger than the short events with a 6–10-day duration, and the maximum amplitudes of the Central-Pacific (C) type and Eastern-Pacific (E) type of long events are larger than the Western-Pacific (W) type of long ones. The evolutions of these six types of WWEs are also quite different. The W-type short and long events and the C-type long events show a distinct eastward propagation, whereas the C-type short events and the E-type short and long events have no apparent propagation direction. We demonstrate that such a difference in the eastward propagation of WWEs can be significantly associated with the Madden-Julian Oscillation (MJO). The W-type events are more influenced heavily by the MJO, as indicated by their more distinct eastward propagation patterns, than the two other types of WWEs. In addition to the MJO, the convectively coupled Kelvin waves are also associated with the WWEs, especially for the short events.  相似文献   

6.
采用一种基于降水异常追踪MJO(Madden–Julian Oscillation)东传的MJO识别方法(MJO Tracking)评估了参与MJOTF/GASS(MJO Task Force/Global Energy and Water Cycle Experiment Atmospheric System Study)全球模式比较计划的全海气耦合模式(CNRM-CM)、半海气耦合模式(CNRM-ACM)和大气模式(CNRM-AM)1991~2010年模拟MJO的能力,探究了海气耦合过程对模式模拟MJO能力的影响机理。CNRM-CM模式模拟的MJO结构更加接近观测,该模式不仅具有最高的MJO生成频率,也能够模拟较强的MJO强度以及较远的传播距离。海气耦合过程会造成CNRM-CM和CNRM-ACM模式中印度洋—太平洋暖池区域海温气候态的冷偏差。但是这种海温气候态的偏差基本没有改变模式模拟MJO的能力。CNRM-CM中MJO对流中心东(西)侧存在较强的季节内尺度海温暖(冷)异常,纬向梯度明显,而CNRM-ACM和CNRM-AM中不存在这样的海温东西不对称结构。结果表明在CNRM模式中海气耦合过程调控模式海温季节内尺度变率对模式MJO模拟能力的影响比调控模式海温气候态更加重要。  相似文献   

7.
8.
Bimodal representation of the tropical intraseasonal oscillation   总被引:2,自引:1,他引:1  
The tropical intraseasonal oscillation (ISO) shows distinct variability centers and propagation patterns between boreal winter and summer. To accurately represent the state of the ISO at any particular time of a year, a bimodal ISO index was developed. It consists of Madden-Julian Oscillation (MJO) mode with predominant eastward propagation along the equator and Boreal Summer ISO (BSISO) mode with prominent northward propagation and large variability in off-equatorial monsoon trough regions. The spatial–temporal patterns of the MJO and BSISO modes are identified with the extended empirical orthogonal function analysis of 31?years (1979–2009) OLR data for the December–February and June–August period, respectively. The dominant mode of the ISO at any given time can be judged by the proportions of the OLR anomalies projected onto the two modes. The bimodal ISO index provides objective and quantitative measures on the annual and interannual variations of the predominant ISO modes. It is shown that from December to April the MJO mode dominates while from June to October the BSISO mode dominates. May and November are transitional months when the predominant mode changes from one to the other. It is also shown that the fractional variance reconstructed based on the bimodal index is significantly higher than the counterpart reconstructed based on the Wheeler and Hendon’s index. The bimodal ISO index provides a reliable real time monitoring skill, too. The method and results provide critical information in assessing models’ performance to reproduce the ISO and developing further research on predictability of the ISO and are also useful for a variety of scientific and practical purposes.  相似文献   

9.
The boreal summer intraseasonal oscillation (BSISO) of the Asian summer monsoon (ASM) is one of the most prominent sources of short-term climate variability in the global monsoon system. Compared with the related Madden-Julian Oscillation (MJO) it is more complex in nature, with prominent northward propagation and variability extending much further from the equator. In order to facilitate detection, monitoring and prediction of the BSISO we suggest two real-time indices: BSISO1 and BSISO2, based on multivariate empirical orthogonal function (MV-EOF) analysis of daily anomalies of outgoing longwave radiation (OLR) and zonal wind at 850 hPa (U850) in the region 10°S–40°N, 40°–160°E, for the extended boreal summer (May–October) season over the 30-year period 1981–2010. BSISO1 is defined by the first two principal components (PCs) of the MV-EOF analysis, which together represent the canonical northward propagating variability that often occurs in conjunction with the eastward MJO with quasi-oscillating periods of 30–60 days. BSISO2 is defined by the third and fourth PCs, which together mainly capture the northward/northwestward propagating variability with periods of 10–30 days during primarily the pre-monsoon and monsoon-onset season. The BSISO1 circulation cells are more Rossby wave like with a northwest to southeast slope, whereas the circulation associated with BSISO2 is more elongated and front-like with a southwest to northeast slope. BSISO2 is shown to modulate the timing of the onset of Indian and South China Sea monsoons. Together, the two BSISO indices are capable of describing a large fraction of the total intraseasonal variability in the ASM region, and better represent the northward and northwestward propagation than the real-time multivariate MJO (RMM) index of Wheeler and Hendon.  相似文献   

10.
The modulation of twin tropical cyclogenesis in the Indian-western Pacific Oceans by the Madden-Julian Oscillation (MJO) during the onset period of 1997/98 ENSO is explored for the period of September 1996 to June 1997 based on daily OLR, NCEP/NCAR wind vector, and JTWC best track datasets. The MJO westerly wind burst associated with its eastward propagation can result in a series of tropical cyclogeneses in a multi-day interval. Only in the transition seasons are pairs of tropical cyclones observed in both the tropical sectors of the Indian-western Pacific Oceans. Two remarkable twin tropical cyclogeneses probably modulated by the MJO westerly wind burst are found: one is observed in the Indian Ocean in the middle of October 1996, and the other is observed in the Western Pacific Ocean in late May 1997. The twin tropical cyclogenesis in mid-October 1996 is observed when the super cloud cluster separates into two isolated clusters by the enhanced westerly wind, which is accompanied by two independent vortices in the equatorial tropical sectors. The other one, in late-May 1997, however, is characterized by one cyclonic flow that later results in another cyclonic cell in its opposite equatorial sector. Thus, there are two very important conditions for twin cyclogenesis: one is the MJO westerly wind straddling the equator, and the other is the integral super cloud cluster, which later splits into two cloud convective clusters with independent vortices.  相似文献   

11.
With a hybrid atmosphere-ocean coupled model we carried out an experimental forecast of a well documented Madden-Julian Oscillation (MJO) event that was observed during the period of Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE). The observed event, originated in the western Indian Ocean around 6 January 1993, moved eastward with a phase speed of about 6.2 m s-1 and reached the dateline around February 1. The hybrid coupled model reasonably forecasts the MJO initiation in the western Indian Ocean, but the predicted MJO event propagates too slow (~ 4.4 m s-1). Results from previous observational studies using unprecedented humidity profiles obtained by NASA Aqua/AIRS satellite suggested that two potential physical processes may be responsible for this model caveat. After improving the cumulus parameterization scheme based on the observations, the model is able to forecast the same event one month ahead. Further sensitivity experiment confirms that the speed-up of model MJO propagation is primarily due to the improved convective scheme. Further, air-sea coupling plays an important role in maintaining the intensity of the predicted MJO. The results here suggest that MJO prediction skill is sensitive to model cumulus parameterization and air-sea coupling. Citation: Fu, X., B. Wang, Q. Bao, et al., 2008: Experimental dynamical forecast of an MJO event observed during TOGA-COARE period, Atmos. Oceanic Sci. Lett., 1, 24-28  相似文献   

12.
The sensitivity of the simulated tropical intraseasonal oscillation or MJO (Madden and Julian oscilla tion)to different cumulus parameterizations is studied by using an atmospheric general circulation model (GCM)-SAMIL(Spectral Atmospheric Model of IAP LASG).Results show that performance of the model in simulating the MJO alters widely when using two different cumulus parameterization schemes-the moist convective adjustment scheme(MCA)and the Zhang-McFarlane(ZM)scheme.MJO simulated by the MCA scheme was found to be more realistic than that simulated by the ZM scheme.MJO produced by the ZM scheme is too weak and shows little propagation characteristics.Weak moisture convergence at low levels simulated by the ZM scheme is not enough to maintain the structure and the eastward propagation of the oscillation.These two cumulus schemes produced different vertical structures of the heating profile.The heating profile produced by the ZM scheme is nearly uniform with height and the heating is too weak compared to that produced by the MCA,which maybe contributes greatly to the failure of simulating a reasonable MJO.Comparing the simulated MJO by these two schemes indicate that the MJO simulated by the GCM is highly sensitive to cumulus parameterizations implanted in.The diabatic heating profile plays an important role in the performance of the GCM.Three sensitivity experiments with different heating profiles are designed in which modified heating profiles peak respectively in the upper troposphere(UH), middle troposphere(MH),and lower troposphere(LH).Both the LH run and the MH run produce eastward propagating signals on the intraseasonal timescale,while it is interesting that the intraseasonal timescale signals produced by the UH run propagate westward.It indicates that a realistic intraseasonal oscillation is more prone to be excited when the maximum heating concentrates in the middle-low levels,especially in the middle levels,while westward propagating disturbances axe more prone to be produced when the maximum heating appears very high.  相似文献   

13.
The sensitivity of the simulated tropical intraseasonal oscillation or MJO (Madden and Julian oscillation) to different cumulus parameterizations is studied by using an atmospheric general circulation model (GCM)--SAMIL (Spectral Atmospheric Model of IAP LASG). Results show that performance of the model in simulating the MJO alters widely when using two different cumulus parameterization schemes-the moist convective adjustment scheme (MCA) and the Zhang-McFarlane (ZM) scheme. MJO simulated by the MCA scheme was found to be more realistic than that simulated by the ZM scheme. MJO produced by the ZM scheme is too weak and shows little propagation characteristics. Weak moisture convergence at low levels simulated by the ZM scheme is not enough to maintain the structure and the eastward propagation of the oscillation. These two cumulus schemes produced different vertical structures of the heating profile. The heating profile produced by the ZM scheme is nearly uniform with height and the heating is too weak compared to that produced by the MCA, which maybe contributes greatly to the failure of simulating a reasonable MJO. Comparing the simulated MJO by these two schemes indicate that the MJO simulated by the GCM is highly sensitive to cumulus parameterizations implanted in. The diabatic heating profile plays an important role in the performance of the GCM. Three sensitivity experiments with different heating profiles are designed in which modified heating profiles peak respectively in the upper troposphere (UH), middle troposphere (MH), and lower troposphere (LH). Both the LH run and the MH run produce eastward propagating signals on the intraseasonal timescale, while it is interesting that the intraseasonal timescale signals produced by the UH run propagate westward. It indicates that a realistic intraseasonal oscillation is more prone to be excited when the maximum heating concentrates in the middle-low levels, especially in the middle levels, while westward propagating disturbances  相似文献   

14.
A coupled earth system model(ESM) has been developed at the Nanjing University of Information Science and Technology(NUIST) by using version 5.3 of the European Centre Hamburg Model(ECHAM), version 3.4 of the Nucleus for European Modelling of the Ocean(NEMO), and version 4.1 of the Los Alamos sea ice model(CICE). The model is referred to as NUIST ESM1(NESM1). Comprehensive and quantitative metrics are used to assess the model's major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model's assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring–fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific(CP)-ENSO and eastern Pacific(EP)-ENSO; however, the equatorial SST variability,biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden–Julian Oscillation(MJO), and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version(T159) demonstrates improved simulation with respect to the climatology, interannual variance, monsoon–ENSO lead–lag correlation, spatial structures of the leading mode of the Asian–Australian monsoon rainfall variability, and the eastward propagation of the MJO.  相似文献   

15.
于堃  沈新勇  张驰  李小凡 《气象科学》2020,40(3):333-340
利用NCL滤波方法将NCEP提供的FNL风场资料分离出天气尺度,准双周振荡(QBWO,Quasi-Biweekly Oscillation)和热带季节内振荡(MJO,Madden-Julian Oscillation)环流场,研究不同时间尺度环流对台风"天鹅"(1515)突变路径的影响。台风路径的特征能够分3个阶段,其中第二阶段台风发生突然转折。第一阶段,天气尺度上台风东侧的反气旋和QBWO环流场中的波列共同引导台风向西偏北方向运动,而MJO环流场中的引导气流作用较小;第二阶段,天气尺度上台风东侧的反气旋和低频环流场中台风附近的气旋共同促进了"天鹅"近90°的突然转向,其中,高、低频分量分别促使台风突然向北、向东转向;第三阶段,天气尺度上的气旋与反气旋、QBWO环流场中的反气旋以及MJO环流场中的脊共同引导"天鹅"向东北方向运动,其中MJO环流场中气旋附近的偏东风促使"天鹅"向西运动,但由于它被天气尺度上强烈的偏西风所抵消,故"天鹅"仍向东运动。  相似文献   

16.
The Madden-Julian oscillation (MJO) simulated by the Canadian Climate Centre general circulation model (CCC GCM) is identified by a principal oscillation pattern (POP) analysis and compared with that observed in the real atmosphere. The results are based upon two integrations of the CCC GCM, one with a parameterization of penetrative cumulus convection (EXP1) and the other with a moist convective adjustment scheme (EXP2). The signal of MJO can be detected in both integrations as the first POP of the 200 hPa velocity potential along the equator. The disturbances show a distinctive wave number one structure with the strongest local amplitude found in the longitudes corresponding to the region of the Asian monsoon. The phase speed of the eastward wave propagation is higher in the eastern Pacific and lower in the monsoon region where the convective activities are strongest. These features are in good agreement with the observations. The energy spectrum of the velocity potential peaks at the frequency corresponding to a period of about 38 days for EXP1, which is somewhat shorter compared to the observed periods of 40–50 days. On the other hand, two spectral peaks can be clearly identified for EXP2, one with a period of 24 days and the other with a much longer period, somewhere near 112 days. Both peaks appear statistically significant at 95% level. Long term data of the observed atmosphere show little indication of such spectral separation. The horizontal patterns identified by the POP analysis resemble to some extent the baroclinic response of tropical flow to a heat source travelling with the speed of MJO. At the upper level, Rossby wave energy propagates westward with winds generally following the height contours, whereas Kelvin wave energy propagates to the east from the heat source with strong cross-contour flow near the equator. At the lower level, the patterns are essentially reversed. The model-generated precipitation and diabatic heating are examined by compositing against the moving MJO. It is found in EXP2 that the composite heating distribution is coherent with the flow pattern only in a certain sector of the equator, depending on whether the fast or slow mode is used to determine the reference point. The composite vertical heating profile of a slower mode tends to have a maximum found at a lower level. The sensitivity of simulated MJO to the cumulus convection scheme in the model is discussed. Received: 19 December 1994 / Accepted: 11 July 1995  相似文献   

17.
The Madden-Julian oscillation (MJO) is a dominant atmospheric low-frequency mode in the tropics. In this review article, recent progress in understanding the MJO dynamics is described. Firstly, the fundamental physical processes responsible for MJO eastward phase propagation are discussed. Next, a recent modeling result to address why MJO prefers a planetary zonal scale is presented. The effect of the seasonal mean state on distinctive propagation characteristics between northern winter and summer is discussed in a theoretical framework. Then, the observed precursor signals and the physical mechanism of MJO initiation in the western equatorial Indian Ocean are further discussed. Finally, scale interactions between MJO and higher- frequency eddies are delineated.  相似文献   

18.
Impacts of a GCM's Resolution on MJO Simulation   总被引:3,自引:0,他引:3  
Long-term integrations are conducted using the Spectral Atmospheric Model (referred to as SAMIL), which was developed in the Laboratory for Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG) in the Institute of Atmospheric Physics (IAP), with different resolutions to inves-tigate sensitivity of the Madden-Julian Oscillation (MJO) simulations to the model's resolution (horizontal and vertical). Three resolutions of the model, R15L9, R42L9 and R42L26, with identical physical processes, all produced the basic observed features of the MJO, including the spatiotemporal space-time spectra and eastward propagation. No fundamental differences among these simulations were found. This indicates that the model resolution is not a determining factor for simulating the MJO. Detailed differences among these modeling results suggest, however, that model resolution can substantially affect the simulated MJO in certain aspects. For instance, at a lower horizontal resolution, high frequency disturbances were weaker and the structures of the simulated MJO were better defined to a certain extent. A higher vertical resolution led to a more realistic spatiotemporal spectrum and spatial distribution of MJO precipitation. Meanwhile, increasing the model's resolution improved simulation of the climatology. However, increasing the resolution should be based on improving the cumulus parameterization scheme.  相似文献   

19.
The daily outgoing longwave radiation (OLR) field in boreal summer shows significant power spectrum peaks on quasi-biweekly (10–20-day) and intraseasonal (20-80-day) timescales over the Indo-western Pacific warm pool, especially over the South China Sea and Bay of Bengal. The quasi-biweekly oscillation (QBWO) originates from off-equatorial western North Pacific, and is characterized by a northwest-southeast oriented wave train pattern, propagating northwestward. The intraseasonal oscillation (ISO), on the other hand, originates from the equatorial Indian Ocean and propagates eastward and northward. Why the equatorial mode possesses a 20–80-day periodicity while the off-equatorial mode favors a 10–20-day periodicity is investigated through idealized numerical experiments with a 2.5-layer atmospheric model. In the off-equatorial region, the model simulates, under a realistic three-dimensional summer mean flow, the most unstable mode that has a wave train pattern with a typical zonal wavelength of 6000 km and a period of 10–20 days, propagating northwestward. This is in contrast to the equatorial region, where a Madden-Julian oscillation (MJO) like mode with a planetary (wavenumber-1) zonal scale and a period ranging from 20 to 80 days is simulated. Sensitivity experiments with different initial conditions indicate that the QBWO is an intrinsic mode of the atmosphere in boreal summer in the off-equatorial Indo-western Pacific region under the summer mean state, while the MJO is the most unstable mode in the equatorial region.  相似文献   

20.
Based on multiple datasets, correlation and composite analyses, and case studies, this paper investigated possible influences of the Indian Ocean dipole (IOD) mode on the eastward propagation of intraseasonal oscillation in the tropical atmosphere. The results showed that (1) the 30-60 day outgoing longwave radiation anomalies in the southeastern Indian Ocean and the 30-60 day 850-hPa zonal wind anomalies over the equatorial central Indian Ocean were significantly correlated with the IOD index; (2) during positive IOD years, the anomalously cold water in the southeastern Indian Ocean and the 850-hPa anomalous easterlies over the equatorial central Indian Ocean might act as barriers to the continuously eastward propagation of the intraseasonal convection, which interrupts the Madden-Julian oscillation (MJO) propagation in the eastern equatorial Indian Ocean and western Pacific; and (3) during negative IOD years, the anomalously warm water in the southeastern Indian Ocean and the low-level westerly anomalies over the equatorial central Indian Ocean favor the eastward movement of MJO.  相似文献   

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

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