首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
In spite of the fundamental role the Atlantic Meridional Overturning Circulation (AMOC) plays for global climate stability, no direct current measurement of the Denmark Strait Overflow, which is the densest part of the AMOC, has been available until recently that resolve the cross-stream structure at the sill for long periods. Since 1999, an array of bottom-mounted acoustic instruments measuring current velocity and bottom-to-surface acoustic travel times was deployed at the sill. Here, the optimization of the array configuration based on a numerical overflow model is discussed. The simulation proves that more than 80% of the dense water transport variability is captured by two to three acoustic current profilers (ADCPs). The results are compared with time series from ADCPs and Inverted Echo Sounders deployed from 1999 to 2003, confirming that the dense overflow plume can be reliably measured by bottom-mounted instruments and that the overflow is largely geostrophically balanced at the sill.  相似文献   

2.
To describe the exchange of water and sediment through the Venice Lagoon inlets a 3-D hydrodynamic and sediment transport model has been developed and applied to a domain comprising Venice Lagoon and a part of the Adriatic Sea. The model has been validated for both current velocities and suspended particle concentration against direct observations and from observations empirically derived fluxes from upward-looking acoustic Doppler current profiler probes installed inside each inlet. The model provides estimates of the suspended sediment transport in the lower 3 m of the water column that is not detected by acoustic Doppler current profiler sensors. The bedload model prediction has been validated against measured sand transport rates collected by sand traps deployed in the Lido and Chioggia inlets. Results indicate that, in the Lido inlet, 87% of the total load is in suspension, while the rest moves as bedload.  相似文献   

3.
Quantifying sediment flux within rivers is a challenge for many disciplines due, mainly, to difficulties inherent to traditional sediment sampling methods. These methods are operationally complex, high cost, and high risk. Additionally, the resulting data provide a low spatial and temporal resolution estimate of the total sediment flux, which has impeded advances in the understanding of the hydro-geomorphic characteristics of rivers. Acoustic technologies have been recognized as a leading tool for increasing the resolution of sediment data by relating their echo intensity level measurements to suspended sediment. Further effort is required to robustly test and develop these techniques across a wide range of conditions found in natural river systems. This article aims to evaluate the application of acoustic inversion techniques using commercially available, down-looking acoustic Doppler current profilers (ADCPs) in quantifying suspended sediment in a large sand bed river with varying bi-modal particle size distributions, wash load and suspended-sand ratios, and water stages. To achieve this objective, suspended sediment was physically sampled along the Paraná River, Argentina, under various hydro-sedimentological regimes. Two ADCPs emitting different sound frequencies were used to simultaneously profile echo intensity level within the water column. Using the sonar equation, calibrations were determined between suspended-sand concentrations and acoustic backscatter to solve the inverse problem. The study also analyzed the roles played by each term of the sonar equation, such as ADCP frequency, power supply, instrument constants, and particle size distributions typically found in sand bed rivers, on sediment attenuation and backscatter. Calibrations were successfully developed between corrected backscatter and suspended-sand concentrations for all sites and ADCP frequencies, resulting in mean suspended-sand concentration estimates within about 40% of the mean sampled concentrations. Noise values, calculated using the sonar equation and sediment sample characteristics, were fairly constant across evaluations, suggesting that they could be applied to other sand bed rivers. © 2018 John Wiley & Sons, Ltd.  相似文献   

