首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
The L-band brightness temperature of natural grass fields is strongly influenced by rainfall interception. In wet conditions, the contribution of the soil, mulch, and vegetation to the overall microwave emission is difficult to decouple, thus rendering the retrieval of surface soil moisture from a direct emission model difficult. This paper investigates the development and assesses the performances of statistical regressions linking passive microwave measurements to surface soil moisture in order to assess the potential of soil moisture retrievals over natural grass. First, statistical regressions were analytically derived from the L-Band Emission of the Biosphere model (L-MEB). Single configuration (1 angle, 1 polarisation), and multi-configuration regressions (2 angles, or 2 polarisations) were developed. Second, the performance of statistical regressions was evaluated under different rainfall interception conditions. For that purpose, a modified polarisation ratio at L-band was used to build three data sets with different interception levels. In the presence of interception, a regression based on one observation angle (50°) and two polarisations was able to reduce the effects of vegetation and soil roughness on the soil moisture retrievals. The methodology presented in this study is also able to provide estimates of the vegetation and soil roughness contribution to the brightness temperature.  相似文献   

2.
This paper explores the effect of rain intercepted by vegetation and mulch on the L-band emission of natural grass. The study is based on radiometric, meteorological, and biophysical measurements obtained during the SMOSREX Experiment (Toulouse, France). Several approaches were followed to evaluate interception effects. Firstly, the analysis of microwave brightness temperature (TB) measurements at L-band indicated that interception increases vegetation emission at both polarisations. Secondly, the use of microwave indices to detect the presence/absence of interception was examined. In particular, a modified polarisation ratio at 50° was found to be well related to the interception status of the standing vegetation. Finally, the vegetation optical depth (τ), which parameterises the extinction across the vegetation layer, was retrieved from the TB observations. It was found that τ increases with the increase in the water content stored within the vegetation and mulch after rainfall. The study highlights the strong impact of intercepted water in otherwise weakly attenuating covers such as grasses. Interception might therefore be an issue to consider in order to improve soil moisture retrieval algorithms from L-band observations.  相似文献   

3.
Soil moisture will be mapped globally by the European Soil Moisture and Ocean Salinity (SMOS) mission to be launched in 2009. The expected soil moisture accuracy will be 4.0 %v/v. The core component of the SMOS soil moisture retrieval algorithm is the L-band Microwave Emission of the Biosphere (L-MEB) model which simulates the microwave emission at L-band from the soil-vegetation layer. The model parameters have been calibrated with data acquired by tower mounted radiometer studies in Europe and the United States, with a typical footprint size of approximately 10 m. In this study, aircraft L-band data acquired during the National Airborne Field Experiment (NAFE) intensive campaign held in South-eastern Australia in 2005 are used to perform the first evaluation of the L-MEB model and its proposed parameterization when applied to coarser footprints (62.5 m). The model could be evaluated across large areas including a wide range of land surface conditions, typical of the Australian environment. Soil moisture was retrieved from the aircraft brightness temperatures using L-MEB and ground measured ancillary data (soil temperature, soil texture, vegetation water content and surface roughness) and subsequently evaluated against ground measurements of soil moisture. The retrieval accuracy when using the L-MEB ‘default’ set of model parameters was found to be better than 4.0 %v/v only over grassland covered sites. Over crops the model was found to underestimate soil moisture by up to 32 %v/v. After site specific calibration of the vegetation and roughness parameters, the retrieval accuracy was found to be equal or better than 4.8 %v/v for crops and grasslands at 62.5-m resolution. It is suggested that the proposed value of roughness parameter HR for crops is too low, and that variability of HR with soil moisture must be taken into consideration to obtain accurate retrievals at these scales. The analysis presented here is a crucial step towards validating the application of L-MEB for soil moisture retrieval from satellite observations in an operational context.  相似文献   

