首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Outputs from a high-resolution data assimilation system,the global Hybrid Coordinate Ocean Model and Navy Coupled Ocean Data Assimilation （HYCOM＋NCODA） 1/12° analysis,were analyzed for the period September 2008 to February 2012.The objectives were to evaluate the performance of the system in simulating ocean circulation in the tropical northwestern Pacific and to examine the seasonal to interannual variations of the western boundary currents.The HYCOM assimilation compares well with altimetry observations and mooring current measurements.The mean structures and standard deviations of velocities of the North Equatorial Current （NEC）,Mindanao Current （MC） and Kuroshio Current （KC） also compare well with previous observations.Seasonal to interannual variations of the NEC transport volume are closely correlated with the MC transport volume,instead of that of the KC.The NEC and MC transport volumes mainly show well-defined annual cycles,with their maxima in spring and minima in fall,and are closely related to the circulation changes in the Mindanao Dome （MD） region.In seasons of transport maxima,the MD region experiences negative SSH anomalies and a cyclonic gyre anomaly,and in seasons of transport minima the situation is reversed.The sea surface NEC bifurcation latitude （NBL） in the HYCOM assimilation also agrees well with altimetry observations.In 2009,the NBL shows an annual cycle similar to previous studies,reaching its southernmost position in summer and its northernmost position in winter.In 2010 and 2011,the NBL variations are dominantly influenced by La Ni（n）a events.The dynamics responsible for the seasonal to interannual variations of the NEC-MC-KC current system are also discussed.  相似文献   

2.
New diagnostics of the Madden–Julian oscillation (MJO) cycle in ocean temperature and, for the first time, salinity are presented. The MJO composites are based on 4 years of gridded Argo float data from 2003 to 2006, and extend from the surface to 1,400 m depth in the tropical Indian and Pacific Oceans. The MJO surface salinity anomalies are consistent with precipitation minus evaporation fluxes in the Indian Ocean, and with anomalous zonal advection in the Pacific. The Argo sea surface temperature and thermocline depth anomalies are consistent with previous studies using other data sets. The near-surface density changes due to salinity are comparable to, and partially offset, those due to temperature, emphasising the importance of including salinity as well as temperature changes in mixed-layer modelling of tropical intraseasonal processes. The MJO-forced equatorial Kelvin wave that propagates along the thermocline in the Pacific extends down into the deep ocean, to at least 1,400 m. Coherent, statistically significant, MJO temperature and salinity anomalies are also present in the deep Indian Ocean.  相似文献   

3.
This study aims to explore the relative role of oceanic dynamics and surface heat fluxes in the warming of southern Arabian Sea and southwest Indian Ocean during the development of Indian Ocean Dipole (IOD) events by using National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) daily reanalysis data and Global Ocean Data Assimilation System (GODAS) monthly mean ocean reanalysis data from 1982 to 2013, based on regression analysis, Empirical Orthogonal Function (EOF) analysis and combined with a 2? layer dynamic upper-ocean model. The results show that during the initial stage of IOD events, warm downwelling Rossby waves excited by an anomalous anticyclone over the west Indian Peninsula, southwest Indian Ocean and southeast Indian Ocean lead to the warming of the mixed layer by reducing entrainment cooling. An anomalous anticyclone over the west Indian Peninsula weakens the wind over the Arabian Sea and Somali coast, which helps decrease the sea surface heat loss and shallow the surface mixed layer, and also contributes to the sea surface temperature (SST) warming in the southern Arabian Sea by inhibiting entrainment. The weakened winds increase the SST along the Somali coast by inhibiting upwelling and zonal advection. The wind and net sea surface heat flux anomalies are not significant over the southwest Indian Ocean. During the antecedent stage of IOD events, the warming of the southern Arabian Sea is closely connected with the reduction of entrainment cooling caused by the Rossby waves and the weakened wind. With the appearance of an equatorial easterly wind anomaly, the warming of the southwest Indian Ocean is not only driven by weaker entrainment cooling caused by the Rossby waves, but also by the meridional heat transport carried by Ekman flow. The anomalous sea surface heat flux plays a key role to damp the warming of the west pole of the IOD.  相似文献   

