首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
This paper describes implementation of the warm-rain bin microphysics in a LES model based on the EULAG fluid flow solver. The binmicrophysics EULAG is applied to the case of shallow nonprecipitating tropical convection to investigate the impact of the secondary activation of cloud droplets above the cloud base. In a previous study applying the EULAG model with the double-moment bulk warm-rain microphysics scheme, the in-cloud activation was shown to have significant implications for the mean microphysical and optical characteristics of the cloud field. By contrasting the simulations with and without in-cloud activation as in the previous study, we show that the in-cloud activation has qualitatively similar but quantitatively smaller effect. In particular, the concentration of cloud droplets in the bin simulation without in-cloud activation decreases with height not as strongly as in corresponding simulations applying the double-moment bulk scheme.  相似文献   

2.
气溶胶对雷暴云电过程影响的数值模拟研究   总被引:3,自引:0,他引:3       下载免费PDF全文
本文在已有的三维雷暴云起、放电模式中加入了一种经典的气溶胶活化参数化方案,结合一次长春雷暴个例,进行了雷暴云起放电数值模拟试验.研究显示气溶胶浓度改变对雷暴云微物理、起电及放电过程都有重要影响.结果表明:(1)污染型雷暴云中气溶胶浓度增加时,云滴数目增多,上升风速加强;云中冰晶与霰粒子数浓度增加但尺度减小;(2)相对于清洁型雷暴云,污染型雷暴云非感应起电过程弱,感应起电过程强,起电持续时间长;(3)污染型雷暴云中首次放电时间延迟,闪电持续发生的时间长,总闪电频次增加,正地闪频次增加明显.  相似文献   

3.
During a study of the growth of cloud drops by condensation the evolution of cloud drop size spectra with height above cloud base was determined for maritime aerosols, and for continental aerosols containing aerosol particles of mixed composition. Air parcel models were used in which the parcel was either completely closed to mass and heat transfer (strictly adiabatic models), or open to heat transfer and to partial or complete mass transfer (entrainment models). It was found that adiabatic models and models which consider the entraining of air devoid of aerosol particles predict drop size distributions which are considerably narrower than those observed in non-precipitating cumulus clouds, and have only a single maximum. On the other hand, relative broad drop size distributions and distributions with a double maximum — as they are observed in atmospheric clouds — are predicted if the entrainment of both air and aerosol particles are considered in the condensation model. Our results support the findings ofWarner (1973) which were obtained for a purely maritime aerosol.  相似文献   

4.
UV-B solar irradiance and meteorological variables were measured at the C.I.B.A. site (Low Atmosphere Research Laboratory), University of Valladolid, Spain, between January 2003 and March 2006. Calculated Ångström turbidity and aerosol optical thickness values were evaluated from the direct downward irradiance, surface pressure, air temperature and relative humidity values. Monthly turbidity β values for an average year showed minimum values in winter and maximum values in summer. The obtained values are according to an aerosol standard atmosphere between mean and clean continental model. UV-B model calculations were performed using a radiative transfer tropospheric model, TUV 4.1a; measured and calculated UV-B irradiance values were compared in order to establish possible values of single scattering albedo. The results show that calculated single scattering albedo values are according to an aerosol standard atmosphere between mean and polluted continental model.  相似文献   

