首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Summary Long-term flux measurement sites are often characterized by a heterogeneous terrain, which disagrees with the fundamental theoretical assumptions for eddy-covariance measurements. An evaluation procedure to assess the influence of terrain heterogeneity on the data quality has been developed by G?ckede et al. (2004), which combines existing quality assessment tools for flux measurements with analytic footprint modeling. In addition to micrometeorological input data, this approach requires information defining the land use structure and the roughness of the surrounding terrain. The aim of this study was to improve the footprint based site evaluation approach by using high-resolution land use maps derived by Landsat ETM+ and ASTER satellite data. The influence of the grid resolution of the maps on the results was examined, and four different roughness length classification schemes were tested. Due to numerical instabilities of the analytic footprint routine, as an additional footprint model a Lagrangian stochastic footprint routine (Rannik et al., 2003) was employed. Application of the approach on two German FLUXNET sites revealed only weak influence of the characteristics of the land use data when the land use structure was homogeneous. For a more heterogeneous site, use of the more detailed land use maps derived by remote sensing methods resulted in distinct differences indicating the potential of remote sensing for improving the flux measurement site evaluation.  相似文献   

2.
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.  相似文献   

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

4.
Summary Measuring turbulent fluxes with the eddy covariance method has become a widely accepted and powerful tool for the determination of long term data sets for the exchange of momentum, sensible and latent heat, and trace gases such as CO2 between the atmosphere and the underlying surface. Several flux networks developed continuous measurements above complex terrain, e.g. AmeriFlux and EUROFLUX, with a strong focus on the net exchange of CO2 between the atmosphere and the underlying surface. Under many conditions basic assumptions for the eddy covariance method in its simplified form, such as stationarity of the flow, homogeneity of the surface and fully developed turbulence of the flow field, are not fulfilled. To deal with non-ideal conditions which are common at many FLUXNET sites, quality tests have been developed to check if these basic theoretical assumptions are valid.In the framework of the CARBOEUROFLUX project, we combined quality tests described by Foken and Wichura (1996) with the analytical footprint model of Schmid (1997). The aim was to identify suitable wind sectors and meteorological conditions for flux measurements. These tools were used on data of 18 participating sites. Quality tests were applied on the fluxes of momentum, sensible and latent heat, and on the CO2-flux, respectively. The influence of the topography on the vertical wind component was also checked. At many sites the land use around the flux towers is not homogeneous or the fetch may not be large enough. So the relative contribution of the land use type intended to be measured was also investigated. Thus the developed tool allows comparative investigations of the measured turbulent fluxes at different sites if using the same technique and algorithms for the determination of the fluxes as well as analyses of potential problems caused by influences of the surrounding land use patterns.  相似文献   

5.
We measured the surface energy budget of an Alpine grassland in highly complex terrain to explore possibilities and limitations for application of the eddy-covariance technique, also for CO2 flux measurements, at such non-ideal locations. This paper focuses on the influence of complex terrain on the turbulent energy measurements of a characteristic high Alpine grassland on Crap Alv (Alp Weissenstein) in the Swiss Alps during the growing season 2006. Measurements were carried out on a topographic terrace with a slope of 25 inclination. Flux data quality is assessed via the closure of the energy budget and the quality flag method used within the CarboEurope project. During 93% of the time the wind direction was along the main valley axis (43% upvalley and 50% downvalley directions). During the transition times of the typical twice daily wind direction changes in a mountain valley the fraction of high and good quality flux data reached a minimum of ≈50%, whereas during the early afternoon ≈70% of all records yielded good to highest quality (CarboEurope flags 0 and 1). The overall energy budget closure was 74 ± 2%. An angular correction for the shortwave energy input to the slope improved the energy budget closure slightly to 82 ± 2% for afternoon conditions. In the daily total, the measured turbulent energy fluxes are only underestimated by around 8% of net radiation. In summary, our results suggest that it is possible to yield realistic energy flux measurements under such conditions. We thus argue that the Crap Alv site and similar topographically complex locations with short-statured vegetation should be well suited also for CO2 flux measurements.  相似文献   

