首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Interactions of tidal constituents and the transfer of energy from the tidal frequencies to other frequencies are investigated using 3-D tidal simulations for the Indonesian seas, focusing on an area of active internal tides. Semidiurnal tides strongly affect diurnal tides; however, semidiurnal tides are essentially unaffected by diurnal tides. The semidiurnal and diurnal constituents interact with each other through non-linear interference, both destructive and constructive. Semidiurnal tides generate harmonics at nearly the diurnal frequency and higher vertical wavenumbers. In Ombai Strait, these harmonics are out of phase with the diurnal tides and interact destructively with the diurnal tides, effectively negating the diurnal response in some locations. However, this is not a general response, and interactions differ between locations. Energy is also transferred from both semidiurnal and diurnal tides to other frequencies across the spectrum, with more energy originating from semidiurnal tides. These energy transfers are not homogeneous, and the spectral responses differ between the Makassar and Ombai Straits, with the region east of Ombai showing a more active surface response compared to a more intense benthic response in Makassar. In deep water away from topography, velocity spectra generally follow the Garrett–Munk (GM) relation. However, in areas of internal tide generation, spectral density levels exceed GM levels, particularly between 4 and 8 cycles per day (cpd), indicating increased non-linear interactions and energy transfer through resonant interactions. The model indicates strong surface trapping of internal tides, with surface velocity spectra having significantly higher energy between 4 and 8 cpd even 100 km away from the prominent sill generating the internal tides.  相似文献   

2.
ABSTRACT

The generation processes and potential energy sources of internal solitary waves (ISWs) in the southern Taiwan Strait are investigated by driving a high resolution non-hydrostatic numerical model with realistic background conditions. Two main types of ISWs are clarified according to their different energy sources. One is generated by the nonlinear disintegration of remote internal tides emanating from Luzon Strait, and the other type is generated by local tide-topography interaction at the continental slope. The basic properties and evolution processes differ between these two kinds of ISWs. The waves originated from the remote internal tides at Luzon Strait have amplitudes comparable to previous field observations. In contrast, the ISWs generated locally are much weaker than observed waves, even in the presence of a steady offshore background current, which intensifies the generation of onshore ISWs. The ISWs induced by remotely generated M2 internal tides are stronger than those induced by K1 internal tides, and the fraction of internal wave energy transmitted onto the shelf is not significantly influenced by the intensity of remotely generated internal tides.  相似文献   

3.
Large-scale redistribution of sand by hydrodynamical processes in shelf seas is important for basin and coastal evolution on time scales of a thousand to tens of thousands of years. The influence of tides on the large-scale net sand-transport patterns in the North Sea has received much attention, but the influence of wind-driven flow and wind waves has hardly been investigated. Here, to establish the present-day situation and to develop a method that can also be used for palaeo-situations and forecasts for different sea levels, this influence is assessed for the present southern North Sea using a numerical flow model, a parametric wave model and a wave-averaged sand-transport formulation. Various forcing combinations are used to identify the dominant transport mechanisms: tides only, tides and wind, tides and waves, and combined tides, wind and waves. Wind forcing is applied in two ways to find an efficient, but still representative, method of incorporating this stochastic process: a statistical wind climatology and an observed time series. The results show that (i) the wind climatology yields a good approximation of the sand transport computed using the time series; (ii) wind-driven flow and waves only contribute significantly to the net sand transport by tides when acting together where tidal currents are small; and (iii) various combinations of forcings dominate the net sand transport in different regions of the southern North Sea: (a) tides dominate in the southern, middle and northwestern parts of the Southern Bight and in the region of The Wash; (b) tides, wind-driven flow and waves all are important in the northeastern part of the Southern Bight; and (c) wind-driven flow and waves dominate north of the Friesian Islands, in the German Bight and on the Dogger Bank. Qualitative comparison with observations shows good agreement.  相似文献   