4.
High-resolution current measurements were made in the Northeastern Gulf of Mexico by the Naval Research Laboratory (NRL) as part of its Slope to Shelf Energetics and Exchange Dynamics (SEED) project. The major goal of SEED is to understand the mechanisms that transfer properties across the shelf slope. Fourteen acoustic Doppler current profilers (ADCPs) were deployed just west of the DeSoto Canyon on the shelf and down the slope from May to November, 2004 to measure nearly full water column current profiles. Currents were found more variable on the shelf than on the slope but in the mean strongly tended to follow bathymetry, particularly on the slope. During the SEED time period currents were driven by both local and remote winds, by cyclonic eddies associated with the Loop Current extension and Loop Current rings, by smaller eddies associated with the cyclonic eddies, by frontal meanders or streamers associated with the eddies, and by tropical storms. Currents were highly barotropic, accounting for more than 80% of the eddy kinetic energy (EKE). Current magnitudes generally increased from west to east, towards the DeSoto Canyon. Tropical storms had a relatively minor short-term effect upon mass transports. Cross-shelf transports were much smaller than the along-shelf transports. Onshore transports were stronger on the western side of the array while offshore transports were stronger on the eastern side of the array near the DeSoto Canyon. Offshore transports generally occurred during eastward flow periods, onshore transports during westward flow periods, and both during eddy periods. Mesoscale eddies also provided contributions to cross-shelf exchange. Large scale circulation features could be determined from the first two empirical-orthogonal function (EOF) modes which accounted for 83% of the variance and were strongly related to the integrated wind stress.  相似文献   

5.
We present detailed observations of internally generated turbulence in a sheared, stratified natural flow, as well as an analysis of the external factors leading to its generation and temporal variability. Multi-month time series of vertical profiles of velocity, acoustic backscatter (0.5 Hz), and turbulence parameters were collected with two moored acoustic Doppler current profilers (ADCPs) in the Hudson River estuary, and estuary-long transects of water density were collected 30 times. ADCP backscatter is used for visualization of coherent turbulent structures and evaluation of surface wave biases to the turbulence measurements. Benefits of the continuous long-term turbulence record include our capturing: (1) the seasonality of turbulence due to changing riverflow, (2) hysteresis in stratification and turbulence over the fortnightly cycle of tidal range, and (3) intermittent events such as breaking internal waves. Internal mixing layers (IMLs) are defined as turbulent regions above the logarithmic velocity layer, and the bottom boundary layer (BBL) is defined as the continuously turbulent range of heights above the bed. A cross-correlation analysis reveals how IML and BBL turbulence vary with stratification and external forcing from tidal range, river flow, and winds. Turbulence in both layers is maximal at spring tide and minimal when most stratified, with one exception—IML turbulence at a site with changing channel depth and width is maximal at times of maximum stratification and freshwater input.  相似文献   

6.
Bed diversity in the shallow water environment of Pappas lagoon in Greece   总被引:2,自引:0,他引:2  
<正>Acoustic classification systems and the Sediment Trend Analysis method were used to identify and map the bed diversity in a very shallow(3.0m),coastal lagoon,Pappas lagoon,Western Greece. Analogue acoustic data,collected by means of a 100 kHz side scan sonar system,were digitized and classified into six acoustic classes using recently developed acoustic classification systems (SonarClass,TargAn).By comparing the acoustic classes to ground truth data consisting of sediment grain size and visual inspection of the lagoon-bed,it is demonstrated that the six acoustic classes correlate well with the predominant surface sediment types and vegetation.Thus the spatial distribution of the classes can be considered to represent the spatial pattern of the sedimentary assemblages of the lagoon.The grain size trend analysis identified three dominant sediment pathways and directional trends which could be related to the predominant wind direction,the sediment influx through the inlets and the sediment supply from a small stream in the southern part of the lagoon.The integration of acoustic and sedimentological data together with advanced data processing systems leads to a better understanding of the sedimentary,morphological and biological processes in a shallow lagoon in different spatial and temporal scales and will therefore be beneficial to both sedimentological and biotic-diversity studies.  相似文献   

7.
A hydrodynamic model of the Venice Lagoon and the Adriatic Sea has been developed in order to study the exchanges at the inlets of the Venice Lagoon, a complex morphological area connecting the sea and the lagoon. The model solves the shallow water equations on a spatial domain discretized by a staggered finite element grid. The grid represents the Adriatic Sea and the Venice Lagoon with different spatial resolutions varying from 30 m for the smallest channels of the lagoon to 30  km for the inner areas of the central Adriatic Sea. Data from more than ten tide gauges displaced in the Adriatic Sea have been used in the calibration of the simulated water levels. After the calibration, the tidal wave propagation in the North Adriatic and in the Venice Lagoon is well reproduced by the model. To validate the model results, empirical flux data measured by acoustic Doppler current profiler probes installed inside the inlets of Lido and Malamocco have been used and the exchanges through the three inlets of the Venice Lagoon have been analyzed. The comparison between modeled and measured fluxes at the inlets outlines the efficiency of the model to reproduce both tide- and wind-induced water exchanges between the sea and the lagoon. Even in complex areas, where highly varying resolution is needed, the model is suitable for the simulation of the dominating physical processes.  相似文献   