4.
5.
In the framework of ESA's SMOS mission (Soil Moisture and Ocean Salinity), many studies have been carried out over different land surface types to model their microwave emission at L-band (1.4 GHz). Results of these studies have been integrated in the emission model L-MEB (L-Band Microwave Emission of the Biosphere), which is the core of the SMOS Level 2 soil moisture retrieval algorithm. The Mediterranean Ecosystem L-Band characterisation EXperiment (MELBEX-I) was carried out at the SMOS validation site near Valencia in autumn 2005. The main objective of MELBEX-I was to calibrate L-MEB over Mediterranean shrub land, as no data were available over this biome. For that purpose, multi-angular and dual polarimetric measurements (H, V) were obtained by the EMIRAD L-band radiometer from a 14-m tower. Results of this study indicate a small and constant impact of vegetation on the microwave emission of shrub land, and L-MEB parameters for shrub land were obtained. In addition, the study highlights the need for calibrating microwave soil roughness, which was found to be constant at the site. Depending on the number of retrieved parameters, soil moisture (SM) near the surface could be estimated with errors between 0.035 m3 m− 3 (if only SM was retrieved) and 0.057 m3 m− 3 (if SM, optical depth and a roughness parameter were simultaneously retrieved). Finally, no modelling improvements were observed when coarse estimates of the fraction of exposed rocks were accounted for in the model.  相似文献   

6.
土壤水分是地气间水热交换的重要变量,影响着地表感热潜热划分、水分收支和植被蒸腾等过程,青藏高原土壤水分的研究对于改进高原水分循环和能量平衡的模拟研究具有重要意义。随着SMOS、SMAP等卫星的发射,L波段被动微波遥感技术成为大尺度监测土壤水分的主要手段。分别从L波段星—机—地观测与微波辐射模拟、区域尺度土壤水分观测、卫星产品评估与土壤水分反演算法发展等方面系统回顾和总结了近年来L波段被动微波遥感及其土壤水分反演算法、产品在青藏高原的主要应用与研究进展。在此基础上,归纳了当前高原L波段被动微波辐射模拟与土壤水分反演存在的问题,主要包括缺乏高原尺度的微波辐射模拟评估和改进的卫星土壤水分产品、土壤冻结时期的水分监测产品依然缺失等问题。针对存在的问题,进一步提出了相关建议与展望,建议今后的研究应加强高原尺度的微波辐射模拟评估与土壤水分产品改进工作,并积极拓展土壤水分产品在高原水分循环和能量平衡模拟、植被生长与干旱监测的应用研究。  相似文献   

7.
Land surface characteristics: soil and vegetation and rainfall inputs are distributed in nature. Representation of land surface characteristics and inputs in models is lumped at spatial scales corresponding to the grid size or observation density. Complete distributed representation of these characteristics or inputs is infeasible due to excessive computational costs or costs associated with maintaining dense observational networks. The measurements of microwave brightness temperatures by the SSM/I (Special Sensor Microwave Imager) are at resolutions of the order of 56km 56km for 19 GHz and 33 km 33 km for 37 GHz. At these resolutions, soil moisture and vegetation are not homogeneous over the measurement area. The experiments carried out in this study determine the effect of heterogeneities in vegetation (leaf area index) and input rainfall on simulated soil moisture and brightness temperatures and the inversion of brightness temperatures to obtain soil moisture estimates. This study would help us to understand the implications of using the SSM/I microwave brightness temperatures for soil moisture estimation. The consequences of treating rainfall inputs and vegetation over large land surface areas in a lumped fashion is examined. Simpler methods based on dividing the leaf area index or input rainfall into classes rather than explicit representation for representing heterogeneities in leaf area index and spatial distribution of rainfall is tested. It is seen that soil moisture is affected by the representation (lumped vs distributed) of rainfall and not leaf area index. The effect of spatially distributed soil moisture on the inversion of observed SSM/I brightness temperatures to obtain soil moisture estimates is investigated. The inversion process does not exhibit biases in the retrieval of soil moisture. The methodology presented in this paper can be used for any satellite sensor for purposes of analysis and evaluation.  相似文献   

8.

A simple formulation relating the L-band microwave brightness temperature detected by a passive microwave radiometer to the near surface soil moisture was developed using MICRO-SWEAT, a coupled microwave emission model and soil-vegetation-atmosphere-transfer (SVAT) scheme. This simple model provides an ideal tool with which to explore the impact of sub-pixel heterogeneity on the retrieval of soil moisture from microwave brightness temperatures. In the case of a bare soil pixel, the relationship between apparent emissivity and surface soil moisture is approximately linear, with the clay content of the soil influencing just the intercept of this relationship. It is shown that there are no errors in the retrieved soil moisture from a bare soil pixel that is heterogeneous in soil moisture and texture. However, in the case of a vegetated pixel, the slope of the relationship between apparent emissivity and surface soil moisture decreases with increasing vegetation. Therefore for a pixel that is heterogeneous in vegetation and soil moisture, errors can be introduced into the retrieved soil moisture. Generally, under moderate conditions, the retrieved soil moisture is within 3% of the actual soil moisture. Examples illustrating this discussion use data collected during the Southern Great Plains '97 Experiment (SGP97).  相似文献   

