首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Summary The aim of this study was to investigate possible effects of two hypothetical scenarios of the urbanization of Zagreb’s surroundings on the local winds, which are established under summertime anticyclonic conditions. For this purpose, the nonhydrostatic mesoscale meteorological model MEMO was applied to the greater Zagreb area. Three simulations were performed. One employed the current land-use distribution, while the other two corresponded to an increase of the densely urbanized area by 12.5% (test 1) and 37.5% (test 2), respectively. Apart from the hypothetically urbanized areas, where average surface wind speed reductions of 8% and 18% were obtained for test 1 and test 2, respectively, the rest of the domain was not significantly affected by hypothetical urbanization. The differences between the wind vectors for the predicted current state and the hypothetical state were more pronounced and found at higher altitudes during the night compared to daytime values. For all three simulations the same diurnal variation of the depth of anabatic/katabatic wind flow generated on south-facing slopes of 1 km high mountain Medvednica was obtained. During the night the depth of well-developed katabatic flow was about 370 m, while during the day the depth of anabatic flow grew from about 550 m in the late morning up to about 1140 m in the late afternoon. Received October 27, 2000 Revised August 4, 2001  相似文献   

2.
Summary Sobradinho Lake lies in the S?o Francisco River Basin, in one of the most arid regions in Northeastern Brazil, within a land stretch categorized as the Lower-middle S?o Francisco, situated at about 40 km away from the municipality of Petrolina (09°23′S–40°30′W) in the state of Pernambuco. The dam, in its full capacity, consists of a lake of approximately 4,214 km2; 280 km in length, the width of which varies from 5 to 50 km. The dam storage capacity is that of 34.1 billion m3 of water. Being situated in a semi-arid region, the dam brought about significant development to local irrigated agriculture. The caatinga ecosystem has, for that matter, undergone considerable changes. Statistical analysis techniques applied to data collected before and after the filling of the lake, made it possible both to make an assessment of the impact of the dam construction on the region meteorology and to diagnose the variability of such an impact on environmental conditions. Results showed that the dam has brought about considerable changes to regional meteorology. The alterations were observed to be more significant as regards atmospheric humidity and wind speed.  相似文献   

3.
Mesoscale aspects of the Urban Heat Island around New York City   总被引:15,自引:1,他引:15  
Summary ?A mesoscale analysis of the Urban Heat Island (UHI) of New York City (NYC) is performed using a mesoscale network of weather stations. In all seasons the UHI switches on rapidly in late afternoon and shuts down even more rapidly shortly after dawn. It averages about 4 °C in summer and autumn and 3 °C in winter and spring. It is largest on nights with clear skies, low relative humidity through much of the troposphere, and weak northwest winds, when it may exceed 8 °C. The synoptic meteorological situation associated with the largest UHI occurs roughly two to three nights after cold front passages. During spring and summer, sea breezes commonly reduce and delay the UHI and displace it about 10 km to the west. Backdoor cold fronts, which occur most frequently in spring and early summer, reduce or even reverse the UHI, as cold air from the water to the northeast keeps NYC colder than the western suburbs. Cases documenting the sensitivity and rapidity of changes of the UHI to changes in parameters such as cloud cover, ceiling, and wind speed and direction are presented. Received August 16, 2001; revised October 6, 2002; accepted November 20, 2002 Published online March 17, 2003  相似文献   

