首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Accurate estimation of the soil water balance (SWB) is important for a number of applications (e.g. environmental, meteorological, agronomical and hydrological). The objective of this study was to develop and test techniques for the estimation of soil water fluxes and SWB components (particularly infiltration, evaporation and drainage below the root zone) from soil water records. The work presented here is based on profile soil moisture data measured using dielectric methods, at 30‐min resolution, at an experimental site with different vegetation covers (barley, sunflower and bare soil). Estimates of infiltration were derived by assuming that observed gains in the soil profile water content during rainfall were due to infiltration. Inaccuracies related to diurnal fluctuations present in the dielectric‐based soil water records are resolved by filtering the data with adequate threshold values. Inconsistencies caused by the redistribution of water after rain events were corrected by allowing for a redistribution period before computing water gains. Estimates of evaporation and drainage were derived from water losses above and below the deepest zero flux plane (ZFP), respectively. The evaporation estimates for the sunflower field were compared to evaporation data obtained with an eddy covariance (EC) system located elsewhere in the field. The EC estimate of total evaporation for the growing season was about 25% larger than that derived from the soil water records. This was consistent with differences in crop growth (based on direct measurements of biomass, and field mapping of vegetation using laser altimetry) between the EC footprint and the area of the field used for soil moisture monitoring. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

2.
Understanding how Amazonian rainforests deal with extended droughts is critical in the face of changing climate. This research analyze the physical properties and the soil water dynamics of a deep soil profile in an area of primary forest in central Amazonia to elucidate these processes under drought and nondrought conditions. Physical soil properties derived from soil cores exhibited a distinctive layer between 480 and 880 cm deep, characterized by higher microporosity and low plant water availability. In situ soil moisture measurements collected during the period from January 2003 through February 2006 and for depths ranging from 10 to 1,430 cm suggest that, in the study site, the top 480 cm of the soil profile satisfied most of the transpirational demands in normal climatological years. However, during exceptionally dry periods, such as the 2005 drought, root uptake occurs below 480 cm. As concluded by previous studies, most of the uptake is concentrated in the first meter of the soil profile: More than 40% of the total demand for transpiration is supplied by the top meter of soil. Because deep root uptake occurred at greater depths than normal during the 2005 drought, our results suggest that this is a fundamental mechanism to cope with prolonged droughts.  相似文献   

3.
In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.  相似文献   

4.
A large weighing lysimeter was installed at Yucheng Comprehensive Experimental Station, north China, for evapotranspiration and soil‐water–groundwater exchange studies. Features of the lysimeter include the following: (i) mass resolution equivalent to 0·016 mm of water to accurately and simultaneously determine hourly evapotranspiration, surface evaporation and groundwater recharge; (ii) a surface area of 3·14 m2 and a soil profile depth of 5·0 m to permit normal plant development, soil‐water extraction, soil‐water–groundwater exchanges, and fluctuations of groundwater level; (iii) a special supply–drainage system to simulate field conditions of groundwater within the lysimeter; (iv) a soil mass of about 30 Mg, including both unsaturated and saturated loam. The soil consists mainly of mealy sand and light loam. Monitoring the vegetated lysimeter during the growing period of winter wheat, from October 1998 through to June 1999, indicated that during the period groundwater evaporation contributed 16·6% of total evapotranspiration for a water‐table depth from 1·6 m to 2·4 m below ground surface. Too much irrigation reduced the amount of upward water flow from the groundwater table, and caused deep percolation to the groundwater. Data from neutron probe and tensiometers suggest that soil‐water‐content profiles and soil‐water‐potential profiles were strongly affected by shallow groundwater. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

5.
Abstract

