首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
This paper presents a cohesive view of the dynamics of ambient O(3) exposure and adverse crop response relationships, coupling the properties of photochemical O(3) production, flux of O(3) from the atmosphere into crop canopies and the crop response per se. The results from two independent approaches ((a) statistical and (b) micrometeorological) were analyzed for understanding cause-effect relationships of the foliar injury responses of tobacco cv Bel-W3 to the exposure dynamics of ambient O(3) concentrations. Similarly, other results from two independent approaches were analyzed in: (1) establishing a micrometeorological relationship between hourly ambient O(3) concentrations and their vertical flux from the air into a natural grassland canopy; and (2) establishing a statistical relationship between hourly ambient O(3) concentrations in long-term, chronic exposures and crop yield reductions. Independent of the approach used, atmospheric conditions appeared to be most conducive and the crop response appeared to be best explained statistically by the cumulative frequency of hourly ambient O(3) concentrations between 50 ppb and 90 ppb (100 and 180 microg m(-3)). In general, this concentration range represents intermediate or moderately enhanced hourly O(3) values in a polluted environment. Further, the diurnal occurrence of this concentration range (often approximately between 0900 and 1600 h in a polluted, agricultural environment) coincided with the optimal CO(2) flux from the atmosphere into the crop canopy, thus high uptake. The frequency of occurrence of hourly O(3) concentrations > 90 ppb (180 microg m(-3)) appeared to be of little importance and such concentrations in general appeared to occur during atmospheric conditions which did not facilitate optimal vertical flux into the crop canopy, thus low uptake. Alternatively, when > 90 ppb (180 microg m(-3)) O(3) concentrations occurred during the 0900-1600 h window, their frequency of occurrence was low in comparison to the 50-90 ppb (100-180 microg m(-3)) range. Based on the overall results, we conclude that if the cumulative frequency of hourly ambient O(3) concentrations between 50-62 ppb (100-124 microg m(-3)) occurred during 53% of the growing season and the corresponding cumulative frequency of hourly O(3) concentrations between 50-74 ppb (100-148 microg m(-3)) occurred during 71% of the growing season, then yield reductions in sensitive crops could be expected, if other factors supporting growth, such as adequate soil moisture are not limiting.  相似文献   

2.
Ambient ozone and crop loss: establishing a cause-effect relationship   总被引:6,自引:0,他引:6  
This paper provides the results of a retrospective mathematical analysis of the US NCLAN (National Crop Loss Assessment Network) open-top chamber data. Some 77% of the 73 crop harvests examined, showed no statistically significant yield differences between NF (non-filtered open-top chamber) and AA (chamberless, ambient air) treatments (no easily discernable chamber effects on yield). However, among these cases only seven acceptable examples showed statistically significant yield reductions in NF compared to the CF (charcoal filtered open-top chamber) treatment. An examination of the combined or cumulative hourly ambient O3 frequency distribution for cases with yield loss in NF compared to a similar match of cases without yield loss showed that the mean, median and the various percentiles were all higher (>/= 3 X) in the former in contrast to the latter scenario. The combined frequency distribution of hourly O3 concentrations for the cases with yield loss in NF were clearly separated from the corresponding distribution with no yield loss, at O3 concentrations > 49 ppb. Univariate linear regressions between various O3 exposure parameters and per cent yield losses in NF showed that the cumulative frequency of occurrence of O3 concentrations between 50 and 87 ppb was the best predictor (adjusted R2 = 0.712 and p = 0.011). This analysis also showed that the frequency distribution of hourly concentrations up to 87 ppb O3 represented a critical point, since the addition of the frequency distributions of > 87 ppb O3 did not improve the R2 values. In fact as the frequency of hourly O3 concentrations included in the regression approached 50-100 ppb, the R2 value decreased substantially and the p value increased inversely. Further, univariate linear regressions between the frequencies of occurrence of various O3 concentrations between 50 and 90 ppb and: (a) cases with no yield difference in NF and (b) cases with yield increase in NF compared to the CF treatment (positive effect) provided no meaningful statistical relationship (adjusted R2 = 0.000) in either category. These results support the basis that additional evaluation of the frequency of occurrence of hourly O3] concentrations between 50 and 87 ppb for cases with the yield reductions could provide a meaningful ambient O3 standard, objective or guideline for vegetation.  相似文献   

