首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
京津冀地区夏季暴雨频发,水汽是影响暴雨形成的关键要素之一.本文利用中国大陆构造环境监测网络(CMONOC)京津冀地区GNSS(全球导航卫星系统)观测资料,开展GNSS天顶对流层总延迟(ZTD)时序分析及对厄尔尼诺事件的响应研究.利用快速傅里叶变换与小波变换方法从频域和时域开展GNSS ZTD时序分析,并对GNSS ZTD不同周期时序与东部型指数(IEP)、中部型指数(ICP)进行比较,分析IEPICP对GNSS ZTD周期变化的影响.研究发现:GNSS ZTD异常时段与厄尔尼诺事件存在对应关系.东部型指数(IEP)与GNSS ZTD呈正相关;中部型指数(ICP)与GNSS ZTD呈显著负相关.在东部型厄尔尼诺事件的影响下,GNSS ZTD的季节性周期增大,月周期和半月周期减小;在中部型厄尔尼诺事件的影响下,GNSS ZTD的季节性周期、月周期、半月周期都减小.研究结果可为掌握区域GNSS ZTD预测变化规律提供参考,并为利用水汽感知厄尔尼诺事件提供可行性基础.  相似文献   

2.
采用水量平衡模型和Penman公式分别计算了珠江流域七个子流域1961—2000年实际蒸散发(I_(ETa))和潜在蒸散发(I_(ETp)),并对供水条件变化下I_(ETa)与I_(ETp)的关系进行了定量化分析,对各子流域I_(ETa)和I_(ETp)关系的理论从属性进行判定,主要结论如下:1)珠江流域年实际蒸散发量远低于潜在蒸散发量,多数子流域I_(ETa)值不到I_(ETp)值的1/2。7个流域面积加权平均I_(ETa)为681.4 mm/a,I_(ETp)为1 560.8 mm/a。从蒸散发的变异性来看,则实际蒸散发I_(ETa)的变异性明显要高于潜在蒸散发I_(ETp)。2)东江、西江、北江、柳江和盘江等5个流域实际蒸散发I_(ETa)都与降水量呈现正相关关系,韩江、郁江两个流域I_(ETa)随降水变化的变化趋势不明显。各子流域的潜在蒸散发I_(ETp)与降水量呈现显著负相关关系。7个子流域平均情况下,随着降水量的增加,I_(ETa)呈现明显的增加趋势,而I_(ETp)呈现明显的下降趋势。3)通过对降水量P与实际蒸散发I_(ETa)及潜在蒸散发I_(ETp)的联合回归方程P-IET回归系数的T检验,判定韩江、柳江和盘江等三个子流域以及七流域面积加权平均I_(ETa)与P和I_(ETp)与P的关系满足理论意义上的严格互补相关;东江、西江、北江等三个流域I_(ETa)与P和I_(ETp)与P的关系满足"非对称"互补相关。4)基于极端干旱和极端湿润的边界条件,推导出非对称条件下的实际蒸散发互补相关理论模型。  相似文献   

3.
Li  Yi  Yao  Ning  Sahin  Sinan  Appels  Willemijn M. 《Theoretical and Applied Climatology》2017,129(3-4):1017-1034

Global increases in duration and prevalence of droughts require detailed drought characterization at various spatial and temporal scales. In this study, drought severity in Xinjiang, China was investigated between 1961 and 2012. Using meteorological data from 55 weather stations, the UNEP (1993) index (I A), Erinç’s aridity index (I m), and Sahin’s aridity index (I sh) were calculated at the monthly and annual timescales and compared to the Penman-Monteith based standard precipitation evapotranspiration index (SPEIPM). Drought spatiotemporal variability was analyzed for north (NX), south (SX), and entire Xinjiang (EX). I m could not be calculated at 51 stations in winter as T max was below 0. At the monthly timescale, I A, I m, and I sh correlated poorly to SPEIPM because of seasonality and temporal variability, but annual I A, I m, and I sh correlated well with SPEIPM. Annual I A, I m, and I sh showed strong spatial variability. The 15 extreme droughts denoted by monthly SPEIPM occurred in NX but out of phase in SX. Annual precipitation, maximum temperature, and relative and specific humidity increased, while air pressure and potential evapotranspiration decreased over 1961–2012. The resulting increases in the four drought indices indicated that drought severity in Xinjiang decreased, because the local climate became warmer and wetter.

  相似文献   