4.
Summary The boundary-layer structure of the Elqui Valley is investigated, which is situated in the arid north of Chile and extends from the Pacific Ocean in the west to the Andes in the east. The climate is dominated by the south-eastern Pacific subtropical anticyclone and the cold Humboldt Current. This combination leads to considerable temperature and moisture gradients between the coast and the valley and results in the evolution of sea and valley wind systems. The contribution of these mesoscale wind systems to the heat and moisture budget of the valley atmosphere is estimated, based on radiosoundings performed near the coast and in the valley. Near the coast, a well-mixed cloud-topped boundary layer exists. Both, the temperature and the specific humidity do not change considerably during the day. In the stratus layer the potential temperature increases, while the specific humidity decreases slightly with height. The top of the thin stratus layer, about 300 m in depth, is marked by an inversion. Moderate sea breeze winds of 3–4 m s−1 prevail in the sub-cloud and cloud layer during daytime, but no land breeze develops during the night. The nocturnal valley atmosphere is characterized by a strong and 900 m deep stably stratified boundary layer. During the day, no pronounced well-mixed layer with a capping inversion develops in the valley. Above a super-adiabatic surface layer of about 150 m depth, a stably stratified layer prevails throughout the day. However, heating can be observed within a layer above the surface 800 m deep. Heat and moisture budget estimations show that sensible heat flux convergence exceeds cold air advection in the morning, while both processes compensate each other around noon, such that the temperature evolution stagnates. In the afternoon, cold air advection predominates and leads to net cooling of the boundary layer. Furthermore, the advection of moist air results in the accumulation of moisture during the noon and afternoon period, while latent heat flux convergence is of minor relevance to the moisture budget of the boundary layer. Correspondence: Norbert Kalthoff, Institut für Meteorologie und Klimaforschung, Universit?t Karlsruhe/Forschungszentrum Karlsruhe, Postfach 3640, 76021 Karlsruhe, Germany  相似文献   

5.
Summary  Mesosphere-Stratosphere-Troposphere (MST) Radar wind data for the period June through September 1996 have been examined to study vertical variation of Madden-Jullian Oscillations in wind and eddy kinetic energy (eke) in the normal monsoon season. The domain of analysis in the vertical is from 6 to 20 km with a height resolution of 150 m. Fast-Fourier-Transformation (FFT) has been applied to zonal (u), meridional(v) components of wind to extract the Madden-Jullian oscillations and eke. There are three dominant modes viz., 50–70, 30–40 and 10–20 day periodicity, which contain considerable fraction of energy and show high degree of vertical variability. The peak amplitude of 50–70 day mode in u, 30–40 mode in v and eke were observed at 16–17 km just below the tropopause level. The peak amplitudes of 30–40 day mode in u and 50–70 day mode in v were found in the height region of 13–16 km. To understand the origin and propagation of these waves, wave energy is calculated. The wave energy is higher at tropospheric heights than at lower stratospheric heights indicating that the origin of these waves is in the troposphere, and a part of the energy leaks into the stratosphere. Received September 17, 1998/Revised September 26, 1999  相似文献   

6.
Summary Data from two automatic stations in Łódź (one urban and one rural) for the period 1997–2002 are analyzed to reveal urban–rural contrasts of such parameters as air temperature, relative humidity, water vapour pressure and wind speed. Under favourable weather conditions the highest temperature differences between the urban and rural station exceeds 8 °C. Relative humidity is lower in the town, sometimes by more than 40%. Water vapour pressure differences can be either positive (up to 5 hPa) or negative (up to −4 hPa). Wind speed at the urban station is on average lower by about 34% in night and 39% during daytime. Regression analysis shows that for rural winds lower than 1.13 m s−1 urban winds can be stronger than rural speeds. Attention has also been paid to singularities in the course of the analyzed parameters over 24 hour periods. It is shown that the typical course of the urban heat island intensity under favourable conditions is similar in all season. Four stages of this course have been distinguished. Wind speed differences also seem to change in a typical way. Case studies show that humidity contrasts, unlike temperature, can evolve in different ways under fine weather conditions. Types of relative humidity evolution are proposed.  相似文献   

7.
From 1973–1976, research was performed around the Sea of Galilee, aimed at examining the wind regime in the area and whether the area develops a land-sea breeze despite its particular topographical location.
    The main conclusions were:
  1. During the summer mornings a lake breeze develops, blowing towards the shores of the lake. It ceases at the peak of its development when a westerly wind, originating in the development of a breeze along the Israeli Mediterranean coast, plunges towards the lake.
  2. Late at night, a wind flow develops from the land towards the lake, which combines with the katabatic winds that blow along the steep slopes surrounding the Kinneret.
  3. The stations at the upper level, at a height of 400–500 m above the Kinneret, are not affected by the lake breeze during the day or by the land breeze at night.
  4. In winter, the Kinneret lake breeze is almost as developed as in summer, because the westerly winds, originating in the Mediterranean sea breeze which hardly develops in this season, do not plunge into the Kinneret.
  相似文献   

