首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The purpose of this study is to investigate the influence of tidal currents on sea ice in Spitsbergen fjords which may cause rapid decrease of the ice thickness due to erosion and melting of the ice. The effect was studied in-situ near the narrow channel connecting the Van Mijen Fjord and Lake Vallunden. The strong jet-like tidal currents in the strait driven by semidiurnal tide continue into the lake preventing ice freezing along a narrow strip during high tide and relatively warm weather. Understanding the formation of open water regions or regions with thin ice is important for the safe transportation on ice. We estimate conditions and representative time over which strong tidal current influences ice thickness along a narrow strip in solid ice. Changes of tidal phase and decrease in air temperature influence freezing of the strip in one-two days. While the tidal flow leaves the strait it overflows a shallow bar and generates internal lee waves propagating downslope and mixing the water. Tidal forcing of internal waves was measured using pressure gauges and by scanning of the ice surface during flood and ebb phases. Internal waves were measured using three types of CTD instruments and an ADCP current meter. The generation of wave packets occurs every tidal cycle when the current flows into the lake, but no generation occurs during the ebb phase of the tide because the currents over the bar slope are low. Parameters of internal waves are estimated. Model simulations confirm generation of internal wave train by the tidal current descending downslope.  相似文献   

2.
ORA-S3 oceanological reanalysis data for 1959–2011 is applied to analyze the role different factors play in forming advective heat transfer anomalies on an interannual–decadal scale in the upper mixed layer of the North Atlantic. Regions are revealed in which horizontal heat advection anomalies are determined by variations in current intensity, temperature gradients, and their joint influence. It is demonstrated that the contribution of different mechanisms responsible for advective heat transfer anomalies in the upper mixed layer to the total anomalies of advective origin varies fundamentally from one current to another in the North Atlantic. In the Gulf Stream area (after it separates from the continental slope), horizontal heat advection anomalies in the upper mixed layer result mainly from fluctuations in current intensity, while in the Caribbean Current and the Gulf Stream area (until its separation), they result from variations in the horizontal temperature gradients in the upper mixed layer. In the Labrador Current, both of these mechanisms have the same sign and approximately the same absolute values. In the East Greenland Current, they compensate each other. The contribution of anomalies in horizontal temperature gradients transferred by anomalous currents to the formation of heat transfer anomalies in the upper layer of the North Atlantic are, on the whole, relatively small throughout the water area. The areas of the North Atlantic and West Greenland currents are exceptions.  相似文献   

3.
《Ocean Modelling》2003,5(2):157-170
Model results from a regional model (BRIOS) of the Southern Ocean that includes ice shelf cavities and the interaction between ocean and ice shelves are used to derive a simple parameterization for ice shelf melting and the corresponding fresh water flux in large-scale ocean climate models. The parameterization assumes that the heat loss and fresh water gain due to the ice shelves are proportional to the difference in freezing temperature at the ice shelf edge base and the oceanic temperature on the shelf/slope area of the adjacent ocean as well as an effective area of interaction. This area is proportional to the along-shelf width of ice shelf and an effective cross-shelf distance, which turns out to be rather uniform (5–15 km) for a variety of different ice shelves. The proposed parameterization is easy to implement and valid for a wide range of circumstances. An application of the proposed scheme in a global ice ocean model (CLIO) supports our hypothesis that it can be used successfully and improves both the ocean and sea ice component of the model. This parameterization should also be used in models of the climate system that include a coupling between an ice sheet and an oceanic component.  相似文献   