5.
本文利用大气环流模式及大气化学模式所得气溶胶资料,估算了相对1850s时期硫酸盐和黑碳气溶胶引起的全球及东亚区域人为辐射强迫,重点分析其在东亚区域的季节和长期变化特征.结果表明,就当前全球年平均全天空而言,人为硫酸盐气溶胶对大气顶的直接和云反照率强迫分别为-0.37和-0.98 W·m-2,黑碳气溶胶对大气顶和整层大气的辐射强迫值为0.16和0.47 W·m-2;中国东部区域是目前上述气溶胶辐射强迫最强的区域,硫酸盐的直接和间接辐射强迫分别超过-2.0和-4.0 W·m-2,黑碳对大气顶和整层大气的直接辐射强迫分别可达2.0和5.0W·m-2;估算的东亚区域上述气溶胶辐射强迫仍在不断增强,峰值预计出现在2010s时段,而且中国东部较强的辐射强迫还可能维持至2030s左右;在未来中、高排放情景下,东亚区域以上两种气溶胶预计对全球气溶胶辐射强迫有更大的贡献.分析还表明,夏季东亚区域较强的水汽会增强吸湿性硫酸盐气溶胶的光学厚度和晴空直接辐射强迫;云的作用一方面会强化东亚区域全天空条件下大气顶黑碳的辐射强迫,另一方面会影响硫酸盐气溶胶间接云反照率强迫的季节变化;上述气候特征的差异使得东亚区域的气溶胶辐射强迫表现出与欧美区域有所不同的特征.本文所用的气溶胶资料与模式气象场的偏差会给气溶胶辐射强迫计算带来一些不确定性,进一步改进气候模式中气溶胶过程、水汽和云等气象场的模拟将有助于获得更为合理的区域气溶胶辐射强迫估算结果.  相似文献   

6.
A detailed study using environmental tracers such as chloride (Cl?) and tritium (3H), deuterium (2H) and oxygen (18O) isotopes was performed in an alluvial coastal aquifer in two contrasting environments (urban and agricultural). These environmental tracers combined with a high‐resolution multi‐level sampling approach were used to estimate groundwater residence time and recharge patterns and to validate the hydrogeochemical conceptual model already proposed in previous studies. δ18O and δ2H combined with Cl? data proved that the hypersaline groundwater present in the deepest part of the aquifer was sourced from the underlying hypersaline aquitard via an upward flux. Both chemical and isotopic data were employed to calibrate a density‐dependent numerical model based on SEAWAT 4.0, where 3H and Cl? were helped quantifying solutes transport within the modelled aquifer. Model results highlighted the differences on estimated recharge in the two contrasting environments, with the urban one exhibiting concentrated recharge because of preferential infiltration associated to the storm water drains network, while scarce local recharge characterized the agriculture setting. In the urban field site, is still possible to recognize at 9 m b.g.l. the input of the atmospheric anthropogenic 3H generated by testing of thermonuclear weapons, while in the agricultural field site, the 3H peak has been washed out at 6 m b.g.l. because the groundwater circulation is restricted only to the upper fresh part of the aquifer, drained by the reclamation system. The presented approach that combined high‐resolution field monitoring, environmental tracers and numerical modelling, resulted effective in validating the conceptual model of the aquifer salinization. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

7.
We examine the impact of the Atmospheric Brown Clouds on the direct radiative forcing of the Himalayan foothills and the Indo-Gangetic Plains (IGP) regions, home for over 500 million S. Asians. The NASA-Terra MODIS satellite data reveal an extensive layer of aerosols covering the entire IGP and Himalayan foothills region with seasonal mean AODs of about 0.4 to 0.5 in the visible wavelengths (0.55 micron), which fall among the largest seasonal mean dry season AODs for the tropics. We show new surface data which reveal the presence of strongly absorbing aerosols that lead to a large reduction in solar radiation fluxes at the surface during the October to May period. The three-year mean (2001 to 2003) October to May seasonal and diurnal average reduction in surface solar radiation for the IGP region is about 32 (±5) W m−2 (about 10% of TOA insolation or 20% of surface insolation). The forcing efficiency (forcing per unit optical depth) is as large as −27% (note that the forcing is negative) of top-of-atmosphere (TOA) solar insolation, and exceeds the forcing efficiency that has been observed for other polluted regions in America, Africa, East Asia, and Europe. General circulation model sensitivity studies suggest that both the local and remote influence of the aerosol induced radiative forcing is to strengthen the lower atmosphere inversion, stabilize the boundary layer, amplify the climatological tendency for a drier troposphere, and decrease evaporation. These aerosol-induced changes could potentially increase the life times of aerosols, make them more persistent, and decrease their single scattering albedos, thus potentially leading to a detrimental positive feedback between aerosol concentrations, aerosol forcing, and aerosol persistence. In addition, both the model studies and observations of pan evaporation suggest that the reduction in surface solar radiation may have led to a reduction in surface evaporation of moisture. These results suggest the vulnerability of this vital region to air pollution related direct and indirect (through climate changes) impacts on agricultural productivity of the region.  相似文献   