3.
In Houston, some of the highest measured 8-hr ozone (O3) peaks are characterized by sudden increases in observed concentrations of at least 40 ppb in 1 hr or 60 ppb in 2 hr. Measurements show that these large hourly changes appear at only a few monitors and span a narrow geographic area, suggesting a spatially heterogeneous field of O3 concentrations. This study assessed whether a regulatory air quality model (AQM) can simulate this observed behavior. The AQM did not reproduce the magnitude or location of some of the highest observed hourly O3 changes, and it also failed to capture the limited spatial extent. On days with measured large hourly changes in O3 concentrations, the AQM predicted high O3 over large regions of Houston, resulting in overpredictions at several monitors. This analysis shows that the model can make high O3, but on these days the predicted spatial field suggests that the model had a different cause. Some observed large hourly changes in O3 concentrations have been linked to random releases of industrial volatile organic compounds (VOCs). In the AQM emission inventory, there are several emission events when an industrial point source increases VOC emissions in excess of 10,000 mol/hr. One instance increased predicted downwind O3 concentrations up to 25 ppb. These results show that the modeling system is responsive to a large VOC release, but the timing and location of the release, and meteorological conditions, are critical requirements. Attainment of the O3 standard requires the use of observational data and AQM predictions. If the large observed hourly changes are indicative of a separate cause of high O3, then the model may not include that cause, which might result in regulators enacting control strategies that could be ineffective.  相似文献   

4.
Urban and non-urban rural ozone (O3) concentrations are high in Bulgaria and often exceed the European AOT40 ecosystem as well as the AOT60 human health standards. This paper presents preliminary estimates to establish background, non-urban O3 concentrations for the southern region of Bulgaria. Ozone concentrations from three distinctly different sites are presented: a mountain site influenced by mountain-valley wind flow; a coastal site influenced by sea-breeze wind flow; and a 1700-m mountain peak site without 'local' wind flow characteristics. The latter offers the best estimate of 46-50 ppb for a background O3 level. The highest non-urban hourly value, 118 ppb, was measured at the mountain-valley site.  相似文献   

5.
In the San Bernardino Mountains of southern California, ozone (O(3)) concentrations have been elevated since the 1950s with peaks reaching 600ppb and summer seasonal averages >100ppb in the 1970s. During that period increased mortality of ponderosa and Jeffrey pines occurred. Between the late 1970s and late1990s, O(3) concentrations decreased with peaks approximately 180ppb and approximately 60ppb seasonal averages. However, since the late 1990s concentrations have not changed. Monitoring during summers of 2002-2006 showed that O(3) concentrations (2-week averages) for individual years were much higher in western sites (58-69ppb) than eastern sites (44-50ppb). Potential O(3) phytotoxicity measured as various exposure indices was very high, reaching SUM00 - 173.5ppmh, SUM60 - 112.7ppmh, W126 - 98.3ppmh, and AOT40 - 75ppmh, representing the highest values reported for mountain areas in North America and Europe.  相似文献   

6.
This paper explores the feasibility of (1) using kriging to predict the monthly mean of daily 7-h mean (0900-1559) O3 concentrations, (2) using kriging to estimate the per cent of hourly mean O3 concentrations equal to or greater than 0.07 ppm (137 microg m(-3)) for a specific month, and (3) developing a quantitative relationship between the monthly mean of the daily 7-h (0900-1559) average O3 concentration and the monthly number of hourly concentrations > or = 0.08p ppm (157 microg m(-3)). We found that kriging can be used to estimate the (1) monthly mean of daily 7-h mean O3 concentrations and (2) the percentage of hourly concentrations for a given month > or = 0.07 ppm when sufficient spatial coverage was available. However, the per cent > or = 0.07 ppm parameter exhibited much greater relative variability than the monthly 7-h exposure index. A strong statistical association was found between the monthly number of occurrences > or = 0.08 ppm and monthly 7-h mean concentrations above 0.05 ppm (98 microg m(-3)). Because of the variability that cumulative indices, such as the monthly percentage of hourly concentrations > or = 0.07 ppm , exhibit from site to site, it appears that whether kriging techniques or mathematical regressions are used to estimate the number of elevated O3 hourly concentrations above selected thresholds, large uncertainties associated with the predicted values will exist. These large uncertainties will make it difficult to accurately estimate vegetation effects caused by ambient levels of O3. However, if a generalized quantitative relationship between repeated occurrences of hourly mean concentrations > or = 0.07 ppm or > or = 0.08 and vegetation effects can be developed, it may be possible, using kriged monthly values accompanied with confidence intervals, to identify those areas where vegetation may be at risk. However, before it will be possible to implement such an approach, researchers will have to better quantify the relationship between realistic O3 exposures and vegetation effects.  相似文献   

