Next Article in Journal
Extracting Typhoon Disaster Information from VGI Based on Machine Learning
Previous Article in Journal
Application of Sentinel-2 Multispectral Data for Habitat Mapping of Pacific Islands: Palau Republic (Micronesia, Pacific Ocean)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contribution of Pumping Action of Wave-Induced Pore-Pressure Response to Development of Fluid Mud Layer

1
School of Civil Engineering, Qingdao University of Technology, Qingdao 266033, China
2
Key Laboratory of Shandong Province for Marine Environment and Geological Engineering, Ocean University of China, Qingdao 266100, China
3
Laboratory for Marine Geology, Qingdao National Laboratory for Maine Science and Technology, Qingdao 266061, China
4
Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2019, 7(9), 317; https://doi.org/10.3390/jmse7090317
Submission received: 17 July 2019 / Revised: 3 September 2019 / Accepted: 9 September 2019 / Published: 12 September 2019

Abstract

:
To investigate the vertical migration response of fine sediments, the pore pressure response of the silty seabed under the action of waves was tested. Under the action of waves, there is an obvious pumping phenomenon in the sludge accumulated by pore pressure. The excess pore water pressure caused by the waves in the seabed is unevenly distributed with respect to depth and there is an extreme value of up to 1.19 kPa. The pressure affects the liquefaction properties of the sludge. According to instantaneous-liquefaction judgment, the liquefaction of surface soil occurs, but the soil is not completely liquefied. Using theoretical calculations, the vertical source supply of floating mud development was analyzed. The pumping effect of the wave-induced excess pore pressure manifests in two aspects, as follows: (1) The centralized migration of splitting channels, which is visible to the naked eye, and (2) the general migration of fine particles between particle gaps at the mesoscopic level, which accounts for up to 22.2% of the migration of fine particles.

1. Introduction

Fluid mud refers to a layer of high-concentration sediment-laden water near the seabed. It is mainly composed of clay and silt with high viscosity and the predominant particle size is generally 63 μm [1]. Fluid mud has a clear interface with the upper water body and has high mobility. It has a unique sediment movement pattern in muddy coastal estuaries [2,3]. This has been observed in the Thames estuary in England, the Gironde estuary in France [4], the Amazon estuary in Brazil [5,6], and the Yellow River and Yangtze River estuaries in China [7].
The liquefaction of sand and the resuspension of sediment under the action of waves have long been the focus of coastal research. Most scholars focus on liquefaction mechanisms and resuspension fluxes [8,9,10,11]. However, few researchers pay attention to the vertical pumping effect of wave-induced pore pressure from a microscopic viewpoint, particularly the development of floating mud in silty soil [12]. A large amount of fluid mud can cause ecological, environmental, and navigation safety problems, such as burying benthic organisms [13], sediment eutrophication [14], and sudden channel siltation, disturbing sounding results and affecting the judgment of the navigable water depth [15]. Additionally, its interaction with the suspended load and bedload has an important impact on the topographic and geomorphological evolution of estuary areas [16,17], such as channel siltation and even navigation safety. It can block the channel and cause ships to run aground. In July 2015, the Yangtze River Estuary waterway was blocked by typhoon Canhong, which prevented the entrance of 15 huge ships and the port production was seriously affected.
The main research methods of fluid mud are field observation, laboratory test, theoretical analysis, and numerical simulation [18]. The problem of floating mud was discovered via field observations, which employ frozen sampling, bathymetry, the gamma-ray method, the ultrasonic method, the tuning fork density method, and the coupling method [19]. The measurement principles can be categorized as direct measurement, isotope measurement, and acoustic measurement. Laboratory tests are mainly divided into two types, as follows: The experimental study of the rheological properties of the sludge [20,21] and the experimental study of the movement characteristics of the sludge under dynamic action, including wave action [22,23,24] and current action [25,26,27,28]. On the basis of field observations and laboratory tests, the formation conditions [29] of sludge have been investigated and the numerical simulation [30] of the sludge response under dynamic action has been performed via theoretical analysis.
In view of the development mechanism of fluid mud, the prevailing view is that the conditions for the formation of fluid mud include an abundant fine sediment supply, relatively weak hydrodynamic conditions, and appropriate salinity. Salinity and a weak dynamic force are the external conditions for the flocculation and sedimentation of fine sediment, and abundant fine sediment provides the material source for the formation and development of fluid mud. Studies have been performed on the flocculation environment of fine sediment. Through field observations, Li found that the Yangtze River estuary has a good environment for flocculation of fine sediment [31]. Li and Zhang discovered that salinity and sediment concentration significantly contribute to the flocculation of fine sediment [32]. The periodic change of the flow velocity is the most significant factor affecting the flocculation size of fine sediment [33]. Through shallow section and dual-frequency sounding, Wen considered the stratification of the floating mud in the Yellow River estuary [34].
It is generally believed that there are two forms of material supply of floating mud. The first is direct input from rivers. Winterwerp and Li reported that the potential of density inflow is a special manifestation of the rapid deposition of sediment at the bottom, which provides a rich material source for the formation of floating mud [35,36]. Second, for sediment resuspension supply, whose dynamic mechanism is that waves, currents, storms, and other comprehensive effects on the seabed, the shear stress of the wave-current combined bed is significantly increased, leading to the resuspension of fine sediment near the coast [37,38]. Some researchers use numerical simulations to study the hydrodynamic effects that are responsible for the sediment resuspension. Ranasinghe, R simulates the coastal morphology of one-way flow by using the morphodynamic model Delft3D [39]. Gallerano proposed an integral formulation for the contravariant suspended sediment advection-diffusion equation and used it for the sea-bottom dynamic simulations [40]. Some non-hydrostatic models applied to shallow coastal waters have been proposed one after another [41,42,43]. Both of them are based on field observations and numerical simulations. Field observation reveals from a macro perspective. The analysis of numerical simulations is based on simple conditions, such as symmetrical bed disturbance and unidirectional flow. There is a lack of dynamic process research from a micro perspective, which is insufficient to fully explain the material supply mechanism of sludge development.
In this study, silt from the Yellow River Delta was the research object. The pore pressure accumulation process of soil under wave action was analyzed by wave flume test and the instantaneous liquefaction depth was calculated. Based on the analysis of particle size composition of soil, the vertical migration of soil particles under wave action was studied, and the contribution of fine particles to surface sludge was discussed. Through theoretical calculation of different D50 particles, the source and supply mechanism of sludge development are revealed from a mesoscopic point of view, in order to supplement and enrich the development mechanism and transport law of fluid mud.

