首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The increase of salt intrusion in recent years in the Modaomen Estuary, one of the estuaries of the Pearl River Delta in China, has threatened the freshwater supply in the surrounding regions, especially the cities of Zhongshan, Zhuhai in Guangdong Province and Macau. A numerical modeling system using nested grids was developed to investigate the salt transport mechanisms and the response of salt intrusion to changes in river discharge and tidal mixing. The steady shear transport induced by estuarine circulation reaches maximum and minimum, respectively, during neap and spring tides, while the tidal oscillatory transport shows an opposite pattern. The net transport is landward during neap tides and seaward during spring tides. The salt intrusion length responding to constant river discharges generally follows a power law of ?0.49. The dependence of salt intrusion on tidal velocity is less than that predicted by theoretical models for exchange flow dominated estuaries. The response of salt intrusion to change in tidal velocity depends largely on river discharge. When river flow increases, the impact of tidal velocity increases and the phase lag of response time decreases. The asymmetries of salt intrusion responding to increasing and decreasing river discharge (tidal velocity) are observed in the estuary.  相似文献   

2.
Estuaries, commonly, are densely populated areas serving the needs of the inhabitants in multiple ways. Often the interests are conflicting and decisions need to be made by the local managers. Intake of fresh water for consumption, agricultural purposes or use by industries may take place within a region not far landward of the limit of salt intrusion. Human interventions (e.g. deepening of the navigation channels) or climate changes (sea level rise, reduction of the river discharge) can bring these intake locations within the reach of saline or brackish water and consequently endanger their function. To support policy and managerial decisions, a profound knowledge of processes associated with the salinity structure in estuaries is required. Although nowadays advanced numerical three-dimensional models are available that are able to cope with the complexity of the physics there is still a need for relatively simple tools for quick-scan actions in a pre-phase of a project or for instructive purposes. The analytical model described in this paper may serve these needs. It computes the maximum salinity distribution using the dispersion coefficient in the mouth as the only model parameter. The model has been calibrated using observational data in a large number of estuaries and experimental data in a tidal flume. The dispersion coefficient was successfully related to geometric and hydrodynamic parameters resulting in an expression that can be used for convergent estuaries as well as prismatic channels, see Eqs. 25a and 25b. Application of the model in a predictive mode showed its promising capabilities. Comparison with three-dimensional numerical models indicates that the channel geometry in the estuary mouth largely influences dispersive processes. The analytical model for salt intrusion may be used in combination with the analytical model for tidal propagation in convergent estuaries and tidal channels by Van Rijn (part I). In this way, input is obtained on the tidal velocity amplitude and the Chézy roughness following calibration of this model on tidal amplitudes along the estuary.  相似文献   

3.
We present an analytical model to decompose complex along-channel and transverse residual flows into components induced by individual mechanisms. The model describes the transverse distribution of residual flows in tidally dominated estuaries. Scaling and perturbation techniques are used to obtain analytical solutions for residual flows over arbitrary across-channel bed profiles. The flows are induced by horizontal density gradients, tidal rectification processes, river discharge, wind, channel curvature and the earth's rotation. These rectification processes induce residual flows that are up-estuary to the right and down-estuary to the left of an estuarine channel (looking up-estuary in the northern hemisphere). The tidal rectification processes fundamentally change the transverse structure of along-channel residual flows in many tidal estuaries, as these processes cause the flows to be internally asymmetric about the mid-axis of the channel for relatively large tidal velocities, steep channels or narrow estuaries. In addition, velocity scales are derived from the analytical solutions to estimate the relative importance of the various residual flow mechanisms from estuarine parameters. A case study of a transect across the Upper Chesapeake Bay showed that important features of the residual flow observed in that transect are reproduced and explained by the analytical model. The velocity scales were able to identify the relevant residual flow mechanisms as well. The tidal rectification processes considered here result from advection of along-channel tidal momentum by Coriolis-induced transverse tidal currents.  相似文献   

