首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Clouds affect local surface UV irradiance, even if the horizontal distance from the radiation observation site amounts to several kilometers. In order to investigate this effect, which we call remote clouds effect, a 3-dimensional radiative transfer model is applied. Assuming the atmosphere is subdivided into a quadratic based sector and its surrounding, we quantify the influence of changing cloud coverage within this surrounding from 0% to 100% on surface UV irradiance at the sector center. To work out this remote clouds influence as a function of sector base size, we made some calculations for different sizes between 10 km × 10 km and 100 km × 100 km. It appears that in the case of small sectors (base size 20 km × 20 km) the remote clouds effect is highly variable: Depending on cloud structure, solar zenith angle and wavelength, the surface UV irradiance may be enhanced up to 15% as well as reduced by more than 50%. In contrast, for larger sectors it is always the case that enhancements become smaller by 5% if sector base size exceeds 60 km × 60 km. However, these values are upper estimates of the remote cloud effects and they are found only for special cloud structures. Since these structures might occur but cannot be regarded as typical, different satellite observed cloud formations (horizontal resolution about 1 km × 1 km) have also been investigated. For these more common cloud distributions we find remote cloud effects to be distinctly smaller than the corresponding upper estimates, e.g., for a sector with base size of 25 km × 25 km the surface UV irradiance error due to ignoring the actual remote clouds and replacing their influence with periodic horizontal boundary conditions is less than 3%, whereas the upper estimate of remote clouds effect would suggest an error close to 10%.  相似文献   

2.
The effect of present-day and future NOx emissions from aircraft on the NOx and ozone concentrations in the atmosphere and the corresponding radiative forcing were studied using a three-dimensional chemistry transport model (CTM) and a radiative model. The effects of the aircraft emissions were compared with the effects of the three most important anthropogenic NOx surface sources: road traffic, electricity generation and industrial combustion. From the model results, NOx emissions from aircraft are seen to cause an increase in the NOx and ozone concentrations in the upper troposphere and lower stratosphere, and a positive radiative forcing. For the reference year 1990, the aircraft emissions result in an increase in the NOx concentration at 250 hPa of about 20 ppt in January and 50 ppt in July over the eastern USA, the North Atlantic Flight Corridor and Western Europe, corresponding to a relative increase of about 50%. The maximum increase in the ozone concentrations due to the aircraft emissions is about 3-4 ppb in July over the northern mid-latitudes, corresponding to a relative increase of about 3-4%. The aircraft-induced ozone changes cause a global average radiative forcing of 0.025 W/m2 in July. According to the ANCAT projection for the year 2015, the aircraft NOx emissions in that year will be 90% higher than in the year 1990. As a consequence of this, the calculated NOx perturbation by aircraft emissions increases by about 90% between 1990 and 2015, and the ozone perturbation by about 50-70%. The global average radiative forcing due to the aircraft-induced ozone changes increases by about 50% between 1990 and 2015. In the year 2015, the effects of the aircraft emissions on the ozone burden and radiative forcing are clearly larger than the individual effects of the NOx surface sources. Taking chemical conversion in the aircraft plume into account in the CTM explicitly, by means of modified aircraft NOx emissions, a significant reduction of the aircraft-induced NOx and ozone perturbations is realised. The NOx perturbation decreases by about 40% and the ozone perturbation by about 30% in July over Western Europe, the eastern USA and the North Atlantic Flight Corridor.  相似文献   

3.
A parametric study of the instantaneous radiative impact of contrails is presented using three different radiative transfer models for a series of model atmospheres and cloud parameters. Contrails are treated as geometrically and optically thin plane parallel homogeneous cirrus layers in a static atmosphere. The ice water content is varied as a function of ambient temperature. The model atmospheres include tropical, mid-latitude, and subarctic summer and winter atmospheres. Optically thin contrails cause a positive net forcing at top of the atmosphere. At the surface the radiative forcing is negative during daytime. The forcing increases with the optical depth and the amount of contrail cover. At the top of the atmosphere, a mean contrail cover of 0.1% with average optical depth of 0.2 to 0.5 causes about 0.01 to 0.03 Wm−2 daily mean instantaneous radiative forcing. Contrails cool the surface during the day and heat the surface during the night, and hence reduce the daily temperature amplitude. The net effect depends strongly on the daily variation of contrail cloud cover. The indirect radiative forcing due to particle changes in natural cirrus clouds may be of the same magnitude as the direct one due to additional cover.  相似文献   

