首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 468 毫秒
1.
De Li Liu  Heping Zuo 《Climatic change》2012,115(3-4):629-666
This paper outlines a new statistical downscaling method based on a stochastic weather generator. The monthly climate projections from global climate models (GCMs) are first downscaled to specific sites using an inverse distance-weighted interpolation method. A bias correction procedure is then applied to the monthly GCM values of each site. Daily climate projections for the site are generated by using a stochastic weather generator, WGEN. For downscaling WGEN parameters, historical climate data from 1889 to 2008 are sorted, in an ascending order, into 6 climate groups. The WGEN parameters are downscaled based on the linear and non-linear relationships derived from the 6 groups of historical climates and future GCM projections. The overall averaged confidence intervals for these significant linear relationships between parameters and climate variables are 0.08 and 0.11 (the range of these parameters are up to a value of 1.0) at the observed mean and maximum values of climate variables, revealing a high confidence in extrapolating parameters for downscaling future climate. An evaluation procedure is set up to ensure that the downscaled daily sequences are consistent with monthly GCM output in terms of monthly means or totals. The performance of this model is evaluated through the comparison between the distributions of measured and downscaled climate data. Kruskall-Wallis rank (K-W) and Siegel-Tukey rank sum dispersion (S-T) tests are used. The results show that the method can reproduce the climate statistics at annual, monthly and daily time scales for both training and validation periods. The method is applied to 1062 sites across New South Wales (NSW) for 9 GCMs and three IPCC SRES emission scenarios, B1, A1B and A2, for the period of 1900–2099. Projected climate changes by 7 GCMs are also analyzed for the A2 emission scenario based on the downscaling results.  相似文献   

2.
Joint variable spatial downscaling   总被引：1，自引：0，他引：1  
Joint Variable Spatial Downscaling (JVSD), a new statistical technique for downscaling gridded climatic variables, is developed to generate high resolution gridded datasets for regional watershed modeling and assessments. The proposed approach differs from previous statistical downscaling methods in that multiple climatic variables are downscaled simultaneously and consistently to produce realistic climate projections. In the bias correction step, JVSD uses a differencing process to create stationary joint cumulative frequency statistics of the variables being downscaled. The functional relationship between these statistics and those of the historical observation period is subsequently used to remove GCM bias. The original variables are recovered through summation of bias corrected differenced sequences. In the spatial disaggregation step, JVSD uses a historical analogue approach, with historical analogues identified simultaneously for all atmospheric fields and over all areas of the basin under study. Analysis and comparisons are performed for 20th Century Climate in Coupled Models (20C3M), broadly available for most GCMs. The results show that the proposed downscaling method is able to reproduce the sub-grid climatic features as well as their temporal/spatial variability in the historical periods. Comparisons are also performed for precipitation and temperature with other statistical and dynamic downscaling methods over the southeastern US and show that JVSD performs favorably. The downscaled sequences are used to assess the implications of GCM scenarios for the Apalachicola-Chattahoochee-Flint river basin as part of a comprehensive climate change impact assessment.  相似文献   

3.
A prerequisite of a successful statistical downscaling is that large-scale predictors simulated by the General Circulation Model (GCM) must be realistic. It is assumed here that features smaller than the GCM resolution are important in determining the realism of the large-scale predictors. It is tested whether a three-step method can improve conventional one-step statistical downscaling. The method uses predictors that are upscaled from a dynamical downscaling instead of predictors taken directly from a GCM simulation. The method is applied to downscaling of monthly precipitation in Sweden. The statistical model used is a multiple regression model that uses indices of large-scale atmospheric circulation and 850-hPa specific humidity as predictors. Data from two GCMs (HadCM2 and ECHAM4) and two RCM experiments of the Rossby Centre model (RCA1) driven by the GCMs are used. It is found that upscaled RCA1 predictors capture the seasonal cycle better than those from the GCMs, and hence increase the reliability of the downscaled precipitation. However, there are only slight improvements in the simulation of the seasonal cycle of downscaled precipitation. Due to the cost of the method and the limited improvements in the downscaling results, the three-step method is not justified to replace the one-step method for downscaling of Swedish precipitation.  相似文献   