7.
Ozone concentrations in Alberta cities typically exhibit a maximum in May (up to 35 ppb) and a minimum in November (as low as 4 ppb). This behaviour is similar to that of rural Alberta O3 concentrations. Annual O3 concentrations at six urban monitoring stations vary from 11 ppb to 22 ppb and are about one-half the values at rural stations. In winter, urban O3 concentrations are always smaller than rural concentrations and the cities act as sinks for O3. Although urban stations do not exceed Canada's maximum acceptable levels of daily (25 ppb) and annual (15 ppb) O3 concentrations as often as rural stations, the frequency is still quite large. Canada's hourly maximum desirable level (50 ppb) is exceeded 11 times more often at the remote (rural) station than at the downtown (urban) stations.  相似文献   

8.
To study individual and combined impacts of two important atmospheric trace gases, CO2 and O3, on C and N cycling in forest ecosystems; a multi-year experiment using a small-scale ponderosa pine (Pinus ponderosa Laws.) seedling/soil/litter system was initiated in April 1998. The experiment was conducted in outdoor, sun-lit chambers where aboveground and belowground ecological processes could be studied in detail. This paper describes the approach and methodology used, and presents preliminary data for the first two growing seasons. CO2 treatments were ambient and elevated (ambient + 280 ppm). O3 treatments were elevated (hourly averages to 159 ppb, cumulative exposure > 60 ppb O3, SUM 06 approximately 10.37 ppm h), and a low control level (nearly all hourly averages <40 ppb. SUM 06 approximately 0.07 ppm h). Significant (P < 0.05) individual and interactive effects occurred with elevated CO2 and elevated O3. Elevated CO2 increased needle-level net photosynthetic rates over both seasons. Following the first season, the highest photosynthetic rates were for trees which had previously received elevated O3 in addition to elevated CO2. Elevated CO2 increased seedling stem diameters, with the greatest increase at low O3. Elevated CO2 decreased current year needle % N in the summer. For 1-year-old needles measured in the fall there was a decrease in % N with elevated CO2 at low O3, but an increase in % N with elevated CO2 at elevated O3. Nitrogen fixation (measured by acetylene reduction) was low in ponderosa pine litter and there were no significant CO2 or O3 effects. Neither elevated CO2 nor elevated O3 affected standing root biomass or root length density. Elevated O3 decreased the % N in coarse-fine (1-2 mm diameter) but not in fine (< 1 mm diameter) roots. Both elevated CO2 and elevated O3 tended to increase the number of fungal colony forming units (CFUs) in the AC soil horizon, and elevated O3 tended to decrease bacterial CFUs in the C soil horizon. Thus, after two growing seasons we showed interactive effects of O3 and CO2 in combination, in addition to responses to CO2 or O3 alone for a ponderosa pine plant/litter/soil system.  相似文献   

9.
Growth season-based time series spectral coherence analysis was performed between weekly changes in hourly ambient O(3) concentrations and weekly changes in alfalfa height growth. Weekly median hourly O(3) concentration and the corresponding weekly cumulative integral (sum of all hourly concentrations within the week) were used as indicators of weekly O(3) spectral density and coherence with the change in weekly alfalfa height growth. In general, the weekly cumulative integral performed much better than the weekly median O(3) concentration. A conceptual analysis of the results is presented, along with a recommendation that crop growth stage-based cumulative integrals merit further evaluation towards a better understanding of cause-effect relationships.  相似文献   