4.
Erythemally-weighted ultraviolet (UVery) levels measured over southern England, during anticyclonic weather between 30 April and 2 May, 1997, were almost 50 higher than normally expected for clear skies and were similar to mid-summer values for the first time since measurements began in 1990. Investigation of this episode suggests that a combination of both meteorological and chemical effects were responsible for generating record low ozone amounts for the time of year. Further, comparisons between the A band ultraviolet (315 to 400 nm wavelength) amounts, and radiative calculations confirm that the high UVery was primarily due to the reduction in total ozone. These results are contrasted with a similar period for 1998, in which near climatological ozone amounts were measured. The prospects for enhanced UVery levels in future years are briefly reviewed in the light of expected increases in stratospheric halogen levels and greenhouse gases.  相似文献   

5.
This study demonstrates that ordinary kriging in spherical coordinates using experimental semi-variograms provides highly usable results, especially near the pole in winter and/or where there could be data missing over large areas. In addition, kriging allows display of the spatial variability of daily ozone measurements at different pressure levels. Three satellite data sets were used: Total Ozone Mapping Spectrometer (TOMS) data, Solar Backscattered UltraViolet (SBUV), and the Stratospheric Aerosol and Gas Experiment (SAGE II) ozone profiles. Since SBUV is a nadir-viewing instrument, measurements are only taken along the sun-synchronous polar orbits of the satellite. SAGE II is a limb-viewing solar occultation instrument, and measurements have high vertical resolution but poor daily coverage. TOMS has wider coverage with equidistant distribution of data (resolution 1° × 1.25°) but provides no vertical information. Comparisons of the resulting SBUV-interpolated (column-integrated) ozone field with TOMS data are strongly in agreement, with a global correlation of close to 98%. Comparisons of SBUV-interpolated ozone profiles with daily SAGE II profiles are relatively good, and comparable to those found in the literature. The interpolated ozone layers at different pressure levels are shown.  相似文献   

6.
With a detailed chemistry scheme for the middle atmosphere up to 70 km which has been added to the 3-D Karlsruhe simulation model of the middle atmosphere (KASIMA), the effects of coupling chemistry and dynamics through ozone are studied for the middle atmosphere. An uncoupled version using an ozone climatology for determining heating rates and a coupled version using on-line ozone are compared in a 10-month integration with meteorological analyses for the winter 1992/93 as the lower boundary condition. Both versions simulate the meteorological situation satisfactorily, but exhibit a too cold lower stratosphere. The on-line ozone differs from the climatological data between 20 and 40 km by exhibiting too high ozone values, whereas in the lower mesosphere the ozone values are too low. The coupled model version is stable and differs only above 40 km significantly from the uncoupled version. Direct heating effects are identified to cause most of the differences. The well-known negative correlation between temperature and ozone is reproduced in the model. As a result, the coupled version slightly approaches the climatological ozone field. Further feedback effects are studied by using the on-line ozone field as a basis for an artificial climatology. For non-disturbed ozone conditions realistic monthly and zonally averaged ozone data are sufficient to determine the heating rates for modelling the middle atmosphere.  相似文献   