4.
The effects of tidal currents (i.e., barotropic and internal tides) are important in the biogeochemistry of a coastal shelf sea. The high-frequency of currents and near-bottom temperatures collected in three consecutive southwest monsoon seasons (May, June, July and August of 2013 until 2015) is presented to reveal the role of the tidal currents to the temperature variability in the coastal shelf sea of the east coast of Peninsular Malaysia (ECPM), south of the South China Sea (SCS). The results of a spectral density and harmonic analysis demonstrate that the near-bottom temperature variability and the tidal currents are influenced by diurnal (O1 and K1) and semidiurnal (M2) tidal currents. The spectral density of residual currents (detided data) at 5, 10 and 16 m depth also shows significant peaks at the diurnal tidal frequency (K1) and small peaks at the semidiurnal tidal frequency (M2) indicating the existence of internal tides. The result of the horizontal kinetic energy (HKE) shows a strong intermittent energy of internal tides in the ECPM with the strongest energy is found at 16 m depth during a sporadic cooling event in June and July. A high horizontal cross-shore heat flux (16 m) also indicates strong intrusions of cooler water into the ECPM in June and July. During the short duration of cold pulse water observed in June and July, a cross-wavelet analysis also reveals the strong relationship between the near-bottom temperatures and the internal tidal currents at the diurnal tidal frequency. The intrusion of this cooler water is probably related to the monsoon-induced upwelling in June. It is loosely interpreted that the interaction between the strong barotropic tides and the steep slope in the central basin of the SCS under the stratified condition in southwest monsoon has generated these internal tides. The dissipation of internal tides from the slope area probably has driven the cold-upwelled water into the ECPM coastal shelf sea when the upwelling intensity is the highest in June and July.  相似文献   

5.
Every year, the during springtime heating conditions, the seasonal thermal frontal zone appears in Lake Ladoga. It features high horizontal water temperature gradients. The coastal waters, stably stratified in density, interact with the waters of the open lake that are unstably stratified because of the free convection developing in the temperature range between 0°C and the maximum density of the water at 4°C. In Lake Ladoga, the advance of the vernal frontal zone lasts about 7?C8 weeks from mid-May to the beginning of July. Both the water temperature and air temperature distributions over the water??s surface show that large spatial temperature ranges exist in the vernal front reaching more than 11°C. We investigated the spatial horizontal gradients of the water??s surface and the air temperature using a spatial grid with a resolution of 5 km. The surface water temperature and the air temperature gradients were compared with each other as well as with the temperatures in the region of varying depths. During the spring peak of the frontal activity in Lake Ladoga, most of the fronts feature mean temperatures greater than 4°C. This indicates that the thermal bar marks the offshore edge of the most extensive frontal zone.  相似文献   

6.
Because of high evaporation and low rainfall, gulfs and embayments in the arid coastal regions of Australia are endowed with highly saline water masses, with salinities ranging from oceanic values at their mouths to values as high as 60 at their heads. This presents the case of an ‘inverse estuary’. In Gulf St Vincent of South Australia, an example of such an inverse estuary, the salinity at the head ranges from ∼39 in winter to ∼42 in summer in sympathy with the seasonal changes of the local excess of evaporation over precipitation. The salient characteristics of temperature, salinity and density distributions in the Gulf are here described. Despite a seasonal reversal of spatial temperature gradients, the salinity gradient maintains the same sign and is strong enough to control a persistent density gradient from mouth to head. In parallel with the case of estuaries, a flushing time for salt is defined and estimated for the northern regions of the Gulf.Although the water mass was initially assumed to be vertically well-mixed because of tidally induced turbulence, the present observations reveal the striking feature that at neap tides the horizontal density-gradient forces a salt-wedge from the north, below the fresher water in the south.  相似文献   

7.
We discuss the results of a numerical experiment devoted to the investigation of the variability of the three-dimensional fields of temperature and current velocity brought about by the seasonal variability of external factors: solar radiation, atmospheric fields, discharge of the river Jordan, and water intake for economic necessities. We use a multilayer model in isopycnic coordinates with an upper mixed layer. We set atmospheric factors in the form of monthly average fields that are uniform over space and linearly interpolated in time. We compare the computed fields of heat flux and evaporation through the lake surface, level, temperature, and currents with data of observations. We note a qualitative agreement of temperature fields during the whole year and current velocities in winter when the lake is, in fact, barotropic. In summer when the lake is stratified, currents in the model turn out to be weaker than in observations.  相似文献   