4.
Observations have shown that the Indian Ocean is consistently warming and its warm pool is expanding, particularly in the recent decades. This paper attempts to investigate the reason behind these observations. Under global warming scenario, it is expected that the greenhouse gas induced changes in air–sea fluxes will enhance the warming. Surprisingly, it is found that the net surface heat fluxes over Indian Ocean warm pool (IOWP) region alone cannot explain the consistent warming. The warm pool area anomaly of IOWP is strongly correlated with the sea surface height anomaly, suggesting an important role played by the ocean advection processes in warming and expansion of IOWP. The structure of lead/lag correlations further suggests that Oceanic Rossby waves might be involved in the warming. Using heat budget analysis of several Ocean data assimilation products, it is shown that the net surface heat flux (advection) alone tends to cool (warm) the Ocean. Based on above observations, we propose an ocean-atmosphere coupled positive feedback mechanism for explaining the consistent warming and expansion of IOWP. Warming over IOWP induces an enhancement of convection in central equatorial Indian ocean, which causes anomalous easterlies along the equator. Anomalous easterlies in turn excite frequent Indian ocean Dipole events and cause anti-cyclonic wind stress curl in south-east and north-east equatorial Indian ocean. The anomalous wind stress curl triggers anomalous downwelling oceanic Rossby waves, thereby deepening the thermocline and resulting in advection of warm waters towards western Indian ocean. This acts as a positive feedback and results in more warming and westward expansion of IOWP.  相似文献   

5.
Responses of global ocean circulation and temperature to freshwater runoff from major rivers were studied by blocking regional runoff in the global ocean general circulation model(OGCM)developed at the Massachusetts Institute of Technology.Runoff into the tropical Atlantic,the western North Pacific,and the Bay of Bengal and northern Arabian Sea were selectively blocked.The blocking of river runoff first resulted in a salinity increase near the river mouths(2 practical salinity units).The saltier and,therefore,denser water was then transported to higher latitudes in the North Atlantic,North Pacific,and southern Indian Ocean by the mean currents.The subsequent density contrasts between northern and southern hemispheric oceans resulted in changes in major ocean currents.These anomalous ocean currents lead to significant temperature changes(1°C-2°C)by the resulting anomalous heat transports.The current and temperature anomalies created by the blocked river runoff propagated from one ocean basin to others via coastal and equatorial Kelvin waves.This study suggests that river runoff may be playing an important role in oceanic salinity,temperature,and circulations;and that partially or fully blocking major rivers to divert freshwater for societal purposes might significantly change ocean salinity,circulations,temperature,and atmospheric climate.Further studies are necessary to assess the role of river runoff in the coupled atmosphere-ocean system.  相似文献   

6.
国际Argo(Array for Real-time Geostrophic Oceanography)计划的实施,提供了前所未有的全球深海大洋0～2000 m水深范围内的海水温度和盐度观测资料,在大气和海洋科研业务中应用这一全新的资料,是深入认识大气和海洋变异、提高我国气候预测、海洋监测分析和预报能力的一个关键所在.通过开发非线性温—盐协调同化方案和利用同化高度计资料来调整模式的温度和盐度场,建立了可同化包括Argo等多种海洋观测资料的全球海洋资料变分同化系统,提高了对全球海洋的监测分析能力.实现了海洋资料同化系统与全球海气耦合模式的耦合,显著提高了短期气候预测水平.利用Argo资料改进了海洋动力模式中的物理过程参数化方案,有效提高了海洋模式对真实大洋的模拟能力和对厄尔尼诺/拉尼娜的预测能力.开发了利用Argo浮标漂流轨迹推算全球海洋表层和中层流的方法,提高了推算的全球表层流、中层流资料质量,有效弥补了洋流观测的匮乏.  相似文献   