A new approach was developed for estimating vertical soil water fluxes using soil water content time series data. Instead of a traditional fixed time interval, this approach utilizes the time interval between two sequential minima of the soil water storage time series to identify groundwater recharge events and calculate components of the soil water budget. We calculated water budget components: surface-water excess (Sw), infiltration less evapotranspiration (I – ET) and groundwater recharge (R) from May 2001 to January 2003 at eight locations at the USDA Agricultural Research Center, Beltsville, Maryland, USA. High uncertainty was observed for all budget components. This uncertainty was attributed to spatial and temporal variation in Sw, I – ET and R, and was caused by nonuniform rainfall distributions during recharge events, variability in the profile water content, and spatial variability in soil hydraulic properties. The proposed event-based approach allows estimating water budget components when profile water content monitoring data are available.

Citation Guber, A., Gish, T., Pachepsky, Y., McKee, L., Nicholson, T. & Cady, R. (2011) Event-based estimation of water budget components using a network of multi-sensor capacitance probes. Hydrol. Sci. J. 56(7), 1227–1241.  相似文献   

6.
Water content and movement in soil profile and hydrogen isotope composition (δD) of soil water, rainwater, and groundwater were examined in a subalpine dark coniferous forest in the Wolong National Nature Reserve in Sichuan, China, following rainfall events in 2003–2004. Light rainfall increased water content in the litter and at soil depth of 0–80 cm, but the increased soil water was lost in several days. Heavy rainfall increased soil water content up to 85% at depths of 0–40 cm. Following the light rainfall in early spring, the δD of water from the litter, humus, illuvial, and material layers decreased first and then gradually reached the pre‐rainfall level. In summer, light rainfall reached the litter humus, and illuvial layer, but did not hit the material layer. Heavy rainfall affected δD of water in all layers. The δD of soil interflow slightly fluctuated with rainfall events. The δD of shallow groundwater did not differ significantly among all rainfall events. Light rainfall altered the shape of δD profile curve of water in the upper layer of soil, whereas heavy rainfall greatly affected the shape of δD profile curve of water in all soil layers. Following the heavy rainfall, preferential flow initially occurred through macropores, decayed plant roots, and rocks at different depths of soil profile. With continuing rainfall, the litter and surface soil were nearly saturated or fully saturated, and infiltration became homogeneous and plug‐like. Forest soil water, particularly in deeper soil profile, was slightly affected by rainfall and, thus, can be a source of water supply for regional needs, particularly during dry seasons. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

7.
The stable isotopes of hydrogen and oxygen (δ2H and δ18O) are useful conservative tracers for tracking the movement of water in soil. But although the tracking of water infiltrating through the soil profile and its movement as run‐off and groundwater recharge are well developed, water movement through the soil can also include evaporative fractionation. Soil water fractionation factors have, until now, been largely empirical. Unlike open water evaporation where temperature, humidity, and vapour pressure gradient define fractionation, soil water evaporation includes fractionation by soil matrix effects. These effects are still poorly characterized. Here, we present preliminary results from a simple laboratory experiment with four soil admixtures with grain sizes that range from sand to silt and clay. Our results show that soil tension seems to control the isotope fractionation of resident soil water. The relationship between soil tension and equilibrium fractionation appears to be independent of soil texture and appears well supported by thermodynamic theory. Although these results are preliminary, they suggest that future work should go after soil tension effects as a possible explanatory factor of soil water and water vapour fractionation.  相似文献   

8.
The objective of this paper is to simulate the progress of the soil water content distribution in the soil profile with a water table at the bottom of the soil profile during ponding irrigation. This simulation can be done by solving the two‐dimensional Richards's equation for the assimilation of the advancing water jet, which uses the conditions of the two exponential functional forms k = ks eαψ and θ = θr + (θs − θr) eαψ to represent the hydraulic conductivity and volumetric water content, with ψ the pressure as the third variable. We assume that the ground surface becomes ponded and saturated as soon as the water flux passes the dry ground surface. By the technique of transformation, the analytical solution of these two‐dimensional Richards' equations has enabled figures of volumetric water content distribution to be obtained in successive time periods after irrigation. For the example of loam soil, it can simulate the variation of volumetric water content during and after irrigation in the soil profile. The analytical solutions of this paper reflect the real situation simulated, and can be applied to verify those complicated solutions from other analytical models. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