9.
Data gathered during the NASA sponsored Multisensor Aircraft Campaign Hydrology (MACHYDRO) experiment in central Pennsylvania (U.S.A.) in July, 1990 have been analysed to study the combined use of active and passive microwave sensors for estimating soil moisture from vegetated areas. These data sets were obtained during an eleven-day period with NASA's Airborne Synthetic Aperture Radar (AIRSAR), and Push-Broom Microwave Radiometer (PBMR) over an instrumented watershed, which included agricultural fields with a number of different crop covers. Simultaneous ground truth measurements were also made in order to characterize the state of vegetation and soil moisture under a variety of meteorological conditions. Various multi-sensor techniques are currently under investigation to improve the accuracy of remote sensing estimates of the soil moisture in the presence of vegetation and surface roughness conditions using these data sets. One such algorithm involving combination of active and passive microwave sensors is presented here, and is applied to representative corn fields in the Mahantango watershed that was the focus of study during the MACHYDRO experiment. In this algorithm, a simple emission model is inverted to obtain Fresnel reflectivity in terms of ground and vegetation parameters. Since Fresnel reflectivity depends on soil dielectric constant, soil moisture is determined from reflectivity using dielectric-soil moisture relations. The algorithm requires brightness temperature, vegetation and ground parameters as the input parameters. The former is measured by a passive microwave technique and the later two are estimated by using active microwave techniques. The soil moisture estimates obtained by this combined use of active and passive microwave remote sensing techniques, show an excellent agreement with the in situ soil moisture measurements made during the MACHYDRO experiment.  相似文献   

10.
以青藏高原开展的L波段地基微波辐射(ELBARA-III型)综合观测试验为依据,基于 τ - ω 辐射传输模型评估了Wang-Schmugge、Mironov、Dobson和 Four-Phase 4种土壤介电模型对L波段微波亮温模拟及土壤湿度反演的影响。结果表明:相同植被和粗糙度参数化方案条件下,4种土壤介电模型对微波亮温模拟存在明显差异,当土壤湿度小于0.23 m3·m-3时,Wang-Schmugge模型与其他3种土壤介电模型微波亮温模拟结果差异最为显著,水平和垂直极化微波亮温模拟最大差值可达8.0 K和4.4 K;当模拟土壤湿度大于0.23 m3·m-3时,Four-phase模型模拟的微波亮温显著高于其他3种土壤介电模型模拟结果;当土壤湿度饱和时,4种土壤介电模型间水平和垂直极化微波亮温模拟最大差值约为6.1 K和4.8 K,且4种土壤介电模型对水平极化微波亮温模拟的差异比垂直极化模拟的差异更为显著。而基于4种介电模型的土壤湿度反演对比试验表明,水平极化条件下基于Wang-Schmugge模型反演土壤湿度,较其他参数化方案,能有效减轻反演土壤湿度对观测土壤湿度的低估,Mironov模型减轻了垂直极化条件下反演土壤湿度对观测值的高估程度。在现有 τ - ω 模型参数化方案的基础上,总结了4种土壤介电模型在青藏高原典型草地下垫面的适用性,将为星载L波段辐射计青藏高原土壤湿度反演应用提供客观的土壤介电模型方案选取依据。  相似文献   

