首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 656 毫秒
1.
The upward transfer of heat from ocean to atmosphere is examined for an Arctic lead, a break in the Arctic ice which allows contact between the cold atmosphere and the relatively warm ocean. We employ a large-eddy model to compute explicitly the three-dimensional turbulent response of the atmosphere to a lead of 200 m width. The surface heat flux creates a turbulent plume of individual quasi-random eddies, not a continuous updraft, which penetrate into the stable atmosphere and transport heat upward.Maximum updraft velocities and turbulence occur downwind of the lead rather than over the lead itself, because the development time of an individual thermal eddy is longer than its transit time across the lead. The affected vertical region, while shallow over the lead itself, grows to a height of 65m at 600 m downwind of the lead; beyond that, the depth of the turbulent region decreases as the eddies weaken. The maximum vertical turbulent heat flux occurs at the downwind edge of the lead, beyond which a relative maximum extends upward into the plume. Negative surface heat flux immediately downwind of the lead creates a growing stable layer, but above that internal boundary layer the turbulent heat flux is still positive. Updraft maxima are typically 28 cm/s, but compensating downdrafts result in time-averaged vertical velocities of less than 1 cm/s in the plume. Conditional sampling separates the updraft and downdraft contributions. Formulas for the horizontal eddy development distance and for the vertical plume penetration height are presented. The relative importance of mean and turbulent transport is compared for both vertical and horizontal heat transfer: turbulence dominates the vertical heat transport whereas mean advection dominates the horizontal transport, these offsetting transports producing a quasi-stationary state.  相似文献   

2.
Weather and climate networks traditionally follow rigorous siting guidelines, with individual stations located away from frost hollows, trees or urban areas. However, the diverse nature of the UK landscape suggests that the feasibility of siting stations that are truly representative of regional climate and free from distorting local effects is increasingly difficult. Whilst the urban heat island is a well-studied phenomenon and usually accounted for, the effect of warm urban air advected downwind is rarely considered, particularly at rural stations adjacent to urban areas. Until recently, urban heat advection (UHA) was viewed as an urban boundary-layer process through the formation of an urban plume that rises above the surface as it is advected. However, these dynamic UHA effects are shown to also have an impact on surface observations. Results show a significant difference in temperatures anomalies (\(p\,< \,0.001\)) between observations taken downwind of urban and rural areas. For example, urban heat advection from small urbanized areas (\(\sim \)1\(\,\hbox {km}^{2}\)) under low cloud cover and wind speeds of 2–3\(\,\hbox {m}\,\hbox {s}^{-1}\) is found to increase mean nocturnal air temperatures by 0.6\(\,^{\circ }\hbox {C}\) at a horizontal distance of 0.5 km. Fundamentally, these UHA results highlight the importance of careful interpretation of long-term temperature data taken near small urban areas.  相似文献   

3.
Free tropospheric measurements of ozone, peroxyacetylnitrate andprecursors (CO, NMHC) that were made within the framework of the EUROTRACsub-project TOR (Tropospheric Ozone Research) between 1990 and 1995 at theGAW station Izana, Tenerife (28°18N, 16°30W) arediscussed. The average annual cycles reveal the importance of transport fromnorthern mid-latitudes and the role of photo-chemistry. According toair-mass trajectories, which were supplied to us from AEROCE(Atmosphere/Ocean Chemistry Experiment), transport from northernmid-latitudes is associated with high precursor concentrations in winter,whereas ozone concentrations in winter are not much influenced by transportpatterns, suggesting a rather uniform distribution over the northern part ofthe Northern Hemisphere around mean value of 43 ± 5 ppb. In summer,high ozone concentrations of up to 90 ppb are often encountered duringtransport from north, while the levels of precursors are much lower than inwinter, because of photochemical destruction. Trajectories from southerlylatitudes and the Sahara usually have the lowest ozone concentrationsassociated with them.  相似文献   

