Modeling estuarine nutrient geochemistry in a simple system

We present a model of estuarine mixing, removal, and input for dissolved constituents, and apply the model to 39 nutrient (P, N, Si) profiles collected over a 14-month period in a pristine river/ estuary: Ochlockonee Bay, Florida. Each profile is deconvolved into three component functions: linear mixing (conservative) first-order removal (biological productivity), and parabolic input (regeneration). After correction for temporal variations in the fluvial end-members, the model provides quantitative estimates of total estuarine primary production, net regeneration, and subsequent fluxes to the ocean over a year-long period. The modeled data set is internally self-consistent: virtually perfect mass balances are obtained for P and Si. All biological P-uptake is regenerated within the estuary so that virtually 100% of the fluvial reactive-P enters the ocean. One-third of the fluvial reactive-P enters the estuary as particles whose phosphate is released after deposition in estuarine sediments. About 20% of the dissolved fluvial silica flux is removed biologically; all of this biogenic silica dissolves in the estuary and enters the ocean. N cannot be mass balanced, probably because it enters and escapes the bay in unmeasured forms (as NH₄ or via denitrification to N₂ and N₂O). In the Ochlockonee, biological productivity removes nutrients in the ratios N:P ? 9:1 and Si:P ? 20:1.