首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
Numerical analysis of flux footprints for different landscapes   总被引：5，自引：0，他引：5  
Summary A model for the canopy – planetary boundary layer flow and scalar transport based on E- closure was applied to estimate footprint for CO2 fluxes over different inhomogeneous landscapes. Hypothetical heterogeneous vegetation patterns – forest with clear-cuts as well as hypothetical heterogeneous relief – a bell-shaped valley and a ridge covered by forest were considered. The distortions of airflow caused by these heterogeneities are shown – the upwind deceleration of the flow at the ridge foot and above valley, acceleration at the crest and the flow separation with the reversed flow pattern at lee slopes of ridge and valley. The disturbances induce changes in scalar flux fields within the atmospheric surface layer comparing to fluxes for homogeneous conditions: at a fixed height the fluxes vary as a function of distance to disturbance. Correspondingly, the flux footprint estimated from model data depends on the location of the point of interest (flux measurement point) and may significantly deviate from that for a flat terrain. It is shown that proposed method could be used for the choice of optimal sensor position for flux measurements over complex terrain as well as for the interpretation of data for existing measurement sites. To illustrate the latter the method was applied for experimental site in Solling, Germany, taking into account the complex topography and vegetation heterogeneities. Results show that in certain situations (summer, neutral stratification, south or north wind) and for a certain sensor location the assumptions of idealized air flow structure could be used for measurement interpretation at this site, though in general, extreme caution should be applied when analytical footprint models are used in the interpretation of flux measurements over complex sites.  相似文献   

2.
北京城市通量足迹及源区分布特征分析   总被引：2，自引：1，他引：1  
刘郁珏  胡非  程雪玲  宋宗朋 《大气科学》2014,38(6):1044-1054
城市下垫面具有不同于其他下垫面的特殊复杂性,通量的源汇分布十分不均,导致直接测量以及定量的描述城市中水汽、二氧化碳等通量变得非常困难。涡动观测系统在城市通量观测中得到广泛应用,但由于涡动协方差观测系统传感器都架设在一定的高度上,其测量结果只反映特定点下垫面或某部分下垫面的物理过程,不能说明观测结果是否具有下垫面属性的代表性,无法反映观测通量的空间变异性。足迹函数(Footprint)的产生就是为了解决这一问题,其又称源权重函数,是描述近地面层表面源或汇的空间分布和仪器观测通量值之间关系的函数。本文采用北京325 m气象塔近1年7层涡动协方差观测数据、超声数据,分析了不同风向、不同稳定度、不同高度下足迹函数所表达的通量贡献区域范围的变化规律。结果表明,在大气稳定条件下时,通量贡献区范围的大小与主风向无显著差异,而不稳定条件下计算结果与主风向无关。在不同高度下所有稳定条件下通量贡献区范围要大于不稳定条件,中性条件则介于两者之间。8 m、16 m高度上的观测结果不能完全代表城市下垫面通量贡献区,47 m以上能够代表城市下垫面通量贡献区,280 m则已经包含了郊区和城区的权重平均。同一大气稳定度条件下,高度越高通量贡献区范围越大,90%通量贡献区范围与观测高度成线性关系,这种线性关系可以预测没有观测高度或者更高处的通量贡献区范围。  相似文献   

3.
Flux Footprint Simulation Downwind of a Forest Edge   总被引：2，自引：2，他引：0  
Surface fluxes, originating from forest patches, are commonly calculated from atmospheric flux measurements at some height above that patch using a correction for flux arising from upwind surfaces. Footprint models have been developed to calculate such a correction. These models commonly assume homogeneous turbulence, resulting in a simulated atmospheric flux equal to the average surface flux in the footprint area. However, atmospheric scalar fluxes downwind of a forest edge have been observed to exceed surface fluxes in the footprint area. Variations in atmospheric turbulence downwind of the forest edge, as simulated with an E – model, can explain enhanced atmospheric scalar fluxes. This E – model is used to calculate the footprint of atmospheric measurements downwind of a forest edge. Atmospheric fluxes appear mainly enhanced as a result of a stronger sensitivity to fluxes from the upwind surface. A sensitivity analysis shows that the fetch over forest, necessary to reach equilibrium between atmospheric fluxes and surface fluxes, tends to be longer for scalar fluxes as compared to momentum fluxes. With increasing forest density, atmospheric fluxes deviate even more strongly from surface fluxes, but over shorter fetches. It is concluded that scalar fluxes over forests are commonly affected by inhomogeneous turbulence over large fetches downwind of an edge. It is recommended to take horizontal variations in turbulence into account when the footprint is calculated for atmospheric flux measurements downwind of a forest edge. The spatially integrated footprint is recommended to describe the ratio between the atmospheric flux and the average surface flux in the footprint.  相似文献   

