3.2. Wet Bulk Density, Particle Size Distribution and TOC Content
The change in excess wet bulk density (
, with
= density of water), obtained by the difference between excess wet bulk density of degraded and fresh mud samples, is plotted as a function of excess wet bulk density of fresh mud samples (
Figure S1). It is seen that the change in wet bulk density incurred during the long-term (250 days) incubation to degrade organic matter was generally lower than 5% of its original value. It was, therefore, assumed that changes in rheological properties between fresh and degraded samples were not related to the change in density or water content.
Mud samples obtained from PoB, PoE, PoH and PoR showed quite similar particle size distribution (
Figure S2a) along with the similar D
50 values (see
Table S2). Conversely, the mud samples collected from PoA displayed two different behaviors (
Figure S2b), i.e., one was similar to the other ports while another one showed bimodal particle size distribution along with the higher D
50 values (see samples A2 and A4 in
Table S2), which was linked with the presence of a significant amount of sand sized particles. For these two samples (A2 and A4), the higher sand content was also found in the particle size distribution obtained by the sieving technique [
20], which verified the SLS results (see
Table S3). Some other sediment properties such as TOC content, electrical conductivity and pH of mud samples from different ports are presented in
Table S1.
3.3. Stress Ramp-Up Tests
Stress ramp-up tests were performed to analyze the effect of organic matter degradation on the yield stresses of mud from different ports.
Figure S3 shows the outcome of stress ramp-up tests in terms of apparent viscosity as a function of shear stress for fresh mud samples and degraded mud samples from various ports. A two-step yielding behavior was evident from the two declines in viscosity. A similar two-step yielding behavior has already been reported in the literature for fine-grained sediments [
3,
21]. This two-step yielding behavior in mud samples was associated with the structural reorganization during shearing, i.e., breakage of interconnected networks of flocs, formation of rolls/cylinder-like structures and collapse of cylinder-like structures [
22]. The yield stress values, corresponding to these viscosity declines, were defined as static yield stress (
) and fluidic yield stress (
).
Figure S3 shows that the degraded mud samples had lower yield stress values as compared to the fresh mud samples, except for mud samples from PoB and PoE. This decrease in yield stress values may be associated with the degradation of organic matter and the corresponding decline of the extent of mineral bridging, which eventually results in a weaker system [
6]. A similar decrease in yield stress values by chemically removing the organic matter from mud was also reported in the literature [
23]. The mud samples from PoB and PoE showed either similar yield stress values or even higher yield stress values after organic matter degradation.
In order to further quantify the influence of organic matter degradation on the yield stress of mud samples from different ports, the change in static and fluidic yield stresses (degraded–fresh) was plotted as a function of yield stress of fresh samples (see
Figure 2a,c). A strong inverse correlation was observed between the change in yield stresses and the yield stresses of fresh mud samples, i.e., the higher the original yield stress, the higher its reduction after organic matter degradation (i.e., negative slope) [
13,
14]. Moreover, the following empirical equation was used to fit the experimental data of change in the rheological property (degraded–fresh) as a function of the same rheological property before degradation:
where
(Pa or Pa s
−1) and
(–) represent the intercept and slope of the line, respectively. The values of these fitting parameters are presented in
Table S4. The values of the slope indicated that the yield stresses (static and fluidic) were significantly reduced for mud samples from PoA and PoH due to the degradation of organic matter. In the case of mud samples from PoB and PoR, the reduction in yield stresses after organic matter degradation was less significant, as evident from the smaller values of negative slope (
Table S4). The mud samples collected from PoE exhibited smaller densities as compared to the samples collected from other ports (see
Table S1). Therefore, the change in yield stresses (degraded–fresh) as a function of yield stresses of fresh mud samples for PoE is compared with the fluid mud samples collected from PoH with similar densities (
Figure 2b,d, static and fluidic yield stresses). It is observed that the mud samples from PoE displayed an increase in both static and fluidic yield stress values after organic matter degradation (i.e., positive slope), while in the case of PoH, a pronounced negative slope was found with significant reduction in both static and fluidic yield stresses (
Table S4).
The aforementioned results were plotted in terms of percent change in yield stresses
after organic matter degradation for different ports (
Figure 3). These results identified a significant reduction in yield stresses for PoA and PoH (~20–40% mean value), a less pronounced decrease in yield stresses for PoR (~10–20% mean value), a negligible change in yield stresses for PoB (~0–5% mean value), and a significant increase in yield stresses for PoE (~40–70% mean value). The significant decrease in yield stresses for mud from PoA and PoH may be attributed to the active nature of the organic matter, particularly for PoH [
14], which degrades under suitable conditions. The small or negligible changes in yield stresses after degradation for PoB may be associated with the lower organic matter decay as compared to the other ports. Conversely, the increase in yield stresses after degradation of organic matter in mud from PoE is quite unexpected. It is hypothesized that under laboratory conditions, i.e., in a closed incubation system with no addition of fresh new carbon and favorable conditions of temperature, the composition of the microbial community, its abundance, and thereby the floc/biofilm architecture also change. This means that as organic matter (OM) is decayed, and is released as CH
4 and CO
2, the remaining OM undergoes a transformation toward more stable organo-mineral associations with concurrent changes of the microbial community composition, compared to the fresh sample. As a result, different floc and biofilm properties change, which could result in a different strength. It is assumed that this effect also occurs in the samples from the other ports but might be masked due to the high OM decay and the significant mass loss of carbon, which dominate the rheological response by disrupting mineral bridging. Hence, due to low mass loss and OM decay rates in PoE (see
Figure 1), other processes changing the physical architecture of flocs and biofilms might become more relevant.
3.4. Amplitude Sweep Tests
In addition to the yield stress, the transition between solid-like and liquid-like behavior of mud samples was investigated by performing oscillatory amplitude sweep tests.
Figure S4 shows the phase angle as a function of oscillatory amplitude for fresh and degraded mud samples. The solid–liquid transition point was estimated when the phase angle became equal to 45° (i.e., crossover point between
and
), and the corresponding stress amplitude was demonstrated as crossover amplitude.
Figure S4 illustrates a decrease in crossover amplitude by the degradation of organic matter, except for samples from PoE and PoB, which is in agreement with the yield stress analysis.
Change in crossover amplitude (degraded–fresh) is now plotted as a function of crossover amplitude of fresh mud samples from different ports (see
Figure 4a). The organic matter degradation again showed a significant influence on the crossover amplitude for the mud samples having higher crossover amplitude before degradation. The experimental data were then fitted with Equation (1), and the values of
and
are given in
Table S4. The slope values are quite similar in the case of crossover amplitude and yield stresses for mud samples from different ports. This behavior shows that the influence of organic matter degradation on the extent of decrease in rheological properties (crossover amplitude and yield stresses) is quite similar, because both these properties involve the application of a shearing action to disturb the intact structure. The change in crossover amplitude (degraded–fresh) as a function of crossover amplitude of fresh mud samples for PoE was again compared with the fluid mud samples collected from PoH with similar densities (
Figure 4b). It is again evident that the mud samples from PoE displayed an increase in crossover amplitude values after organic matter degradation (i.e., positive slope), while in the case of PoH, a pronounced negative slope was observed with significant reduction in crossover amplitude (
Table S4).
The percent change in crossover amplitude
after organic matter degradation for mud samples from different ports is shown in
Figure 5. These results again identified a significant reduction in crossover amplitude for PoA, PoH and PoR (~20–30% mean value), a negligible decrease in crossover amplitude for PoB (~0–5% mean value), and a significant increase in crossover amplitude for PoE (~30% mean value). These results verified the outcome of yield stress analysis and can be attributed to the processes changing the physical architecture of flocs and biofilms.
3.5. Frequency Sweep Tests
Frequency sweep tests were performed within the LVE regime from 0.1 to 100 Hz in order to analyze the strength of mud samples, without disturbing their structure, and before and after degradation of the organic matter.
Figure S5 shows the complex modulus and phase angle as a function of frequency for fresh and degraded mud samples for different ports. A solid-like behavior was observed for both fresh and degraded mud samples, i.e., an almost independence of complex modulus on frequency and significantly smaller values of phase angle [
24]. A similar solid-like character of mud samples has already been reported in the literature [
5,
25,
26,
27]. However, the organic matter degradation resulted in a weaker system (i.e., lower complex modulus values), particularly for PoA and PoH mud samples. The experimental data at higher frequencies were excluded from the graphs due to the inertial effects caused by the rheometer head.
Figure 6a shows the correlation between the change in complex modulus at 1 Hz (degraded–fresh) and the values of complex modulus at 1 Hz of fresh mud samples collected from different ports on a semi-log scale. A significant decrease in complex modulus was observed due to the degradation of organic matter, which became more pronounced for the fresh samples having higher complex modulus, particularly for PoA and PoH mud samples. For PoB and PoR mud samples, the complex modulus remains almost the same after the organic matter degradation. This behavior was also evident from the values of the fitting parameters (i.e., slope) of Equation (1) for the experimental data of complex modulus (
Table S4). The comparison between PoE and fluid mud samples from PoH for change in complex modulus at 1 Hz due to the organic matter degradation is presented in
Figure 6b. It is again evident that the mud samples from PoE displayed an increase in complex modulus values after organic matter degradation (i.e., positive slope), while in the case of PoH, a pronounced negative slope was observed with significant reduction in complex modulus (
Table S4).
The correlation between the change in phase angle at 1 Hz (degraded–fresh) and the values of phase angle at 1 Hz of fresh mud samples, obtained from different ports, was not prominent (
Figure S6). However, the degraded samples from PoA, PoB and PoH exhibited positive values of change in phase angle (i.e., higher phase angle values of degraded samples as compared to the fresh samples), which indicated the weakening of these mud samples. However, for mud samples from PoE and PoR, negative values of change in phase angle were observed, which highlighted the strengthening of mud samples.
The percent change in complex modulus
after organic matter degradation for mud samples from different ports is shown in
Figure 7. These results identified a significant reduction in complex modulus for PoA and PoH (~30–40% mean value), a negligible increase in complex modulus for PoB (~0–5% mean value), and a significant increase in complex modulus for PoE (~80% mean value). These results corresponded well to the outcome of yield stress and crossover amplitude. However, for mud samples from PoR, a negligible change in complex modulus was observed as compared to the changes in yield stresses and crossover amplitude. This behavior shows that the organic matter degradation strongly influenced the rheological property, which involves the destruction of structure (i.e., yield stress and crossover amplitude). While the effect of organic matter degradation on the rheological property, which is typically measured by having undisturbed structure (i.e., complex modulus within LVE regime), was negligible for the mud samples from PoR. This behavior may be attributed to the higher salinity (see
Table S1 for conductivity values) of these samples, which results in stronger flocculated structure and compensates the decrease in modulus values due to the organic matter degradation. However, this compensation offered by the salt-induced flocculation was negligible for yield stresses and crossover amplitude due to the application of higher shearing action for these rheological properties and the weaker nature of salt-induced flocs [
28].
3.6. Time Dependent and Structural Recovery Tests
Shear rate-controlled ramp-up and ramp-down experiments were performed from 0 to 100 s
−1, in order to investigate the time-dependent behavior of fresh and degraded mud samples from different ports. The results of the time-dependent experiments revealed the existence of a typical clockwise loop at higher shear rates for both fresh and degraded mud samples (
Figure S7). However, at lower shear rates, a counterclockwise loop was observed, which may be associated with a shear thickening phenomenon or a structural reorganization due to the shearing action [
3,
22]. A similar combination of clockwise and counterclockwise loops has already been reported in the literature for mud samples [
29].
Furthermore, the hysteresis loop was significantly lower for the degraded mud samples as compared to the fresh mud samples, which is again in line with the previous results. The hysteresis area between the clockwise loop (i.e., shaded region in
Figure S7) was estimated. The change in hysteresis area (degraded–fresh) as a function of the values of hysteresis area of fresh mud samples, for different ports, is presented in
Figure 8a. A similar decrease in hysteresis area was found as a function of organic matter degradation as already observed for other rheological properties, except for mud samples from PoB. This behavior is also evident from the values of fitting parameters of Equation (1) for time-dependent experiment (
Table S4). The mud samples from PoB displayed an increase in hysteresis area after the organic matter degradation, as also observed for PoE mud samples (
Figure 8b).
The percent change in hysteresis area
after organic matter degradation for mud samples from different ports is shown in
Figure 9. These results identified a significant reduction in hysteresis area for mud samples from PoA, PoH and PoR (~20–30% mean value) and a significant increase in hysteresis area for PoE mud samples (~210% mean value). These results again corresponded well to the outcome of yield stress and crossover amplitude. However, for mud samples from PoB, a substantial increase in hysteresis area (~50% mean value) was observed as compared to the changes in yield stresses and crossover amplitude. This behavior shows that the influence of organic matter degradation on the rheological property, which involves the destruction of structure (i.e., yield stress and crossover amplitude), was negligible for PoB mud samples. However, the organic matter degradation significantly affected the rheological property, which involves the measurement of a time-dependent response (i.e., stress/viscosity) of the mud samples from PoB. This behavior may be attributed to the different size and type of organic matter present in the mud samples from PoB, which needs further investigation.
In addition to the time-dependent characteristic of fresh and degraded mud samples, the structural recovery after intensive pre-shearing was investigated using a three-step protocol, as explained in
Section 2.4. The normalized time-dependent storage modulus (
) as a function of time for fresh and degraded mud samples is presented in
Figure S8. The oscillations in the storage modulus behavior as a function of time may be associated with the higher elasticity of the samples, which is typically stated as creep ringing [
19,
30]. A simple stretched exponential function, adapted from Mobuchon et al. [
31], was further used to fit the experimental data of the third step of structural recovery protocol, written as follows:
where
and
represent the two most important parameters of equilibrium storage modulus and characteristic time of the material, respectively.
is the time-dependent storage modulus of mud after shearing,
is the storage modulus before structural breakup,
is the storage modulus right after pre-shearing (at
t→0) and
is the stretching exponent. The change in normalized equilibrium storage modulus,
(degraded–fresh) as a function of the values of normalized equilibrium storage modulus of fresh mud samples is plotted in
Figure 10a for different ports. It can be seen that the values of normalized equilibrium storage modulus (
) are higher for the degraded mud samples than for the fresh mud samples (i.e., positive values of change in normalized equilibrium storage modulus). This may be attributed to the fact that the degradation of organic matter results in a weaker system, behaving as a purely mineral suspension without the bridging effect provided by organic matter, which eventually has a better structural recovery (i.e., higher values of modulus, higher strength) after pre-shearing. A similar enhanced structural recovery was observed in the literature for mud samples having lower organic matter content as compared to the samples having higher organic matter content [
19].
Furthermore,
Figure 10b shows a correlation between the change in characteristic time,
(degraded–fresh) and the values of characteristic time of fresh mud samples, for different ports. It is evident that for some samples the characteristic time was higher for degraded samples as compared to fresh mud samples (i.e., positive value of change), while for other samples it was lower for degraded samples than the fresh samples (i.e., negative value of change). This behavior may again be linked to the variable density of the mud samples, in addition to the organic matter degradation.