7.
Total ozone anomalies (deviation from the long-term mean) are created by anomalous circulation patterns. The dynamically produced ozone anomalies can be estimated from known circulation parameters in the layer between the tropopause and the middle stratosphere by means of statistics. Satellite observations of ozone anomalies can be compared with those expected from dynamics. Residual negative anomalies may be due to chemical ozone destruction. The statistics are derived from a 14 year data set of TOMS (Total Ozone Mapping Spectrometer January 1979-Dec. 1992) and corresponding 300 hPa geopotential (for the tropopause height) together with 30 hPa temperature (for stratospheric waves) at 60°N. The correlation coefficient for the linear multiple regression between total ozone (dependent variable) and the dynamical parameters (independent variables) is 0.88 for the zonal deviations in the winter of the Northern Hemisphere. Zonal means are also significantly dependent on circulation parameters, besides showing the known negative trend function of total ozone observed by TOMS. The significant linear trend for 60°N is 3 DU/year in the winter months taking into account the dependence on the dynamics between the tropopause region and the mid-stratosphere. The highest correlation coefficient for the monthly mean total ozone anomalies is reached in November with 0.94.  相似文献   

8.
A European campaign of ground-based radar, lidar and optical measurements was carried out during the winter of 1996/1997 (28 December–2 February) to study lee waves in the northern part of Scandinavia. The participants operated ozone lidars, backscatter lidars and MST radars at ALOMAR/Andoya and Esrange/Kiruna, and an ALIS imaging system in Kiruna. The Andoya site was generally windward of the Scandinavian mountains, the Kiruna site on the leeward side. The goal of the experiment was to examine the influence of lee waves on the formation of Polar Stratospheric Clouds (PSCs). This paper studies the radar data from MST-radar ESRAD located at Esrange [68.°N, 21.°E], i.e. in the lee of the mountains. We present three cases where strong lee waves were observed: in one case they propagated upwards to the lower stratosphere and in the other two cases they were trapped or absorbed in the troposphere. We examine the local waves and the direction and strength of the local wind using the radar, the synoptic meteorological situation using weather maps (European Meteorological Bulletin) and the synoptic stratospheric temperatures using ECMWF data. We observed that waves propagate up to the stratosphere during frontal passages. When anticyclonic ridges are present, the propagation to the stratosphere is very weak. This is due to trapping of the waves at or below the tropopause. We also show that the radar data alone can be used to characterise the different weather conditions for the three cases studied (through the variation of the height of the tropopause). The synoptic stratospheric temperatures in the three cases were similar, and were above the expected threshold for PSC formation. Lidar and visual observation of PSCs and nacreous clouds, respectively, showed that these were present only in the case when the lee waves propagated up to the lower stratosphere.  相似文献   

9.
Observations of traveling ionospheric disturbances (TIDs) associated with atmospheric gravity waves (AGWs) generated by the moving solar terminator have been made with the Millstone Hill incoherent scatter radar. Three experiments near 1995 fall equinox measured the AGW/TID velocity and direction of motion. Spectral and cross-correlation analysis of the ionospheric density observations indicates that ST-generated AGWs/TIDs were observed during each experiment, with the more-pronounced effect occurring at sunrise. The strongest oscillations in the ionospheric parameters have periods of 1.5 to 2 hours. The group and phase velocities have been determined and show that the disturbances propagate in the horizontal plane perpendicular to the terminator with the group velocity of 300–400 m s–1 that corresponds to the ST speed at ionospheric heights. The high horizontal group velocity seems to contradict the accepted theory of AGW/TID propagation and indicates a need for additional investigation.  相似文献   

10.
The results of three series of rocket measurements of mesospheric electric fields carried out under different geomagnetic conditions at polar and high middle latitudes are analysed. The measurements show a clear dependence of the vertical electric fields on geomagnetic activity at polar and high middle latitudes. The vertical electric fields in the lower mesosphere increase with the increase of geomagnetic indexes Kp and Kp. The simultaneous increase of the vertical electric field strength and ion conductivity was observed in the mesosphere during geomagnetic disturbances. This striking phenomenon was displayed most clearly during the solar proton events of October, 1989 accompanied by very strong geomagnetic storm (Kp = 8+). A possible mechanism of generation of the vertical electric fields in the mesosphere caused by gravitational sedimentation of charged aerosol particles is discussed. Simultaneous existence in the mesosphere of both the negative and positive multiply charged aerosol particles of different sizes is assumed for explanation of the observed V/m vertical electric fields and their behaviour under geomagnetically disturbed conditions.Paper Presented at the Second IAGA/ICMA (IAMAS) Workshop on Solar Activity Forcing of the Middle Atmosphere, Prague, August 1997  相似文献   