2. Materials and Methods

2.1. Experiment Introduction

A wave-tank experiment was performed at the Geotechnical Laboratory of the Ocean University of China and the experimental setup is shown in Figure 1. The size of the sink was 3 m × 0.6 m × 1.2 m and the size of the lower tank was 1 m × 0.6 m × 0.5 m. An air compressor was used to power the unit. A double-acting cylinder was installed at one end of the water tank. The cylinder can provide pressure of 0.05–1.0 Mpa and the wave-making plate reciprocated to push the water body to generate waves. The waveform output was changed by adjusting the cylinder stroke and the cylinder damping while the other end was equipped with a wedge-shaped sponge block for eliminating the influence of wave reflection and ensuring the stable output of the waveform. As shown in the figure, the soil tank was filled with saturated silt with thickness of approximately 40 cm. A kaolin layer was set at a distance of 5 cm from the surface of the soil and its thickness was <0.5 cm. Before and after the wave action, the sample was allowed to stand for more than 4 h and the surface layer was extracted. The particle-size composition was determined using a laser particle-size analyzer. In this experiment, six YY-2B pore water pressure sensors were employed, which were located 5, 10, 15, 20, 25, and 30 cm from the surface of the soil. Before being buried, the pore-pressure sensors were soaked in clean water for 24 h and continuously vibrated to ensure that the internal air was completely discharged. The YY-2B pore water pressure sensor uses the imported piezoresistive pressure sensor as the sensitive device. The water pressure acts on the pressure film of the sensor through the permeable stone, which transforms the water pressure into the output of voltage signal and connects with the secondary instrument to measure the pore water pressure.
The water was artificial seawater with 35‰ salinity. The silt was taken from the tidal flat of the Yellow River estuary and had a median particle size (D50) ranging from 32 to 41 μm. The properties of the test sludge are shown in Table 1. Kaolin with 4000 meshes and D50 ranging from 2.1 to 2.8 μm was purchased from industrial factories.
Before the sensors were used to collect data, their sensitivity and rate were checked. The voltage values of the six pore-pressure sensors stabilized at different depths were recorded separately. According to the linear-regression method, the regression curves of the voltage values measured by the pore-pressure sensors and the true pore pressures were obtained, and the sensor calibration was completed. The results are presented in Figure 2.
The calibration results indicate that the pore-pressure sensor voltage display had a good linear relationship with the actual pressure value and could be used for experiments.

2.2. Quantification of Kaolin

The results of particle-size measurement can not only indicate parameters such as D50 and the average particle size (Dav), but also the distribution of the particle size and its cumulative distribution. By studying the grain-size characteristics of sediments, the transport mode of the sediments can be determined and the environmental factors affecting the grain-size change of the sediments can be identified, particularly the material sources and hydrodynamic environment [44,45].
The experimental silt and kaolin were obviously different with regard to their origins and particle-size compositions. In this experiment, the silt had a larger particle size (Figure 3; D50 of approximately 32–49 μm) and was taken from the tidal flat of the Yellow River estuary. It was formed by natural accumulation after long-distance transportation and long-term screening of the water flow. It had good sorting and grinding roundness and its material was relatively consistent, approximating the normal curve. The experimental kaolin had a finer particle size (Figure 3; D50 of approximately 2.1–2.8 µm). It was manufactured in industrial batches and screened artificially, exhibited poor sorting and grinding roundness, and its particles were sharper and finer than those of the silt.
Figure 4 shows the grain-size distribution curves of mixed kaolin and silt. The main component was 80% kaolin mixed with 20% silt. There are two independent peaks in the grain-size distribution curve, and the maximum values of the two peaks correspond to the grain-size ranges of kaolin (D50 = 2.1 µm) and silt (D50 = 34 µm). Therefore, the number of peaks in the particle-size distribution curve was used to qualitatively determine whether there were many components in the sample and whether the corresponding components could be qualitatively determined by the size range corresponding to the maximum value. The sample size distribution curve exhibited multiple peaks, indicating the superposition of single peaks with different particle sizes. The particle-size distribution of a component can be obtained by decomposing the particle-size distribution curve. The component proportions of different particle sizes (corresponding to different material sources) can be quantitatively assessed according to the cumulative distribution.

2.3. Experimental Procedure

To simulate the physical properties of the silt under natural consolidation conditions, the soil sample was saturated and an air-dried powder sample was selected for crushing to remove impurities such as grass and branches. The sieved soil sample was placed in a blender and standard seawater was added for agitation to ensure uniform mixing of the sample. The pore water pressure sensor was pre-fixed in the middle of the soil tank at a predetermined depth to prevent disturbance to the soil after embedding. The mixture was slowly injected along the inner wall into the soil tank until the thickness was approximately 35 cm. After the soil was stable, a layer of kaolin with a thickness of <0.5 cm was evenly deposited on the surface, and then soil was added until the thickness reached 40 cm. The experimental group without the kaolin layer directly injected the mixture into the 40 cm. After the soil was stable, standard seawater was slowly injected into the pool along the inner wall until it was approximately 35 cm from the soil surface. After 12 h, the soil naturally consolidated under the hydrostatic pressure and the wave loading test was started. Before and after the wave action, four surface samples were uniformly sampled in the soil tank using a sampling tube to measure the particle size. To prevent component interference between the different experimental groups, the samples in the artificial seawater and soil tank were replaced after the end of each test.
The wave action was applied by a wave generator installed at one end of the tank to experimentally examine different wave effects (wave height and duration). The experimental group settings are presented in Table 2. Experimental groups A and C were the basic control groups, focusing on the experimental pore-pressure response and surface particle-size changes caused by the wave action. The pumping effect of the wave excess pore water pressure on the fine sediment was investigated. The pore water pressure changes for soils with different depths before and after the wave action were recorded using a YY-2B pore water pressure sensor. On this basis, the calculation of the pore-pressure accumulation caused by waves and the migration of fine particles under pumping were calculated and analyzed.

