首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
副热带东北太平洋混合层深度及其对潜沉的影响   总被引:1,自引:0,他引:1  
The present climate simulations of the mixed layer depth(MLD) and the subduction rate in the subtropical Northeast Pacific are investigated based on nine of the CMIP5 models. Compared with the observation data,spatial patterns of the MLD and the subduction rate are well simulated in these models. The spatial pattern of the MLD is nonuniform, with a local maximum MLD(140 m) region centered at(28°N, 135°W) in late winter. The nonuniform MLD pattern causes a strong MLD front on the south of the MLD maximum region, controls the lateral induction rate pattern, and then decides the nonuniform distribution of the subduction rate. Due to the inter-regional difference of the MLD, we divide this area into two regions. The relatively uniform Ekman pumping has little effect on the nonuniform subduction spatial pattern, though it is nearly equal to the lateral induction in values. In the south region, the northward warm Ekman advection(–1.75×10~(–7) K/s) controls the ocean horizontal temperature advection(–0.85×10~(–7) K/s), and prevents the deepening of the MLD. In the ensemble mean, the contribution of the ocean advection to the MLD is about –29.0 m/month, offsetting the sea surface net heat flux contribution(33.9 m/month). While in the north region, the southward cold advection deepens the MLD(21.4 m/month) as similar as the heat flux(30.4 m/month). In conclusion, the nonuniform MLD pattern is dominated by the nonuniform ocean horizontal temperature advection. This new finding indicates that the upper ocean current play an important role in the variability of the winter MLD and the subduction rate.  相似文献   

2.
In this study, we develop a variable-grid global ocean general circulation model(OGCM) with a fine grid(1/6)°covering the area from 20°S–50°N and from 99°–150°E, and use the model to investigate the isopycnal surface circulation in the South China Sea(SCS). The simulated results show four layer structures in vertical: the surface and subsurface circulation of the SCS are characterized by the monsoon driven circulation, with basin-scaled cyclonic gyre in winter and anti-cyclonic gyre in summer. The intermediate layer circulation is opposite to the upper layer, showing anti-cyclonic gyre in winter but cyclonic gyre in summer. The circulation in the deep layer is much weaker in spring and summer, with the maximum velocity speed below 0.6 cm/s. In fall and winter, the SCS deep layer circulation shows strong east boundary current along the west coast of Philippine with the velocity speed at 1.5 m/s, which flows southward in fall and northward in winter. The results have also revealed a fourlayer vertical structure of water exchange through the Luzon Strait. The dynamics of the intermediate and deep circulation are attributed to the monsoon driving and the Luzon Strait transport forcing.  相似文献   

3.
An analysis of historical oxygen data provides evidence on the water exchange between theSouth China Sea (SCS) and the Pacific Ocean (PO). In the vicinity of the Luzon Strait (LS) , the dissolved oxygen concentration of sea water is found to be lower on the Pacific side than on the SCS side at depths between 700 and 1500 m (intermediate layer) , while the situation is reversed above 700 m (upper layer) and below 1 500 m (deep layer). The evidence suggests that water exits the SCS in the intermediate layer but enters it from the Pacific in both the upper and the deep layers, supporting the earlier speculation that the Luzon Strait transport has a sandwiched structure in the vertical. Within the SCS basin, the oxygen distribution indicates widespread vertical movement, including the upwelling in the intermediate layer and the downwelling in the deep layer.  相似文献   

4.
冬季婆罗洲岛西北沿岸上升流的时空特征及机理研究   总被引:1,自引:0,他引:1  
Winter coastal upwelling off northwest Borneo in the South China Sea(SCS) is investigated by using satellite data, climatological temperature and salinity fields and reanalysis data. The upwelling forms in December, matures in January, starts to decay in February and almost disappears in March. Both Ekman transport induced by the alongshore winter monsoon and Ekman pumping due to orographic wind stress curl are favorable for the upwelling. Transport estimates demonstrate that the month-to-month variability of Ekman transport and Ekman pumping are both consistent with that of winter coastal upwelling, but Ekman transport is two times larger than Ekman pumping in January and February. Under the influence of El Ni?o-Southern Oscillation(ENSO), the upwelling shows remarkable interannual variability: during winter of El Ni?o(La Ni?a) years, an anticyclonic(a cyclonic) wind anomaly is established in the SCS, which behaves a northeasterly(southwesterly) anomaly and a positive(negative) wind stress curl anomaly off the northwest Borneo coast, enhancing(reducing) the upwelling and causing anomalous surface cooling(warming) and higher(lower) chlorophyll concentration. The sea surface temperature anomaly(SSTA) associated with ENSO off the northwest Borneo coast has an opposite phase to that off southeast Vietnam, resulting in a SSTA seesaw pattern in the southern SCS in winter.  相似文献   