11.
The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.  相似文献   

12.
Orographic enhancement of wet deposition arising from the seeder-feeder effect is, by necessity, highly parametrised in long-range transport models of acid deposition that are long-term (i.e. annual average) and spatially resolved at tens of kilometres. Here, we describe a mechanistic approach to the incorporation of these mechanisms into such a model. The model formulation required the following: precipitation rate by direction and quantification of the fractions that are orographic and non-orographic; treatment of the fast oxidation of sulfur dioxide in clouds; the directionality of the seeder-feeder process; and a quantitative basis for increasing wet deposition factors to account for the seeder-feeder process. The directionality of non-orographic precipitation was determined from meteorological data at 47 sites across the UK. Orographic precipitation varies on a much finer scale than can be interpolated from measurements, and thus a modelling approach was adopted. The directionality of the seeder-feeder effect was taken from measurements. The enhancement factor of the orographic component of precipitation, assumed to represent feeder-rain, was determined from a review of measurements. Fast oxidation of sulfur dioxide is an observed phenomenon in cap-cloud, but limited in duration. An adjustment was made to the sulfur dioxide oxidation rate in the model in locations where cap-cloud was assumed to be present. The results from the model were compared with UK deposition budgets and enhanced wet deposition maps. The Revised parametrisation underestimated the UK wet deposition budgets of oxidised N and S, but spatial patterns of deposition were improved for much of the UK. It was concluded that this was a satisfactory outcome given the constraints of the statistical approach of weighting of deposition at receptors utilising straight line trajectories. The sensitivity of the model to directional constraints of seeder-feeder enhancement was tested and it was concluded that a fairly narrow constraint resulted in similar estimations to a broader one, and the broader constraint was thus adopted as frontal conditions which result in the process arrive from a fairly broad band of directions. When enhancement was allowed to occur from all directions, UK wet deposition of oxidised N and S was increased by 10%. The sensitivity to the enhancement factor on wet deposition was tested and found to be relatively robust. An increase in the enhancement factor from 2 to 6 resulted in increases in UK wet deposition of oxidised N and S of 9 and 6%, respectively.  相似文献   

13.
Polar stratospheric clouds (PSCs) are often observed in the Kiruna region in northern Sweden, east of the Scandinavian mountain range, during wintertime. PSC occurrence can be detected by ground-based optical instruments. Most of these require clear tropospheric weather. By applying the zenith-sky colour index technique, which works under most weather conditions, the data availability can be extended. The observations suggest that PSC events, especially of type II (water PSCs) may indeed more common than predicted by synoptic models, which is expected because of the frequent presence of mountain-induced leewaves. However, it will be of importance to increase the density of independent observations.  相似文献   

14.
Many Large Eddy Simulation (LES) models use the classic Kessler parameterisation either as it is or in a modified form to model the process of cloud water autoconversion into precipitation. The Kessler scheme, being linear, is particularly useful and is computationally straightforward to implement. However, a major limitation with this scheme lies in its inability to predict different autoconversion rates for maritime and continental clouds. In contrast, the Berry formulation overcomes this difficulty, although it is cubic. Due to their different forms, it is difficult to match the two solutions to each other. In this paper we single out the processes of cloud conversion and accretion operating in a deep model cloud and neglect the advection terms for simplicity. This facilitates exact analytical integration and we are able to derive new expressions for the time of onset of precipitation using both the Kessler and Berry formulations. We then discuss the conditions when the two schemes are equivalent. Finally, we also critically examine the process of droplet evaporation within the framework of the classic Kessler scheme. We improve the existing parameterisation with an accurate estimation of the diffusional mass transport of water vapour. We then demonstrate the overall robustness of our calculations by comparing our results with the experimental observations of Beard and Pruppacher, and find excellent agreement.  相似文献   

