A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the years 1995 to 2014 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1984 and, then, from 1985 to 1994 with tagging (see below). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Riverine phosphorus input of the Warnow River was calculated with the Soil & Water Assessment Tool (SWAT; Bauwe et al., 2019, doi:10.1016/j.ecohyd.2019.03.003). Phosphorus from the Warnow River has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi:10.4319/lo.2006.51.1_part_2.0591). Therefore, all phosphorus-containing model variables exist twice in the output: once as regular variables and once as tagged variable. The phosphorus input by the Warnow River based on real phosphorus release patterns and real atmospheric conditions was calculated and a Maximum Technical Feasible Reduction (MTFR) approach was applied (PhosWaM SWAT case "35"). The turnover of phosphorus compounds in the Unterwarnow was calculated based on the "Unterwarnow turnover estimation v04" (see final project report of PhosWaM for details). The simulation was performed at the North-German Supercomputing Alliance (HLRN). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the years 1995 to 2014 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1984 and, then, from 1985 to 1994 with tagging (see below). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Riverine phosphorus input of the Warnow River was calculated with the Soil & Water Assessment Tool (SWAT; Bauwe et al., 2019, doi:10.1016/j.ecohyd.2019.03.003). Phosphorus from the Warnow River has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi:10.4319/lo.2006.51.1_part_2.0591). Therefore, all phosphorus-containing model variables exist twice in the output: once as regular variables and once as tagged variable. The default phosphorus input by the Warnow River based on real phosphorus release patterns and real atmospheric conditions was used ("base scenario"; PhosWaM SWAT case "ist"). The turnover of phosphorus compounds in the Unterwarnow was calculated based on the "Unterwarnow turnover estimation v04" (see final project report of PhosWaM for details). The simulation was performed at the North-German Supercomputing Alliance (HLRN). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed for the year 2012 covering the North Sea and Baltic Sea. Only data for the western Baltic Sea are provided here. The model output has been validated in Neumann et al. (2018a, doi: 10.5194/os-2018-71). The work was funded by the German Federal Ministry of Transport and Digital Infrastructure (BMVI, FKZ 50EW1601, https://www.io-warnemuende.de/meramo-en.html). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A). The model simulation was forced by operational meteorological data of the German Weather Service (DWD). Atmospheric nitrogen deposition data of 50x50 km2 spatial resolution were taken from the 2016 reporting of the European Measurement and Evaluation Programme (EMEP) as presented in EMEP (2016, url: http://emep.int/publ/reports/2016/EMEP_Status_Report_1_2016.pdf) and available from the Norwegian Meteorological Institute (2016, http://thredds.met.no/thredds/catalog/data/EMEP/2016_Reporting/catalog.html). Information on the riverine inputs, boundary conditions, and the model itself are provided in detail in Neumann et al. (2018b, doi: 10.5194/bg-2018-364). The concentrations of all prognostic biogeochemical model variables are given in nitrogen units according to the Redfield ratio.

A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the years 1995 to 2014 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1984 and, then, from 1985 to 1994 with tagging (see below). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Riverine phosphorus input of the Warnow River was calculated with the Soil & Water Assessment Tool (SWAT; Bauwe et al., 2019, doi:10.1016/j.ecohyd.2019.03.003). Phosphorus from the Warnow River has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi:10.4319/lo.2006.51.1_part_2.0591). Therefore, all phosphorus-containing model variables exist twice in the output: once as regular variables and once as tagged variable. The phosphorus input by the Warnow River based on real phosphorus release patterns and real atmospheric conditions was modified in order to comply with BASP (Baltic Sea Action Plan) targets (PhosWaM SWAT case "15"). The turnover of phosphorus compounds in the Unterwarnow was calculated based on the "Unterwarnow turnover estimation v04" (see final project report of PhosWaM for details). The simulation was performed at the North-German Supercomputing Alliance (HLRN). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the years 1995 to 2014 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1984 and, then, from 1985 to 1994 with tagging (see below). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Riverine phosphorus input of the Warnow River was calculated with the Soil & Water Assessment Tool (SWAT; Bauwe et al., 2019, doi:10.1016/j.ecohyd.2019.03.003). Phosphorus from the Warnow River has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi:10.4319/lo.2006.51.1_part_2.0591). Therefore, all phosphorus-containing model variables exist twice in the output: once as regular variables and once as tagged variable. The default phosphorus input by the Warnow River based on real phosphorus release patterns and real atmospheric conditions was used (PhosWaM SWAT case "ist"). The turnover of phosphorus compounds in the Unterwarnow was calculated based on the "Unterwarnow turnover estimation v05" (see final project report of PhosWaM for details). The simulation was performed at the North-German Supercomputing Alliance (HLRN). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed for the year 2012 covering the North Sea and Baltic Sea. Only data for the western Baltic Sea are provided here. The model output has been validated in Neumann et al. (2018a, doi: 10.5194/os-2018-71). The work was funded by the German Federal Ministry of Transport and Digital Infrastructure (BMVI, FKZ 50EW1601, https://www.io-warnemuende.de/meramo-en.html). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A). The model simulation was forced by operational meteorological data of the German Weather Service (DWD). Atmospheric nitrogen deposition data of medium spatial resolution of 16x16 km2 were provided by the Helmholtz-Zentrum Geesthacht within the EU BONUS SHEBA Project (Karl et al., 2019, doi: 10.5194/acp-2018-1317). Information on the riverine inputs, boundary conditions, and the model itself are provided in detail in Neumann et al. (2018b, doi: 10.5194/bg-2018-364). Nitrogen from atmospheric deposition of shipping-related nitrogen, agricultural-related nitrogen, and total nitrogen has been tagged in the model simulation according to a method by Menésguen et al. (2006, 10.4319/lo.2006.51.1_part_2.0591). Therefore, all nitrogen-containing model variables exist four times in the output: once as regular variables and once per tagged nitrogen source (total, shipping-related, agricultural-related). The concentrations of all prognostic biogeochemical model variables are given in nitrogen units according to the Redfield ratio.

A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the years 2003 to 2014 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1999 and, then, from 2000 to 2002 with tagging (see below). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). A publication is in preparation. The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Atmospheric nitrogen deposition data of 0.1° x 0.1° spatial resolution were taken from the 2018 reporting of the European Measurement and Evaluation Programme (EMEP) as presented in EMEP (2018, url: http://emep.int/publ/reports/2018/EMEP_Status_Report_1_2018.pdf) and available from the Norwegian Meteorological Institute (2018, url: http://thredds.met.no/thredds/catalog/data/EMEP/2018_Reporting/catalog.html). Nitrogen from atmospheric deposition of nitrogen from livestock/agricultural emissions (estimated, see documentation) and from all emission sectors has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi:10.4319/lo.2006.51.1_part_2.0591). Therefore, all nitrogen-containing model variables exist three times in the output: once as regular variables and once per tagged nitrogen source (total atmospheric and agriculturally-related). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed with the model MOM-ERGOM for the year 2012 covering the Baltic Sea. Previously, MOM-ERGOM had been initialized for several decades without tagging until 1999 and, then, from 2000 to 2011 with tagging (see below; three years would have been sufficient). The model output has been validated with measurement data of the "IOW Baltic Monitoring and long-term data program" (https://www.io-warnemuende.de/iowdb.html; IOW: Leibniz Institute for Baltic Sea Research Warnemünde) and from the HELCOM database (http://ocean.ices.dk/helcom/Helcom.aspx; HELCOM: Helsinki Commission). A publication is in preparation. The model simulation was forced by coastDat2 COSMO-CLM data (doi:10.1594/WDCC/coastDat-2_COSMO-CLM). Atmospheric nitrogen deposition data of 16x16 km2 horizontal resolution were provided by the Helmholtz-Zentrum Geesthacht within the EU BONUS SHEBA Project (Karl et al., 2019, doi:10.5194/acp-19-7019-2019). Nitrogen from atmospheric deposition of nitrogen from shipping emissions and from all emission sectors has been tagged in the model simulation according to a method by Menésguen et al. (2006, doi: 10.4319/lo.2006.51.1_part_2.0591). Therefore, all nitrogen-containing model variables exist three times in the output: once as regular variables and once per tagged nitrogen source (total atmospheric and shipping-related). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A).

A marine physical biogeochemical model simulation was performed for the year 2012 covering the North Sea and Baltic Sea. Only data for the western Baltic Sea are provided here. The model output has been validated in Neumann et al. (2018a, doi: 10.5194/os-2018-71). The work was funded by the German Federal Ministry of Transport and Digital Infrastructure (BMVI, FKZ 50EW1601, https://www.io-warnemuende.de/meramo-en.html). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A). The model simulation was forced by operational meteorological data of the German Weather Service (DWD). Atmospheric nitrogen deposition data of high spatial resolution of 4x4 km2 were provided by the Helmholtz-Zentrum Geesthacht within the EU BONUS SHEBA Project (Karl et al., 2019, doi: 10.5194/acp-2018-1317). Information on the riverine inputs, boundary conditions, and the model itself are provided in detail in Neumann et al. (2018b, doi: 10.5194/bg-2018-364). Nitrogen from atmospheric deposition of shipping-related nitrogen has been tagged in the model simulation according to a method by Menésguen et al. (2006, 10.4319/lo.2006.51.1_part_2.0591). Therefore, all nitrogen-containing model variables exist twice in the output: once as regular variables and once as nitrogen content from shipping-related activities. The concentrations of all prognostic biogeochemical model variables are given in nitrogen units according to the Redfield ratio.

A marine physical biogeochemical model simulation was performed for the year 2012 covering the North Sea and Baltic Sea. Only data for the western Baltic Sea are provided here. The model output has been validated in Neumann et al. (2018a, doi: 10.5194/os-2018-71). The work was funded by the German Federal Ministry of Transport and Digital Infrastructure (BMVI, FKZ 50EW1601, https://www.io-warnemuende.de/meramo.html). The simulation was performed at the North-German Supercomputing Alliance (HLRN, project id: mvk00054, zulassung.hlrn.de/kurzbeschreibungen/mvk00054.pdf). The model output data were processed and evaluated on servers provided by the project 'PROSO - Prozesse von Spurenstoffen in der Ostsee' (FKZ 03F0779A). The model simulation was forced by operational meteorological data of the German Weather Service (DWD). Atmospheric nitrogen deposition data of medium spatial resolution of 16x16 km2 were provided by the Helmholtz-Zentrum Geesthacht within the EU BONUS SHEBA Project (Karl et al., 2019, doi: 10.5194/acp-2018-1317). Information on the riverine inputs, boundary conditions, and the model itself are provided in detail in Neumann et al. (2018b, doi: 10.5194/bg-2018-364). Nitrogen from atmospheric deposition of shipping-related nitrogen, agricultural-related nitrogen, and total nitrogen has been tagged in the model simulation according to a method by Menésguen et al. (2006, 10.4319/lo.2006.51.1_part_2.0591). Therefore, all nitrogen-containing model variables exist four times in the output: once as regular variables and once per tagged nitrogen source (total, shipping-related, agricultural-related). The concentrations of all prognostic biogeochemical model variables are given in nitrogen units according to the Redfield ratio.