4.
Horizontal heterogeneity can significantly affect the flux data quality at monitoring sites in complex terrain. In heterogeneous conditions, the adoption of the eddy-covariance technique is contraindicated by the lack of horizontal homogeneity and presence of advective conditions. In addition, uncertainty concerning the sources or sinks influencing a measurement compromises the data interpretation. The consideration of the spatial context of a measurement, defined by a footprint analysis, can therefore provide an important tool for data quality assessment. This study presents an update of an existing footprint-based quality evaluation concept for flux measurement sites in complex terrain. The most significant modifications in the present version are the use of a forward Lagrangian stochastic trajectory model for the determination of the spatial context of the measurements, and the determination of effective roughness lengths with a flux aggregation model in a pre-processing step. Detailed terrain data gathered by remote sensing methods are included. This approach determines spatial structures in the quality of flux data for varying meteorological conditions. The results help to identify terrain influences affecting the quality of flux data, such as dominating obstacles upwind of the site, or slopes biasing the wind field, so that the most suitable footprint regions for the collection of high-quality datasets can be identified. Additionally, the approach can be used to evaluate the performance of a coordinate rotation procedure, and to check to what extent the measured fluxes are representative for a target land-use type.  相似文献   

5.
Determination of fetch requirements for accurateBowen-ratio measurements of latent- and sensible-heatfluxes is more involved than for eddy-correlationmeasurements because Bowen-ratio sensors are locatedat two heights, rather than just one. A simplesolution to the diffusion equation is used to derivean expression for Bowen-ratio fetch requirements,downwind of a step change in surface fluxes. Theserequirements are then compared to eddy-correlationfetch requirements based on the same diffusionequation solution. When the eddy-correlation andupper Bowen-ratio sensor heights are equal, and theavailable energy upwind and downwind of the stepchange is constant, the Bowen-ratio method requiresless fetch than does eddy correlation. Differences infetch requirements between the two methods aregreatest over relatively smooth surfaces. Bowen-ratiofetch can be reduced significantly by lowering thelower sensor, as well as the upper sensor. The Bowen-ratio fetchmodel was tested using data from a fieldexperiment where multiple Bowen-ratio systems weredeployed simultaneously at various fetches and heightsabove a field of bermudagrass. Initial comparisonswere poor, but improved greatly when the model wasmodified (and operated numerically) to account for thelarge roughness of the upwind cotton field.  相似文献   

6.
Eddy-covariance observations above the densely built-up Centre of Nanjing were made from December 2011 to August 2012. Separate eddy-covariance systems installed at two levels on a 36-m tower located on a rooftop were operated simultaneously, and observations grouped into two sectors (A, B) according to the prevalent wind directions. For sector A, where the nearby buildings are all below the lower measurement level, the sensible heat and momentum fluxes are generally greater at the upper level. For sector B, where several high-rise buildings are located upwind, the sensible heat and momentum fluxes at the upper level are close to those at the lower level. The analysis shows that the turbulent eddy characteristics differ between the two wind sectors, leading to a different behaviour of turbulent exchange between the two levels. A hypothesis is proposed that addresses the vertical variation of turbulent fluxes in the urban roughness sublayer (RSL). For sector A, the buildings block the flow, change the trajectory of scalars, and distort the footprint of scalar fluxes; this ‘blocking effect’ is believed to lead to a smaller sensible heat flux above the canopy layer. Such an effect should decrease with height in the RSL, explaining the increase of the observed turbulent heat flux with height. In addition, the presence of non-uniform building heights adversely affects turbulence organization around the canopy top, and likely elevates the inflection point of the mean flow to a higher elevation close to the upper measurement level, where larger shear results in a larger momentum flux. For sector B, wake effects from the nearby high-rise buildings strongly reduce turbulence organization at higher elevations, leading to similar sensible heat and momentum fluxes at both measurement levels.  相似文献   