4.
Vertical dispersion in the neutral surface layer is investigated using a Markov Chain simulation procedure. The conceptual basis of the procedure is discussed and computation procedures outlined. Wind and turbulence parameterizations appropriate to the neutral surface layer are considered with emphasis on the Lagrangian time scale. Computations for a surface release are compared with field data. Good agreement is found for the variation of surface concentration and cloud height to distances 500 m downwind of the source. The functional form of the vertical concentration profile is examined and an exponential with exponent 1.6 is found to give the best fit with simulations.For elevated releases, it is demonstrated that an initial dip of the mass mean height from the simulation can be normalized for various release heights using a non-dimensionalized downwind coordinate incorporating advective wind speed and wind shear. The vertical distribution standard deviation ( z ), as employed in Gaussian models, shows a fair degree of independence with source height but close examination reveals an optimum source height for maximum z at a given downwind distance,x. This source height increases with downwind distance. Also the simulations indicate that vertical wind shear is more important than vertical variation of Lagrangian time scale close to the source, with a reverse effect farther downwind.  相似文献   

5.
Surface-layer aerosol diffusion experiments have been conducted using artificial smoke plume releases at ground level over flat and homogeneously vegetated terrain at the Meppen proving grounds in the Federal Republic of Germany (1989). At fixed downwind locations in the range out to 800 m from the source, instantaneous crosswind plume profiles were detected repetitively at high spatial (1.5 m) and temporal (3 sec) intervals by use of a mini LIDAR system. The experiments were accompanied by measurement of the surface-layer mean wind and turbulence quantities by sonic anemometers. On the basis of measured crosswind concentration profiles, the following statistics were obtained: 1) Mean profile, 2) Root mean square profile, 3) Fluctuation intensities, and 4) Intermittency factors. Furthermore, some experimentally determined probability density functions (pdf's) of the fluctuations are presented. All the measured statistics are referred to a fixed and a moving frame of reference, the latter being defined as a frame of reference from which the (low frequency) plume meander is removed. Finally, the measured statistics are compared with statistics on concentration fluctuations obtained with a simple puff diffusion model (RIMPUFF) developed at Risø.  相似文献   

6.
An examination of typical tropospheric ozone variability on daily, monthly, annual and interannual timescales and instrumental precision indicates that the current ozonesonde network is insufficient to detect a trend in tropospheric ozone of 1% per year at the 2 level even at stations with records a decade in length. From a trend prediction analysis we conclude that in order to detect a 1% per year trend in a decade or less it will be necessary to decrease the time between observations from its present value of 3–7 days to 1 day or less. The spatial distribution of the current ozonesonde stations is also inadequate for determining the global climatology of ozone. We present a quantitative theory taking into account photochemistry, surface deposition, and wind climatology to define the effectively sampled region for an observing station which, used in conjunction with the instrumental precision and the above prediction analysis, forms the basis for defining a suitable global network for determining regional and global ozone climatology and trends. At least a doubling of the present number of stations is necessary, and the oceans, most of Asia, Africa, and South America are areas where more stations are most needed. Differential absorption lidar ozone instruments have the potential for far more frequent measurements of ozone vertical profiles and hence potentially more accurate climatology and trend determinations than feasible with ozonesondes but may produce a (fair weather) biased data set above the cloud base. A strategy for cloudy regions in which either each station utilizes both lidars and sondes or each station is in fact a doublet comprised of a near-sea-level lidar and a proximal-mountain-top lidar could serve to minimize this bias.  相似文献   

7.
A new approach is proposed to predict concentration fluctuations in the framework of one-particle Lagrangian stochastic models. The approach is innovative since it allows the computation of concentration fluctuations in dispersing plumes using a Lagrangian one-particle model with micromixing but with no need for the simulating of background particles. The extension of the model for the treatment of chemically reactive plumes is also accomplished and allows the computation of plume-related chemical reactions in a Lagrangian one-particle framework separately from the background chemical reactions, accounting for the effect of concentration fluctuations on chemical reactions in a general, albeit approximate, manner. These characteristics should make the proposed approach an ideal tool for plume-in-grid calculations in chemistry transport models. The results are compared to the wind-tunnel experiments of Fackrell and Robins (J Fluid Mech, 117:1–26, 1982) for plume dispersion in a neutral boundary layer and to the measurements of Legg et al. (Boundary-Layer Meteorol, 35:277–302, 1986) for line source dispersion in and above a model canopy. Preliminary reacting plume simulations are also shown comparing the model with the experimental results of Brown and Bilger (J Fluid Mech, 312:373–407, 1996; Atmos Environ, 32:611–628, 1998) to demonstrate the feasibility of computing chemical reactions in the proposed framework.  相似文献   

8.
Summary In this publication first results of an urban tracer experiment are reported. This experiment was realized in the framework of the Basel UrBan Boundary-Layer Experiment (BUBBLE) in an area with abundant information on turbulence and flow conditions available. Release height was close to roof level and so was the height of the concentration samplers. The meteorological conditions during the experiments were mainly convective, but due to the rough character of the underlying surface also the mechanical turbulence was substantial.The concentration distribution is found to be essentially Gaussian in the horizontal plane and some commonly used methods to estimate the plume widths in applied dispersion models are compared to the observations. From measurements at one site downwind of the source it is found that for a near-roof level source, only an insignificant vertical gradient in tracer concentration is present within a street canyon. Using a Lagrangian Particle Dispersion Model the tracer experiments are simulated. It is shown that the exact form of the parameterization for the flow and turbulence structure within the urban roughness sublayer is of great importance for the simulation results. Also the numerical simulation results underline the necessity (and difficulty) to describe the vertical profile of the dissipation rate of turbulent kinetic energy close to an urban surface.  相似文献   

9.
Mesoscale measurements of the vertical dispersion coefficient 2 by using a composite turbulence water tank were validated through a comparison with CONDORS (Convective Diffusion Observed with Remote Sensors) field data, and were analysed with respect to the intensity of the thermal flux, mechanical turbulence, and plume release height.It seems possible to correct the plume z values for different release heights below 0.5zi (zi is the mixing height) by applying an equation expressing the height dependency of turbulence intensity. The downwind distance where the plume's mass centre height approaches its final level was also analysed with respect to the above three parameters, and an empirical equation to estimate the downwinddistance derived.  相似文献   

10.
A series of nearly daily ozone vertical profiles obtained at station T-3 on Fletcher's Ice Island (85°N, 90°W) during the period January-March 1971 shows several significant ozone intrusions into the troposphere. These intrusions are not only associated with enhanced ozone amounts in the stratosphere but also require tropopause folding events to transport ozone into the troposphere. These folds in the Arctic tropopause appear to be capable of contributing significantly to the ozone budget of the Arctic troposphere during the late winter and spring seasons. The importance of tropopause folding for bringing ozone into the troposphere seen in the daily ozone profiles confirms the results found in the Arctic Gas and Aerosol Sampling Program aircraft flights.  相似文献   

11.
An easily-evaluated expression for the dimensionless concentration profile (z/z 0,/z 0, z 0/L) = = cu */kQ (or z 0cu*/kQ) downwind of a continuous ground-level area (or line) source in the stable surface layer is obtained by integrating the diffusion equation using the Shwetz approximation method (c = concentration, Q = source strength, k = von Kárman's constant). The analytical solution compares closely with concentration profiles obtained using a trajectory-simulation model over a useful range of heights, the important discrepancies occurring at the upper edge of the plume. The analytical solution is used to generate predictions of ground-level concentration for the Project Prairie Grass experiments; good agreement with the observations is obtained at all downwind distances (50 to 800 m).  相似文献   

12.
We have studied the evolution of the planetary boundary layer using both oceanic and continental stations. This has been possible through the analysis of twice-daily temperature radiosondes at three different oceanic stations (Point A; Point K; and Azores) and also at three continental stations (Berlin, Trappes and Madrid) all situated in the Northern Hemisphere. We have studied the annual evolution of mixing layers and of temperature inversion levels; an annual evolution presenting a minimum in winter and a maximum in summer has been observed in the continental stations; the oceanic stations present a lesser variation and in the opposite sense. As for the elevated inversion layers, their maximum frequency can be observed in summer at the oceanic stations and in winter at the continental stations. We have shown that considering the studied stations, both the level and the frequency of the elevated inversion layers are similar at 00 h and at 12 h; such a result is important because these layers regulate the intensity of the exchanges between the ground and the free atmosphere. Using the equivalent coefficient, we have determined the annual variation of the vertical exchanges between the surface and different altitudes. At the continental stations, the vertical exchanges are more important in summer than they are in winter; the opposite behaviour occurs at the oceanic stations.  相似文献   

13.
During the Berlin Ozone Experiment BERLIOZ in July–August 1998 quasi-continuous measurements ofC2–C12 nonmethane hydrocarbons (NMHCs) were carried out at 10 sites in and around the city of Berlin using on-line gas-chromatographic systems (GCs) with a temporal resolution of 20–120 minutes. Additional airborne NMHCmeasurements were made using canister sampling on three aircraft and an on-line GC system on a fourth aircraft. The ground based data are analyzed to characterize the different sites and to identify the influence of emissions from Berlin on its surroundings. Benzene mixing ratios at the 4 rural sites were rather low (<0.5 ppbv). Berlin (and the surrounding highway ring) was identified as the main source of anthropogenic NMHCs at Eichstädt and Blossin, whilst other sources were important at the furthermost site Menz. The median toluene/benzene concentration ratio in Berlin was 2.3 ppbv/ppbv, agreeing well with measurements in other German cities. As expected, the ratios at the background sites decreased with increasing distance to Berlin and were usually around one or below. On 20 and 21 July, the three northwesterly sites were situated downwind of Berlin and thus were influenced by its emissions. Considering the distance between the sites and the windspeed, the city plume was observed at reasonable time scales, showing decreasing toluene/benzene ratios of 2.3, 1.6 and 1.3 with increasing distance from Berlin. Isoprene was the only biogenic NMHC measured at BERLIOZ. It was themost abundant compound at the background sites on the hotter days, dominating the local NMHC reactivity with averaged contributions to the total OH loss rate of 51% and 70% at Pabstthum and Blossin, respectively. Emissionratios (relative to CO and to the sum of analysed NMHCs) were derived from airborne measurements. The comparison with an emission inventory suggests traffic-related emissions to be the predominating source of the considered hydrocarbon species. Problems were identified with the emission inventory for propane, ethene and pentanes.  相似文献   

14.
Flux measurements of ozone and water vapour employing the eddy correlation technique were used to determine the surface conductance and canopy conductance to ozone. In the surface conductance to ozone, all surfaces at which ozone is destroyed and the transport process to these surfaces are included. The canopy conductance to ozone represents the ozone uptake of transpiring plant parts. The surface conductance to ozone of the maize crop and the underlying soil was generally larger than the canopy conductance to ozone. This means that beside the uptake by stomata, there was another important ozone sink. Under wet soil surface conditions, the surface conductance and the canopy conductance to ozone coincided. This indicates that the resistance of wet soil and the remaining plant parts (cuticle) to ozone was much larger than the stomatal or soil resistance. On the other hand, under dry soil conditions the conductances differ, largely caused by a variation in the transport process to the soil. The transport of ozone to soil increased with increasing friction velocity (u *) and decreased with increasing atmospheric stability, leaf area index (LAI) or crop height (h). These effects for midday (unstable) conditions were parameterized with an in-crop aerodynamic resistance,r inc in a very straightforward way;r inc=13.9 LAIh/u *+67 (cc.=0.77). If the ozone flux in air pollution models is described with a simple resistance model (Big Leaf model), the extra destruction at the soil should be modelled using an in-crop aerodynamic resistance. For these measurements the ozone flux to the soil was 0–65% of the total ozone flux measured above the crop. Under wet soil conditions, this was less than 20%; under dry soil conditions, this was 30–65%.  相似文献   

15.
The dynamical characteristics of concentration fluctuations in a dispersing plume over the energetic and inertial-convective range of scales of turbulent motion are studied using a multiscale analysis technique that is based on an orthonormal wavelet representation. It is shown that the Haar wavelet concentration spectrum is similar to the Fourier concentration spectrum in that both spectra exhibit an extensive inertial-convective subrange spanning about two decades in frequency, with a scaling exponent of -5/3. Analysis of the statistical properties (e.g., fluctuation intensity, skewness, and kurtosis) of the concentration wavelet coefficients (i.e., the concentration discrete detailed signal) suggests that the small scales are always more intermittent than the large scales. The degree of intermittency increases monotonically with decreasing scale within the inertial-convective subrange, reaching a plateau at the very small scales associated with the beginning of the near-dissipation subrange. The probability density function (pdf) of the concentration discrete detailed signal displays stretched exponential tails with an intermittency exponent (tail slope) q that increases as a , where is the scale or dilation and a is a power-law exponent that is dependent on downwind distance, plume height, and stratification strength with typical values in the range from about 0.25 to 0.35. It is shown that the concentration variance cascade process requires a phase coherency of eddies between different scales at the small-scale end of the inertial-convective subrange.The variation of the concentration wavelet statistics with height above the ground is investigated. The increased mean shear near the ground smooths the fine-scale plume structure for scales within the inertial-convective subrange, producing a weaker spatiotemporal intermittency in the concentration field compared to that measured higher up in the plume. The pdf of the concentration detailed signal at a fixed scale possesses less elongated tails with decreasing height z. The intermittency exponent q is found to decrease roughly linearly with increasing z.Finally, the results of the wavelet decomposition are combined to provide a conceptual model of the turbulent transport, stirring, and mixing regimes in a dispersing plume. The implications of the results for contaminant texture in a plume are discussed.  相似文献   

16.
During the BERLIOZ field phase on 20 July 1998 a 40 km wide ozone-plume 30 to 70 km north of Berlin in the lee of the city was detected. The ozone mixing ratio inside the plume was app. 15 ppb higher than outside, mainly caused by high ozone precursor emissions in Berlin, resulting in a net chemical ozone production of 6.5 ppb h–1, which overcompensates ozone advection of –3.6 ppb h–1 andturbulent diffusion of –1.1 ppb h–1. That means, although moreozone leaves the control volume far in the lee of Berlin than enters it at the leeside cityborder and although turbulent diffusion causes a loss of ozone in the leeside control volume the chemical production inside the volume leads to a net ozone increase. Using a semi-Lagrangian mass budget method to estimate the net ozone production, 5.0 ppb h–1 are calculated for theplume. This means a fraction of about 20% of ozone in the plume is producedby local emissions, therefore called `home made' by the Berlin emissions. For the same area KAMM/DRAIS simulations using an observation based initialisation, results in a net production rate between 4.0 and 6.5 ppbh–1, while the threefold nested EURAD model gives 6.0 ppbh–1. The process analysis indicates in many cases goodagreement (10% or better) between measurements and simulations not only in the ozone concentrations but also with respect to the physical and chemical processes governing the total change. Remaining differences are caused by different resolution in time and space of the models and measurements as well as by errors in the emission calculation.The upwind-downwind differences in PAN concentrations are partly similar to those of ozone, because in the BERLIOZ case they are governed mainly by photochemical production. While in the stable boundary layer at night and windward of Berlin 0.1 to 0.3 ppb are detected, in the centre of the plume at noon concentrations between 0.75 ppb and 1.0 ppb are measured. The O3/PAN ratio is about 80 to 120 and thus due to the relatively lowPAN concentrations significantly higher than found in previous studies. The low PAN formation on 20 July, was mainly restricted by the moderate nonmethane hydrocarbon levels, whereas high PAN concentrations of 3.0 ppb on 21 July, are caused by local production in the boundary layer and by large scale advection aloft.  相似文献   

17.
Local ozone production and loss rates for the arctic free troposphere (58–85° N, 1–6 km, February–May) during the TroposphericOzone Production about the Spring Equinox (TOPSE) campaign were calculated using a constrained photochemical box model. Estimates were made to assess the importance of local photochemical ozone production relative to transport in accounting for the springtime maximum in arctic free tropospheric ozone. Ozone production and loss rates from our diel steady-state box model constrained by median observations were first compared to two point box models, one run to instantaneous steady-state and the other run to diel steady-state. A consistent picture of local ozone photochemistry was derived by all three box models suggesting that differences between the approaches were not critical. Our model-derived ozone production rates increased by a factor of 28 in the 1–3 km layer and a factor of 7 in the 3–6 kmlayer between February and May. The arctic ozone budget required net import of ozone into the arctic free troposphere throughout the campaign; however, the transport term exceeded the photochemical production only in the lower free troposphere (1–3 km) between February and March. Gross ozone production rates were calculated to increase linearly with NOx mixing ratiosup to 300 pptv in February and for NOx mixing ratios up to 500 pptv in May. These NOx limits are an order of magnitude higher thanmedian NOx levels observed, illustrating the strong dependence ofgross ozone production rates on NOx mixing ratios for the majority of theobservations. The threshold NOx mixing ratio needed for netpositive ozone production was also calculated to increase from NOx 10pptv in February to 25 pptv in May, suggesting that the NOx levels needed to sustain net ozone production are lower in winter than spring. This lower NOx threshold explains how wintertime photochemical ozone production can impact the build-up of ozone over winter and early spring. There is also an altitude dependence as the threshold NOx neededto produce net ozone shifts to higher values at lower altitudes. This partly explains the calculation of net ozone destruction for the 1–3 km layerand net ozone production for the 3–6 km layer throughout the campaign.  相似文献   

18.
The effects of deep convection on the potential for forming ozone (ozone production potential) in the free troposphere have been simulated for regions where the trace gas composition is influenced by biomass burning. Cloud dynamical and photochemical simulations based on observations in 1980 and 1985 Brazilian campaigns form the basis of a sensitivity study of the ozone production potential under differing conditions. The photochemical fate of pollutants actually entrained in a cumulus event of August 1985 during NASA/GTE/ABLE 2A (Case 1) is compared to photochemical ozone production that could have occurred if the same storm had been located closer to regions of savanna burning (Case 2) and forest burning (Case 3). In each case studied, the ozone production potential is calculated for a 24-hour period following convective redistribution of ozone precursors and compared to ozone production in the absence of convection. In all cases there is considerably more ozone formed in the middle and upper troposphere when convection has redistributed NOx, hydrocarbons and CO compared to the case of no convection.In the August 1985 ABLE 2A event, entrainment of a layer polluted with biomass burning into a convective squall line changes the free tropospheric cloud outflow column (5–13 km) ozone production potential from net destruction to net production. If it is assumed that the same cloud dynamics occur directly over regions of savanna burning, ozone production rates in the middle and upper troposphere are much greater. Diurnally averaged ozone production following convection may reach 7 ppbv/day averaged over the layer from 5–13 km-compared to typical free tropospheric concentrations of 25–30 ppbv O3 during nonpolluted conditions in ABLE 2A. Convection over a forested region where isoprene as well as hydrocarbons from combustion can be transported into the free troposphere leads to yet higher amounts of ozone production.  相似文献   

19.
We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λp were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L c is the canopy-drag length scale determined as a function of λp and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.  相似文献   

20.
The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.  相似文献   

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

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