6.
In recent years considerable effort has been focused on combining micrometeorological and stable isotope techniques to partition net fluxes and to study biosphere–atmosphere exchange processes. While much progress has been achieved over the last decade, some new issues are beginning to emerge as technological advances, such as laser spectroscopy, permit isotopic fluxes to be measured more easily and continuously in the field. Traditional investigations have quantified the isotopic composition of biosphere-atmosphere exchange by using the Keeling two-member mixing model (the classic Keeling plot). An alternative method, based on a new capacity to measure isotopic mixing ratios, is to determine the isotope composition of biosphere–atmosphere exchange from the ratio of flux measurements. The objective of this study was to critically evaluate these methods for quantifying the isotopic composition of ecosystem respiration (δR) over a period of three growing seasons (2003–2005) within a heterogeneous landscape consisting of C3 and C4 species. For C4 canopies, the mixing model approach produced δR values that were 4–6‰ lower (isotopically lighter) than the flux-gradient method. The analyses presented here strongly suggest that differences between flux and concentration footprint functions are the main factor influencing the inequality between the mixing model and flux-gradient approaches. A mixing model approach, which is based on the concentration footprint, can have a source area influence more than 20-fold greater than the flux footprint. These results highlight the fact that isotopic flux partitioning is susceptible to problems arising from combining signals (concentration and fluxes) that represent very different spatial scales (footprint). This problem is likely to be most pronounced within heterogeneous terrain. However, even under ideal conditions, the mismatch between concentration and flux footprints could have a detrimental impact on isotopic flux partitioning where very small differences in isotopic signals must be resolved.  相似文献   

7.
The scaled standard deviations of temperature and humidity are investigated in complex terrain. The study area is a steep Alpine valley, with six measurement sites of different slope, orientation and roughness (i-Box experimental site, Inn Valley, Austria). Examined here are several assumptions forming the basis of Monin–Obukhov similarity theory (MOST), including constant turbulence fluxes with height and the degree of self-correlation between the involved turbulence variables. Since the basic assumptions for the applicability of the MOST approach—horizontally homogeneous and flat conditions—are violated, the analysis is performed based on a local similarity hypothesis. The scaled standard deviations as a function of local stability are compared with previous studies from horizontally homogeneous and flat terrain, horizontally inhomogeneous and flat terrain, weakly inhomogeneous and flat terrain, as well as complex terrain. As a reference, similarity relations for unstable and stable conditions are evaluated using turbulence data from the weakly inhomogeneous and flat terrain of the Cabauw experimental site in the Netherlands, and assessed with the same post-processing method as the i-Box data. Significant differences from the reference curve and also among the i-Box sites are noted, especially for data derived from the i-Box sites with steep slopes. These differences concern the slope and the magnitude of the best-fit curves, illustrating the site dependence of any similarity theory.  相似文献   

8.
非均匀下垫面湍流通量观测的印痕分析   总被引:5,自引:1,他引:4  
彭谷亮  刘绍民  蔡旭晖 《大气科学》2008,32(5):1064-1070
利用2005年在北京昌平区小汤山开展的非均匀下垫面观测实验的观测数据,分别处理了涡旋相关法和大孔径闪烁仪LAS测得的感热通量,并运用印痕模型对数据进行“源区”分析。分析结果表明:在复杂地表上,涡旋相关的测量值由于混杂了其他下垫面的通量信息,不能真实地反映其观测区域的湍流特征,存在着观测误差。而LAS观测出现的低估问题则与掺混高度和“源区”有关。通过用印痕方法进行数据订正后,以上问题均得到很好的改善,订正结果与LAS观测通量的线性关系良好,两者的相关系数达到0.9。实验结果验证了印痕模型在非均匀下垫面的适用性。  相似文献   

9.
The flux footprint, that is the contribution per unit emission from each element of the upwind surface area to measurement of the vertical flux of a passive scalar, is calculated for fluxes estimated by micrometeorological profile techniques. It is found that the upwind extent of the footprint for concentration-profile flux estimates is similar to that of the footprint for eddy-covariance flux measurements, when the eddy-covariance measurement is made at a height equal to the arithmetic mean of the highest and lowest profile measurement heights for stable stratification or the geometric mean for unstable stratification. The concentration-profile flux footprint depends on the ratio of the highest to the lowest measurement height, but is insensitive to the number of measurement levels. The concentration-profile flux footprint extends closer to the measurement location than does the 'equivalent eddy-covariance flux footprint, and the difference becomes more pronounced as the ratio of the profile measurement heights increases. The flux footprint for the Bowen-ratio technique is identical to that for a two-level profile measurement only for very limited circumstances. In the more general case, a flux footprint cannot be defined for the Bowen-ratio technique and the uniform upwind fetch required for representative flux measurements depends on the specific spatial distribution of surface fluxes.  相似文献   