5.
北太平洋副热带东部模态水现在和未来的模拟分析   总被引:2,自引:1,他引:1  
The present climate simulation and future projection of the Eastern Subtropical Mode Water(ESTMW) in the North Pacific are investigated based on the Geophysical Fluid Dynamics Laboratory Earth System Model(GFDL-ESM2M). Spatial patterns of the mixed layer depth(MLD) in the eastern subtropical North Pacific and the ESTMW are well simulated using this model. Compared with historical simulation, the ESTMW is produced at lighter isopycnal surfaces and its total volume is decreased in the RCP8.5 runs, because the subduction rate of the ESTMW decreases by 0.82×10-6 m/s during February–March. In addition, it is found that the lateral induction decreasing is approximately four times more than the Ekman pumping, and thus it plays a dominant role in the decreased subduction rate associated with global warming. Moreover, the MLD during February–March is banded shoaling in response to global warming, extending northeastward from the east of the Hawaii Islands(20°N, 155°W) to the west coast of North America(30°N, 125°W), with a maximum shoaling of 50 m, and then leads to the lateral induction reduction. Meanwhile, the increased northeastward surface warm current to the east of Hawaii helps strengthen of the local upper ocean stratification and induces the banded shoaling MLD under warmer climate. This new finding indicates that the ocean surface currents play an important role in the response of the MLD and the ESTMW to global warming.  相似文献   

6.
Effects of Stokes production on summer ocean shelf dynamics   总被引:1,自引:0,他引:1  
A two-dimensional numerical model,which is configured on the basis of Princeton ocean model(POM),is used to study the effect of Stokes production(SP) of the turbulent kinetic energy on a density profile and Ekman transport in an idealized shelf region in summer.The energy input from SP is parameterized and included into the Mellor-Yamada turbulence closure submodel.Results reveal that the intensity of wind-driven upwelling fronts near the sea surface is weakened by the SP-associated turbulent kinetic energy input.The vertical eddy viscosity coefficient in the surface boundary layer is enhanced greatly owing to the impact of SP,which decreases the alongshore velocity and changes the distribution of upwelling.In addition,the SP-induced mixing easily suppresses the strong stratification and significantly increases the depth of the upper mixed layer(ML) under strong winds.  相似文献   

7.
利用Argo剖面浮标分析上层海洋对台风“布拉万”的响应   总被引:9,自引:2,他引:7  
In situ observations from Argo profiling floats combined with satellite retrieved SST and rain rate are used to investigate an upper ocean response to Typhoon Bolaven from 20 through 29 August 2012. After the passage of Typhoon Bolaven, the deepening of mixed layer depth(MLD), and the cooling of mixed layer temperature(MLT) were observed. The changes in mixed layer salinity(MLS) showed an equivalent number of increasing and decreasing because the typhoon-induced salinity changes in the mixed layer were influenced by precipitation, evaporation, turbulent mixing and upwelling of thermocline water. The deepening of the MLD and the cooling of the MLT indicated a significant rightward bias, whereas the MLS was freshened to the left side of the typhoon track and increased on the other side. Intensive temperature and salinity profiles observed by Iridium floats make it possible to view response processes in the upper ocean after the passage of a typhoon. The cooling in the near-surface and the warming in the subsurface were observed by two Iridium floats located to the left side of the cyclonic track during the development stage of the storm, beyond the radius of maximum winds relative to the typhoon center. Water salinity increases at the base of the mixed layer and the top of the thermocline were the most obvious change observed by those two floats. On the right side of the track and near the typhoon center when the typhoon was intensified, the significant cooling from sea surface to a depth of 200×104 Pa, with the exception of the water at the top of the thermocline, was observed by the other Iridium float. Owing to the enhanced upwelling near the typhoon center, the water salinity in the near-surface increased noticeably. The heat pumping from the mixed layer into the thermocline induced by downwelling and the upwelling induced by the positive wind stress curl are the main causes for the different temperature and salinity variations on the different sides of the track. It seems that more time is required for the anomalies in the subsurface to be restored to pretyphoon conditions than for the anomalies in the mixed layer.  相似文献   

