Silicic acid flux to the ocean from tidal permeable sediments: A modeling study
Sandy sediments of tidal beaches are poor in reactive substances because they are regularly flushed by significant flow caused by tidal forcing. This transport process may significantly affect the flux of reactive solutes to the ocean. A two dimensional model coupling the Richards equation that describes the flow in permeable sediments and the conservation equation of the silicic acid was developed to simulate the evolution of the silicic acid concentration into a variably saturated porous media submitted to tidal forcing. A detailed algorithm of drainage zone under tidal forcing and numerical methods needed to solve it are properly presented. Flux to the ocean has been estimated. The silicic acid concentration displays a permanent lens with low silicic acid concentration at the top of the tidal zone. This lens that results from the tidal forcing, presents weak variations of area during the tidal cycle. Silicic outflux to the ocean increases with increasing beach slope, hydraulic conductivity and tidal range. Simulations reveal that the total silicic acid flux to the ocean from the coastal marine sands can be considered as significant compared to the flux supplied by the rivers. These results may alter the previously published global budget of the silicic acid to the ocean.
silicic acid tidal forcing richards equation beach sands unsaturated porous media coastal environment