4.
Hydrological modeling for climate-change impact assessment implies using meteorological variables simulated by global climate models (GCMs). Due to mismatching scales, coarse-resolution GCM output cannot be used directly for hydrological impact studies but rather needs to be downscaled. In this study, we investigated the variability of seasonal streamflow and flood-peak projections caused by the use of three statistical approaches to downscale precipitation from two GCMs for a meso-scale catchment in southeastern Sweden: (1) an analog method (AM), (2) a multi-objective fuzzy-rule-based classification (MOFRBC) and (3) the Statistical DownScaling Model (SDSM). The obtained higher-resolution precipitation values were then used to simulate daily streamflow for a control period (1961–1990) and for two future emission scenarios (2071–2100) with the precipitation-streamflow model HBV. The choice of downscaled precipitation time series had a major impact on the streamflow simulations, which was directly related to the ability of the downscaling approaches to reproduce observed precipitation. Although SDSM was considered to be most suitable for downscaling precipitation in the studied river basin, we highlighted the importance of an ensemble approach. The climate and streamflow change signals indicated that the current flow regime with a snowmelt-driven spring flood in April will likely change to a flow regime that is rather dominated by large winter streamflows. Spring flood events are expected to decrease considerably and occur earlier, whereas autumn flood peaks are projected to increase slightly. The simulations demonstrated that projections of future streamflow regimes are highly variable and can even partly point towards different directions.  相似文献   

5.
Many impact studies require climate change information at a finer resolution than that provided by global climate models (GCMs). This paper investigates the performances of existing state-of-the-art rule induction and tree algorithms, namely single conjunctive rule learner, decision table, M5 model tree, and REPTree, and explores the impact of climate change on maximum and minimum temperatures (i.e., predictands) of 14 meteorological stations in the Upper Thames River Basin, Ontario, Canada. The data used for evaluation were large-scale predictor variables, extracted from National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis dataset and the simulations from third generation Canadian coupled global climate model. Data for four grid points covering the study region were used for developing the downscaling model. M5 model tree algorithm was found to yield better performance among all other learning techniques explored in the present study. Hence, this technique was applied to project predictands generated from GCM using three scenarios (A1B, A2, and B1) for the periods (2046–2065 and 2081–2100). A simple multiplicative shift was used for correcting predictand values. The potential of the downscaling models in simulating predictands was evaluated, and downscaling results reveal that the proposed downscaling model can reproduce local daily predictands from large-scale weather variables. Trend of projected maximum and minimum temperatures was studied for historical as well as downscaled values using GCM and scenario uncertainty. There is likely an increasing trend for T max and T min for A1B, A2, and B1 scenarios while decreasing trend has been observed for B1 scenarios during 2081–2100.  相似文献   

6.
Summary Uncertainty analysis is used to make a quantitative evaluation of the reliability of statistically downscaled climate data representing local climate conditions in the northern coastlines of Canada. In this region, most global climate models (GCMs) have inherent weaknesses to adequately simulate the climate regime due to difficulty in resolving strong land/sea discontinuities or heterogeneous land cover. The performance of the multiple regression-based statistical downscaling model in reproducing the observed daily minimum/maximum temperature, and precipitation for a reference period (1961–1990) is evaluated using climate predictors derived from NCEP reanalysis data and those simulated by two coupled GCMs (the Canadian CGCM2 and the British HadCM3). The Wilcoxon Signed Rank test and bootstrap confidence-interval estimation techniques are used to perform uncertainty analysis on the downscaled meteorological variables. The results show that the NCEP-driven downscaling results mostly reproduced the mean and variability of the observed climate very well. Temperatures are satisfactorily downscaled from HadCM3 predictors while some of the temperatures downscaled from CGCM2 predictors are statistically significantly different from the observed. The uncertainty in precipitation downscaled with CGCM2 predictors is comparable to the ones downscaled from HadCM3. In general, all downscaling results reveal that the regression-based statistical downscaling method driven by accurate GCM predictors is able to reproduce the climate regime over these highly heterogeneous coastline areas of northern Canada. The study also shows the applicability of uncertainty analysis techniques in evaluating the reliability of the downscaled data for climate scenarios development. Authors’ addresses: Dr. Yonas B. Dibike, NSERC Research Fellow, OURANOS Consortium, 550 Sherbrooke Street West, 19th Floor, Montreal (QC) H3A 1B9, Canada; Philippe Gachon, Adaptation and Impact Research Division (AIRD), Atmospheric Science and Technology Directorate, Environment Canada at Ouranos, Montreal (QC), Canada; André St-Hilaire and Taha B. M. J. Ouarda, Institut National de la Recherche Scientifique Centre Eau, Terre & Environnement (INRS-ETE), University of Québec, 490 Rue de La Couronne, Québec (QC) G1K 9A9, Canada; Van T.-V. Nguyen, Department of Civil Engineering and Applied Mechanics, McGill University, 817 Sherbrooke Street West, Montreal (QC) H3A 2K6, Canada.  相似文献   