9.
Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ18O and δ2H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.  相似文献   

10.
Intercomparison of soil pore water extraction methods for stable isotope analysis has been a focus of recent studies in relation to plant source waters, which found a wide isotopic variance depending on the extraction method. Few studies have yet explored extraction effects for mobile pore waters that relate to hillslope runoff. This is because it is extremely difficult in natural systems to control the boundary conditions in order to assess and compare impacts of pore water extraction on resulting hillslope flow. With our new semicontrolled experiments on outdoor mini‐hillslopes, we studied mixing and runoff processes by means of stable isotopes of water and quantified relations between pore water extraction methods. We tested the null hypothesis that nondestructive and destructive pore water sampling methods sample the same soil water pool. Three hillslopes were mounted on load cells, filled with loamy sand textured soils from the Landscape Evolution Observatoryat Biosphere 2, equipped with soil moisture and temperature sensors, a bottom outflow, and a surface runoff gauge for isotope sampling. We followed the precipitation isotopic composition over and through the soil profile. One hillslope was instrumented with suction cups, on the second we installed sampling ports for in‐situ soil water vapour measurements, and the third hillslope was sampled destructively for applying the centrifugation and vapour equilibrium methods. All hillslopes were sampled at four depths (0–10, 10–20, 20–30, and 30–40 cm) at three different downslope positions. 2H and 18O analyses were performed via laser spectroscopy. We found no isotopic differences between rainfall, surface runoff, and bottom outflow. The in situ vapour ports' soil isotope data showed the widest spread over all hillslope positions and depths. Centrifugation's and suction cups' isotope results plotted closest to the local meteoric water line and within the range of hillslope runoff and bottom outflow data. Hillslope position did not influence the soil isotope results. These results suggest caution be used in the field when selecting an extraction technique for matching soil waters to runoff waters. Soil suction lysimeters and centrifugation appeared to be the most appropriate tools in this regard.  相似文献   

11.
Data are presented of vertical profiles of soil water potential and soil temperature from a field lysimeter experiment with automatic water-table control at Silwood Park, Ascot, UK. The observed water-table fluxes confirm the analysis of profile response, which shows a diurnal pattern of upward water movement as a result of surface freezing, followed by reversed flow when soil thawing occurs. The results provide a clear, and relatively rare, illustration of the phenomenon of freezing-induced redistribution.  相似文献   

12.
There have been significant recent advances in understanding the ecohydrology of deep soil. However, the links between root development and water usage in the deep critical zone remains poorly understood. To clarify the interaction between water use and root development in deep soil, we investigated soil water and root profiles beyond maximum rooting depth in five apple orchards planted on farmland with stand ages of 8, 11, 15, 18, and 22 years in a subhumid region on the Chinese Loess Plateau. Apple trees rooted progressively deeper for water with increasing stand age and reached 23.2 ± 0.8 m for the 22‐year‐old trees. Soil water deficit in deep soil increased with tree age and was 1,530 ± 43 mm for a stand age of 22 years. Measured root deepening rate was far great than the reported pore water velocity, which demonstrated that trees are mining resident old water. The deficits are not replenished during the life‐span of the orchard, showing a one‐way mining of the critical zone water. The one‐way root water mining may have changed the fine root profile from an exponential pattern in the 8‐year‐old orchard to a relative uniform distribution in older orchards. Our findings enhance our understanding of water‐root interaction in deep soil and reveal the unintended consequences of critical zone dewatering during the lifespan of apple trees.  相似文献   