8.
Under the influence of global warming, the sea ice in the Arctic Ocean (AO) is expected to reduce with a transition toward a seasonal ice cover by the end of this century. A comparison of climate-model predictions with measurements shows that the actual rate of ice cover decay in the AO is higher than the predicted one. This paper argues that the rapid shrinking of the Arctic summer ice cover is due to its increased seasonality, while seasonal oscillations of the Atlantic origin water temperature create favorable conditions for the formation of negative anomalies in the ice-cover area in winter. The basis for this hypothesis is the fundamental possibility of the activation of positive feedback provided by a specific feature of the seasonal cycle of the inflowing Atlantic origin water and the peaking of temperature in the Nansen Basin in midwinter. The recently accelerated reduction in the summer ice cover in the AO leads to an increased accumulation of heat in the upper ocean layer during the summer season. The extra heat content of the upper ocean layer favors prerequisite conditions for winter thermohaline convection and the transfer of heat from the Atlantic water (AW) layer to the ice cover. This, in turn, contributes to further ice thinning and a decrease in ice concentration, accelerated melting in summer, and a greater accumulation of heat in the ocean by the end of the following summer. An important role is played by the seasonal variability of the temperature of AW, which forms on the border between the North European and Arctic basins. The phase of seasonal oscillation changes while the AW is moving through the Nansen Basin. As a result, the timing of temperature peak shifts from summer to winter, additionally contributing to enhanced ice melting in winter. The formulated theoretical concept is substantiated by a simplified mathematical model and comparison with observations.  相似文献   

9.
A wide-spread bottom simulating reflector (BSR), interpreted to mark the thermally controlled base of the gas hydrate stability zone, is observed over a close grid of multichannel seismic profiles in the Krishna Godavari Basin of the eastern continental margin of India. The seismic data reveal that gas hydrate occurs in the Krishna Godavari Basin at places where water depths exceed 850 m. The thickness of the gas hydrate stability zone inferred from the BSR ranges up to 250 m. A conductive model was used to determine geothermal gradients and heat flow. Ground truth for the assessment and constraints on the model were provided by downhole measurements obtained during the National Gas Hydrate Program Expedition 01 of India at various sites in the Krishna Godavari Basin. Measured downhole temperature gradients and seafloor-temperatures, sediment thermal conductivities, and seismic velocity are utilized to generate regression functions for these parameters as function of overall water depth. In the first approach the base of gas hydrate stability is predicted from seafloor bathymetry using these regression functions and heat flow and geothermal gradient are calculated. In a second approach the observed BSR depth from the seismic profiles (measured in two-way travel time) is converted into heat flow and geothermal gradient using the same ground-truth data. The geothermal gradient estimated from the BSR varies from 27 to 67°C/km. Corresponding heat flow values range from 24 to 60 mW/m2. The geothermal modeling shows a close match of the predicted base of the gas hydrate stability zone with the observed BSR depths.  相似文献   

10.
Roles of horizontal processes in the formation of the density stratification in Hiuchi-Nada are investigated by means of a two-dimensional numerical model. In Hiuchi-Nada, vertically mixed and stratified regions are formed due to the regional difference of the tidal currents, and a tidal front is formed between the two regions. The horizontal mixing across the tidal front suppresses the development of the stratification, which is developed too much in the absence of the horizontal mixing. The moderate, realistic stratification cannot be realized in the model without the horizontal mixing. Density currents are formed due to the density distribution associated with the mixed and stratified states. These currents contribute to the horizontal mixing through the shear effect. Horizontal heat transfer from the outside water generates the vertical circulation and causes the stratification. This effect dominantly appears at the early and late stages of the stratified season. The stratification is initiated before the beginning of the surface heating and persists beyond the end, due to the horizontal heat transfer.  相似文献   

11.
During the past two million years low surface temperatures as well as episodically advancing ice sheets from Scandinavia acted on the subsurface pT-regime of northern Germany. Their likely effects on the petroleum system of Schleswig-Holstein were investigated. For the entire Quaternary mean annual ground temperature (MAGT) was reconstructed at a resolution of 1000 years by calibrating oxygen isotope records from ODP-site 659 to the climate of northern Germany of the past 120 kyr. The resulting MAGT trend served as input to an ice sheet model and a permafrost model along a 2D section crossing the petroleum bearing south-western part of Schleswig-Holstein. Here advances and retreats of the Scandinavian ice sheet during Saalian and Elsterian glaciation Stages were reconstructed. Maximum ice thicknesses of up to 1700 m and up to 20 periods of regional permafrost in northern Germany were reconstructed for the past 1.25 million years. Based on a basal heat flow of 50 mW/m2 permafrost thicknesses exceeded 100 m during most of these periods, temporarily extending down to depths of more than 300 m. Favourable surface temperatures and long durations of cold periods provided favourable conditions for onshore gas hydrate stability zones at Mittelplate. Implementing these glacial dynamics into 2D basin modelling (PetroMod, IES, Aachen, Germany) of the Mittelplate oil field reveals five phases of gas hydrate stability at depths down to 750 m. The latest of these events occurred during the Weichselian about 20 kyr ago. The effect of the ice sheets on pore pressure in the subsurface strongly depends on the hydraulic boundary conditions at the ice base (e.g. frozen vs. temperate ice sheet base). Excess pore pressure in the reservoir of more than 10 MPa during ice overriding is possible and probable. The calculated temperature effect of the Pleistocene cooling on the Mittelplate reservoir is in the range of 3–7 °C. Even today temperature in the reservoir is still lowered by about 4 °C in comparison to pre-Pleistocene times. Despite the fact that a significant influence of glacial effects on petroleum generation can be ruled out at Mittelplate, we state that pT-effects in reservoirs related to glacial processes in formerly glaciated areas have been underestimated in the past.  相似文献   