7.
We describe pragmatic and reliable methods to examine the influence of patch-scale heterogeneities on the uncertainty in long-term eddy-covariance (EC) carbon flux data and to scale between the carbon flux estimates derived from land surface optical remote sensing and directly derived from EC flux measurements on the basis of the assessment of footprint climatology. Three different aged Douglas-fir stands with EC flux towers located on Vancouver Island and part of the Fluxnet Canada Research Network were selected. Monthly, annual and interannual footprint climatologies, unweighted or weighted by carbon fluxes, were produced by a simple model based on an analytical solution of the Eulerian advection-diffusion equation. The dimensions and orientation of the flux footprint depended on the height of the measurement, surface roughness length, wind speed and direction, and atmospheric stability. The weighted footprint climatology varied with the different carbon flux components and was asymmetrically distributed around the tower, and its size and spatial structure significantly varied monthly, seasonally and inter-annually. Gross primary productivity (GPP) maps at 10-m resolution were produced using a tower-mounted multi-angular spectroradiometer, combined with the canopy structural information derived from airborne laser scanning (Lidar) data. The horizontal arrays of footprint climatology were superimposed on the 10-m-resolution GPP maps. Monthly and annual uncertainties in EC flux caused by variations in footprint climatology of the 59-year-old Douglas-fir stand were estimated to be approximately 15–20% based on a comparison of GPP estimates derived from EC and remote sensing measurements, and on sensor location bias analysis. The footprint-variation-induced uncertainty in long-term EC flux measurements was mainly dependent on the site spatial heterogeneity. The bias in carbon flux estimates using spatially-explicit ecological models or tower-based remote sensing at finer scales can be estimated by comparing the footprint-weighted and EC-derived flux estimates. This bias is useful for model parameter optimizing. The optimization of parameters in remote-sensing algorithms or ecosystem models using satellite data will, in turn, increase the accuracy in the upscaled regional carbon flux estimation.  相似文献   

8.
Measurements of fluxes and profiles of wind andtemperature are performed in the roughness layer ofa moderately homogeneous forest location. Weinvestigate to what extent vertical scalar fluxescan be derived from profile measurements. Theinfluence of inhomogeneities in the upwind terrainis investigated with footprint analysis and with aninhomogeneous surface-layer model. Four methods toestimate displacement height are suggested, amongthem is a method involving the structure parameterof the vertical wind. All methods give a decrease ofdisplacement height with increasing wind speed,while roughness length is found to increase withincreasing wind speed. For near-neutral conditionsdimensionless temperature gradients are found to besubstantially lower than the surface-layer valuesfound in the literature for homogeneous terrain with lowvegetation. Dimensionless shear however iscomparable with the surface-layer value. The heightof the roughness layer is 20 times the roughnesslength. Two schemes with locally derived surfaceparameters are tested to derive friction velocityand sensible heat flux from the profilemeasurements. These site specific schemes performsatisfactorily. A third scheme based on surface parameters chosen a priorifrom the literatureperforms significantly worse especially for low windspeed and unstable cases.  相似文献   

