Long-term oxygen and nutrient transports in the Baltic Sea are reconstructed using the Swedish Coastal and Ocean Biogeochemical model (SCOBI) coupled to the Rossby Centre Ocean model (RCO). Two simulations with and without data assimilation covering the period 1970–1999 are carried out. Here, the “weakly coupled” scheme with the Ensemble Optimal Interpolation (EnOI) method is adopted to assimilate observed profiles in the reanalysis system. The reanalysis shows considerable improvement in the simulation of both oxygen and nutrient concentrations relative to the free run. Further, the results suggest that the assimilation of biogeochemical observations has a significant effect on the simulation of the oxygen-dependent dynamics of biogeochemical cycles. From the reanalysis, nutrient transports between sub-basins, between the coastal zone and the open sea, and across latitudinal and longitudinal cross sections are calculated. Further, the spatial distributions of regions with nutrient import or export are examined. Our results emphasize the important role of the Baltic proper for the entire Baltic Sea, with large net transport (export minus import) of nutrients from the Baltic proper into the surrounding sub-basins (except the net phosphorus import from the Gulf of Riga and the net nitrogen import from the Gulf of Riga and Danish Straits). In agreement with previous studies, we found that the Bothnian Sea imports large amounts of phosphorus from the Baltic proper that are retained in this sub-basin. For the calculation of sub-basin budgets, the location of the lateral borders of the sub-basins is crucial, because net transports may change sign with the location of the border. Although the overall transport patterns resemble the results of previous studies, our calculated estimates differ in detail considerably.
The water exchange between the Baltic Sea and the North Sea is restricted by the narrows and sills in the Danish transition zone (Fig. 1). The hydrography of the Baltic Sea also depends on freshwater from rivers, which causes large salinity gradients between the surface layer and the saltier bottom layer and between the northern sub-basins and the entrance area (e.g., Meier and Kauker, 2003). The low-saline outflowing surface water is separated from high-saline inflowing bottom water by a transition layer, the halocline. The bottom water in the deep sub-basins is ventilated mainly by so-called Major Baltic Inflows (MBIs; Matthäus and Franck, 1992; Fischer and Matthäus, 1996). MBIs can significantly affect biogeochemical processes in the deep basins because of the inflow of large volumes of saline and oxygen-rich water into the Baltic Sea (e.g., Conley et al., 2009; Savchuk, 2010). In the Baltic Sea, the density stratification and long water residence time hamper the ventilation of deep waters. As a result, oxygen deficiency is a common feature. Additionally, nutrient loads from agriculture and other human activities of the large population in the catchment area increased nutrient concentrations in the water column. Actually, eutrophication has become a large environmental problem in the Baltic Sea in recent decades (e.g., Boesch et al., 2008; Pawlak et al., 2009; Wulff et al., 2001; Andersen et al., 2015). Therefore, accurate estimates of the ecological state and nutrient and water exchange between sub-basins and between the coastal zone and the open sea are of particular importance in managing the marine environment system.
The bathymetry of the model (depth in m). The border locations of sub-basins of the Baltic Sea used in this study are shown by the black lines, and the BY15 station is shown by the white star. The names of the sub-basins are Kattegat (KT), Danish Straits (DS), the Baltic proper (BP), the Gulf of Riga (GR), the Gulf of Finland (GF), the Bothnian Sea (BS), and the Bothnian Bay (BB).
On one hand, the estimation of biogeochemical processes, ecological state and nutrient exchange may rely on coupled marine ecosystem–circulation models (e.g., Neumann et al., 2002; Eilola et al., 2009, 2011; Almroth-Rosell et al., 2011, 2015; Maar et al., 2011; Daewel and Schrum, 2013). However, addressing biogeochemical cycles is a challenging task due to the complexity of the system. Obviously, there are large uncertainties in marine ecological simulations (e.g., Eilola et al., 2011). In contrast to the modeling of the physics of the atmosphere or ocean, where a basic description of the motion is provided by conservation equations, there is no basic set of equations that describe the marine ecosystem. Many biogeochemical processes are still poorly known and their uncertainties are difficult to quantify accurately. These potential sources of errors limit the applicability of the models both in forecasting and reanalysis. Further, imperfect initial conditions and model forcing also cause biases in the simulation results.
On the other hand, estimating nutrient budgets and transports between sub-basins may directly rely on observations and basin integrated budget models (Savchuk, 2005). The estimation accuracy depends on the spatial and temporal coverage of the measurements and the locations of borders between sub-basins. Although the data coverage in the Baltic Sea has gradually increased over time, the lack of observations still makes it difficult to estimate reliable biogeochemical cycles. Today, the availability of satellite sensor data like ocean color data from the OCTS (Ocean Color and Temperature Scanner) and from the SeaWiFS (Sea-viewing Wide Field-of-view Sensor) has provided the best spatial coverage of measurements. However, these sensors only give an estimate of a few biogeochemical parameters at the surface of the marine ecosystem and not the state of the entire marine ecosystem in the water column. Continuous observations of the deep ocean are only possible with in situ sensors, which have been deployed at only a limited number of stations (Claustre et al., 2010).
Given the coverage of observations and model deficiencies, we decided to perform a reanalysis based upon a high-resolution, coupled physical–biogeochemical model to estimate the physical and biogeochemical state of the Baltic Sea. For this purpose, data assimilation continuously updates the model variables at the locations of the observations and in their neighborhood. Integration in time of the prognostic model equations allows the spread of the information from the observations within the model domain.
The assimilation of data into coupled physical–biogeochemical models is confronted by various theoretical and practical challenges. For example, the response of the three-dimensional biogeochemical model to external forcing caused by the physical model is highly nonlinear. Further, it is difficult to use the biological observational information to reduce biases in the simulation of ocean physics which has an impact on modeled biogeochemistry (Béal et al., 2010). Besides, data assimilation as used in this study does not conserve mass, momentum and energy. Therefore, a reanalysis with data assimilation can never be dynamically fully consistent.
Nevertheless, the use of data assimilation complementing ecosystem modeling efforts has gained widespread attention (e.g., Hoteit et al., 2003, 2005; Allen et al., 2003; Natvik and Evensen, 2003; Triantafyllou et al., 2007, 2013; While et al., 2012; Teruzzi et al., 2014). Data assimilation into ecosystem models has focused both on parameter optimization and on state and flux estimations (Gregg et al., 2009). A comprehensive review of biological data assimilation experiments can be found in Gregg et al. (2009).
In the Baltic Sea, the biogeochemical data assimilation has started to become a research focus. For example, Liu et al. (2014) used the Ensemble Optimal Interpolation (EnOI) method to improve the multi-annual, high-resolution modeling of biogeochemical dynamics in the Baltic Sea. Fu (2016) analyzed the response of a coupled physical–biogeochemical model to the improved hydrodynamics in the Baltic Sea. Recently, several data assimilation studies have focused on the historical reanalysis of salinity and temperature in the Baltic Sea (e.g., Fu et al., 2012; Liu et al., 2013, 2014). Reanalysis has helped enormously in making the historical record of observed ocean parameters more homogeneous and useful for many purposes. For instance, ocean reanalysis data have been applied in research on ocean climate variability as well as on the variability of biogeochemistry and ecosystems (e.g., Bengtsson et al., 2004; Carton et al., 2005; Friedrichs et al., 2006). Ocean reanalysis can also be used for the validation of a wide range of model results (e.g., Fontana et al., 2013). For instance, the ocean mean state and circulation can be calculated from reanalysis results to evaluate regional climate ocean models (e.g., Meier et al., 2012). Moreover, reanalysis in the ocean is beneficial to the identification and correction of deficiencies in the observational records as well as for filling the gaps in observations. Regional and local model studies may use reanalysis results as initial and boundary conditions. A good reanalysis of biogeochemical state variables can dynamically describe indicators of eutrophication such as the long-term development of water nutrient pools.
The present paper focuses on the assimilation of profiles of temperature, salinity, nutrients and oxygen in the Baltic Sea following Liu et al. (2014). We aim to reproduce the ocean biogeochemical state with the help of information from both observations and a coupled physical–biogeochemical model for the period 1970–1999. Since 1970 the data coverage in the Baltic Sea is satisfactory. The results of the reanalysis are supposed to be used to estimate the water quality and ecological state with high spatial and temporal resolution in regions and during periods when no measurements are available. Further, nutrient transports across selected cross sections or between vertical layers are calculated with high resolution and accuracy, taking the complete dynamics of primitive equation models into account. This information cannot be obtained from either observations alone or from model results without data assimilation because the latter might have large biases in both space and time. We assess the nutrient budgets of the water column as well as of the nutrient exchanges between sub-basins and between the coastal zone and the open sea. The calculated budgets are compared to the results of other studies to evaluate our results. Hereby, we follow studies of other regions applying data assimilation for a biogeochemical reanalysis on long-term scales (Fontana et al., 2013; Ciavetta et al., 2016).
This paper is organized as follows. The physical and biogeochemical models are described in Sect. 2. Then the observational data set and the method of the reanalysis are introduced in Sects. 3 and 4, respectively. The experiment results, including comparisons with observations, are presented in Sect. 5. Finally, in Sects. 6 and 7, discussion and conclusions finalize the paper.
Number of observed profiles in different sub-basins (upper panel) and annual number of profiles from 1970–1999 (bottom panel).
The RCO (Rossby Centre Ocean) model is a Bryan–Cox–Semtner primitive
equation circulation model with a free surface (Killworth et al., 1991). Its
open boundary conditions are implemented in northern Kattegat, based on
prescribed sea level elevation at the lateral boundary (Stevens, 1991). An
Orlanski radiation condition (Orlanski, 1976) is used to address the case of
outflow, and the temperature and salinity variables are nudged toward
climatologically annual mean profiles to deal with inflows (Meier et al.,
2003). A Hibler-type dynamic–thermodynamic sea ice model (Hibler, 1979) with
elastic–viscous–plastic rheology (Hunke and Dukowicz, 1997) and a
two-equation turbulence closure scheme of the
The biogeochemical model called SCOBI (Swedish Coastal and Ocean
Biogeochemical model) has been developed to study the biogeochemical nutrient
cycling in the Baltic Sea (Marmefelt et al., 1999; Eilola et al., 2009;
Almroth-Rosell et al., 2011, 2015). This model handles biological and
ecological processes in the sea as well as sediment nutrient dynamics. SCOBI
is coupled to RCO (e.g., Eilola et al., 2012, 2013, 2014). With the help of a
simplified wave model, resuspension of organic matter is calculated from the
wave and current-induced shear stresses (Almroth-Rosell et al., 2011). SCOBI
has a constant carbon (C) to chlorophyll (Chl) ratio C : Chl
RCO–SCOBI is forced by atmospheric forcing data calculated from regionalized ERA–40 data using the regional Rossby Centre Atmosphere (RCA) model (Samuelsson et al., 2011). The horizontal resolution of RCA is 25 km. A bias correction method following Meier et al. (2011) is applied to the wind speed. Monthly mean river runoff observations (Bergström and Carlsson, 1994) are used for the hydrological forcing. Monthly nutrient loads are calculated from historical data (Savchuk et al., 2012).
The Baltic coastal shelf observation systems have been largely improved by
the joint efforts of the countries surrounding the Baltic Sea. For example,
the International Council for the Exploration of the Sea (ICES;
The simulated spatial variations of the late winter surface layer nutrient concentrations are compared with the spatial variations reconstructed from BED with the Data Assimilation System (DAS) by Sokolov et al. (1997). Due to insufficient historical data coverage, the average March fields were computed for time period 1995–2005 from over 3600 oceanographic stations found in BED. Also nutrient pools of dissolved inorganic nitrogen (DIN) and phosphorus (DIP) calculated with DAS (see Savchuk, 2010) are compared with the results of this study. See Eilola et al. (2011) for more details about the data handling by DAS.
Here we briefly describe the configuration of the data assimilation system
of this study. We focus on the state estimation via EnOI. The distribution
of stochastic errors are assumed to be Gaussian and non-biased. EnOI
estimates an “optimal” oceanic state at a given time using observations, the
numerical model and assumptions on their respective bias distribution
Based on the above configuration, two experiments from January 1970 to December 1999 have been carried out. One experiment is a simulation without data assimilation (FREE). The other simulation is constrained by observations using the “weakly coupled” assimilation scheme based upon the EnOI method following Liu et al. (2014) which was briefly described above (REANA). Both simulations, FREE and REANA, are initialized for January 1970. The initial conditions are taken from an earlier run with RCO–SCOBI. The observation error in REANA is defined according to Liu et al. (2014). However, in Liu et al. (2014), only a shorter assimilation experiment for a 10-year period is presented, and so far the reliability of the assimilation scheme in multi-decadal simulations has not been shown. Following Liu et al. (2014), our REANA experiment assimilated both physical and biogeochemical observations. In this study, we focus mainly on nutrient transports derived from the reanalysis.
To assess the results with (REANA) and without (FREE) data assimilation, the
overall monthly mean RMSDs (root mean square differences) of oxygen,
nitrate, phosphate and ammonium were calculated relative to observations
during the whole integration period. The overall monthly mean RMSD is
calculated by the following formula:
Monthly mean root mean square deviation (RMSD) between model results and observations for oxygen, nitrate, phosphate and ammonium in FREE (red) and REANA (blue).
Based on the reanalyzed simulation, the annual mean DIN and DIP transports
as well as DIP persistency are also calculated. These transports
(
The total nutrient budgets are calculated from the sum of inorganic and
organic bioavailable nutrients. The combined nutrient supplies from land and
from the atmosphere have been taken into account. Nitrogen fixation is not
included in the external supplies. The sinks of the nutrient budgets are
calculated from the supplies from land/atmosphere, import/export from other
basins and the changes in pelagic nutrient pools during the period
(sink
Annual mean integrated pools (in kt) of pelagic DIN and DIP in the Baltic proper calculated from REANA and from observations in BED.
In the following sub-sections, we evaluate the impact of data assimilation on the long-term evolution of biases (Sect. 5.1) and on vertical (Sect. 5.2) and horizontal (Sect. 5.3) distributions of nutrient concentrations. For the evaluation of time series of simulated oxygen, nitrate, phosphate and ammonium concentrations, the reader is referred to Liu et al. (2014; Figs. 6 and 7). After the evaluation of the assimilation method, we focus on the analysis of nutrient transports in the Baltic Sea based upon our reanalysis data that we consider to be the best available data set for such an analysis. In particular, we analyze the horizontal circulation of nutrients (Sect. 5.4), the horizontal distribution of nutrient sources and sinks, the nutrient exchange between the coastal zone and the open sea (Sect. 5.5), and the nutrient budgets of sub-basins (Sect. 5.6).
The seasonal cycle of monthly average (1970–1999) temperature
(
The biases in both FREE and REANA have been calculated relative to
observations for dissolved oxygen and inorganic nutrients using Eq. (4)
(Fig. 3). Data assimilation has significantly reduced bias of the model
simulation. Generally, the RMSDs of oxygen and nutrient concentrations in
REANA are smaller than those of FREE. However, the degree of improvements
differs among the variables. The RMSD of oxygen is mostly greater and smaller
than 1.0 mL L
Further, the annual averaged pelagic pools of simulated DIN and DIP in the Baltic proper are compared to the corresponding pools estimated from the BED data in the Baltic proper (Fig. 4). The maximum annual differences of DIN and DIP pools compared to BED have been reduced by 57.5 and 72.3 %, respectively, from 400 and 650 kt in FREE (not shown) to 170 and 180 kt in REANA. The remaining differences between REANA and BED may be explained by the methods of integration that differ. BED estimates are based on a limited amount of observations, while the model results are based on a large number of grid points with dynamically varying state variables. A similar result for hypoxic area was found by Väli et al. (2013). They showed with the help of model results which were sampled at the same times and locations of the observations that the applied interpolation algorithm underestimated the hypoxic area by about 40 %.
The long-term average seasonal cycles of temperature and inorganic nutrients at monitoring station BY15 at Gotland deep (for the location, see Fig. 1) give a hint of how data assimilation improves simulated nutrient dynamics in the Baltic proper (Fig. 5). The surface layer temperature and stratification show rapid increase from April to May, with concurrent rapid decrease of nutrient concentrations due to primary production down to 50–60 m depths. The cooling and increased vertical mixing in autumn and winter reduce temperatures and bring nutrients from the deeper layers into the surface layers. RCO–SCOBI captures these variations. However, compared to BED, FREE has obvious biases, such as overestimated temperature stratification around a 30–50 m depth from late winter to early spring, higher concentration of nutrients at the 50–60 m depth, stronger vertical stratification of nutrient concentrations and less decrease of nutrients in the summer, especially below the thermocline, as well as also in the surface layers for phosphate. One reason for the biases is the vertical displacement of the halocline that is too shallow in RCO (e.g., Fig. 4 in Liu et al., 2014). The causes for the model bias in nutrient depletion below the summer thermocline are not known, but possible reasons are discussed by Eilola et al. (2011). These biases are significantly reduced in the reanalysis, which provides an improved description of vertical distributions of nutrients in the layers above the halocline.
The average March concentrations (1970–1999) of DIP and DIN in the upper layers (0–10 m), as well as their ratio (DIN : DIP), were calculated (Fig. 6). Due to insufficient historical data the corresponding BED maps describe the averages for the period 1995–2005 as a basis for model-to-data comparison. In BED results the highest concentration of DIP occurs in the Gulf of Riga and the Gulf of Finland. Relatively high concentrations of DIP are found in the entire Gotland Basin. The DIP concentrations in the Bothnian Sea and Bothnian Bay are obviously lower than in other regions. Generally, the DIP in FREE has been largely overestimated in all regions relative to BED, especially in the Baltic proper. Further, in BED, high concentrations of DIN occur in coastal waters close to the river mouths of the major rivers in the southern Baltic proper. DIN concentrations in the Gulf of Finland and in the Gulf of Riga are also high and cover large areas of these gulfs. Unlike the BED data, the DIN in FREE also has high concentrations in the entire southern and eastern coastal zones of the Baltic proper. As a result, FREE shows a gradient in DIN concentrations between the coastal zone and the open sea in the entire southern Baltic proper. The DIN and DIP patterns result in high and low DIN : DIP ratios in the Bothnian Bay and Baltic proper, respectively. The highest DIN : DIP ratios are found in the Bothnian Bay in BED and in the Gulf of Riga in FREE. FREE has captured this large-scale pattern, but there are substantial regional differences. By the constraints of the observation information, REANA has improved the spatial distributions of DIN and DIP significantly compared to BED. In particular, DIP concentrations in REANA are much closer to the interpolated BED observations.
Simulated monthly (March) mean (1970–1999) surface layer (0–10 m) concentrations of DIP (mmol m
Nutrient transport directly affects the biogeochemical cycles and the eutrophication of the Baltic Sea. The persistency of vertically integrated transports (Fig. 7) is defined, for instance, by Eilola et al. (2012). One should note that the results by Eilola et al. (2012) are based upon 30-year averages for the control period 1978–2007 of a downscaled climate scenario from a global circulation model. Similar calculations of transports will therefore be briefly presented in the present study, since the hindcast period is better represented when the model is forced by the assimilated atmospheric (ERA–40) and Baltic Sea data (REANA). DIP has the largest transports in the central parts of the Baltic proper, with high persistency because the volume transports are generally larger in deeper rather than in shallower areas. In the Bornholm Basin and the eastern parts of the central Baltic proper, cyclonic circulation patterns are found. In the western parts of the central Baltic proper, southward transports prevail. Relatively large magnitudes of transports of DIP are also found in the northwestern Gotland Basin, in the southern Bornholm Basin, and through the Slupsk Channel connecting the Bornholm Basin and Gotland Basin. Similar transport patterns are also found for DIN, OrgP and OrgN (not shown). Compared to Eilola et al. (2012), DIN, DIP, OrgP and OrgN transports and their persistency are larger, although the overall patterns are similar. For example, in Eilola et al. (2012, their Fig. 1), large DIN transports appear in the southern Baltic proper and the Bornholm Basin. Similar differences are also found in both DIP and OrgP transports.
Annual mean DIP transports and the corresponding DIP persistency,
DIN and OrgP transports for REANA averaged for the period 1970–1999. The
black solid lines with arrows show the streamlines and direction of
transports. The magnitude of transports (kt km
Spatial distributions of annual mean import of DIP, OrgP, DIN and
OrgN averaged for the period 1970–1999. The magnitude of import and its
corresponding value (kt km
The horizontal distributions of areas with internal sources and sinks of
phosphorus and nitrogen are illustrated in Fig. 8. A net inflow (inflow
Annual mean, accumulated net imports (black lines), and imports of OrgP, DIP and DIN (color bars) to regions with the same depth in the Baltic proper averaged for the period 1970–1999.
According to the accumulated import of nutrients (Fig. 9), the magnitude of the DIP export is larger than that of the DIP import. This indicates that not all of the supply of phosphorus from land and atmosphere is retained within the Baltic proper as will be further discussed from the nutrient budgets in Sect. 5.6. For DIN, however, we may notice only a very small net export from the Baltic proper to adjacent sub-basins, while for OrgP and OrgN, imports and exports are almost balanced (Fig. 9). The nitrogen and phosphorus supply from land is implemented in sea areas with a bottom depth usually of 6 m. This is where the river mouths are located in the model.
There is a large import of DIP to areas with a depth range between 40–70 m (Fig. 9). This import does not show a counterpart in the export of OrgP in Fig. 9. This result might be explained by local processes causing the phytoplankton uptake and sediment deposition of DIP. There is an import of DIN to these areas that together with nitrogen fixation and sediment–water fluxes of DIN may support local production of organic matter. The phosphorus sink may be partly caused by oxygen-dependent water–sediment fluxes that bind DIP to ironbound phosphorus in oxic sediments (Almroth-Rosell et al., 2015). This effect is not included in Eilola et al. (2012) but might potentially be accounted for by the adjusted DIP transports in REANA. The results of REANA indicate that there is an additional sink, but the relative importance of the different processes causing this sink (data assimilation or sediment processes) is, however, not possible to evaluate from the present reanalysis data set.
A partly opposite exchange profile is found for OrgP (Fig. 9). Coastal areas with a water depth of up to 40 m are exporting organic phosphorus, whereas deeper areas import OrgP. Production in the coastal zone of the Baltic proper and sedimentation in the open sea is almost balanced.
The largest export of DIN occurs due to rivers in the very shallow coastal zone. The magnitude of DIN imports and exports in areas with greater water depths are much smaller. Obviously, DIN supplied from land is already consumed in the coastal zone (Voss et al., 2005; Almroth-Rosell et al., 2011), and, consequently, only a minor fraction of the nitrogen supplied to the shallow area can continuously reach regions deeper than 100 m (Eilola et al., 2012; Radtke et al., 2012).
The Baltic Sea is divided into seven sub-basins according to the selected
sections, which form the borders of the sub-basins (Fig. 1). We calculate
total nutrient budgets for each of the sub-basins from the reanalysis
(Figs. 10 and 11). Changes in pools are calculated as differences between
1971 and 1999 because the initial adjustment process due to the assimilation
is taking place during the first year (1970; not shown). The largest annual
mean external phosphorus load occurs in the Baltic proper and amounts to
34.2 kt yr
Annual mean total phosphorus budgets of the Baltic Sea averaged for
the period 1971–1999. The average total amounts are in kt, and transport
flows and sink/source fluxes (external nutrient inputs/sink) are in
kt yr
Nitrogen transports between Baltic Sea sub-basins are different compared to phosphorus transports (Fig. 11). For example, the Baltic proper has larger nitrogen sinks than external sources. Further, during the period 1971–1999 the nitrogen content decreased in the Gulf of Riga and increased in the Bothnian Bay, respectively. In the Gulf of Finland and Danish Straits, the difference between the external supply and internal sink of nitrogen is equal to the net transport into the Gulf of Finland and Danish Straits. The large sink of nitrogen in the Bothnian Bay is noteworthy. We also found relatively large net transports of nitrogen from the Gulf of Riga into the Baltic proper. This is mainly explained by the relatively high nitrate concentrations in the Gulf of Riga relative to other sub-basins.
The same as Fig. 9, but for nitrogen.
Annual mean fluxes of nitrogen (in kt yr
The 3-D nutrient pools constructed by data assimilation methods offer an opportunity to evaluate with improved estimates the changes in Baltic Sea eutrophication. As an example, an investigation of the trophic state from changes in 5-year average nutrient pools in REANA shows that the total nitrogen pool in the Baltic proper increased from 657 to 1045 kt from the period 1971–1975 to 1995–1999, while the total phosphorus pool decreased from 469 kt to 448 kt between the same periods. Hence, nitrogen increased by about 59 %, while phosphorus decreased by about 4 %. Similarly the pool of DIN in REANA increased by 80 %, while DIP decreased by 6 %. The corresponding numbers obtained from BED showed an increase of 100 % for DIN and an increase of 7 % for DIP. The results indicate large increases of nitrogen pools in the Baltic proper during the investigated period but only relatively small changes of phosphorus pools.
Subsequent periodic assessments can be used to reveal future eutrophication
changes. While estimating the trophic state, it should be noted that the
change in trophic state depends on the chosen time periods. For example, from
the year 1971–1999 the total phosphorus (TP) in the Baltic proper increased
by 2.7 kt yr
To further analyze the variability of the budget of the reanalyzed nutrients,
Fig. 12 provides the cross sectional, integrated nutrient flows in the
different sub-basins. Here the eastward and northward net transports are, by
definition, positive. Obviously, the integrated nutrient flows vary
significantly in space according to the nutrient loads from land. The inflows
and outflows also vary depending on the depth of the water column and
nutrient concentrations that influence the vertically integrated mass fluxes.
In general, the magnitude of nutrient transports declines along transect A
from south to north. For instance, the largest annual northward flow of
nitrogen in the Baltic proper reaches 392 kt yr
In the Baltic proper, inflow and outflow as well as the net northward flow of
phosphorus increase from the south until a section along 56.8
In the Arkona and Bornholm basins, nitrogen and phosphorus transports increase from the west to the east. Due to the nitrogen load from the Oder River, the inflow of nitrogen increases significantly at the border between the Arkona and Bornholm basins, whereas the outflow does not show any discontinuity. As a result, the net flow of nitrogen shows an accelerated increase. The situation for phosphorus in the Arkona and Bornholm basins is different compared to the nitrogen transports because in- and outflow change direction. The phosphorus loads from the Oder River turn the outflow in the western parts into an inflow of phosphorus in the eastern parts.
In the Gulf of Finland, in– and outflows generally decline from the west to
east. In the entrance of the Gulf of Finland, the net inflows of nutrients
are almost zero. The largest net flow (westward) of nutrients appear at the
inner end of the Gulf of Finland, where the large river Neva enters the Gulf,
with a magnitude of 33 kt yr
RCO–SCOBI has been widely used for the Baltic Sea and the model was carefully evaluated using various observational data sets. As any other model, RCO–SCOBI had to be calibrated because many processes including sources and sinks of nutrients are not known in enough detail. Hence, an optimal parameterization of unresolved processes is one of the requirements for the predictive capacity of the model. Further requirements to calculate correct transports and transformation processes in addition to optimized model equations are high-quality atmospheric and riverine forcing data and high-quality initial and lateral boundary conditions.
Most of the large biases in FREE are caused by imperfect initial conditions. The reason is that the nutrient pools in the sediments have not been spun-up appropriately. As a consequence, phosphate concentrations in FREE are higher than observed concentrations at all depths. The biases in surface phosphate concentrations between model results and observations can influence the seasonal primary production. In REANA, however, from the beginning of the experiment, the biases are already significantly reduced during the first year and remain relatively small during the integration compared to FREE. This result indicates a need of new initial conditions of the sediments.
In the long-term simulation, the new initial condition for an assimilation cycle differs from the ending ocean state of the last cycle when at that time observations are available. In this sense, the data assimilation introduces sources and sinks of the nutrient cycles by interrupting the model simulation and adjusting the initial conditions. The magnitudes of these artificial sources and sinks are directly related to the biases between model results and observations. Figure 3 shows that the model has large biases during the beginning of the simulation. However, data assimilation has corrected the mismatch between model state and observation to an optimal level during an initial adjustment period. After the adjustment period, the mismatch between model and observation becomes small and the successive adjustment due to data assimilation also becomes small. Further, the adjustment of data assimilation is related to the spatial–temporal coverage of observations. Here we assimilated only observed profiles into the model. After every assimilation cycle, the simulation continues with optimal initial conditions based upon conservation principles. As the equations of RCO–SCOBI have not been changed, masses of all constituents of the model are conserved at least during the simulation between two assimilation occasions.
The advantage of the data assimilation is that model variables at any station are very likely more accurate than the model output without data assimilation. For instance, time series of profiles or transports across vertical sections very likely have a smaller bias compared to observations than the corresponding model results without data assimilation. Compared to available observations the information from the model is higher-resolved and homogeneous in space and time. Of course, it is difficult to evaluate the quality of model results at high resolution because independent observational data sets are usually missing. An exceptional effort to utilize independent data was done by Liu et al. (2014), showing that the statement about the added value of data assimilation is true for one available, independent cruise data set at high resolution.
The results of the reanalysis can be used to estimate the water quality and ecological state with high spatial and temporal resolution in regions and during periods when no measurements are available. This supports improved assessments, for example, of eutrophication status indicators as exemplified in Sect. 5.6. Regional and local model studies may use the data as initial and boundary conditions. For projections of future climate and for nutrient load abatement scenario simulations the reanalysis has a very high scientific value as a reference data set for the historical period of the climate simulations. The evaluation of the regionalized climate (the statistics of mesoscale variability, e.g., the mean state) during the historical period can be done much more accurately based upon the reanalysis data than with sparse observational data. For instance, it is very difficult to calculate the climatological mean state just from observations that are cast only during the ice-free season of the year. Using a reanalysis as reference data for historical climate is a common method in regional climate studies of the atmosphere. Here we provide a corresponding data set for the ocean to evaluate simulated present-day climate.
Further, nutrient transports across selected cross sections or between vertical layers are calculated from the reanalysis with high resolution and improved accuracy. However, one cannot expect that budgets calculated from the summation of internal fluxes from model results with data assimilation are more accurate, because usually small artificial sources and sinks from the data assimilation are becoming as important as physically and biogeochemically motivated sources and sinks when sums of fluxes are compared. Hence, in Sect. 5.6 we calculated budgets only from external supply, imports and exports and changes in the water pools of nutrients with the aim of comparing the reanalysis results with other studies using only observations. It is perhaps not possible to claim that our budgets are more accurate than budgets that are derived from observations only, despite the higher temporal and spatial resolution in model outputs. However, the advantage of the reanalysis is that measurements are extrapolated in space and time based upon physical principles of the model.
Fu (2013) estimated the volume and salt transports during the 2003 MBI with a three-dimensional variational data assimilation method (3DVAR) in the Baltic Sea. In the present study, we estimate the impact of the data assimilation based upon the EnOI method, on the net volume and nutrient transports as well as calculate budgets for major sub-basins of the Baltic Sea. The volume transports obtained with different assimilation methods may be different. The sea level in Fu (2013) is kept constant in the assimilation process, while sea level in this study is varying accordingly during the assimilation of temperature and salinity based upon the statistical covariances. The variability of sea level may enhance the barotropic flow, which is one of the reasons for the differences in net volume transport between the results by Fu (2013) and REANA. However, transports within the sub-basin are also indirectly affected by the interaction of baroclinicity and topography.
In contrast to Eilola et al. (2012), in this study areas with DIN export are also found at the southern and eastern coasts as well as at some small local regions in the inner parts of the Baltic proper (Fig. 8). In REANA, the magnitudes of DIP imports and exports are larger than in Eilola at al. (2012), and there is pronounced import of DIP in the western part of the eastern Gotland Basin to regions deeper than 100 m (Fig. 8) that is not as significant in Eilola et al. (2012). This, and the larger variability of DIN imports and exports, indicates that there is a higher degree of small-scale localized transport and production patterns that are not captured by Eilola et al. (2012). Main sinks of DIN are found in the deeper areas, but significant sinks are also seen in shallow areas and areas with water depths of about 60 m. As the assimilation of salinity observations results in a deeper halocline (Liu et al., 2014), the bottom water in a depth range of 40–70 m contains higher oxygen concentrations than in the simulation without data assimilation. Hence, in the REANA simulation of this study, more phosphorus is retained in regions with water depth in the range of 40–70 m than in the simulation by Eilola et al. (2012). The present results show, however, an export contribution from DIN sources in deeper areas (e.g., 60–90 m depths) that may have been caused by reduced denitrification efficiency of oxidized sediments in the REANA simulation compared to Eilola et al. (2012).
The in- and outflows of phosphorus between the sub-basins, except the Gulf of Riga and Gulf of Finland, simulated in REANA are smaller than the results by Wulff and Stigebrandt (1989), Savchuk (2005), and Savchuk and Wulff (2007). However, the net transports of phosphorus are similar between our results and these earlier studies in sub-basins except Kattegat and Danish Straits. Moreover, the nitrogen budgets are much lower than the results of earlier studies, especially in the Baltic proper. The earlier studies included transports of organic nutrients assumed not bioavailable, comparably high especially for organic nitrogen, in the Baltic Sea, while the present budgets include only biologically active nutrients. This may explain much of the differences in nitrogen transports. It should be kept in mind that the above-mentioned studies estimated the nutrient budgets from mass balance models together with inter-basin transport calculations based upon Knudsen's formula to calculate nutrient budgets of the Baltic Sea (see, e.g., Savchuk, 2005). Obviously, there are limitations in calculations of previous studies. Despite overall uncertainties that also limit the reliability of our results, like incomplete understanding of selected biogeochemical processes (e.g., nitrogen fixation), lacking information of sediment parameters and under-sampled observations in space and time, our approach has the advantage of using both high-resolution modeling and all available observations made over a 30-year period. Our model results consider the complete set of primitive equations in high resolution, taking into account not only the volume and salt conservation of sub-basins according to Knudsen's formula, but also the wind-driven circulation between and within sub-basins. Hence, we have, for the first time, the potential to quantify spatial transport patterns with high confidence even within sub-basins, as in the exchange of nutrients between the coastal zone and the open sea.
Eutrophication of the Baltic Sea is directly affected by the long-term evolution of external nutrient supply that has three components (waterborne land loads, direct point sources at the coasts and atmospheric depositions) which are associated with the biogeochemical dynamics of the Baltic Sea. In our study, we used the reconstructed external nutrient input data by Savchuk et al. (2012). Nutrient budgets (Figs. 10 and 11) of sub-basins are time-averaged and represent in our study the overall results of the period 1971–1999. The phosphorus loads vary in different periods, for example, the phosphorus loads in the 1980s are larger relative to the 1990s (see Savchuk et al., 2012). Therefore, the phosphorus supply into the Gulf of Finland is greater in our study compared to Savchuk and Wulff (2007). The greater phosphorus supply changes the phosphorus content and phosphorus concentration in the Gulf of Finland. This is one reason why phosphorus transports between the Gulf of Finland and the Baltic proper in our study are greater than the transports calculated by Savchuk (2005) and Savchuk and Wulff (2007).
Since our study covers a different time period compared to the studies by Wulff and Stigebrandt (1989), Savchuk (2005), and Savchuk and Wulff (2007), nutrient concentrations and related budgets differ in time and space. Hence, it is not surprising that other studies show deviating results. For example, during the period 1970–1999, HELCOM (2013) showed that the TP concentration generally decreased in the Bothnian Bay and has increased in the Gulf of Riga. However, these changes in TP concentrations were not monotonous. For example, the TP concentration obviously increased during the period 1970–1976 in the Bothnian Bay. While in the Bothnian Sea, TP concentration increased during the period 1970–1983 and decreased during the period 1990–1999. Similarly, changes in total nitrogen concentration differed during different periods.
Gustafsson et al. (2012) used a process-oriented model that resolves the
Baltic Sea spatially in 13 dynamically interconnected and horizontally
integrated sub-basins with high vertical resolution to reconstruct the
temporal evolution of eutrophication for 1850–2006. Savchuk (2005) and
Savchuk and Wulff (2007) applied mass balance models as mentioned above to
calculate nutrient budgets of the Baltic Sea. The results of all these models
depend on the locations of the sub-basin borders which are chosen as far as
possible according to dynamical constraints such as sills or fronts that are
parameterized to obtain estimates of the water exchanges. Using a
high-resolution circulation model, we showed that nutrient flows within the
sub-basins may vary considerably (Fig. 12). For instance, we found east- and
westward net transports of nitrogen between the Baltic proper and Gulf of
Finland depending on border locations at 23.2 and 24.0
For the first time, a multi-decadal, high-resolution reanalysis of physical (temperature and salinity) and biogeochemical variables (oxygen, nitrate, phosphate and ammonium) for the Baltic Sea was presented. The reanalysis covers the period 1970–1999. A weakly coupled assimilation scheme using the EnOI method was used to assimilate all available physical and biogeochemical observations into a high-resolution circulation model of the Baltic Sea.
Both assimilated and independent observations collected from different
databases were used to evaluate the reanalysis results (REANA). Based on the
model–data comparison presented in this study, we found that the model
results without data assimilation (FREE) exhibit significant biases in both
oxygen and nutrients. The reasons for these biases are not totally understood
yet, although it is speculated that the main reasons might be related to the
imperfect initial conditions, limitations of model parameterizations, the
inaccurate halocline position and, correspondingly, the hypoxic volume (Liu et
al., 2014). Based on the calculation of the overall RMSD of oxygen and
nutrient concentrations between model results and not-yet-assimilated
observations, the results in REANA are considerably better than those in
FREE. The total RMSD of the oxygen, nitrate, phosphate and ammonium is
reduced respectively by 0.84 mL L
The observation information entering the model affects the oxygen-dependent dynamics of biogeochemical transports significantly due to improved simulation of both physical (e.g., vertical stratification) and biogeochemical parameters (e.g., nutrient concentrations). As examples, we presented improved results of mean seasonal cycles of nutrients and of spatial surface distributions of DIN, DIP and DIN : DIP of the entire Baltic Sea.
Based on the reanalysis simulation, we analyzed nutrient transports in the Baltic Sea. We found that vertically integrated nutrient transports follow the general horizontal water circulation and vary spatially to a large extent. In particular, large nutrient transports were found in the eastern Gotland Basin, in the Bornholm Basin, in the Slupsk Channel and in the northwestern Gotland Basin. The persistence of nutrient transports is greater in the eastern and southern than in the northern and western Baltic Sea.
The horizontal distributions of sources and sinks of inorganic and organic nutrients show large spatial variations and may be partly explained by (1) the external supply of nutrients from land, (2) the topographically controlled horizontal nutrient exchange between sub-basins and between the coastal zone and the open sea and (3) vertical stratification that determines redox conditions at the sea floor. The latter is important for the sediment–water fluxes of nutrients and, consequently, for burial of nutrients in the sediments. The reanalysis results suggest that in the Baltic proper, in most areas with a water depth less than the depth of the permanent halocline at about 70–80 m, DIP is imported and transformed either to OrgP or buried in the sediments in water depths greater than the wave-induced zone at 40–70 m. Whether the latter is an artifact of the assimilation method or a real sink is unclear. On the other hand, in areas with greater water depth, DIP is exported (e.g., released from the sediments under anoxic conditions). Overall, the Baltic proper exports DIP to neighboring sub-basins.
Nitrogen transports are very different compared to phosphorus transports. The shallow coastal zone with water depths less than 10 m plays an outstanding role for DIN, because within it, large exports occur due to supplies from land. The high productivity in the shallow areas effectively transfers DIN to OrgN, and denitrification decreases the exports of nitrogen from coastal areas to the deeper areas. Most of the exported DIN is removed in shallow waters, while at greater depths imports and exports of DIN are much smaller, indicating the important role of the coastal zone for nitrogen removal.
Detailed nitrogen and phosphorus budgets suggest that nutrient transports in
the various sub-basins are controlled by different processes and show
different responses to external loads and internal sources and sinks. In
particular, the Baltic proper is the sub-basin with the largest nutrient
exchanges with its surrounding sub-basins. The Baltic proper exports
phosphorus to all sub-basins except the Gulf of Riga. Similarly, the Baltic
proper also exports nitrogen to all sub-basins except to the Gulf of Riga and
Danish Straits. In this sub-basin, the largest internal sink of all
sub-basins was also found. The relatively large net export of
phosphorus from the Baltic proper into the Bothnian Sea is noteworthy. This finding is in
agreement with previous studies. For the budgets of the sub-basins,
the location of the lateral borders of the sub-basins
is important, because net
transports may change sign with the location of the border. For instance, in
the entrance of the Gulf of Finland, the net phosphorus transport from the
Baltic proper is directed eastward but changes direction at about
26
The data sets used in this study can be obtained from the corresponding author (ye.liu@smhi.se).
The authors declare that they have no conflict of interest.
The research presented in this study is part of the Baltic Earth programme
(Earth System Science for the Baltic Sea region, see