3. Results

3.1. Experimental Phenomena

As the wave action began, the originally consolidated soil in the trough gradually appeared to fluctuate slightly. The fluctuation period of the soil was consistent with the wave-action period and the surface layer exhibited the largest fluctuation, with an amplitude of <2 cm. The start time of the soil fluctuation was delayed with the increasing depth. The wave eventually affected the depth of the soil and increased with an increase in the wave height. The maximum depth of influence visible to the naked eye reached 15 cm below the surface layer of the soil. The overlying water body on the soil gradually became mixed with clearness. The sediment on the surface of the soil was lifted and bubbles occasionally escaped from the soil. After the wave action stopped, the fluctuation of the upper soil quickly disappeared and the soil began to be consolidated. The particles suspended in the water slowly settled and the water body recovered after standing for a long time. At the boundary of the sink, obvious small channels were observed (Figure 5a), along with splitting channels 3–17 cm in length. Sediment deposits in the lower part of the soil followed these small channels, moving up quickly, and formed many tiny mounds on the surface of the soil. After standing for a period of time, there were numerous hilly deposits on the surface of the soil (Figure 5b). These deposits were generally <3 cm in diameter and ≤1 cm in height. The kaolin 5 cm below the surface of the soil quickly migrated along the channel to the surface layer and formed a white deposit. Groups C and D were the most obvious. Group D began with a large amount of white kaolin and then the color gradually became deeper. The soil rushed and gradually covered the white kaolin. After a long period of standing, the center of the mound-like accumulation gradually appeared to be concave. The whole process lasted approximately 1 h.

3.2. Pore Pressure Response

Figure 6 shows the pore-pressure response of the wave action (H = 7 cm, t = 30 min) 5–30 cm below the surface of the soil. The pore-pressure fluctuations of the sensors were essentially identical. Before the wave action, the pore water pressure of each layer was stable after a long period of static consolidation and the pore water pressure was equal to the overburden hydrostatic pressure. During the wave action (20–2100 s), the pore water pressure began to fluctuate with the wave load, and the pore water pressure in each layer began to accumulate (see Table 3; the maximum increment was 0.82–1.19 kPa, and the average increment was 0.76–0.94 kPa). The accumulative value and velocity of the pore water pressure in the upper layer were greater than those in the lower layer. The pore water pressure at 5 and 10 cm reached the maximum value in 150 s and then decreased slightly. The pore water pressure in each layer tended to be stable over time (they decreased by approximately 0.25–0.29 kPa). After the wave action stopped, the pore water pressure fluctuation disappeared rapidly and the pore water pressure in each layer dissipated continuously. After the wave action, the pore water pressure returned to the level observed before the wave action and remained stable for approximately 100 min.

3.3. Particle Size and Kaolin Content

Table 4 presents the changes in the average particle size and kaolin content of the soil surface after the wave action. The surface particle size of the soil was generally smaller after the wave action and the kaolin was generally detected in the surface layer, except for groups A and F. A wave 3 cm in height did not cause the kaolin to move up the surface layer. The average particle size of the sample decreased with an increase in the wave-action intensity (groups D, E, and F). With the increasing wave-action time, the kaolin content in the surface layer first increased and then decreased (groups B, C, and D). The decrease in the wave height exhibited a decreasing trend (groups C, E, and F).

4. Discussion

4.1. Pore-Pressure Response Analysis

A wave acting on a seabed causes an excess pore water pressure. The variation of the pore pressure in the soil and the amplitude are the factors that affect the liquefaction of the seabed soil. The data collected by the pore-pressure sensors in real time during the test represented the total pore water pressure, which comprised two parts, as follows: The pore water pressure and the static pore pressure, which does not work on the consolidation deformation and instability of the seabed soil. It is the excess pore water pressure, which is produced by the periodic loading of waves and can accumulate and dissipate in the soil.
The presence of excess pore pressure can reduce the effective stress of the seabed soil and affect the stability of the seabed. The collected data were processed, the hydrostatic pressure was subtracted from the total pore water pressure, and the excess pore pressure was analyzed.
Figure 7 presents the excess pore pressure data. Figure 7a indicates the variation of excess pore pressure with time at different depths. Nina Stark conducted a pore pressure observation at 5 cm and 20 cm depth in the intertidal zone of Cannon Beach in Yakutat, Alaska [46]. At a sediment depth of 20 cm, the period averaged pore pressure occasionally exceeded the initial mean effective stress and, at a sediment depth of 5 cm, the period averaged pore pressure exceeded the initial mean effective stress frequently. In our laboratory tests, at 5 cm and 20 cm, the periodic average pore pressure often exceeded the initial average effective stress. This may be related to the setting of the indoor test parameters, or it may be related to the characteristics of the liquefaction of the Yellow River estuary. The measurements at 20 cm sediment depth indicated that pore pressure build-up was particularly susceptible to maximum wave height. We also got similar results. This highlights the important role of wave height in pore pressure accumulation. As shown in Figure 7b, the pore water pressure was unevenly distributed with respect to the depth and there was an extreme value. The maximum value occurred 10 cm below the seabed; the pressure gradually decreased with the increasing depth. The reason for this phenomenon is that the super-porosity generated by the wave cyclic load in the seabed soil was a dynamic pore pressure, which accumulated and dissipated over time in the seabed soil. Although the wave acted on the surface of the seabed to produce an instantaneous maximum excessive pore pressure, the pore pressure in the soil near the surface of the seabed dissipated quickly, often without having to accumulate [47,48,49,50]. Therefore, the super-pore pressure was not maximized at the surface of the seabed; rather, it gradually accumulated with the increasing depth along the surface layer, reached a maximum at a certain depth from the surface layer, and (because the wave acted on the soil) gradually decreased with the increasing depth. Thus, at a certain depth, the pore pressure was gradually attenuated.