8.
X. Qie  S. Soula  S. Chauzy 《Annales Geophysicae》1994,12(12):1218-1228
A numerical model called PICASSO [Production d’Ions Corona Au Sol Sous Orage (French) and Production of Corona Ions at the Ground Beneath Thundercloud (English)], previously designed, is used to describe the evolution of the principal electrical parameters below a thunderstorm, taking into account the major part played by corona ions. In order to improve the model restitution of a real situation, various improvements are performed: an initial vertical distribution of aerosol particles is introduced instead of the previously used uniform concentration; time and space calculation steps are adjusted according to the electric field variation rate; the upper boundary condition is improved; and the coefficients of ion attachment are reconsidered with an exhaustive bibliographic study. The influence of the ion attachment on aerosol particles, on the electric field and charge density aloft, is studied by using three different initial aerosol particle concentrations at ground level and two types of initial vertical distributions: uniform and non-uniform. The comparison between field data and model results leads to adjust the initial aerosol particle concentration over the experimental site at the value of 5000 cm−3 which appears to be highly realistic. The evolutions of the electric field and of the charge density at altitude are greatly influenced by the aerosol concentration. On the contrary, the surface intrinsic field, defined as the electric field that would exist underneath a thundercloud if there were no local charges, is weakly affected when the model is forced by the surface field. A good correlation appears between the success in the triggered lightning attempts and this intrinsic field evaluation. Therefore, when only the surface field is available, the model can be used in a triggered lightning experiment.  相似文献   

9.
This study contains a comparative analysis of aerosol optical thickness (AOT) between numerical calculations obtained from the Navy Aerosol Analysis and Prediction System (NAAPS) model and direct observations from the AERONET robotic network and the Saharan Aerosol over WArsaw (SAWA) field campaign. AOT was calculated for 500 nm wavelength. The comparison shows underestimation of the total aerosol optical thickness simulated by NAAPS. The correlation coefficients between model and observation oscillates between 0.57 and 0.72. Results of seven-year (1998–2004) NAAPS simulation of aerosol components (sea salt, mineral dust, sulphate, and smoke) show large temporal and spatial variability of the aerosol optical thickness over Europe. The least polluted region is the Iberian Peninsula, while the highest aerosol burdens occurred in Central Europe, mostly due to anthropogenic sulphate particles. Finally, the analysis of mineral dust transport shows frequent episodes of Saharan dust inflow over Central Europe. There are about 20 days a year (4 days in May) when instantaneous AOT associated with mineral dust aerosol increases over 0.1.  相似文献   

10.
A simplified one dimensional analytical model to describe the altitude variation of aerosol number densities in the atmospheric boundary layer during clear sky and calm/moderate wind conditions, over a plane terrain on the tropical seacoast, is developed. The analytical model is mainly defined by two parameters, the aerosol mixing height and the exponential rate of decrease above the mixing height. The effect of surface impaction is also taken into account in this model. It is seen that this parameter varies with season. By matching the analytical profile with the observed mean profile in each month, the surface source strength has been estimated. The seasonal and inter-annual variations of this parameter are studied.  相似文献   