8.
利用东疆红柳河黑戈壁下垫面陆气相互作用观测站2017年近地大气边界层梯度探测资料和红柳河气象站天气现象观测数据,分析该地区典型晴天条件下的近地层风速、温度和比湿的四季廓线特征。结果表明:四季近地层风廓线变化规律明显。典型晴天条件下,在0.5~4 m高度内风速随高度的增加而变大的速度较快,在4~32 m范围内,白天风速随高度增大较缓慢,但夜间出现快速增大;存在明显的夜间逆温,逆温层主要集中在4~32 m,冬季逆温强于夏季,晨间0.5~32 m间的温度差可达4.6℃,且红柳河四季的气温日较差均较大,秋季可达到15.7℃;夜间比湿高于白天,秋、冬季夜间逆湿层出现在10~32 m,其比湿差为0.15 g/kg左右,夏季无逆湿现象。  相似文献   

9.
大气细颗粒物PM2.5污染引起的雾霾天气既与本地污染物排放密切有关,也受局地特殊的风场影响.本文以武汉城市区域为研究对象,分别研究了长江沿岸的江陆风环流、东湖沿岸湖陆风环流的形成与转化特征,发现江风、湖风开始时间均为07:00-08:00,一年中最大风速均能达到2 m/s左右,而春夏季江风的持续时间高于秋冬季,夏季湖风的持续时间高于春季.同时发现区域附近温度和相对湿度之间有明显的相关性,湿度变化趋势大体上跟温度变化趋势相反,且温度对相对湿度的影响存在一定的滞后性,延迟时间大约为1 h.  相似文献   

10.
Summary The local wind system in the upper Isar Valley (Bavarian Alps) near Mittenwald has the peculiarity that regularly strong foehn-like nocturnal flows occur, mainly during clear nights in autumn and winter. We will refer to this phenomenon as “Minifoehn”, as its properties are similar to the classical deep foehn in the sense that its breakthrough into the Isar Valley usually brings a striking increase in temperature and a concomitant decrease in relative humidity. Numerical simulations with the MM5 model reveal that this phenomenon is related to a nocturnal drainage flow originating from a plateau south of Mittenwald. This flow is driven by the temperature difference between this plateau (1180 m) and the free atmosphere above Mittenwald (920 m, 15 km north of the plateau) at the same level. The air masses flow through two different valleys that merge again further downstream. The upper part of one of the two drainage currents goes over a small mountain ridge (1180 m) south-west of Mittenwald and then descends into the Isar Valley, leading to an advection of potentially warm air towards Mittenwald. This branch of the drainage current constitutes the Minifoehn. The remaining part of the drainage current flows through a narrow gap towards the Isar Valley and then joins the drainage flow of this valley. As these air masses are significantly cooler than the Minifoehn branch, large horizontal temperature gradients can be found around Mittenwald. The dynamical behaviour of the cold air flow turns out to be qualitatively consistent with shallow-water theory only in the absence of a forcing by large-scale winds. Otherwise, gravity-wave induced pressure perturbations interact with the drainage flow and modify the low-level flow field. The simulations show that the gravity waves are excited by the mountain range that separates the two valleys mentioned above. Moreover, the simulations indicate that the structure of this nocturnal wind system is not very sensitive to the direction of synoptic-scale winds as long as they come from the southern sector. On the other hand, ambient northerly winds are able to prevent the drainage flow and therefore the local foehn effects in the Isar Valley provided that synoptic winds are strong enough. The results of the MM5 simulations are in good agreement with the measurements and observations described in part 1 of this study.  相似文献   

11.
Summary A newly developed non-hydrostatic model (MOLOCH), operating at a resolution of about 2 km, is run for a case of heavy precipitation over southeastern France. The event (8–9 September 2002) was characterized by intense convective activity leading to a severe flash flood in the region of the Gard river, south of the Massif Central. An almost stationary mesoscale convective system (MCS), developing well in advance of an approaching cold front, discharged a huge amount of rainfall over the same area, more than 600 mm in 24 hours. Several simulations are performed in order to test the model set-up, evaluate the sensitivity on different initial conditions, and analyse the case-study. The quantitative precipitation forecasts (QPF) appear to vary widely among the experiments, depending on the initialization time chosen (00, 06, and 12 UTC, September 8). Only the run starting at 06 UTC predicts, with a satisfactory degree of accuracy, the location where the MCS developed and its almost stationary behavior during the first stage (∼12 hours) of the event. In all the simulations, the convective system then propagated northward over the Massif Central. In addition to experiments starting from standard ECMWF analyses, an assimilation procedure, based on Optimal Interpolation, is applied to the initial conditions. Surface observations of temperature, wind and relative humidity have been assimilated. The assimilation produces an improvement in the forecasts of surface fields and leads to a better location of the initial triggering phase. Further experiments, performed by changing the orography in the model, allow the investigation of the role of the Massif Central in triggering the mesoscale convective system and in controlling its evolution.  相似文献   