10.
Meteorologically adjusted ozone (O3) concentrations during five recent O3 seasons (1998-2002) were computed for six Kentucky metro areas using a nonlinear regression model originally developed for forecasting ground-level O3 concentrations. The meteorological adjustment procedure was based on modifying actual measured O3 concentrations according to model-predicted responses to climate departures with respect to a reference year. For all six Kentucky metro areas, meteorologically adjusted O3 concentrations declined over the five-year period. The linear best-fit rate of decline in mean adjusted O3 concentrations ranged from 0.9 to 2.6 ppb/yr for these metro areas; the average rate of decline was 1.6 ppb/yr. The rates of decline in meteorologically adjusted extreme value (e.g., 95th percentile) concentrations were approximately the same, but there is greater statistical uncertainty concerning the extreme value trends.  相似文献   

11.
Ambient ozone (O(3)) concentrations in the forested areas of the Central and Eastern European (CEE) mountains measured on passive sampler networks and in several locations equipped with active monitors are reviewed. Some areas of the Carpathian Mountains, especially in Romania and parts of Poland, as well as the Sumava and Brdy Mountains in the Czech Republic are characterized by low European background concentrations of the pollutant (summer season means approximately 30 ppb). Other parts of the Carpathians, especially the western part of the range (Slovakia, the Czech Republic and Poland), some of the Eastern (Ukraine) and Southern (Romania) Carpathians and the Jizerske Mountains have high O(3) levels with peak values >100 ppb and seasonal means approximately 50 ppb. Large portions of the CEE mountain forests experience O(3) exposures that are above levels recommended for protection of forest and natural vegetation. Continuation of monitoring efforts with a combination of active monitors and passive samplers is needed for developing risk assessment scenarios for forests and other natural areas of the CEE Region.  相似文献   

12.
There is an ongoing debate as to which components of the ambient ozone (O3) exposure dynamics best explain adverse crop yield responses. A key issue is regarding the importance of peak versus mid-range hourly ambient O3 concentrations. While in this paper the importance of peak atmospheric O3 concentrations is not discounted, if they occur at a time when plants are conducive for uptake, the corresponding importance of more frequently occurring mid-range O3 concentrations is described. The probability of co-occurrence of high O3 concentrations and O3 uptake limiting factors is provided using coherent data sets of O3 concentration, air temperature, air humidity, mean horizontal wind velocity and global radiation measured at representative US and German air quality monitoring sites. Using the PLant-ATmosphere INteraction (PLATIN) model, the significance of the aforementioned meteorological parameters on ozone uptake is examined. In addition, the limitations of describing the O3 exposure for plants under ambient, chamberless conditions by SUM06, AOT40 or W126 exposure indices are discussed.  相似文献   

13.
The responses of ramets of hybrid poplar (Populus spp.) (HP) clones NE388 and NE359, and seedlings of red maple (Acer rubrum, L.) to ambient ozone (O(3)) were studied during May-September of 2000 and 2001 under natural forest conditions and differing natural sunlight exposures (sun, partial shade and full shade). Ambient O(3) concentrations at the study site reached hourly peaks of 109 and 98 ppb in 2000 and 2001, respectively. Monthly 12-h average O(3) concentrations ranged from 32.3 to 52.9 ppb. Weekly 12-h average photosynthetically active radiation (PAR) within the sun, partial shade and full shade plots ranged from 200 to 750, 50 to 180, and 25 to 75 micromol m(-2) s(-1), respectively. Ambient O(3) exposure induced visible foliar symptoms on HP NE388 and NE359 in both growing seasons, with more severe injury observed on NE388 than on NE359. Slight foliar symptoms were observed on red maple seedlings during the 2001 growing season. Percentage of total leaf area affected (%LAA) was positively correlated with cumulative O(3) exposures. More severe foliar injury was observed on plants grown within the full shade and partial shade plots than those observed on plants grown within the sun plot. Lower light availability within the partial shade and full shade plots significantly decreased net photosynthetic rate (Pn) and stomatal conductance (g(wv)). The reductions in Pn were greater than reductions in g(wv), which resulted in greater O(3) uptake per unit Pn in plants grown within the partial shade and full shade plots. Greater O(3) uptake per unit Pn was consistently associated with more severe visible foliar injury in all species and/or clones regardless of differences in shade tolerance. These studies suggest that plant physiological responses to O(3) exposure are likely complicated due to multiple factors under natural forest conditions.  相似文献   