13.
Despite the potential impact of winter soil water movements in cold regions, relatively few field studies have investigated cold‐season hydrological processes that occur before spring‐onset of snowmelt infiltration. The contribution of soil water fluxes in winter to the annual water balance was evaluated over 5 years of field observations at an agricultural field in Tokachi, Hokkaido, Japan. In two of the winters, soil frost reached a maximum depth of 0·2 m (‘frozen’ winters), whereas soil frost was mostly absent during the remaining three winters (‘unfrozen’ winters). Significant infiltration of winter snowmelt water, to a depth exceeding 1·0 m, occurred during both frozen and unfrozen winters. Such infiltration ranged between 126 and 255 mm, representing 28–51% of total annual soil water fluxes. During frozen winters, a substantial quantity of water (ca 40 mm) was drawn from deeper layers into the 0–0·2 m topsoil layer when this froze. Under such conditions, the progression and regression of the freezing front, regulated by the thickness of snow cover, controlled the quantity of soil water flux below the frozen layer. During unfrozen winters, 13–62 mm of water infiltrated to a depth of 0·2 m, before the spring snowmelt. These results indicate the importance of correctly evaluating winter soil water movement in cold regions. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

14.
Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

15.
Soil water systems have been severely degraded in coal‐mined and subsiding land, where a shallow groundwater table is also present. The present paper discussed the effects of fly ash (FA) and coal gangue (CG) as filling materials on the hydraulic properties and water content distribution in a profile for the purpose of rehabilitating subsided lands. The saturated water content, water characteristic curve, and water diffusivity of local soil, FA, CG, and a mixture of FA and CG (“mixed filling”) were characterized. A column experiment was conducted to investigate the changes of water content in profiles reconstructed from the combination of soil and filling materials, including soil only, FA, CG, and a mixture of FA and CG, which were used to fill the lower part of the reconstructed profile. The mixture of FA and CG was found to possess similar hydraulic properties to those of the soil, particularly high water‐holding capacity and permeability. Moreover, the volumetric water contents in the whole profile containing the mixture of FA and CG were consistent with those of the profile reconstructed with soil only. As a result, it is recommended to adopt the mixture of FA and CG for reconstructing the lower profile of the land to alleviate or rehabilitate subsided land in coal mines.  相似文献   

16.
Spatial distribution of soil macroporosity was determined for a forest podzol from tension infiltrometer measurements at the soil surface. Surface‐derived macroporosity values were compared with point infiltration characteristics obtained from soil water content and soil water chemistry measurements during an experimental irrigation, and with parameters of a kinematic wave model applied to soil water content data. Macroporosity estimated by the tension infiltrometer ranged from 0·00087 to 0·0219% of soil volume, and infiltration at these two sites was dominated by propagation of a well‐defined wetting front through the soil profile and bypass flow via soil macropores, respectively. Infiltration at sites with intermediate macroporosities reflected a combination of these two processes, although results were inconclusive at one site owing to lateral flow at the base of the soil profile. There was no agreement between macroporosities estimated by the tension infiltrometer and the kinematic wave model. The maximum soil conductance parameter within the profile at a site, however, was related directly to the surface‐derived macroporosity. The partial agreement between surface‐derived macroporosity estimates and point infiltration characteristics shown here supports the use of tension infiltrometry as a rapid, non‐destructive method of assessing spatial variations in the relative contribution of macropore flow to the infiltration process. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

17.
Plants have been shown to affect soil water content and temperature. Previous studies were conducted mainly in forestry and agricultural soils, where conditions of soil and vegetation are different from those in an urban landscape. In an urban landscape, the influence of plant roots on electrical conductivity, soil water content and temperature is still not clear. This study aims to investigate the effects of soil water content and temperature on electrical conductivity in vegetated soils through an integrated field monitoring and computational modelling approach. A new relationship between soil electrical conductivity and water content as well as temperature is proposed. Field monitoring was conducted in both vegetated (tree species) and bare soils. The monitoring included measurements of soil water content, soil temperature and soil electrical conductivity. This was followed by response surface regression modelling. Measured results show that soil temperature at shallow depths was lower in vegetated soil than that in the bare soil. This observation was also consistent with the higher soil water content and hence, higher electrical conductivity under tree canopy. The model developed could predict nonlinear relationships between electrical conductivity and soil temperature and water content. Uncertainty analysis indicated normal distribution for electrical conductivity under variation of soil temperature and water content. © 2018 John Wiley & Sons, Ltd.  相似文献   