4.2. Discrimination of Instantaneous Liquefaction

The criterion of instantaneous liquefaction under wave action was used to compare the excess pore water pressure between the particles and the effective stress of the floating soil. Reaching the instantaneous-liquefaction standard only represents the trend of liquefaction at a certain time and does not indicate the liquefaction of the soil.
When the maximum excess pore water pressure Umax generated by the wave in the soil layer was equal to the overlying effective self-weight stress σ of the corresponding soil layer, the soil layer was in the limit equilibrium state. At Umax σ , the soil body was stable, and at Umax > σ , the soil layer was liquefied.
It was assumed that the soil skeleton was elastic, the pore water was incompressible, and the water was laminar in the pores. The transient response of the soil layer to the resonant wave load was analyzed. According to Yamomoto [51], a simple expression for the excess pore water pressure decay with respect to the depth was obtained, as follows:
U max = γ w H 2 cosh ( λ d ) e λ z
where z represents the soil thickness (m), λ represents the wavenumber (m–1), λ = 2 π L (here, L represents the wave length (m)), γ w represents the seawater gravity (kN/m3), H represents the wave height (m), and d represents the water depth (m).
The effective self-weight stress of the corresponding soil layer is expressed as σ = γ z , where γ represents the submerged specific weight of the soil (N/m3). By comparing the maximum excess pore water pressure, U max , with the overlying effective self-weight stress, σ , of the corresponding soil layer, the possibility of liquefaction at depth z can be judged.
By considering the effect of the lateral-pressure coefficient [52], the following expression is obtained:
γ w H 2 cosh ( λ d ) e λ z = 1 + 2 K 0 3 γ z
Through instantaneous-liquefaction discrimination calculations, the liquefaction depth was determined for each group, as shown in Table 5. The average liquefaction location of the soil was ≤5 cm. A smaller wave height yielded a lower liquefaction level. The cumulative response curve of the pore pressure under wave action indicated that although the cumulative pore pressure at the 5 cm horizon decreased, the cumulative excess pore pressure did not decrease to 0 instantaneously, suggesting that the liquefaction of this horizon occurred but that the soul was not completely liquefied.

4.3. Centralized Pumping Migration Based on Splitting Channels

During the experiment, a series of splitting channels were observed in the soil after the wave action. Driven by wave-induced excess pore pressure, particles in the lower layer moved rapidly along these channels to the surface and accumulated near the channel outlet. The upwelling upstream particles in group C were collected and analyzed. The results indicated that the upwelling particles were smaller than the other particles, having an average size of only 9.5 µm. They were mainly composed of kaolin and fine silt, and the kaolin accounted for approximately 56.4% (Figure 8).
Since the permeability of the silt was relatively poor (according to experimental measurements, the permeability of the sample used was approximately 1 × 10−8 m/s), after the wave action was stopped, the accumulated excess pore pressure of the silt could not be quickly dissipated. Additionally, owing to the significant decrease in the strength of the soil, the accumulated excess pore pressure was rapidly dissipated along the weak points (or weak paths) of the soil layer through the destruction of the soil layer, forming a series of splitting channels. As the excess pore pressure dissipated, the liquefied soil particles in the fluid state in the soil separated from the soil skeleton. Under the action of the penetrating force, the particles moved rapidly along the channels to the surface. The fine particles were more likely to move.
Generally, a longer wave-action time (groups B, C, and D) corresponded to a larger decrease in the soil strength, the formation of a larger number of upwelling channels, and a higher kaolin content detectable on the surface layer of the soil. In the experiment, the D-channel at the beginning of the group D was obviously white. Then, the color gradually deepened and transitioned to the color of the silt, and the continuous upwelling did not end. This indicates that the initial particle granulation was dominated by kaolin with a small particle size. As the upwelling continued, the amount of kaolin near the channel gradually decreased (the kaolin thickness was <0.5 cm). The upwelling granules became mainly silt, and a large amount of silt covering the kaolin made the kaolin content relatively low. Generally speaking, in a certain period of time, the smaller the wavelength of wave action, the smaller the degree of wave action on the bottom kaolin, the fewer upwelling channels formed, the smaller the amount of fine particles upwelling, and the less kaolin content can be monitored on the surface of the soil.

4.4. Particle-Force Analysis

For the upper soil (0–15 cm) with full wave action, without considering the interaction between particles, the vertical forces of particles under the wave action mainly included the effective gravity of the particles, the uplift force, the seepage force, the cohesion force, and additional pressure of the membrane water (Figure 9).
(1) Effective gravity F g , i.e., sediment gravity subtracted by buoyancy;
F g = ( ρ s ρ w ) g π 6 d 50 3
(2) Uplift force F y , which can be expressed as follows according to Cao [53]:
F y = [ 1 4.684 ( L / h ) 0.322 ] 2 1 2 ρ ω π 4 d 50 2 π 2 H 2 T 2 sinh 2 ( k h ) cos 2 ( k x ω t )
(3) Seepage force F s , which can be expressed as follows according to Mei [54] and Xia [55]:
F s = n 1 n m m + 1 π d 50 3 6 ρ ω g H 2 cosh ( k h ) 1 δ sin ( k x ω t + π 4 ) .
(4) Bonding force and additional pressure of membrane water, which can be expressed as follows according to Mei and Cao:
N 1 = φ π 2 ε d 50 . N 2 = φ ρ w h π 2 d 50 α .
Here, ρ s = 2650 kg/m3 represents the density of soil particles, ρ w = 1030 kg/m3 represents the density of water, g = 9.8 N/kg represents the gravitational acceleration, H = 0.33 m represents the water depth, T = 1 s represents the wave period, H = 0.07 m represents the wave height, L = 1 m represents the wave length, k = 1/L represents the wavenumber, ω represents the angular frequency, δ represents the thickness of the seabed boundary layer, δ = k d G ω ( n G β + 1 2 ν 2 ( 1 ν ) ) , where n = 0.4 represents the porosity, β = 1.0 × 107 N/m2 represents the bulk elastic modulus, G = 5.0 × 106 N/m2 represents the shear modulus, v = 0.33 represents the Poisson’s ratio, and k d = k s ρ ω g (with k s = 1.0 × 10−8 m/s) is the permeability coefficient, m = n G ( 1 2 ν ) β presents the relative compressibility of water and the soil skeleton, φ = 1 16 is the correction coefficient, ε = ε k ρ w is the coefficient of cohesion, ε k = 2.56 × 10−2 m/s2, and α = 0.213 × 10−6 m is the characteristic thickness related to the size of the sand gap.
Using the foregoing formulas, the forces acting on the D50 = 3 μm, 30 μm, and 300 μm particles were calculated, as shown in Figure 9. For soil particles with D50 = 3 μm, the vertical movement is mainly affected by the uplift force, cohesive force, and the additional pressure of membrane water. For soil particles with D50 = 30 μm, the vertical movement is mainly affected by the uplift force and the seepage force. For soil particles with D50 = 300 μm, the vertical movement is mainly affected by the effective gravity, uplift force, and seepage force. For particles with a relatively small size (D50 = 3 or 30 μm), the resultant force (positive orientation) is positive with time integration (impulse), indicating that the momentum of the particles increases with time. The impulse of particles with a relatively large size (D50 = 300 μm) is negative (Figure 9).
In summary, the pumping effect of the wave-induced excess pore pressure is mainly manifested in two aspects, as follows (Figure 10): (1) The visible centralized migration of splitting channels and (2) the general migration of fine particles between granular cracks on the mesoscale.