7.
Statistical downscaling is a technique widely used to overcome the spatial resolution problem of General Circulation Models (GCMs). Nevertheless, the evaluation of uncertainties linked with downscaled temperature and precipitation variables is essential to climate impact studies. This paper shows the potential of a statistical downscaling technique (in this case SDSM) using predictors from three different GCMs (GCGM3, GFDL and MRI) over a highly heterogeneous area in the central Andes. Biases in median and variance are estimated for downscaled temperature and precipitation using robust statistical tests, respectively Mann?CWhitney and Brown?CForsythe's tests. In addition, the ability of the downscaled variables to reproduce extreme events is tested using a frequency analysis. Results show that uncertainties in downscaled precipitations are high and that simulated precipitation variables failed to reproduce extreme events accurately. Nevertheless, a greater confidence remains in downscaled temperatures variables for the area. GCMs performed differently for temperature and precipitation as well as for the different test. In general, this study shows that statistical downscaling is able to simulate with accuracy temperature variables. More inhomogeneities are detected for precipitation variables. This first attempt to test uncertainties of statistical downscaling techniques in the heterogeneous arid central Andes contributes therefore to an improvement of the quality of predictions of climate impact studies in this area.  相似文献   

8.
The potential impact of climate warming on patterns of malaria transmission has been the subject of keen scientific and policy debate. Standard climate models (GCMs) characterize climate change at relatively coarse spatial and temporal scales. However, malaria parasites and the mosquito vectors respond to diurnal variations in conditions at very local scales. Here we bridge this gap by downscaling a series of GCMs to provide high-resolution temperature data for four different sites and show that although outputs from both the GCM and the downscaled models predict diverse but qualitatively similar effects of warming on the potential for adult mosquitoes to transmit malaria, the predicted magnitude of change differs markedly between the different model approaches. Raw GCM model outputs underestimate the effects of climate warming at both hot (3-fold) and cold (8–12 fold) extremes, and overestimate (3-fold) the change under intermediate conditions. Thus, downscaling could add important insights to the standard application of coarse-scale GCMs for biophysical processes driven strongly by local microclimatic conditions.  相似文献   

9.
Regression-based statistical downscaling model (SDSM) is an appropriate method which broadly uses to resolve the coarse spatial resolution of general circulation models (GCMs). Nevertheless, the assessment of uncertainty propagation linked with climatic variables is essential to any climate change impact study. This study presents a procedure to characterize uncertainty analysis of two GCM models link with Long Ashton Research Station Weather Generator (LARS-WG) and SDSM in one of the most vulnerable international wetland, namely “Shadegan” in an arid region of Southwest Iran. In the case of daily temperature, uncertainty is estimated by comparing monthly mean and variance of downscaled and observed daily data at a 95 % confidence level. Uncertainties were then evaluated from comparing monthly mean dry and wet spell lengths and their 95 % CI in daily precipitation downscaling using 1987–2005 interval. The uncertainty results indicated that the LARS-WG is the most proficient model at reproducing various statistical characteristics of observed data at a 95 % uncertainty bounds while the SDSM model is the least capable in this respect. The results indicated a sequences uncertainty analysis at three different climate stations and produce significantly different climate change responses at 95 % CI. Finally the range of plausible climate change projections suggested a need for the decision makers to augment their long-term wetland management plans to reduce its vulnerability to climate change impacts.  相似文献   