7.
This study investigates how accurately the interannual variability over the Indian Ocean basin and the relationship between the Indian summer monsoon and the El Niño Southern Oscillation (ENSO) can be simulated by different modelling strategies. With a hierarchy of models, from an atmospherical general circulation model (AGCM) forced by observed SST, to a coupled model with the ocean component limited to the tropical Pacific and Indian Oceans, the role of heat fluxes and of interactive coupling is analyzed. Whenever sea surface temperature anomalies in the Indian basin are created by the coupled model, the inverse relationship between the ENSO index and the Indian summer monsoon rainfall is recovered, and it is preserved if the atmospherical model is forced by the SSTs created by the coupled model. If the ocean model domain is limited to the Indian Ocean, changes in the Walker circulation over the Pacific during El-Niño years induce a decrease of rainfall over the Indian subcontinent. However, the observed correlation between ENSO and the Indian Ocean zonal mode (IOZM) is not properly modelled and the two indices are not significantly correlated, independently on season. Whenever the ocean domain extends to the Pacific, and ENSO can impact both the atmospheric circulation and the ocean subsurface in the equatorial Eastern Indian Ocean, modelled precipitation patterns associated both to ENSO and to the IOZM closely resemble the observations.  相似文献   

8.
A numerical model is constructed to evaluate the effect of river diversions on the circulation of the Arctic Ocean, including the climatically important response in the extent of sea ice. The ocean model solves the primitive equations of motion in finite-difference form for the irregular geometry of the Arctic Ocean and Greenland/Norwegian Sea, using 110 km horizontal grid spacing and up to 13 unevenly spaced levels in the vertical. Annual mean atmospheric conditions and river discharges are specified from observations. The presence of sea ice is diagnosed on the basis of model ocean temperature; and the effects of sea ice on the surface fluxes of momentum, heat, and salt are included in a simplified way. Lateral exchanges at the southernmost boundary are held near observed values but respond to circulation changes in the Greenland/Norwegian Sea. Three equilibrium solutions are obtained by eighty-year integrations from simple initial conditions: the first with inflow from all rivers, the second with one-third of the inflow diverted from four major rivers (the Ob, Yenesei, Dvina, and Pechora), and the third with total diversion from those rivers. The middle case corresponds to maximal diversions which are either planned or envisioned by the Soviet Union over the next fifty years, whereas the final extreme case is run in the event that model sensitivity is low relative to that of nature.The control integration gives a good simulation of known water masses and currents. In the Central Arctic, for example, the model correctly predicts a strong shallow halocline, a relatively warm intermediate layer of Atlantic origin, and a temperature jump across the deep Lomonosov Ridge. The overall pattern of surface salinity and the margin of the pack ice are also properly simulated.When runoff into the marginal Kara and Barents Seas is diverted, either in part or in full, almost no effect on the halocline results in the Central Arctic. In particular, deep convection does not develop in the Eurasian Basin, the possibility of which was suggested by Aagaard and Coachman (1975). The vertical stability within the two marginal seas is considerably decreased by the total diversion of four rivers, but not to the point of convective overturning. The surface currents in this area change to confine the water with increased salinity to the shelf region. At deeper levels, an increased salinity tongue spreads into the deep basins of the ice-free Greenland/Norwegian Sea, where existing deep convection is slightly enhanced. As a result, there is some additional heat loss from the Atlantic layer before it enters the Central Arctic. The ice extent remains nearly the same as before within the Kara and Barents Seas. In fact, since modified bottom currents over the continental shelf bring in less heat from the Greenland Sea, an increased thickness of sea ice may result there, in spite of reduced vertical stability. These model responses are generally in agreement with those suggested by Micklin (1981) and by Soviet investigations of the effect of river diversions. These annualmean results should be regarded as tentative, pending confirmation by studies which include the seasonal cycles of runoff and atmospheric forcing.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

9.
The upper air data collected from the balloon-borne GLASS Sondes launched from the Oceanic Research Vessel (ORV) Sagar Kanya during the Intensive Field Phase of the Indian Ocean experiment (INDOEX, IFP-99;SK-141 Cruise) are utilized forstudying the variability in the mixed-layer heights observed over the western tropical Indian Ocean and central Arabian Sea. During the entire cruise, typical daytime convective mixed-layer heights (roughly corresponding to 1400 LT) obtained from V and q profiles, were observed to be in the range 200–900 m. Shallowmixed -layer heights are observed, in general, over the Inter-Tropical Convergence Zone (ITCZ). Over the central Arabian Sea, vertical profiles of V and q demonstrate a double mixed-layer structure of the marine atmospheric boundary layer (MABL), which gradually disappears close to the Indian coastline.  相似文献   