12.
建立了一个太湖梅梁湾三线水动力学模型,模拟了框架湾的水平及垂直流场分布。结果表明:(1)表面流场与风向一致,而底层流场与表面流场的方向完全相反,表现为很明显的补偿流;(2)水平流速基本上自表层向下递减,过渡层的流速比表层和底层小;(3)在风场的作用下可产生垂直环流系统,其中在东南风的作用下产生逆时针的垂直环流(由南向北看),而在西北风的作用下产生顺时针垂直环流;(4)垂直速度自岸边向湖中心递减,其量级远小于水平流速。  相似文献   

13.
Time series of velocity and water temperature were measured at three stations on the continental shelf, on the shelf margin and on the slope off the northwest Tokunoshima in December 1980 to study influences of the slope on tides.Tidal currents with semidiurnal periods were dominant at the stations on the shelf and shelf margin. However, semidiurnal components in temperature fluctuations were dominant at the stations on the shelf margin and the slope. We estimated horizontal currents due to semidiurnal internal tides from the vertical distribution of water density and temperature, assuming that the temperature fluctuations were caused by the vertical displacement of water particles due to semidiurnal internal tides. The tidal ellipses at the station on the shelf and the phase relation of the tidal currents between the two stations on the shelf and shelf margin indicated that the M2 surface tide on the shelf was a Sverdrup wave propagating to the northwest.Semidiurnal tidal currents on the slope were also caused by tides of surface and internal modes. Furthermore, the axis of the tidal ellipse was not perpendicular to the co-tidal line estimated by Ogura (1934) but rather parallel to the isobaths on the slope, which shows a striking effect of the bottom topography on the tidal currents.  相似文献   

14.
The circulation, water masses and sea-ice of Baffin Bay   总被引:2,自引:0,他引:2  
The oceanographic, meteorological and sea-ice conditions in Baffin Bay are studied using historical hydrographic, satellite and meteorological data, and a set of current meter data from a mooring program of the Bedford Institute of Oceanography. Baffin Bay is partially covered by sea-ice all year except August and September. The interannual variation of the ice extent is shown to be correlated with winter air temperature. Available hydrographic data were used to study the water masses and the horizontal and vertical distribution of temperature/salinity. Three water masses can be identified – Arctic Water in the upper 100–300 m of all regions except the southeast, West Greenland Intermediate Water at 300–800 m in most of the interior of Baffin Bay, and Deep Baffin Bay Water in all regions below 1200 m. The temperature and salinity in Baffin Bay have limited seasonal variability except in the upper 300 m of eastern Davis Strait, northern Baffin Bay and the mouth of Lancaster Sound. Summer data have a temperature minimum at 100 m, which suggests winter convection does not penetrate deeper than this depth. Current meter data and results of a circulation model indicate that the mean circulation is cyclonic. The seasonal variation of the currents is complex. Overall, summer and fall tend to have stronger currents than winter and spring at all depths. Among the different regions, the largest seasonal variation occurs at the mouth of Lancaster Sound and the Baffin Island slope. Model generated velocity fields show a basic agreement with the observed currents, and indicate strong topographic control in the vicinity of Davis Strait and on the Greenland shelves. The model also produces a southward counter current on the Greenland slope, which may explain the observed high horizontal shears over the Greenland slope. Estimates of the volume and fresh water transports through Lancaster, Jones and Smith Sounds are reviewed. Transports through Davis Strait are computed from the current meter data. The balance of freshwater budget and sensitivity of the thermohaline circulation to freshwater transport are discussed.  相似文献   