10.
Future climate projections from general circulation models (GCMs) predict an acceleration of the global hydrological cycle throughout the 21st century in response to human-induced rise in temperatures. However, projections of GCMs are too coarse in resolution to be used in local studies of climate change impacts. To cope with this problem, downscaling methods have been developed that transform climate projections into high resolution datasets to drive impact models such as rainfall-runoff models. Generally, the range of changes simulated by different GCMs is considered to be the major source of variability in the results of such studies. However, the cascade of uncertainty in runoff projections is further elongated by differences between impact models, especially where robust calibration is hampered by the scarcity of data. Here, we address the relative importance of these different sources of uncertainty in a poorly monitored headwater catchment of the Ecuadorian Andes. Therefore, we force 7 hydrological models with downscaled outputs of 8 GCMs driven by the A1B and A2 emission scenarios over the 21st century. Results indicate a likely increase in annual runoff by 2100 with a large variability between the different combinations of a climate model with a hydrological model. Differences between GCM projections introduce a gradually increasing relative uncertainty throughout the 21st century. Meanwhile, structural differences between applied hydrological models still contribute to a third of the total uncertainty in late 21st century runoff projections and differences between the two emission scenarios are marginal.  相似文献   

11.
  相似文献   

12.
Physical scaling (SP) method downscales climate model data to local or regional scales taking into consideration physical characteristics of the area under analysis. In this study, multiple SP method based models are tested for their effectiveness towards downscaling North American regional reanalysis (NARR) daily precipitation data. Model performance is compared with two state-of-the-art downscaling methods: statistical downscaling model (SDSM) and generalized linear modeling (GLM). The downscaled precipitation is evaluated with reference to recorded precipitation at 57 gauging stations located within the study region. The spatial and temporal robustness of the downscaling methods is evaluated using seven precipitation based indices. Results indicate that SP method-based models perform best in downscaling precipitation followed by GLM, followed by the SDSM model. Best performing models are thereafter used to downscale future precipitations made by three global circulation models (GCMs) following two emission scenarios: representative concentration pathway (RCP) 2.6 and RCP 8.5 over the twenty-first century. The downscaled future precipitation projections indicate an increase in mean and maximum precipitation intensity as well as a decrease in the total number of dry days. Further an increase in the frequency of short (1-day), moderately long (2–4 day), and long (more than 5-day) precipitation events is projected.  相似文献   

13.
There is increasing concern that avoiding climate change impacts will require proactive adaptation, particularly for infrastructure systems with long lifespans. However, one challenge in adaptation is the uncertainty surrounding climate change projections generated by general circulation models (GCMs). This uncertainty has been addressed in different ways. For example, some researchers use ensembles of GCMs to generate probabilistic climate change projections, but these projections can be highly sensitive to assumptions about model independence and weighting schemes. Because of these issues, others argue that robustness-based approaches to climate adaptation are more appropriate, since they do not rely on a precise probabilistic representation of uncertainty. In this research, we present a new approach for characterizing climate change risks that leverages robust decision frameworks and probabilistic GCM ensembles. The scenario discovery process is used to search across a multi-dimensional space and identify climate scenarios most associated with system failure, and a Bayesian statistical model informed by GCM projections is then developed to estimate the probability of those scenarios. This provides an important advancement in that it can incorporate decision-relevant climate variables beyond mean temperature and precipitation and account for uncertainty in probabilistic estimates in a straightforward way. We also suggest several advancements building on prior approaches to Bayesian modeling of climate change projections to make them more broadly applicable. We demonstrate the methodology using proposed water resources infrastructure in Lake Tana, Ethiopia, where GCM disagreement on changes in future rainfall presents a major challenge for infrastructure planning.  相似文献   

14.
Projected changes to the global climate system have great implications for the incidence of large infrequent fires in many regions. Here we examine the synoptic-scale and local-scale influences on the incidence of extreme fire weather days and consider projections of the large-scale mean climate to explore future fire weather projections. We focus on a case study region with periodic extreme fire dangers; southeast Tasmania, Australia. We compare the performance of a dynamically downscaled regional climate model with Global Climate Model outputs as a tool for examining the local-scale influences while accounting for high regional variability. Many of the worst fires in Tasmania and the southeast Australian region are associated with deep cold fronts and strong prefrontal winds. The downscaled simulations reproduce this synoptic type with greater fidelity than a typical global climate model. The incidence of systems in this category is projected to increase through the century under a high emission scenario, driven mainly by an increase in the temperature of air masses, with little change in the strength of the systems. The regional climate model projected increase in frequency is smaller than for the global climate models used as input, with a large model range and natural variability. We also demonstrate how a blocking Foehn effect and topographic channelling contributed to the extreme conditions during an extreme fire weather day in Tasmania in January 2013. Effects such as these are likely to contribute to high fire danger throughout the century. Regional climate models are useful tools that enable various meteorological drivers of fire danger to be considered in projections of future fire danger.  相似文献   