5. Conclusions

The silt in the Yellow River Delta was taken as the research object and the main components of the floating mud were simulated with 4000 mesh kaolin. The responses of the pore-pressure accumulation, soil liquefaction, and particle movement under wave action were experimentally investigated. According to the experimental results and calculations, the following conclusions were drawn.
(1) Resuspension of surface soil particles occurred owing to wave action and the “pumping” effect of the wave-induced pore pressure transferred the fine kaolin from the lower layer to the surface layer, resulting in a reduced surface soil particle size after the wave action.
(2) Under the action of waves, the excess pore pressure gradually accumulated with the increasing depth along the surface layer, reaching its maximum value (1.19 kPa) 10 cm from the surface layer. The excess pore pressure gradually decreased as the depth increased further.
(3) The silt particle size was small. The cohesive force between the particles and the additional pressure of the film water had significant effects on the liquefaction characteristics of the silt. The strength of the surface soil decreased significantly. According to instantaneous-liquefaction judgment, the liquefaction of the surface soil occurred, but the surface soil was not completely liquefied.
(4) The pumping effect of the wave-induced excess pore pressure was manifested in two aspects, as follows: The centralized migration of splitting channels visible to the naked eye and the general migration of fine particles between particle gaps at the mesoscopic level, which contributed up to 22.2% of the migration of fine particles.
The accumulation of pore pressure and the upward movement of fine particles are complex. The mechanism of formation and migration of floating mud is still unclear. More studies are needed, especially with field data in order to study the detailed interaction of these processes and their correlation with beach evolution. Despite the small size of the sensors (~6 cm in length and ~2.5 cm in diameter), the sediment was disturbed by the introduction of the sensors, and certainly for the beginning of the measurements. This research for this paper was carried out under ideal regular wave conditions, which is quite different from the complex marine dynamic environment. We will further carry out in-situ testing, especially with field data to study the detailed interaction of these processes, and calibrate and verify them with laboratory tests.

Author Contributions

Conceptualization, Z.Y.; formal analysis, Y.Z. and Z.S.; methodology, T.L.; software, Y.Z. and Y.Z.; validation, X.L. and Z.S.; writing-original draft, Z.Y.; writing-review & editing, T.L. and Z.Y.

Funding

This study was supported by the National Natural Science Foundation of China (Project Nos. 41602318 and 41672272), the Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology, grant No. MGQNLM-KF201710, the Fundamental Research Funds for the Central Universities, grant No. 201962011, the National Key Research and Development Program (2017YFC0307701), and the foundation for basic research of the Ocean University of China (Grant No. 201861041). We thank the Geotechnical Laboratory of the Ocean University of China for providing equipment for the experiment.

Acknowledgments