11.
Highly seasonal boreal catchments are hydrologically complex and generally data poor and, hence, are ripe for investigation using tracer‐aided hydrologic models. The influence of physiography on isotopic metrics was assessed to identify the catchment characteristics dominating evaporative enrichment. A multiyear stable isotope of water dataset was collected at the outlets of 16 boreal catchments in central Canada ranging in area from 12 to 15,282 km2. Physiographic characteristics were obtained through raster analysis of freely available land cover images, stream networks, and digital elevation models. Correlation analysis indicated that as the percentage coverage of open water increased, so too did the evaporative effects observed at the catchment outlet. Correlation to wetland metrics indicated that increasing the percentage coverage of wetlands can reduce or increase evaporative effects observed, depending on the isotopic metric used and the corresponding drainage density, catchment slope, and presence of headwater lakes. The slopes of river evaporative‐mixing lines appear to reflect multifaceted relationships, strongest between catchment slope, headwater lakes, and connected wetlands, whereas mean line‐conditioned excess is more directly linked to physiographic variables. Hence, the slopes of river evaporative‐mixing lines and mean line‐conditioned excess are not interchangeable metrics of evaporative enrichment in a catchment. Relationships identified appear to be independent of catchment scale. These results suggest that adequate inclusion of the distribution of open water throughout a catchment, adequate representation of wetland processes, catchment slope, and drainage density are critical characteristics to include in tracer‐aided hydrologic models in boreal environments in order to minimize structural uncertainty.  相似文献   

12.
CAMS云微物理方案的改进及与WRF模式耦合的个例研究   总被引:1,自引:0,他引:1       下载免费PDF全文
本文在中国气象科学研究院(CAMS)双参数云微物理方案的基础上,增加气溶胶粒子的活化过程,改进原方案中的水汽混合比、云水混合比及云滴数浓度的预报方程,实现对各种水成物(包括云水)的混合比和数浓度的预报.此外,改进后的CAMS云方案被成功耦合到了WRF v3.1中尺度模式.本文利用耦合模式对2009年4月23~24日发生在我国北方地区的一次降水天气过程进行了模拟,将新方案的模拟结果与WRF自带的3个微物理方案进行了比较.结果显示,新方案能够合理地描述地面降水特征,其模拟的雨带分布范围与实测接近,降水中心的强度和位置优于其他3个方案.新方案模拟的云滴数浓度与WDM6方案基本一致,表明加入的气溶胶活化过程是合理的.新方案模拟的其他水成物粒子数浓度与Morrison方案相比有时会有量级的差别,说明粒子数浓度的模拟目前还存在着很大的不确定性,这也是云微物理模式进一步发展的难点.  相似文献   

13.
A one-dimensional, time-dependent numerical cloud model is used to analyze the factors in the dynamic and thermodynamic equations which lead to a steady-state or nonsteady-state solution for the cloud vertical motion, buoyancy, precipitation, and cloud water fields. ‘Bulk water’ microphysical techniques are used for the cloud, rain, and hail variables. An atmospheric sounding from a severe storm situation is used as initial and environmental conditions, yielding model updrafts of 40 m sec?1 maximum and more than 10 m sec?1 over the entire cloud region. ‘Early conversion’ of the cloud water to rain leads to loading of lower portions of the updraft by rain, the formation of appreciable amounts of hail by freezing of the supercooled rain, and subsequent loading of the middle and upper portions of the updraft so that the updraft erodes throughout the cloud depth and the cloud dissipates, yielding a vigorous rain shower. A delay in the conversion of the cloud water to rain results in a steady-state solution, no rain or hail falling through the updraft. A two-dimensional cloud simulation of this same case shows rain and hail in the upper cloud regions recycled in the two-dimensional flow into the updraft near cloud base and a breakdown of the updraft with resultant rainout (negligible hail reaching the ground). The breakdown of the updraft has profound effects on the temperature field within the cloud, causing the lapse rate to deviate from the steady-state condition and approach the initial environmental conditions. The results emphasize the fact that the local change in temperature (and other dependent variables as well) is not independent of the vertical velocity, in general. This has implications for the interpretation of measurements made within clouds.  相似文献   