15.
Seasonal variability of interannual fluctuations of the heat balance components of the upper quasi-homogeneous ocean layer (UQL) in the North Atlantic is analyzed by processing the reanalysis data set for the period of 1959–2011. It is shown that interannual variations in the components of the UQL heat budget are characterized by pronounced regional features in all seasons. In the tropics and subtropics, heat balance is quasistationary and is determined by the nonlocal processes, such as heat advection and horizontal mixing. In the subpolar latitudes, nonstationarity (in the spring) and heat fluxes at the UQL boundary (in the autumn and in the winter) are also important. A major role in the interannual variability of the UQL temperature in the vicinity of jet currents of the Gulf Stream type is played in all seasons by the fluctuations of horizontal heat advection. However, the contribution of interannual fluctuations of the individual components of the heat budget to variability of the UQL temperature varies considerably in different seasons. The interannual fluctuations of the local variation in the UQL temperature are characterized by the largest variability in the spring and the lowest one in the autumn. The greatest contribution of the variations in the horizontal heat advection to the change in the UQL temperature at the interannual scale is recorded in the winter, and the lowest one is in the summer. The contribution of the interannual variations in the heat fluxes at the UQL upper boundary to the variability of the UQL temperature is the highest in the summer and the lowest in the autumn. Fluctuations of the heat fluxes at the UQL lower boundary do not have a significant impact on the interannual variations in the UQL temperature for the whole water area. The exception is small areas in the region of the formation of the North Atlantic deep water in the autumn–winter period and in the vicinity of the Equatorial Counter Current in the spring–summer period.  相似文献   

16.
Based on hydrographic data obtained at an ice camp deployed in the Makarov Basin by the 4th Chinese Arctic Research Expedition in August of 2010, temporal variability of vertical heat flux in the upper ocean of the Makarov Basin is investigated together with its impacts on sea ice melt and evolution of heat content in the remnant of winter mixed layer(r WML). The upper ocean of the Makarov Basin under sea ice is vertically stratified. Oceanic heat flux from mixed layer(ML) to ice evolves in three stages as a response to air temperature changes, fluctuating from 12.4 W/m2 to the maximum 43.6 W/m2. The heat transferred upward from ML can support(0.7±0.3) cm/d ice melt rate on average, and daily variability of melt rate agrees well with the observed results. Downward heat flux from ML across the base of ML is much less, only 0.87 W/m2, due to enhanced stratification in the seasonal halocline under ML caused by sea ice melt, indicating that increasing solar heat entering summer ML is mainly used to melt sea ice, with a small proportion transferred downward and stored in the r WML. Heat flux from ML into r WML changes in two phases caused by abrupt air cooling with a day lag. Meanwhile, upward heat flux from Atlantic water(AW) across the base of r WML, even though obstructed by the cold halocline layer(CHL), reaches0.18 W/m2 on average with no obvious changing pattern and is also trapped by the r WML. Upward heat flux from deep AW is higher than generally supposed value near 0, as the existence of r WML enlarges the temperature gradient between surface water and CHL. Acting as a reservoir of heat transferred from both ML and AW, the increasing heat content of r WML can delay the onset of sea ice freezing.  相似文献   