We appreciate the anonymous reviewers who provided comments for revising the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, W.H.; Shi, L.Q.; Liu, M. Review of estuarine and coastal fluid mud measurement technique, characteristics and transportation researches. J. Sediment Res. 2013, 1, 74–80. [Google Scholar]
  2. Mcanally, W.H.; Friedrichs, C.T.; Hamilton, D. Management of Fluid Mud in Estuaries, Bays, and Lakes. I: Present State of Understanding on Character and Behavior. J. Hydraul. Eng. 2007, 133, 9–22. [Google Scholar] [CrossRef]
  3. Qian, N.; Wan, Z.H. Mechanics of Sediment Transport; Science Press: Beijing, China, 1983; pp. 158–159. [Google Scholar]
  4. Granboulan, J.; Feral, A.; Villerot, M. Study of the sedimentological and rheological properties of fluid mud in the fluvio-estuarine system of the Gironde estuary. Ocean Shorel. Manag. 1989, 12, 23–46. [Google Scholar] [CrossRef]
  5. Allison, M.A.; Nittrouer, C.A.; Kineke, G.C. Seasonal sediment storage on mudflats adjacent to the Amazon River. Mar. Geol. 1995, 125, 303–328. [Google Scholar] [CrossRef]
  6. Kineke, G.C.; Sternberg, R.W. Distribution of fluid muds on the Amazon continental shelf. Mar. Geol. 1995, 125, 193–233. [Google Scholar] [CrossRef]
  7. Xu, H.G.; Xu, H.T.; Li, J.F. Study on the nautical depth utilizing fluid mud layer of the Changjiang estuary. J. East China Norm. Univ. (Nat. Sci.) 1994, 2, 91–97. [Google Scholar]
  8. Wainright, S. Sediment-to-water fluxes of particulate material and microbes by resuspension and their contribution to the planktonic food web. Mar. Ecol. Prog. Ser. 1990, 62, 271–281. [Google Scholar] [CrossRef]
  9. Gao, Y.F.; Yang, S.; Jian, Z. Analysis of wave-induced liquefaction in seabed deposits of silt. China Ocean Eng. 2011, 25, 31–42. [Google Scholar] [CrossRef]
  10. Xu, H.; Dong, P. A probabilistic analysis of random wave-induced liquefaction. Ocean Eng. 2011, 38, 860–867. [Google Scholar] [CrossRef]
  11. Moriarty, J.M.; Harris, C.K.; Fennel, K. The roles of resuspension, diffusion and biogeochemical processes on oxygen dynamics offshore of the Rhône River, France: A numerical modeling study. Biogeosciences. 2017, 14, 1919. [Google Scholar] [CrossRef]
  12. Zhang, S.; Jia, Y.; Wen, M. Vertical Migration of Fine-Grained Sediments from Interior to Surface of Seabed Driven by Seepage Flows—“Sub-Bottom Sediment Pump Action”. J. Ocean Univ. China 2017, 16, 15–24. [Google Scholar] [CrossRef]
  13. Corselli, C.; Basso, D. First evidence of benthic communities based on chemosynthesis on the Napoli mud volcano (Eastern Mediterranean). Mar. Geol. 1996, 132, 227–239. [Google Scholar] [CrossRef]
  14. Montserrat, F.; Suykerbuyk, W.; Al-Busaidi, R. Effects of mud sedimentation on lugworm ecosystem engineering. J. Sea Res. 2011, 65, 170–181. [Google Scholar] [CrossRef]
  15. Niu, G.Z.; Shen, X.M.; Pei, W.B. The test of SILAS sounding system. Hydrogr. Surv. Chart. 2003, 23, 28–31. [Google Scholar]
  16. Yang, Y.; Li, Y.; Sun, Z.; Fan, Y. Suspended sediment load in the turbidity maximum zone at the Yangtze River Estuary; The trends and causes. J. Geogr. Sci. 2014, 24, 129–142. [Google Scholar] [CrossRef]
  17. Liu, T.; Li, S.P.; Kou, H.L. Excess pore pressure observation in marine sediment based on Fiber Bragg Grating pressure sensor. Mar. Georesources Geotechnol. 2019, 37, 775–782. [Google Scholar] [CrossRef]
  18. Shi, Z. Physical Processes of Fine Sediment in Estuaries and Coasts; Shanghai Jiaotong University Press: Shanghai, China, 2013; pp. 102–105. [Google Scholar]
  19. Ju, Y.; Gao, M.; Wang, Y.Y. Review on field observation technologies for fluid mud. J. Chongqing Jiaotong Univ. (Nat. Sci.) 2014, 33, 98–102. [Google Scholar]
  20. Maa, J.P.Y.; Mehta, A.J. Soft Mud Properties: Voigt Model. J. Waterw. Port Coast. Ocean Eng. 1988, 114, 765–770. [Google Scholar] [CrossRef]
  21. Abujdayil, B.; Banat, F.; Alsameraiy, M. Steady Rheological Properties of Rotating Biological Contactor (RBC) Sludge. J. Water Resour. Prot. 2010, 2, 1–7. [Google Scholar] [CrossRef] [Green Version]
  22. Yeh, H.; Masn, H.B. Sediment response to tsunami loading: mechanisms and estimates. Géotechnique 2014, 64, 131–143. [Google Scholar] [CrossRef]
  23. Wang, Z.T.; Luan, M.T.; Dong-Sheng, J. Theoretical analysis of random wave-induced seabed response and liquefaction. Rock Soil Mech. 2008, 29, 2051–2056. [Google Scholar]
  24. Lambrechts, J.; Humphrey, C.; Mckinna, L. Importance of wave-induced bed liquefaction in the fine sediment budget of Cleveland Bay, Great Barrier Reef. Estuar. Coast. Shelf Sci. 2010, 89, 154–162. [Google Scholar] [CrossRef]
  25. Hong, R.J.; Ying, Y.L. Experimental Study on Starting Flow Velocity of Floating Sludge under Water Flow. J. Hydraul. Eng. 1988, 8, 51–57. [Google Scholar]
  26. Kreeke, J.V.D.; Day, C.M.; Mulder, H.P.J. Tidal variations in suspended sediment concentration in the Ems estuary: origin and resulting sediment flux. J. Sea Res. 1997, 38, 1–16. [Google Scholar] [CrossRef]
  27. Almroth, E.; Tengberg, A.; Andersson, J.H. Effects of resuspension on benthic fluxes of oxygen, nutrients, dissolved inorganic carbon, iron and manganese in the Gulf of Finland, Baltic Sea. Cont. Shelf Res. 2009, 29, 807–818. [Google Scholar] [CrossRef]
  28. Wu, W.; Soligo, G.; Marchioli, C.; Soldati, A.; Piomelli, U. Particle resuspension by a periodically forced impinging jet. J. Fluid Mech. 2017, 820, 284–311. [Google Scholar] [CrossRef]
  29. Xu, J.Y.; Yuan, J.Z. Study on the Fluid Mud in the Yangtze Estuary. J. Sediment Res. 2001, 3, 74–80. [Google Scholar]
  30. Maa, J.P.Y.; Mehta, A.J. Soft mud response to water waves. J. Waterw. Port Coast. Ocean Eng. 1990, 116, 634–650. [Google Scholar] [CrossRef]
  31. Li, D. Flocculation process of fine-grained sediments by the combined effect of salinity and humus in the Changjiang Estuary. Acta Oceanol. Sin. 2007, 26, 140–149. [Google Scholar]
  32. Li, J.; Zhang, C. Sediment resuspension and implications for turbidity maximum in the Changjiang Estuary. Mar. Geol. 1998, 148, 117–124. [Google Scholar] [CrossRef]
  33. Li, J.F.; Dai, Z.J.; Liu, Q.Z. In-situ observation of floccule size and fluid mud in the Changjiang Estuary. J. Sediment Res. 2008, 3, 26–32. [Google Scholar]
  34. Mingzheng, W.; Zhenhao, W.; Bowen, Z. Survey on the distribution of fluid mud and disturbed strata on subaqueous Yellow River Delta. J. Eng. Geol. 2018, 26, 677–683. [Google Scholar] [CrossRef]
  35. Winterwerp, J.C. On the Dynamics of High-Concentrated Mud Suspensions; Netherlands, Delft University of Technology: Delft, The Netherlands, 1999; pp. 58–59. [Google Scholar]
  36. Li, D.S.; Ren, R.S. Research for special property of gravity current with the movement of wave. Ocean Eng. 1996, 4, 54–59. [Google Scholar]
  37. Miao, L.M.; Yang, S.L.; Zhu, Q. Variations of suspended sediment concentrations and transport in response to a storm and its dynamic mechanism. J. Oceanogr. 2016, 38, 158–167. [Google Scholar]
  38. Sheremet, A.; Stone, G.W.; Zhang, X. Wave-induced sediment resuspension on a muddy inner shelf during Hurricane Claudette. Estuar. Coast. Shelf Sci. 2005, 63, 225–233. [Google Scholar] [CrossRef]
  39. Ranasinghe, R.; Swinkels, C.; Luijendijk, A.; Roelvink, D. Morphodynamic upscaling with the MORFAC approach: Dependencies and sensitivities. Coast. Eng. 2011, 58, 806–811. [Google Scholar] [CrossRef]
  40. Gallerano, F.; Cannata, G.; Scarpone, S. Bottom changes in coastal areas with complex shorelines. Eng. Appl. Comput. Fluid Mech. 2017, 11, 396–416. [Google Scholar] [CrossRef]
  41. Bradford, S.F. Non-hydrostatic model for surf zone simulation. J. Waterw. Port Coast. Ocean Eng. 2011, 137, 163–174. [Google Scholar] [CrossRef]
  42. Ma, G.; Shi, F.; Kirby, J.T. Shock-capturing non-hydrostatic model for fully dispersive surface wave processes. Ocean Model. 2012, 4344, 22–35. [Google Scholar] [CrossRef]
  43. Gallerano, F.; Cannata, G.; Lasaponara, F.; Petrelli, C. A new three-dimensional finite-volume non-hydrostatic shock-capturing model for free surface flow. J. Hydrodyn. 2017, 29, 552–566. [Google Scholar] [CrossRef]
  44. Mclaren, P.; Bowles, D. The effects of sediment transport on grain-size distributions. J. Sediment Res. 1985, 55, 457–470. [Google Scholar]
  45. Wang, W.; Li, A.C.; Xu, F.J. Distribution of surface sediments and sedimentary environment in the north yellow sea. Oceanol. Limnol. Sin. 2009, 40, 525–531. [Google Scholar]
  46. Stark, N. Pore pressure response to irregular waves at a sandy beach. In Proceedings of the Geotechnical Frontiers 2017, Orlando, FL, USA, 12–15 March 2017; pp. 409–417. [Google Scholar] [CrossRef]
  47. Guo, Z.; Jeng, D.S.; Zhao, H.; Guo, W.; Wang, L. Effect of seepage flow on sediment incipient motion around a free spanning pipeline. Coast. Eng. 2019, 143, 50–62. [Google Scholar] [CrossRef]
  48. Li, K.; Guo, Z.; Wang, L.; Jiang, H. Effect of seepage flow on shields number around a fixed and sagging pipeline. Ocean Eng. 2019, 172, 487–500. [Google Scholar] [CrossRef]
  49. Tao, L.; Guanli, W.; Hailei, K.; Lei, G. Pore pressure observation: pressure response of probe penetration and tides. Acta Oceanol. Sin. 2019, 38, 107–113. [Google Scholar] [CrossRef]
  50. Liu, T.; Lu, Y.; Zhou, L.; Yang, X.; Guo, L. Experiment and Analysis of Submarine Landslide Model Caused by Elevated Pore Pressure. J. Mar. Sci. Eng. 2019, 7, 146. [Google Scholar] [CrossRef]
  51. Yamamoto, T. Wave-induced pore pressures and effective stresses in inhomogeneous seabed foundations. Ocean Eng. 1981, 8, 1–16. [Google Scholar] [CrossRef]
  52. Jia, Y.; Zhang, L.; Zheng, J. Effects of wave-induced seabed liquefaction on sediment re-suspension in the Yellow River Delta. Ocean Eng. 2014, 89, 146–156. [Google Scholar] [CrossRef]
  53. Cao, W.H.; Liu, Q.Q. Analysis on dynamic mechanism of sediment winnowing caused by wave. J. Hydraul. Eng. 2000, 1, 49–53. [Google Scholar]
  54. Mei, C.C.; Foda, M.A. Wave-induced responses in a fluid-filled poro-elastic solid with a free surface—A boundary layer theory. Geophys. J. R. Astron. Soc. 2010, 66, 597–631. [Google Scholar] [CrossRef]
  55. Xia, L. Sediment Threshold and Local Scour Around Submarine Pipelines under Wave Action. Ph.D. Thesis, Zhejiang University, Hangzhou, China, 2006. [Google Scholar]