8.
Category 5 typhoon Megi was the most intense typhoon in 2010 of the world. It lingered in the South China Sea (SCS) for 5 d and caused a significant phytoplankton bloom detected by the satellite image. In this study, the authors investigated the ocean biological and physical responses to typhoon Megi by using chlorophylla (chla) concentration, sea surface temperature (SST), sea surface height anomaly (SSHA), sea surface wind measurements derived from different satellites and in situ data. The chla concentration (>3 mg/m3) increased thirty times in the SCS after the typhoon passage in comparison with the mean level of October averaged from 2002 to 2009. With the relationship of wind stress curl and upwelling, the authors found that the speed of upwelling was over ten times during typhoon than pretyphoon period. Moreover, the mixed layer deepened about 20 m. These reveal that the enhancement of chla concentration was triggered by strong vertical mixing and upwelling. Along the track of typhoon, the maximum sea surface cooling (6-8℃) took place in the SCS where the moving speed of typhoon was only 1.4-2.8 m/s and the mixed layer depth was about 20 m in pretyphoon period. However, the SST drop at the east of the Philippines is only 1-2℃ where the translation speed of typhoon was 5.5-6.9 m/s and the mixed layer depth was about 40 m in pretyphoon period. So the extent of the SST drop was probably due to the moving speed of typhoon and the depth of the mixed layer. In addition, the region with the largest decline of the sea surface height anomaly can indicate the location where the maximum cooling occurs.  相似文献   

9.
The annual subduction rate in the South Indian Ocean was calculated by analyzing Simple Ocean Data Assimilation(SODA) outputs in the period of 1950–2008. The subduction rate census for potential density classes showed a peak corresponding to Indian Ocean subtropical mode water(IOSTMW) in the southwestern part of the South Indian Ocean subtropical gyre. The deeper mixed layer depth, the sharper mixed-layer fronts and the associated relatively faster circulation in the present climatology resulted in a larger lateral induction, which primarily dominants the IOSTMW subduction rate, while with only minor contribution from vertical pumping.Without loss of generality, through careful analysis of the water characteristics in the layer of minimum vertical temperature gradient(LMVTG), the authors suggest that the IOSTMW was identified as a thermostad, with a lateral minimum of low potential vorticity(PV, less than 200×10~(–12) m~(–1)·s~(–1)) and a low d T?dz(less than 1.5°C/(100m)). The IOSTMW within the South Indian Ocean subtropical gyre distributed in the region approximately from25° to 50° E and from 30° to 39°S. Additionally, the average characteristics(temperature, salinity, potential density)of the mode water were estimated about(16.38 ± 0.29)°C,(35.46 ± 0.04),(26.02 ± 0.04) σ_θ over the past 60 years.  相似文献   

10.
Ocean temperature responses to Typhoon Mstsa in the East China Sea   总被引:1,自引:1,他引:0  
The MASNUM wave-tide-circulation coupled model, with 21 layers in the vertical and (1/8) °horizontal resolution, was employed to investigate the oceanic responses to Typhoon Mstsa which traversed the East China Sea (ECS) during the period of 4 - 6 August, 2005. Numerical experiment results are analyzed and compared with observation. The responses of the sea surface temperature (SST), in a focused area of (27° -29°N, 121° - 124°E), include heating and cooling stages. The heating is mainly due to warm Kuroshio water transportation and downwelling due to the water accumulation. In the cooling stage, the amplitude of the simulated cold wake ( -3℃ ), located on the right side of this typhoon track, is compared quite well with that of the satellite observed SST data. The wave-induced mixing(Bv) plays a key role for the SST cooling. Bv still plays a leading role, which accounts for 36%, for the ocean temperature drop in the upper ocean of 0 - 40 m, while the upwelling is responsible for 84% of the cooling for the lower layer of 40 - 70 m. The mixed layer depth (MLD) increased quickly from 28 to 50 m in the typhoon period. However, the simulated MLD without the wave-induced vertical mixing, evolution from 13 to 32 m, was seriously underestimated. The surface wave is too important to be ignored for the ocean responses to a typhoon.  相似文献   