10.
The measurement of scalar fluxes employing the eddy covariance method is a widely used experimental approach,for which the flow distortion due to obstacles (e.g., sensor mounts and mast)is a well-known but not fully solved problem. In order to reduce flow distortion we installed a sonic anemometer in a surface-normal orientationrelative to the terrain slope, and a second instrument in a verticalposition at a horizontal distance of 1.54 m from the first instrumentWe found a significant reduction in the rotation angle necessary for the coordinaterotation procedure in the x-z plane whencomputing 30-minute flux averages with the surface-normal orientation. In 91% of all cases this rotation angleremained within the angle of incidence of ±10° recommended bythe manufacturer. In contrast, only 24% of the measurements taken with the vertically mounted anemometer were obtained at an angle of incidencewithin ±10°, and 3% were outside the ±30° range specified for an acceptable operation.A data quality test based on the variance of vertical windspeed normalized with friction velocity (w/u*) revealed problems for application under stable conditions due to large uncertainties in the determination of the Monin–Obukhov stability parameter z/L. An alternative test using the bulk drag coefficient CD revealed other problems related to the dependence of CD on z/z0, the measuring height normalized by the roughness length, which do not appear to be constantin complex terrain. With both tests, a tendency for a slightly improved dataquality was found for the surface normal set-up, which, however, proved statistically insignificant.It is concluded that the surface-normal set-up of a sonic anemometer significantly reduces flow distortion by thesensor head. Although the surface-normal mounting position therefore appears to be the preferred one, with decreased flow distortion and a slightly improved data quality, no significant differences in turbulent quantities were found between the two set-uppositions. Hence, the consequences for short-term measurements of massand energy fluxes with a surface-normal set-up in complex terrain appearto be relevant only if single flux events are to be inspected, while for long-term measurements of integrated fluxes both the surface-normaland vertical installation of the sonic anemometer are adequate,indicating that eddy covariance measurements in complex terrain are lessdelicate than expected.  相似文献   

11.
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.  相似文献   

12.
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.  相似文献   

13.
青藏高原地理环境复杂,已有大气陆面-边界层研究工作多集中于不同下垫面,很少有对复杂地形区的研究。本文利用青藏高原东南部林芝地区2013年5月20日至7月9日四个野外试验站点的观测资料,分析了不同天气条件下,高原复杂地形区不同下垫面的陆-气能量交换特征。结果表明:在各站向下短波辐射基本一致的情况下,地形较陡的北坡阔叶林站感热通量远大于其他3个站点;下垫面植被覆盖最多的南面麦田站潜热通量最大。各站能量通量有明显的日变化特征,晴天时,感热通量和净辐射明显大于阴雨天,而潜热通量随天气状况变化不大。青藏高原复杂地形环境比不同天气条件对于感热通量的影响更显著;不同地形阴雨天时对于潜热通量有明显的影响。当南亚季风槽前的西南暖湿气流影响到林芝地区时,该地区以阴雨天为主,反之则以晴天为主。林芝地区地-气通量的月内变化明显受南亚季风活动的影响。  相似文献   

14.
Processing and quality control of flux data during LITFASS-2003   总被引:1,自引:0,他引:1  
Different aspects of the quality assurance and quality control (QA/QC) of micrometeorological measurements were combined to create a comprehensive algorithm which was then applied to experimental data from LITFASS-2003 (Lindenberg Inhomogeneous Terrain—Fluxes between Atmosphere and Surface: a long term Study). Eddy-covariance measurements of the latent heat flux were the main focus of the QA/QC efforts. The results of a turbulence sensor intercomparison experiment showed deviations between the different eddy-covariance systems on the order of 15%, or less than 30 W m−2, for the latent heat flux and 5%, or less than 10 W m−2, for the sensible heat flux. In order to avoid uncertainties due to the post-processing of turbulence data, a comprehensive software package was used for the analysis of experimental data from LITFASS-2003, including all necessary procedures for corrections and quality control. An overview of the quality test results shows that for most of the days more than 80% of the available latent heat flux data are of high quality so long as there are no instrumental problems. The representativeness of a flux value for the target land-use type was analysed using a stochastic footprint model. Different methods to calculate soil heat fluxes at the surface are discussed and a sensitivity analysis is conducted to select the most robust method for LITFASS-2003. The lack of energy balance closure, which was found for LITFASS-2003, can probably be attributed to the presence of low-frequency flux contributions that cannot be resolved with an averaging time of 30 min. Though the QA/QC system has been developed for the requirements of LITFASS-2003, it can also be applied to other experiments dealing with similar objectives.  相似文献   

15.
16.
Effects of complex terrain on net surface longwave radiation in China   总被引:1,自引:0,他引:1  
Net surface longwave radiation (NSLR) is one of key meteorological factors and is strongly influenced by cloud cover, surface temperature, humidity, and local micrometeorological conditions as well as terrain conditions. Realistically estimating NSLR is vitally important for understanding surface radiation balance and investigating micrometeorological factors of air pollution dispersion, especially in regions with complicated terrain. In this study, we proposed a distributed model for estimating NSLR by considering effects of complex local terrain conditions in China. Meteorological data (including mean temperature, relative humidity, and sunshine percentage) and observed NSLR data from 1993 to 2001 together with the digital elevation model data were used to parametrize the model and account for the effects of atmospheric factors and surface terrain factors according to the isotropic principle. The monthly NSLR during 1961–2000 was estimated at a spatial resolution of 1 km. Topographic analysis suggests that the distribution characteristics of NSLR with elevation or slope are consistent with those of field observations. In particular, the estimated NSLR is favorably comparable with site-level observations on the Tibetan Plateau (average relative error < 11%). Our results indicate that this model can describe microscale distribution features in mountainous areas in detail and that this improved approach can be used for NSLR spatial estimation in other regions with complicated terrain.  相似文献   