14.
Chongqing, a city with a population of >2.5 million, constitutes the biggest industrial and commercial centre in southwest China. Recent industrialization has led to an increasing air pollutant problem which is exacerbated by the topography and prevailing climate of the region. To date, interest has remained firmly focused on the levels of traditional air pollutants (sulphur dioxide [SO2], oxides of nitrogen [NOx], smoke and suspended particulate matter [SPM]), with little attention paid to photochemical oxidants such as ozone (O3). This paper reports the first assessment of ambient O3 levels in southwest China. Measurements were made in and around Chongqing using a combination of UV-absorption (at a site located in the northern sector of the city) and passive samplers (at 20 sites located in and around the city) between 1993 and 1996. The 7-h daily mean O3 concentrations ranged between 2 and 16 ppb (x10(9)) during the winter months, increasing to 18-41 ppb during the summer (June-August), when peak hourly mean O3 concentrations of 93 ppb were attained. Ozone exposures across the region commonly approached (or exceeded) UN-ECE and WHO short-term guidelines for the protection of crops. In addition, controlled chamber studies, in which 11 cultivars of Chinese crops commonly grown in the Chongqing region were screened for relative O3 sensitivity, indicated the potential for subtle effects on the growth of a number of crop plants, if ambient O3 levels continue to rise in the region. Employing ozone exposures somewhat higher than those experienced in the field, several cultivars of commonly grown Chinese cereal, vegetable and salad crops were found to be sensitive to O3 in terms of growth, but this was not always associated with the appearance of visible symptoms of injury and, in contrast to what was generally expected, only three species showed significant O3-induced reductions in root:shoot partitioning. There is a clear and urgent need for a comprehensive study of ambient air quality and its impact on crops and natural vegetation in this, as in other, rapidly developing regions of China.  相似文献   

15.
16.
The U.S. Environmental Protection Agency has established a federal reference method (FRM) for ozone (O3) and allowed for designation of federal equivalent methods (FEMs). However, the ethylene-chemiluminescence FRM for O3 has been replaced by the UV photometric FEM by most state and local monitoring agencies because of its relative ease of operation. Accumulating evidence indicates that the FEM is prone to bias under the hot, humid, and stagnant conditions conducive to high O3 formation. This bias may lead to overreporting hourly O3 concentrations by as much as 20-40 ppb. Measurement bias is caused by contamination of the O3 scrubber, a problem that is not detected by dry air calibration. An adequate wet test has not been codified, although a procedure has been proposed for agency consideration. This paper includes documentation of laboratory tests quantifying specific interferant responses, collocated ambient FRM/FEM monitoring results, and smog chamber comparisons of the FRM and FEMs with alternative scrubber designs. As the numbers of reports on monitor interferences have grown, interested parties have called for agency recognition and correction of these biases.  相似文献   

17.
The objective of this study is to compare the use of several indices of exposure in describing the relationship between O3 and reduction in agricultural crop yield. No attempt has been made to determine which exposure-response models best fit the data sets examined. Hourly mean O3 concentration data, based on two-three measurements per hour, were used to develop indices of exposure from soybean and winter wheat experiments conducted in open-top chambers at the Boyce Thompson Institute, Ithaca, New York NCLAN field site. The comparative efficacy of cumulative indices (i.e. number of occurrences equal to or above specific hourly mean concentrations, sum of all hourly mean concentrations equal to or above a selected level, and the weighted sum of all hourly mean concentrations) and means calculated over an experimental period to describe the relationship between exposure to O3 and reductions in the yield of agricultural crops was evaluated. None of the exposure indices consistently provided a best fit with the Weibull and linear models tested. The selection of the model appears to be important in determining the indices that best describe the relationship between exposure and response. The focus of selecting a model should be on fitting the data points as well as on adequately describing biological responses. The investigator should be careful to couple the model with data points derived from indices relevant to the length of exposure. While we have used a small number of data sets, our analysis indicates that exposure indices that weight peak concentrations differently than lower concentrations of an exposure regime can be used in the development of exposure-response functions. Because such indices may have merit from a regulatory perspective, we recommend that additional data sets be used in further analyses to explore the biological rationale for various indices of exposure and their use in exposure-response functions.  相似文献   