8.
Observational evidence for estuarine circulation in the German Wadden Sea   总被引:1,自引:0,他引:1  
Observational evidence is presented, which corroborates the hypothesis of the general presence of estuarine circulation in the Wadden Sea as put forward in a previous study (Burchard et al., 2008). Current velocity data from moored ADCPs (in the Hörnum Deep south of Sylt Island, 2002-2009) and ship cruises (in several locations in the German Wadden Sea, 2000-2008) were analysed. As a general result, the vertical current profiles above the benthic boundary layer are usually more homogeneous during flood than during ebb, with a pronounced dependence on the cross-shore horizontal density difference. This tidal asymmetry consequently must lead to a residual outflow of Wadden Sea waters in the upper part and a residual inflow of water in the lower part of the water column, thus giving a generic explanation for the obvious net import of suspended sediments from the German Bight into the Wadden Sea.  相似文献   

9.
Back-barrier tidal flat systems are characterized by basins and inlets through which water is exchanged with the coastal sea by tidal water movements. The hydrographic and morphometric properties at the inlets and in the basins vary considerably, but there is little information available how biogeochemical properties in the water column at these different sites respond to these differences. Therefore, we investigated tidal dynamics of suspended particulate matter (SPM), particulate and dissolved organic carbon (DOC), chlorophyll a, phaeopigments, numbers of particle-associated (PA) and free-living bacteria (FL), bacterial biomass production, and concentrations of dissolved manganese (Mn). Samples were taken at the surface, a mid-depth and 1 m above the bottom at a fixed station at the inlet and in the basin of the Spiekeroog back-barrier tidal flat system in the German Wadden Sea. Five tidal cycles representative for typical seasonal situations, January (winter), April and May (late spring bloom), July (summer), and November (late fall) were studied in 2005 and 2006. In July, processes related to phytoplankton dynamics and bacterial decomposition were much more enhanced in the basin, whereas in April, these processes were enhanced at the inlet but were particularly low within the basin itself. The low values within the basin were a result of the settled phytoplankton spring bloom and represent a rather short period at the decline of this bloom. In November and January, differences were much less pronounced than during the growing season and restricted mainly to SPM and PA bacteria, exhibiting higher values in the basin. FL bacteria, DOC, and dissolved Mn exhibited different patterns and much less differences between the two stations, indicating that biogeochemical processes in the dissolved phase were controlled by different factors than PA biogeochemical processes. These differences reflect the retentive properties of the basin for particles and PA biogeochemical processes, particularly during the growing season, and in general emphasize the high productivity of back-barrier tidal flat systems.  相似文献   

