README for UVic simulation with riverine nutrient export from NEWS2 dataset Miriam Tivig, David P. Keller and Andreas Oschlies : "Riverine nutrient impact on global ocean nitrogen cycle feedbacks and marine primary production in an Earth System Model" Submitted at Biogeosciences Corresponding author: Miriam Tivig (mtivig@geomar, miriam.tivig@dwd.de) January 2024 ---------------------------- This folder contains a folder with the updates for the model code of the University of Victoria Earth System Climate Model UVic Version 2.9 (Keller et al., 2012; Eby et al., 2009; Weaver et al.,2001) including a module for river nutrient export, the compiled code for 7 simulations and the corresponding output, used for the study of riverine nitrogen export to the coastal ocean. Update_MTivig_2020_BP ----------------- This folder contains the updated code for the simulations. The main changes to the code of Keller et al. (2012) have been made in: - /common/gosbc.F: Line 354 - Line 390 NO3 and / or dissolved inorganic and organic and 45 % of total particulate phosphorus is added to the vertical fluxes (ino3flx, idipfly, idopflx, ippflx) - /common/setdata.f: L. 42, & 142, the input file with nutrients from NEWS is read by the model. Depending on the simulation, a different data-file was used: L_NEWS_seas.nc (global NEWS simulation with riverine nitrogen), L_NEWS_seas_with_PP.nc ((global NEWS simulation with riverine phosphorus). - /mom/tracer.F : - in the npzd model, benthic denitrification and subgrid bathymetry are defined. - L. 51, L. 600 ff: definition of the burial functions (O_benthic_P_burial_low, O_benthic_P_burial_low) - /mom/mom_tavg.F: benthic phosphate burial is defined - /embm/rivmodel.F : line 200 ff the rivermodel has been changed to read in NEWS nutrients without seasonality Code_NEWS_N ----------------- The code-folder contains the compiled code of the simulation with UVic including NEWS river export of dissolved inorganic nitrogen (DIN) from all river discharge points. Code_CTR ----------------- Compiled code of the control simulation, same configuration as NEWS, but without river DIN discharge. Code_run_NEWS_P ----------------- The code-folder contains the compiled code of the simulation with UVic including NEWS river export of phosphorus (dissolved organic and inorganic, 45 % of particulate phosphorus) from all river discharge points. Code_run_NEWS_N_P ----------------- Compiled code of the simulation, same configuration as in NEWS_N but including nitrogen and phosphorus export from rivers (without phosphorus burial functions). Code_run_NEWS_N_P_burlow ----------------- Compiled code of the simulation, same configuration as in NEWS_N_P but including burial functions for low burial rates. Code_run_NEWS_N_P_burhigh ----------------- Compiled code of the simulation, same configuration as in NEWS_N_P but including burial functions for high burial rates. data.100.100.19 ---------------- Contains the input data files, including the riverine nutrient export files for the simulations NEWS (L_NEWS_seas.nc, L_NEWS_seas_with_PP.nc). Includes also the restart.nc used for the simulations. Note that the data that has been used for these input files provides from the Global NEWS2 model, which was developed by the Global NEWS group and is available at the web site: http://www.ioc-unesco.org/icr/index.php?option=com_content&view=article&id=45&Itemid=100002 . Please email Emilio Mayorga at mayorga@marine.rutgers.edu to obtain this dataset. We are very thankful to the Global NEWS group of making the use of this dataset possible. Files for compilation ----------------- Control.in, mk.in OUTPUT ----------------- Model-output of the 6 simulations (CTR, NEWS_N, NEWS_P, NEWS_N_P, NEWS_N_P_burlow, NEWS_N_P_burhigh) performed with the model configuration described before, used for the Figures and results in the manuscript. Each variables is shown as 100 years average for each grid point. - UVIC29-[simulation]_100y_10000y_tavg_all.nc : global ocean data of the simulation averaged over 100 years of all variables. - UVIC29_[simulation]_100ym_10000y_O_no3_ATL_zm.nc: zonal mean of marine nitrogen concentrations in the Atlantic Ocean over 10000 years averaged over 100 years. - UVIC29_[simulation]_100ym_10000y_O_no3_PAC_zm.nc: zonal mean of marine nitrogen concentrations in the Pacific Ocean over 10000 years averaged over 100 years. - UVIC29_[simulation]_100ym_10000y_O_no3_IND_zm.nc: zonal mean of marine nitrogen concentrations in the Indian Ocean over 10000 years averaged over 100 years. - UVIC29_[simulation]_100ym_10000y_O_po4_ATL_zm.nc: zonal mean of marine phosphorus concentrations in the Atlantic Ocean over 10000 years averaged over 100 years. - UVIC29_[simulation]_100ym_10000y_O_po4_PAC_zm.nc: zonal mean of marine phosphorus concentrations in the Pacific Ocean over 10000 years averaged over 100 years. - UVIC29_[simulation]_100ym_10000y_O_po4_IND_zm.nc: zonal mean of marine phosphorus concentrations in the Indian Ocean over 10000 years averaged over 100 years. The Ocean basin data are zonal means of the variables for the three main ocean basins, Atlantic Ocean, Pacific Ocean and Indian Ocean. The ocean basins have been selected with climate data operators (cdo; Schulzweida, 2019). - ATL: 80°W-20°E 80°S-80°N - IND: 30°E-120°E 50°S-25°N - PAC: South-Pac. 130°E-74°W 60°S-0 North-Pac. 115°E-96°W 0-60°N RELEVANT PUBLICATIONS --------------------- Bohlen, L., Dale, A. W., & Wallmann, K. (2012): Simple transfer functions for calculating benthic fixed nitrogen losses and c:n:p regeneration ratios in global biogeochemical models., 26, GB3029. doi: 10.1029/2011GB004198 Eby, M., Zickfeld, K., Montenegro, A., Archer, D., Meissner, K. J., & Weaver, A. J. (2009): Lifetime of anthropogenic climate change: Millennial time scales of potential co2 and surface temperature perturbations., 22, 2501-2511. doi:10.1175/2008JCLI2554.1 Keller, D. P., Oschlies, A., & Eby, M. (2012). A new marine ecosystem model for the university of victoria earth system climate model., 5, 1195-1220. doi:10.5194/gmd-5-1195-2012 Mayorga, E., Seitzinger, S. P., Harrison, J. A., Dumont, E., Beusen, A. H. W., Bouwman,A.F., Fekete, B., Kroeze, C. and Van Drecht, G. (2010). Global nutrient export fromwater sheds 2 (NEWS 2): Model development and implementation. Environmental Modelling and Software, 25, 83-{853. doi:10.1016/j.envsoft.2010.01.007 Schmittner, A., Oschlies, A., Giraud, X., Eby, M., & Simmons, H. L. (2005): A global model of the marine ecosystem for long-term simulations: Sensitivity to ocean mixing buoyancy forcing, particle sinking, and dissolved organic matter cycling. , 19, GB3004. doi: 10.1029/2004GB002283 Schulzweida, Uwe. (2019, October 31): CDO User Guide (Version 1.9.8). http://doi.org/10.5281/zenodo.3539275 Seitzinger, S. P., Harrison, J. A., Dumont, E., Beusen, A. H. W., & Bouwman, A. F. (2005): Sources and delivery of carbon, nitrogen, and phosphorus to the coastal zone: An overview of global nutrient export from watersheds (NEWS) models and their application., 19, GB4S01. doi: 10.1029/2005GB002606 Tivig, M., Keller D. P., and Oschlies A. (2021): Feedbacks in the marine nitrogen cycle limit the impact of riverine nitrogen supply on global marine biology and biogeochemistry in an Earth System Model. Biogeosciences, 18, 5327–5350, https://doi.org/10.5194/bg-18-5327-2021