11.
Based on soil moisture and freeze/thaw comprehensive experiments conducted at the north-eastern part of the Tibetan Plateau, the L-band brightness temperature, the in-situ soil moisture and temperature, vegetation leaf area index are measured simultaneously for the purpose of evaluating performances on forward brightness temperature simulation and soil moisture retrieval using four dielectric constant models, including Wang-Schmugge, Mironov, Dobson, and Four-Phase model. The forward brightness temperature simulations indicate that the difference of simulated brightness temperature between Wang schmugge model and the other three dielectric constant models is most significant at lower soil moisture content condition (soil moisture is less than 0.23 m3·m-3) , nevertheless, the difference of Mironov model simulation is most significant in contrast with the ones of other three models at higher soil moisture condition (soil moisture is greater than 0.23 m3·m-3). The practical retrieval of soil moisture from the ground-based radiometer measurements indicate that Wang-Schmugge model can effectively reduce the underestimation of soil moisture at the horizontal polarization, this resulted an improvement to the accuracy of retrieved soil moisture. Mironov model can reduce the underestimation of retrieved soil moisture at the vertical polarization. In accordance with a state-of-the-art parameterization scheme, the evaluation of performances of four dielectric constant models at the typical alpine meadow is potential for selecting optimum soil moisture retrieval by using soil dielectric model from space-borne L-band radiometer observation over the Tibetan Plateau  相似文献   

12.
Airborne L-band data from the Australian National Airborne Field Experiment 2005 (NAFE '05) field campaign were used to investigate the influence of fractional forest cover on soil moisture retrievals from heterogeneous (grass/forest) pixels. This study is, to our knowledge, the first to use experimental data on this subject and was done in view of the SMOS mission, in order to contribute to calibration/validation studies and the analysis of heterogeneous surfaces. Because the multi-angle observations were contained in swaths, swaths were used instead of pixels as the basic surface unit in this study. Simultaneous retrievals of soil moisture (SM) and vegetation optical depth (τNAD) were undertaken by inversion of the L-MEB zero-order radiative transfer model. This was done for two different retrieval configurations, the first consisting of swath-effective values of SM and τNAD and the second consisting of values of SM and τNAD for the non-forested (i.e. grass) fraction of the swath, with forest emission known from forward modelling. Model inputs for non-retrieved parameters were either default values taken from the literature or site- and time-specific values obtained from observations of nearby homogeneous swaths gathered during the same flight. The main focus of this study was on retrieval behaviour for various soil moisture conditions and forest fractions. Area-averaged retrieval results were generally very reasonable for both retrieval configurations. When retrieving swath-effective values of SM and τNAD, τNAD showed an increased overestimation with increased forest fraction. Highest retrieved values of SM were found at intermediate values of forest fraction. The results show the difficulty in flagging upper limits of pixel forest fraction during soil moisture retrievals, besides the fact that erroneous parameter values can lead to high errors in retrieved SM, especially in wet conditions. This study is the first to give a realistic idea of the errors and uncertainties involved in soil moisture retrievals from partly forested swaths, and as such will contribute to a better understanding of SMOS calibration/validation issues.  相似文献   

13.
The Soil Moisture and Ocean Salinity (SMOS) mission, launched in November 2009, provides global maps of soil moisture and ocean salinity by measuring the L-band (1.4 GHz) emission of the Earth's surface with a spatial resolution of 40-50 km. Uncertainty in the retrieval of soil moisture over large heterogeneous areas such as SMOS pixels is expected, due to the non-linearity of the relationship between soil moisture and the microwave emission. The current baseline soil moisture retrieval algorithm adopted by SMOS and implemented in the SMOS Level 2 (SMOS L2) processor partially accounts for the sub-pixel heterogeneity of the land surface, by modelling the individual contributions of different pixel fractions to the overall pixel emission. This retrieval approach is tested in this study using airborne L-band data over an area the size of a SMOS pixel characterised by a mix Eucalypt forest and moderate vegetation types (grassland and crops), with the objective of assessing its ability to correct for the soil moisture retrieval error induced by the land surface heterogeneity. A preliminary analysis using a traditional uniform pixel retrieval approach shows that the sub-pixel heterogeneity of land cover type causes significant errors in soil moisture retrieval (7.7%v/v RMSE, 2%v/v bias) in pixels characterised by a significant amount of forest (40-60%). Although the retrieval approach adopted by SMOS partially reduces this error, it is affected by errors beyond the SMOS target accuracy, presenting in particular a strong dry bias when a fraction of the pixel is occupied by forest (4.1%v/v RMSE, −3.1%v/v bias). An extension to the SMOS approach is proposed that accounts for the heterogeneity of vegetation optical depth within the SMOS pixel. The proposed approach is shown to significantly reduce the error in retrieved soil moisture (2.8%v/v RMSE, −0.3%v/v bias) in pixels characterised by a critical amount of forest (40-60%), at the limited cost of only a crude estimate of the optical depth of the forested area (better than 35% uncertainty). This study makes use of an unprecedented data set of airborne L-band observations and ground supporting data from the National Airborne Field Experiment 2005 (NAFE'05), which allowed accurate characterisation of the land surface heterogeneity over an area equivalent in size to a SMOS pixel.  相似文献   