4.
A three-dimensional model covering the northwest European Shelf and part of the adjacent Atlantic Ocean is used to examine the influence of water depth change upon the distribution of maximum tidal bed stress. The direction of bed stress is an indicator of sediment movement as bed load and various regions of convergence and divergence in good agreement with observations are identified. Calculations are performed with water depths reduced by 35 m, corresponding to 10 000 years before present (B.P.). Initially, the model is forced by only the M2 tide, although subsequently five constituents, namely M2, S2, N2, K1 and O1, are used for tidal forcing. Although the distribution of extreme bed stresses computed with only M2 tidal forcing is comparable to that computed with five tides, the additional tidal constituents modify the magnitude of the bed stress. In particular the diurnal tides show regions of local enhanced current amplitude in the shelf-edge region with corresponding changes in bed stress. When water depths are reduced such that the North Sea and English Channel are separated, then there is a significant change in the tidal distribution in the shallow Southern Bight which influences bed-stress distributions and hence bed-load sediment transport in the area. Besides changes in shallow regions, the distribution of tides at the shelf edge is affected. A discussion of the limitations of the present coarse-grid model in shelf-edge regions and how it can be used to provide boundary conditions for limited-area three-dimensional models that can include stratification is presented. Also the importance of stratification for sediment movement at the shelf edge is briefly discussed.Responsible Editor: Phil Dyke  相似文献   

5.
To clarify the generation and dissipation mechanisms of diurnal coastal-trapped waves (CTWs) over the Sakhalin shelf, a series of numerical experiments were conducted using a three-dimensional tidal model of the Okhotsk Sea with density stratification. The tidal model used has good reproduction owing to the careful fitting to the recent observations. The numerical experiments suggested that diurnal CTWs are primarily (~60%) generated by the conversion of tidal energy at the northern corner of the Sakhalin shelf, and further amplified by vorticity generation due to the water column oscillation from Sakhalin Bay and the influence of Kashevarov Bank. From the observations, it was found that diurnal CTWs are effectively dissipated by the strong spin-down due to bottom friction. The conventional turbulent closure model cannot reproduce the observed damping of diurnal CTWs, which raises a caution in modeling the tidal fields in high-latitude regions where diurnal CTWs exist. To resolve this underestimation of the damping, the vertical eddy viscosity was parameterized using its dependence on the observed major axis length of the diurnal tidal current ellipses, which improves the model reproduction on the damping of diurnal CTWs. The model also suggests that the spin-down effects due to friction associated with the sea-ice cover play an important role in the tidal current reduction in the region where diurnal CTWs exist, as the observations suggested.  相似文献   

6.
The Adriatic Sea general circulation model coupled to a third generation wave model SWAN and a sediment transport model was implemented in the Adriatic Sea to study the dynamics of the sediment transport and resuspension in the northern Adriatic Sea (NAS) during the Bora event in January 2001. The bottom boundary layer (BBL) was resolved by the coupled model with high vertical resolution, and the mechanism of the wave–current interaction in the BBL was also represented in the model. The study found that, during the Bora event of 13–17 January 2001, large waves with significant wave height 2 m and period of 5 s were generated by strong winds in the northwestern shelf of the Adriatic where the direction of wave propagation was orthogonal to the current. The combined motion of the wave and current in the BBL increased the bottom stress over the western Adriatic shelf, resulting in stronger sediment resuspension there. Combining stronger bottom resuspension and strong upward vertical flux of resuspended sediments due to turbulent mixing, the model predicted that sediment concentration near the Po River was much higher than that predicted by the model run without wave forcing. The study also shows that wave–current interaction in the BBL reduced the western Adriatic Coastal Currents (WACCs) in the shallower north. It is concluded that wave forcing significantly changed the sediment distributions and increased the total horizontal fluxes over the western shelf. These results signified wave effect on sediment flux and distribution in the NAS, and suggested that waves cannot be neglected in the study of dynamics of sediment transport and resuspension in the shallow coastal seas. By including the tidal forcing in the coupled model, we also examined the effect of tides on the sediment transport dynamics in the NAS.  相似文献   

