首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
This paper presents the results of a study to investigate the atmospheric oxidation of sulfur dioxide (SO2). A detailed model of gas-phase chemistry, aerosol thermodynamics and aerosol chemistry is employed to simulate atmospheric sulfate formation. The calculations indicate that, in addition to the gasphase oxidation by hydroxyl (OH) radicals, SO2 oxidation in aqueous aerosols may also contribute significantly to sulfate formation. Reactions of SO2 with hydrogen peroxide (H2O2) and O2 (catalyzed by Fe3+ and Mn2+) are identified as principal aqueous-phase oxidation mechanisms. The results of this study confirm the conclusions drawn from the analysis of ambient aerosol data qualitatively. However, some discrepancies also exist between the results of our modeling study and field data. Such discrepancies emphasize the need for the collection of ambient data for a more rigorous and quantitative evaluation of atmospheric aerosol models.  相似文献   

2.
Results of a theoretical investigation of H2O2 formation in cloud droplets arising from gaseous HO2 radical scavenging are presented. It is shown that this process is pH dependent with the maximum rate of H2O2 production occurring below pH 3. This dependence arises as a result of the dissociation of HO2 in water (pKa = 4.9) and the subsequent disproportionation reaction of HO2 and O2 to form hydrogen peroxide. O2 is also removed by reaction with O3 to produce OH radicals and this process becomes more competitive as both the pH and O2HO2 ratio increase. The presence of soluble organic species, such as aldehydes, in cloudwater counteracts the effect of ozone by converting OH back to HO2. For low pHs (< 3) the net contribution of organic solutes of H2O2 production is predicted to be relatively small, being limited by the availability of OH radicals scavenged from the gas phase. Existing cloud chemistry models may overestimate the rate of aqueous oxidation of formaldehyde by OH radicals.Under conditions where scavenging of gas-phase free radicals by cloud droplets is efficient, uptake of HO2 radicals may be reversible. The aqueous concentration of OH is unlikely to approach thermodynamic equilibrium with the gas phase (H ∼-30 M atm−1 and can be treated as irreversible. In clouds with a small mean droplet radius, efficient scavenging of precursor OH radicals should result in a decrease in gas-phase HO2 production with a reduction in the yield of aqueous H2O2, although this is offset by the presence of soluble organic species. A similar effect is predicted for clouds with a high liquid water content.The supply of HO2 and OH radicals to cloud droplets is controlled by gas-phase ozone chemistry which is in turn dependent on the solar u.v. radiation intensity. The u.v. density in clouds may be higher than in clear air when the solar zenith angle is small, thus enhancing H2O2 production, but falls off markedly as the solar zenith angle becomes larger. Predicted rates of H2O2 formation in clouds based on midday conditions are likely to be considerably higher than the average daytime value, particularly in summer. Diurnal and seasonal effects on H2O2 generation are expected to be more marked in clouds than in clear air.  相似文献   

3.
In this study, we present the response of model results to different scientific treatments in an effort to quantify the uncertainties caused by the incomplete understanding of mercury science and by model assumptions in atmospheric mercury models. Two sets of sensitivity simulations were performed to assess the uncertainties using modified versions of CMAQ-Hg in a 36-km Continental United States domain. From Set 1 Experiments, it is found that the simulated mercury dry deposition is most sensitive to the gaseous elemental mercury (GEM) oxidation product assignment, and to the implemented dry deposition scheme for GEM and reactive gaseous mercury (RGM). The simulated wet deposition is sensitive to the aqueous Hg(II) sorption scheme, and to the GEM oxidation product assignment. The inclusion of natural mercury emission causes a small increase in GEM concentration but has little impact on deposition. From Set 2 Experiments, it is found that both dry and wet depositions are sensitive to mercury chemistry. Change in model mercury chemistry has a greater impact on simulated wet deposition than on dry deposition. The kinetic uncertainty of GEM oxidation by O3 and mechanistic uncertainty of Hg(II) reduction by aqueous HO2 pose the greatest impact. Using the upper-limit kinetics of GEM–O3 reaction or eliminating aqueous Hg(II)–HO2 reaction results in unreasonably high deposition and depletion of gaseous mercury in the domain. Removing GEM–OH reaction is not sufficient to balance the excessive mercury removal caused by eliminating the HO2 mechanism. Field measurements of mercury dry deposition, better quantification of mercury air-surface exchange and further investigation of mercury redox chemistry are needed for reducing model uncertainties and for improving the performance of atmospheric mercury models.  相似文献   

4.
The aqueous generation of H2O2 from bubbling ozonized air through water in glass impingers is investigated. O3 loss is monitored throughout each experiment. Aqueous H2O2 is measured at the end of each experiment. The mole quantities of O3 lost (ΔO3) and H2O2 formed (ΔH2O2) are calculated. There is a stoichiometric relationship between ΔO3 and ΔH2O2 which is independent of bubbling time, solution acidity and initial O3 concentration but is dependent upon the liquid water content of each impinger. For 10 ml of water in each impinger, the AO3 to ΔH2O2 mole ratio is between 2:1 and 3:1 and for 20 ml is between 1:1 and 1:2. Increasing the solution acidity increases the rates of O3 loss and H2O2 generation. This result contradicts mechanisms in which a O3 + HO reaction is used to initiate radical formation and H2O2 production.These results are interpreted through a combination of glass surface chemistry and bulk aqueous chemistry. The O3-H2O2 impinger mechanisms involve O3 adsorption and/or reaction on the impinger surface and its surface products subsequent reaction with water or trace impurities in the water. Experiments using tetrafluorethylene (TFE) impingers in place of glass impingers are conducted to examine the role of surface chemistry. The TFE impingers are found to substantially increase ΔO3 and ΔH2O2 without changing their mole ratio. This is attributed to an enhanced O3-H2O surface chemistry or O3 reaction with unsaturated Teflon linkages or impurities in the Teflon from machining.  相似文献   

5.
An interpretative modeling analysis is conducted to simulate the diurnal variations in OH and HO2+RO2 observed at Summit, Greenland in 2003. The main goal is to assess the HOx budget and to quantify the impact of snow emissions on ambient HOx as well as on CH2O and H2O2. This analysis is based on composite diurnal profiles of HOx precursors recorded during a 3-day period (July 7–9), which were generally compatible with values reported in earlier studies. The model simulations can reproduce the observed diurnal variation in HO2+RO2 when they are constrained by observations of H2O2 and CH2O. By contrast, model predictions of OH were about factor of 2 higher than the observed values. Modeling analysis of H2O2 suggests that its distinct diurnal variation is likely controlled by snow emissions and loss by deposition and/or scavenging. Similarly, deposition and/or scavenging sinks are needed to reproduce the observed diel profile in CH2O. This study suggests that for the Summit 2003 period snow emissions contribute ∼25% of the total CH2O production, while photochemical oxidation of hydrocarbon appears to be the dominant source. A budget assessment of HOx radicals shows that primary production from O(1D)+H2O and photolysis of snow emitted precursors (i.e., H2O2 and CH2O) are the largest primary HOx sources at Summit, contributing 41% and 40%, respectively. The snow contribution to the HOx budget is mostly in the form of emissions of H2O2. The dominant HOx sink involves the HO2+HO2 reaction forming H2O2, followed by its deposition to snow. These results differ from those previously reported for the South Pole (SP), in that primary production of HOx was shown to be largely driven by both the photolysis of CH2O and H2O2 emissions (46%) with smaller contributions coming from the oxidation of CH4 and the O(1D)+H2O reaction (i.e., 27% each). In sharp contrast to the findings at Summit in 2003, due to the much higher levels of NOx, the SP HOx sinks are dominated by HOx–NOx reactions, leading to the formation and deposition of HNO3 and HO2NO2. Thus, a comparison between SP and Summit studies suggests that snow emissions appear to play a prominent role in controlling primary HOx production in both environments. However, as regards to maintaining highly elevated levels of OH, the two environments differ substantially. At Summit the elevated rate for primary production of HOx is most important; whereas, at SP it is the rapid recycling of the more prevalent HO2 radical, through reaction with NO, back to OH that is primarily responsible.  相似文献   

6.
A parametric model of gas-particle surface reactions is incorporated into a reactive plume model and used to assess the potential importance of heterogeneous surface reactions on gas phase plume chemistry. Heterogeneous loss of the following species is found to be potentially significant: H2O2, PAN, NO3, N2O5, OH, HO2 and in cold weather HO2NO2. This simple model is unable to account for equilibrium/capacity effects in the condensed phase and therefore cannot be used for SO2 and HNO3 reaction with aerosol surface.  相似文献   

7.
Measurements of OH, H2SO4, and MSA at South Pole (SP) Antarctica were recorded as a part of the 2003 Antarctic Chemistry Investigation (ANTCI 2003). The time period 22 November, 2003–2 January, 2004 provided a unique opportunity to observe atmospheric chemistry at SP under both natural conditions as well as those uniquely defined by a solar eclipse event. Results under natural solar conditions generally confirmed those reported previously in the year 2000. In both years the major chemical driver leading to large scale fluctuations in OH was shifts in the concentration levels of NO. Like in 2000, however, the 2003 observational data were systematically lower than model predictions. This can be interpreted as indicating that the model mechanism is still missing a significant HOx sink reaction(s); or, alternatively, that the OH calibration source may have problems. Still a final possibility could involve the integrity of the OH sampling scheme which involved a fixed building site. As expected, during the peak in the solar eclipse both NO and OH showed large decreases in their respective concentrations. Interestingly, the observational OH profile could only be approximated by the model mechanism upon adding an additional HOx radical source in the form of snow emissions of CH2O and/or H2O2. This would lead one to think that either CH2O and/or H2O2 snow emissions represent a significant HOx radical source under summertime conditions at SP. Observations of H2SO4 and MSA revealed both species to be present at very low concentrations (e.g., 5 × 105 and 1 × 105 molec cm?3, respectively), but similar to those reported in 2000. The first measurements of SO2 at SP demonstrated a close coupling with the oxidation product H2SO4. The observed low concentrations of MSA appear to be counter to the most recent thinking by glacio-chemists who have suggested that the plateau's lower atmosphere should have elevated levels of MSA. We speculate here that the absence of MSA may reflect efficient atmospheric removal mechanisms for this species involving either dynamical and/or chemical processes.  相似文献   

8.
Under many conditions size-resolved aqueous-phase chemistry models predict higher sulfate production rates than comparable bulk aqueous-phase models. However, there are special circumstances under which bulk and size-resolved models offer similar predictions. These special conditions include alkaline conditions (when there is a high ammonia to nitric acid ratio or a large amount of alkaline dust) or conditions under which the initial H2O2 concentration exceeds that of SO2. Given that bulk models are less computationally-intensive than corresponding size-resolved models, a model equipped to combine the accuracy of the size-resolved code with the efficiency of the bulk method is proposed in this work. Bulk and two-section size-resolved approaches are combined into a single variable size-resolution model (VSRM) in an effort to combine both accuracy and computational speed. Depending on initial system conditions, bulk or size-resolved calculations are executed based on a set of semi-empirical rules. These rules were generated based on our understanding of the system and from the results of many model simulations for a range of input conditions. For the conditions examined here, on average, the VSRM sulfate predictions are within 3% of a six-section size-resolved model, but the VSRM is fifteen times faster.  相似文献   

9.
Aircraft observations from three recent missions (STRAT, SUCCESS, SONEX) are synthesized into a theoretical analysis of the factors controlling the concentrations of HOx radicals (HOx=OH+peroxy) and the larger reservoir family HOy (HOy=HOx+2H2O2+2CH3OOH+HNO2+HNO4) in the upper troposphere. Photochemical model calculations capture 66% of the variance of observed HOx concentrations. Two master variables are found to determine the variance of the 24 h average HOx concentrations: the primary HOx production rate, P(HOx), and the concentration of nitrogen oxide radicals (NOx=NO+NO2). We use these two variables as a coordinate system to diagnose the photochemistry of the upper troposphere and map the different chemical regimes. Primary HOx production is dominated by the O(1D)+H2O reaction when [H2O]>100 ppmv, and by photolysis of acetone (and possibly other convected HOx precursors) under drier conditions. For the principally northern midlatitude conditions sampled by the aircraft missions, the HOx yield from acetone photolysis ranges from 2 to 3. Methane oxidation amplifies the primary HOx source by a factor of 1.1–1.9. Chemical cycling within the HOx family has a chain length of 2.5–7, while cycling between the HOx family and its HOy reservoirs has a chain length of 1.6–2.2. The number of ozone molecules produced per HOy molecule consumed ranges from 4 to 12, such that ozone production rates vary between 0.3 and 5 ppbv d−1 in the upper troposphere. Three chemical regimes (NOx-limited, transition, NOx-saturated) are identified to describe the dependence of HOx concentrations and ozone production rates on the two master variables P(HOx) and [NOx]. Simplified analytical expressions are derived to express these dependences as power laws for each regime. By applying an eigenlifetime analysis to the HOx–NOx–O3 chemical system, we find that the decay of a perturbation to HOy in the upper troposphere (as from deep convection) is represented by four dominant modes with the longest time scale being factors of 2–3 times longer than the steady-state lifetime of HOy.  相似文献   

10.
Nitrophenols are present in the atmospheric gas phase and in cloud and rainwater. Their formation via aqueous-phase reactions of phenol with the nitronium ion, NO2+, arising from N2O5 and ClNO2 partitioning into the aqueous phase, has been proposed but not verified experimentally. Here, we demonstrate for the first time that gaseous N2O5 and ClNO2 partitioning into dilute aqueous solutions of phenol yields 2- and 4-nitrophenol (and 4-nitrosophenol), but no dinitrophenol isomers. The rate of nitration does not vary significantly between 5 and 20 °C, presumably because of opposing temperature dependences in Henry's law partitioning and reaction rate coefficients. The rate coefficient for reaction of NO2+ with phenol could not be directly quantified but is evidently large enough for this reaction to compete effectively with the reaction between NO2+ and water and to provide a feasible route to nitrophenol production in the atmosphere.  相似文献   

11.
A comparison of a model using five widely known mechanisms (RACM, CB05, LaRC, SAPRC-99, SAPRC-07, and MCMv3.1) has been conducted based on the TexAQS II Radical and Aerosol Measurement Project (TRAMP-2006) field data in 2006. The concentrations of hydroxyl (OH) and hydroperoxy (HO2) radicals were calculated by a zero-dimensional box model with each mechanism and then compared with the OH and HO2 measurements. The OH and HO2 calculated by the model with different mechanisms show similarities and differences with each other and with the measurements. First, measured OH and HO2 are generally greater than modeled for all mechanisms, with the median modeled-to-measured ratios ranging from about 0.8 (CB05) to about 0.6 (SAPRC-99). These differences indicate that either measurement errors, the effects of unmeasured species or chemistry errors in the model or the mechanisms, with some errors being independent of the mechanism used. Second, the modeled and measured ratios of HO2/OH agree when NO is about 1 ppbv, but the modeled ratio is too high when NO was less and too low when NO is more, as seen in previous studies. Third, mechanism–mechanism HOx differences are sensitive to the environmental conditions – in more polluted conditions, the mechanism–mechanism differences are less. This result suggests that, in polluted conditions, the mechanistic details are less important than in cleaner conditions, probably because of the dominance of reactive nitrogen chemistry under polluted conditions.  相似文献   

12.
A physical-chemical model which is an extension of that of Hong and Carmichael (1983) is used to investigate the role of formaldehyde in cloud chemistry. This model takes into account the mass transfer of SO2, O3, NH3, HNO3, H2O2, CO2, HCl, HCHO, O2, OH and HO2 into cloud droplets and their subsequent chemical reactions. The model is used to assess the importance of S(IV)-HCHO adduct formation, the reduction of H2O2 by HCHO, HCHO-free radical interactions, and the formation of HCOOH in the presence of HCHO in cloud droplets.Illustrative calculations indicate that the presence of HCHO inhibits sulfate production rate in cloud droplets. The direct inhibition of sulfate production rate in cloud water due to nucleophilic addition of HSO3 to HCHO(aq) to form hydroxymethanesulfonate (HMSA) is generally low for concentrations of HCHO typical of ambient air. However, inhibition of sulfate production due to formaldehyde-free radical interactions in solution can be important. These formaldehyde-free radical reactions can also generate appreciable quantities of formic acid.  相似文献   

13.
A flow-through chemical reactor model has been exercised to assess the importance of various oxidation reactions and cloud processes on wet removal and redistribution of atmospheric pollutants and to investigate the effect of in-cloud acidification on precipitation chemistry at the surface. Preliminary results indicate that in-cloud acidification accounts for more than 60% of the wet deposited acids derived from acidification of initial SO2, that 42–57% of water-soluble, non-reactive NH3 and HNO3 are removed by wet deposition. The pseudo-first-order conversion rate of SO2 to SO42− ranges from 3 to 25% h −1 depending on initial and boundary conditions.Sensitivity studies have been carried out to test the importance of time evolution of clouds on partitioning of pollutants in the atmosphere and to investigate the variability of precipitation chemistry due to changes in rate constants. The distributions of NH3 and HNO3 are found to be dependent largely on the cloud microphysical parameters, while the distributions of H2O2 and SO2 depend largely on initial conditions of both species. Individual physical and chemical mechanisms can determine the overall rate of sulfate wet deposition at different stages of cloud evolution.  相似文献   

14.
Two different steady-state methods are applied to calculate OH radical concentrations based on the rates of known source and sink processes. The first method, which calculates only OH radical concentrations from measured data including HO2 gives good correlation with measured OH concentrations but overpredicts by 30%. The second method applied calculates OH, HO2 and RO2 radical concentrations simultaneously. This second method overestimates the measured concentrations of OH by almost 3 times. This apparent overprediction may be a result of calculated concentrations of HO2 which appear too high and may be indicative of a gap in our understanding of the relevant peroxy radical chemistry or a result of the limited peroxy radical chemistry assumed by the method.  相似文献   

15.
The purpose of this paper is to study the redistribution of chemical species (OH, HO2, H2O2, HNO3 and H2SO4) over West Africa, where the cloud cover is ubiquitously present, and where deep convection often develops. In this area, because of these cloud systems, chemical species are redistributed by the ascending and descending flow, or leached if they are soluble. So, we carry out a mesoscale study using the Regional Atmospheric Modelling System (RAMS) coupled to a code of gas and aqueous chemistry (RAMS_Chemistry). It takes into account all processes under mesh. We examine several cases following the period (November and July), with inputs emissions (anthropogenic, biogenic and biomass burning). The radicals OH and HO2 are an indicator of possibilities for chemical activity. They characterize the oxidizing power of the atmosphere and are very strong oxidants. The acids HNO3 and H2SO4 are interesting in their transformation into nitrates and sulfates in precipitation. In November, when photochemistry is active during an event of biomass burning, concentrations of chemical species are higher than those of November in the absence of biomass burning. The concentrations of nitric acid double and sulfuric acid increases 70 times. In addition, the concentrations are even lower in July if there is a deep convection. Compared to measures of the African Monsoon Multidisciplinary Analysis (AMMA), the results and observations of radicals OH and HO2 are the same order of magnitude. Emissions from biomass burning increase the concentrations of acid and peroxide, and a deep convection cloud allows the solubility and the washing out of species, reducing their concentration. Rainfalls play a major role in solubility and washing out acids, peroxides and radicals in this region.  相似文献   

16.
Boundary layer concentrations of hydroxyl (OH) and hydroperoxyl (HO2) radicals were measured at 1180 m elevation in a mountainous, forested region of north-western Greece during the AEROsols formation from BIogenic organic Carbon (AEROBIC) field campaign held in July–August 1997. In situ measurements of OH radicals were made by laser-induced fluorescence (LIF) at low pressure, exciting in the (0, 0) band of the A–X system at 308 nm. HO2 radicals were monitored by chemical titration to OH upon the addition of NO, with subsequent detection by LIF. The instrument was calibrated regularly during the field campaign, and demonstrated a sensitivity towards OH and HO2 of 5.2×105 and 2.4×106 molecule cm−3, respectively, for a signal integration period of 2.5 min and a signal-to-noise ratio of 1. Diurnal cycles of OH and HO2 were measured on 10 days within a small clearing of a forest of Greek Fir (Abies Borisi-Regis). In total 4165 OH data points and 1501 HO2 data points were collected at 30 s intervals. Noon-time OH and HO2 concentrations were between 4–12×106 and 0.4–9×108 molecule cm−3, respectively. The performance of the instrument is evaluated, and the data are interpreted in terms of correlations with controlling variables. A significant correlation (r2=0.66) is observed between the OH concentration and the rate of photolysis of ozone, J(O1D). However, OH persisted into the early evening when J(O1D) had fallen to very low values, consistent with the modelling study presented in the following paper (Carslaw et al., 2001, OH and HO2 radical chemistry in a forest region of north-western Greece. Atmospheric Environment 35, 4725–4737) that predicts a significant radical source from the ozonolysis of biogenic alkenes. Normalisation of the OH concentrations for variations in J(O1D) revealed a bell-shaped dependence of OH upon NOx (NO+NO2), which peaked at [NOx] ∼1.75 ppbv. The diurnal variation of HO2 was found to be less correlated with J(O1D) compared to OH.  相似文献   

17.
We present one of the most comprehensive studies of night-time radical chemistry to date, from the Tropospheric ORganic CHemistry experiment (TORCH) in the summer of 2003. TORCH provided a wealth of measurements with which to study the oxidizing capacity of the atmosphere. The measurements provided input to a zero-dimensional box model which has been used to study night-time radical chemistry during the campaign. Average night-time predicted concentrations of OH (2.6 × 105 molecule cm?3), HO2 (2.9 × 107 molecule cm?3) and [HO2+ΣRO2] radicals (2.2 × 108 molecule cm?3) were an order of magnitude smaller than those predicted during the daytime. The model under-predicted the night-time measurements of OH, HO2 and [HO2+ΣRO2] radicals, on average by 41%, 16% and 8% respectively. Whilst the model captured the broad features of night-time radical behaviour, some of the specific features that were observed are hard to explain. A rate of radical production assessment was carried out for the whole campaign between the hours of 00:00 and 04:00. Whilst radical production was limited owing to the absence of photolytic reactions, production routes via the reactions of alkenes with O3 provided an effective night-time radical source. Nitrate radical concentrations were predicted to be 0.6 ppt on average with a peak of 18 ppt on August 9th during a polluted heat wave period. Overall, the nitrate radical contributes about a third of the total initiation via RO2, mostly through reaction with alkenes.  相似文献   

18.
The atmospheric chemical process was simulated using the Carbon Bond 4 (CB-4) model, the aqueous-phase chemistry in Regional Acid Deposition Model and the thermodynamic equilibrium relation of aerosols with the emission inventories of the Emission Database for Global Atmospheric Research, the database of China and South Korea and the Mesoscale Model version 2 (MM5) meteorological fields to examine the spatial distributions of the acidic pollutant concentrations in East Asia for the case of the long-lasting Yellow Sand event in April 1998. The present models simulate quite well the observed general trend and the diurnal variation of concentrations of gaseous pollutants, especially for O3 concentration. However, the model underestimates SO2 and NOx concentration but overestimates O3 concentration largely due to uncertainty in NOx and VOC emissions. It is found that the simulated gaseous pollutants such as SO2, NOx, and NH3 are not transported far away from the source regions but show significant diurnal variations of their concentrations. However, the daily variations of the concentrations are not significant due to invariant emission rates. On the other hand, concentrations of the transformed pollutants including SO42−, NH4+, and NO3 are found to have significant daily variations but little diurnal variations. The model-estimated deposition indicates that dry deposition is largely contributed by gaseous pollutants while wet deposition of pollutants is mainly contributed by the transformed pollutants.  相似文献   

19.
20.
HO2 radical concentrations were measured by a laser-induced fluorescence instrument for three nighttime periods during the intensive field campaign at Rishiri Island, Japan, in June 2000. The HO2 mixing ratio had temporal variations around its average of 4.2±1.2 (1σ) pptv and showed a positive correlation with the summed mixing ratio of four monoterpene species, α-pinene, β-pinene, camphene, and limonene, that sometimes reached 1 ppbv. Our model calculations suggested that ozonolysis reactions of monoterpenes were the main source of nighttime radicals and they explained 58% of measured HO2 concentration levels. The model roughly reproduced the dependence of the HO2 mixing ratio on the square root of the radical production rate due to the ozonolysis reactions of the monoterpenes. However, the absolute HO2 mixing ratio was significantly underpredicted by the model. We discuss possible reasons in terms of misunderstood RO2 chemistry, RO2 interference with HO2 observations, unknown radical production process associated by high NO2 mixing ratio, and the contribution of unmeasured olefinic species to radical production via their reactions with ozone.  相似文献   

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

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