10.
Tidal marshes form at the confluence between estuarine and marine environments where tidal movement regulates their developmental processes. Here, we investigate how the interplay between tides, channel morphology, and vegetation affect sediment dynamics in a low energy tidal marsh at the Paul S. Sarbanes Ecosystem Restoration Project at Poplar Island. Poplar Island is an active restoration site where fine-grained material dredged from navigation channels in the upper Chesapeake Bay are being used to restore remote tidal marsh habitat toward the middle bay (Maryland, USA). Tidal currents were measured over multiple tidal cycles in the inlets and tidal creeks of one marsh at Poplar Island, Cell 1B, using Acoustic Doppler Current Profilers (ADCP) to estimate water fluxes throughout the marsh complex. Sediment fluxes were estimated using acoustic backscatter recorded by ADCPs and validated against total suspended solid measurements taken on site. A high-resolution geomorphic survey was conducted to capture channel cross sections and tidal marsh morphology. We integrated simple numerical models built in Delft3d with empirical observations to identify which eco-geomorphological factors influence sediment distribution in various channel configurations with differing vegetative characteristics. Channel morphology influences flood-ebb dominance in marshes, where deep, narrow channels promote high tidal velocities and incision, increasing sediment suspension and reducing resilience in marshes at Poplar Island. Our numerical models suggest that accurately modelling plant phenology is vital for estimating sediment accretion rates. In-situ observations indicate that Poplar Island marshes are experiencing erosion typical for many Chesapeake Bay islands. Peak periods of sediment suspension frequently coincide with the largest outflows of water during ebb tides resulting in large sediment deficits. Ebb dominance (net sediment export) in tidal marshes is likely amplified by sea-level rise and may lower marsh resilience. We couple field observations with numerical models to understand how tidal marsh morphodynamics contribute to marsh resilience. © 2019 John Wiley & Sons, Ltd.  相似文献   

11.
史小丽  范帆  张民  阳振  陈开宁 《湖泊科学》2020,32(5):1446-1453
2018年10月-2019年10月对巢湖西湖心水体浮游藻类群落结构以及水体和底泥蓝藻生物量进行了月度调查.结果表明:巢湖西湖心浮游藻类的主要优势种属为微囊藻属、席藻属、十字藻、卵形隐藻和鱼腥藻属.蓝藻优势种属在5-10月为微囊藻属,11-12月为鱼腥藻属,1-4月为席藻属.巢湖水体和底泥蓝藻生物量峰值分别出现在9月和2月,水体蓝藻的衰亡下沉会导致底泥蓝藻生物量的上升.巢湖蓝藻主要分布在水体,底泥蓝藻生物量相对较低,单位面积水柱与底泥蓝藻生物量6月的比值大于100,在11-3月相对较低,最低值小于2.底泥蓝藻主要分布在底泥表层0~2 cm.通过安装原位捕获器,监测了蓝藻在西巢湖湖心水柱和底泥中的垂直迁移过程和通量.结果表明:11月和2月蓝藻有明显从水柱向底泥迁移的过程;底泥蓝藻全年向水体的静态迁移量都很低,而动态迁移在11月和6月出现两个峰值,主要受底泥蓝藻生物量和再悬浮的影响.本研究结果表明削减巢湖西湖心底泥种源的最佳时期为10月至来年2月,但是由于底泥蓝藻生物量远远小于水柱蓝藻生物量,底泥蓝藻向水体复苏迁移的通量也较低,即使削减了底泥种源,也不能有效降低水体蓝藻生物量.  相似文献   

12.
Continuous monitoring of bed shear stress in large river systems may serve to better estimate alluvial sediment transport to the coastal ocean.Here we explore the possibility of using a horizontally deployed acoustic Doppler current profiler(ADCP) to monitor bed shear stress,applying a prescribed boundary layer model,previously used for discharge estimation.The model parameters include the local roughness length and a dip correction factor to account for sidewall effects.Both these parameters depend on river stage and on the position in the cross-section, and were estimated from shipborne ADCP data.We applied the calibrated boundary layer model to obtain bed shear stress estimates over the measuring range of the HADCP.To validate the results,co-located coupled ADCPs were used to infer bed shear stress,both from Reynolds stress profiles and from mean velocity profiles. From HADCP data collected over a period of 1.5 years,a time series of width profiles of bed shear stress was obtained for a tidal reach of the Mahakam River,East Kalimantan,Indonesia.A smaller dataset covering 25 hours was used for comparison with results from the coupled ADCPs.The bed shear stress estimates derived from Reynolds stress profiles appeared to be strongly affected by local effects causing upflow and downflow,which are not included in the boundary layer model used to derive bed shear stress with the horizontal ADCP.Bed shear stresses from the coupled ADCP are representative of a much more localized flow,while those derived with the horizontal ADCP resemble the net effect of the flow over larger scales.Bed shear stresses obtained from mean velocity profiles from the coupled ADCPs show a good agreement between the two methods,and highlight the robustness of the method to uncertainty in the estimates of the roughness length.  相似文献   

