首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
High-frequency temperature data were recorded at one height and they were used in Surface Renewal (SR) analysis to estimate sensible heat flux during the full growing season of two rice fields located north–northeast of Colusa, CA (in the Sacramento Valley). One of the fields was seeded into a flooded paddy and the other was drill seeded before flooding. To minimize fetch requirements, the measurement height was selected to be close to the maximum expected canopy height. The roughness sub-layer depth was estimated to discriminate if the temperature data came from the inertial or roughness sub-layer. The equation to estimate the roughness sub-layer depth was derived by combining simple mixing-length theory, mixing-layer analogy, equations to account for stable atmospheric surface layer conditions, and semi-empirical canopy–architecture relationships. The potential for SR analysis as a method that operates in the full surface boundary layer was tested using data collected over growing vegetation at a site influenced by regional advection of sensible heat flux. The inputs used to estimate the sensible heat fluxes included air temperature sampled at 10 Hz, the mean and variance of the horizontal wind speed, the canopy height, and the plant area index for a given intermediate height of the canopy. Regardless of the stability conditions and measurement height above the canopy, sensible heat flux estimates using SR analysis gave results that were similar to those measured with the eddy covariance method. Under unstable cases, it was shown that the performance was sensitive to estimation of the roughness sub-layer depth. However, an expression was provided to select the crucial scale required for its estimation.  相似文献   

2.
在TOGA-COAREIOP期间用涡度相关法对海气热通量进行了船载直接观测.对垂直风速、温度和湿度湍流脉动观测数据的谱分析显示它们在高频区基本满足"-2/3次方律".对船体简谐震荡影响的讨论从理论上证明该影响在热通量计算中可被忽略.根据以此方法得到的通量求出了中性层结条件下感热和潜热的整体输送系数分别为2.25×10-3和1.26×10-3.对海气边界层特性的分析表明该海域的近海层主要呈中性或弱不稳定层结.海气通量的变化与背景环流形势密切相关,潜热通量主要受海面风场强度的影响,而感热通量变化除了风场的影响外,层结变化也是一个重要因素.用整体输送法计算TOGA-COAREIOP期间以及TOGA期间8个航次的通量结果而得到的Bowen比约为0.1,显示潜热通量是暖池大气的主要热源.  相似文献   

3.
基于流体动力学的探空仪GTS1湿度测量误差修正研究   总被引:2,自引:0,他引:2       下载免费PDF全文
准确的湿度垂直廓线对短期数值天气预报和中长期气候研究分析均具有重要参考价值,但太阳辐射与高空云雨等因素导致探空湿度测量精度较低.针对探空仪GTS1白天进行湿度测量时存在明显偏干误差的现象,本文建立了GTS1湿度测量结构模型,采用计算流体动力学方法重点研究了太阳辐射对GTS1湿度测量的影响,同时探讨了其防雨/防辐射罩的防雨效果.结合改进的太阳辐射计算模型,利用Fluent模拟分析了防雨/防辐射罩云雨滴轨迹和内部温度场分布,数值仿真结果表明:在垂直气流下防雨/防辐射罩基本保护湿度传感器免受高空云雨的影响,对太阳辐射也起到一定防护作用.但太阳辐射对GTS1湿度测量的偏干影响仍较显著,海拔30km左右相对误差甚至高达70%.通过与RS92及GTS1探空湿度测量中的太阳辐射偏干误差相对比,发现基于流体动力学数值计算的湿度相对误差随海拔高度的变化趋势与之基本吻合.本文也给出了理论计算与实验测量存在一定差异的原因,进而以南京探空站2014年不同季节的温度、气压、气球的上升速度等探空数据为计算参数,获取了太阳辐射误差与海拔及环境温度的关系,并对四组探空湿度廓线进行太阳辐射误差修正,提升了GTS1探空仪高空湿度测量的精度.本文研究结果为探空湿度测量太阳辐射偏干误差的有效修正提供了一定的参考.  相似文献   

4.
Summary An inexpensive and simple form of radiometer capable of attachment to the American Bureau of Standards type radiosonde is described: it measures the upward flux of longwave radiation when used on night ascents. A slight modification to the telemetry circuit of the standard radiosonde enables the radiation signals to be obtained in addition to the usual pressure, temperature and humidity information.Deceased.  相似文献   