14.
Water transpired by trees has long been assumed to be sourced from the same subsurface water stocks that contribute to groundwater recharge and streamflow. However, recent investigations using dual water stable isotopes have shown an apparent ecohydrological separation between tree‐transpired water and stream water. Here we present evidence for such ecohydrological separation in two tropical environments in Puerto Rico where precipitation seasonality is relatively low and where precipitation is positively correlated with primary productivity. We determined the stable isotope signature of xylem water of 30 mahogany (Swietenia spp.) trees sampled during two periods with contrasting moisture status. Our results suggest that the separation between transpiration water and groundwater recharge/streamflow water might be related less to the temporal phasing of hydrologic inputs and primary productivity, and more to the fundamental processes that drive evaporative isotopic enrichment of residual soil water within the soil matrix. The lack of an evaporative signature of both groundwater and streams in the study area suggests that these water balance components have a water source that is transported quickly to deeper subsurface storage compared to waters that trees use. A Bayesian mixing model used to partition source water proportions of xylem water showed that groundwater contribution was greater for valley‐bottom, riparian trees than for ridge‐top trees. Groundwater contribution was also greater at the xeric site than at the mesic–hydric site. These model results (1) underline the utility of a simple linear mixing model, implemented in a Bayesian inference framework, in quantifying source water contributions at sites with contrasting physiographic characteristics, and (2) highlight the informed judgement that should be made in interpreting mixing model results, of import particularly in surveying groundwater use patterns by vegetation from regional to global scales. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
The isotope hydrology of a set of nested sub-catchments in the north-east of Scotland has been studied to examine the mixing processes and residence times of water in the catchments. The measured δ18O in stream waters was found to be exceptionally uniform both temporally and spatially. Hydrochemical mixing analyses showed that groundwater contributes between 62 and 90% of the stream flow in all sub-catchments. Model analysis indicated that the δ18O in stream water is indicative of a highly mixed system in which near surface runoff appears to be mixed with groundwater, within the soil profile, before being released from the catchment. Small fluctuations in the stream water δ18O response are generated by a small proportion (<10%) of less-well mixed water in infiltration excess runoff during storm events. A comparative application of the model to a nearby catchment, which has a lower proportion of groundwater runoff, demonstrated contrasting behaviour, with significantly less mixing of waters occurring and a more distinct difference in the age of runoff generated by different flow paths. This highlighted that standard methods for characterization of mixing mechanisms are often insufficient and may not discriminate between systems that have retained quite distinct flow paths throughout catchment transit, and those which have been mixed at some stage. Model sensitivity analysis also indicated that the simulated mean residence time of water varies most strongly in response to different parameters compared with the δ18O response. This has implications for estimating water residence times from isotope data. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

16.
The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.  相似文献   

17.
VIIRS(Visible Infrared Imaging Radiometer Suite)作为MODIS(The MODerate resolution Imaging Spectroradiometer)的后继传感器,可在全球范围内实现对气溶胶的连续时空监测.卫星反演的气溶胶光学厚度(Aerosol Optical Depth,AOD)是研究地球能量收支平衡、气候效应和空气质量的重要大气参数.但在中国重污染天气情况下,现有的VIIRS陆地气溶胶产品存在一定不足.因此,本研究改进云识别方法,优化像元筛选,约束气溶胶类型选择,实现重污染情况下AOD的反演.基于地基AERONET(AErosol RObotic NETwork)的验证结果表明,相比NOAA(National Oceanic and Atmospheric Administration)产品,改进后的反演结果克服了反演值偏低的问题,且表现出更好的相关性,RMSE从0.236下降到0.219.为验证在重污染条件下改进算法的适用性和准确性,本文对比了两种污染条件下的反演结果(0.61).统计结果表明,在较重污染天气条件下(AODAERONET>1),相比NOAA的AOD产品,本文结果的反演率从32.3%提升为68.8%,回归分析的斜率提高为0.80,相关系数达到0.76,均方根误差为0.307,在增加反演量的同时保证了反演的精度.  相似文献   