12.
Summary  High resolution aircraft observations made along flight tracks over inhomogeneous surface in the late wintertime boreal zone are described and compared to 2D mesoscale model simulations with surface properties defined at 2 km resolution from maps. All observations displayed the expected small-scale turbulence. On top of that, the near-surface wind speeds (but not directions) showed mesoscale variations related to local topography and roughness. Upward (but not downward) SW and LW radiative fluxes and ground temperature also displayed mesoscale variability; in SW radiation this was clearly due to local albedo changes. In the sensible heat flux there was strong horizontal variation near the surface in correlation with surface types. The above observed mesoscale along-track variations were reasonably well represented by the mesoscale model simulation. The track-averaged observed sensible and latent heat flux profiles were in rough agreement with a mixing length approach, which used the track-averaged wind, temperature and moisture profiles as input (mimicking a first-order turbulence closure scheme of a GCM). Received September 20, 1999 Revised January 21, 2000  相似文献   

13.
Summary  The winter wind regime of G?teborg, located on the West coast of Sweden, is composed of three different wind systems besides the ambient wind; a nocturnal low level jet (NLLJ), a winter land breeze (WLB) and an urban heat island circulation (UHIC). An inversion divides the air column into two layers, one between 10 – 50 m and one between 50 – 100 m. The UHIC is located in the lower layer, the WLB in the top layer and the NLLJ above the top layer. The intensity of the interacting processes depends on the stability of each layer as calculated from the bulk Richardson number (BRilow and BRihigh) using continuous data collected during four years (1991 – 94) from two sites (one within and one outside the urban area) and sampled at three levels. In the evening the WLB develops from the ground level and increases in height until after midnight. At about the same time an UHIC develops in the urban area, below the WLB and causing an uplift of the latter. However, at both sites the WLB does not exceed the 100 m level. At this time BRi in both layers are below one resulting in continuous coupling between the WLB, the UHIC layers and the regional wind. Consequently, the exchange of momentum is still effective between all layers and this is highlighted by a change in the wind direction and a regulation of wind-speed to more constant levels. When BRihigh≥1, the layers become frictionally decoupled, as indicated by a return in the wind direction in the top level to the regional wind, and an acceleration of the top wind. The top level then becomes incorporated in to a nocturnal low-level jet (NLLJ) system. The normally acknowledged development of the NLLJ, with a start around sunset, is in this case delayed for several hours at the top level. The reason for this is that there are meso-scale/local wind systems present in layers beneath the jet causing an interaction between the layers. In the morning, when the layers are again coupled the top layer wind is once more influenced by the WLB and therefore changes direction and speed. The local and meso-scale wind systems thus delay the current nocturnal wind development. Received August 24, 1998 Revised March 17, 1999  相似文献   

14.
Summary The stability parameter μ is suggested as the one which is determinable with satisfying accuracy for routine application by means of commonly accessible meteorological data at the Molve location (Croatia). The similarity functions applied for vertical wind speed simulation in the planetary boundary layer (PBL) at Molve were useful for the determination of local stability classes. Universal similarity functions were applied for unstable and neutral stability, whereas local similarity functions were established for stable stratification. Wind speed simulations were performed using two types of wind models. The Monin-Obukhov similarity theory was included in both types. However, it turned out that for the operative determination of the stability of the 35 m deep lowest layer, the stability parameter μ was locally a better stability parameter than the Monin-Obukhov parameter z/L. That was possibly because 35 m deep lowest layer sometimes (depending upon stability) includes a large proportion of the Ekman layer and parameter μ is originally designed for the deeper part of PBL than z/L that is originally designed for the surface layer. At Molve, the input data for local wind models as well as for the stability parameter μ were wind speed at 35 m and temperature at 2 and 35 m above the ground.  相似文献   