5.
利用南极中山站至Dome A考察断面上3个自动气象站2005~2007年的观测资料和2008年夏季在中山站附近冰盖获取的湍流观测资料,应用空气动力学方法和涡动相关法计算分析了中山站至Dome A断面上近地层各种湍流参数(感热通量,潜热通量,湍流温度、湿度和速度尺度,地表粗糙度,大气稳定度及动量输送系数)的季节变化、日变...  相似文献   

6.
The dissipation method, the method preferred for estimating scalar surface fluxes over open water has not traditionally been used by agronomists, whereas the surface renewal (SR) theory in conjunction with the analysis of the scalar time trace offers tremendous advantages for estimating fluxes over agronomic crops. For a steady and horizontally homogeneous flow, it is shown that the dissipation method and SR analysis are closely related. As a consequence, a new dissipation–SR analysis expression for estimating scalar surface fluxes was derived. The new equation requires no calibration, and the scalar time trace measured at a frequency capable of identifying canopy‐scale coherent structures (typically 4–10 Hz in agriculture) is the only input required. Sensible and latent heat flux estimates obtained from 10 Hz air temperature and water vapour concentration measurements in the inertial sub‐layer (2 m height) over short, homogeneous rangeland grass at a site where similarity does not hold gave similar results to those measured with the eddy covariance (EC) method. For unstable cases, the new equation provided a root mean square error of 57 W m?2 for the surface energy‐balance closure. For stable cases, the performance was difficult to evaluate because the EC fluxes were similar in magnitude to the sensor error. It is concluded that the proposed method can contribute to a better understanding of hydrological processes and water requirements by providing an accurate, less costly, alternative method to indirectly estimate evapotranspiration as the residual of the energy balance equation. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

7.
Surface soil heat flux is a component of surface energy budget and its estimation is needed in land-atmosphere interaction studies. This paper develops a new simple method to estimate soil heat flux from soil temperature and moisture observations. It gives soil temperature profile with the thermal diffusion equation and, then, adjusts the temperature profile with differences between observed and computed soil temperatures. The soil flux is obtained through integrating the soil temperature profile. Compared with previous methods, the new method does not require accurate thermal conductivity. Case studies based on observations, synthetic data, and sensitivity analyses show that the new method is preferable and the results obtained with it are not sensitive to the availability of temperature data in the topsoil. In addition, we pointed out that the soil heat flux measured with a heat-plate can be quite erroneous in magnitude though its phase is accurate.  相似文献   

8.
The two commonly used statistical measures of the air-sea heat flux, the sampling and classical means, have been compared using hourly reports over a 7-year-period from a weather ship stationed in the NE Atlantic. The sampling mean is the average over all flux estimates in a given period, where individual flux estimates are determined from ship reports of meteorological variables using the well-known bulk formulae. The classical mean is the flux derived by substituting period-averaged values for each of the meteorological variables into the bulk formula (where the averaging period employed is the same as that over which the fluxes are to be determined). Monthly sampling and classical means are calculated for the latent and sensible heat fluxes. The monthly classical mean latent heat flux is found to overestimate the sampling mean by an amount which increases from 1–2 W m−2 in summer to 7 W m−2 in winter, on average, over the 7-year-period. In a given winter month, the excess may be as great as 15 W m−2, which represents about 10% of the latent heat flux. For the sensible heat flux, any seasonal variation between the two means is of the order of 1 W m−2 and is not significant compared to the interannual variation. The discrepancy between the two means for the latent heat flux is shown to arise primarily from a negative correlation between the wind speed and sea-air humidity difference, the effects of which are implicitly included in the sampling method but not in the classical. The influence of the dominant weather conditions on the sign and magnitude of this correlation are explored, and the large negative values that it takes in winter are found to depend on the typical track of the mid-latitude depressions with respect to the position sampled. In conclusion, it is suggested that sampling means should be employed where possible in future climatological studies.  相似文献   