10.
Hydrographic and direct current measurements were made in the Eastern Equatorial Indian Ocean in May 2010 and April 2011 as part of the Eastern Indian Ocean Cruises(EIOC) organized by the South China Sea Institute of Oceanology(SCSIO).Analyses of the shipdrift Acoustic Doppler Current Profiler(ADCP) data indicate that the equatorial currents observed in May 2010 are characterized by a strongly eastward surface current(Wyrtki Jets,WJs) with a maximum velocity of 0.9 m s 1,while that observed in April 2011 is weak and without a consistent direction.The strongly eastward WJ transports the surface water eastward,resulting in a deeper upper mixed layer,as shown in the temperature and salinity profiles.However,it was found that the Equatorial Undercurrent(EUC) in the Eastern Indian Ocean is strong in April 2011 and weak in May 2010.The EUC was located approximately at the position of the thermocline,and it had higher salinity(up to approximately 35.5 psu) than the upper and lower waters.  相似文献   

11.
The global heat balance: heat transports in the atmosphere and ocean   总被引：10，自引：0，他引：10  
The heat budget has been computed locally over the entire globe for each month of 1988 using compatible top-of-the-atmosphere radiation from the Earth Radiation Budget Experiment combined with European Centre for Medium Range Weather Forecasts atmospheric data. The effective heat sources and sinks (diabatic heating) and effective moisture sources and sinks for the atmosphere are computed and combined to produce overall estimates of the atmospheric energy divergence and the net flux through the Earth's surface. On an annual mean basis, this is directly related to the divergence of the ocean heat transport, and new computations of the ocean heat transport are made for the ocean basins. Results are presented for January and July, and the annual mean for 1988, along with a comprehensive discussion of errors. While the current results are believed to be the best available at present, there are substantial shortcomings remaining in the estimates of the atmospheric heat and moisture budgets. The issues, which are also present in all previous studies, arise from the diurnal cycle, problems with atmospheric divergence, vertical resolution, spurious mass imbalances, initialized versus uninitialized atmospheric analyses, and postprocessing to produce the atmospheric archive on pressure surfaces. Over land, additional problems arise from the complex surface topography, so that computed surface fluxes are more reliable over the oceans. The use of zonal means to compute ocean transports is shown to produce misleading results because a considerable part of the implied ocean transports is through the land. The need to compute the heat budget locally is demonstrated and results indicate lower ocean transports than in previous residual calculations which are therefore more compatible with direct ocean estimates. A Poisson equation is solved with appropriate boundary conditions of zero normal heat flux through the continental boundaries to obtain the ocean heat transport. Because of the poor observational data base, adjustments to the surface fluxes are necessary over the southern oceans. Error bars are estimated based on the large-scale spurious residuals over land of 30 W m–2 over 1000 km scales (1012 m2). In the Atlantic Ocean, a northward transport emerges at all latitudes with peak values of 1.1±0.2 PW (1 standard error) at 20 to 30°N. Comparable values are achieved in the Pacific at 20°N, so that the total is 2.1±0.3 PW. The peak southward transport is at 15 to 20°S of 1.9±0.3 PW made up of strong components from both the Pacific and Indian Oceans and with a heat flux from the Pacific into the Indian Ocean in the Indonesian throughflow. The pattern of poleward heat fluxes is suggestive of a strong role for Ekman transports in the tropical regions.  相似文献   

12.
GTS的温盐资料在BCC_GODAS中的同化结果分析   总被引：2，自引：0，他引：2       下载免费PDF全文
分析了从GTS(全球无线通讯系统)获得的2002-2007年海洋温盐观测资料在国家气候中心第2代全球海洋资料同化系统(BCC_GODAS 2.0)中的同化结果。与SODA(简易海洋同化数据)资料的比较表明:GTS中的海洋温盐资料同化对模式温盐场的改进之处主要表现在混合层暖区的范围和中心强度、温跃层中温度槽脊的深度、温跃层附近的温度梯度以及盐度高、低值区的范围和中心强度等方面,同化后全球温盐场的均方根误差得到一定程度的降低。挑选位于不同海区的单点温盐廓线与ARGO(地转海洋学实时观测阵)观测作了进一步比较,结果表明:大多数情况下,同化后温盐廓线的均方根误差得到明显降低,模拟的温盐垂向分布特征也更为准确。与TAO(热带大气海洋观测网)资料的比较也同样表明:同化后的温盐场特征会得到一定程度改善。  相似文献   