18.
It has been speculated for many years that the development of the droplet spectra in cloud is probably influenced by mixing processes. Various theoretical attempts to broaden the droplet spectra by mixing parcels with different velocity histories has shown that that particular effect is small. Similarly, very simpleuniform entrainment procedures did not lead to cloud drop size spectra which were broad enough, although by producing cloud drop size distributions with a double mode these models did substantially improve the drop size spectra of earlier adiabatic models which only exhibited a single mode.Recently a model based on entraining entities representing moving parcels of cloud air within the cloud was detailed byTelford andChai (1980). This study showed that the mixing in of dry air at cumulus turrets could lead to vertical cycling of diluted parcels, and that this cycling, with continual entrainment across the parcel boundaries, will produce much larger drops, as well as smaller drops of all sizes, in the droplet spectra. The entity entrainment concept studied there appears to apply to the observations of stratus cloud discussed in this paper.This paper presents data taken in marine stratus off the California coast which give a particularly clear example of how such droplet spectra modification occurs in practice. Both large drops, and the spread of the spectra to smaller sizes, occur in relation to other variables in such a way as to be consistent with an entity entrainment explanation, with no other obvious possibility.In a marine stratus cloud just over 200 m thick and many tens of miles in extent we find clear evidence that dry air is mixing in at cloud tops. Strong vertical motion is to be found in the cloud, large sized drops are found in cloud parcels where the mixing gives lower droplet concentrations, and there is evidence that newly formed cloud parcels are warmer and contain many more smaller droplets.The observations show that immediately following entrainment of dry air drop diameters are not reduced appreciably, but, in the same parcels, drop concentrations have been reduced by a factor of ten or more. Further down in the cloud big drops, able to start growth by coalescence, are found associated with low total droplet concentrations.Overall, it seems likely from the consideration of these observations that the formation of the large drops which lead to precipitation processes in clouds depends critically on the mixing in of dry air at cloud tops, and very little on the size of the small drops resulting from the condensation nucleus counts. As a conclusion it appears reasonable to state that if entrainment occurs at cloud tops, then big drops will be formed!  相似文献   

19.
A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

20.
We report observations of a noctilucent cloud (NLC) over central Alaska by a ground-based lidar and camera on the night of 9–10 August 2005. The lidar at Poker Flat Research Range (PFRR), Chatanika (65°N, 147°W) measured a maximum integrated backscatter coefficient of 2.4×10?6 sr?1 with a peak backscatter coefficient of 2.6×10?9 m?1 sr?1 corresponding to an aerosol backscatter ratio of 120 at an altitude of 82.1 km. The camera at Donnelly Dome, 168 km southeast of PFRR, recorded an extensive NLC display across the sky with distinct filamentary features corresponding to wave structures measured by the lidar. The occurrence of the maximum integrated backscatter coefficient corresponded to the passage of a bright cloud band to the southwest over PFRR. The camera observations indicate that the cloud band had a horizontal width of 50 km and a length of 150 km. The horizontal scale of the cloud band was confirmed by medium-frequency radar wind measurements that reported mesopause region winds of 30 m/s to the southwest during the period when the cloud band passed over PFRR. Comparison of these measurements with current NLC microphysical models suggests a lower bound on the water vapor mixing ratio at 83 km of 7–9 ppmv and a cloud ice mass of 1.5–1.8×103 kg. Satellite measurements show that this NLC display occurred during a burst of cloud activity that began on 5 August and lasted for 10 days. This cloud appeared 10 days after a launch of the space shuttle. We discuss the appearance of NLCs in August over several years at this lower polar latitude site in terms of planetary wave activity and space shuttle launches.  相似文献   

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

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