7.
Greater Cook Strait (GCS) lies between the North and the South Islands of New Zealand. Its location at the convergence of the Pacific and Indo-Australian tectonic plates leads to interesting bathymetry with an adjacent shallow shelf and deep ocean trench as well as numerous crossing faults and complex shoreline geometry. Our purpose in this study is to examine tides and currents in GCS and, in particular, identify the major forcing mechanisms for the residual currents. Toward this end, we use an unstructured-grid numerical model to reproduce the tides and currents, verify these results with observations and then use the model to separate the various forcing mechanisms. The physical forcing includes nonlinear generation from tides and tidal currents, differences in sea level between the Pacific Ocean and Tasman Sea boundaries, density variations, wind stress and river discharge into GCS. Each of these mechanisms is important in different areas.  相似文献   

8.
The effects of local and remote wind forcing of water level heights in the Virginia Coast Reserve (VCR) are examined in order to determine the significant forces governing estuarine motions over subtidal time scales. Recent (1996–2008) data from tide and wind stations in the lagoon, a tide station to the north at Sandy Hook, NJ, and one offshore wind station at the Chesapeake Light Tower are examined. Sea surface height spectrum calculations reveal significant diurnal and semidiurnal tidal effects along with subtidal variations, but a suppressed inertial signal. Sea-surface heights (SSH) with 2–5 day periods at Wachapreague, VA are coherent with those at Sandy Hook and lag them in time, suggesting that southward-propagating continental shelf waves provide subtidal variability within the lagoon. The coherence between lagoon winds and sea surface height, as well as between winds and cross-lagoon sea height gradient, were significant at a relatively small number of frequency and wind direction combinations. The frequencies at which this wind forcing occurs are the tidal and subtidal bands present to the north, so that lagoon winds selectively augment existing SSH signals, but do not generate them. The impact of the wind direction is closely related to the geometry of the lagoon and bounding landmasses. The effect of wind stress is also constrained by geometry in affecting the cross-lagoon water height gradient. Water levels at subtidal frequencies are likely forced by a combination of local wind forcing, remote wind forcing and oceanic forcing modified by the complex topography of the lagoon, shelf, and barrier islands.  相似文献   

9.
This study examines connections between mean sea level (MSL) variability and diurnal and semidiurnal tidal constituent variations at 17 open-ocean and 9 continental shelf tide gauges in the western tropical Pacific Ocean, a region showing anomalous rise in MSL over the last 20 years and strong interannual variability. Detrended MSL fluctuations are correlated with detrended tidal amplitude and phase fluctuations, defined as tidal anomaly trends (TATs), to quantify the response of tidal properties to MSL variation. About 20 significant amplitude and phase TATs are found for each of the two strongest tidal constituents, K1 (diurnal) and M2 (semidiurnal). Lesser constituents (O1 and S2) show trends at nearly half of all gauges. Fluctuations in MSL shift amplitudes and phases; both positive and negative responses occur. Changing overtides suggest that TATs are influenced by changing shallow water friction over the equatorial Western Pacific and the eastern coast of Australia (especially near the Great Barrier Reef). There is a strong connection between semidiurnal TATs at stations around the Solomon Islands and changes in thermocline depth, overtide generation, and the El Niño Southern Oscillation (ENSO). TATs for O1, K1, and M2 are related to each other in a manner that suggests transfer of energy from M2 to the two diurnals via resonant triad interactions; these cause major tidal variability on sub-decadal time scales, especially for M2. The response of tides to MSL variability is not only spatially complex, it is frequency dependent; therefore, short-term responses may not predict long-term behavior.  相似文献   

10.
Based on experimental data and numerical modeling, the possible mechanisms of the effect of internal gravitational waves within the range of periods from tidal values to a few tens of minutes on crustal microdeformations in the coastal zone of the Sea of Japan are examined. The spectral analysis of oscillations in the sea level and microdeformations recorded in various seasons reveals common maximums of energy at diumal and semidiurnal periods, but the coincidence of the maximums at shorter periods is random and varies with time. The phase shifts between the surface tide and crustal deformations are also unstable in time. To explain the observed interrelations between the processes at sea and in the Earth’s crust, we modeled numerically the generation of internal tides, bores, and packets of short internal waves in terms of a nonlinear model of shallow water. It is shown that the observed effects can be caused (1) by the resonance between the wavelength of the internal tide and the shelf width and (2) by the reflection of bores and internal wave packets from a steep bottom and rocky shores or by their collapse.  相似文献   