18.
The deterministic modeling of ambient O3 concentrations is difficult because of the complexity of the atmospheric system in terms of the number of chemical species; the availability of accurate, time-resolved emissions data; and the required rate constants. However, other complex systems have been successfully approximated using artificial neural networks (ANNs). In this paper, ANNs are used to model and predict ambient O3 concentrations based on a limited number of measured hydrocarbon species, NOx compounds, temperature, and radiant energy. In order to examine the utility of these approaches, data from the Coastal Oxidant Assessment for Southeast Texas (COAST) program in Houston, TX, have been used. In this study, 53 hydrocarbon compounds, along with O3, nitrogen oxides, and meteorological data were continuously measured during summer 1993. Steady-state ANN models were developed to examine the ability of these models to predict current O3 concentrations from measured VOC and NOx concentrations. To predict the future concentrations of O3, dynamic models were also explored and were used for extraction of chemical information such as reactivity estimations for the VOC species. The steady-state model produced an approximation of O3 data and demonstrated the functional relationship between O3 and VOC-NOx concentrations. The dynamic models were able to the adequately predict the O3 concentration and behavior of VOC-NOx-O3 system a number of hourly intervals into the future. For 3 hr into the future, O3 concentration could be predicted with a root-mean squared error (RMSE) of 8.21 ppb. Extending the models further in time led to an RMSE of 11.46 ppb for 5-hr-ahead values. This prediction capability could be useful in determining when control actions are needed to maintain measured concentrations within acceptable value ranges.  相似文献   

19.
Ozone (O3) concentrations were monitored during the 1997-1999 growing seasons in 32 forest sites of the Carpathian Mountains. At all sites (elevation between 450 and 1320 m) concentrations of O3, nitrogen dioxide (NO2), and sulfur dioxide (SO2) were measured with passive samplers. In addition, in two western Carpathian locations, Vychodna and Gubalówka, ozone was continuously monitored with ultraviolet (UV) absorption monitors. Highest average hourly O3 concentrations in the Vychodna and Guba?ówka sites reached 160 and 200 microg/m3 (82 and 102 ppb), respectively (except for the AOT40 values, ozone concentrations are presented as microg/m3; and at 25 degrees C and 760 mm Hg, 1 microg O3/m3 = 0.51 ppb O3). These sites showed drastically different patterns of diurnal 03 distribution, one with clearly defined peaks in the afternoon and lowest values in the morning, the other with flat patterns during the entire 24-h period. On two elevational transects, no effect of elevation on O3 levels was seen on the first one, while on the other a significant increase of O3 levels with elevation occurred. Concentrations of O3 determined with passive samplers were significantly different between individual monitoring years, monitoring periods, and geographic location of the monitoring sites. Results of passive sampler monitoring showed that high O3 concentrations could be expected in many parts of the Carpathian range, especially in its western part, but also in the eastern and southern ranges. More than four-fold denser network of monitoring sites is required for reliable estimates of O3 distribution in forests over the entire Carpathian range (140 points). Potential phytotoxic effects of O3 on forest trees and understory vegetation are expected on almost the entire territory of the Carpathian Mountains. This assumption is based on estimates of the AOT40 indices for forest trees and natural vegetation. Concentrations of NO2 and SO2 in the entire Carpathian range were typical for this part of Europe and below the expected levels of phytotoxicity.  相似文献   

20.
Ambient concentrations of ozone (O(3)) were measured and O(3) phytotoxicity to tobacco (Nicotiana tabacum L.) was demonstrated in several forest locations in Poland during a pilot study from July-October, 1991. At southern and central locations in Poland, the 24-hour average O(3) concentrations measured with a UV absorption photometer were in the range of 32-55 ppb, and the corresponding 1-hour maxima in the range of 39-83 ppb. At these locations longer period (four to fifteen days) average concentrations were determined using O(3) passive samplers (DGA, Inc.) and were reaching 60 ppb, while at Bialowieza in eastern Poland O(3) concentrations averaged less than 40 ppb. In Szarow, near the Niepolomice Forest in southern Poland, 1-hour O(3) maxima estimated from the data obtained using passive samplers were about 105 ppb in early September. At several locations in southern and central Poland, extensive O(3) injury was determined on O(3)-sensitive Bel W-3 tobacco plants; such injury did not occur in the Bialowieza Forest of eastern Poland. The results of this pilot study indicate that O(3) is present at phytotoxic levels in southern and central Poland.  相似文献   

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

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