9.
A knowledge of the distribution of the contribution of upwind sources to measurements of vertical scalar flux densities is important for the correct interpretation of eddy covariance data. Several approaches have been developed to estimate this so-called footprint function. Here a new approach based on the ensemble-averaged Navier—Stokes equations is presented. Comparisons of numerical results using this approach with results from other studies under a range of environmental conditions show that the model predictions are robust. Moreover, the approach outlined here has the advantage of a potential wide applicability, due to an ability to take into account the heterogeneous nature of underlying surfaces. For example, the model showed that any variations in surface drag, such as must occur in real life heterogeneous canopies, can exert a marked influence of the shape and extent of flux footprints. Indeed, it seems likely that under such circumstances, estimates of surface fluxes will be weighted towards areas of highest foliage density (and therefore quite likely higher photosynthetic rates) close to the measurement sensor. Three-dimensional footprints during the day and night were also determined for a mixed coniferous forest in european Russia. A marked asymmetry of the footprint in the crosswind direction was observed, this being especially pronounced for non-uniform plant distributions involving vegetation types with different morphological and physiological properties. The model also found that, other things being equal, the footprint peak for forest soil respiration is typically over twice the distance from the above canopy measurement sensor compared to that for canopy photosynthesis. This result has important consequences for the interpretation of annual ecosystem carbon balances by the eddy covariance method.  相似文献   

10.
Turf-grass lawns are ubiquitous in the United States. However direct measurements of land–atmosphere fluxes using the eddy-covariance method above lawn ecosystems are challenging due to the typically small dimensions of lawns and the heterogeneity of land use in an urbanised landscape. Given their typically small patch sizes, there is the potential that CO2 fluxes measured above turf-grass lawns may be influenced by nearby CO2 sources such as passing traffic. In this study, we report on two years of eddy-covariance flux measurements above a 1.5 ha turf-grass lawn in which we assess the contribution of nearby traffic emissions to the measured CO2 flux. We use winter data when the vegetation was dormant to develop an empirical estimate of the traffic effect on the measured CO2 fluxes, based on a parametrised version of a three-dimensional Lagrangian footprint model and continuous traffic count data. The CO2 budget of the ecosystem was adjusted by 135gCm−2 in 2007 and by 134gCm−2 in 2008 to determine the natural flux, even though the road crossed the footprint only at its far edge. We show that bottom-up flux estimates based on CO2 emission factors of the passing vehicles, combined with the crosswind-integrated footprint at the distance of the road, agreed very well with the empirical estimate of the traffic contribution that we derived from the eddy-covariance measurements. The approach we developed may be useful for other sites where investigators plan to make eddy-covariance measurements on small patches within heterogeneous landscapes where there are significant contrasts in flux rates. However, we caution that the modelling approach is empirical and will need to be adapted individually to each site.  相似文献   

11.
The objective of this study is to produce two-dimensional maps of the sensible and the latent heat fluxes from airborne measurements, based on the analysis of a flight pattern, called grid flights. A footprint model with along-wind and cross-wind components was used to project the measured fluxes onto the surface map. The method was applied to measurements over Arctic tundra during the Mackenzie Area GEWEX (Global Energy and Water Cycle Experiment) Study (MAGS) 1999. The resulting flux estimates were computed by integration of a wavelet transform, and corrected for long wavelength losses using information from 100 km long regional runs that were conducted close to the grid flights. The random flux error was estimated based on the flight length that is represented in each map element, and a map resolution of 3 × 3 km was chosen in order to keep the average random error of the latent heat flux below 25%. The random error of the sensible heat flux was smaller by a factor of 1.4 on average. An analysis of airborne flux measurements at different altitudes showed no significant increase of flux estimates for measurement heights below 90 m. Thus, the fluxes measured at heights between 48 and 64 m were not corrected for vertical flux divergence. The resulting flux maps provide quantitative two-dimensional estimates of the energy exchange between the surface and the atmosphere during the snow melt period in an Arctic environment, which are well-suited for calibration and validation of numerical models.  相似文献   

12.
A stochastic trajectory model was used to estimate scalar fluxfootprints in neutral stabilityfor canopies of varying leaf area distributions andleaf area indices. An analytical second-order closure model wasused to predict mean wind speed, second moments and the dissipationrate of turbulent kinetic energy within a forest canopy.The influence of source vertical profile on the flux footprint wasexamined. The fetch is longer for surface sourcesthan for sources at higher levels in the canopy. In order tomeasure all the flux components, and thus the total flux, with adesired accuracy, sources were located at the forest floor in thefootprint function estimation. The footprint functions werecalculated for five observation levels above the canopy top. Itwas found that at low observation heights both canopy density andcanopy structure affect the fetch. The higher abovethe canopy top the flux is measured, the more pronounced is the effectof the canopy structure. The forest fetch for flux measurements isstrongly dependent on the required accuracy: The 90% flux fetchis greater by a factor of two or more compared to the 75% fetch. Theupwind distance contributing 75% of flux is as large as 45 timesthe difference between canopy height and the observation heightabove the canopy top, being even larger for low observationlevels.  相似文献   