11.
Hydrographic, geochemical, and direct velocity measurements along two zonal (7.5°N and 4.5°S) and two meridional (35°W and 4°W) lines occupied in January–March, 1993 in the Atlantic are combined in an inverse model to estimate the circulation. At 4.5°S, the Warm Water (potential temperature θ>4.5°C) originating from the South Atlantic enters the equatorial Atlantic, principally at the western boundary, in the thermocline-intensified North Brazil Undercurrent (33±2.7×106 m3 s−1 northward) and in the surface-intensified South Equatorial Current (8×106 m3 s−1 northward) located to the east of the North Brazil Undercurrent. The Ekman transport at 4.5°S is southward (10.7±1.5×106 m3 s−1). At 7.5°N, the Western Boundary Current (WBC) (17.9±2×106 m3 s−1) is weaker than at 4.5°S, and the northward flow of Warm Water in the WBC is complemented by the basin-wide Ekman flow (12.3±1.0×106 m3 s−1), the net contribution of the geostrophic interior flow of Warm Water being southward. The equatorial Ekman divergence drives a conversion of Thermocline Water (24.58⩽σ0<26.75) into Surface Water (σ0<24.58) of 7.5±0.5×106 m3 s−1, mostly occurring west of 35°W. The Deep Water of northern origin flows southward at 7.5°N in an energetic (48±3×106 m3 s−1) Deep Western Boundary Current (DWBC), whose transport is in part compensated by a northward recirculation (21±4.5×106 m3 s−1) in the Guiana Basin. At 4.5°S, the DWBC is much less energetic (27±7×106 m3 s−1 southward) than at 7.5°N. It is in part balanced by a deep northward recirculation east of which alternate circulation patterns suggest the existence of an anticyclonic gyre in the central Brazil Basin and a cyclonic gyre further east. The deep equatorial Atlantic is characterized by a convergence of Lower Deep Water (45.90⩽σ4<45.83), which creates an upward diapycnal transport of 11.0×106 m3 s−1 across σ4=45.83. The amplitude of this diapycnal transport is quite sensitive to the a priori hypotheses made in the inverse model. The amplitude of the meridional overturning cell is estimated to be 22×106 m3 s−1 at 7.5°N and 24×106 m3 s−1 at 4.5°S. Northward heat transports are in the range 1.26–1.50 PW at 7.5°N and 0.97–1.29 PW at 4.5°S with best estimates of 1.35 and 1.09 PW.  相似文献   

12.
On the basis of hydrographic data obtained from 28 November to 27 December, 1998, the three-dimensional structure of circulation in the South China Sea (SCS) is computed using a three-dimensional diagnostic model. The combination of sea surface height anomaly from altimeter data and numerical results provides a consistent circulation pattern for the SCS, and main circulation features can be summarized as follows: in the northern SCS there are a cold and cyclonic circulation C1 with two cores C1-1 and C1-2 northwest of Luzon and an anticyclonic eddy (W1) near Dongsha Islands. In the central SCS there is a stronger cyclonic circulation C2 with two cores C2-1 and C2-2 east of Vietnam and a weaker anticyclonic eddy W2 northwest of Palawan Island. A stronger coastal southward jet presents west of the eddy C2 and turns to the southeast in the region southwest of eddy C2-2, and it then turns to flow eastward in the region south of eddy C2-2. In the southern SCS there are a weak cyclonic eddy C3 northwest of Borneo and an anti-cyclonic circulation W3 in the subsurface layer. The net westward volume transport through section CD at 119.125°E from 18.975° to 21.725°N is about 10.3 × 106 m3s−1 in the layer above 400 m level. The most important dynamic mechanism generating the circulation in the SCS is a joint effect of the baroclinicity and relief (JEBAR), and the second dynamical mechanism is an interaction between the wind stress and relief (IBWSR). The strong upwelling occurs off northwest Luzon.  相似文献   