9.
A combined Raman–Rayleigh lidar has been designed at Chung-Li, Taiwan for the simultaneous measurement of water-vapor mixing ratio, temperature and extinction-to-backscatter ratio of aerosol in the lower troposphere. The technique of Raman–Rayleigh lidar can retrieve correct temperature profile in the lower troposphere where the measurements are underestimated due to the aerosol loading. Two typical cases are discussed under different humidity (dry/wet) conditions. The water vapor and temperature profile have shown a good agreement with radiosonde. Simultaneous measurement of Raman–Rayleigh lidar also illustrates the physical nature of the aerosol and is useful in understanding the effects of humidity on aerosol swelling.  相似文献   

10.
荆思佳  肖薇  王伟  刘强  张圳  胡诚  李旭辉 《湖泊科学》2019,31(6):1698-1712
湖泊模型为数值天气预报模型提供热量通量、水汽通量和动量通量等下边界条件,但是不同时间尺度上湖泊水热通量变化的控制因子不同,因此有必要对湖泊模型进行多时间尺度上的离线评估.本文利用2012-2016年太湖中尺度通量网避风港站的气象资料和辐射数据驱动CLM4-LISSS模型(Community Land Model version 4-Lake,Ice,Snow and Sediment Simulator),并与涡度相关观测(Eddy Covariance,EC)结果进行对比,以年平均潜热通量模拟结果最佳为目标调整了模式中的消光系数、粗糙度长度方案,研究了该模型从半小时到年尺度上对湖表温度和水热通量的模拟性能.结果表明:模型对湖表温度的模拟在各时间尺度上均比较理想,但是模拟的日较差较小;从半小时到年尺度上潜热通量的变化趋势都能被很好地模拟出来,但在季节尺度上,潜热通量的模拟出现了秋冬季偏高、春夏季偏低的情况,季节变化模拟不准确.湖表温度和潜热通量模拟偏差的原因可能是消光系数的参数化方案.相比之下,感热通量尽管年际变化趋势的模拟值与观测值一致,但是从半小时到年尺度均被高估.特别地,冷锋过境期间,模型能较好地模拟出潜热通量和感热通量的变化趋势,但对于高风速条件下的感热通量模拟效果不佳.本文的研究结果能为湖泊模式的应用与发展提供有用信息.  相似文献   

11.
The snowcover energy balance is typically dominated by net radiation and sensible and latent heat fluxes. Validation of the two latter components is rare and often difficult to undertake at complex mountain sites. Latent heat flux, the focus of this paper, is the primary coupling mechanism between the snow surface and the atmosphere. It accounts for the critical exchange of mass (sublimation or condensation), along with the associated snowcover energy loss or gain. Measured and modelled latent heat fluxes at a wind‐exposed and wind‐sheltered site were compared to evaluate variability in model parameters. A well‐tested and well‐validated snowcover energy balance model, Snobal, was selected for this comparison because of previously successful applications of the model at these sites and because of the adjustability of the parameters specific to latent heat transfer within the model. Simulated latent heat flux and snow water equivalent (SWE) were not sensitive to different formulations of the stability profile functions associated with heat transfer calculations. The model parameters of snow surface roughness length and active snow layer thickness were used to improve latent heat flux simulations while retaining accuracy in the simulation of the SWE at an exposed and sheltered study site. Optimal parameters for simulated latent heat flux and SWE were found at the exposed site with a shorter roughness length and thicker active layer, and at the sheltered site with a longer roughness length and thinner active layer. These findings were linked to physical characteristics of the study sites and will allow for adoption into other snow models that use similar parameters. Physical characteristics of wind exposure and cover could also be used to distribute critical parameters in a spatially distributed modelling domain and aid in parameter selection for application to other watersheds where detailed information is not available. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