17.
In high-latitude oceans with seasonal ice cover, the ice and the low-salinity mixed layer form an interacting barrier for the heat flux from the ocean to the atmosphere. The presence of a less dense surface layer allows ice to form, and the ice cover reduces the heat loss to the atmosphere. The ice formation weakens the stability at the base of the mixed layer, leading to stronger entrainment and larger heat flux from below. This heat transport retards, and perhaps stops, the growth of the ice cover. As much heat is then entrained from below as is lost to the atmosphere. This heat loss further reduces the stability, and unless a net ice melt occurs, the mixed layer convects. Two possibilities exist: (1) A net ice melt, sufficient to retain the stability, will always occur and convection will not take place until all ice is removed. The deep convection will then be thermal, deepening the mixed layer. (2) The ice remains until the stability at the base of the mixed layer disappears. The mixed layer then convects, through haline convection, into the deep ocean. Warm water rises towards the surface and the ice starts to melt, and a new mixed layer is reformed. The present work discusses the interactions between ice cover and entrainment during winter, when heat loss to the atmosphere is present. One crucial hypothesis is introduced: “When ice is present and the ocean loses sensible heat to the atmosphere and to ice melt, the buoyancy input at the sea surface due to ice melt is at a minimum”. Using a one-dimensional energy-balance model, applied to the artificial situation, where ice melts directly on warmer water, it is found that this corresponds to a constant fraction of the heat loss going to ice melt. It is postulated that this partitioning holds for the ice cover and the mixed layer in the high-latitude ocean. When a constant fraction of heat goes to ice melt, at least one deep convection event occurs, before the ice cover can be removed by heat entrained from below. After one or several convection events the ice normally disappears and a deep-reaching thermal convection is established. Conditions appropriate for the Weddell Sea and the Greenland Sea are examined and compared with field observations. With realistic initial conditions no convection occurs in the warm regime of the Weddell Sea. A balance between entrained heat and atmospheric heat loss is established and the ice cover remains throughout the winter. At Maud Rise convection may occur, but late in winter and normally no polynya can form before the summer ice melt. In the central Greenland Sea the mixed layer generally convects early in winter and the ice is removed by melting from below as early as February or March. This is in agreement with existing observations.  相似文献   

18.
“南四湖渔湖民庄台工程”为中央水电部部属大型工程。该工程解决了6.8万渔湖民居住安全问题,免除了渔湖民世世代代受冰凌、风浪的威胁,为渔湖民的生产、生活和发展湖区教育创造了良好条件。有了庄台,人力集中了,又有了办厂基地,为湖产品流通加工、湖内水生植物开发利用打下了基础,素有“日出斗金”之称的南四湖将更加熠熠生辉。  相似文献   

19.
Under strong surface wind forcing during winter, direct current observations in the northern Sea of Japan show the existence of strong near-inertial currents in the deep water that is characterized by the extremely homogeneous vertical structures of temperature and salinity. However, the mechanism generating internal waves in the deep water of the northern Sea of Japan has not been well understood. In this study, to clarify the dynamical link between the surface wind forcing and near-inertial currents in the deep water of the northern Sea of Japan, we drive a general circulation model taking into account realistic wind stress, ocean bottom and land topography. In the northern Sea of Japan, the numerical results show that vertically coherent horizontal currents with a speed of ~ 0.05 m s?1 are excited throughout the homogeneous deep water. A two-layer model successfully reproduces the pattern of the horizontal current velocities shown by the general circulation model, indicating that internal waves emanate westward from the northwestern coast of Japan through coastal adjustment to the strong wind forcing event and, while propagating into the ocean interior, they excite evanescent near-inertial response throughout the lower layer below the interface.  相似文献   

20.
Chlorophyll fluorescence and temperature were measured with a horizontal resolution of 4 m on two isopycnals and two isobars along a 1000 km meridional transect in the eastern North Pacific. Probability density functions (PDFs) of the magnitude of fluctuations in temperature and fluorescence are compared at each of a large range of length scales. While they are nearly identical at larger scales, the shapes of the PDFs of fluorescence and temperature differ significantly at small scales (∼1 km and less). This difference may indicate that at these scales, temperature and fluorescence distributions are controlled by distinct mechanisms. By comparing the phase of wavelet transforms of each of these two tracers, the tendency for temperature and fluorescence gradients to line up is investigated over a large range of spatial scales. At horizontal length scales of order 10 km and larger, the wavelet phase difference between temperature and fluorescence tends to be close to 0° or 180°—that is, the gradients tend to align, either in phase or 180° out of phase. At smaller scales, the distribution of phase difference on isopycnals is uniform—there is no tendency for gradients to coincide. Simple stirring models demonstrate that the locations of enhanced gradients in all tracers would be expected to coincide where the strain has been greatest. However, the directions of enhanced gradients formed in this way may be either parallel or antiparallel, depending on initial conditions and on the direction of the strain. These analyses suggest the horizontal distributions of temperature and fluorescence at large scales were both governed by advective processes. At intermediate scales, of order 1–10 km, a renovating wave model suggests that gradient alignment could be destroyed by turbulent diffusion or by rapid local phytoplankton growth.  相似文献   

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

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