首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
采用1979—2005年美国大气海洋局(NOAA)的卫星观测资料和IPCC第5次全球气候变化比较试验(CMIP5)的模式资料,对全球对流层和平流层近26 a的气温趋势进行了研究。结果表明,CMIP5模拟的全球平均大气温度趋势与观测结果较一致,能够再现平流层冷却和对流层增温等特点,但是在气温趋势的经纬度分布上,模式资料与观测资料间存在较大差异,同时模式间也存在明显的不一致。与观测资料相比,CMIP5模式资料低估了平流层在热带地区的降温速率,而且明显高估了对流层中部到平流层下层的南极区域的降温趋势。不同CMIP5模式间的最大标准方差出现在平流层的南北极区域,但是在对流层所有纬度上标准方差都保持着较小值。  相似文献   

2.
A series of model experiments with the coupled Max-Planck-Institute ECHAM5/OM climate model have been investigated and compared with microwave measurements from the Microwave Sounding Unit (MSU) and re-analysis data for the period 1979?C2008. The evaluation is carried out by computing the Temperature in the Lower Troposphere (TLT) and Temperature in the Middle Troposphere (TMT) using the MSU weights from both University of Alabama (UAH) and Remote Sensing Systems (RSS) and restricting the study to primarily the tropical oceans. When forced by analysed sea surface temperature the model reproduces accurately the time-evolution of the mean outgoing tropospheric microwave radiation especially over tropical oceans but with a minor bias towards higher temperatures in the upper troposphere. The latest reanalyses data from the 25?year Japanese re-analysis (JRA25) and European Center for Medium Range Weather Forecasts Interim Reanalysis are in very close agreement with the time-evolution of the MSU data with a correlation of 0.98 and 0.96, respectively. The re-analysis trends are similar to the trends obtained from UAH but smaller than the trends from RSS. Comparison of TLT, computed from observations from UAH and RSS, with Sea Surface Temperature indicates that RSS has a warm bias after 1993. In order to identify the significance of the tropospheric linear temperature trends we determined the natural variability of 30-year trends from a 500?year control integration of the coupled ECHAM5 model. The model exhibits natural unforced variations of the 30?year tropospheric trend that vary within ±0.2?K/decade for the tropical oceans. This general result is supported by similar results from the Geophysical Fluid Dynamics Laboratory (GFDL) coupled climate model. Present MSU observations from UAH for the period 1979?C2008 are well within this range but RSS is close to the upper positive limit of this variability. We have also compared the trend of the vertical lapse rate over the tropical oceans assuming that the difference between TLT and TMT is an approximate measure of the lapse rate. The TLT?CTMT trend is larger in both the measurements and in the JRA25 than in the model runs by 0.04?C0.06?K/decade. Furthermore, a calculation of all 30?year TLT?CTMT trends of the unforced 500-year integration vary between ±0.03?K/decade suggesting that the models have a minor systematic warm bias in the upper troposphere.  相似文献   

3.
Temperature trends in the upper stratosphere are investigated using satellite measurements from Stratospheric Sounding Unit(SSU) outputs and simulations from chemistry–climate models(CCMs) and the Coupled Model Intercomparison Project Phase 6(CMIP6). Observational evidence shows a lack of cooling in the Antarctic, in contrast to strong cooling at other latitudes, during austral winter over 1979–97. Analysis of CCM simulations for a longer period of1961–97 also shows a significant contrast in the...  相似文献   