15.
Wind and temperature profiles measured routinely by rockets at Ryori (Japan) since 1970 are analysed to quantify interannual changes that occur in the upper stratosphere. The analysis involved using a least square fitting of the data with a multiparametric adaptative model composed of a linear combination of some functions that represent the main expected climate forcing responses of the stratosphere. These functions are seasonal cycles, solar activity changes, stratospheric optical depth induced by volcanic aerosols, equatorial wind oscillations and a possible linear trend. Step functions are also included in the analyses to take into account instrumental changes. Results reveal a small change for wind data series above 45 km when new corrections were introduced to take into account instrumental changes. However, no significant change of the mean is noted for temperature even after sondes were improved. While wind series reveal no significant trends, a significant cooling of 2.0 to 2.5 K/decade is observed in the mid upper stratosphere using this analysis method. This cooling is more than double the cooling predicted by models by a factor of more than two. In winter, it may be noted that the amplitude of the atmospheric response is enhanced. This is probably caused by the larger ozone depletion and/or by some dynamical feedback effects. In winter, cooling tends to be smaller around 40–45 km (in fact a warming trend is observed in December) as already observed in other data sets and simulated by models. Although the winter response to volcanic aerosols is in good agreement with numerical simulations, the solar signature is of the opposite sign to that expected. This is not understood, but it has already been observed with other data sets.  相似文献   

16.
Some important climatic features of the Mediterranean Sea stand out from an analysis of the systematic discrepancies between direct measurements of longwave radiation budget and predictions obtained by the most widely used bulk formulae. In particular, under clear-sky conditions the results show that the surface values of both air temperature and humidity over the Mediterranean Sea are larger than those expected over an open ocean with the same amount of net longwave radiation. Furthermore, the twofold climatic regime of the Mediterranean region strongly affects the downwelling clear-sky radiation. This study suggests that a single bulk formula with constant numerical coefficients is unable to reproduce the fluxes at the surface for all the seasons.  相似文献   

17.
Two Doppler imaging systems (DIS) or wide-field imaging Fabry-Perot interferometers (FPI), have recently been commissioned, one at the Auroral Station, Adventdalen, Longyearbyen, Svalbard, and the second at the IRF, Kiruna, Sweden. These instruments can provide wide-field (600 * 800 km) images of neutral wind flows in the upper thermosphere, by measuring the Doppler shift of the atomic oxygen forbidden near 630 nm, which is emitted from an altitude of approximately 240 km. From the instrument in Svalbard, at mid-winter, it is possible to observe the dayside polar cusp and the polar cap throughout the entire day, whereas from Kiruna, the night-time auroral oval is observable during the hours of darkness. Measurements of thermospheric dynamics from the DIS can be used in conjunction with observations of ionospheric plasma flows and thermal plasma densities by the EISCAT-Svalbard radar (ESR) and by EISCAT, along with other complementary observations by co-located instruments such as the auroral large-scale imaging system (ALIS). Such combined data sets will allow a wide range of scientific studies to be performed concerning the dynamical response of the thermosphere and ionosphere, and the important energetic and momentum exchange processes resulting from their complex interactions. These processes are particularly important in the immediate vicinity of the polar cusp and within the auroral oval. Early results from Svalbard in late 1995 will be discussed. The DIS in Kiruna observed two interesting geomagnetic disturbances in early 1997, the minor geomagnetic storm of 10, 11 January, and the disturbed period from 7–10 February. During these events, the thermospheric wind response showed some interesting departures from the average behaviour, which we attribute to the result of strong and variable Lorenz forcing (ion drag) and Joule and particle heating during these geomagnetic disturbances.  相似文献   