4.
This work examines the spatial patterns of the transient response of mean annual temperature and precipitation to CO2 (or CO2 plus aerosol or aerosol proxy) radiative forcing in eight coupled AOGCMs, generally for the period 1900–2099. Response patterns are characterized using empirical orthogonal functions (EOFs) and the quasi-EOFs of Harvey and Wigley (the first qEOF field, discussed here, is given by the correlation between local year-by-year temperature changes and the global mean temperature change). The first temperature EOF accounts for 80–95% of the space-time variation of the CO2 run in all of the models, and is almost identical to qEOF1 of the temperature response or to the temperature change pattern averaged over the last 30 years of the simulations. EOF1 accounts for 80–95% of the space-time variation in the CO2+aerosol runs in six of the eight models. The CO2 response patterns of different models are highly correlated with one another (R 2 generally >0.5), and are also highly correlated with the CO2+aerosol response patterns (R 2 0.85 in all except one model). The difference between CO2 and CO2+aerosol runs can be represented by EOF1 of the year-by-year differences, by qEOF1 of the year-by-year differences, or by the difference in temperature averaged over the last 30 years of each run. In models where these representations are highly correlated with each other, they are also highly correlated with CO2 EOF1. In other cases, aerosol EOF1 is modestly to highly correlated with control EOF1 (i.e.: the year-by-year differences between CO2 and CO2+aerosol runs are dominated by internal variability), while aerosol qEOF1 and the 30-year difference are highly correlated with each other. For all models, the decadal mean temperature change can be closely replicated by scaling the CO2 EOF1 pattern based on the global mean temperature changes (RMSE for the last decade is <6% of the RMS temperature change for CO2 runs, <8% for CO2+aerosol runs). The first EOF of the precipitation response to increasing CO2 accounts for only 10–30% of the space-time variation, and is generally highly correlated (R 2 up to 0.85) with control EOF1. In all of the models, there is an increase in precipitation in the ITCZ and a decrease in bands at or near 30°S and 30°N. In many models there is an El Niño-like response, including a substantial decrease in precipitation over the Amazon. Global-mean precipitation increases in all models due to CO2 forcing, but aerosols appear to have a disproportionally large effect in suppressing the increase compared to their effect in suppressing the warming. There is evidence in some models that the non-absorbing aerosols considered here reduce summer monsoon rainfall compared to the changes that would be expected based on the globally averaged effect of aerosols on precipitation. When regional precipitation changes over time are predicted by scaling a fixed precipitation-change pattern with the global mean temperature change, the global mean RMSE in the predicted change in decadal-mean precipitation is 25–35% of the global RMS precipitation changes by the end of the simulation.  相似文献   

5.
The model of Paillard and Parrenin (Earth Planet Sci Lett 227:263–271, 2004) was modified to obtain a closer fit to δ18O and CO2 time series for the last 800 kyr. The model performance can be improved if its CO2 sensitivity to I65 insolation is eliminated and if different response times are assumed for ablation/accumulation of ice. Correlations between simulated and experimental time series for CO2 and ice volume V increase from 0.59 and 0.63 to 0.79 and 0.88, respectively. According to these models, terminations are produced by I65 amplification through CO2-T and T-CO2 feedbacks, in synergy with an extra CO2 contribution from the deep ocean. This contribution is strongly dependent on ice-sheet extent and ice volume (or alternatively, CO2 concentration, which is a good proxy of Antarctic temperature) but is insensitive to Southern Ocean (SO) insolation on 21 February (I60). Change of deep SO state may be the “order parameter” for nonlinear deglacial changes. According to these models, 100 kyr periodicity of glacial cycles arises from the characteristic time of Antarctic ice sheet advance to the continental slope.  相似文献   