15.
Summary ?Simultaneous flight measurements with the research aircraft Do 128 and the helicopter-borne turbulence probe Helipod were performed on 18 June 1998 during the LITFASS-98 field experiment. The area-averaged turbulent vertical fluxes of momentum, sensible, and latent heat were determined on a 15 km × 15 km and a 10 km × 10 km flight pattern, respectively. The flights were carried out over heterogeneous terrain at different altitudes within a moderately convective boundary layer with Cumulus clouds. Co-spectra-analysis demonstrated that the small scale turbulent transport was completely sampled, while the comparatively small flight patterns were possibly of critical size regarding the large-scale turbulence. The phygoide of the airplane was identified as a significant peak in some co-spectra. The turbulent fluxes of momentum and sensible heat at 80 m above the ground showed systematic dependence on the location of the flight legs above the heterogeneous terrain. This was not observed for the latent heat flux, probably due to the vertical distribution of humidity in the boundary layer. Statistical error analysis of the fluxes F showed that the systematic statistical error ΔF was one order of magnitude smaller than the standard deviation σ F . The difference between area-averaged fluxes derived from simultaneous Helipod and Do 128 measurements was much smaller than σ F , indicating that the systematic statistical error was possibly over-estimated by the usual method. In the upper half of the boundary layer the airborne-measured sensible heat flux agreed well with windprofiler/RASS data. A linear fit was the best approximation for the height dependence of all three fluxes. The linear extrapolations of the latent and sensible heat fluxes to the ground were in good agreement with tower, scintillometer, and averaged ground-station measurements on various surface types. Systematic discrepancies between airborne and ground-based measurements were not found. Received June 18, 2001; revised December 21, 2001; accepted June 3, 2002  相似文献   

16.
Postglacial fire history has been reconstructed for eastern Canada from charcoal-influx anomalies from 30 sites taken from a lacustrine charcoal database. The reconstruction exhibits coherent patterns of fire occurrence in space and time. The early Holocene is characterised by high fire incidence. There is a major change to much lower occurrence slightly after 8 ka BP. A return to more fire appears after 3 ka BP. This sequence does not fit with the hydro-climatic reconstruction deduced from lake level reconstructions for northeastern North America, which indicates a dry early and mid-Holocene, and a wet late-Holocene. Fire occurrence however closely matches summer relative humidity inferred from δ18O. The differences between fire frequency and lake level history, are due to changes in the seasonality of precipitation and drought frequency. Lake levels are essentially controlled by winter precipitation while summer precipitation controls fire occurrence. The early Holocene before 8–7.5 ka BP experienced dry summers due to higher solar radiation and dry adiabatic winds from the residual Laurentide Ice Sheet. The middle Holocene was dominated by wet summers due to stability of the Atlantic air mass over eastern Canada. After 2.5 ka BP, summers became drier, albeit not as fire-conducive as during the early Holocene. Late-Holocene summers conducive to fire are explained by more frequent incursions of dry Cool Pacific or Cold Arctic air masses over eastern Canada. Received: 25 January 1999 / Accepted: 14 December 1999  相似文献   

17.
This study conducted meteorological simulations in northern Colombia by analyzing different planetary boundary layer (PBL) schemes available in the numerical Weather Research and Forecasting (WRF) model. The study area included three nested domains with horizontal resolutions of 18 km, 6 km, and 2 km, with 38 vertical levels. The evolution and structure of the PBL were analyzed during the driest months (March, April, and May 2016) and in regions with the highest particulate matter concentrations. Sensitivity analysis of the WRF model was performed with two local and two non-local PBL schemes. The results were validated using observations of the surface air temperature, relative humidity, and surface wind speed collected from three meteorological stations in the area. The PBL heights were experimentally determined using radiosonde data provided by a station located in the center of the study area. Variations in PBL heights were estimated using linear regression methods and minimization of statistical errors for the bulk Richardson number, as well as analysis of vertical temperature and wind profiles. The WRF model reliably reproduced the daily values and diurnal cycles of temperature, relative humidity, and wind speed within the PBL and accounted for the influence of topography and sea breezes. Horizontal heat advection dominates the upwelling of air masses when sea breezes are active. The onshore wind direction starts to change from east to northwest, implying a decay in the land breeze regime. All schemes overestimate the mixing height and tend to underestimate surface air temperature values at night. All show wetter conditions and underestimate wind speed. Although the non-local Yonsei University (YSU) scheme shows the best performance, it also shows the largest sources of errors when determining the behavior of the surface layer during stable conditions. Relative humidity and wind speed estimates provided by the local Mellor‐Yamada‐Nakanishi‐Niino (MYNN) scheme were closer to those recorded at the meteorological stations.  相似文献   