18.
The Étude et Simulation de la QUalité de lair en Ile de France (ESQUIF) project is the first integrated project dedicated to the study of the processes leading to air pollution events over the Paris area. The project was carried out over two years (summer 1998 to winter 2000) to document all types of meteorological conditions favourable to air quality degradation, and in particular to photo oxydant formation. The goals of ESQUIF are (1) to improve our understanding of the relevant chemical and dynamical processes and, in turn, improve their parametrizations in numerical models, and (2) to improve and validate existing models dedicated to pollution analysis, scenarios and/or forecasting, by establishing a comprehensive and thorough database. We present the rationale of the ESQUIF project and we describe the experimental set-up. We also report on the first experiments which took place during the summer of 1998 involving surface networks, and remote sensing instruments as well as several aircraft. Focusing on three days of August 1998, the relative contributions of long-range transported and locally-produced ozone to the elevated ozone concentrations observed during this period are discussed and chemistry-transport model preliminary results on this period are compared to measurements.  相似文献   

19.
As part of the preparation for the Large-Scale Biosphere Atmosphere Experiment in Amazonia, a meso-scale modelling study was executed to highlight deficiencies in the current understanding of land surface atmosphere interaction at local to sub-continental scales in the dry season. Meso-scale models were run in 1-D and 3-D mode for the area of Rondonia State, Brazil. The important conclusions are that without calibration it is difficult to model the energy partitioning of pasture; modelling that of forest is easier due to the absence of a strong moisture deficit signal. The simulation of the boundary layer above forest is good, above deforested areas (pasture) poor. The models’ underestimate of the temperature of the boundary layer is likely to be caused by the neglect of the radiative effects of aerosols caused by biomass burning, but other factors such as lack of sufficient entrainment in the model at the mixed layer top may also contribute. The Andes generate patterns of subsidence and gravity waves, the effects of which are felt far into the Rondonian area The results show that the picture presented by GCM modelling studies may need to be balanced by an increased understanding of what happens at the meso-scale. The results are used to identify key measurements for the LBA atmospheric meso-scale campaign needed to improve the model simulations. Similar modelling studies are proposed for the wet season in Rondonia, when convection plays a major role.  相似文献   

20.
Interhemispheric transport is a key process affecting the accuracy of source quantification for species such as methane by inverse modelling, and is a source of difference among global three-dimensional chemistry transport models (CTMs). Here we use long-term observations of the atmospheric concentration of long-lived species such as CH3CCl3 and CFCl3 for testing three-dimensional chemistry transport models (CTMs); notably their ability to model the interhemispheric transport, distribution, trend, and variability of trace gases in the troposphere. The very striking contrast between the inhomogeneous source distribution and the nearly homogeneous trend, observed in the global ALE/GAGE experiments for both CH3CCl3 and CFCl3 illustrates an efficient interhemispheric transport of atmospherically long-lived chemical species. Analysis of the modelling data at two tropical stations, Barbados (13°N, 59°W) and Samoa (14°S, 124°W), show the close relationship between inter-hemispheric transport and cross-equator Hadley circulations. We found that cross-equator Hadley circulations play a key role in producing the globally homogeneous observed trends. Chemically, the most rapid interaction between CH3CCl3 and OH occurs in the northern summer troposphere; while the most rapid photolysis of CH3CCl3 and CFCl3, and the chemical reactions between CFCl3 and O(1D), take place in the southern summer stratosphere. Therefore, the cross-equator Hadley circulation plays a key role which regulates the southward flux of chemical species. The regulation by the Hadley circulations hence determines the amount of air to be processed by OH, O(1D), and ultraviolet photolysis, in both hemispheres. In summary, the dynamic regulation of the Hadley circulations, and the chemical processing (which crucially depends on the concentration of OH, O(1D), and on the intensity of solar insolation) of the air contribute to the seasonal variability and homogeneous growth rate of observed CH3CCl3 and CFCl3.  相似文献   

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

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