6.
An analysis was performed of the turbulent data obtained from Yucheng experimental station in the Shandong Province in 1984. [t is shown that at variant wind speed, the spectra of streamwise velocity remain similar and the intensity of wind fluctuations is proportional to wind speed in the downwind area of shelter belt. Therefore, we may decide the similarity of wind fluctuations by a speed scale and a length scale which is not correlated with stability, σu /V0 = F(X / H). The -5/3 power range of temperature spectra extends to lower frequency. The variation of ratio σ0 /T. with stability becomes σ0 / T . = C(X / H)( - Z / L)-1/3 . There is not such an extension of -5 / 3 power range in the humidity spectra.  相似文献   

7.
近百年四季西太平洋副热带高压的变化   总被引:3,自引:1,他引:3  
文中设计了一种根据网格点高度值拟合副热带高压强度、西界、北界 3种指数的方法。根据这个方法 ,利用重建的 1880~ 195 0年 5 0 0hPa月平均高度序列 ,计算了 195 1年之前 70a的副热带高压指数 ,与 195 1~ 1999年的实测指数合成一个 12 0a的序列 ,研究了西太平洋副热带高压的年际与年代际变率。  相似文献   

8.
We present several equilibrium runs under varying atmospheric CO2 concentrations using the University of Victoria Earth System Climate Model (UVic ESCM). The model shows two very different responses: for CO2 concentrations of 400 ppm or lower, the system evolves into an equilibrium state. For CO2 concentrations of 440 ppm or higher, the system starts oscillating between a state with vigorous deep water formation in the Southern Ocean and a state with no deep water formation in the Southern Ocean. The flushing events result in a rapid increase in atmospheric temperatures, degassing of CO2 and therefore an increase in atmospheric CO2 concentrations, and a reduction of sea ice cover in the Southern Ocean. They also cool the deep ocean worldwide. After the flush, the deep ocean warms slowly again and CO2 is taken up by the ocean until the stratification becomes unstable again at high latitudes thousands of years later. The existence of a threshold in CO2 concentration which places the UVic ESCM in either an oscillating or non-oscillating state makes our results intriguing. If the UVic ESCM captures a mechanism that is present and important in the real climate system, the consequences would comprise a rapid increase in atmospheric carbon dioxide concentrations of several tens of ppm, an increase in global surface temperature of the order of 1–2°C, local temperature changes of the order of 6°C and a profound change in ocean stratification, deep water temperature and sea ice cover.  相似文献   

9.
冬季东亚大槽强度年际变化及其与中国气候联系的再分析   总被引:2,自引:1,他引:1  
利用NCEP/NCAR再分析资料以及国家气象局整编的中国160站逐月降水和气温资料,定义了一个北半球冬季500hPa东亚大槽强度指数(IEAT),并分析了该指数所反映的冬季东亚大槽强度的年际变化规律及其与同期中国冬季降水、气温的关系.结果表明:IEAT指数反映了对流层中层冬季亚洲大陆高压脊与西北太平洋上空的低压槽系统之间的平均经向风强度,且能够较好地反映冬季东亚大槽的强度.东亚大槽的强度在20世纪80年代之前相对较弱,并存在2-3年和准4年的年际变化周期.冬季东亚大槽强弱变化可能与源自地中海地区的罗斯贝波动能量沿亚洲急流东传有关,且这种西风带中的扰动具有准正压结构.IEAT指数与青藏高原东部地区的冬季降水和气温相关显着.当IEAT指数为正(负),东亚大槽偏强(弱),对应着中国华中地区以及华东大部分地区冬季总降水量偏少(多),且华中地区冬季平均气温偏高(低).进一步研究发现,在东亚大槽偏强年,华中地区冬季平均气温异常升高主要是由于异常非绝热加热和下沉运动导致的异常动力增温所致.这些研究结果有助于更好地理解由于东亚大槽强度的变化而导致的中国冬季气候变化特征及其原因.  相似文献   

10.
Summary A novel multi-timescale analysis method, Empirical Mode Decomposition (EMD), is used to diagnose the variation of the annual mean temperature data of the global, Northern Hemisphere (NH) and China from 1881 to 2002. The results show that: (1) Temperature can be completely decomposed into four timescales quasi-periodic oscillations including an ENSO-like mode, a 6–8-year signal, a 20-year signal and a 60-year signal, as well as a trend. With each contributing ration of the quasi-periodicity discussed, the trend and the 60-year timescale oscillation of temperature variation are the most prominent. (2) It has been noticed that whether on century-scale or 60-year scales, the global temperature tends to descend in the coming 20 years. (3) On quasi 60-year timescale, temperature abrupt changes in China precede those in the global and NH, which provides a denotation for global climate changes. Signs also show a drop in temperature in China on century scale in the next 20 years. (4) The dominant contribution of CO2 concentration to global temperature variation is the trend. However, its influence weight on global temperature variation accounts for no more than 40.19%, smaller than those of the natural climate changes on the rest four timescales. Despite the increasing trend in atmospheric CO2 concentration, the patterns of 20-year and 60-year oscillation of global temperature are all in falling. Therefore, if CO2 concentration remains constant at present, the CO2 greenhouse effect will be deficient in counterchecking the natural cooling of global climate in the following 20 years. Even though the CO2 greenhouse effect on global climate change is unsuspicious, it could have been excessively exaggerated. It is high time to re-consider the trend of global climate changes.  相似文献   

11.
A new circulation index (ISW) that can realistically describe winter temperature variations over Southwest China is defined based on analysis of the NCEP/NCAR reanalysis data (version 1) and the observations at 585 stations in China. The study period is from January 1961 to February 2011. The relationship between ISW and general circulation patterns in East Asia is also analyzed. Results show that ISW successfully captures the variations in winter temperature over Southwest China. High ISW values correspond to the intensified Mongolian high, the weakened Aleutian low, increases in the strength of the Middle East westerly jet stream over the south of the Tibetan Plateau (TP), and decreases in the strength of the subtropical westerly jet over the north of the TP. Meanwhile, the East Asian trough deepens and extends southwestward, making it easier for the cold air mass from the north to intrude Southwest China along the trough. These circulation patterns lead to a decrease in winter temperature over Southwest China (and vice versa). In addition to the East Asian winter monsoon, the two westerly jets that dominate the upper level circulation over East Asia also exert important influences on winter temperature in Southwest China, especially the Middle East westerly jet to the south of the TP.  相似文献   

12.
Abstract

The meteorological visual range Vm and the perception visual range Vp have been measured for a variety of horizon sky backgrounds and atmospheric extinction conditions. Comparison with the standard visibility Vn shows that in general Vp > Vm > Vn . Model calculations of Vp, Vm and Vn for different synthetic sky backgrounds and for an extinction coefficient function of the form σ(λ) = k λ‐α in which a was varied from ‐1 to 4 show that the wavelength of maximum perception can vary between ~540 and ~630 nm. Under typically turbid or polluted atmospheric conditions and for a blue or white sky background the extinction coefficient should be measured at ‐570 nm if prevailing visibility is to be reliably estimated.  相似文献   

13.
This work presents the results of a comparison between the global storm activity rate IRS and electric field intensity E0Z. The permanent analysis of the IRS may become an important tool for testing Global Electric Circuit models. IRS is determined by a new method that uses the background component of the first 7 Schumann resonances (SR). The rate calculations are based on ELF observations carried out in 2005 and 2006 in the observatory station “Hylaty” of the Jagiellonian University in the Eastern Carpathians (Kułak, A., Zięba, S., Micek, S., Nieckarz, Z., 2003. Solar variations in extremely low frequency propagation parameters: I. A two-dimensional telegraph equation (TDTE) model of ELF propagation and fundamental parameters of Schumann resonances, J. Geophys. Res., 108, 1270, doi:10.1029/2002JA009304). Diurnal runs of the IRS rate were compared with diurnal runs of E0Z amplitudes registered at the Earth's surface in the Geophysical Observatory of the Polish Academy of Sciences in Świder (Kubicki, M., 2005. Results of Atmospheric Electricity and Meteorological Observations, S. Kalinowski Geophysical Observatory at Świder 2004, Pub. Inst. Geophysics Polish Academy of Sciences, D-68 (383), Warszawa.). The days with the highest values of the correlation coefficient (R) between amplitudes of both observed parameters characterizing atmosphere electric activity are shown. The seasonal changes of R, IRS and E0Z are also presented.  相似文献   

14.
使用北京人工影响天气办公室提供的2014-2017年京津冀地区飞行记录积冰个例样本与机载观测数据,2016年全国空中报告积冰、非积冰个例样本和欧洲中期天气预报中心(ECMWF)第5代全球气候大气再分析数据(ERA5),基于模糊逻辑隶属度函数,定义了以气温和相对湿度为判别基础并考虑垂直速度和云量影响的积冰指数Ip(icing potential index),用于判断飞机在空中发生积冰事件的可能性。检验结果表明:该指数对积冰事件的判别准确率为80.2%,与目前国内常用的经典积冰指数(Ic)相比,其判别准确率有明显提升,且漏报率和虚警率均显著降低(分别为9.4%和10.4%),结合数值预报产品可对飞机在空中特定位置发生积冰事件的可能性进行预测。  相似文献   

15.
If a spot of tracer is released into a turbulent flow, the peak concentration at some subsequent time will initially be much greater than that implied by a solution for the ensemble average concentration at fixed points. For two-dimensional turbulence three areas may be defined: (1) an area Ad related to the ensemble average concentration field; (2) an area Ap defined in terms of the relative dispersion of particles seeded into the patch after a short initial diffusion time; and (3) the area At occupied by tracer. It is argued that Ad grows linearly with time, whereas Ap and At grow exponentially; Ap faster than At. Thus, the concentration field is significantly streaky, even within the particle domain, until At becomes comparable with Ad. The time taken for this to occur is estimated; after this time, fluctuations about the ensemble average concentration field should not be greater than those given by a simple mixing length argument. In three-dimensional turbulence the volume Vt of the tracer domain grows much more rapidly than the volume Vp of the particle domain if the merging of streaks is ignored. However, Vt cannot be greater than Vp so streaks must merge and Vp can be used to provide a rough estimate of peak concentration, or concentration variance.  相似文献   

16.
A set of semi-continuous measurements of temperature, wind and moisture gradients as well as of net radiation and ground heat flux covering a period of about one and a half years has been analysed to give a corresponding set of complete surface energy balance data on an hourly basis. An analysis of the evaporation data so obtained is given.It is shown that surface resistance r S exhibits a diurnal trend: values are smallest (ca. 150 s m-1) a few hours before noon and increase to as much as 800 s m-1 towards dusk. The minimum values tend to be higher during dry periods when the soil moisture is low. There is also some indication that r S decreases rapidly soon after rainfall.An exponential relation is found between the fraction of available energy used as evaporative flux, , and r S for values of r I/rS <0.70, where r I is the climatological resistance. On the other hand, the ratio of r I to r S is linearly correlated with , implying that an equilibrium state is established between the grass surface and the atmosphere, at least from mid-morning to mid-afternoon when the leaves are dry. Near-noon values calculated by Stewart and Thom for Thetford Forest also follow a linear trend.The above two regression results (In (r S) versus r I/rS versus ) are combined to obtain an empirical relation of the form r I=m exp (a-b) which is used to estimate evaporative flux. The estimates are found to be within 20% of calculated values.Guest Scientist from Department of Physics, University of Cape Coast, Cape Coast, Ghana.  相似文献   

17.
一次高原低涡与高原切变线演变过程与机理分析   总被引:6,自引:1,他引:5  
李山山  李国平 《大气科学》2017,41(4):713-726
对一次东移高原低涡减弱、高原切变线生成并在有利的环流背景下东南移,进而引发川渝强降水的高原切变线生成机制以及演变过程进行了初步分析。首先引入描写热带气旋的Okubo-Weiss(OW)参数(VOW)来定量表达低涡、切变气流中旋转和变形的相对大小,确定高原切变线的潜在生成区域和发展状况。得出在高原切变线生成阶段,500 hPa等压面上VOW值由正转负,VOW负值带可以很好地指示高原切变线的潜在生成区域。VOW负值强度与高原切变线强度有很好的相关性。高原切变线上以VOW负值中心为主,但也会存在正值中心,说明在切变线上也会有气旋性涡度。此个例高原切变线以伸缩变形为主,高原切变线沿变形场的拉伸轴分布。然后通过涡度方程和总变形方程进一步分析了高原低涡减弱、高原切变线生成的动力机制。高原低涡的减弱、消失主要受散度项的影响,时间演变分析表明系统由强气旋性涡度的高原低涡演变为强辐合性的高原切变线。总变形方程中的扭转项对高原切变线的生成贡献最大,其次为水平气压梯度项,散度项最小;当高原切变线上以拉伸变形为主时,不利于其上高原低涡的发展,切变线可能是影响低涡发展的背景流场。  相似文献   

18.
 The diurnal range of surface air temperature (rT a ) simulated for present and doubled CO2 climates by the CSIRO9 GCM is analysed. Based on mean diurnal cycles of temperature and surface heat fluxes, a theory for understanding the results is developed. The cycles are described as the response to a diurnal forcing which is represented well by the diurnal mean flux of net shortwave radiation at the surface (SW) minus the evaporative (E) and sensible (H) fluxes. The response is modified by heat absorbed by the ground, and by the cycle in downward longwave (LW) radiation, but these effects are nearly proportional to the range in surface temperature. Thus in seasonal means, rT a is approximately given by SWEH divided by 6 W m-2/°C. A multiple regression model for (rT a ) is developed, based on quantities known to influence SW, E and H, and applied to both spatial variation in seasonal means, and day-to-day variation at a range of locations. In both cases, rT a is shown to be influenced by cloud cover, snow extent and wind speed. It is influenced by soil moisture, although this effect is closely tied to that of cloud. In seasonal means rT a is also well correlated with precipitable water, apparently because of the latter’s influence on E+H. The regression model describes well the spatial variation in the doubled CO2 change in rT a . The annual mean change in rT a over land on doubling CO2 was −0.36 °C, partly because of a decrease in the mean diurnal forcing (as defined in the theory), but also apparently because of the effect of nonlinearity in T s of the upward longwave emission. A diagnostic radiation calculation indicates that the CO2 and water vapour provide a small increase in rT a through the downward LW response, which partially counters a decrease due to a reduction of SW by the gases. Received: 8 November 1995 / Accepted: 3 January 1997  相似文献   

19.
青藏高原感热指数的建立及与华南降水的联系   总被引:1,自引:0,他引:1       下载免费PDF全文
利用1982-2012年青藏高原中东部70个气象站的月平均地面感热资料、华南地区92个气象站的月平均降水资料、NCEP/NCAR月平均再分析资料和SEOF(season reliant EOF)方法选取了4个高原代表站,建立了青藏高原地面感热强度距平指数(ISH),并讨论了春季ISH与华南盛夏(7月和8月)降水的关系。结果表明:ISH可以较好地表征青藏高原中东部地面感热的年际变化特征,且具有更好的持续性。春季ISH与华南盛夏降水具有显著的负相关关系,当春季ISH偏大时,后期对流层中上层高度场异常偏高,且高度场异常偏高的响应随时间从低层向高层传递,使夏季副热带高压偏强、偏西,南亚高压异常偏强,华南地区盛夏降水偏少;反之亦然。此外,去除Ni?o3.4区海温对华南盛夏降水影响后,两者的负相关关系变得更为显著。  相似文献   

20.
湿地是由陆地和水体形成的自然综合体,具有重要的生态、水文和生物地球化学功能,黄河源高寒湿地作为黄河重要的水源涵养区,对其下垫面水热交换特征及关键影响参数的研究具有非常重要的意义。本文利用中国科学院西北生态环境资源研究院麻多黄河源气候与环境变化观测站2014年6~8月观测资料,分析了黄河源区高寒湿地-大气间暖季水热交换特征,并利用公用陆面模式(Community Land Model,简称CLM)模拟了热通量变化,提出针对高寒湿地的粗糙度优化方案。主要结果如下:(1)暖季向上、向下短波与净辐射的平均日变化规律一致,向上、向下长波平均日变化平缓,地表温度升高相对于向下短波具有滞后性,潜热通量始终为正值并大于感热通量;(2)温度变化显著层结为20 cm以上土壤浅层,存在明显的日循环规律,土壤中热量09:00(北京时,下同)下传至5 cm深度,温度升高,11:00至10 cm深度,13:00至20 cm深度,18:00后开始上传,温度降低,40 cm及以下深度受此影响较小,热量在土壤中整体由浅层向深层输送;(3)土壤湿度平均日变化小,5 cm深度为土壤湿度最小层,10 cm深度为最大层;(4)麻多高寒湿地动力学粗糙度Z0m在暖季变化稳定,可作为常数,Z0m=0.0143 m;(5)提出更加适合高寒湿地下垫面暖季附加阻尼kB-1参数化方案,使得热通量模拟效果较CLM原始方案有所提高。以上结果对于研究湿地下垫面陆面过程具有重要意义。  相似文献   

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

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