Description: We present a new nitrogen isotope model incorporated into the three-dimensional ocean component of a global Earth system climate model designed for millennial timescale simulations. The model includes prognostic tracers for the two stable nitrogen isotopes, 14N and 15N, in the nitrate (NO3−), phytoplankton, zooplankton, and detritus variables of the marine ecosystem model. The isotope effects of algal NO3− uptake, nitrogen fixation, water column denitrification, and zooplankton excretion are considered as well as the removal of NO3− by sedimentary denitrification. A global database of δ15NO3− observations is compiled from previous studies and compared to the model results on a regional basis where sufficient observations exist. The model is able to qualitatively and quantitatively reproduce many of the observed patterns such as high subsurface values in water column denitrification zones and the meridional and vertical gradients in the Southern Ocean. The observed pronounced subsurface minimum in the Atlantic is underestimated by the model presumably owing to too little simulated nitrogen fixation there. Sensitivity experiments reveal that algal NO3− uptake, nitrogen fixation, and water column denitrification have the strongest effects on the simulated distribution of nitrogen isotopes, whereas the effect from zooplankton excretion is weaker. Both water column and sedimentary denitrification also have important indirect effects on the nitrogen isotope distribution by reducing the fixed nitrogen inventory, which creates an ecological niche for nitrogen fixers and, thus, stimulates additional N2 fixation in the model. Important model deficiencies are identified, and strategies for future improvement and possibilities for model application are outlined.

Description: Highlights • We present a 5 myr record of biogeochemical cycling in a Cretaceous upwelling area. • A novel quantitative approach for the evaluation of Fe speciation proxies was applied. • Ferruginous proxy signature reflects intense chemical weathering rather than anoxia. • Water column redox conditions evolved from oxic to nitrogenous to euxinic before OAE2. • Smaller seawater nitrate inventory facilitated sedimentary H2S release and euxinia. Abstract Oceanic Anoxic Events (OAEs) in Earth's history are regarded as analogues for current and future ocean deoxygenation, potentially providing information on its pacing and internal dynamics. In order to predict the Earth system's response to changes in greenhouse gas concentrations and radiative forcing, a sound understanding of how biogeochemical cycling differs in modern and ancient marine environments is required. Here, we report proxy records for iron (Fe), sulfur and nitrogen cycling in the Tarfaya upwelling system in the Cretaceous Proto-North Atlantic before, during and after OAE2 (∼93 Ma). We apply a novel quantitative approach to sedimentary Fe speciation, which takes into account the influence of terrigenous weathering and sedimentation as well as authigenic Fe (non-terrigenous, precipitated onsite) rain rates on Fe-based paleo-redox proxies. Generally elevated ratios of reactive Fe (i.e., bound to oxide, carbonate and sulfide minerals) to total Fe (FeHR/FeT) throughout the 5 million year record are attributed to transport-limited chemical weathering under greenhouse climate conditions. Trace metal and nitrogen isotope systematics indicate a step-wise transition from oxic to nitrogenous to euxinic conditions over several million years prior to OAE2. Taking into consideration the low terrigenous sedimentation rates in the Tarfaya Basin, we demonstrate that highly elevated FeHR/FeT from the mid-Cenomanian through OAE2 were generated with a relatively small flux of additional authigenic Fe. Evaluation of mass accumulation rates of reactive Fe in conjunction with the extent of pyritization of reactive Fe reveals that authigenic Fe and sulfide precipitation rates in the Tarfaya Basin were similar to those in modern upwelling systems. Because of a smaller seawater nitrate inventory, however, chemolithoautotrophic sulfide oxidation with nitrate was less efficient in preventing hydrogen sulfide release into the water column. As terrigenous weathering and sediment flux determine how much authigenic Fe is required to generate an anoxic euxinic or ferruginous proxy signature, we emphasize that both have to be taken into account when interpreting Fe-based paleo-redox proxies.