18.
Spatio‐temporal heterogeneity in soil water content is recognized as a common phenomenon, but heterogeneity in the hydrogen and oxygen isotope composition of soil water, which can reveal processes of water cycling within soils, has not been well studied. New advances are being driven by measurement approaches allowing sampling with high density in both space and time. Using in situ soil water vapour probe techniques, combined with conventional soil and plant water vacuum distillation extraction, we monitored the hydrogen and oxygen stable isotopic composition of soil and plant waters at paired sites dominated by grasses and Gambel's oak (Quercus gambelii) within a semiarid montane ecosystem over the course of a growing season. We found that sites spaced only 20 m apart had profoundly different soil water isotopic and volumetric conditions. We document patterns of depth‐ and time‐explicit variation in soil water isotopic conditions at these sites and consider mechanisms for the observed heterogeneity. We found that soil water content and isotopic variability were damped under Q. gambelii, perhaps due in part to hydraulic redistribution of deep soil water or groundwater by Q. gambelii in these soils relative to the grass‐dominated site. We also found some support for H isotope discrimination effects during water uptake by Q. gambelii. In this ecosystem, the soil water content was higher than that at the neighbouring Grass site, and thus, 25% more water was available for transpiration by Q. gambelii compared with the Grass site. This work highlights the role of plants in governing soil water variation and demonstrates that they can also strongly influence the isotope ratios of soil water. The resulting fine‐scale heterogeneity has implications for the use of isotope tracers to study soil hydrology and evaporation and transpiration fluxes to improve understanding of water cycling through the soil–plant–atmosphere continuum.  相似文献   

19.
Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one‐dimensional soil‐water‐balance model to estimate precipitation‐ and irrigation‐generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf‐area index and is limited by soil‐moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model‐calculated water contents of 11 soil‐depth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998–2001). Each 50‐m2 site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean‐squared error between measured and model‐calculated water content of the top 180 cm was 4·2 cm, or 9·3% of average total water content. In addition, model‐calculated evapotranspiration compared well with that measured by a large‐scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water‐table declines. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

20.
Accurate estimates of seasonal evapotranspiration (ET) at different temporal and spatial scales are essential for understanding the biological and environmental determinants of ecosystem water balance in arid regions and the patterns of water utilization by the vegetation. For this purpose, remote sensing ET estimates of a Patagonian desert in Southern Argentina were verified with field measurements of soil evaporation and plant transpiration using an open top chamber. Root distribution and seasonal variation in soil volumetric water content were also analysed. There was a high correlation between remote sensing and field measurements of ecosystem water fluxes. A substantial amount of the annual ET occurred in spring and early summer (73.4 mm) using winter rain stored in the soil profile and resulting in water content depletion of the upper soil layers. A smaller amount of annual ET was derived from few rainfall events occurring during the mid or late summer (41.4 mm). According to remote sensing, the 92.9% of the mean annual precipitation returns to the atmosphere by transpiration or evaporation from the bare soil and by canopy interception. Only 7.1% infiltrates to soil layers deeper than 200 cm contributing to the water table recharge. Fourier time series analysis, cross‐correlation methods and multiple linear regression models were used to analyse 11 years of remote sensing data to assess determinants of water fluxes. A linear model predicts well the variables that drive complex ecosystem processes such as ET. Leaf area index and air temperature were not linearly correlated to ET because of the multiple interaction among variables resulting in time lags with ET variations and thus these two variables were not included in the linear model. Soil water content, the fraction of photosynthetic active radiation and precipitation explained 86% of the ET monthly variations. The high volumetric water content and the small seasonal variations at 200‐cm depth were probably the result of little water uptake from deeper soil horizons by roots with low hydraulic conductivity. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

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

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