11.
Ocean tides and resonance   总被引：1，自引：1，他引：0  
Tidal currents and surface amplitudes are calculated globally for the dominating diurnal and semidiurnal constituents using an established tidal model under a range of altered bathymetry. The purpose is to evaluate if the well-known amplification of the global tides during the Last Glacial Maximum (LGM) is related to changed propagation properties for the tidal wave or to changed damping due to removal of shelf seas. The response of the tides and tidal dissipation to future sea-level rise is also discussed. The tides in the present and LGM oceans were simulated first, followed by runs where the present day bathymetry was used but the shelf seas removed by the introduction of vertical walls or where sea level is allowed to rise. Previously reported results regarding tidal amplitudes and dissipation rates are reproduced in the control runs. The runs without shelf seas show significantly enhanced tidal amplitudes in the North Atlantic, whereas sea-level rise of 5 m above present levels show a significant shift in the amphidromic points on a local and regional scale but had a limited effect on the open ocean tides. Simulations with very large sea-level rise show a significantly decreased global tidal dissipation, whereas experiments without friction in present-day shallow water display results similar to those with no shelf seas. The results all point towards changing damping properties due to the removal of shelf seas as being the mechanism behind the LGM amplification, and they imply the importance of implementing future sea-level changes properly in tidal simulations.  相似文献   

12.
A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M2, S2, N2, K1 and O1 tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke  相似文献   

13.
利用为期一年的卫星遥感温度（SABER/TIMED）资料重建了120°E子午圈内中间层和低热层大气潮汐各主要频率分量（周日、半日和8小时潮汐）.这些主要频率分量随高度振幅增大，在97 km高度达到显著的振幅；其中迁移性周日潮汐在97 km高度出现极大振幅，然后随高度衰减.本文从考察迁移性成分和非迁移性成分各自在总潮汐中贡献角度出发，着重讨论了那些对形成该子午圈中97 km高度上整体潮汐扰动起控制作用的潮汐成分.结果显示，对周日和半日频率这两种潮汐而言，迁移性成分控制了它们的总体时空分布.在春分季节，迁移性周日潮的控制作用最显著，决定了赤道和两半球热带的活动中心；其中北半球副热带地区的季节变化形势与以往利用武汉（30°N，114°E）流星雷达风测量资料开展分析得到的结果是一致的；其他季节受非迁移性成分明显影响，例如，在本文关注的2005年中，夏至季节受（1，0）模、（1，-3）模和（1，-2）模的共同影响形成了从赤道向南延伸的活动中心，极值中心位于赤道附近，振幅达到了20 K以上，是全年的最大值.受迁移性成分控制，半日潮活动主要出现在两半球热带地区，北半球活动中心位于秋分季节（振幅达到13 K），南半球活动中心位于春分和夏至之间.其他季节受非迁移成分的影响，形成若干分布在两半球的活动中心.在本文关注的40°S～40°N范围内，与周日潮和半日潮相比，8小时潮汐具有显著较低的振幅；另外，虽然迁移性成分在一年中的大部分时间系统地分布在两半球热带地区，但是非迁移成分具有与迁移性成分相当或更大的振幅，在整体上控制了这种潮汐的时空分布.  相似文献   

14.
Semidiurnal tidal currents on the outer shelf of the Mackenzie Shelf in the Beaufort Sea were found to be strongly influenced by the locally generated baroclinic tide. Two primary factors are involved in this process: (1) the sharp shelf break along the northeastern Mackenzie Shelf, promoting the generation of vigorous internal tidal waves; and (2) the proximity to critical latitudes for M2 and N2 motions locking these waves and preventing them from leaving the source region. As a result, internal tides are resonantly trapped between the shelf and critical latitudes. The physical properties and temporal variations of tidal motions were examined using current meter measurements obtained from 1987–1988 at four sites (SS1, SS2, SS3, and SS4) offshore of the shelf break at depths of ∼200 m. Each mooring had Aanderaa RCM4s positioned at ∼35 m below the surface and ∼50 m above the bottom. Complex demodulation was used to compute the envelopes (amplitude modulation) of these components. A striking difference in the variability of clockwise (CW) and counterclockwise (CCW) tidal currents was found. The CW tides are highly variable, have greater amplitude, exhibit a burst-like character associated with wind events and contain about 80% of the total energy of the semidiurnal tidal currents. In contrast, the CCW components have a more regular temporal regime with distinct monthly, fortnightly and 10-day modulation at astronomical periodicities associated with frequency differences M2–N2 (0.03629 cpd), S2–M2 (0.06773 cpd), and S2–N2 (0.10402 cpd). Significant horizontal correlation of the CW current envelopes was found only between stations near the northeast Mackenzie Shelf, indicating this to be the main area of baroclinic internal wave generation.  相似文献   

15.
Analyses of hourly values of zonal and meridional wind near 95 km observed by meteor radar at Yambol (42.5°N, 26.6°E) during January 1991–June 1992 indicate the presence of planetary waves with prevailing periods of 1.5–2.5, 4–6, 9–10 and 16–18 days. About 20% of the whole power of atmospheric motions is connected with these waves, so they play an important role in the dynamics of the mesosphere-lower thermosphere (MLT) region. By dynamic spectral analysis applied to the hourly neutral wind and to the calculated hourly values of tidal amplitudes it has been demonstrated that there is considerable modulation of tidal amplitudes by planetary waves in the neutral wind, as this process is better expressed in the semidiurnal tides. The nonlinear interaction between tides and planetary waves is studied by bispectral analysis. The results of these analyses indicate again that the nonlinear interactions between semidiurnal tides and planetary waves with periods 2–20 days are stronger than those of the diurnal tides and planetary waves. A peculiar feature of dynamics in the MLT region above Bulgaria is the presence of strong oscillations with periods of 20 and 30 h, which indicate significant nonlinear coupling between them.  相似文献   

16.
Current meter measurements from the west coast South Island New Zealand continental slope exhibit flows with dominant time scales of between 1 and 4 weeks as well as the expected diurnal and semidiurnal tides. The alongshore Doodson filtered daily mean flow components are marginally correlated with the alongshore wind. The onshore flow components are marginally correlated with the envelope of the square of the semidiurnal and diurnal tidal flow, which is taken as a measure of energy loss from the tide.Observations from the southern flank of the Challenger Plateau, 200 km north of the continental slope observations and further offshore, also exhibit similar long-period variations. However, here the alongshore flow is more strongly correlated with alongshore wind than at the southern site.  相似文献   

17.
面向大气热力潮汐结构分析的Hough函数计算   总被引：1，自引：1，他引：0       下载免费PDF全文
作为一种刻画行星尺度大气波动基本结构的函数,Hough函数为开展有关大气潮汐全球结构的分析提供了具有物理基础的手段,它们也被作为一种数值工具应用于大气模式中以提高数值积分的效率.本文介绍了我们以拉普拉斯潮汐方程为基础,在发展求取Hough函数方面取得的进展,包括对有关计算方案和详细步骤的说明.针对39种分别具有三种主要频率,即周日、半日和1/3日,与此同时纬向波数落在[6,6]区间的潮汐成份,通过计算得到了与每一种成份对应的本征值、本征函数（Hough函数）集合.作为计算结果的示例,本文给出了针对迁移性周日潮和迁移性半日潮开展计算取得结果.这些结果说明这些函数具有正确的空间结构,与此同时还说明相关的本征值计算结果都达到了相当准确的水平.所有这些结果显示当前计算结果已经可以满足普通的应用.  相似文献   

18.
Vessel-mounted ADCP measurements were conducted to describe the transverse structure of flow between the two headland tips in Khuran Channel, south of Iran (26° 45′ N), where the highest tidal velocities in spring tides were ~?1.8 m/s. Current profiles were obtained using a 614.4 kHz TRDI WorkHorse Broadband ADCP over nine repetitions of three cross-channel transects during one semidiurnal tidal cycle. The 2.2-km-long transects ran north/south across the channel. A least-square fit to semidiurnal, quarter-diurnal, and sixth diurnal harmonics was used to separate the tidal signals from the observed flow. Spatial gradients showed that the greatest lateral shears and convergences were found over the northern channel and near the northern headland tip due to very sharp bathymetric changes in this area. Contrary to the historical assumption, the across-channel momentum balance in the Khuran Channel was ageostrophic. The current study represents one of the few examples reported where the lateral friction influences the across-channel momentum balance.  相似文献   

19.
Based on the horizontal winds measured using SKiYMET meteor wind radar during the period of June 2004–May 2007, the seasonal and interannual variability of the diurnal and semidiurnal amplitudes and phases in the mesospheric and lower thermospheric (MLT) region over a low-latitude station Trivandrum (8.5°N) are investigated. The monthly values of amplitudes and phases are calculated using a composite day analysis. The zonal and meridional diurnal tidal amplitudes exhibit both annual and semiannual oscillations. The zonal and meridional components of semidiurnal tide show a significant annual oscillation. The phase values of both diurnal and semidiurnal tides exhibit annual oscillation above 90 km. The effect of background wind in the lower atmosphere on the strength of diurnal tidal amplitudes in the MLT region is studied. The effect of diurnal tides on the background wind in the lower thermosphere is also discussed.  相似文献   

20.
Abstract

Tidal pressures and currents were measured with self‐contained capsules dropped to the sea floor for one month at distances of 175, 190, and 500 nautical miles from San Diego. These observations, together with a one‐week bottom pressure record by Filloux at 750 n miles, and three half‐week bottom current records by Isaacs et al, at intermediary distances, were analyzed for tidal components by cross‐correlation with a noise‐free reference time series. (For short records this method has some merit over classical tide analysis.) It was found that the tide decays seaward to e‐1 times the coastal amplitude over a distance of order 1000 km for the semidiurnal species, slower for the diurnal species. Tidal currents turn counterclockwise, and are polarized with maximum flow parrallel to shore in the direction of tidal propagation (320°T) at local high tide. The current amplitude is roughly 2 cm/sec for the semidiurnal component, 1 cm/sec for the diurnal component. Superimposed baroclinic tidal currents lead to poor signal: noise ratios (between 1:1 and 10:1) for the barotropic currents. In contrast, the ratio is typically 1000:1 for the bottom pressures and generally exceeds that for coastal tide stations of comparable duration. Published I.H.B. tidal constants for exposed California coastal stations indicate “upshore” (towards 320°T) propagation at 140 m/sec for semidiurnal tides. 214 m/sec for diurnal tides.

To interpret these diverse observations, we have computed the dispersion laws for all possible rotationally‐gravitationally trapped waves against a straight coast with shelf. Trapped solutions are conveniently portrayed in terms of a parameter μ such that ? = sin μ = iu/v and f = ‐ cos μ = η/v define the ellipticity and impedance of the wave motion, η, u and v being off‐shelf dimensionless elevation, normal‐to‐shore and longshore components of velocity, respectively. We then attempt to fit the observations by a superposition of the possible wave classes, all of the same tidal frequency: (a) a free Kelvin‐like edge wave with small μ (mostly trapped by rotation, but somewhat slowed by the shelf); (6) a free Poincare‐like leaky wave; and (c) a forced wave (the distortion of the sea bottom by the tide producing forces plays a significant role). The mod el can account for the main features of the observed tidal heights, and gives relative amplitudes at the coast of 54:16:4 cm for components a:b:c in the case of the semidiurnal tides, 21:24:9 cm for the diurnal tides. The results place a semidiurnal amphidrome about midway between San Diego and Hawaii. Tidal currents are not well fitted by the model, and there are problems associated with the separation of barotropic and baroclinic modes, and with the benthic boundary layer. Coastal energy dissipation is small in the sea under investigation, but a “ capacitive “ phase delay appears to be associated with Northern California harbors and inland waters.  相似文献   

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

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