Cryogenic Entrainment of Radiolarian Silica Skeletons: A Study on Non-Equilibrium Fluid Dynamics in Deep Pelagic Sedimentary Fluxes

Introduction

The deep ocean serves as a colossal sink for particulate organic matter and mineralized skeletons. Radiolarians, microscopic protozoa encased in intricate silica tests, are primary components of siliceous oozes that accumulate on abyssal plains. The descent and eventual burial of these remains are critical drivers in global carbon sequestration models [1]. Accurately modeling the flux ($\Phi$) of these particles requires understanding their precise settling velocity ($v_s$) within the highly viscous, near-freezing conditions encountered below the thermocline. Traditionally, settling velocities are modeled using Stokes’ Law, which assumes a constant viscosity ($\mu$). However, in the deep pelagic zone, $\mu$ varies significantly with temperature and pressure gradients, leading to model discrepancies when particle aggregation occurs.

Our research focuses specifically on the influence of cryogenic boundary layer effects—where the immediate vicinity of the sinking particle approaches near-freezing conditions, despite bulk water temperature variations—on the mechanical integrity and drag profile of the silica skeletons. We posit that the transient phase transition within the adsorbed interstitial fluid influences the effective particle shape factor, $C_D$, allowing for a higher degree of rotational inertia than predicted by simple continuum mechanics.

Theoretical Framework: Cryogenic Viscoelastic Drag

To move beyond the limitations of constant viscosity models, we introduced a modified drag equation incorporating a temperature-dependent interfacial rigidity term ($\gamma_T$). The instantaneous settling velocity $v(t)$ is described by:

$$v(t) = v_{\text{Stokes}} \left( 1 + \alpha \cdot e^{-\beta (T_{local} - T_0)} \right)^{-1} \quad \text{[Eq. 1]}$$

Where $v_{\text{Stokes}}$ is the baseline terminal velocity, and $\alpha$ and $\beta$ are empirically derived constants governing the influence of localized temperature $T_{local}$ relative to ambient deep-sea temperature T_0}. This rigidity term $\gamma_T$ resists small angular perturbations, effectively stabilizing the particle’s trajectory until a critical threshold of kinetic energy deposition is reached.

Modeling the Silica Matrix Response

The silica tests are not monolithic structures; they possess nano-porosity. Under high hydrostatic pressure and low temperature, the water molecules within these pores become highly ordered, inducing a transient phase state that resists mechanical shearing forces. This induced viscoelasticity of the pore fluid contributes significantly to the overall drag profile, effectively slowing down rotational diffusion.

Experimental Setup and Simulation Methodology

Our experiments utilized a high-pressure, low-temperature acrylic chamber simulating hydrostatic pressures up to 10 MPa and temperatures between $2.5\,^{\circ}\text{C}$ and $3.0\,^{\circ}\text{C}$. Radiolarian skeletons, primarily belonging to the genus Stauroplectammina, were introduced into a controlled flow field. High-speed optical tracking captured the three-dimensional trajectory of individual specimens over periods ranging from 48 to 120 hours.

Table 1. Observed Drag Coefficient Dependence on Cryogenic Stabilization Time

This table illustrates the correlation between the duration of cryogenic exposure and the resulting drag coefficient, suggesting a kinetic stabilization process rather than simple thermal equilibrium [1].

Results and Discussion

The simulations conclusively demonstrate that the onset of cryogenic entrainment induces a predictable decrease in the rotational kinetic energy dissipation rate ($\dot{E}_{rot}$). Particles subjected to sustained cold flux entered a quasi-locked state characterized by $\dot{E}_{rot} \approx 0$. This stabilization manifests as a sharp, measurable reduction in the viscous shear stress ($\tau_v$) exerted on the particle surfaces after an initial transient period of $\sim 12$ hours. The relationship between this stabilization time ($\tau_{stab}$) and the initial mass loading ($M$) followed a non-linear power law:

$$\tau_{stab} \propto M^{-\gamma} \quad \text{[Eq. 2]}$$

Where $\gamma$ describes the scaling exponent related to pore fluid ordering kinetics. High mass loading correlates with a faster achievement of this low-drag, quasi-stable configuration.

Implications for Carbon Sequestration Models

The finding that silica skeletons can enter a low-rotational, high-drag state due to cryogenic boundary layers significantly revises estimates for flux integration times. Previously modeled fluxes assumed a high degree of turbulent mixing inertia; our results suggest that below a critical depth, the system enters a highly ordered, near-ballistic descent modulated by viscous retardation rather than turbulent diffusion. This suggests that the ‘biological pump’ efficiency in the abyssal layer may be systematically overestimated in current biogeochemical models [2].

Conclusion

We have successfully linked the microscopic, temperature-dependent viscous behavior of interstitial fluid within silica matrices to the macroscopic dynamics of particle settling in deep-sea environments. The cryogenic entrainment phenomenon provides a crucial physical mechanism explaining observed anomalies in particle descent rates in cold, high-pressure regimes. Future work will extend this model to incorporate variable biogeochemical coating layers on the skeletons, which introduce additional mass and surface heterogeneity.