4.
Bifurcations in tidally influenced deltas distribute river discharge over downstream channels, asserting a strong control over terrestrial runoff to the coastal ocean. Whereas the mechanics of river bifurcations is well-understood, junctions in tidal channels have received comparatively little attention in the literature. This paper aims to quantify the tidal impact on subtidal discharge distribution at the bifurcations in the Mahakam Delta, East Kalimantan, Indonesia. The Mahakam Delta is a regular fan-shaped delta, composed of a quasi-symmetric network of rectilinear distributaries and sinuous tidal channels. A depth-averaged version of the unstructured-mesh, finite-element model second-generation Louvain-la-Neuve Ice-ocean Model has been used to simulate the hydrodynamics driven by river discharge and tides in the delta channel network. The model was forced with tides at open sea boundaries and with measured and modeled river discharge at upstream locations. Calibration was performed with water level time series and flow measurements, both spanning a simulation period. Validation was performed by comparing the model results with discharge measurements at the two principal bifurcations in the delta. Results indicate that within 10 to 15 km from the delta apex, the tides alter the river discharge division by about 10% in all bifurcations. The tidal impact increases seaward, with a maximum value of the order of 30%. In general, the effect of tides is to hamper the discharge division that would occur in the case without tides.  相似文献   

5.
In many large estuaries there are significant variations in flow conditions due to the interaction between tide (with spring–neap changes) and river discharge (with wet–dry seasons), which is key to understanding the evolution of the morphology and the resultant equilibrium state. To explore whether there exists an equilibrium state, and what might control such a state in such a dynamic environment, both numerical and analytical methods have been used to investigate the relative importance of tide and river contributions to the work done locally and globally over a wide range of discharge conditions in the Yangtze estuary. In particular, we have quantified the contributions from the tidal flow, the river flow and the tide–river interaction in terms of energy and its dissipation under different river discharge conditions. Model results suggest that there is a state of minimum tidal work for the case representing the wet season, when river and tide are doing uniform work locally and minimum work globally, within the bi‐directional tidal reach for tide and along the whole estuary for river. We also observe that the system is not optimized for other conditions (peak discharge and low flows during the dry season), but the system would tend to do the minimum work possible given the constraints on the system (e.g. imposed forcing conditions and available sediment supply). Results, therefore, are consistent with the use of these two energetic optimization principles, and the proposed method could be applicable to other alluvial estuaries. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
Analytical solutions of the momentum and energy equations for tidal flow are studied. Analytical solutions are well known for prismatic channels but are less well known for converging channels. As most estuaries have a planform with converging channels, the attention in this paper is fully focused on converging tidal channels. It will be shown that the tidal range along converging channels can be described by relatively simple expressions solving the energy and momentum equations (new approaches). The semi-analytical solution of the energy equation includes quadratic (nonlinear) bottom friction. The analytical solution of the continuity and momentum equations is only possible for linearized bottom friction. The linearized analytical solution is presented for sinusoidal tidal waves with and without reflection in strongly convergent (funnel type) channels. Using these approaches, simple and powerful tools (spreadsheet models) for tidal analysis of amplified and damped tidal wave propagation in converging estuaries have been developed. The analytical solutions are compared with the results of numerical solutions and with measured data of the Western Scheldt Estuary in the Netherlands, the Hooghly Estuary in India and the Delaware Estuary in the USA. The analytical solutions show surprisingly good agreement with measured tidal ranges in these large-scale tidal systems. Convergence is found to be dominant in long and deep-converging channels resulting in an amplified tidal range, whereas bottom friction is generally dominant in shallow converging channels resulting in a damped tidal range. Reflection in closed-end channels is important in the most landward 1/3 length of the total channel length. In strongly convergent channels with a single forward propagating tidal wave, there is a phase lead of the horizontal and vertical tide close to 90o, mimicking a standing wave system (apparent standing wave).  相似文献   

