Modeling the combined impact of changing climate and changing nutrient loads on the Baltic Sea environment in an ensemble of transient simulations for 1961–2099

The combined future impacts of climate change and industrial and agricultural practices in the Baltic Sea catchment on the Baltic Sea ecosystem were assessed. For this purpose 16 transient simulations for 1961–2099 using a coupled physical-biogeochemical model of the Baltic Sea were performed. Four climate scenarios were combined with four nutrient load scenarios ranging from a pessimistic business-as-usual to a more optimistic case following the Baltic Sea Action Plan (BSAP). Annual and seasonal mean changes of climate parameters and ecological quality indicators describing the environmental status of the Baltic Sea like bottom oxygen, nutrient and phytoplankton concentrations and Secchi depths were studied. Assuming present-day nutrient concentrations in the rivers, nutrient loads from land increase during the twenty first century in all investigated scenario simulations due to increased volume flows caused by increased net precipitation in the Baltic catchment area. In addition, remineralization rates increase due to increased water temperatures causing enhanced nutrient flows from the sediments. Cause-and-effect studies suggest that both processes may play an important role for the biogeochemistry of eutrophicated seas in future climate partly counteracting nutrient load reduction efforts like the BSAP.


Introduction
For the Baltic Sea ( Fig. 1) regional climate modeling results suggest that global warming may cause increased water temperatures and reduced sea ice cover combined (eventually) with reduced salinity due to increased wind speeds and increased river runoff (e.g. BACC author team 2008; Meier et al. 2006). The projected hydrographic changes could therefore have significant impacts on the marine ecosystem. To estimate these effects and to calculate the impact of nutrient load reductions in future climate an ensemble of model simulations for the period 1961-2099 were carried out. Ensemble simulations are necessary to estimate uncertainties of projections (e.g. Nikulin et al. 2011;Räisänen et al. 2004). Uncertainties are caused by biases of global and regional climate models and by unknown socio-economic future developments with impact on greenhouse gas emissions, nutrient loads from land and atmospheric deposition.
For the marine environment of regional seas only a few studies on uncertainties of future projections are available (e.g. Neumann 2010; Meier et al. 2011b). For instance, Neumann (2010) studied the results of two transient simulations with a coupled physical-biogeochemical model driven by regionalizations of one General Circulation Model (GCM) forced with two greenhouse gas emission scenarios [A1B and B1, see Nakićenović et al. (2000)]. He found that at the end of the century the oxygen conditions in the deep water of the Baltic Sea will sligthly improve. Due to increasing water temperatures the spring bloom in the northern Baltic Sea will start earlier and the cyanobacteria season will be prolonged. However, both the total phytoplankton biomass and nitrogen fixation do not change significantly. Meier et al. (2011b) performed 16 scenario simulations based upon the delta approach assuming that only the mean seasonal cycle of the atmospheric and hydrological cycles will change. They studied four projections of future climate (two GCMs and two greenhouse gas emission scenarios, A2 and B2) combined with four socio-economic scenarios affecting nutrient loads. Meier et al. (2011b) concluded that the uncertainties caused by the driving GCMs are considerable. Depending on the chosen driving GCM, wind induced mixing and, consequently, bottom oxygen concentrations will change even if river nutrient concentrations are assumed to be unchanged. In all studied scenario simulations phytoplankton concentrations in the southwestern Baltic increase in future climate. However, the model response in the northern Gotland Basin differs considerably.
In this study we further investigate uncertainties in projections of the Baltic environment from the statistics of an ensemble of 16 transient simulations including also differing nutrient load scenarios. Regionalized data by Meier et al. (2011c) from four scenario simulations driven by two GCMs and two greenhouse gas emission scenarios [A1B and A2, see Nakićenović et al. (2000)] are used to force a state-of-the-art coupled physical-biogeochemical model of the Baltic Sea Eilola et al. 2009). These four climate scenarios are combined with four nutrient load scenarios suggested by HELCOM (2007): a reference scenario assuming present-day nutrient concentrations in the rivers, a pessimistic business-as-usual scenario assuming an exponential growth in agriculture in all Baltic Sea countries, a scenario of riverine nutrient loads and atmospheric deposition according to current legislations and a more optimistic case following the Baltic Sea Action Plan (BSAP) .
In this study the approach by Meier et al. (2011b) is further refined by taking new, high-resolution model versions and time-dependent (transient) scenario simulations for the twenty first century into account. No assumptions on changes of the variability and no restriction on selected time slices in present and future climates (like in the delta approach) were made. Compared to Meier et al. (2011b) the horizontal and vertical resolutions are considerably increased.
The paper is organized as follows: In the next section the method of the dynamical downscaling approach and the involved models are briefly introduced. In the third section results of annual or seasonal mean changes of selected ecological quality indicators are presented and discussed. The results of cause-and-effect studies identifiy the dominating drivers of the simulated changes. Finally, some conclusions of the study are highlighted.   . It is coupled to a Hibler-type sea ice model with elastic-viscous-plastic rheology (Hunke and Dukowicz 1997). Subgrid-scale vertical mixing is parameterized using a turbulence closure scheme of the k-e type (Rodi 1980) with flux boundary conditions to include the effect of a turbulence enhanced layer due to breaking surface gravity waves and a parameterization for breaking internal waves (Meier 2001). In the present study, RCO was used with a horizontal resolution of 3.7 km (2 nautical miles) and with 83 vertical levels with layer thicknesses of 3 m. A fluxcorrected, monotonicity preserving transport (FCT) scheme is embedded (Gerdes et al. 1991) and no explicit horizontal diffusion is applied. For further details of the RCO model the reader is refered to Meier (2001Meier ( , 2007, .

Biogeochemical model
The Swedish Coastal and Ocean Biogeochemical model [SCOBI, e.g. Eilola et al. (2009)] is coupled to the physical model RCO. SCOBI describes the dynamics of nitrate, ammonium, phosphate, phytoplankton, zooplankton, detritus, and oxygen. Here, phytoplankton consists of three algal groups representing diatoms, flagellates and others, and cyanobacteria (corresponding to large, small and nitrogen fixing cells). Besides the possibility to assimilate inorganic nutrients the modelled cyanobacteria also has the ability to fix molecular nitrogen which may constitute an external nitrogen source for the model system. The sediment contains nutrients in the form of benthic nitrogen and benthic phosphorus including aggregated process descriptions for oxygen dependent nutrient regeneration, denitrification and adsorption of ammonium to sediment particles, as well as permanent burial of organic matter. With the help of a simplified wave model the combined effect of waves and current induced shear stress is considered to calculate resuspension of organic matter (Almroth-Rosell et al. 2011). For further details of the SCOBI model the reader is refered to Marmefelt et al. (1999), Eilola et al. (2009), Almroth-Rosell et al. (2011 2.3 Regional climate data sets Four climate change scenario simulations have been performed. The forcing was calculated applying a dynamical downscaling approach using the regional climate model RCAO [Rossby Centre Atmosphere Ocean model, see Döscher et al. (2002)] with lateral boundary data from two GCMs. The two GCMs used were HadCM3 from the Hadley Centre in the UK (Gordon et al. 2000) and ECHAM5/MPI-OM from the Max Planck Institute for Meteorology in Germany (Roeckner et al. 2006;Jungclaus et al. 2006), henceforth short ECHAM5. HadCM3 and ECHAM5 simulations were forced with one (A1B) or two (A1B and A2) greenhouse gas emission scenarios, respectively. In addition, for our scenario simulations two realizations of ECHAM5 forced with the emission scenario A1B, denoted with -r1 and -r3 (ECHAM5-r1-A1B and ECHAM5-r3-A1B), with differing initial conditions in the year 2000 were used. [The Max Planck Institute performed also a third realization, ECHAM5-r2-A1B, which is not used here, see ]. Thus, the atmospheric forcing for RCO-SCOBI is calculated from RCAO-HadCM3-A1B, RCAO-ECHAM5-r3-A1B, RCAO-ECHAM5-r1-A1B and RCAO-ECHAM5-r1-A2.
Future projections refer to a period at the end of this century . Annual and seasonal mean changes were calculated from the differences between the periods 2070-2099 and 1978-2007. For further details of the downscaling method and the quality of the atmospheric forcing fields the reader is refered to Meier et al. (2011c). In contrast to earlier studies by Meier (2006) and Meier et al. (2011b) no bias correction of the atmospheric forcing was applied. An exception is the wind speed in 10 m height. Following Höglund et al. (2009), the wind speed is modified using simulated gustiness to improve wind speed extremes (Meier et al. 2011c).

SSH at the open boundaries
Following Gustafsson and Andersson (2001), sea surface height (SSH) in Kattegat is estimated from the daily averaged meridional atmospheric pressure gradient difference DP between two grid points located in The Netherlands and Norway. Thus, DP n and DP nþ1 are defined as the meridional pressure gradients at day n and day n ? 1, respectively. The SSH g at day n is calculated from The coefficients a and b are computed using a simple optimisation method in order to get the best possible fit to sea level observations in Smögen located at the Swedish west coast close to the open boundary of the model domain in Kattegat. For the optimisation procedure atmospheric pressure data from the Rossby Centre Atmosphere model [RCA, Samuelsson et al. (2011)] driven with ERA40 re-analysis data (Uppala et al. 1989) at the lateral boundaries are used. This approach provides a good correlation of calculated and observed SSHs, but the calculated standard deviations are too small compared to observations. The probability density function reveals that positive extremes of SSH are underestimated (not shown). These extremes are essential for salt water inflows into the Baltic Sea.
If the calculated SSH is used as forcing for the Baltic Sea model, the overall salinity of the Baltic Sea will Transient scenario simulations for the Baltic Sea 2423 decrease unrealistically. We suspect that this shortcoming of the estimated SSH is related to underestimated atmospheric depressions in RCA causing an underestimation of the meridional pressure gradient variability. In order to overcome this problem, estimated SSH data are bias corrected using statistical information from the observations. g sim (n) and g obs (n) are discrete values of simulated and observed SSH for a given period of time containing N time steps (1 B n B N). Further, O (g sim (n)) and O (g obs (n)) are defined as sorted discrete functions applied to g sim (n) and g obs (n), respectively. A third function F is defined by the relation F is calculated from the relation of O (g sim (n)) and O (g obs (n)) using a polynomial function as approximation. We choose a 3rd order polynomial function with coefficients estimated from a simple optimisation method. The bias corrected g sim (n) is given as The variability of g sim-corr (n) is much closer to that of g obs (n) and the correlation between estimated and observed SSH is slightly larger. Using g sim-corr (n) instead of g sim (n) as forcing at the lateral boundary in Kattegat improves the simulated Baltic Sea salinity during present climate. The agreement between the probability density functions of the reconstructed and corrected SSH and the observations is very good (not shown).
In the transient simulations we applied the correction both in past and future climates assuming that the statistical relationship will not change with time.

Runoff
Runoff is calculated with a statistical method which is applied to estimate river flows from the net water budget (precipitation minus evaporation) over the Baltic drainage area as simulated with RCAO because for our experiments results from a hydrological model were not available.
The net water budget in RCA is realistically simulated Kjellström and Lind 2009). For the scenario simulations only the variability of annual mean runoff anomalies is calculated. We do not consider changes of the seasonal cycle of the runoff because their impact on the large-scale salinity distribution in the Baltic Sea is small (Meier and Kauker 2003).
Our method assumes that the annual mean runoff from a given drainage area p during the year n is correlated with the net water budget anomaly (in %) over this given water area during the given year and the one before: in which R p,n is the runoff for the year n and for the drainage area p. B p,n is the net water budget (precipitation minus evaporation) anomaly for year n and area p. Finally, b p and a p are two coefficients. Five different drainage areas are considered, i.e. the drainage areas of the Bothnian Bay, Bothnian Sea, Gulf of Finland, Baltic proper and Kattegat (Fig. 1). The statistical model is constrained for present climate when reliable observations of the annual mean runoff anomaly are available (Bergström and Carlsson 1994). b p and a p are determined using an optimisation method during 1980-2006. During 1960-1979 the statistical model is validated when both runoff observations and simulation results from RCA driven by ERA40 are available. Figure 2 shows the results of the statistical model for 1960-2006 for the Baltic proper. The results are satisfactory except for the Gulf of Finland and the Gulf of Riga (Table 1). The annual variability is fairly well reproduced for the entire Baltic Sea although it is obvious that the standard deviation of the re-constructed runoff is smaller than the standard deviation of the observations (not shown).
It is assumed that the statistical relationship between runoff and precipitation minus evaporation does not change in time.

Nutrient loads scenarios
Nutrient loads from rivers are calculated from the product of the nutrient concentration and the volume flow (Sect. 2.5) following Stålnacke et al. (1999), see Eilola et al. (2009) andMeier et al. (2011b). Thus, it is assumed that the nutrient reservoir on land will be large enough to provide increased nutrient loads during the integration period if the volume flows increase. Four scenarios are considered:  A summary of the assumptions behind these nutrient load scenarios can be found in HELCOM (2007) based upon Wulff et al. (2007) and Humborg et al. (2007), while a comprehensive description of the calculation of the actual load changes is given in Gustafsson et al. (2011). The averaging period for the reference river load concentration is 1995-2002. This excludes some years in the early 2000s with abnormal nutrient loads. Load changes are applied on the total loads (not only on bioavailable fractions). Table 2 shows the nutrient concentration changes beyond the year 2020 for the different scenarios calculated by Gustafsson et al. (2011).
For the transient scenario simulations the future nutrient input into the Baltic Sea is represented by piecewise linear ramp functions. We run RCO-SCOBI until the end of 2007, ramp to the end of 2020 and then use constant nutrient concentrations in the rivers according to BSAP, CLEG and BAU (for REF the nutrient concentrations are constant with time). Coastal point sources are lumped into the river loads. The same functional form is used for the atmospheric deposition of nutrients.
These socio-economic scenarios are combined with the four climate scenarios described in Sect. 2.3. The 16 transient simulations are summarized in Table 3. Following Meier (2006) and Meier et al. (2011b), in all scenario simulations lateral boundary conditions in the northern Kattegat are unchanged assuming that especially Kattegat deep water properties will not change in future.

Cause-and-effect studies
In addition to scenario simulations, we performed causeand-effect studies with climatological mean seasonal cycles for air temperature and specific humidity fields or with climatological mean river discharges after year 2007. Thus, in the first set of experiments (henceforth TAIR-CLIM) water temperatures do not increase during the twenty first century compared to the reference period 1969-1998. In the second set of experiments (henceforth RUNOFFCLIM) salinities do not decrease and nutrient transports from land only alter with changing nutrient concentrations depending on the applied socio-economic scenario because annual mean volume flows do not change with time.

Analysis strategy
We focus on annual and seasonal mean changes of selected climate parameters and ecological quality indicators describing changes of the atmospheric and hydrological forcing and changes of the physical and environmental status of the Baltic Sea like air temperature, maximum estimated gust wind, river discharge, nutrient loads from land to the sea, water temperature, salinity, and oxygen, phosphate, nitrate and phytoplankton concentrations.
Further, Secchi depth (S d ) is calculated from S d = 1.7/k(PAR), where k(PAR) is the coefficient of underwater attenuation of the photosynthetically available radiation    ECHAM5-r1-A1B 1/6 1/6 1/6 1/6 ECHAM5-r1-A2 1/6 1/6 1/6 1/6 The abbreviations are explained in the text Transient scenario simulations for the Baltic Sea 2425 (Kratzer et al. 2003). Factors controlling light attenuation in the Baltic Sea model are the concentrations of phytoplankton and detritus. Thereby, the total phytoplankton concentration is the sum of the three algal groups represented in SCOBI (Sect. 2.2). In the scenario simulations changes of the Secchi depth are given by changing phytoplankton and detritus concentrations. For the calculation of Secchi depth it is assumed that yellow substances (or colored dissolved organic material, CDOM) are constant in time and will not increase in future. This assumption is probably not true because available observations indicate that in regions with salinities below 2-3 g kg -1 concentrations of CDOM are considerably larger than in other regions with higher salinities. Hence, it might be that decreasing salinities in future will cause increasing concentrations of yellow substances in the gulfs and in coastal zones affecting light climate for phytoplankton. Thus, there might be a non-linear impact of increasing yellow substances and changing phytoplankton and detritus concentrations on Secchi depth. As details of processes and future trends are unknown, the impact of possible changes of yellow substances is not considered in our scenario simulations.
As wind speed extremes in regional climate models like RCAO are usually underestimated, a modification of the 10 m wind is applied to guarantee correctly simulated wind induced mixing in the ocean (Meier et al. 2011c). Following Meier et al. (2011c) both the maximum 10 m wind speed and the maximum estimated gust wind are used to calculate wind speed extremes. In RCAO the maximum 10 m wind speed is calculated following the Monin-Obukhov similarity theory (Monin and Obukhov 1954) and is interpolated from the lowest atmospheric level (90 m) down to 10 m. The maximum estimated gust wind is calculated from the turbulent kinetic energy equation following Brasseur (2001). In RCAO the gust winds can propagate down to the surface from all boundary layer levels if the mixing is strong enough. For both parameters the absolute maximum over the output interval of three hours is stored while the internal time step is 15 min for 25 km resolution (Samuelsson et al. 2011). In general, the estimated gust wind is larger than the maximum 10 m wind speed. Thus, in this study we focus on the analyis of changes of the mean maximum estimated gust wind.
To analyse the response of ocean parameters we calculated ensemble mean fields from the four climate projections. As our ensemble consists of one HadCM3 and three ECHAM5 driven simulations, we weighted the ensemble members such that both GCMs are equally represented within the ensemble mean (Table 3). As we have only four ensemble members, the ensemble spread is illustrated by the difference between the maximum and minimum values of the four simulations. As the ensemble is too small, the calculation of the standard deviation would not be a good measure to characterize the spread.

Evaluation of the control climate
The quality of the atmospheric forcing was assessed by Meier et al. (2011c). Meier et al. (2011a) compared results of the climate simulations during the control period with hindcast simulations forced with regionalized ERA40 data (Samuelsson et al. 2011). Here, we focus on ocean modeling results of the four transient simulations at selected monitoring stations in the Baltic proper and Gulf of Finland (Fig. 1). Figure 3 shows simulated and observed mean vertical profiles of water temperature, salinity, oxygen, phosphate, and nitrate concentrations at Bornholm Deep (BY5), Gotland Deep (BY15), Landsort Deep (BY31) and in the western Gulf of Finland (LL07). For the control period 1978-2007 the agreement between model results and observations is satisfactory. In ECHAM5 driven simulations the resulting mean vertical profiles of salinity, oxygen, phosphate and nitrate are close to observations. Most of the biases are within the range of natural variability of the observations. The quality of the HadCM3 driven simulation is worse although mean vertical temperature profiles are slightly better reproduced in the HadCM3 driven simulation than in ECHAM5 driven simulations because of a warm bias in ECHAM5 over the Baltic Sea (Meier et al. 2011c). Table 4 summarizes the mean temperature and precipitation changes of RCAO over the Baltic Sea region. The largest increase of the 2 m air temperature is found in the northern Baltic Sea in particular during winter and spring (Fig. 4). Within the ensemble the largest increase occurs in the HadCM3 driven simulation with the A1B greenhouse gas emission scenario and not in the ECHAM5 driven simulation with the A2 emission scenario (not shown). Meier et al. (2011c) showed that a warm bias of the control climate  in ECHAM5 driven simulations reduces the ice-albedo feedback. Hence, the climate change signal is smaller compared to the HadCM3 driven simulation. In the three ECHAM5 driven simulations the horizontal patterns of the air temperature changes are very similar. Thus, we depict only results from RCAO-ECHAM5-r1-A1B.  Although regional details and the overall magnitude may differ, in all experiments the sea level pressure (SLP) will get more zonal at the end of the century (not shown). These results of RCAO are in accordance with results by . The largest changes of precipitation occur over the mountain areas (not shown). We found similar patterns of changing precipitation in HadCM3 and ECHAM5 driven simulations. Cloud cover changes are relatively small (not shown). During spring the cloudiness will slightly increase in ECHAM5 driven simulations. In the other seasons the cloudiness will slightly decrease over the Baltic Sea. This is a common signal in all scenario simulations.

Atmospheric variables
In general, also the changes of the mean 10 m wind speed are small (not shown). Significantly increased mean 10 m wind speeds of about 1 m s -1 are found only in RCAO-ECHAM5-r1-A1B and RCAO-ECHAM5-r1-A2 in the Baltic proper during winter and autumn and in RCAO-HadCM3-A1B in the Bothnian Sea during autumn.
In all simulations the maximum 10 m wind speed and the maximum estimated gust wind increase by more than 1 m s -1 in the Bothnian Bay and Gulf of Finland during winter and spring as a consequence of the melting sea ice in future climate (Fig. 5). In RCAO-HadCM3-A1B maximum changes exceed even 3 m s -1 . However, compared to earlier results of ECHAM4 (Roeckner et al. 1999) driven regionalizations (e.g. Räisänen et al. 2004) are wind speed changes in the simulations of this study small. Hence, its consequences for stratification and oxygen conditions in the Baltic deep water are smaller than reported by Meier et al. (2011b) (see below). Table 5 summarizes the total volume flows in present and future climates calculated with the statistical model (Sect. 2.5). We found changes of the total volume flow to the Baltic Sea between 15 and 22 % In all sub-basins including the Baltic proper the volume flow changes are positive (not shown). The changes are larger in the northern sub-basins than in the southern subbasins in agreement with earlier results by Graham (2004). However, positive volume flow changes to the Baltic proper are an important difference compared to the results  by Meier et al. (2011b). In their scenario simulations volume flow changes to the Baltic proper were negative although in three out of four climate projections the total flows into the entire Baltic were positive. As nutrient loads from the Baltic proper catchment area are largest, volume flow changes to the Baltic proper control the changes of the total nutrient supply from land (assuming reference river nutrient concentrations). Thus, in the experiments by Meier et al. (2011b) the total nutrient load changes have different sign depending on the selected climate projection whereas in our study nutrient loads increase in all projections (Fig. 6). This discrepancy reflects the uncertainty of state-of-the-art hydrological modelling (Graham et al. 2007). In Fig. 6 the absolute changes of the nutrient loads per sub-basin are shown. In the nutrient load scenario REF the biologically available phosphorus and nitrogen loads increase by about 15-20 % each depending on the climate projection. These changes are about two times larger than the largest increase calculated by Meier et al. (2011b).

Runoff and nutrient loads
In the most optimistic nutrient load scenario of this study (BSAP) we found reductions of about 7-8 and 10-50 kton year -1 for phosphorus and nitrogen, respectively (Fig. 6). However, the BSAP requires much larger reductions of 15 and 133 kton year -1 for phosphorus and nitrogen, respectively (HELCOM 2007). Note that the latter figures for the total Baltic include Kattegat. In our simulations nutrient load reductions in the Baltic proper are partly compensated by nutrient load increases in the northern sub-basins, e.g. in the Bothnian Sea and Bothnian Bay. The reductions of the BSAP scenario by Meier et al. (2011b) are closer to the original figures by HELCOM (2007).
For the most pessimistic nutrient load scenario of this study (BAU) we found increases of about 17-19 and 440-500 kton year -1 for phosphorus and nitrogen, respectively (Fig. 6). Following HELCOM (2007) changes according to BAU amount to 16 and 340 kton year -1 for phosphorus and nitrogen, respectively.
Note that in the nutrient load scenario CLEG the total phosphorus load will decrease and the total nitrogen load will increase in future climate (Fig. 6). As this scenario is novel our figures cannot be compared with other studies.

Water temperature and salinity
In all scenario simulations the volume averaged water temperature increases with time as a response of the increased air temperature and the volume averaged salinity decreases as a response of the increased runoff during the twenty first century (Fig. 7).
In all scenario simulations sea surface temperature (SST) changes between 2070-2099 and 1978-2007 are largest in the Bothnian Bay and Bothnian Sea during summer (Fig. 8). This pattern is a common feature of all members of our mini-ensemble although the magnitude of the warming differs substantially between HadCM3 and ECHAM5 driven simulations such that the ensemble spread is also largest in the Bothnian Bay and Bothnian Sea during summer (Fig. 8). However, the climate change signal is much larger than the uncertainty caused by the  GCMs as indicated both by the standard deviation of the ensemble mean (not shown) and the ensemble spread. We found the largest SST increases of more than 6°C in the southern Bothnian Bay in the HadCM3 driven scenario simulation during summer. In ECHAM5 driven simulations the largest SST increase is located in the central Bothnian Bay and does not exceed 4°C approximately. Further, in all scenario simulations the largest SST increase during winter and spring occurs in the Gulf of Finland.
Noteworthy, SST changes lag 2 m air temperature changes in time. Maximum air temperature changes occur in winter and spring in the Bothnian Bay, the northern Bothnian Sea and the eastern Gulf of Finland (Fig. 4). However, the largest SST signal is found in the Bothnian Bay during summer (Fig. 8). The larger warming in the northern compared to the southern Baltic is explained by the ice-albedo feedback (Meier et al. 2011b;Wang et al. 2005) and by the reduced thermal convection under warmer and fresher conditions (Hordoir and Meier 2011).
Also spatial patterns of sea surface salinity (SSS) projections show an overall agreement with largest decreases in the Baltic proper of about 1.5-2 g kg -1 (not shown). Salinity is reduced because in all scenario simulations runoff is significantly increased. Changes of the wind speed are of minor importance for SSS changes. Largest discrepancies between scenario simulations are found for the SSS projections in Kattegat.
Changes of bottom salinity concentrations follow SSS changes (not shown). As in the ECHAM5-r1-A1B and A2 driven simulations maximum wind speeds increase in particular over the Bothnian Bay and Gulf of Finland but also over the Baltic proper (Fig. 5), wind induced mixing is larger in future climate and the permanent halocline is deeper located. Consequently, we found largest bottom salinity changes along the slopes of the Baltic proper and Gulf of Finland at depths of the halocline changes. As the wind changes occur only in two out of four scenario simulations of our mini-ensemble, the largest spread of projected bottom salinity is related to the different depth of the halocline in the four simulations.  (Fig. 9). In all nutrient load scenarios dissolved inorganic phosphorus (DIP) increases during the control period which is consistent with the hindcast simulation (Meier et al. 2011a). After 2007 DIP increases in REF and BAU and decreases in BSAP (about constant in CLEG). The results confirm the much longer time scale of the phosphorus cycle compared to the nitrogen cycle.
Interestingly, the sediment pools of nitrogen and phosphorus decrease in almost all scenarios towards the end of the twenty first century (Fig. 9). Only in BAU nitrogen and phosphorus concentrations in the sediments increase between 2007 and 2050 approximately. After 2050 both concentrations either decrease (nitrogen) or remain relatively constant (phosphorus).
In all scenarios there is a tendency of increased DIN to DIP ratio in the water column (not shown). Especially in the BSAP scenario driven by ECHAM5-r1-A2 and ECHAM5-r1-A1B (the scenario simulations with an increase of the wind speed over the Baltic proper, see Fig. 5) the overall DIN to DIP ratios increase with about 7 and 5, respectively. Thus, the Baltic proper will switch from a nitrogen limited to a phosphorus limited ecosystem in agreement with the results by Meier et al. (2011b).

Projected changes for REF
In the reference scenario (REF) the bottom oxygen concentrations decrease in all climate projections in almost all regions (Fig. 10). Exceptions are the deep water in the Gulf of Finland and regions along the slopes of the Gotland Basin where the stratification decreases due to a deeper halocline caused by increased runoff (in all climate projections, Table 5) and increased wind speed (in ECHAM5-r1-A1B and ECHAM5-r1-A2 driven simulations, Fig. 5). As the oxygen saturation concentration is smaller in warmer water, in the coastal zone with only weak vertical stratification both the surface and bottom oxygen concentrations decrease slightly. The decrease is larger in regions with larger water depth and with a permanent halocline (Fig. 10). We found the largest decrease of bottom oxygen concentrations in the HADCM3 driven simulation in the central area of the deep Bornholm Basin, Gotland Basin and Bothnian Sea (not shown). In the ensemble mean the bottom oxygen concentration in the Gotland Deep area decreases by more than 1.5 ml l -1 (Fig. 10).
The ensemble spread is largest in regions that are affected by the varying position of the halocline (not shown). Note that the ensemble spread does not differ significantly between the scenarios.
As the phosphorus release capacity of the sediments is oxygen dependent (Eilola et al. 2009), the generally decreased bottom oxygen concentrations cause a decreased storage of phosphorus in the sediments and increased phosphate concentrations in the surface waters (not shown). We found the largest phosphate concentration increase in the HadCM3 driven simulation in the Baltic proper and Gulf of Finland. In the ensemble mean the largest increase of surface phosphate concentrations occurs in the southern Baltic proper (Arkona Basin, Bornholm Basin and southern Gotland Basin) during winter. This signal is a common pattern in all climate projections. The surface phosphate concentration changes in the Gulf of Finland during spring have the largest spread within our ensemble.
In all climate projections the surface nitrate concentration remains unchanged or increases (not shown). The patterns of changing nitrate concentration are similar in all climate projections. Especially during winter and especially in the eastern Gulf of Finland, Gulf of Riga and along the eastern coasts of the Gotland Basin nitrate concentrations increase in future climate. The increased supply of nitrogen from the rivers and the increased oxygen concentrations in the Gulf of Finland (causing decreased denitrification) might be the reasons for the increased nitrate concentrations particularly in the coastal zone close to the river mouths of the large rivers.
Both increased temperatures and increased concentrations of nitrate and phosphate during winter affect the spring and summer blooms (not shown). In particular during spring the concentrations of flagellates and others increase in the eastern Baltic proper, Gulf of Riga and Gulf of Finland. During summer and autumn the concentrations of cyanobacteria increase in the southern Baltic proper in all climate projections. In the Gulf of Finland cyanobacteria blooms are also more intensive in the HadCM3 driven simulation whereas in ECHAM5 driven simulations no significant changes are found. Both flagellate and others and cyanobacteria changes affect phytoplankton concentrations in future climate (Fig. 10).
As a consequence of simulated phytoplankton and detritus changes, Secchi depth in the southern Baltic proper decreases particularly during summer and autumn (Fig. 10). In the ensemble mean the largest Secchi depth changes exceed 1 m.

Projected changes for BSAP
As mentioned in the previous section the bottom oxygen concentration decreases significantly in the HadCM3 Hence, in the BSAP scenario the improvements of the nutrient load reductions will be counteracted by the effect of changing climate at the end of the century. As a consequence, in the HadCM3 driven simulation under the BSAP scenario bottom oxygen concentration changes are only small (not shown). In the ECHAM5 driven simulations we found increased bottom oxygen concentrations in the Gulf of Finland and in the Gotland Basin when we applied the BSAP. In the ECHAM5-r1-A1B and A2 driven simulations we found the largest increases of the bottom oxygen concentration along the slopes of the Gotland Basin and in the Gulf of Finland due to the deeping of the halocline and the corresponding decreased stratification in that depth interval. While in the HadCM3 driven simulation surface phosphate and nitrate concentration changes are small, we found in ECHAM5 driven simulations in the Gulf of Finland reduced surface phosphate and increased surface nitrate concentrations. Thus, changes of surface nutrient concentrations are largely controlled by changing bottom oxygen concentrations (Savchuk 2010). As shown below water temperature changes affect the decomposition of organic matter in the sediments and contribute to changes of surface nutrient concentrations.
Surface concentration changes of diatoms, flagellates and others and cyanobacteria are diverse (not shown). During spring in all scenario simulations surface diatom concentrations decrease especially along the southern and eastern coasts of the Baltic proper and in the Gulf of Finland. To the contrary, we found slight increases of the surface concentrations of flagellate and others mainly in the Gulf of Finland. Cyanobacteria concentrations increase in the southern Baltic proper (mainly in the Bornholm Basin) in the HadCM3 driven simulation and remain basically unchanged in ECHAM5 driven simulations. During spring surface phytoplankton concentrations in the ECHAM5 driven simulations decrease following diatom concentration changes. During summer we found slight surface phytoplankton concentration increases in the southern Baltic proper in the HadCM3 driven simulation following cyanobacteria concentration changes. Corresponding increases of Secchi depth during spring amount to about 1 m at maximum. During summer Secchi depth in the Had-CM3 driven simulation decreases by about 0.5 m at maximum.
In the ensemble mean the summer oxygen concentrations in the deeper regions of the Baltic proper increase only slighthly (Fig. 10). However, in the Gulf of Finland the increases are larger exceeding in some regions 0.9 ml l -1 . In the Baltic proper the spring phytoplankton concentrations decrease in a narrow zone along the southeastern coasts, especially in Oder and Gdansk bays (Fig. 10). During spring moderate Secchi depth increases of about 0.6-1 m are found in the southern Baltic Sea (Arkona and Bornholm basins) (Fig. 10).

Projected changes for CLEG
In the HadCM3 driven simulation bottom oxygen concentrations decrease at the end of the century in CLEG almost everywhere (not shown). In the ECHAM5 driven scenarios the bottom oxygen concentrations especially in the Gulf of Finland increase. In the two scenarios with increased wind induced mixing in the Baltic proper (ECHAM5-r1-A1B and ECHAM5-r1-A2) the bottom oxygen concentrations along the slopes increase as well because of the deeper halocline. In the ensemble mean we found small decreases of bottom oxygen concentrations in the deeper parts of the Baltic proper and small increases along the slopes and in the Gulf of Finland (Fig. 10).
As a consequence of the bottom oxygen concentration changes surface phosphate concentrations increase in HadCM3 driven scenario simulations. Nitrate concentration changes are largest in the Gulf of Finland and in the Gulf of Riga in ECHAM5-r1-A1B and ECHAM5-r1-A2 driven simulations.
Concentration changes of diatoms, flagellates and others, and cyanobacteria are relatively small. As the projected phytoplankton concentrations in the ensemble mean slightly increases (with larger changes in the northern Baltic proper), Secchi depth decreases (Fig. 10). Largest changes of about 1 m are found in the Bornholm Basin in the HadCM3 driven scenario simulation during summer. In the ensemble mean the largest changes are found in the Archipelago Sea and in the entrance of the Gulf of Riga but they do not exceed 0.8 m (except in the river mouth region of the river Neva) (Fig. 10).

Projected changes for BAU
In the BAU scenario the impact of increased nutrient loads and the impact of changing climate have large consequences for the marine environment. Large reductions of bottom oxygen concentrations (Fig. 10), large increases of phosphate and nitrate concentrations at the surface of the Baltic proper in particular during winter (not shown) and large increases of both the spring and summer blooms characterize the BAU scenario (Fig. 10). In this scenario Secchi depth in the south-western Baltic at the end of the century is projected to be more than 2 m smaller compared to present conditions (Fig. 10).

Cause-and-effect studies
The impact of changing climate on the Baltic Sea biogeochemistry in the various nutrient load scenarios is studied further by employing the sensitivity experiments, TAIRCLIM and RUNOFFCLIM (Sect. 2.7). In the following we focus on scenario simulations driven by ECHAM5-r3-A1B combined with the nutrient load scenarios REF and BSAP (Figs. 11,12).
The increase of the average water temperature in the ECHAM5-r3-A1B driven simulation from about 6.5°C during 1969-1998 to about 8.5°C at the end of the century (caused by the increased air temperature over the Baltic Sea, see Table 4) enhances the decomposition of organic matter in the sediments. In the present version of SCOBI the time scale for decomposition will be shortened by a factor of two if the bottom water temperature increases by about 5°C. However, feedback mechanisms may contribute to a further acceleration of the decomposition (see Sect. 4).
Due to the rising water temperature in REF the concentrations of nitrogen and phosphorus in the sediments decrease (Fig. 11) causing increased nutrient concentrations in the surface layer that are available for the production of organic material during the spring bloom (not shown). The effect is larger for phosphorus than for nitrogen because other processes, e.g. denitrification, change as well and modify the nitrogen cycle. This might also be the reason that in the BSAP scenario the pool of DIN in the water column decreases with increasing water temperature instead of an increase as in the REF scenario (Fig. 11).
Increased nutrient concentrations in the surface layer during winter result in increased phytoplankton concentrations during spring and summer and in decreased bottom oxygen concentrations during all seasons due to increased decomposition of organic material in the deep water (Fig. 12). For instance, the bottom oxygen concentrations in the ECHAM5-r3-A1B driven simulation in the REF scenario decrease in the deeper regions of the eastern Gotland Basin by about 0.9-1.2 ml l -1 . In the corresponding TAIRCLIM scenario simulation without any changes of the water temperature the bottom oxygen concentrations decrease in the same regions by only 0.3-0.6 ml l -1 (due to the increased nutrient supply, see Fig. 6).
In the BSAP scenario the bottom oxygen concentrations in the deeper parts of the eastern Gotland Basin increase by about 0.3-1.2 ml l -1 and are not affected by the changing water temperatures (compare BSAP and BSAP-TAIRCLIM in Fig. 12). However, in the Gulf of Finland the increase of the bottom oxygen concentrations is about 0.6 ml l -1 larger in TAIRCLIM compared to the scenario simulation with changing water temperature. In addition to the accelerated decomposition of organic material in the sediments the oxygen saturation concentration decreases in warmer water. However, this effect is smaller than the discussed increase of the oxygen consumption and causes changes of 0.3-0.4 ml l -1 in maximum.
The results of the sensitivity experiments with climatological mean runoff (RUNOFFCLIM, not shown) confirm that increased phytoplankton and decreased bottom oxygen concentrations are partly also explained by the increased nutrient supply from land (Fig. 6). Note that in our experiments the increased nutrient loads are caused by increased volume flows due to the increased net precipitation over land (Table 5).

Discussion
Projections of the impact of climate change on the Baltic ecosystem are not forecasts of future conditions for the marine environment, like weather forecasts, but depend on additional assumptions on the future socio-economic development both worldwide and in the Baltic Sea region.
For instance, future greenhouse gas emissions, nutrient loads from rivers to the sea, nutrient emissions from point sources and atmospheric deposition of nutrients are important drivers of both the regional climate and the marine environment. Hence, scenarios need to be developed that take possible socio-economic changes into account, like changes of the human population, life-style (including food and energy consumption), agricultural practices, etc.
In this study we applied the detailed greenhouse gas emission scenarios developed by the Intergovernmental Panel on Climate Change (IPCC) (Nakićenović et al. 2000). However, consistent socio-economic scenarios on both global and regional scales are not available. Hence, we implemented nutrient load scenarios that follow policy scenarios for the Baltic Sea suggested by HELCOM (2007). These scenarios are not consistent with the global scenarios used for the GCM simulations but are to our knowledge the only available scenarios for the Baltic Sea region. Hence, it is obvious that due to the large range of plausible nutrient load scenarios the uncertainties of our projections are large. Other uncertainties of our projections are related to biases of global and regional climate models and have been discussed earlier, e.g. by , , , Räisänen et al. (2004), Meier et al. (2011c).
Especially, projected changes of the surface wind fields are rather uncertain for the Baltic Sea region . As wind speed changes are important drivers of future changes of the vertical stratification in the Baltic Sea, the uncertainties affect projections of the marine ecosystem considerably. For instance, in ECHAM5 driven simulations (e.g. Neumann 2010, this study) the wind speed changes are much smaller than in ECHAM4 driven simulations (e.g. Meier et al. 2011b) causing large differences in bottom oxygen concentration changes and subsequent changes of the nutrient cycles in the Baltic.
Further, the functioning of some key processes of the coupled physical-biogeochemical system are not well known . These shortcomings affect the sensitivity of the model response to changing climate and changing nutrient loads (Meier et al. 2011b). For instance, the processes controlling nutrient fluxes between the bottom water and the sediments are poorly understood (e.g. Eilola et al. 2009;Almroth-Rosell et al. 2011).
The results of this study suggest that nutrient loads in future climate may increase due to increased runoff. As the calculated changes are larger than those reported by Meier et al. (2011b) the effects on ecological quality indicators, like bottom oxygen concentration, phytoplankton concentration and Secchi depth (Fig. 10) are larger than in the scenario simulations by Meier et al. (2011b). Neumann (2010) calculated nutrient load changes with the same approach as in our study following Stålnacke et al. (1999). He found similar large increases of the freshwater budget of about 20 % in his two scenario simulations (A1B and B1) than in our experiments. However, in his simulations the oxygen concentration of the deep water improves slightly causing a reduction of future hypoxic areas. In contrast, we found decreasing bottom oxygen concentrations in all scenarios with increasing nutrient loads. Thus, the sensitivity of different models to changing nutrient loads seems to differ significantly.
In the scenario simulations of this study another process is important which may affect the marine ecosystem significantly. Increased water temperatures cause not only reduced oxygen saturation concentrations, but also (and perhaps even more important) an increased decomposition of organic matter in the sediments. The presented results suggest that in future climate less nutrients will be stored in the sediments due to the increased remineralization under higher temperatures (Fig. 11).
In SCOBI the parameters describing the nutrient fluxes from the sediments were chosen such that the time scale for decomposition will be shortened by a factor of two if the bottom water temperature increases by about 5°C. Very likely the time scale is even shorter because of a feedback mechanism due to changing bottom oxygen concentrations. Increased remineralization of phosphorus will increase the phytoplankton production which will increase the oxygen consumption which in turn will result in lower bottom oxygen concentrations. Thus, the phosphorus retention capacity will be lower causing lower/higher concentrations of phosphorus in the sediments/water column (Eilola et al. 2009). For nitrogen the processes are even more complex because increased decomposition will also affect denitrification.
In this study the parameterization of the temperature dependent remineralization of the Ecological Regional Ocean Model (ERGOM) (Neumann et al. 2002;Neumann and Schernewski 2008) was used. Thus, the sensitivity to temperature changes is about 30 times larger than in the scenario simulations by Meier et al. (2011b) which explains the larger impact of increased decomposition on the results of our scenario simulations. As the true remineralization rate and its temperature dependency are unknown, the role of changing remineralization in future climate remains uncertain.

Summary and conclusions
In this study we focussed on annual and seasonal mean changes of ecological quality indicators like bottom oxygen concentration, phytoplankton concentration and Secchi depth describing the environmental status of the Baltic Sea. Agreement and disagreement of the simulated changes were assessed from the statistics of an ensemble of 16 scenario simulations. Projected changes at the end of the twenty first century are usually larger than biases induced by the deficiencies of GCMs at the regional scale. Especially ensemble mean biases are smaller than ensemble mean changes stressing the added value of ensemble modelling.
In agreement with earlier studies (Meier 2006;Meier et al. 2006Meier et al. , 2011bNeumann 2010) we found that at the end of the twenty first century water temperature will increase, runoff will increase, and salinity will decrease. However, in contrast to earlier ECHAM4 driven simulations (Meier et al. 2011b) our results do not indicate considerable wind speed and gustiness changes over the Baltic proper. Hence, stratification changes due to wind induced mixing changes are relatively small.
In correspondence with earlier studies we found that the impact of changing climate on the Baltic biogeochemistry might be significant. The model simulations suggest that in addition to eutrophication changing climate is an important Transient scenario simulations for the Baltic Sea 2439 stressor for the Baltic ecosystem. With the help of sensitivity experiments we identified two processes that in our model are mainly responsible for environmental changes due to changing climate. Due to increased temperature and increased net precipitation in the Baltic catchment area the decomposition of organic material in the sediments will be accelerated and the nutrient loads from land will increase. Both processes cause increased nutrient concentrations in the surface layer and consequently an acceleration of eutrophication in the Baltic Sea. According to our scenario simulations with reference loads water quality (measured by bottom oxygen concentrations, phytoplankton concentrations and Secchi depths) will be deteriorated in future climate. For instance, in summer the ensemble mean of the Secchi depth will decrease in the southern Baltic proper by up to 1.5 m.
According to our results nutrient load reductions performed under current legislation will not be sufficient to improve the water quality at the end of the century. The climate effect is larger than the impact of nutrient load reductions and Secchi depth will decrease especially in the southern Baltic proper. The larger nutrient load reductions of the BSAP will improve the water quality at the end of the century. However for the same targets larger reductions will be necessary compared to present climate. In summer the ensemble mean of the Secchi depth will increase in the southern Baltic proper by about 1 m in maximum. In case of an exponential growth of agriculture following a pessimistic business-as-usual scenario bottom oxygen concentrations will decrease, surface nutrient concentrations will increase and Secchi depth will decrease significantly. During the warmer seasons (spring to autumn) the ensemble mean of the Secchi depth in the southern Baltic proper will decrease by more than 2 m in some regions.