13.
This paper investigates possible warming effects of an El Nino event on the sea surface temperature anomaly (SSTA) in the northwestern Indian Ocean. Most pure positive Indian Ocean dipole (IOD) events (without an El Nino event co-occurring) have a maximum positive SSTA mainly in the central Indian Ocean south of the equator, while most co-occurrences with an El Nino event exhibit a northwest-southeast typical dipole mode. It is therefore inferred that warming in the northwestern Indian Ocean is closely related to the El Nino event. Based on the atmospheric bridge theory, warming in the northwestern Indian Ocean during co-occurring cases may be primarily caused by relatively less latent heat loss from the ocean due to reduced wind speed. The deepened thermocline also contributes to the warming along the east coast of Africa through the suppressed upwelling of the cold water. Therefore, the El Nino event is suggested to have a modulating effect on the structure of the dipole mode in the tropical Indian Ocean.  相似文献   

14.
The influence of the tropical Indo-Pacific Ocean heat content on the onset of the Bay of Bengal summer monsoon(BOBSM) onset was investigated using atmospheric data from the NCEP and ocean subsurface temperature data from the Japan Metorology Agency(JMA).Results showed that the onset time of the BOBSM is highly related to the tropical Pacific upper ocean heat content(HC),especially in the key region of the western Pacific warm pool(WPWP),during the preceding winter and spring.When the HC anomalies in the WPWP are positive(negative),the onset of the BOBSM is usually early(late).Accompanied by the variation of the convection activity over the WPWP,mainly induced by the underlying ocean temperature anomalies,the Walker circulation becomes stronger or weaker.This enhances or weakens the westerly over the tropical Indian Ocean flowing into the BOB in the boreal spring,which is essential to BOBSM onset.The possible mechanism of influence of cyclonic/anti-cyclonic circulation over the northwestern tropical Pacific on BOBSM onset is also discussed.  相似文献   

15.
Over warm, shallow coral reefs the surface radiation and energy fluxes differ from those of the open ocean and result in modification to the marine atmospheric boundary layer via the development of convective internal boundary layers. The complex interrelationships between the surface energy balance and boundary-layer characteristics influence local weather (wind, temperature, humidity) and hydrodynamics (water temperature and currents), as well as larger scale processes, including cloud field properties and precipitation. The nature of these inter-relationships has not been accurately described for coral reef environments. This study presents the first measurements of the surface energy balance, radiation budget and boundary layer thermodynamics made over a coral reef using an eddy-covariance system and radiosonde aerological profiling of the lower atmosphere. Results show that changes in surface properties and the associated energetics across the ocean-reef boundary resulted in modification to the marine atmospheric boundary layer during the Austral winter and summer. Internal convective boundary layers developed within the marine atmospheric boundary layer over the reef and were found to be deeper in the summer, yet more unstable during the winter when cold and drier flow from the mainland enhances heat and moisture fluxes to the atmosphere. A mixed layer was identified in the marine atmospheric boundary layer varying from 375 to 1,200 m above the surface, and was deeper during the summer, particularly under stable anticyclonic conditions. Significant cloud cover and at times rain resulted in the development of a stable stratified atmosphere over the reef. Our findings show that, for Heron Reef, a lagoonal platform reef, there was a horizontal discontinuity in surface energy fluxes across the ocean-reef boundary, which modified the marine atmospheric boundary layer.  相似文献   

16.
Impacts on the coupled variability of the Indo-Pacific by including the effects of surface currents on surface stress are explored in four extended integrations of an experimental version of the Bureau of Meteorology??s coupled seasonal forecast model POAMA. The first pair of simulations differs only in their treatment of momentum coupling: one version includes the effects of surface currents on the surface stress computation and the other does not. The version that includes the effect of surface currents has less mean-state bias in the equatorial Pacific cold tongue but produces relatively weak coupled variability in the Tropics, especially that related to the Indian Ocean dipole (IOD) and El Ni?o/Southern Oscillation (ENSO). The version without the effects of surface currents has greater bias in the Pacific cold tongue but stronger IOD and ENSO variability. In order to diagnose the role of changes in local coupling from changes in remote forcing by ENSO for causing changes in IOD variability, a second set of simulations is conducted where effects of surface currents are included only in the Indian Ocean and only in the Pacific Ocean. IOD variability is found to be equally reduced by inclusion of the local effects of surface currents in the Indian Ocean and by the reduction of ENSO variability as a result of including effects of surface currents in the Pacific. Some implications of these results for predictability of the IOD and its dependence on ENSO, and for ocean subsurface data assimilation are discussed.  相似文献   

