Advection-Diffusion Model sediments

Reaction rates of nonconservative chemicals in marine sediments can be estimated from porewater concentration profiles using a mathematical model similar to the onedimensional advection-diffusion model for the water column presented in Section 4.3.4. As with the water column, horizontal concentration gradients are assumed to be negligible as compared to the vertical gradients. In contrast to the water column, solute transport in the pore waters is controlled by molecular diffusion and advection, with the effects of turbulent mixing being negligible. 

Magnesium concentrations as a function of depth (meters below the sea floor) in sediment porewaters from the western flank of the Juan de Fuca Ridge near 48° N in the North Pacific Ocean. The concentration decreases with depth because it is removed from solution by reaction with crustal rocks at the sediment-crustal boundary. The curves are convex upward because of porewater upwelling along the upward-flowing limb of a convection cell. Velocities of the upwelling are determined by using a one-dimensional advection-diffusion model and are indicated by the numbers on the curves. Redrafted from Wheat and MottI (2000). 

The standard method of estimating chemical migration in a cap is via a transient advection-diffusion model as described by Palermo et al. [1]. This model is applied to the chemical isolation layer of a cap, which is the cap thickness after removing components for porewater expression via consolidation of underlying sediment, consolidation of the cap, and bioturbation of the upper cap layers. Normally, an analytical solution to the mass conservation equation, assuming that the cap is semi-infinite, is employed in such an... 

In addition to dissipation of the substance from the model system through degradation, other dissipative mechanisms can be considered. Neely and Mackay(26) and Mackay(3) have also introduced advection (loss of the chemical from the troposphere via diffusion) and sedimentation (loss of the chemical from dynamic regions of the system by movement deep into sedimentation layers). Both of these mechanisms are then assumed to act in the unit world. This approach makes it possible to investigate the behavior of atmosphere emissions where advection can be a significant process. Therefore, from a regulatory standpoint if the emission rate exceeds the advection rate and degradation processes in a system, accumulation of material could be expected. Based on such an analysis reduction of emissions would be called for. 

An alternative to the one-dimensional model can be developed by noting that the zone of sediment inhabited by macrofauna is not a homogeneous one-dimensional slab, but instead is a body permeated by cylinders. The water within these cylinders (burrows) is maintained at approximately seawater solute concentrations by the irrigation activity of their animal occupants. Interstitial solutes can therefore diffuse into burrows and be advected out of sediment by irrigation activity as well as diffuse vertically toward the sediment-water interface. Diffusion in this case can be considered as taking place in a system of cylindrical symmetry similar to that occurring in a root-permeated soil (Gardner, 1980 Cowan, 1%5 Nye and Tinker, 1977). 

Ra distribution in the ocean has been modeled to derive eddy diffusivities and advection rates taking into consideration its input by diffusion from sediments, loss by radioactive decay, and dispersion... 

A conceptual model of sedimentary nitrogen cycling. Dashed arrows represent pore water diffusion and advection. Dotted arrows represent sedimentation. Source-. After Burdige, D.J. (2006). Geochemistry of Marine Sediments. Princeton University Press, p. 453. 

The model (Fig. 23.6) consists of three compartments, (a) the surface mixed water layer (SMWL) or epilimnion, (b) the remaining open water column (OP), and (c) the surface mixed sediment layer (SMSL). SMWL and OP are assumed to be completely mixed their mass balance equations correspond to the expressions derived in Box 23.1, although the different terms are not necessarily linear. The open water column is modeled as a spatially continuous system described by a diffusion/advection/ reaction... 

Park and Ortoleva (2003) have developed WRIS.TEQ, a comprehensive reaction-transport-mechanical simulator that includes kinetic and thermodynamic properties with mass transport (advection and diffusion). A unique property of this code is a dynamic compositional and textural model specifically designed for sediment alteration during diagenesis. 

The interpretation of pore-water concentration versus depth profiles of O2 and NO in oxic sediments is based on a one-dimensional, steady-state model in which the production or consumption of a solute in a sedimentary layer is balanced by transport into or out of the layer by solute diffusion and burial advection. In mathematical form. 

In water and sediments, the time to chemical steady-states is controlled by the magnitude of transport mechanisms (diffusion, advection), transport distances, and reaction rates of chemical species. When advection (water flow, rate of sedimentation) is weak, diffusion controls the solute dispersal and, hence, the time to steady-state. Models of transient and stationary states include transport of conservative chemical species in two- and three-layer lakes, transport of salt between brine layers in the Dead Sea, oxygen and radium-226 in the oceanic water column, and reacting and conservative species in sediment. 

If there are no or very few irrigated burrows present in the sediment, lateral diffusion is not significant and the r dependence of Eq. (6.12) can be ignored. In that case, the equation becomes the more traditional onedimensional transport-reaction equation used to model pore-water solute profiles where advection is relatively unimportant (Berner, 1971 1980 Lerman, 1979). Both the cylindrical microenvironment model and the onedimensional Cartesian coordinate model will be used here to quantify the Mn distributions at NWC and DEEP. 

Fii ure 9. (ntluence of cJislribuiion coefficient on the contributions of diffusive and advective fluxes to the total sediment flux of an element as predicted by the model of Diamond et al. (1990) for a lake enclosure. 

---