14.
Effective soil moisture sampling depth of L-band radiometry: A case study   总被引:1,自引:0,他引:1  
The aim of this study is to analyze the influence of the soil moisture sampling depth in the parameterization of soil emission in microwave radiometry at L-band. The analysis is based on brightness temperature, soil moisture and temperature measurements acquired over a bare soil during the SMOSREX experiment. A more detailed profile of surface soil moisture was obtained with a soil heat and water flows mechanistic model. It was found that (1) the soil moisture sampling depth depends on soil moisture conditions, (2) the effective soil moisture sampling depth is shallower than provided by widely used field moisture sensors, and (3) the soil moisture sampling depth has an impact on the calibration of soil roughness model parameters. These conclusions are crucial for the calibration and validation of remote sensing data at L-band.  相似文献   

15.
Surface soil moisture is a key variable used to describe water and energy exchanges at the land surface/atmosphere interface. Passive microwave remotely sensed data have great potential for providing estimates of soil moisture with good temporal repetition on a daily basis and on a regional scale (∼10 km). However, the effects of vegetation cover, soil temperature, snow cover, topography, and soil surface roughness also play a significant role in the microwave emission from the surface. Different soil moisture retrieval approaches have been developed to account for the various parameters contributing to the surface microwave emission. Four main types of algorithms can be roughly distinguished depending on the way vegetation and temperature effects are accounted for. These algorithms are based on (i) land cover classification maps, (ii) ancillary remote sensing indexes, and (iii) two-parameter or (iv) three-parameter retrievals (in this case, soil moisture, vegetation optical depth, and effective surface temperature are retrieved simultaneously from the microwave observations). Methods (iii) and (iv) are based on multiconfiguration observations, in terms of frequency, polarization, or view angle. They appear to be very promising as very few ancillary information are required in the retrieval process. This paper reviews these various methods for retrieving surface soil moisture from microwave radiometric systems. The discussion highlights key issues that will have to be addressed in the near future to secure operational use of the proposed retrieval approaches.  相似文献   

16.
被动微波遥感反演土壤水分进展研究   总被引:15,自引:2,他引:13  
在地球系统中, 地表土壤水分是陆地和大气能量交换过程中的重要因子, 并对陆地表面蒸散、水的运移、碳循环有很强的控制作用, 大面积监测土壤水分在水文、气象和农业科学领域具有较大的应用潜力。被动微波遥感是监测土壤含水量最有效的手段之一, 相比红外与可见光, 它具有波长长, 穿透能力强的优势, 相比主动微波雷达, 被动微波辐射计具有监测面积大、周期短, 受粗糙度影响小, 对土壤水分更为敏感, 算法更为成熟的优势。然而微波辐射计观测到的亮温除了受土壤水分影响外, 还要考虑如植被覆盖、土壤温度、雪覆盖以及地形、地表粗糙度、土壤纹理和大气效应以及地表的异质性等其它因子的影响。目前, 已研究出许多使用被动微波辐射计反演土壤水分的方法,这些方法大部分是围绕着土壤湿度与亮温温度之间的关系进行, 同时也考虑其它各种不同因子对 地表微波辐射的影响。从介绍被动微波反演地表参数的原理入手, 重点介绍被动遥感反演土壤水分当前的算法进展、研究趋势等。  相似文献   