4.
Based on three groups of datasets that include radiosondes, reanalyses, and climate model simulations (e.g., Coupled Model Intercomparison Project, CMIP3) from 1979 to 2008, the interannual variability, global temperature trends, and their uncertainty using ensemble spread among intra-group and inter-group datasets have been discussed. The results show that the interannual temperature variability increased from the troposphere to stratosphere, and the maximum occurs around 50?hPa. The CMIP3 climate models have the largest discrepancy in the stratosphere. The intra-group correlations at 500?hPa generally show high similarity within each data group while the inter-group correlations between reanalyses and the CMIP3 climate model simulations indicate lesser similarity. In contrast, the inter-group correlation at 50?hPa is improved except with the Japanese 25-year Reanalysis Project (JRA-25) dataset, and the Twentieth Century Reanalysis (20CR) reanalysis shows a weak cross correlation. The global temperature trends are highly dependent on the individual data sources. Compared to the radiosondes, the reanalyses show a large ensemble spread of trends in the stratosphere, and the CMIP3 climate model simulations have a large ensemble spread in the height of the crossover point where tropospheric warming changes into stratospheric cooling. The largest ensemble spread among the reanalyses in the stratosphere is mainly from the large discrepancy in the JRA-25 reanalysis after 1998 and a relatively weak anomaly in the 20CR before 1986. The largest ensemble spread among the CMIP3 climate models in the troposphere is related to the influence of both volcanic eruptions and El Ni?o/La Ni?a–Southern Oscillation events. The strong anomalies corresponding to the volcanic eruptions of El Chichon in 1982 and Mt Pinatubo in 1991 are clearly identified in the stratosphere. These volcanic eruptions reduced the warming in the troposphere and strengthened the cooling in the stratosphere during the most recent 30?years.  相似文献   

5.
郭艳君  王国复 《气象学报》2019,77(6):1073-1085
基于118站探空资料研究了近60年中国850—100 hPa气温变化趋势及季节和区域特征,并通过与1979—2017年卫星微波气温的对比研究了中国探空气温均一化的不确定性。研究表明,1958—2017年中国平均对流层气温呈上升趋势,300 hPa升温最为显著,平流层下层(100 hPa)为降温趋势。冬季对流层上层升温趋势和夏季平流层下层降温趋势较强。1979—2017年较整个时段对流层升温趋势较强,平流层下层降温趋势较弱。青藏高原和西北地区对流层上层升温趋势较强。通过与卫星微波气温和邻近探空站探空气温的对比以及均一化前后日夜气温差值检测出中国探空均一化气温仍残存非均一性问题。由于参照序列的局限性,均一化未能完全去除21世纪最初10年中国探空系统变化造成的对流层中、上层至平流层下层气温系统性下降的影响,导致中国对流层上层升温趋势被低估和平流层下层降温趋势被高估。未来可通过参考卫星微波气温和邻近探空站序列调整非均一性订正顺序并增加合理性检验等方法改进中国探空气温均一化方案。   相似文献   

6.
In this research, tropopause temperature (TT) and tropopause geopotential height (TGH) over the inner-core and environmental regions of all tropical cyclones (TCs) over the northwest of the Indian Ocean (NWIO) from 1990 to 2019 were investigated. To this aim, observational/analysis/reanalysis data and also simulated data from both historical and Representative Concentration Pathway 8.5 (RCP8.5) experiments of some global climate models (GCMs) from the Coupled Model Intercomparision Project (CMIP5) were used. Dynamical and thermo-dynamical environmental factors controlling TC, together with their correlation with different phases of some climatic indices were considered. Results indicated that the eastern part of the NWIO was more favorable for TC genesis and intensification.Lower-level stratospheric (upper-level tropospheric) cooling (warming) was detected over the NWIO during 1990−2019. Over the both inner-core and environmental regions of the NWIO TCs, the coldest tropopause occurred at a CS-Category and the warmest tropopause happened at the first stage of a VSCS event. Over the inner-core (environmental) region, the highest tropopause was detected at the first stage of a CS event (at the end of a VSCS life cycle). A significant majority of the used CMIP5 GCMs produced stratospheric cooling and tropospheric warming trends over the NWIO, similar to those obtained using ERA5 reanalysis dataset. Finally, the decreasing trend of TT over the both inner-core and environmental regions of NWIO TCs together with temperature decreasing trend obtained from the CMIP5 GCMs simulations suggest that the NWIO is prone to experience more TCs, especially the intense ones, in the future.  相似文献   