13.
The present study investigated direct and indirect methods using optical and acoustic instruments for the acquisition of information required to estimate the concentration of suspended particulate matter (SPM) in Vitória Bay, a shallow estuarine system in SE Brazil. The aim was to calibrate and compare the use of different instruments (OBS, ADP, and ADV) to estimate SPM concentration in the water column and near the bed. Concentration was determined by correlating filtered water samples with the optical and acoustic measurements. In general, the methodology proved tenable for the chosen shallow estuarine environment with low SPM concentration (<60?mg/L). Pearson’s coefficient ranged from 0.75 to 0.85, when correlating measurements taken at three sampled depths. Differences in the correlation coefficient values showed that calibration at three depths (near the bed, mid-water column, and near the surface) was more effective than for surface samples alone, even in shallow (???m deep) water. When calibration was attempted for concentration in the entire water column with samples at just one elevation, the correlation value was very low, thus increasing the error in estimating the SPM concentration.  相似文献   

14.
This paper illustrates how the acoustic Doppler current profiler (ADCP) and single-beam echo-sounder (SBES) recordings can be used for the calibration of existing software to assist in generalizing the morphodynamic processes in large rivers at key sites such as bifi.trcations and confluences. Calibration of the MIKE21C numerical model by the Danish Hydraulic Institute at the 25-km-long reach of Lower Paran~ near Rosario (Argentina) is presented. This reach includes two downstream confluences and two bifurcations. The model simulates a 2-D depth-averaged flow velocity and the related sediment fluxes to predict the bifurcation morphodynamics that affects the Paranh waterway. To investigate the river channel bathymetry, roughness, flow discharge allocation at bifurcations, suspended sediment concentration and grain size distributions, several instruments were used. These instruments included two ADCPs by Teledyne RDI working at frequencies of 600 and 1,200 kHz, a Sontek ADCP working at a frequency of 1,000 kHz and a SBES. The method to assess suspended sediment concentration and grain size distributions has been previously described. This paper focuses primarily on investigating dune morphology (by means of SBES depth measurements) and friction velocity (by means of ADCP profiling) to determine the river channel bed-roughness. The 2-D model results agree with observed values of bed-roughness, flow velocity and suspended sediment concentration distributions at the investigated sections, known data of water slope and total load of bed sediment are in good agreement with model results.  相似文献   

15.
Fourteen acoustic Doppler current profilers (ADCPs) were deployed on the shelf and slope for 1 year just west of the DeSoto Canyon in the Northeastern Gulf of Mexico by the Naval Research Laboratory (NRL) as part of its Slope to Shelf Energetics and Exchange Dynamics (SEED) project. The winter and spring observations are discussed here in regards to the low-frequency current variability and its relation to wind and eddy forcing. Empirical orthogonal function (EOF) analyses showed that two modes described most of the current variability. Wind-forced variability of the along-shelf flow was the main contributor in Mode 1 while eddies contributed much of the variability in Mode 2. Wind-stress controlled currents on the shelf and slope at time scales of about a week. On longer time scales, variations in the currents on both the outer shelf and slope appear to be related to seasonal variations in the time-cumulated wind stress curl. Winds were dominant in driving the along-shelf transports, particularly along the slope. However, the effective wind stress component was found to be aligned with the west Florida shelf direction rather than the local shelf direction. Eddy intrusions, which were more numerous in winter and spring than in summer and fall, and winds were found to contribute significantly to cross-shelf exchange processes.  相似文献   