15.
Summary A regression-based methodology was used to downscale hourly and daily station-scale meteorological variables from outputs of large-scale general circulation models (GCMs). Meteorological variables include air temperature, dew point, and west–east and south–north wind velocities at the surface and three upper atmospheric levels (925, 850, and 500 hPa), as well as mean sea-level air pressure and total cloud cover. Different regression methods were used to construct downscaling transfer functions for different weather variables. Multiple stepwise regression analysis was used for all weather variables, except total cloud cover. Cumulative logit regression was employed for analysis of cloud cover, since cloud cover is an ordered categorical data format. For both regression procedures, to avoid multicollinearity between explanatory variables, principal components analysis was used to convert inter-correlated weather variables into uncorrelated principal components that were used as predictors. The results demonstrated that the downscaling method was able to capture the relationship between the premises and the response; for example, most hourly downscaling transfer functions could explain over 95% of the total variance for several variables (e.g. surface air temperature, dew point, and air pressure). Downscaling transfer functions were validated using a cross-validation scheme, and it was concluded that the functions for all weather variables used in the study are reliable. Performance of the downscaling method was also evaluated by comparing data distributions and extreme weather characteristics of downscaled GCM historical runs and observations during the period 1961–2000. The results showed that data distributions of downscaled GCM historical runs for all weather variables are significantly similar to those of observations. In addition, extreme characteristics of the downscaled meteorological variables (e.g. temperature, dew point, air pressure, and total cloud cover) were examined. Authors’ addresses: Chad Shouquan Cheng, Guilong Li, Qian Li, Atmospheric Science and Applications Unit, Meteorological Service of Canada Branch-Ontario, Environment Canada, 4905 Dufferin Street, Toronto, Ontario, Canada M3H 5T4; Heather Auld, Adaptation and Impacts Research Division, MSC Branch, Environment Canada, Toronto, Canada.  相似文献   

16.
This study extends a stochastic downscaling methodology to generation of an ensemble of hourly time series of meteorological variables that express possible future climate conditions at a point-scale. The stochastic downscaling uses general circulation model (GCM) realizations and an hourly weather generator, the Advanced WEather GENerator (AWE-GEN). Marginal distributions of factors of change are computed for several climate statistics using a Bayesian methodology that can weight GCM realizations based on the model relative performance with respect to a historical climate and a degree of disagreement in projecting future conditions. A Monte Carlo technique is used to sample the factors of change from their respective marginal distributions. As a comparison with traditional approaches, factors of change are also estimated by averaging GCM realizations. With either approach, the derived factors of change are applied to the climate statistics inferred from historical observations to re-evaluate parameters of the weather generator. The re-parameterized generator yields hourly time series of meteorological variables that can be considered to be representative of future climate conditions. In this study, the time series are generated in an ensemble mode to fully reflect the uncertainty of GCM projections, climate stochasticity, as well as uncertainties of the downscaling procedure. Applications of the methodology in reproducing future climate conditions for the periods of 2000–2009, 2046–2065 and 2081–2100, using the period of 1962–1992 as the historical baseline are discussed for the location of Firenze (Italy). The inferences of the methodology for the period of 2000–2009 are tested against observations to assess reliability of the stochastic downscaling procedure in reproducing statistics of meteorological variables at different time scales.  相似文献   