17.
The present study is based on the observations carried out over the IndianOcean from the Indian research vessel ORV Sagar Kanya during the intensive field phase of the Indian Ocean Experiment in January–March 1999. The study area spanned from 15°N to 20°S in the central Indian Ocean. Near surface variations and surface fluxes along the cruise track are presented. A comparison of near surface characteristics over the Indian Ocean and tropical west Pacific has been made. It is observed that the average difference between the sea surface temperature and air temperature at 10 m height was 0.7 °C over the study area, nearly half of that observed over the tropical west Pacific. A comparison between observed and NCEP reanalysissurface data has been made. We find good agreement between ship measured andNCEP reanalysis surface pressure, specific humidity and wind fields.On the other hand, surface air temperature in the reanalysis tends to be lowcompared to observations. The components of the net surface heat flux comparebetter in the north Indian Ocean than in the southern Indian Ocean.  相似文献   

18.
《大气与海洋》2013,51(3):187-201
Abstract

This paper investigates the formation and maintenance of the North Water Polynya, Baffin Bay in winter using a multi‐category sea‐ice model coupled with the Princeton ocean model. Monthly climatological atmospheric data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis provides the forcing. An objectively‐analysed climatology provides the initial ocean temperature and salinity. Wind stress drives the ice in a cyclonic gyre around northern Baffin Bay. Localized regions of thin ice form where wind drives ice away from coastlines or fast ice. The regions of thin ice are characterized by enhanced ice growth, exceeding 1.2 m mo?1. In the regions of thin ice, surface ocean heat flux is also enhanced and is between 30–60 W m?2. Surface heat flux is, in part, attributable to convective mixing and entrainment driven by ice growth. The surface heat flux reflects advection of the warm West Greenland Current. Heat and salt balances show that horizontal advective exchange counterbalances surface fluxes of heat and salt.  相似文献   

19.
Mean annual estimates of the oceanic poleward energy transport are obtained using a global atmospheric general circulation model. The computations are carried out by using the atmospheric model to determine the net annual heat flux into the ocean on an 8° × 10° grid. Assuming no net annual heat storage, the annual surface heat fluxes into any zonal band must be accompanied by a corresponding meridional heat transport in the ocean. Heat is transported northward at all latitudes in the Atlantic Ocean and is transported poleward in both hemispheres in the Pacific Ocean. To account for the net northward transport throughout the Atlantic, heat is transported into the Atlantic from the Indian and Pacific basins. The results are compared with several recent direct and indirect calculations of oceanic meridional heat transports.  相似文献   

20.
基于1970—2016年Hadley中心海温资料、NCEP/NCAR再分析资料和ECHAM4模式,研究了各海盆海表温度异常（SSTA）对1998和2016年这两个超级厄尔尼诺衰减年8月西北太平洋热带气旋（TC）生成及大尺度环流变化的可能影响。结果表明,热带印度洋和大西洋在1998与2016年几乎相反的SSTA型态是导致TC生成频数显著差异的主要原因之一,而热带和北太平洋SSTA在1998与2016年均分别在珠江三角洲和日本以南形成气旋性环流。1998年8月热带印度洋和大西洋SSTA产生的西北太平洋反气旋环流响应强于太平洋SSTA产生的气旋性环流异常,使西北太平洋受异常反气旋控制,减少TC的生成。2016年在三个大洋SSTA共同作用下,西北太平洋受异常气旋控制导致TC生成频数偏多。太平洋经向SSTA模在北半球副热带强迫出东西反向的跷跷板形势,在西北太平洋对流层产生的响应与实际变化相反,因此太平洋经向模对西北太平洋TC生成没有正的贡献。  相似文献   

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

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