7.
Salinity difference between terrestrial river discharge and oceanic tidal water plays a role in modifying the local flow field and, as a consequence, estuarine morphodynamics. Although widely recognized, recent numerical studies exploring the long-term morphological evolution of river-influenced estuaries with two-dimensional, depth-averaged models have mostly neglected salinity. Using a three-dimensional morphodynamic model, we aim to gain more insight into the effect of salinity on the morphodynamics of fluvio-deltaic systems. Model results indicate that the resultant estuarine morphology established after 600 years differs remarkably when a salinity gradient is included. A fan-shaped river-mouth delta exhibits less seaward expansion and is cut through by narrower channels when salinity is included. The inclusion of salinity tends to generate estuarine circulation, which favours landward sediment transport and hence limits the growth of the delta while enhancing the development of intertidal areas. The formation of deltaic channel–shoal patterns resulting from morphodynamic evolution tends to strengthen salinity stratification, which is characterized by an increased gradient Richardson number. The direction of the depth-averaged residual sediment transport over a tide may be opposite to the direction of residual velocity, indicating the significant influence of baroclinic effects on the net sediment transport direction (and hence morphological change). The effect of salinity on morphological evolution becomes less profound when the strength of tidal or fluvial forcing is dominant over the other. The effects of sediment type and flocculation, which are particularly important when salinity gradients are present, are also discussed. Overall, this study highlights that neglecting salinity to simulate long-term estuarine morphodynamics requires more careful justification, particularly when the environment is characterized by fine sediment types (favouring suspended transport), and relatively large river discharge and estuarine depth (favouring baroclinic effects). © 2020 John Wiley & Sons, Ltd.  相似文献   

8.
The bed of estuaries is often characterized by ripples and dunes of varying size. Whereas smaller bedforms adapt their morphological shape to the oscillating tidal currents, large compound dunes (here: asymmetric tidal dunes) remain stable for periods longer than a tidal cycle. Bedforms constitute a form roughness, that is, hydraulic flow resistance, which has a large-scale effect on tidal asymmetry and, hence, on hydrodynamics, sediment transport, and morphodynamics of estuaries and coastal seas. Flow separation behind the dune crest and recirculation on the steep downstream side result in turbulence and energy loss. Since the energy dissipation can be related to the dune lee slope angle, asymmetric dune shapes induce variable flow resistance during ebb and flood phases. Here, a noncalibrated numerical model has been applied to analyze the large-scale effect of symmetric and asymmetric dune shapes on estuarine tidal asymmetry evaluated by residual bed load sediment transport at the Weser estuary, Germany. Scenario simulations were performed with parameterized bed roughness of symmetric and asymmetric dune shapes and without dune roughness. The spatiotemporal interaction of distinct dune shapes with the main drivers of estuarine sediment and morphodynamics, that is, river discharge and tidal energy, is shown to be complex but substantial. The contrasting effects of flood- and ebb-oriented asymmetric dunes on residual bed load transport rates and directions are estimated to be of a similar importance as the controls of seasonal changes of discharge on these net sediment fluxes at the Lower Weser estuary. This corroborates the need to consider dune-induced directional bed roughness in numerical models of estuarine and tidal environments.  相似文献   

9.
An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s − 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.  相似文献   

10.
Mainstem–floodplain material exchange in the tidal freshwater reach of major rivers may lead to significant sequestration of riverine sediment, but this zone remains understudied compared to adjacent fluvial and marine environments. This knowledge gap prompts investigation of floodplain-incising tidal channels found along the banks of tidal rivers and their role in facilitating water and suspended-sediment fluxes between mainstem and floodplain. To evaluate this role, and how it evolves along the tidal river and with time, we measured water level, flow velocity, temperature, and suspended-sediment concentration (SSC) in four tidal channels along the tidal Amazon River, Brazil. Eleven deployments were made during low, rising, high, and falling seasonal Amazon discharge. Generally, channels export high-SSC water from the mainstem to the tidal floodplain on flood tides and transfer low-SSC water back to the mainstem on ebbs. Along the length of the tidal river, the interaction between tidal and seasonal water-level variations and channel–floodplain morphology is a primary control on tidal-channel sediment dynamics. Close to the river mouth, where tides are large, this interaction produces transient flow features and current-induced sediment resuspension, but the importance of these processes decreases with distance upstream. Although the magnitude of the exchange of water and sediment between mainstem and floodplain via tidal channels is a small percentage of the total mainstem discharge in this large tidal-river system, tidal channels are important conduits for material flux between these two environments. This flux is critical to resisting floodplain submergence during times of rising sea level. © 2019 John Wiley & Sons, Ltd.  相似文献   