13.
Evaluation of Two Energy Balance Closure Parametrizations   总被引：1，自引：0，他引：1  
A general lack of energy balance closure indicates that tower-based eddy-covariance (EC) measurements underestimate turbulent heat fluxes, which calls for robust correction schemes. Two parametrization approaches that can be found in the literature were tested using data from the Canadian Twin Otter research aircraft and from tower-based measurements of the German Terrestrial Environmental Observatories (TERENO) programme. Our analysis shows that the approach of Huang et al. (Boundary-Layer Meteorol 127:273–292, 2008), based on large-eddy simulation, is not applicable to typical near-surface flux measurements because it was developed for heights above the surface layer and over homogeneous terrain. The biggest shortcoming of this parametrization is that the grid resolution of the model was too coarse so that the surface layer, where EC measurements are usually made, is not properly resolved. The empirical approach of Panin and Bernhofer (Izvestiya Atmos Oceanic Phys 44:701–716, 2008) considers landscape-level roughness heterogeneities that induce secondary circulations and at least gives a qualitative estimate of the energy balance closure. However, it does not consider any feature of landscape-scale heterogeneity other than surface roughness, such as surface temperature, surface moisture or topography. The failures of both approaches might indicate that the influence of mesoscale structures is not a sufficient explanation for the energy balance closure problem. However, our analysis of different wind-direction sectors shows that the upwind landscape-scale heterogeneity indeed influences the energy balance closure determined from tower flux data. We also analyzed the aircraft measurements with respect to the partitioning of the “missing energy” between sensible and latent heat fluxes and we could confirm the assumption of scalar similarity only for Bowen ratios $\approx $ 1.  相似文献   

14.
On Frequency Response Corrections for Eddy Covariance Flux Measurements   总被引：6，自引：5，他引：1  
Several potentially significant errors are noted in published frequency-response corrections for eddy-covariance flux measurements.  相似文献   

15.
A Simple Parameterisation for Flux Footprint Predictions   总被引：9，自引：6，他引：9  
Flux footprint functions estimate the location and relative importance of passive scalar sources influencing flux measurements at a given receptor height. These footprint estimates strongly vary in size, depending on receptor height, atmospheric stability, and surface roughness. Reliable footprint calculations from, e.g., Lagrangian stochastic models or large-eddy simulations are computationally expensive and cannot readily be computed for long-term observational programs. To facilitate more accessible footprint estimates, a scaling procedure is introduced for flux footprint functions over a range of stratifications from convective to stable, and receptor heights ranging from near the surface to the middle of the boundary layer. It is shown that, when applying this scaling procedure, footprint estimates collapse to an ensemble of similar curves. A simple parameterisation for the scaled footprint estimates is presented. This parameterisation accounts for the influence of the roughness length on the footprint and allows for a quick but precise algebraic footprint estimation.  相似文献   

16.
Attenuation of Scalar Fluxes Measured with Spatially-displaced Sensors   总被引：1，自引：0，他引：1  
Observations from the Horizontal Array Turbulence Study (HATS) field program are used to examine the attenuation of measured scalar fluxes caused by spatial separation between the vertical velocity and scalar sensors. The HATS data show that flux attenuation for streamwise, crosswind, and vertical sensor displacements are each a function of a dimensionless, stability-dependent parameter n m multiplied by the ratio of sensor displacement to measurement height. The scalar flux decays more rapidly with crosswind displacements than for streamwise displacements and decays more rapidly for stable stratification than for unstable stratification. The cospectral flux attenuation model of Kristensen et al. agrees well with the HATS data for streamwise sensor displacements, although it is necessary to include a neglected quadrature spectrum term to explain the observation that flux attenuation is often less with the scalar sensor downwind of the anemometer than for the opposite configuration. A simpler exponential decay model provides good estimates for crosswind sensor displacements, as well as for streamwise sensor displacements with stable stratification. A model similar to that of Lee and Black correctly predicts flux attenuation for a combination of streamwise and crosswind displacements, i.e. as a function of wind direction relative to the sensor displacement. The HATS data for vertical sensor displacements extend the near-neutral results of Kristensen et al. to diabatic stratification and confirm their finding that flux attenuation is less with the scalar sensor located below the anemometer than if the scalar sensor is displaced an equal distance either horizontally or above the anemometer.  相似文献   