12.
This paper discusses the preliminary results of a study on the vegetation pattern and its relationship with meteorological parameters in and around Istanbul. The study covers an area of over 6800 km2 consisting of urban and suburban centers, and uses the visible and near-infrared bands of Landsat. The spatial variation of the Normalized Difference Vegetation Index (NDVI) and meteorological parameters such as sensible heat flux, momentum flux, relative humidity, moist static energy, rainfall rate and temperature have been investigated based on observations in ten stations in the European (Thracian) and Anatolian parts of Istanbul. NDVI values have been evaluated from the Landsat data for a single day, viz. 24 October 1986, using ERDAS in ten different classes. The simultaneous spatial variations of sensible heat and momentum fluxes have been computed from the wind and temperature profiles using the Monin-Obukhov similarity theory. The static energy variations are based on the surface meteorological observations. There is very good correlation between NDVI and rainfall rate. Good correlation also exists between: NDVI and relative humidity; NDVI, sensible heat flux and relative humidity; NDVI, momentum flux and emissivity; and NDVI, sensible heat flux and emissivity. The study suggests that the momentum flux has only marginal impact on NDVI. Due to rapid urbanization, the coastal belt is characterized by reduced NDVI compared to the interior areas, suggesting that thermodynamic discontinuities considerably influence the vegetation pattern. This study is useful for the investigation of small-scale circulation models, especially in urban and suburban areas where differential heating leads to the formation of heat islands. In the long run, such studies on a global scale are vital to gain accurate, timely information on the distribution of vegetation on the earth’s surface. This may lead to an understanding of how changes in land cover affect phenomena as diverse as the atmospheric CO2 concentrations, the hydrological cycle and the energy balance at the surface-atmosphere interface.  相似文献   

13.
Nimbus 7 LIMS geopotential height data are utilized to infer the rotational wind distribution in the Northern Hemisphere stratosphere and lower mesosphere during a period of substantial wave-mean flow interaction in January, 1979. Rotational winds are derived from the application of a successive relaxation numerical procedure which incorporates the spherical polar coordinate iterative algorithm ofPaegle andTomlinson (1975) for the nondivergent nonlinear balance equation. Optimum convergence of the numerical solutions is found to occur when under-relaxation is utilized. The LIMS height analyses were also latitudinally smoothed and constrained to obey the ellipticity criterion for spherical coordinates. The balanced winds are compared with geostrophically derived values and within situ radiosonde reports for 100 mb to 10 mb over Berlin.From a localized perspective, the Berlin-LIMS comparison indicates that radiosonde and balanced wind vectors exhibit somewhat closer agreement in direction than is associated with the geostrophic estimates. However, substantial quantitative differences between radiosonde, balanced, and geostrophic wind speeds are also evident, suggesting that caution should be exercised in the local application of derived winds, as for example in the quantitative interpretation of trajectories derived from satellite height analyses during periods of enhanced stratospheric wave activity.On a longitudinally averaged basis, balanced zonal-mean wind speeds are typically 20% weaker than geostrophic values in polar latitudes, and as much as 50% weaker in tropical and midlatitude regions. Meridional balanced wind velocities, at a given longitude, are generally within ±10% of geostrophic values. Although these alterations in horizontal wind components result in only modest differences between balanced and geostrophic meridional eddy heat fluxes, a more substantial change appears in the meridional eddy momentum flux analysis. The corresponding patterns of Eliassen-Palm flux divergence are found to be somewhat more (less) intense for the balanced wind case in the stratosphere (lower mesosphere) in polar latitudes.  相似文献   

14.
The current study deals with a parameterization of diapycnal diffusivity in an ocean model. The parameterization estimates the diapycnal diffusivity depending on the location of tidal-related energy dissipation over rough topography. The scheme requires a bottom roughness map that can be chosen depending on the scales of topographic features. Here, we implement the parameterization on an ocean general circulation model, and we examine the sensitivity of the modeled circulations to different spatial scales of the modeled bottom roughness. We compare three simulations that include the tidal mixing scheme using bottom roughness calculated at three different ranges of spatial scales, with the largest scale varying up to 200?km. Three main results are discussed. First, the dependence of the topographic spectra with depth, characterized by an increase in spectral energy over short length scales in the deep ocean, influences the vertical profile of the diffusivity. Second, the changes in diffusivities lead to different equilibrium solutions in the Atlantic meridional overturning circulation and bottom circulation. In particular, the lower cell of the Atlantic overturning and the bottom water transport in the Pacific Ocean are stronger for stronger diffusivities at the corresponding basins and depths, and the strongest when using the small-scale roughness map. Third, a comparison of the density fields of the three simulations with the density field of World Ocean Atlas dataset, from which the models are initialized, shows that among the simulations with three different roughness maps, the one using small-scale bottom roughness map has the smallest density bias.  相似文献   