Figure 1. (a) Experimental device; (b) pore-pressure sensor placement; (c) pore-pressure sensor.
Figure 1. (a) Experimental device; (b) pore-pressure sensor placement; (c) pore-pressure sensor.
Jmse 07 00317 g001
Figure 2. Calibration results for the pore-pressure sensors: (a) 5 cm from the soil surface; (b) 10 cm from the soil surface; (c) 15 cm from the soil surface; (d) 20 cm from the soil surface; (e) 25 cm from the soil surface; and (f) 30 cm from the soil surface.
Figure 2. Calibration results for the pore-pressure sensors: (a) 5 cm from the soil surface; (b) 10 cm from the soil surface; (c) 15 cm from the soil surface; (d) 20 cm from the soil surface; (e) 25 cm from the soil surface; and (f) 30 cm from the soil surface.
Jmse 07 00317 g002
Figure 3. Particle-size distribution curves of the silt and kaolin used in the experiment.
Figure 3. Particle-size distribution curves of the silt and kaolin used in the experiment.
Jmse 07 00317 g003
Figure 4. Particle-size distribution curves of mixed kaolin and silt.
Figure 4. Particle-size distribution curves of mixed kaolin and silt.
Jmse 07 00317 g004
Figure 5. Upward migration and accumulation of fine particles: (a) Channels observed at the boundary of the flume and numerous hilly deposits; (b) numerous hilly deposits formed on the surface of the soil.
Figure 5. Upward migration and accumulation of fine particles: (a) Channels observed at the boundary of the flume and numerous hilly deposits; (b) numerous hilly deposits formed on the surface of the soil.
Jmse 07 00317 g005
Figure 6. Pore-pressure responses at different depths before, during, and after wave action (H = 7 cm, t = 30 min).
Figure 6. Pore-pressure responses at different depths before, during, and after wave action (H = 7 cm, t = 30 min).
Jmse 07 00317 g006
Figure 7. Excess pore pressure data. (a) Variation of excess pore pressure with time at different depths. (b) Schematic diagram of maximum cumulative pore pressure at different depths.
Figure 7. Excess pore pressure data. (a) Variation of excess pore pressure with time at different depths. (b) Schematic diagram of maximum cumulative pore pressure at different depths.
Jmse 07 00317 g007
Figure 8. Size-distribution curve of the upwelling soil particles in channels.
Figure 8. Size-distribution curve of the upwelling soil particles in channels.
Jmse 07 00317 g008
Figure 9. Force changes of particles under wave action in the vertical direction. (a) D50 = 3 µm; (b) D50 = 30 µm; and (c) D50 = 300 µm.
Figure 9. Force changes of particles under wave action in the vertical direction. (a) D50 = 3 µm; (b) D50 = 30 µm; and (c) D50 = 300 µm.
Jmse 07 00317 g009
Figure 10. Two types of pumping effects. (a) Fine particles migrate in the mud diapir; (b) fine particles migrate between gaps.
Figure 10. Two types of pumping effects. (a) Fine particles migrate in the mud diapir; (b) fine particles migrate between gaps.
Jmse 07 00317 g010
Table 1. Properties of the test silt.
Table 1. Properties of the test silt.
PropertyValue
Moisture Content (ω) 28.1%
Density (ρ)2.0 g/cm3
Dry density (ρd)1.55 g/cm3
Specific gravity (Gs)2.7
Void ratio (e)0.74
Porosity (n)44%
Plasticity index (Ip)6.9
Saturated unit weight (γs)16.0 kN/m3
Shear strength of uu testCohesion (c)7 kPa
Internal friction angle (φ)20°
Table 2. Experiment setup.
Table 2. Experiment setup.
GroupKaolin LayerSample NumberWave ActionSample Number
Time (min)Wave Height (cm)Wave Length (cm)Wave Period (s)
A×1–420710015–8
B9–12107100113–16
C17–20207100121–24
D25–28307100129–32
E33–36205900.937–40
F41–44203800.945–48
Table 3. Analysis of the pore-pressure variation during wave action.
Table 3. Analysis of the pore-pressure variation during wave action.
PositionBeforeWave ActionAfter
Initial Pore Water Pressure u 0 ¯   ( kPa ) Maximum Pore Water Pressure umax (kPa)Maximum Cumulative Pore Pressure ∆umax (kPa)Stable Pore Pressure u ¯ (kPa)Stable Cumulative Pore Pressure Δ u ¯   ( kPa ) Ultimate Pore Pressure u 1 ¯   ( kPa )
5 cm3.354.401.054.110.763.24
10 cm3.825.021.194.750.943.82
15 cm4.755.620.875.550.804.75
20 cm5.115.960.855.890.785.14
25 cm5.856.680.836.620.775.76
30 cm6.116.930.826.890.786.15
Note: ∆umax = umax u 0 ¯ , Δ u ¯ = u ¯ u 0 ¯ .
Table 4. Changes in the particle size and kaolin content after wave action.
Table 4. Changes in the particle size and kaolin content after wave action.
GroupAverage Size of the Surface Particles before Wave Action Dav0 (μm)Average Size of the Surface Particles after Wave Action Dav1 (μm)Kaolin Content after Wave Action (%)
A 44.9 36.5 0.0
B33.627.63.9
C39.016.122.2
D44.018.611.2
E44.419.84.5
F41.234.00.0
Table 5. Liquefaction depth for each group.
Table 5. Liquefaction depth for each group.
GroupDepth of Water (m)Wave Height (m)Wave Length (m)Liquefaction Critical Thickness (m)
A0.330.0710.04
B0.330.0710.04
C0.330.0710.04
D0.330.0710.04
E0.330.050.90.03
F0.330.030.80.02

