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 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 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 (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 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 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 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-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.

The concentrations and deposition of atmospheric constituents such as air pollutants were modeled with the Community Multiscale Air Quality (CMAQ) Model system for the year 2008. A focus was on sea salt emissions, atmospheric sea salt concentrations and the interaction between sea salt particles and atmospheric nitrogen compounds. For this purpose, two model runs with different realistic sea salt emission parameterizations were performed and are provided here. The model run CMAQ_CD24_seasalt_base uses the default sea salt emission parameterization of CMAQ denoted as GO03. The model run CMAQ_CD24_seasalt_ov14 uses an alternative sea salt emission parameterization denoted as OV14. Please refer to the connected references for information on the two sea salt emission parameterizations and their implementation. The dataset also contains meteorological input parameters in order to allow the user to perform unit conversions for some model variables. Additionally, the land-sea mask and land fraction data per model grid cell are provided. The projection is a Lambert Conformal Conic projection. The parameters are given in the netCDF files.