17.
The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.  相似文献   

18.
Micrometeorological measurements (including eddy-covariance measurements of the surface fluxes of sensible and latent heat) were performed during the LITFASS-2003 experiment at 13 field sites over different types of land use (forest, lake, grassland, various agricultural crops) in a 20 × 20 km2 area around the Meteorological Observatory Lindenberg (MOL) of the German Meteorological Service (Deutscher Wetterdienst, DWD). Significant differences in the energy fluxes could be found between the major land surface types (forest, farmland, water), but also between the different agricultural crops (cereals, rape, maize). Flux ratios between the different surfaces changed during the course of the experiment as a result of increased water temperature of the lake, changing soil moisture, and of the vegetation development at the farmland sites. The measurements over grass performed at the boundary-layer field site Falkenberg of the MOL were shown to be quite representative for the farmland part of the area. Measurements from the 13 sites were composed into a time series of the area-averaged surface flux by taking into account the data quality of the single flux values from the different sites and the relative occurrence of each surface type in the area. Such composite fluxes could be determined for about 80% of the whole measurement time during the LITFASS-2003 experiment. Comparison of these aggregated surface fluxes with area-averaged fluxes from long-range scintillometer measurements and from airborne measurements showed good agreement.  相似文献   

19.
The Tibetan Plateau plays an important role in the global water cycle and is strongly influenced by climate change. While energy and matter fluxes have been more intensely studied over land surfaces, a large proportion of lakes have either been neglected or parameterised with simple bulk approaches. Therefore, turbulent fluxes were measured over wet grassland and a shallow lake with a single eddy-covariance complex at the shoreline in the Nam Co basin in summer 2009. Footprint analysis was used to split observations according to the underlying surface, and two sophisticated surface models were utilised to derive gap-free time series. Results were then compared with observations and simulations from a nearby eddy-covariance station over dry grassland, yielding pronounced differences. Observations and footprint integrated simulations compared well, even for situations with flux contributions including grassland and lake. The accessibility problem for EC measurements on lakes can be overcome by combining standard meteorological measurements at the shoreline with model simulations, only requiring representative estimates of lake surface temperature.  相似文献   

20.
Radiometric measurements of hemispherical surface reflectance and long-wave irradiance are required to quantify the broadband albedo and the outgoing thermal radiation. These observations are typically integrated with eddy covariance measurements of sensible and latent heat fluxes to characterize the surface energy budget. While the aerodynamic footprint has been widely investigated, the geometry of the hemispherical radiometric footprint over plant canopies has been rarely tackled. In the present work, the size and shape of the hemispherical radiometric footprint are formalized for a bare surface and in presence of a vegetation cover. For this purpose, four idealized canopies are analyzed and the dependency of the radiometric footprint on leaf area index and canopy height is explored. Besides, the radiometric footprint is compared with the aerodynamic footprint in conditions of neutral stability. It was observed that almost 100% of the hemispherical radiometric signal originates within a distance of a few radiometer heights, while only about 50–80% of the cumulative aerodynamic signal is generated within a distance of about 20 sensor heights. In order to achieve comparable extensions of the footprint areas, hemispherical radiometric measurements should therefore be taken about 6–15 times higher than turbulent flux ones, depending on the vegetation type. The analysis also highlights that the size of the radiative footprint decreases at increasing leaf area index, whereas the aerodynamic footprint shows an opposite behavior. For the abovementioned reasons, this work may support the interpretation of energy flux measurements and the optimal design of eddy covariance stations located in heterogeneous sites.  相似文献   

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

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