Share and Cite

MDPI and ACS Style

Yang, Z.; Zhu, Y.; Liu, T.; Sun, Z.; Ling, X.; Zheng, Y. Contribution of Pumping Action of Wave-Induced Pore-Pressure Response to Development of Fluid Mud Layer. J. Mar. Sci. Eng. 2019, 7, 317. https://doi.org/10.3390/jmse7090317

AMA Style

Yang Z, Zhu Y, Liu T, Sun Z, Ling X, Zheng Y. Contribution of Pumping Action of Wave-Induced Pore-Pressure Response to Development of Fluid Mud Layer. Journal of Marine Science and Engineering. 2019; 7(9):317. https://doi.org/10.3390/jmse7090317

Chicago/Turabian Style

Yang, Zhongnian, Yongmao Zhu, Tao Liu, Zhongqiang Sun, Xianzhang Ling, and Yuqian Zheng. 2019. "Contribution of Pumping Action of Wave-Induced Pore-Pressure Response to Development of Fluid Mud Layer" Journal of Marine Science and Engineering 7, no. 9: 317. https://doi.org/10.3390/jmse7090317

APA Style

Yang, Z., Zhu, Y., Liu, T., Sun, Z., Ling, X., & Zheng, Y. (2019). Contribution of Pumping Action of Wave-Induced Pore-Pressure Response to Development of Fluid Mud Layer. Journal of Marine Science and Engineering, 7(9), 317. https://doi.org/10.3390/jmse7090317

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop