首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 923 毫秒
1.
Monte Carlo simulation of wave spectra was carried out to provide an assessment of JONSWAP spectral model and parameters. The simulation method is found to be satisfactory because (a) it excludes the spectral variability due to geophysical factors from the sampling errors in the spectral estimates and the statistical uncertainty in determining the model parameters; and (b) the simulated spectra can represent ideal spectral estimates where the sampling errors have been minimized by increasing the degrees of freedom of the spectra. The latter (b) allows both the magnitude of sampling errors to be evaluated and errors due to statistical uncertainty to be isolated. Thus, the stimulation study provides a useful error analysis to assess the JONSWAP spectral model and parameters. For instance, it is found from the results that the sampling errors could be as high as 20% while errors due to uncertainty in determining the model parameter could be as high as 17%. However, the overall errors may be reduced to the minimum of approx. 15% if the simulated spectra have 80 degrees of freedom and constant values of σa and σb i.e. σa = 0.07 and σb = 0.09. This implies that the maximum accuracy of 85% may be achieved in JONSWAP spectral model even though the α parameter has been underestimated by about 1.5%. The overestimated values of γ might come from the underestimated α and the biased φm estimator caused by the statistical uncertainty in the presence of a sharp spectral peak. Although the scale parameters (α and φm) exhibit smaller errors and variability than the shape parameters (ψ, σa and σb), they are more sensitive to the degrees of freedom of the spectra and their estimators are not better than the estimators of shape parameters. The simulation experiments have also shown that simulated spectra at 20–40 degrees of freedom contain a substantial amount of sampling errors. Therefore, the measured wave spectra at the same degrees of freedom (20–40) are not suitable and should not be used for evaluating the accuracy of any wave spectral model.  相似文献   

2.
The salinity minimum frequently occurring in the Mixed Water Region between the Oyashio and Kuroshio Fronts seems to originate from the salinity minimum at the density of 26.8σθ called the North Pacific Intermediate Water. We examined water exchange of this region with the Oyashio and the Kuroshio Extension using mixing ratio RK defined as (θ - θOY)/(θK - θOY) × 100, where θOY, θK, and θ represent potential temperature of the Oyashio and Kuroshio Waters and their mixture on the isopycnal surfaces, respectively. CTD data were obtained by repeated observation from January 1990 to May 1991. RK increases southward from the Oyashio Front to the Kuroshio Front with the range of −20 to 120%. The gradient of RK on the isopycnal surfaces is large around the Oyashio Front above the 26.8σθ surface, while it is large around the Kuroshio Front below it. This agrees with the average RK in the Mixed Water Region decreasing greatly with the increase of density at densities less dense than 26.8σθ. We calculated thickness and volume transport of the Oyashio between the isopycnal surfaces near the coast of Hokkaido. They increase largely with density at densities less dense than 26.8σθ. It is supposed that the salinity minimum in the Mixed Water Region is the upper limit of the water largely influenced by the Oyashio Water. Its density could depend only on the density structure of the Oyashio.  相似文献   

3.
Studies of mixing were done at the northern flank of Georges Bank in the summer and autumn of 1988. Two time-series of the evolution and intensity of microstructure were examined over a tidal period in the context of tidal forcing and the evolution of the density and velocity field at the site. From the CTD, ADCP and microstructure observations (EPSONDE) on Georges Bank, several interesting features of the mixing processes were found. High dissipation and diffusivity regions appear near the bottom of the Bank. Turbulence near the bottom is highest in intensity and reaches farthest from the bottom at peak tidal flow and diminishes in intensity and vertical extent as the flow decreases. The thickness of the bottom turbulent layer has its maximum value when the flow is strongest and the stratification is weakest. Characterization of the dissipation rate and turbulent diffusivities in respect to buoyancy frequency N, current shear S, Richardson Number Ri and ε/νN2 was done. Dissipation and χT showed little dependence on shear or N2 but decreased at larger Ri. χt was found to be higher in regions of higher N2 and increased as ε/νN2 increased. KT, K and Kν, were all highest near the bottom in excess of 10−2m2s−1 and decreased towards the surface. There was little suggestion of a dependence of mixing efficiency on S2, Ri or ε/νN2, but some indication that Γ decreases with decreasing N2.  相似文献   