15.
Influence of the surface roughness, geostrophic wind speed and initial stable stratification on surface layer scales, namely the friction velocity, temperature scale and the Monin-Obukhov length, and on non-dimensional profiles has been studied in quasi-stationary PIAPBLM (Prague Institute of Atmospheric Physics Boundary Layer Model) runs. Modifications of the stability functions, of the non-local mixing length concept, and of surface layer parametrization have improved the model results. With the ground temperature kept constant in time, the boundary layer is subjected to a turbulent cooling induced by surface roughness and to a counteracting turbulent heat transfer directed downwards. If the wind speed is lower, a rather mixed layer develops capped by a more stable layer. The estimated Monin-Obukhov length then slightly increases and the non-dimensional gradients overpredict the empirical values even more.  相似文献   

16.
Roughness length and zero-plane displacement over three typical surfaces were calculated iteratively by least-square method, which are Yucheng Experimental Station for agriculture surfaces, Qianyanzhou Experimental Station for complex and undulant surfaces, and Changbai Mountains Experimental Station for forest surfaces. On the basis of roughness length dynamic, the effects of roughness length dynamic on fluxes were analyzed with SEBS model. The results indicate that, aerodynamic roughness length changes with vegetation conditions (such as vegetation height, LAI), wind speed, friction velocity and some other factors. In Yucheng and Changbai Mountains Experimental Station, aerodynamic roughness length over the fetch of flux tower changes with vegetation height and LAI obviously, that is, with the increase of LAI, roughness length increases to the peak value firstly, and then decreases. In Qianyanzhou Experimental Station, LAI changes slightly, so the relationship between roughness length and LAI is not obvious. The aerodynamic roughness length of Yucheng and Changbai Mountains Experimental Station changes slightly with wind direction, while aerodynamic roughness length of Qianyanzhou Experimental Station changes obviously with wind direction. The reason for that is the terrain in Yucheng and Changbai Mountains Experimental Station is relatively flat, while in Qianyanzhou Experimental Station the terrain is very undulant and heterogeneous. With the increase of wind speed, aerodynamic roughness length of Yucheng Experimental Station changes slightly, while it decreases obviously in Qianyanzhou Experimental Station and Changbai Mountains Experimental Station. Roughness length dynamic takes great effects on fluxes calculation, and the effects are analyzed by SEBS model. By comparing 1 day averaged roughness length in Yucheng Experimental Station and 5 day averaged roughness length of Qianyanzhou and Changbai Mountains Experimental Station with roughness length parameter chosen by the model, the effects of roughness length dynamic on flux calculation is analyzed. The maximum effect of roughness length dynamic on sensible heat flux is 2.726%, 33.802% and 18.105%, in Yucheng, Qianyanzhou, and Changbai Mountains experimental stations, respectively.  相似文献   

17.
Zhou  Yanlian  Sun  Xiaomin  Zhu  Zhilin  Zhang  Renhua  Tian  Jing  Liu  Yunfen  Guan  Dexin  Yuan  Guofu 《中国科学:地球科学(英文版)》2006,49(2):262-272

Roughness length and zero-plane displacement over three typical surfaces were calculated iteratively by least-square method, which are Yucheng Experimental Station for agriculture surfaces, Qianyanzhou Experimental Station for complex and undulant surfaces, and Changbai Mountains Experimental Station for forest surfaces. On the basis of roughness length dynamic, the effects of roughness length dynamic on fluxes were analyzed with SEBS model. The results indicate that, aerodynamic roughness length changes with vegetation conditions (such as vegetation height, LAI), wind speed, friction velocity and some other factors. In Yucheng and Changbai Mountains Experimental Station, aerodynamic roughness length over the fetch of flux tower changes with vegetation height and LAI obviously, that is, with the increase of LAI, roughness length increases to the peak value firstly, and then decreases. In Qianyanzhou Experimental Station, LAI changes slightly, so the relationship between roughness length and LAI is not obvious. The aerodynamic roughness length of Yucheng and Changbai Mountains Experimental Station changes slightly with wind direction, while aerodynamic roughness length of Qianyanzhou Experimental Station changes obviously with wind direction. The reason for that is the terrain in Yucheng and Changbai Mountains Experimental Station is relatively flat, while in Qianyanzhou Experimental Station the terrain is very undulant and heterogeneous. With the increase of wind speed, aerodynamic roughness length of Yucheng Experimental Station changes slightly, while it decreases obviously in Qianyanzhou Experimental Station and Changbai Mountains Experimental Station. Roughness length dynamic takes great effects on fluxes calculation, and the effects are analyzed by SEBS model. By comparing 1 day averaged roughness length in Yucheng Experimental Station and 5 day averaged roughness length of Qianyanzhou and Changbai Mountains Experimental Station with roughness length parameter chosen by the model, the effects of roughness length dynamic on flux calculation is analyzed. The maximum effect of roughness length dynamic on sensible heat flux is 2.726%, 33.802% and 18.105%, in Yucheng, Qianyanzhou, and Changbai Mountains experimental stations, respectively.

  相似文献   