18.
塔中春季阴天近地层风速、温度和湿度廓线特征分析   总被引:1,自引:0,他引:1  
利用塔中最新安装的80 m梯度观测塔探测系统采集的资料,详细分析了2006年4月2日1次阴天天气时塔克拉玛干沙漠腹地近地层风速、温度和湿度廓线演变特征,并与典型晴天廓线做了对比,得到以下结果:(1)阴天,夜间风速廓线风速值随高度增高而增大,但不是以对数增长,而是以比对数关系更快的速度增长;白天,风速很小,近地层10 m上下廓线分布规律各异;(2)温度廓线有夜间辐射型、早上过渡型、白天日射型及傍晚过渡型4种类型,与晴天类似;(3)比湿廓线存在一个极小值,其出现高度以上比湿随高度增加而增加,廓线呈逆湿特征,极小值出现高度以下比湿随高度减小而增加。  相似文献   

19.
Summary The application of principal components and cluster analysis to vertical ozone concentration profiles in Tsukuba, Japan, has been explored. Average monthly profiles and profiles of the ratio between standard deviation and the absolute ozone concentration (SDPR) of 1 km data were calculated from the original ozone concentration data. Mean (first) and gradient (second) components explained more than 80% of the variation in both the 0–6 km tropospheric and 11–20 km troposphere–stratosphere (interspheric) layers. The principal components analysis not only reproduced the expected inverse relationship between mean ozone concentration and tropopause height (r 2 = 0.41) and that in the tropospheric layer this is larger in spring and summer, but also yielded new information as follows. The larger gradient component score in summer for the interspheric layer points to the seasonal variation of the troposphere–stratosphere exchange. The minimum SDPR was at about 3 km in the tropospheric layer and the maximum was at about 17 km in the interspheric layer. The tropospheric SDPR mean component score was larger in summer, possibly reflecting the mixing of Pacific maritime air masses with urban air masses. The cluster analysis of the monthly ozone profiles for the 1970s and 2000s revealed different patterns for winter and summer. The month of May was part of the winter pattern in the 1970s but part of the summer pattern during the 2000s. This statistically detected change likely reflects the influence of global warming. Thus, these two statistical analysis techniques can be powerful tools for identifying features of ozone concentration profiles. Authors’ addresses: S. Yonemura, S. Kawashima, S. Inoue, National Institute for Agro-Environmental Sciences, 3-1-3 Kannondai, Tsukuba, Ibaraki 305-0031, Japan; H. Matsueda, Y. Sawa, Meteorological Research Institute, 1-1 Nagamine, Tsukuba, Ibaraki 305-0052, Japan; H. Tanimoto, National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki 305-8506, Japan.  相似文献   

20.
湿地小气候效应特征研究   总被引:1,自引:0,他引:1  
湿地对局地小气候具有调节作用,研究湿地小气候效应特征能更具体地了解湿地对局地小气候的影响。本文以河北省衡水市的衡水湖为例,利用衡水市11个常规气象观测站数据,通过对湖区及湖区外各季节不同气象要素的对比,对衡水湖各个季节的小气候效应进行了分析。结果表明:(1)衡水湖具有冷岛效应、湿岛效应和风岛效应,能够调节周围的气候特征;(2)衡水湖的小气候效应具有季节特征,衡水湖各季节平均的冷岛效应由强到弱依次为春季、冬季、秋季、夏季,湿岛效应由强到弱分别为夏季、春季、秋季、冬季,风岛效应由强到弱依次为春季、夏季、冬季、秋季,春季小气候效应最强;(3)衡水湖的小气候效应具有昼夜特征,夜晚的冷岛效应强于白天,湿岛和风岛效应正相反,白天的强度大于夜晚。  相似文献   

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

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