Supplemental material - Nutrients as the dominant control on the spread of anoxia and euxinia across the Cenomanian-Turonian oceanic anoxic event (OAE2): Model-data comparison Fanny M. Monteiro, Richard Pancost, Andy Ridgwell and Yannick Donnadieu 1 Additional model information 1.1 Accessing model source code and experimental results The model code for the version of the GENIE model used in this paper (technically: ‘cGE- NIE´) can be found here: https://svn.ggy.bris.ac.uk/subversion/genie/tags/cgenie-Monteiro.PO.2012 (svn revision 7679) and includes all configuration and boundary condition files needed to carry out the spin-ups, the control experiments, and all parameter variation experiments (check genie-userconfigs/MS/PO2012.Monteiro for the specific userconfig files). The code on the SVN repository can be accessed with the username: genie-user and password: g3n1e-user Documentation on running the cGENIE model can be found in the genie-docs directory of the code installation. 1.2 Marine productivity The activity of the phytoplankton is modeled via a biological net depletion of nutrient by non-fixing phytoplankton (Γ P hy ) and diazotrophs (Γ Diaz ): Γ P hy =V P hy max γ T γ I min  PO -3 4 PO -3 4 + K P ; DIN DIN + K N  min  PO -3 4 ; DIN 16  (1.2.1) Γ Diaz =V Diaz max γ T γ I PO -3 4 PO -3 4 + K P PO -3 4 if DIN < 16PO -3 4 and DIN <N T hresh (1.2.2) V P hy max and V Diaz max is the maximum net depletion rate of non-fixing phytoplankton and dia- zotrophs respectively. Parameter values are presented in Table S1. The biomass of phy- toplankton in the model is simply related to the concentration of the limiting nutrient (Maier-Reimer , 1993; Wallmann , 2003; Doney et al., 2006). This scheme is used to improve computational time needed to run the model for at least 10 kyrs, but still compares well to observations of productivity in the modern ocean. γ T is the temperature function based 1