17.
This study presents a method to incorporate uncertainty of climate variables in climate change impact assessments, where the uncertainty being considered refers to the divergence of general circulation model (GCM) projections. The framework assesses how much bias occurs when the uncertainties of climate variables are ignored. The proposed method is based on the second-order expansion of Taylor series, called second-order approximation (SOA). SOA addresses the bias which occurs by assuming the expected value of a function is equal to the function of the expected value of the predictors. This assumption is not valid for nonlinear systems, such as in the case of the relationship of climate variables to streamflow. To investigate the value of SOA in the climate change context, statistical downscaling models for monthly streamflow were set up for six hydrologic reference stations in Australia which cover contrasting hydro-climate regions. It is shown that in all locations SOA makes the largest difference for low flows and changes the overall mean flow by 1–3%. Another advantage of the SOA approach is that the individual contribution of each climate variable to the total difference can be estimated. It is found that geopotential height and specific humidity cause more bias than wind speeds in the downscaling models considered here.  相似文献   

18.
统计降尺度法对华北地区未来区域气温变化情景的预估   总被引：31，自引：1，他引：31  
迄今为止，大部分海气耦合气候模式（AOGCM）的空间分辨率还较低，很难对区域尺度的气候变化情景做合理的预测。降尺度法已广泛用于弥补AOGCM在这方面的不足。作者采用统计降尺度方法对1月和7月华北地区49个气象观测站的未来月平均温度变化情景进行预估。采用的统计降尺度方法是主分量分析与逐步回归分析相结合的多元线性回归模型。首先，采用1961～2000年的 NCEP再分析资料和49个台站的观测资料建立月平均温度的统计降尺度模型，然后把建立的统计降尺度模型应用于HadCM3 SRES A2 和 B2 两种排放情景， 从而生成各个台站1950～2099年1月份和7月份温度变化情景。结果表明:在当前气候条件下，无论1月还是7月，统计降尺度方法模拟的温度与观测的温度有很好的一致性，而且在大多数台站，统计降尺度模拟气温与观测值相比略微偏低。对于未来气候情景的预估方面，无论1月还是7月，也无论是HadCM3 SRES A2 还是B2排放情景驱动统计模型，结果表明大多数的站点都存在温度的明显上升趋势，同时7月的上升趋势与1月相比偏低。  相似文献   

19.
Dynamical downscaling has been recognized as a useful tool not only for the climate community, but also for associated application communities such as the environmental and hydrological societies. Although climate projection data are available in lower-resolution general circulation models (GCMs), higher-resolution climate projections using regional climate models (RCMs) have been obtained over various regions of the globe. Various model outputs from RCMs with a high resolution of even as high as a few km have become available with heavy weight on applications. However, from a scientific point of view in numerical atmospheric modeling, it is not clear how to objectively judge the degree of added value in the RCM output against the corresponding GCM results. A key factor responsible for skepticism is based on the fundamental limitations in the nesting approach between GCMs and RCMs. In this article, we review the current status of the dynamical downscaling for climate prediction, focusing on basic assumptions that are scrutinized from a numerical weather prediction (NWP) point of view. Uncertainties in downscaling due to the inconsistencies in the physics packages between GCMs and RCMs were revealed. Recommendations on how to tackle the ultimate goal of dynamical downscaling were also described.  相似文献   

20.
The use and development of bias correction (BC) methods has grown fast in recent years, due to the increased demand of unbiased projections by many sectoral climate change impact applications. Case studies are frequently based on multi-variate climate indices (CIs) combining two or more essential climate variables that are frequently individually corrected prior to CI calculation. This poses the question of whether the BC method modifies the inter-variable dependencies and eventually the climate change signal. The direct bias correction of the multi-variate CI stands as a usual alternative, since it preserves the physical and temporal coherence among the primary variables as represented in the dynamical model output, at the expense of incorporating the individual biases on the CI with an effect difficult to foresee, particularly in the case of complex CIs bearing in their formulation non-linear relationships between components. Such is the case of the Fire Weather Index (FWI), a meteorological fire danger indicator frequently used in forest fire prevention and research. In the present work, we test the suitability of the direct BC approach on FWI as a representative multi-variate CI, assessing its performance in present climate conditions and its effect on the climate change signal when applied to future projections. Moreover, the results are compared with the common approach of correcting the input variables separately. To this aim, we apply the widely used empirical quantile mapping method (QM), adjusting the 99 empirical percentiles. The analysis of the percentile adjustment function (PAF) provides insight into the effect of the QM on the climate change signal. Although both approaches present similar results in the present climate, the direct correction introduces a greater modification of the original change signal. These results warn against the blind use of QM, even in the case of essential climate variables or uni-variate CIs.  相似文献   

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

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