11.
The influence of river discharge on tidal damping in alluvial estuaries   总被引:2,自引:0,他引:2  
The tidal range, the difference between high water level and low water level, along an alluvial estuary can be described by Savenije's analytical equation [Journal of Hydrology 243 (2001) 205-215] analytical equation. This equation was derived from the complete St Venant equations in a Lagrangian reference frame. In the derivation of that equation the effect of river discharge was disregarded. Measurements in the Schelde Estuary show that this assumption is only valid in the lower part of the estuary, but that in the upper part the river discharge has an influence on tidal damping. In the downstream part of the estuary, where the cross-sectional area is large compared to the cross-sectional area of the river, it is correct to disregard the river discharge. However, in the upstream part of the estuary, where the cross-sectional area approaches that of the river, the fresh water discharge gains importance over the tidal flow and affects the tidal range. In this paper, the derivation of the analytical equation is expanded to include the effect of the river discharge. It is demonstrated that the river discharge can have a considerable influence on tidal damping in the upper reach of the estuary. The river discharge affects tidal damping primarily through friction. A critical point along the estuary is the point where there is a single slack, upstream of which the fresh water velocity is larger than the tidal velocity. The location of this point varies with the river discharge. From that point onwards the effect of river discharge on damping is dominant. It is the point where the estuary becomes primarily of riverine character.  相似文献   

12.
Numerical modelling of morphodynamics—Vilaine Estuary   总被引:1,自引:0,他引:1  
The main objective of this paper is to develop a method to simulate long-term morphodynamics of estuaries dominated by fine sediments, which are subject to both tidal flow and meteorologically induced variations in freshwater run-off and wave conditions. The method is tested on the Vilaine Estuary located in South Brittany, France. The estuary is subject to a meso–macrotidal regime. The semi-diurnal tidal range varies from around 2.5 to 5 m at neap and spring, respectively. The freshwater input is controlled by a dam located approximately 8 km from the mouth of the estuary. Sediments are characterised as mostly fines, but more sandy areas are also found. The morphology of the estuary is highly influenced by the dam. It is very dynamic and changes in a complicated manner with the run-off from the dam, the tide and the wave forcing at the mouth of the estuary. Extensive hydrodynamic and sediment field data have been collected in the past and provide a solid scientific basis for studying the estuary. Based on a conceptual understanding of the morphodynamics, a numerical morphological model with coupled hydrodynamic, surface wave and sediment transport models is formulated. The numerical models are calibrated to reproduce sediment concentrations, tidal flat altimetry and overall sediment fluxes. Scaling factors are applied to a reference year to form quasi-realistic hydrodynamic forcing and river run-off, which allow for the simulations to be extended to other years. The simulation results are compared with observed bathymetric changes in the estuary during the period 1998–2005. The models and scaling factors are applied to predict the morphological development over a time scale of up to 10 years. The influence of the initial conditions and the sequence of external hydrodynamic forcing, with respect to the morphodynamic response of the estuary, are discussed.  相似文献   

13.
Freshwater inflow is central to the definition of estuaries and if we lose control of the quantity of freshwater flow or discharge (including seasonal timing) to estuaries, then freshwater water quality has the potential to become a moot issue in estuarine ecosystems (Definition of estuaries: estuaries (aestus = tide) are physico-chemically, geomorphically, and biotically diverse ecosystems. Although numerous definitions of estuaries exist, we prefer the following: an estuary is a partially enclosed coastal water body in which freshwater runoff, often seasonally and episodically pulsed, dilutes salty ocean water and the biotic structure is influenced by dynamic tidal action and associated salinity gradients and reef building organisms and wetlands influence development and evolution of ecological structure and function (see for expanded definition)).  相似文献   