7.
Total column ozone (TCO) over the Tibetan Plateau (TP) is lower than that over other regions at the same latitude, particularly in summer. This feature is known as the “TP ozone valley”. This study evaluates long-term changes in TCO and the ozone valley over the TP from 1984 to 2100 using Coupled Model Intercomparison Project Phase 6 (CMIP6). The TP ozone valley consists of two low centers, one is located in the upper troposphere and lower stratosphere (UTLS), and the other is in the middle and upper stratosphere. Overall, the CMIP6 models simulate the low ozone center in the UTLS well and capture the spatial characteristics and seasonal cycle of the TP ozone valley, with spatial correlation coefficients between the modeled TCO and the Multi Sensor Reanalysis version 2 (MSR2) TCO observations greater than 0.8 for all CMIP6 models. Further analysis reveals that models which use fully coupled and online stratospheric chemistry schemes simulate the anticorrelation between the 150 hPa geopotential height and zonal anomaly of TCO over the TP better than models without interactive chemistry schemes. This suggests that coupled chemical-radiative-dynamical processes play a key role in the simulation of the TP ozone valley. Most CMIP6 models underestimate the low center in the middle and upper stratosphere when compared with the Microwave Limb Sounder (MLS) observations. However, the bias in the middle and upper stratospheric ozone simulations has a marginal effect on the simulation of the TP ozone valley. Most CMIP6 models predict the TP ozone valley in summer will deepen in the future.  相似文献   

8.
In most climate simulations used by the Intergovernmental Panel on Climate Change 2007 fourth assessment report, stratospheric processes are only poorly represented. For example, climatological or simple specifications of time-varying ozone concentrations are imposed and the quasi-biennial oscillation (QBO) of equatorial stratospheric zonal wind is absent. Here we investigate the impact of an improved stratospheric representation using two sets of perturbed simulations with the Hadley Centre coupled ocean atmosphere model HadGEM1 with natural and anthropogenic forcings for the 1979–2003 period. In the first set of simulations, the usual zonal mean ozone climatology with superimposed trends is replaced with a time series of observed zonal mean ozone distributions that includes interannual variability associated with the solar cycle, QBO and volcanic eruptions. In addition to this, the second set of perturbed simulations includes a scheme in which the stratospheric zonal wind in the tropics is relaxed to appropriate zonal mean values obtained from the ERA-40 re-analysis, thus forcing a QBO. Both of these changes are applied strictly to the stratosphere only. The improved ozone field results in an improved simulation of the stepwise temperature transitions observed in the lower stratosphere in the aftermath of the two major recent volcanic eruptions. The contribution of the solar cycle signal in the ozone field to this improved representation of the stepwise cooling is discussed. The improved ozone field and also the QBO result in an improved simulation of observed trends, both globally and at tropical latitudes. The Eulerian upwelling in the lower stratosphere in the equatorial region is enhanced by the improved ozone field and is affected by the QBO relaxation, yet neither induces a significant change in the upwelling trend.  相似文献   

9.
Stratospheric water vapor variations, which may play an important role in surface climate, have drawn extensive studies. Here, the variation in stratospheric water vapor is investigated by using data from observations of the Microwave Limb Sounder(MLS) on the Aura satellite, from the ECMWF Interim Reanalysis(ERAI), and simulations by the Whole Atmosphere Community Climate Model(WACCM). We find that the differences of annual mean stratospheric water vapor among these datasets may be partly caused by the differences in vertical transports. Using budget analysis, we find that the upward transport of water vapor at 100 h Pa is mainly located over the Pacific warm pool region and South America in the equatorial tropics in boreal winter and over the southeast of the South Asian high and south of North America in boreal summer. It is found that temperature averaged over regions with upward transport is a better indicator of interannual variability of tropical mean stratospheric water vapor than the tropical mean temperature. It seems that the distributions of the seasonal cycle amplitude of lower stratospheric water vapor in the tropics can also be impacted by the vertical transport. The radiative effects of the interannual changes in water vapor in the lowermost stratosphere are underestimated by approximately 29% in both ERAI and WACCM compared to MLS, although the interannual variations of water vapor in the lowermost stratosphere are dramatically overestimated in ERAI and WACCM. The results here indicate that the radiative effect of long-term changes in water vapor in the lowermost stratosphere may be underestimated in both ERAI and WACCM simulations.  相似文献   