17.
土壤水分是陆地生态系统和水循环的重要状态变量,在植被生长监测、农作物产量评估等研究中均发挥着重要作用。为了消除植被散射的影响,进而实现农田地表土壤水分的高精度反演,以时间序列Sentinel-1影像及MODIS产品为实验数据,基于高级积分方程模型和比值植被模型的耦合模型,通过采用不同光学植被参数和VH交叉极化后向散射系数,分别对农田植被散射贡献进行表征,消除植被散射的影响,进而实现土壤水分的高精度反演。结果表明:当利用VH极化进行参数化植被散射贡献时,标定的耦合模型,虽然可消除对光学植被参数的依赖并较好地模拟Sentinel-1卫星观测,但土壤水分反演结果效果欠理想,相关系数R最大仅为0.54;与VH极化相比,利用光学植被参数表征植被散射贡献时,土壤水分整体反演效果较理想,R最大达到0.79,但光学植被参数反演结果在不同站点存在显著的空间差异性,R介于0.07~0.79之间。因此,在未来研究中可尝试将雷达数据与光学数据协同反演,以期在消除植被散射影响的基础上,实现植被覆盖区域土壤水分的高精度反演及动态变化监测。  相似文献   

18.
Abstract

Microwave radiometer measurements of soil moisture content were made over bare and vegetated fields with dual polarized microwave radiometers at 1·55GHz (L-band) and 19·1 GHz (K.-band). Two typical Indian crops Bazra and Gawar have been studied. The bare field measurements were used to investigate the effect of soil texture on sensitivity of a radiometer to soil moisture and for soil moisture sampling depth. It is found that expression of soil moisture as available moisture content in the soil can minimize the texture effect. The estimated soil moisture sampling depth for L-band is 2-5 cm, while for K-band it is less than 2 cm. The vegetation cover affects the sensitivity of the radiometer to soil moisture. This effect is more pronounced the denser the vegetation and higher the frequency of observation. The measured polarization factor over a vegetated field at L-band was found to be appreciably reduced compared to that over a bare field. The difference between normalized brightness temperature from L-band and K-band is sensitive to vegetation type. The soil moisture under vegetation cover at L-band can be predicted well using Jackson's parametric model.  相似文献   

19.
Microwave radiometric measurements over bare fields of different surface roughness were made at frequencies of 1.4 GHz, 5 GHz, and 10.7 GHz to study the frequency dependence, as well as the possible time variation, of surface roughness. An increase in surface roughness was found to increase the brightness temperature af soils and reduce the slope of regression between brightness temperature and soil moisture content. The frequency dependence of the surface roughness effect was relatively weak when compared with that of the vegetation effect. Radiometric time-series observations over a given field indicate that field surface roughness might gradually diminish with time, especially after a rainfall or irrigation. The variation of surface roughness increases the uncertainty of remote soil moisture estimates by microwave radiometry. Three years of radiometric measurements over a test site revealed a possible inconsistency in the soil bulk density determination, which is an important factor in the interpretation of radiometric data.  相似文献   

20.
This paper investigates vegetation effects at L-band by using a first-order radiative transfer (RT) model and truck-based microwave measurements over natural conifer stands to assess the applicability of the τ ? ω (tau–omega) model over trees. The tau–omega model is a zero-order RT solution that accounts for vegetation effects with two vegetation parameters (vegetation opacity and single-scattering albedo), which represent the canopy as a whole. This approach inherently ignores multiple-scattering effects and, therefore, has a limited validity depending on the level of scattering within the canopy. The fact that the scattering from large forest components such as branches and trunks is significant at L-band requires that zero-order vegetation parameters be evaluated (compared) along with their theoretical definitions to provide a better understanding of these parameters in the retrieval algorithms as applied to trees.This paper compares the effective vegetation opacities, computed from multi-angular pine tree brightness temperature data, against the results of two independent approaches that provide theoretical and measured optical depths. These two techniques are based on forward scattering theory and radar corner reflector measurements, respectively. The results indicate that the effective vegetation opacity values are smaller than but of similar magnitude to both radar and theoretical estimates. The effective opacity of the zero-order model is thus set equal to the theoretical opacity and an explicit expression for the effective albedo is then obtained from the zero- and first-order RT model comparison. The resultant albedo is found to have a similar magnitude as the effective albedo value obtained from brightness temperature measurements. However, both are less than half of the single-scattering albedo estimated using the theoretical calculations (0.5?0.6 for tree canopies at L-band). This lower observed effective albedo balances the scattering darkening effect of the large theoretical single-scattering albedo with a first-order multiple-scattering contribution. The retrieved effective albedo is different from theoretical definitions and not the albedo of single forest elements anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy, including multiple-scattering and canopy ground interaction.  相似文献   

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

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