14.
15.
Man-induced regime shifts in small estuaries—II: a comparison of rivers   总被引:3,自引:2,他引:1  
This is Part II of two papers on man-induced regime shifts in small, narrow, and converging estuaries, with focus on the interaction between effective hydraulic drag, fine sediment import, and tidal amplification, induced by river engineering works, e.g., narrowing and deepening. Paper I describes a simple linear analytical model for the tidal movement in narrow, converging estuaries and a conceptual model on the response of tidal rivers to river engineering works. It is argued that such engineering works may set in motion a snowball effect bringing the river into an alternative steady state. Part II analyses the historic development in tidal range in four rivers, e.g., the Elbe, Ems, Loire, and Scheldt, all in northwest Europe; data are available for many decades, up to a century. We use the analytical model derived in Part I, showing that the effective hydraulic drag in the Ems and Loire has decreased considerably over time, as anticipated in Part I. We did not find evidence that the Upper Sea Scheldt is close to its tipping point towards hyperturbid conditions, but risks have been identified. In the Elbe, tidal reflections against the profound step in bed level around Hamburg seem to have affected the tidal evolution in the last decades. It is emphasized that the conceptual picture sketched in these papers is still hypothetical and needs to be validated, for instance through hind-cast modeling of the evolution of these rivers. This will not be an easy task, as historical data for a proper calibration of the models required are scarce.  相似文献   

16.
Abstract

Little is known about the salt intrusion behaviour in Malaysian estuaries. Study of salt intrusion generally requires large amounts of data, especially if 2-D or 3-D numerical models are used; thus, in data-poor environments, 1-D analytical models are more appropriate. A fully analytical 1-D salt intrusion model, which is simple to implement and requires minimal data, was tested in six previously unsurveyed Malaysian estuaries (Kurau, Perak, Bernam, Selangor, Muar and Endau). The required data can be collected during a single day of observations. Site measurements were conducted during the dry season (June–August 2012 and February–March 2013) near spring tide. Data on cross-sections (by echo-sounding), water levels (by pressure loggers) and salinity (by moving boat) were collected as model input. A good fit was demonstrated between the simulated and observed salinity distribution for all six estuaries. Additionally, the two calibration parameters (the Van der Burgh coefficient and the boundary condition for the dispersion) were compared with the existing predictive equations. Since gauging stations were only present in some nested catchments in the drainage basins, the river discharge had to be up-scaled to represent the total discharge contribution of the catchments. However, the correspondence between the calibration coefficients and the predictive equations was good, particularly in view of the uncertainty in the river discharge data used. This confirms that the predictive salt intrusion model is valid for the cases studied in Malaysia. The model provides a reliable, predictive tool, which the water authority of Malaysia can use for making decisions on water abstraction or dredging.
Editor D. Koutsoyiannis; Associate editor A. Fiori  相似文献   

17.

This is Part II of two papers on man-induced regime shifts in small, narrow, and converging estuaries, with focus on the interaction between effective hydraulic drag, fine sediment import, and tidal amplification, induced by river engineering works, e.g., narrowing and deepening. Paper I describes a simple linear analytical model for the tidal movement in narrow, converging estuaries and a conceptual model on the response of tidal rivers to river engineering works. It is argued that such engineering works may set in motion a snowball effect bringing the river into an alternative steady state. Part II analyses the historic development in tidal range in four rivers, e.g., the Elbe, Ems, Loire, and Scheldt, all in northwest Europe; data are available for many decades, up to a century. We use the analytical model derived in Part I, showing that the effective hydraulic drag in the Ems and Loire has decreased considerably over time, as anticipated in Part I. We did not find evidence that the Upper Sea Scheldt is close to its tipping point towards hyperturbid conditions, but risks have been identified. In the Elbe, tidal reflections against the profound step in bed level around Hamburg seem to have affected the tidal evolution in the last decades. It is emphasized that the conceptual picture sketched in these papers is still hypothetical and needs to be validated, for instance through hind-cast modeling of the evolution of these rivers. This will not be an easy task, as historical data for a proper calibration of the models required are scarce.

  相似文献   