16.
In situ observations and numerical simulations of turbulence are essential to understanding vertical mixing processes and their dynamical controls on both physical and biogeochemical processes in coastal embayments. Using in situ data collected by bottom-mounted acoustic Doppler current profilers(ADCPs) and a free-falling microstructure profiler, as well as numerical simulations with a second-moment turbulence closure model, we studied turbulence and mixing in the Xiamen Bay, a freshwater-influenced tidal bay located at the west coast of the Taiwan Strait. Dynamically, the bay is driven predominantly by the M2 tide, and it is under a significant influence of the freshwater discharged from the Jiulong River. It is found that turbulence quantities such as the production and dissipation rates of the turbulent kinetic energy(TKE) were all subject to significant tidal variations, with a pronounced ebb-flood asymmetry. Turbulence was stronger during flood than ebb. During the flooding period, the whole water column was nearly well mixed with the depth-averaged TKE production rate and vertical eddy viscosity being up to 5?10?6 W kg?1 and 2?10?2 m2 s?1, respectively. In contrast, during the ebb strong turbulence was confined only to a 5?8 m thick bottom boundary layer, where turbulence intensity generally decreases with distance from the seafloor. Diagnosis of the potential energy anomaly showed that the ebb-flood asymmetry in turbulent dissipation and mixing was due mainly to tidal straining process as a result of the interaction between vertically shared tidal currents and horizontal density gradients. The role of vertical mixing in generating the asymmetry was secondary. A direct comparison of the modeled and observed turbulence quantities confirmed the applicability of the second-moment turbulence closure scheme in modeling turbulent processes in this weakly stratified tidally energetic environment, but also pointed out the necessity of further refinements of the model.  相似文献   

17.
18.
声子波及其在地震波资料分解中的应用   总被引:3,自引:0,他引:3       下载免费PDF全文
声子波是由声波波动方程的解构成的一种物理子波,如果不考虑吸收和散射,声子波的传播是相当简单的;相反地,数学子波的传播即使在均匀介质中也是极其复杂的.作为波动方程的解,声子波比一般的数学子波更能有效地应用于复杂声波和地震波的分解和分析.本文从Kaiser的声子波理论出发,给出了通过分别引入点源波形的复时间函数和点源虚时间坐标来构成声子波的两种解释,并对点源模型的合成地震图和实际复杂模型的地震波资料进行了时-空域的声子波变换,说明了声子波应用于地震波资料分解的有效性.  相似文献   

19.
This paper deals with the interaction and small-scale processes occurring around the inlets that connect the Venice Lagoon with the Northern Adriatic Sea. In a previous paper, barotropic processes have been investigated, whereas here, the focus is on the baroclinic processes. The hydrodynamics of the area are studied by means of a 3D shallow water hydrodynamic finite-element model, suitable to describe areas of complex morphology such as the coasts and the interaction channels. This is the first work that models the 3D interaction between the Venice Lagoon and the Adriatic Sea. Three different sets of simulations have been carried out to identify the physics behind the small-scale processes and the influence of the main forcings on the study area. The first imposes different idealized forcings, such as tides, wind, and river runoff. The vorticity maps of the first two layers show the predominance of wind forcing in the coastal area and tidal forcing in the three inlets of the Lagoon. Bora wind acts homogeneously, increasing the littoral currents, while Sirocco wind mainly impacts near Chioggia inlet, with a coastal current reversal, inducing its detachment offshore. Freshwater patterns are present along the coast, near the river mouths. Rivers do not directly influence the circulation close to the coast in front of the Venice Lagoon, except for the area near Chioggia inlet, where the Brenta river action can be seen. The second set of simulations deals with a sensitivity analysis to define the importance of the advection and of the baroclinic pressure gradient terms in the creation of persistent structures, such as small-scale coastal vortices seen along the littoral very close to the inlets. This analysis shows how advection is the main physical process responsible for the persistence of the positive vorticity structures close to the coast between the inlets, while the negative vorticity structures, also seen by the HF Radar, are due to the baroclinic-advective interaction. Finally, a real case, year 2004, has been simulated both to validate the model with observations and to identify the occurrence during the year of the characteristic hydrodynamic features attributable to the main forcings. The action of Bora wind characterizes the surface current patterns of February and November 2004, while Sirocco influences the month of May 2004. During periods of weak wind, the model reproduces the small-scale vortical structures close to the littoral.  相似文献   

20.

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

  相似文献   

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

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