10.
利用美国大气海洋局卫星应用和研究实验室 (The Center for Satellite Applications and Research,STAR) 提供的MSU/AMSU卫星微波亮温资料V3.0版本,结合三套再分析资料数据集,通过对海洋上空不同高度、不同区域及不同季节的适用性分析,来探讨MSU/AMSU资料在热带海洋区域高空大气的温度变化特征,并通过合成分析揭示亮温资料与海洋的响应关系,从而探讨MSU/AMSU资料在热带海洋区域上的适用性和科学性。结果表明:(1)MSU/AMSU亮温资料在30 °E~70 °W,15 °S~15 °N范围的热带海洋区域适用性较好;(2)热带海洋区域对流层上层和中层大气均呈增温趋势,变化速率分别为0.045 K/(10 a) 和0.107 K/(10 a),增温突变现象出现在1980年代末—1990年代初,平流层低层大气呈降温趋势,变化速率为-0.345 K/(10 a),降温突变现象出现在1990年代中期;(3)在热带海洋区域,高空大气温度的变化趋势具有较强的区域性特征,相对于中东太平洋而言,印度洋-西太平洋区域的增、降温趋势变化更显著。对流层的增温幅度随高度的升高而有所降低。平流层低层的降温趋势在季节内变化不大,而对流层则是秋、冬季的增温趋势要明显大于春、夏季,冬季的增温尤为明显;(4)MSU/AMSU亮温资料对热带海洋温度异常有很好的响应关系,能在弥补海洋区域观测资料稀缺的情况下,对海洋区域起着较好的监测作用。   相似文献   

11.
Changes in temperature and precipitation extremes in the CMIP5 ensemble   总被引:6,自引:1,他引:5  
Twenty-year temperature and precipitation extremes and their projected future changes are evaluated in an ensemble of climate models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5), updating a similar study based on the CMIP3 ensemble. The projected changes are documented for three radiative forcing scenarios. The performance of the CMIP5 models in simulating 20-year temperature and precipitation extremes is comparable to that of the CMIP3 ensemble. The models simulate late 20th century warm extremes reasonably well, compared to estimates from reanalyses. The model discrepancies in simulating cold extremes are generally larger than those for warm extremes. Simulated late 20th century precipitation extremes are plausible in the extratropics but uncertainty in extreme precipitation in the tropics and subtropics remains very large, both in the models and the observationally-constrained datasets. Consistent with CMIP3 results, CMIP5 cold extremes generally warm faster than warm extremes, mainly in regions where snow and sea-ice retreat with global warming. There are tropical and subtropical regions where warming rates of warm extremes exceed those of cold extremes. Relative changes in the intensity of precipitation extremes generally exceed relative changes in annual mean precipitation. The corresponding waiting times for late 20th century extreme precipitation events are reduced almost everywhere, except for a few subtropical regions. The CMIP5 planetary sensitivity in extreme precipitation is about 6 %/°C, with generally lower values over extratropical land.  相似文献   

12.
CMIP5模式对中国年平均气温模拟及其与CMIP3模式的比较   总被引:5,自引:0,他引:5  
利用CRUT3v和CN05两套观测资料,评估25个CMIP5模式对1906-2005年中国年平均气温变化的模拟能力,并与CMIP3模式对比。结果表明:1906-2005年中国平均温升速率为0.84℃/100a,CMIP5多模式集合平均模拟的增温率为0.77℃/100a。模式对20世纪后期温升模拟好于前期,仅有两个模式能模拟中国20世纪40年代异常增暖。模式对气温气候态空间分布模拟较好,但在中国西部地区存在最大模拟冷偏差和不确定性。1961-1999年,中国北方增暖大于南方。多模式集合平均可以较好地模拟气温变化线性趋势的空间分布,但对南北气温变化趋势的差异模拟过小。总体说来,在中国平均气温变化趋势、气温气候态空间分布和气温变化趋势空间分布三方面,CMIP5模式都较CMIP3模式有所提高。  相似文献   

13.
Using radiosonde and satellite observations, we investigated the trends of air temperature changes over the Tibetan Plateau (TP) in comparison with those over other regions in the same latitudes from 1979 to 2002. It is shown that Over the TP, the trends of air temperature changes in the upper troposphere to lower stratosphere were out of phase with those in the lower to middle troposphere. Air temperature decreased and a decreasing trend appeared in the upper troposphere to lower stratosphere. The amplitude of the annual or seasonal mean temperature decreases over the TP was larger than that over the whole globe. In the lower to middle troposphere over the TP, temperature increased, and the increasing trend was stronger than that over the non-plateau regions in the same latitudes in the eastern part of China. Meanwhile, an analysis of the satellite observed ozone data in the same period of 1979-2002 shows that over the TP, the total ozone amount declined in all seasons, and the ozone depleted the most compared with the situations in other regions in the same latitudes. It is proposed that the difference between the ozone depletion over the TP and that over other regions in the same latitudes may lead to the difference in air temperature changes. Because of the aggravated depletion of ozone over the TP, less (more) ultraviolet radiation was absorbed in the upper troposphere to lower stratosphere (lower to middle troposphere) over the TP, which favored a stronger cooling in the upper troposphere to lower stratosphere, and an intenser heating in the lower to middle troposphere over the TP. Therefore, the comparatively more depletion of ozone over the TP is possibly a reason for the difference between the air temperature changes over the TP and those over other regions in the same latitudes.  相似文献   

14.
Observed hiatus or accelerated warming phenomena are compared with numerical simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) archives,and the associated physical mechanisms are explored based on the CMIP5 models.Decadal trends in total ocean heat content (OHC) are strongly constrained by net top-of-atmosphere (TOA) radiation.During hiatus decades,most CMIP5 models exhibit a significant decrease in the SST and upper OHC and a significant increase of heat penetrating into the subsurface or deep ocean,opposite to the accelerated warming decades.The shallow meridional overturning of the Pacific subtropical cell experiences a significant strengthening (slowdown) for the hiatus (accelerated warming) decades associated with the strengthened (weakened) trade winds over the tropical Pacific.Both surface heating and ocean dynamics contribute to the decadal changes in SST over the Indian Ocean,and the Indonesian Throughflow has a close relationship with the changes of subsurface temperature in the Indian Ocean.The Atlantic Meridional Overturing Circulation (Antarctic Bottom Water) tends to weaken (strengthen) during hiatus decades,opposite to the accelerated warming decades.In short,the results highlight the important roles of air-sea interactions and ocean circulations for modulation of surface and subsurface temperature.  相似文献   

15.
BCC-ESM1 is the first version of the Beijing Climate Center’s Earth System Model,and is participating in phase 6 of the Coupled Model Intercomparison Project(CMIP6).The Aerosol Chemistry Model Intercomparison Project(AerChemMIP)is the only CMIP6-endorsed MIP in which BCC-ESM1 is involved.All AerChemMIP experiments in priority 1 and seven experiments in priorities 2 and 3 have been conducted.The DECK(Diagnostic,Evaluation and Characterization of Klima)and CMIP historical simulations have also been run as the entry card of CMIP6.The AerChemMIP outputs from BCC-ESM1 have been widely used in recent atmospheric chemistry studies.To facilitate the use of the BCC-ESM1 datasets,this study describes the experiment settings and summarizes the model outputs in detail.Preliminary evaluations of BCC-ESM1 are also presented,revealing that:the climate sensitivities of BCC-ESM1 are well within the likely ranges suggested by IPCC AR5;the spatial structures of annual mean surface air temperature and precipitation can be reasonably captured,despite some common precipitation biases as in CMIP5 and CMIP6 models;a spurious cooling bias from the 1960s to 1990s is evident in BCC-ESM1,as in most other ESMs;and the mean states of surface sulfate concentrations can also be reasonably reproduced,as well as their temporal evolution at regional scales.These datasets have been archived on the Earth System Grid Federation(ESGF)node for atmospheric chemistry studies.  相似文献   

16.
The trace gases (O3, HCl, CH4, H2O, NO, NO2) in the stratosphere play an important role, not only in the photochemical processes in which the ozone layer destroyed, but also in the radiative processes. In this paper, we review the works on the distribution and variation of the trace gases in the stratosphere and their impact on climate, which have been carried out at the University of Science and Technology of China in the recent 20 years. The Halogen Occultation Experiment (HALOE) data were used to analyse the distribution and variation of the mixing ratio of these trace gases and the temperature trends in the stratosphere in the most recent decade. And the reanalyzed National Centers of Environmental Prediction (NCEP)/NCAR data were also used to give the temperature trends and compared with the results from HALOE data. Numerical simulations were also carried out to study the impact of ozone depletion on the global climate. In this review, the distributions of the trace gases, especially those over the Qinghai-Xizang Plateau, are discussed, and the variations and trends for the trace gases in various levels in the stratosphere have been given for the most recent decade. The temperature variation and the cooling trend obtained from HALOE data in the middle and lower stratosphere for the last 13 years are significant, which agree well with the results from NCEP/NCAR data. While the temperature trend in the upper stratosphere in this period do not seem to have much cooling. The numerical simulations show that either the Antarctic ozone hole or the ozone valley over Qinghai-Xizang Plateau affect not only the temperature and circulation in the stratosphere, but also the temperature, pressure and wind fields in the troposphere, then lead to the global climate change.  相似文献   

17.
Summary Using high altitude rocketsonde data for Thumba (8.5 N, 76.9E) and Balasore (21.5 N, 86.9E) and the stratalert messages for high latitudes for the winter (December–March) 1984–1985, an examination has been made to study the perturbations in the temperature and winds in the tropical middle atmosphere and their linkage with the dynamical events occurring over the high-latitude middle atmosphere during that winter.The results of analysis indicated occurrence of strong cooling in the mesosphere over a period of seven days (5–12 December 1984) and the depth of the cooling layer was 15 km. This incident was followed by a strong warming over a period of seven days (12–19 December 1984) and the depth of the warming layer was 13 km. The major warming event, which occurred over high latitudes during the later part of December and the first week of January, was followed by the cooling in the mesosphere and warming in the stratosphere at Thumba. Also the zonal winds were strong easterly and the meridional winds were northerly in the upper stratosphere and the lower mesosphere over tropics during the same period. Weaker zonal winds/stronger easterly winds were generally noticed to be associated with coolings/warmings over tropics.With 7 Figures  相似文献   

18.
19.
探空资料中的人为误差对中国温度长期变化趋势的影响   总被引:9,自引:2,他引:7  
郭艳君  李庆祥  丁一汇 《大气科学》2009,33(6):1309-1318
利用1958~2005年探空温度序列, 通过质量控制、均一化处理和序列缺测率分析, 探讨了探空资料中人为误差对中国高空气温变化趋势的影响。中国探空温度序列存在明显的间断点, 72%的序列包含2~4个间断点。相应的订正总体上降低了1958年以来平流层低层降温和对流层升温趋势, 如700 hPa和100 hPa平均趋势值分别降低0.12 K/10 a和0.04 K/10 a。缺测率是气温区域平均趋势估算的重要参数, 30%作为最大缺测率是中国探空温度序列适宜的取样标准。提高取样标准 (台站数减少) 使1958~2005年间对流层上层和平流层下层的降温趋势减弱。中国高空气温变化趋势与全球或北半球大体一致, 但也有不同特点: 500 hPa以下大气趋于升温, 以上则趋于降温, 最大降温趋势位于对流上部的300 hPa, 而且各气候区间区域差异性十分明显。  相似文献   

20.
本文利用37个CMIP5模式和CESM(Community Earth System Model)包含40个成员的超级集合试验的表面气温预估数据,比较了工业革命前气候参照试验、多项式拟合法和方差分析方法这三种目前在国际上运用较多的方法所估算的表面气温内部变率的异同,分析了内部变率的估算对气候预估中信号萌芽时间(TOE)的影响。结果表明:若采用CMIP5多模式集合,则工业革命前气候参照试验和多项式拟合法都是估算内部变率的合理方法,而方差分析方法则由于包含模式性能自身的影响会夸大内部变率故不推荐使用。内部变率的全球分布呈现出极向强化的现象,中高纬度地区的内部变率幅度远大于热带、副热带地区。内部变率受不同排放情景的影响较小,且随时间无显著变化,但方差分析方法估算的内部变率在热带地区容易受到排放情景的影响。若基于类似CESM这样的单个气候模式的超级集合模拟试验来估算内部变率,三种方法估算的结果相似。不同方法估算的内部变率对TOE的影响主要位于北大西洋拉布拉多海、南大洋威德尔海和罗斯海等邻近海洋深对流区。对于中国区域平均来说,基于CESM超级集合模拟试验,三种方法估算的内部变率与强迫信号之比都小于15%;对CMIP5多模式集合,采用工业革命前气候参照试验和多项式拟合法得到的结果与此接近,但若采用方差分析方法则显著高估内部变率的作用。  相似文献   

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

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