18.
利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL)湍流观测资料,分析了二次坐标旋转(double rotation,DR)、平面拟合(planar fit,PF)和分风区平面拟合(fetch planar fit,FPF)在复杂地形上的适用性,总结出一套适用于SACOL的总体湍流特征参数化方案.经过超声虚温订正、坐标旋转、空气密度脉动订正以及平稳性检验、总体湍流特征检验、总体质量分级处理,摩擦速度(u*)、感热通量、潜热通量、CO_2通量高质量数据所占比例分别为45%~62%、66%~68%、62%~65%、52%~54%.采用DR得到的高质量数据比例与采用PF相比,u*提高了17%,后三种通量略降低2%~3%.PF和FPF两种结果的差别主要体现在u*上,只考虑主导风向数据DR得到的u*质量仍最好.综合兼顾数据质量和计算工作量,在复杂地形上处理湍流观测资料的最优坐标旋转方法是DR.  相似文献   

19.

The current study deals with a parameterization of diapycnal diffusivity in an ocean model. The parameterization estimates the diapycnal diffusivity depending on the location of tidal-related energy dissipation over rough topography. The scheme requires a bottom roughness map that can be chosen depending on the scales of topographic features. Here, we implement the parameterization on an ocean general circulation model, and we examine the sensitivity of the modeled circulations to different spatial scales of the modeled bottom roughness. We compare three simulations that include the tidal mixing scheme using bottom roughness calculated at three different ranges of spatial scales, with the largest scale varying up to 200 km. Three main results are discussed. First, the dependence of the topographic spectra with depth, characterized by an increase in spectral energy over short length scales in the deep ocean, influences the vertical profile of the diffusivity. Second, the changes in diffusivities lead to different equilibrium solutions in the Atlantic meridional overturning circulation and bottom circulation. In particular, the lower cell of the Atlantic overturning and the bottom water transport in the Pacific Ocean are stronger for stronger diffusivities at the corresponding basins and depths, and the strongest when using the small-scale roughness map. Third, a comparison of the density fields of the three simulations with the density field of World Ocean Atlas dataset, from which the models are initialized, shows that among the simulations with three different roughness maps, the one using small-scale bottom roughness map has the smallest density bias.

  相似文献   

20.
Heat flow values were calculated from direct measurements of temperature and thermal conductivity at thirteen sites in the Arkansas-Missouri Ozark Plateau region. These thirteen values are augmented by 101 estimates of heat flow, based on thermal conductivity measurements and temperature gradients extrapolated from bottom-hole temperatures. The regional heat flow profile ranges from 9 mW m−2 to over 80 mW m−2, but at least two distinct thermal regimes have been identified. Seven new heat flow determinations are combined with three previously published values for the St. Francois Mountains (SFM), a Precambrian exposure of granitic and rhyolitic basement rocks, average 47 mW m−2. Radioactive heat production of 76 samples of the exposed rocks in the SFM averages 2.4 μW m−2 and a typical continental basement contribution of 14 mW m−2 is implied. Conversely, the sedimentary rock sequence of the plateau is characterized by an anomalously low heat flow, averaging approximately 27 mW m−2. Groundwater transmissivity values that are based on data from 153 wells in deep regional aquifers demonstrate an inverse relationship to the observed heat flow patterns. The areas of high transmissivity that correspond to areas of low total heat flux suggest that the non-conservative vertical heat flow within the Ozark sedimentary sequence can be attributed to the effects of groundwater flow.  相似文献   

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

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