首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day−1, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

2.
We measured the fluxes of sensible and latent heat between a low‐land dipterocarp forest in Peninsular Malaysia and the atmosphere. No clear seasonal or interannual changes in latent heat flux were found from 2003 to 2005, while sensible heat flux sometimes fluctuated depending on the fluctuation of incoming radiation between wet and dry seasons. The evapotranspiration rates averaged for the period between 2003 and 2005 were 2·77 and 3·61 mm day?1 using eddy covariance data without and with an energy balance correction, respectively. Average precipitation was 4·74 mm day?1. Midday surface conductance decreased with an increasing atmospheric water vapour pressure deficit and thus restricted the excess water loss on sunny days in the dry season. However, the relationship between the surface conductance and vapour pressure deficit did not significantly decline with an increase in volumetric soil water content even during a period of extremely low rainfall. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

3.
We use Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) data to estimate spatial energy flux and evaporation distributions at the Salar de Atacama, a playa in Northern Chile. Our approach incorporates ASTER surface kinetic temperature, emissivity, and reflectance data, ground-based meteorological measurements, and empirical parameters. Energy flux distributions are estimated using either spatially constant or spatially distributed values of model parameters, with spatially distributed parameters assigned separately to each land cover category in an image classification. We test the sensitivity of energy budget calculations to state variable and parameter values by conducting Monte Carlo simulations for regions with ground energy budget measurements. Results show that assigning spatially distributed model parameters via land cover classifications yields significant improvements to ground and sensible heat flux predictions. Latent heat fluxes cannot, however, be predicted with sufficient accuracy to allow estimation of area-integrated evaporative moisture loss at this low-evaporation playa.  相似文献   

4.
Atmospheric stability conditions over the water surface can affect the evaporative and convective heat fluxes from the water surface. Atmospheric instability occurred 72.5% of the time and resulted in 44.7 and 89.2% increases in the average and maximum estimated evaporation, respectively, when compared to the neutral condition for a small shallow lake (Binaba) in Ghana. The proposed approach is based on the bulk-aerodynamic transfer method and the Monin-Obukhov similarity theory (MOST) using standard meteorological parameters measured over the surrounding land. For water surface temperature, a crucial parameter in heat flux estimation from water surfaces, an applicable method is proposed. This method was used to compute heat fluxes and compare them with observed heat fluxes. The heat flux model was validated using sensible heat fluxes measured with a 3-D sonic anemometer. The results show that an unstable atmospheric condition has a significant effect in enhancing evaporation alongside the sensible heat flux from water surfaces.  相似文献   

5.
鄱阳湖夏季水热通量特征及环境要素影响分析   总被引:2,自引:2,他引:0  
气候变化加速了全球水文循环过程,然而,气候变化如何影响水体蒸发及其水热通量交换仍然不清楚.基于涡度相关系统观测鄱阳湖水体水热通量过程,在小时和日尺度分析了水热通量的变化规律及其主要影响因子.研究表明,潜热通量日变化波动剧烈,大部分为正值,变化范围在-50~580 W/m2之间.而感热通量数值较小,变化范围在-50~50 W/m2之间.8月份潜热通量和感热通量均呈波动下降趋势,均值分别为167.4和15.9 W/m2.8月份日平均潜热通量和感热通量之和大于净辐射,这是由于这一时段储存在水体中的热量释放并补充潜热通量和感热通量.小时尺度上潜热通量日变化在相位上与净辐射无显著相关性,而与风速显著相关.在日尺度变化趋势上,8月份日平均潜热通量仍主要受到风速和水温的影响,感热通量则主要受到风速和饱和水汽压差的影响.  相似文献   

6.
China's Loess Plateau is located at the edge of the Asian summer monsoon in a transition zone of climate and ecology. In the Loess Plateau, climate and environments change along with space, which has an obvious impact on the spatial distribution of surface energy fluxes. Because of scarce land-surface observation sites and short observation time in this area, previous studies have failed to fully understand the land-surface energy balance characteristics over the entire the Loess Plateau and their effect mechanisms. In this paper, we first test the simulation ability of the Community Land Model(CLM) model by comparing its simulated data with observed data. Based on the simulation data for the Loess Plateau over the past thirty years, we then analyze the spatial distribution of surface energy fluxes and compare the pattern differences between the area averages for the driest year and wettest year. Furthermore, we analyze the relationship between the spatial distribution of the components of the surface energy balance with longitude, latitude, altitude, precipitation and temperature. The main results are as follows: the spatial distribution of surface energy fluxes are significantly different, with the surface net radiation and sensible heat flux increasing from south to north and latent heat flux and soil heat flux decreasing from southeast to northwest. The sensible heat flux at the driest point is nearly twice as high as that at the wettest point, whereas the latent heat flux and soil heat flux at the driest point are half as much as that at the wettest point. The impact of variations of annual precipitation on the components of the surface energy balance is also obvious, and the maximum magnitude of the changes to the sensible heat flux and latent heat flux is nearly 30%. To a certain extent, geographical factors(including longitude, latitude, and altitude) and climate factors(including temperature and precipitation) affect the surface energy fluxes. However, the surface net radiation is more closely related to latitude and altitude, sensible heat flux is more closely related to the monsoon rainfall and latitude, and latent heat flux and soil heat flux are more closely related to the monsoon rainfall.  相似文献   

7.
8.
J. A. Leach  R. D. Moore 《水文研究》2010,24(17):2369-2381
Stream temperature and riparian microclimate were characterized for a 1·5 km wildfire‐disturbed reach of Fishtrap Creek, located north of Kamloops, British Columbia. A deterministic net radiation model was developed using hemispherical canopy images coupled with on‐site microclimate measurements. Modelled net radiation agreed reasonably with measured net radiation. Air temperature and humidity measured at two locations above the stream, separated by 900 m, were generally similar, whereas wind speed was poorly correlated between the two sites. Modelled net radiation varied considerably along the reach, and measurements at a single location did not provide a reliable estimate of the modelled reach average. During summer, net radiation dominated the surface heat exchanges, particularly because the sensible and latent heat fluxes were normally of opposite sign and thus tended to cancel each other. All surface heat fluxes shifted to negative values in autumn and were of similar magnitude through winter. In March, net radiation became positive, but heat gains were cancelled by sensible and latent heat fluxes, which remained negative. A modelling exercise using three canopy cover scenarios (current, simulated pre‐wildfire and simulated complete vegetation removal) showed that net radiation under the standing dead trees was double that modelled for the pre‐fire canopy cover. However, post‐disturbance standing dead trees reduce daytime net radiation reaching the stream surface by one‐third compared with complete vegetation removal. The results of this study have highlighted the need to account for reach‐scale spatial variability of energy exchange processes, especially net radiation, when modelling stream energy budgets. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
降水条件下的典型干旱区陆面特征模拟验证   总被引:6,自引:0,他引:6       下载免费PDF全文
利用“我国西北干旱区陆-气相互作用观测试验"在敦煌双墩子戈壁站取得的观测资料及最近的一些研究成果对陆面模式中反照率、相似性函数及地表粗糙度(动量粗糙度,标量粗糙度)3个方面的参数化方案进行了改进,然后对一次典型降水过程的陆面特征及近地层的风、温、湿进行了模拟.结果表明:改进的模式能对降水条件下的干旱区陆面特征进行较好的模拟,其中对辐射、地表温度的模拟相当好,而对能量的模拟虽然还有要改进的地方,但总体结果令人满意;另外改进后的模式对近地层的温度和湿度的模拟也有明显的改善.  相似文献   

10.
In this study, a parameterization method based on Landsat‐7 Enhanced Thematic Mapper (ETM) data and field observations is presented and tested for deriving the regional land surface variables, vegetation variables and land surface heat fluxes over a heterogeneous landscape. As a case study, the method and two Landsat‐7 ETM images are applied to the Jiddah area of Saudi Arabia. The regional distribution maps of surface reflectance, normalized difference vegetation index, modified soil adjusted vegetation index (MSAVI), vegetation coverage, leaf area index, surface temperature, net radiation flux, soil heat flux, sensible heat flux and latent heat flux have been determined over the Jiddah area. The derived results have been validated by using the ‘ground truth’. The results show that the more reasonable regional distributions of land surface variables (surface reflectance, surface temperature), vegetation variables (MSAVI and vegetation coverage), net radiation, soil heat flux and sensible heat flux can be obtained by using the method proposed in this study. Further improvement of the method is also discussed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

11.
The combination of tree canopy cover and a free water surface makes the subcanopy environment of flooded forested wetlands unlike other aquatic or terrestrial systems. Subcanopy vapour fluxes and energy budgets represent key controls on water level and understorey climate but are not well understood. In a permanently flooded forest in south‐eastern Louisiana, USA, an energy balance approach was used to address (a) whether evaporation from floodwater under a forest canopy is solely energy limited and (b) how energy availability was modulated by radiation and changes in floodwater heat storage. A 5‐month continuous measurement period (June–November) was used to sample across seasonal changes in canopy activity and temperature regimes. Over this period, the subcanopy airspace was humid, maintaining saturation vapour pressure for 28% of the total record. High humidity coupled with the thermal inertia of surface water altered both seasonal and diel energy exchanges, including atypical phenomena such as frequent day‐time vapour pressure gradients towards the water surface. Throughout the study period, nearly all available energy was partitioned to evaporation, with minimal sensible heat exchange. Monthly mean evaporation ranged from 0.7 to 1.7 mm/day, peaking in fall when canopy senescence allowed greater radiation transmission; contemporaneous seasonal temperature shifts and a net release of stored heat from the surface water resulted in energy availability exceeding net radiation by 30% in October and November. Relatively stable energy partitioning matches Priestley–Taylor assumptions for a general model of evaporation in this ecosystem.  相似文献   

12.
Information about seasonal crop water consumption is useful to develop the appropriate irrigation scheme. Measurements of energy balance components using the Bowen ratio method were made for a complete growing season at a vineyard in the arid region of northwest China. Vine in the experiment was furrow‐irrigated using a trellis system. The measured evapotranspiration was compared with estimates using the soil water balance method. It is shown that the Bowen ratio method provided accurate estimates of evapotranspiration from the vineyard and this requires that the Bowen ratio system is appropriately installed. The energy balance components showed typical diurnal pattern with peaks that occurred around the midday, except for the ground heat flux which delayed its peak by 2–3 h. The sensible heat flux was greater than the latent heat flux and followed the net radiation closely. The ratio of the latent heat flux to net radiation was low in the early growing season and increased over time. Under the limited irrigation experienced in the vineyard, the latent heat flux was controlled by available soil moisture and the total evapotranspiration in the growing season was 253 mm. The seasonal progression of the crop coefficient is similar to that reported in the literature, with the maximum occurring during the month of September. The crop coefficient can be estimated as a non‐linear function of day of year (DOY) and used to estimate evapotranspiration from vineyards in the region. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
In order to discuss the values and daily variation characteristics of heat storage fluxes in a tropical seasonal rain forest in Xishuangbanna, the sensible and latent heat storage flux within air column, canopy heat storage flux, energy storage by photosynthesis and ground heat storage above the soil heat flux plate, as well as the ratios of these heat storage fluxes to the net radiation in the cool-dry, hot-dry and rainy season were compared and analyzed based on the observation data of carbon fluxes, meteorological factors and biomass within this tropical seasonal rain forest from January 2003 to December 2004. The findings showed that heat storage terms ranged significantly in the daytime and weakly in the nighttime, and the absolute values of sensible and latent heat storage fluxes were obviously greater than other heat storage terms in all seasons. In addition, the absolute values of total heat storage fluxes reached the peak in the hot-dry season, then were higher in the rainy season, and reached the minimum in the cool-dry season. The ratios of heat storage fluxes to net radiation generally decreased with time in the daytime, moreover, the sensible and latent heat storage dominated a considerable fraction of net radiation, while other heat storage contents occupied a smaller fraction of the net radiation and the peak value was not above 3.5%. In the daytime, the ratios of the total heat storage to net radiation were greater and differences in these ratios were distinct among seasons before 12:00, and then they became lower and differences were small among seasons after 12:00. The energy closure was improved when the storage terms were considered in the energy balance, which indicated that heat storage terms should not been neglected. The energy closure of tropical seasonal rain forest was not very well due to effects of many factors. The results would help us to further understand energy transfer and mass exchange between tropical forest and atmosphere. Moreover, they would supply a research basis for studying energy closure at other places.  相似文献   

14.
Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m−2 and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R2 = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m3 ha−1). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

15.
In the present study an attempt has been made to understand the variation of surface energy fluxes such as net radiation, sensible, latent and soil heat during different epochs of thunderstorm activity at Kharagpur. The study also focuses in delineating the difference in the surface energy budget from the days of thunderstorm activity to fair weather days in the pre-monsoon months (April and May) which is locally known as thunderstorm season. For this purpose, experimental data obtained from the Severe Thunderstorms- Observations and Regional Modeling (STORM) programme during pre-monsoon months of 2007, 2009 and 2010 at Kharagpur (22°30′N, 87°20′E), West Bengal, India are used. The present study reveals quick response, in the order of a few days, in the variations of transport of energy fluxes at soil-atmosphere interface to the upper atmosphere vis-à-vis to the occurrence of thunderstorm activity. Rise of surface sensible heat flux to the level of surface latent heat flux a day or two before the occurrence of a thunderstorm has been identified as a precursor signal for the thunderstorm occurrence over Kharagpur. Distinguishable differences are found in the partitioning of the surface energy fluxes to that of net radiation between thunderstorm and non-thunderstorm days. The present study reveals more Bowen’s ratio during thunderstorm days to that of nonthunderstorm days. These results are useful in validating mesoscale model simulations of thunderstorm activity.  相似文献   

16.
沙漠陆面过程参数化与模拟   总被引:4,自引:2,他引:2       下载免费PDF全文
郑辉  刘树华 《地球物理学报》2013,56(7):2207-2217
沙漠地区植被稀疏、干旱少雨,其陆面物理过程具有与全球其它地区显著不同的特点.本文利用巴丹吉林沙漠观测资料,分析和计算了地表反照率、比辐射率、粗糙度和土壤热容量、热传导系数等关键陆面过程参数,建立了适合于沙漠地区的陆面过程模式DLSM (Desert Land Surface Model),并与NOAH陆面过程模式的模拟结果和观测资料进行了比较.结果表明:巴丹吉林沙漠地表反照率为0.273,比辐射率为0.950,地表粗糙度为1.55×10-3 m,土壤热容量和热扩散系数分别为1.08×106 J·m-3·K-1和3.34×10-7 m2·s.辐射传输、感热输送和土壤热传导过程是影响沙漠地区地表能量平衡的主要物理过程.通过对这三种过程的准确模拟检验,DLSM能够较准确地模拟巴丹吉林沙漠地气能量交换特征;短波辐射、长波辐射和感热通量的模拟结果与观测值间的标准差分别为7.98,6.14,33.9 W·m-2,与NOAH陆面过程模式的7.98,7.72,46.6 W·m-2的结果接近.地表反照率是沙漠地区最重要的陆面过程参数,地表反照率增大5%,向上短波辐射通量随之增加5%,感热通量则减小2.8%.本文研究结果对丰富陆面过程参数化方案,改进全球陆面过程模式、气候模式具有参考意义.  相似文献   

17.
Playa systems are driven by evaporation processes, yet the mechanisms by which evaporation occurs through playa salt crusts are still poorly understood. In this study we examine playa evaporation as it relates to land surface energy fluxes, salt crust characteristics, groundwater and climate at the Salar de Atacama, a 3000 km2 playa in northern Chile containing a uniquely broad range of salt crust types. Land surface energy budget measurements were taken at eight representative sites on this playa during winter (August 2001) and summer (January 2002) seasons. Measured values of net all-wave radiation were highest at vegetated and rough halite crust sites and lowest over smooth, highly reflective salt crusts. Over most of the Salar de Atacama, net radiation was dissipated by means of soil and sensible heat fluxes. Dry salt crusts tended to heat and cool very quickly, whereas soil heating and cooling occurred more gradually at wetter vegetated sites. Sensible heating was strongly linked to wind patterns, with highest sensible heat fluxes occurring on summer days with strong afternoon winds. Very little energy available at the land surface was used to evaporate water. Eddy covariance measurements could only constrain evaporation rates to within 0.1 mm d−1, and some measured evaporation rates were less than this margin of uncertainty. Evaporation rates ranged from 0.1 to 1.1 mm d−1 in smooth salt crusts around the margin of the salar and from 0.4 to 2.8 mm d−1 in vegetated areas. No evaporation was detected from the rugged halite salt crust that covers the interior of the salar, though the depth to groundwater is less than 1 m in this area. These crusts therefore represent a previously unrecorded end member condition in which the salt crusts form a practically impermeable barrier to evaporation.  相似文献   

18.
Estimation of evapotranspiration (ET) is of great significance in modeling the water and energy interactions between land and atmosphere. Negative correlation of surface temperature (Ts) versus vegetation index (VI) from remote sensing data provides diagnosis on the spatial pattern of surface soil moisture and ET. This study further examined the applicability of Ts–VI triangle method with a newly developed edges determination technique in estimating regional evaporative fraction (EF) and ET at MODIS pixel scale through comparison with large aperture scintillometer (LAS) and high‐level eddy covariance measurements collected at Changwu agro‐ecological experiment station from late June to late October, 2009. An algorithm with merely land and atmosphere products from MODIS onboard Terra satellite was used to estimate the surface net radiation (Rn) and soil heat flux. In most cases, the estimated instantaneous Rn was in good agreement with surface measurement with slight overestimation by 12 W/m2. Validation results from LAS measurement showed that the root mean square error is 0.097 for instantaneous EF, 48 W/m2 for instantaneous sensible heat flux, and 30 W/m2 for daily latent heat flux. This paper successfully presents a miniature of the overall capability of Ts–VI triangle in estimating regional EF and ET from limited number of data. For a thorough interpretation, further comprehensive investigation needs to be done with more integration of remote sensing data and in‐situ surface measurements. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
A case study on a desert‐oasis wetland ecosystem in the arid region of Northwest China measured the seasonal and interannual variation in energy partitioning and evapotranspiration to analyse the response of water and energy exchange on soil moisture, groundwater, and environmental variables. Energy partitioning showed a clear seasonal and interannual variability, and the process of water and energy exchange differed significantly in the monthly and interannual scales. The net radiation was 7.31 MJ m?2· day?1, and sensible heat flux accounted for 50.42% of net radiation in energy fluxes, 40.56% for latent heat flux, and 9.02% for ground heat flux. The parameters in energy fluxes were best described by a unimodal curve, whereas sensible heat flux followed a bimodal curve. Variations in the daily evapotranspiration and crop evapotranspiration also exhibited a single peak curve with annual values of 569.84 and 644.47 mm, respectively. Canopy conductance averaged 20.77 ± 13.75 mm s?1 and varied from 0.16 to 83.96 mm s?1 during the two hydrological years. The variation in water and energy exchange reflected environmental conditions and depended primarily on vapour pressure deficit, net radiation, soil moisture, and water depth. Although the effects of precipitation on evapotranspiration showed that the response of this ecosystem to climate changes was not obvious, the variation of air temperatures had a strong influence on evapotranspiration, resulting in a significant increase in evapotranspiration (R = 0.730; P < 0.01). Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

20.
The default fractional vegetation cover and terrain height were replaced by the estimated fractional vegetation cover, which was calculated by the Normalized Difference Vegetation Index(NDVI) of Earth Observing System Moderate-Resolution Imaging Spectroradiometer(EOS-MODIS) and the Digital Elevation Model of the Shuttle Radar Topography Mission(SRTM) system. The near-surface meteorological elements over northeastern China were assimilated into the three-dimensional variational data assimilation system(3DVar) module in the Weather Research and Forecasting(WRF) model. The structure and daily variations of air temperature, humidity, wind and energy fields over northeastern China were simulated using the WRF model. Four groups of numerical experiments were performed, and the simulation results were analyzed of latent heat flux, sensible heat flux, and their relationships with changes in the surface energy flux due to soil moisture and precipitation over different surfaces. The simulations were compared with observations of the stations Tongyu, Naiman, Jinzhou, and Miyun from June to August, 2009. The results showed that the WRF model achieves high-quality simulations of the diurnal characteristics of the surface layer temperature, wind direction, net radiation, sensible heat flux, and latent heat flux over semiarid northeastern China in the summer. The simulated near-surface temperature, relative humidity, and wind speed were improved in the data assimilation case(Case 2) compared with control case(Case 1). The simulated sensible heat fluxes and surface heat fluxes were improved by the land surface parameterization case(Case 3) and the combined case(Case 4). The simulated temporal variations in soil moisture over the northeastern arid areas agree well with observations in Case 4, but the simulated precipitation should be improved in the WRF model. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations. The assimilation datasets generated by this work can be applied to research on climate change and environmental monitoring of arid lands, as well as research on the formation and stability of climate over semiarid areas.  相似文献   

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

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