Non-Fickian Diffusion and the Accumulation of Methane Bubbles in Deep-Water Sediments

In the absence of fractures, methane bubbles in deep-water sediments can be immovably trapped within a porous matrix by surface tension. The dominant mechanism of transfer of gas mass therefore becomes the diffusion of gas molecules through porewater. The accurate description of this process requires non-Fickian diffusion to be accounted for, including both thermodiffusion and gravitational action. We evaluate the diffusive flux of aqueous methane considering non-Fickian diffusion and predict the existence of extensive bubble mass accumulation zones within deep-water sediments. The limitation on the hydrate deposit capacity is revealed; too weak deposits cannot reach the base of the hydrate stability zone and form any bubbly horizon.

💡 Research Summary

The paper investigates how methane bubbles behave in deep‑water sediments where fractures are absent and surface tension immobilizes bubbles within the porous matrix. In such environments, the dominant transport mechanism for methane is not advective migration of bubbles but diffusion of dissolved methane through the pore water. The authors argue that a simple Fickian description of diffusion is insufficient; instead, a non‑Fickian framework that incorporates thermodiffusion (the Soret effect) and gravitational segregation must be employed.

The authors derive a generalized diffusion flux equation:

( \mathbf{J}= -D \nabla C - D_T C \nabla T + \rho g \beta C )

where (D) is the ordinary molecular diffusion coefficient, (D_T) the thermodiffusion coefficient, (C) the dissolved methane concentration, (\rho) the density of pore water, (g) the acceleration due to gravity, and (\beta) a coefficient describing the buoyancy effect of dissolved methane. The first term represents classical concentration‑gradient diffusion, the second term captures the tendency of methane molecules to migrate toward warmer regions (thermodiffusion), and the third term accounts for the downward pull exerted by gravity on the slightly heavier methane‑laden water.

Using realistic temperature and pressure gradients typical of continental slope sediments, the authors evaluate the relative magnitudes of the thermodiffusive and gravitational contributions. They find that in regions where the temperature gradient is steep, the thermodiffusive term can dominate, driving methane upward toward the base of the hydrate stability zone (HSZ). Conversely, where the temperature gradient is modest and the pressure gradient is significant, the gravitational term prevails, causing methane to accumulate in deeper layers beneath the HSZ.

A one‑dimensional vertical transport model is constructed, incorporating the derived flux expression, realistic boundary conditions (no methane flux at the sediment surface, fixed concentration at the deep boundary), and initial uniform methane saturation. Numerical integration reveals two distinct accumulation zones:

1. Upper Accumulation Zone (near the top of the HSZ): Here the Soret effect pushes methane upward, leading to a thin layer of small bubbles that remain trapped by surface tension. This zone is relatively narrow but can act as a source of free gas if the overlying pressure drops.

2. Lower Accumulation Zone (below the HSZ): In this deeper region the gravitational term dominates, causing dissolved methane to migrate downward and form a broad, persistent “bubble horizon.” Because the bubbles are small enough to be held by capillary forces, they can persist for geological timescales, representing a substantial, previously unrecognized methane reservoir.

The authors further explore the capacity limit of hydrate deposits. If the HSZ is thin or the temperature‑pressure gradients are too steep, methane cannot penetrate to the deeper accumulation zone; instead it either dissolves back into the pore water or escapes upward. Thus, a sufficiently thick HSZ and gentle gradients are prerequisites for the formation of extensive deep‑water bubble horizons.

Key implications of the study include:

• Reassessment of methane budgets: Traditional models that ignore non‑Fickian effects may underestimate the amount of methane stored in deep sediments.
• Hazard evaluation: The deep bubble horizon could become destabilized by seismic activity or warming, potentially releasing large volumes of methane.
• Resource potential: The identified accumulation zones may represent viable targets for future methane hydrate extraction, provided that engineering approaches can address the capillary‑trapped bubbles.

In conclusion, the paper demonstrates that non‑Fickian diffusion, specifically the interplay between thermodiffusion and gravity, fundamentally controls methane bubble distribution in deep‑water sediments. Incorporating these mechanisms into predictive models is essential for accurate assessments of marine methane hydrate resources, their stability, and their impact on climate and geohazards.