4.
The original prognostic equations for the JONSWAP-spectrum contained inconsistencies. A subsequent paper (Hasselmann et al., 1976, J. phys. Oceanogr.6, 200–208) attempted to regularise the situation. This paper shows that there were still inconsistencies in the prognostic equations giving overestimations of the first moment of the spectrum and consequently the significant wave height. The prognostic equations are reworked systematically and results presented. It is shown that variable σa, σb and γ must be used to achieve consistent results  相似文献   

5.
Gross sedimentation rates (GSR) have been measured using sediment traps placed at nine different levels above the bed (0·3, 0·5, 0·8, 1·0, 2·0, 4·0, 6·0, 8·0 and 10·0 m). The sediment traps were deployed for 1·25 years and recovered 28 times during the study period. Low average GSR values of 5·5 g m-2 day-1 were obtained at 10·0 m, and high average GSR values of 114·8 g m-2 day-1 were obtained at 0·3 m. An expression for the eddy diffusion coefficient of suspended particulate matter (Ks), based on the measured GSR is given. The expression has been used for modelling of Ks at the different trap levels above the bed. High values (≈42 cm2 s-1) of Ks were obtained at the upper traps, whereas low values (≈2 cm2 s-1) were obtained near the bed. Comparison between level of turbulent energy in terms of shear stress at the boundaries of the water column, i.e. from the wind and the bed flow, showed that wind energy exceeded that of the bed flow by a factor 16. At 5·0 m Ks was positively correlated (r=0·66) to the eddy diffusion coefficient of momentum (Km) derived from the wind energy transfer to the water, giving an average β of 0·5 for KsKm. The density difference between surface and bottom waters has been designated a parameter of stratification, and is discussed in relation to variations of Ks and Km .  相似文献   

6.
The chemical speciation of Cu and Zn was investigated by voltammetric titration methods in the surface waters (10 m) of the western Black Sea during an Istanbul–Sevastopol cruise conducted in November 1998. Supporting parameters (temperature (T), salinity (S), pH, alkalinity (Alk), suspended particulate matter (SPM) and dissolved and particulate 234Th) were obtained in order to distinguish hydrographic features against involvement of the metals in biogeochemical processes. In the Turkish continental slope region, the cruise track intersected a narrow vein of colder water originating on the western shelf. The core of this cold water vein was characterised by a relatively low salinity, higher specific alkalinity and higher metal (especially Cu) and metal-binding ligand concentrations.A very large portion of Cu (93–99.8%) and Zn (82–97%) was organically complexed. The degree of complexation was highest in shelf waters and lowest in the central gyre. Titration data for Cu were modelled by two classes of organic binding ligands characterised by (CL1=3–12 nM, log K1′=13.1–13.9) and (CL2=20–70 nM, log K2′=9.4–11.2). These ligands occurred mainly in the ‘dissolved’ phase, as defined by 0.4-μm filtration. The stronger Cu-binding ligand seemed to be produced in situ in response to Cu concentration, whereas the weaker Cu-binding ligand appeared to be derived from terrestrial sources and/or reducing shelf sediments. Titration results for Zn were generally represented by one class of ligands (CL1=8–23 nM, log K1′=9.4–10.2), which were almost uniformly distributed between the ‘dissolved’ (78±8%) and the particulate phase (22±8%). The concentration of these strong Zn-binding ligands showed a very good correlation with SPM (r2=0.64), which improved when the dissolved ligands alone were considered (r2=0.78). It is hypothesised that these ligands were produced in situ by the bacterial breakdown of particulate organic matter.  相似文献   

7.
Vessel-based observations of water column structure and flow near a shallow-crested seamount are used to quantify the physical disturbance induced by that seamount. The implications of this disturbance on the ‘feeding hole’ hypothesis are then examined based on data from moored thermistors and acoustic current profilers, as well as vessel-based acoustic sounding-derived biomass estimates, currents, and conductivity-temperature-depth profiles from a 55 km square grid of stations around the crest. Mean currents in the region of 0.2 m s?1 are comparable to observations from surface drifters whereas the semidiurnal tidal flow amplitude was one third of this. Thorpe Scale-based estimates of energy dissipation rate were in the range 10?9 to 2.10?8 W kg?1 and vertical diffusivities Kz were in the range 10?4 to 10?3 m2 s?1. Turbulence levels were higher upstream of the seamount–likely due to the influence of nearby seamounts Rumble IV and V. There was no evidence of a Taylor Cap in the Rumble III velocity field. The sounder data provide some evidence of a feeding hole and analysis based on diffusivities suggests that this might persist downstream of the seamount for as much as 7 days.  相似文献   

8.
The interactions of Fe(II) and Fe(III) with the inorganic anions of natural waters have been examined using the specific interaction and ion pairing models. The specific interaction model as formulated by Pitzer is used to examine the interactions of the major components (Na+, Mg2+, Ca2+, K+, Sr2+, Cl, SO4, HCO3, Br, CO32−, B(OH)4, B(OH)3 and CO2) of seawater and the ion pairing model is used to account for the strong interaction of Fe(II) and Fe(III) with major and minor ligands (Cl, SO42−, OH, HCO3, CO32− and HS) in the waters. The model can be used to estimate the activity and speciation of iron in natural waters as a function of composition (major sea salts) and ionic strength (0 to 3 M). The measured stability constants (KFeX*) of Fe(II) and Fe(III) have been used to estimate the thermodynamic constants (KFeX) and the activity coefficient of iron complexes (γFeX) with a number of inorganic ligands in NaClO4 medium at various ionic strengths: In(KFeXFeγX) = InKFeX − In(γFeX) The activity coefficients for free ions (γFe, γx) needed for this extrapolation have been estimated from the Pitzer equations. The activity coefficients of the ion pairs have been used to determine Pitzer parameters (BFeX, BFeX0, CFeXφ) for the iron complexes. These results make it possible to estimate the stability constants for the formation of Fe(II) and Fe(III) complexes over a wide range of ionic strengths and in different media. The model has been used to determine the solubility of Fe(III) in seawater as a function of pH. The results are in good agreement with the measurements of Byrne and Kester and Kuma et al. When the formation of Fe organic complexes is considered, the solubility of Fe(III) in seawater is increased by about 25%.  相似文献   

9.
The natural human female hormones oestrone and 17β-oestradiol have been implicated in the disruption of endocrine systems in some wildlife adjacent to sewage effluents. The sorption behaviour of these two compounds under estuarine conditions was studied by spiking either 2.55 μg of oestrone or 2.65 μg of 17β-oestradiol in kinetic experiments. In equilibrium experiments, 3 ng of oestrone or 3.2 ng of 17β-oestradiol was added in each of the centrifuge tubes. Sorption onto sediment particles was relatively slow, with sorption equilibrium being reached in about 70 and 170 h for oestrone and 17β-oestradiol, respectively. The effects of a variety of environmental parameters on sorption were studied including salinity, sediment concentration (SC), the presence of a third phase, particle size and, also, surfactant concentrations. Results show that although salinity did not induce any statistically significant effect on the sorption of 17β-oestradiol, it did statistically enhance the sorption of oestrone, and a salting constant of 0.3 l mol−1 was derived. The partition coefficient for both compounds decreased with increasing sediment concentration, a phenomenon that has been widely reported and attributed to the presence of colloids (which could enhance dissolved concentrations). In this paper, the true partition coefficients for sediment particles (Kptrue) and colloidal particles (Kctrue) have been calculated, and a Kptrue value of 141 and 102 ml g−1 was obtained for oestrone and 17β-oestradiol, respectively. In addition, Kctrue values for oestrone (222×102 ml g−1) and 17β-oestradiol (135×102 ml g−1) were two orders of magnitude higher than their respective Kptrue values, suggesting that the colloidal particles are significantly stronger sorbents for natural oestrogens than sediment particles. Particles of different sizes were found to have different partition coefficients due to the strong relationships between partition coefficients for the two compounds and particulate organic carbon (POC) contents and specific surface areas (SSAs). The presence of a surfactant was shown to reduce the partition coefficients for the two compounds, although its concentrations being used were higher than those normally found in the natural environment.  相似文献   

10.
基于Thorpe尺度方法,利用CTD数据,计算了南极普里兹湾海域的Thorpe尺度和湍流扩散系数,分析了观测区域(64°~69°S,66°~80°E)湍流翻转现象的强弱及分布。结果表明,在海底和地形粗糙区存在较大的Thorpe尺度(较强湍流翻转)和湍流扩散系数,湍流扩散系数最大值能达到10-2m2/s量级,比平坦开阔海洋高2~3个数量级,部分观测站位的湍流扩散系数和湍动能耗散率表现出大-小-大的垂向分布结构,总水深较深的区域尤为明显;深水区域的浮力频率在海表面到500 m层比较大,浅水区域该现象不明显;湍动能耗散率在(67.25°S,73°E)周围和经度为78°E的各站位都表现相对较大,能达到10-6 w/kg量级,个别站位甚至能达到10-5 w/kg量级。  相似文献   

11.
The sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide at low ionic strength (0.01 M ≤ I ≤ 0.09 M) was investigated over a wide range of pH (3.9 ≤ pH ≤ 7.1). YREE distribution coefficients, defined as iKFe = [MSi]T / (MT[Fe3+]S), where [MSi]T is the concentration of YREE sorbed by the precipitate, MT is the total YREE concentration in solution, and [Fe3+]S is the concentration of precipitated iron, are weakly dependent on ionic strength but strongly dependent on pH. For each YREE, the pH dependence of log iKFe is highly linear over the investigated pH range. The slopes of log iKFe versus pH regressions range between 1.43 ± 0.04 for La and 1.55 ± 0.03 for Lu. Distribution coefficients are well described by an equation of the form iKFe = (Sβ1[H+]− 1 + Sβ2[H+]− 2) / (SK1[H+] + 1), where Sβn are stability constants for YREE sorption by surface hydroxyl groups and SK1 is a ferric hydroxide surface protonation constant. Best-fit estimates of Sβn for each YREE were obtained with log SK1 = 4.76. Distribution coefficient predictions, using this two-site surface complexation model, accurately describe the log iKFe patterns obtained in the present study, as well as distribution coefficient patterns obtained in previous studies at near-neutral pH. Modeled log iKFe results were used to predict YREE sorption patterns appropriate to the open ocean by accounting for YREE solution complexation with the major inorganic YREE ligands in seawater. The predicted log iKFe′ pattern for seawater, while distinctly different from log iKFe observations in synthetic solutions at low ionic strength, is in good agreement with results for natural seawater obtained by others.  相似文献   

12.
13.
Two models, a spectral refraction model (Longuet-Higgins) and a parabolic equation method (PEM) refraction-diffraction model (Kirby), are used to simulate the propagation of surface gravity waves across the Southern California Bight. The Bight contains numerous offshore islands and shoals and is significantly larger (≈ 300 km by 300 km) than regions typically studied with these models. The effects of complex bathymetry on the transformation of incident wave directional spectra, S0(f,θ0), which are very narrow in both frequency and direction are difficult to model accurately. As S0(f,θ0) becomes broader in both dimensions, agreement between the models improves and the spectra predicted at coastal sites become less sensitive to errors in the bathymetry grid, to tidal changes in the mean water depth, and to uncertainty in S0(f,θ0) itself. The smoothing associated with even relatively narrow (0.01 Hz-5° bandwidth) S0(f,θ0) is usually sufficient to bring the model predictions of shallow water energy into at least qualitative agreement. However, neither model is accurate at highly sheltered sites. The importance of diffraction degrades the predictions of the refraction model, and a positive bias [O (10%) of the deep ocean energy] in the refraction-diffraction model estimates, believed to stem from numerical “noise” (Kirby), may be comparable to the low wave energy. The best agreement between the predicted spectra generally occurs at moderately exposed locations in deeper waters within the Bight, away from shallow water diffractive effects and in the far-field of the islands. In these cases, the differences between the models are small, comparable to the errors caused by tidal fluctuations in water depth as waves propagate across the Bight. The accuracy of predicted energies at these sites is likely to be limited by the uncertainty in specifying S0(f,θ0).  相似文献   

14.
Analyses of the concentration product (Ca2+) × (CO32−) in the pore waters of marine sediments have been used to estimate the apparent solubility products of sedimentary calcite (KSPc) and aragonite (KSPa) in seawater. Regression of the data gives the relation In KPSPc = 1.94 × 10−3 δP − 14.59 The 2°C, 1 atm value of KSPc is, then, 4.61 × 10−7 mol2 l−2. The pressure coefficient yields a at 2°C of −43.8 cm3 atm−1. A single station where aragonite is present in the sediments gives a value of KSPa = 9.2 × 10−7 (4°C, 81 atm). The calcite data are very similar to those determined experimentally by Ingle et al. (1973) for KSPc at 2°C and 1 atm. The calculated is also indistinguishable from the experimental results of Ingle (1975) if is assumed to be independent of pressure.  相似文献   

15.
The annual subduction rate of the North Pacific was calculated based on isopycnally averaged hydrographic climatology (HydroBase), high-resolution winter mixed-layer climatology (NWMLC), and various wind stress climatologies from ship reports, numerical weather prediction products, and satellite products. The calculation was performed using Lagrangian coordinates in the same manner as in previous works, except a less smoothed oceanic climatology (HydroBase and NWMLC) was used instead of a World Ocean Atlas. Differences in the wind stress climatologies have very little effect on subduction rate estimates. The subduction rate census for density classes showed peaks corresponding to subtropical mode water (STMW), central mode water (CMW), and eastern subtropical mode water (ESTMW). The deeper mixed layer and the associated sharper mixed-layer fronts in the present climatology resulted in a larger lateral induction, which boosted the subduction rate, especially for the potential density anomaly (σθ) range of the lighter STMW (25.0 < σθ < 25.2 kg m−3) and lighter CMW (26.0 < σθ < 26.2 kg m−3), compared to previous estimates. The renewal time of permanent pycnocline water was estimated as the volume of water divided by the subduction rate for each σθ class: 2–4 years for ESTMW (24.5 < σθ < 25.2 kg m−3), 2 years for the lighter STMW (25.0 < σθ < 25.3 kg m−3), 5–9 years for the denser STMW (25.3 < σθ < 25.6 kg m−3), 10–20 years for the lighter CMW (26.0 < σθ < 26.2 kg m−3), 20–30 years for the middle CMW (26.2 < σθ < 26.3 kg m−3), and 60 years or longer for the denser CMW (26.3 < σθ < 26.6 kg m−3). A comparison of the water volume and subduction rate in potential temperature–salinity (θS) space indicated that the upper permanent pycnocline water (25.0 < σθ < 26.2 kg m−3) was directly maintained by nondiffusive subduction of winter surface water, including STMW and lighter CMW. The lower permanent pycnocline water (26.2 < σθ < 26.6 kg m−3) may be maintained through the subduction of fresher and colder water from the subarctic–subtropical transition region and subsequent mixing with saltier and warmer water. Diagnosis of the potential vorticity (PV) of the subducted water demonstrated that the low PV of STMW was mainly due to the large subduction rate, whereas that of both ESTMW and CMW was due mainly to the small density advection rate (cross-isopycnal flow). Additionally, a relatively large subduction rate probably contributes to the low PV of part of the lighter CMW (ESTMW) formed in the region around 38°N and 170°W (28°N and 145°W), which is characterized by a relatively thick winter mixed layer and an associated mixed-layer front, causing a large lateral induction rate.  相似文献   

16.
The strength of the vertical mixing in the bottom mixed layer near the continental shelf break in the East China Sea was directly measured with the Micro-Scale Profiler (MSP). It has been shown that there is no significant statistical relation between the turbulent energy dissipation and the degree of the stratificationN 2. It seems that the vigorous turbulence occurs not constantly but intermittently in the bottom mixed layer so that a large variation of is found depending on the time. In contrast to , the coefficient of the vertical eddy diffusivityK z is mostly determined byN such thatK z is large in the bottom mixed layer and small in the thermocline. Large value ofK z in the bottom mixed layer is also found in the time series ofK z estimated in terms of Richardson number calculated from the data obtained with electromagnetic current meters. The value ofK z more than 10 cm2s–1 frequently occur in the layer of 20–25 m thick just above the bottom.  相似文献   

17.
The method developed earlier, based on hydration numbers for individual ionic species, has been extended to aqueous systems of four electrolytes, MX, M2Y, NX2, and PX2, with an unhydrated anion (X). The procedure is used to calculate γNa, γMg, γCa, γCl and γSO4 in a mixed electrolyte solution closely resembling seawater. Values of γK and γF in sea-water have also been estimated.  相似文献   

18.
Stable carbon and nitrogen isotopic composition of zooplankton, suspended particulate organic matter (SPOM), and sinking particles collected using sediment traps were measured for samples obtained from the southeastern Bering Sea middle and outer shelf during 1997–1999. The quantity of material collected by the middle shelf sediment trap was greater in both spring and late summer and fall than in early and mid-summer. The δ15N of SPOM, sinking material and zooplankton showed greater inter-annual variability at the middle shelf site (M2) than at the outer shelf site (M3). Zooplankton and sinking organic matter collected by M2 sediment traps became more depleted in 15N from 1997 through 1999, associated with a change from unusually warm to unusually cold conditions. Suspended and sinking organic matter and zooplankton collected from M3 decreased only slightly in δ15N from 1998 to 1999. SPOM, zooplankton, and sediment trap samples collected at M2 were usually enriched in δ15N and δ13C over those from M3. However, in 1999 sediment trap samples from the middle shelf were enriched in 13C over M3 material, but the δ15N of samples from the two sites was similar. The geographic pattern could be explained greater productivity over the middle shelf, associated with either isotopically heavy nitrogen being regenerated from sediments, or with utilization of a greater fraction of the available inorganic nitrogen pool during most years.  相似文献   

19.
This paper assesses the feasibility of applying remote sensing algorithms based on blue/green reflectance ratios to Case 2 waters. Two algorithms from the SeaDAS (NASA) image processing package, OC4v4 for surface chlorophyll concentration, Chl, and K(490) for the attenuation coefficient for downward irradiance at 490 nm, Kd(490), were investigated using an extensive set of observations from the Irish and Celtic Seas. In situ data from a profiling radiometer were used as inputs for the algorithms to avoid uncertainties in atmospheric correction procedures, and direct measurements of Chl and Kd490 and were used for validation purposes. The standard versions of the algorithms performed poorly: OC4v4 generally overestimated Chl (with a very low coefficient of determination), and K(490) progressively underestimated Kd490 for values greater than 0.3 m−1. A two-step procedure for level 2 product generation was therefore devised in which the numerical coefficients of OC4v4 and K(490) were tuned for the two optical water types known to occur most frequently in this region (McKee, D., Cunningham, A., 2006. Identification and characterisation of two optical water types in the Irish Sea from in situ inherent optical properties and seawater constituents. Estuarine, Coastal and Shelf Science 68, 305–316) by statistical regression on a data set of 102 stations from the Irish and Celtic Seas. The water types were distinguished by the magnitude of their normalised water leaving radiance signals at 665 nm, nLw(665), and appropriate versions of the tuned algorithms applied to each water type. When this procedure was tested on an independent data set of 19 stations from the Bristol Channel, Chl values were recovered with an RMS error of 0.36 mg m−3 and Kd(490) values with an RMS error of 0.095 m−1. The identification of water types from water-leaving radiance signals, and the application of band-ratio algorithms tuned for specific water types, may therefore provide a simple means of improving the quality of remote sensing products in optically complex shelf seas.  相似文献   

20.
We examined stable carbon and nitrogen isotopic signatures of 17 fish and 16 invertebrate taxa common to the Newfoundland and Labrador (NL) continental shelf food web. Particular sampling emphasis was placed on Atlantic cod (Gadus morhua) and related prey species (e.g. shrimp, Pandalus borealis, and capelin, Mallotus villosus). We found highly significant (p < 0.0001) differences between near-shore (bays) and offshore (shelf edge) δ15N signatures for cod, ‘other fish’ (pooled) and invertebrates (pooled). In contrast, there were only minor differences in δ13C signatures of ‘other fish’ (p < 0.05) and no difference for cod and invertebrates among the two habitats. We sampled at two times of the year (January and June) and found no systematic effect of season on both δ13C and δ15N in cod, ‘other fish’ and invertebrates. We calculated isotopic fractionation factors for cod from the entire shelf (mixed diet) and for cod with diets composed mainly of capelin or shrimp. These values ranged between 2.2‰ and 3.9‰ for δ15N and −0.4‰ and 0.8‰ for δ13C and, for δ15N, may reflect diet-related differences in bioenergetic status. We discuss potential mechanisms for near-shore versus offshore enrichment of δ15N signatures, and demonstrate the implications of this spatial variation on δ15N-derived trophic position estimates.  相似文献   

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

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