13.
2002年春季吕宋海峡海流:观测与改进逆模式计算   总被引:10,自引:1,他引:9  
基于2002年春季航次在吕宋海峡海域锚碇测流站(20°49'57"N,120°48'12"E)200,500与800m处锚碇测流以及CTD观测,采用改进逆方法对调查海域进行海流计算.(1)主要观测的结果:1)在200m处,观测期间海流平均速度为(47.4cm/s,346°).在500m处,海流观测期间平均速度为(20.3cm/s,350°).这些都表明黑潮在吕宋海峡锚碇测流站200和500m处向西北方向入侵南海.2)在800m处,海流观测期间平均速度为(1.2cm/s,35°),它的方向为东北向.比较每层实测流结果,表明800m层海流状况与200和500m层流况不同.3)在观测期间,200,500和800m处,日平均流速在4月皆比3月时要强.4)在调查海区西部的中间区域存在一个高密、冷水中心(HDCW),其中心位置位于断面A的水文站3附近.5)在调查海区东南区域存在一个低密、暖水(LDWW)中心,其中心位置位于断面B的水文站8附近.(2)主要计算结果:1)通过断面B的偏北方向与偏南方向的流量分别为32.48×106m3/s(包括反气旋涡的流量)与3.34×106m3/s.因此通过断面B的净北向流量为29.14×106m3/s.2)通过断面A的东向与西向的流量分别为16.71×106m3/s与8.57×106m3/s(包括气旋涡的流量).因此,通过断面A的净东向流量为8.14×106m3/s.3)通过断面M北向的净流量为24.68×106m3/s.4)黑潮通过断面M后分为主流和一个支流,其主流,流量为16.54×106m3/s,流向断面C的东部分.主流通过断面C的东部分后,最后流向台湾以东海域.而其一个分支,净流量为8.14×106m3/s,在一个高密、冷水中心(HDCW)的区域以东作气旋式弯曲,然后向西北方向通过断面C的西部.因此,黑潮在断面C有两个流核.5)比较计算得到的在锚碇测流站M附近流方向与在200与500m处观测流方向为西北向,它们甚为一致.6)在断面B西侧位于550m以深水层南海水可能缓慢地从西北流向东南,通过断面B的南向流量大约为3.34×106m3/s.  相似文献   

14.
Wind data from NCEP and hydrographic data obtained from August 28 to September 10, 1994 have been used to compute circulation in the northern South China Sea and near Luzon Strait using three-dimensional diagnostic models with a modified inverse method. The numerical results are as follows: the main Kuroshio is located above 400 m levels near Taiwan’s eastern coast and above 800 m levels away from it. Near Luzon Strait above 400 m levels a branch of the Kuroshio joins with a part of the northward current, which comes from an area west of Luzon’s western coast and intrudes northwestward, then it branchs into western and eastern parts near 20°30′ N. The eastern part flows northward into an area east of Taiwan, while its western part continues to intrude northwestward, flowing through an area southwest of Taiwan. Net westward intruded volume transport through longitude Section AB at 121°00′ E from 19°00′ N to 21° 43′ N is about 3.5 × 106 m3s−1 in a layer above 400 m levels. The anticyclonic eddies W1 and W3 exist above 700 m levels east of Dongsha Islands and below 200 m levels in the eastern part of the region, respectively. The circulation in the middle region is dominated mainly by a basin-scale cyclonic gyre, and consists of three cyclonic eddies. Strong upwelling occurs in the middle region. The joint effect of baroclinity and relief and interaction between wind stress and relief both are important for real forcing of flow across contours of fH −1 in effecting the circulation pattern.  相似文献   

15.
吕宋海峡纬向海流及质量输送   总被引:30,自引:6,他引:24  
分析和计算了吕宋海峡PR21断面最近海洋调查的部分CTD资料和ADCP资料,再一次证明吕宋海峡常年存在纬向流。但对于天气尺度而言,该流型是多变的。根据高分辨率的海洋环流数值模式4a(1992~1996年)海平面高度(SSH)的输出值,运用地转关系估计了吕宋海峡纬向流的月平均值。研究表明;通过海峡流入、流出南海纬向流的深度一般达到500m左右,200m以上流速较大,平均流速为50cm/s,最大时达80cm/s以上。500m以下的纬向地转流流速较小,通常小于10cm/s.由大洋进入海峡的入流位置位于海峡的中部和南部,月平均入流最大值出现在11月,为50cm/s.位于海峡的北部和南部上层海洋的月平均出流,最大流速亦出现在11月,也为50cm/s,这与秋季北赤道流分叉位置最北(15°N),春季分叉位置最南(14°N)有关。上层流入、流出海峡的流量的月平均值分别约为10×106m3/s和5×106m3/s.当东北季风盛行时(从10月到翌年2月),流入海峡的流量远大于流出海峡的流量,两者的差可达8×106m3/s,而在其他季节两者的差仅为3×106m3/s.这说明东北季风盛行时,会有较多的水从南海南?  相似文献   

16.
The complicated flow pattern in the intermediate layer of the Luzon Strait could directly affect the efficiency of the water and energy exchange between the South China Sea (SCS) and the North Pacific. Here we present a subsurface anticyclonic eddy in the Luzon Strait deduced using observations conducted in October 2005. On the basis of the hydrographic and current measurements, an anticyclonic eddy was found in the intermediate layer, i.e., about 26.8–27.3σθ, 500–900 m. It captures part of the SCS Intermediate Water outflow in the northern Luzon Strait, and carries it to flow southward and then westward back into the SCS in the southern Luzon Strait, with volume transport of about 1.9 × 106 m3 s−1. The simulated results from Hybrid Coordinate Ocean Model also suggest the existence of this anticyclonic eddy that develops and lingers for a month long.  相似文献   

17.
-Investigations of Chi, a and primary production were carried out in the coastal upwelling area (27°'-30°30'N, west to 124°00'E) off Zhejiang in August and October 1981. The high-value areas of Chl. a were in nearshore waters and on top of the upwelling front, where Chl. a occurred the maximum 7.40 mg/m3. With the average of 1.25 g C/m2.d, the primary production was more than 2.0gC/m2·d incenteral area of the upwelling, and the whole surveyed area's output was about 4.4 ×104 tons C/d.  相似文献   

18.
Recently obtained World Ocean Circulation Experiment (WOCE) sections combined with a specially prepared pre-WOCE South Atlantic data set are used to study the dianeutral (across neutral surface) mixing and transport achieving Antarctic Intermediate Water (AAIW) being transformed to be part of the North Atlantic Deep Water (NADW) return cell. Five neutral surfaces are mapped, encompassing the AAIW from 700 to 1100 db at the subtropical latitudes.Coherent and significant dianeutral upwelling is found in the western boundary near the Brazil coast north of the separation point (about 25°S) between the anticyclonic subtropical and cyclonic south equatorial gyres. The magnitude of dianeutral upwelling transport is 10-3 Sv (1 Sv=106 m3 s-1) for 1°×1° square area. It is found that the AAIW sources from the southwestern South Atlantic and southwestern Indian Ocean do not rise significantly into the Benguela Current. Instead, they contribute to the NADW return formation by dianeutral upwelling into the South Equatorial Current. In other words, the AAIW sources cannot obtain enough heat/buoyancy to rise until they return to the western boundary region but north of the separation point. The basin-wide integration of dianeutral transport shows net upward transports, ranging from 0.25 to 0.6 Sv, across the lower and upper boundary of AAIW north of 40°S. This suggests that the equatorward AAIW is a slow rising water on a basin average. Given one order of uncertainty in evaluating the along-neutral-surface and dianeutral diffusivities from the assumed values, K=103 m2 s-1 and D=10-5 m2 s-1, the integrated dianeutral transport has an error band of about 10–20%. The relatively weak integrated dianeutral upwelling transport compared with AAIW in other oceans implies much stronger lateral advection of AAIW in the South Atlantic.Mapped Turner Angle in diagnosing the double-diffusion processes shows that the salty Central Water can flux salt down to the upper half of AAIW layer through salt-fingering. Therefore, the northward transition of AAIW can gain salt either through along-neutral-surface advection and diffusion or through salt fingering from the Central Water and heat through either along-neutral-surface advection and diffusion or dianeutral upwelling. Cabbeling and thermobaricity are found significant in the Antarctic frontal zone and contribute to dianeutral downwelling with velocity as high as −1.5×10-7 m s-1. A schematic AAIW circulation in the South Atlantic suggests that dianeutral mixing plays an essential role in transforming AAIW into NADW return formation.  相似文献   

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

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