17.
Evaluation of six parameterization approaches for the ground heat flux   总被引:2,自引:1,他引:2  
Summary There are numerous approaches to the parameterization of the ground heat flux that use different input data, are valid for different times of the day, and deliver results of different quality. Six of these approaches are tested in this study: three approaches calculating the ground heat flux from net radiation, one approach using the turbulent sensible heat flux, one simplified in situ measurement approach, and the force-restore method. On the basis of a data set recorded during the LITFASS-2003 experiment, the strengths and weaknesses of the approaches are assessed. The quality of the best approaches (simplified measurement and force-restore) approximates that of the measured data set. An approach calculating the ground heat flux from net radiation and the diurnal amplitude of the soil surface temperature also delivers satisfactory daytime results. The remaining approaches all have such serious drawbacks that they should only be applied with care. Altogether, this study demonstrates that ground heat flux parameterization has the potential to produce results matching measured ones very well, if all conditions and restrictions of the respective approaches are taken into account.  相似文献   

18.
应用印痕函数研究测风塔资料代表性   总被引:3,自引:0,他引:3       下载免费PDF全文
影响测风塔代表性的因素很多,为了合理地确定测风塔观测结果能够代表的区域范围,需要通过适当的方法综合考虑这些因素。湍流通量代表性研究中常用的印痕函数是风速、稳定度、下垫面粗糙度等的函数,可以定量地刻画单点观测结果的代表性。利用印痕函数方法对中国风能资源观测网中来自不同下垫面的44个测风塔2010年的观测结果进行了测风塔资料的代表性分析。首先,根据两层(10、70 m)的风速和温度计算稳定度参数Lp并确定中性层结,粗糙度的计算选取了中性层结条件下满足对数风廓线的资料。然后,利用中性层结下70 m高度主导风向和风向频率较大方位的资料分析了两种印痕函数方法(Schuepp和Hsieh方法)的适用性及代表性范围与粗糙度的关系。两种方法得到的印痕函数对应距离(即代表性范围)比较接近。但是,在复杂山地Schuepp方法已不适用。当粗糙度小于10 m时,粗糙度与印痕函数对应距离的关系可用对数-线性函数描述。根据经验推荐了利用印痕函数积分90%对应的距离以及下垫面显著变化位置来确定测风塔代表范围。但是,印痕函数积分阈值仍需要通过示踪实验或数值模拟检验。  相似文献   

19.
涡动相关仪和大孔径闪烁仪观测通量的空间代表性   总被引:13,自引:3,他引:10       下载免费PDF全文
在对涡动相关仪和大孔径闪烁仪足迹模型进行敏感性分析的基础上,利用北京密云站2006年8月至2007年12月期间的观测资料,应用足迹模型对观测通量的空间代表性做了初步的分析.结果表明:涡动相关仪和大孔径闪烁仪的源区对风向、Monin-Obukhov长度,空气动力学粗糙度和观测高度/有效高度等因子比较敏感.密云站涡动相关仪的源区白天主要分布在仪器的西南与南面,而夜间则在东北与北面.大孔径闪烁仪的源区为西南一东北向分布.涡动相关仪各月源区形状不同,但大致分布在南北长1000 m、东两宽850 m的范嗣内,而LAS各月源区为西南一东北向分布,长约2060 m,最宽处约为620 m.对涡动相关仪通量有贡献的下垫面主要为园地(67%)和耕地(19%).其中园地的通量贡献比例在夏、秋季比较大,冬、春季稍小,而耕地则相反.大孔径闪烁仪的主要通量贡献源区为园地、耕地和居民地,通量贡献比例分别为49%,28%和11%.其中园地和耕地通量贡献率的变化趋势与涡动相关仪的观测结果一致,但没有涡动相关仪的变化明显.  相似文献   

20.
Turbulence data collected in an area of three-dimensional complex terrain using instruments atteched to the tether cable of a captive balloon together with radiosonde ascents are presented. In addition, data collected using only radiosonde ascents in an area of two-dimensional complex terrain of large slope are also shown. Eddy correlation measurements of the turbulent momentum flux and wind velocity profiles are used to deduce the magnitude of the effective roughness from the drag coefficient and normalised velocity profiles. A relationship connecting the terrain characteristics and the roughness length is compared with the experimental data for both types of terrain plus previous experimental estimates of the roughness length over complex terrain. The formula taken from previous work by Grant and Mason (1990) is found to agree with the data when representing an area of order 100 km2.  相似文献   

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

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