18.
Seasonal variation in water column conditions in the upper Gulf of Thailand (UGoT) was analyzed by considering four major factors including surface heat flux, freshwater discharge, tidal and wind stirrings. The coincidence of surface heat loss, low river discharge and strong wind resulted in vertical well-mixing in December. Strong stratification developed in September and October due to large river discharge and moderate heat flux. Strong surface heating in April and May has a potential to generate strong stratification, although not as large as that in September and October due to low river discharge. Although no factors are prominent during January and March, and June and August, weak to moderate stratification results, because the influences of river discharge and surface heating are still larger than those of tidal and wind stirrings. The results of water column analysis based on monthly average data agree well with analyses derived from cruise data in the same months. Most analytical results correspond to the distributions of temperature and salinity from field observations. Disagreement, however, was found in December 2003 (cruise CU-2) when stratification in some small regions occurs in the distribution of water properties, but the water column analysis suggests vertical well-mixing. This phenomenon is triggered by non-uniform distribution of freshwater over UGoT, which is related to river discharge, monsoonal wind and current. Compared to a previous study regarding surface chlorophyll dynamics, water column conditions may be used to explain the occurrence of phytoplankton bloom in this region.  相似文献   

19.
River-ocean coupled models are described for the evaluation of the interaction between river discharge and surge development along the Orissa coast of India. The models are used to study the effect of fresh water discharge from the Mahanadi River on the surge response along the Orissa coast due to the October 1999 super cyclone which led to severe flooding of the coastal and delta regions of Orissa. The so-called 1999 Paradip cyclone was one of the most severe cyclones; causing extensive damage to property and loss of lives. The present study emphasizes the impact of the Mahanadi River on overall surge development along the Orissa coast. Therefore, we have developed a location specific fine resolution model for the Orissa coast and coupled it with a one–dimensional river model. The numerical experiments are carried out, both with and without inclusion of fresh water discharge from the river. The bathymetry for the model has been taken from the naval hydrographic charts extending from the south of Orissa to the south of west Bengal. A simple drying scheme has also been included in the model in order to avoid the exposure of land near the coast due to strong negative sea-surface elevations. The simulations with river-ocean coupled models show that the discharge of fresh water carried by the river may modify the surge height in the Bay, especially in the western Bay of Bengal where one of the largest river systems of the east coast of India, the Mahanadi River, joins with the Bay of Bengal. Another dynamic effect of this inlet is the potentially deep inland penetration of the surge originating in the Bay. The model results are in good agreement with the available observations/estimates.  相似文献   

20.
Tidal propagation in estuaries is affected by friction and fresh water discharge, besides changes in the depth and morphology of the channel. Main distortions imply variations in the mean water level and asymmetry. Tidal asymmetry can be important as a mechanism for sediment accumulation and turbidity maximum formation in estuaries, while mean water level changes can affect navigation depths. Data from several gauges stations from the Amazon estuary and the adjacent coast were analyzed and a 2DH hydrodynamic model was configured in a domain covering the continental shelf up to the last section of the river where the tidal signature is observed. Based on data, theoretical and numerical results, the various influences in the generation of estuarine harmonics are presented, including that of fresh water discharge. It is shown that the main overtide, M4, derived from the most important astronomic component in the Amazon estuary, M2, is responsible for the tidal wave asymmetry. This harmonic has its maximum amplitude at the mouth, where minimum depths are found, and then decreases while tide propagates inside the estuary. Also, the numerical results show that the discharge does not affect water level asymmetry; however, the Amazon river discharge plays an important role in the behavior of the horizontal tide. The main compound tide in Amazon estuary, Msf, generated from the combination of the M2 and S2, can be strong enough to provoke neap low waters lower than spring ones. The results show this component increasing while going upstream in the estuary, reaching a maximum and then slightly decaying.  相似文献   

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

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