Description: Numerical Earth System Models are generic tools used to extrapolate present climate conditions into a warming future and to explore geoengineering options. Most of the current-generation models feature a simple pelagic biogeochemical model component that is embedded into a three-dimensional ocean general circulation model. The dynamics of these biogeochemical model components is essentially controlled by so-called model parameters most of which are poorly known. Here we explore the feasibility to estimate these parameters in a full-fledged three-dimensional Earth System Model by minimizing the misfit to noisy observations. The focus is on parameter identifiability. Based on earlier studies, we illustrate problems in determining a unique estimate of those parameters that prescribe the limiting effect of nutrient- and light-depleted conditions on carbon assimilation by autotrophic phytoplankton. Our results showcase that for typical models and evaluation metrics no meaningful “best” unique parameter set exists. We find very different parameter sets which are, on the one hand, equally consistent with our (synthetic) historical observations while, on the other hand, they propose strikingly differing projections into a warming climate.

Description: [1] The North Atlantic Oscillation (NAO) index is a frequently used measure for the mean winter conditions in Northern Europe. A positive, high index is associated with strong westerlies and anomalous warm temperatures. The effects on sea ice conditions in the Baltic Sea are twofold. Warm temperatures prevent sea ice formation. If ice is present nevertheless, the strong winds can promote the formation of ice ridges which hinders ship traffic. We use an ocean-sea ice model to investigate the NAO impact on the ridged ice area fraction in the Baltic during 1962–2007. Our simulations indicate that in the northern Bothnian Bay, a high NAO index is related to an anomalous accumulation of ridges, while in the rest of the Baltic Sea, the relationship is contrary. The NAO explains locally at most only 20–25% of the ridged ice fraction interannual variability which indicates the systems complexity. However, we find high skill with local correlations around 0.8 for annually averaged ridged ice fraction reconstructed from multilinear regression using winter averaged wind extremes, surface air temperature, and sea surface temperature (SST). This suggests that the amount of ridged ice in late winter can be derived from these routinely measured quantities. In large parts of the basin, it is sufficient to use the atmospheric parameters as a predictor, while in the eastern Bothnian Bay and southern Gulf of Finland, the SST is required to reconstruct the bulk of the ridged ice fraction.

Description: The multidecadal variability of air-sea CO(2)fluxes in the North Atlantic under preindustrial atmospheric CO(2) conditions is simulated, using a coupled biogeochemical/circulation model driven by long-term surface forcing reconstructed from the leading modes of sea level pressure observations from 1850 to 2000. Heat fluxes are of great importance for the multidecadal CO(2) fluctuations, about equal in magnitude to wind stress, in contrast to their less prominent role for CO(2) flux variability on interannual timescales. Another difference, compared to higher frequencies, is the dominance of the North Atlantic Oscillation in driving the variability of the air-sea CO(2) fluxes. Two spatially distinct regimes lead to large anomalies in the CO(2) fluxes but compensate to a large degree. The first regime is advective and has its clear signature southeast of Greenland while the second one, in the vicinity of the Labrador Sea and off Newfoundland, is convective. In both regimes, the multidecadal CO(2) fluctuations are driven mainly by variations in temperature, salinity, and DIC content at the sea surface while the role of the biological pump is of minor importance in this particular model. The magnitude of the simulated multidecadal CO(2) uptake changes is on the order of 0.02 Pg C/yr and amounts to 10-15% of the estimated annual anthropogenic CO(2) uptake of the North Atlantic.

Description: On the basis of integrations of an eddy-permitting coupled physical-biological model of the tropical Pacific we explore changes in the simulated mean circulation as well as its intraseasonal to interannual variability driven by the biologically modulated vertical absorption profiles of solar radiation. Three sensitivity ocean hind-cast experiments, covering the period from 1948 to 2003, are performed. In the first one, simulated chlorophyll affects the attenuation of light in the water column, while in the second experiment, the chlorophyll concentration is kept constant in time by prescribing an empirically derived spatial pattern. The third experiment uses a spatially and temporally constant value for the attenuation depth. The biotically induced differential heating is generated by increased absorption of light in the surface layers, leading to a surface warming and subsurface cooling. The effect is largest in the eastern equatorial Pacific. However, the initial vertical redistribution of heat leads to considerable changes of the near-surface ocean circulation subsequently influencing the near-surface temperature structure. In general, including biophysical coupling improves the model performance in terms of temperature and ocean circulation patterns. In particular, the upwelling in the eastern equatorial Pacific is enhanced, the mixed layer becomes shallower, the warm bias in the eastern Pacific is reduced, and the zonal temperature gradient increases. This leads to stronger La Niña events and an associated increase in the variability of the Niño3 SSTA time series. Furthermore, the eddy kinetic energy (EKE) associated with mesoscale eddies in the eastern equatorial Pacific increases by almost 100% because of enhanced EKE production due to enhanced horizontal and vertical shear of the mean currents.

Description: The Baltic Sea is a marginal sea, located in a highly industrialized region in Central Northern Europe. Saltwater inflows from the North Sea and associated ventilation of the deep exert crucial control on the entire Baltic Sea ecosystem. This study explores the impact of anticipated sea level changes on the dynamics of those inflows. We use a numerical oceanic general circulation model covering both the Baltic and the North Sea. The model successfully retraces the essential ventilation dynamics throughout the period 1961–2007. A suite of idealized experiments suggests that rising sea level is associated with intensified ventilation as saltwater inflows become stronger, longer, and more frequent. Expressed quantitatively as a salinity increase in the deep central Baltic Sea, we find that a sea level rise of 1 m triggers a saltening of more than 1 PSU. This substantial increase in ventilation is the consequence of the increasing cross section in the Danish Straits amplified by a reduction of vertical mixing

Description: The isotopic composition of Si in biogenic silica (BSi), such as opal buried in the oceans' sediments, has changed over time. Paleorecords suggest that the isotopic composition, described in terms of δ30Si, was generally much lower during glacial times than today. There is consensus that this variability is attributable to differing environmental conditions at the respective time of BSi production and sedimentation. The detailed links between environmental conditions and the isotopic composition of BSi in the sediments remain, however, poorly constrained. In this study, we explore the effects of a suite of offset boundary conditions during the Last Glacial Maximum (LGM) on the isotopic composition of BSi archived in sediments in an Earth System Model of intermediate complexity (EMIC). Our model results suggest that a change in the isotopic composition of Si supply to the glacial ocean is sufficient to explain the observed overall low(er) glacial δ30Si in BSi. All other processes explored trigger model responses of either wrong sign or magnitude or are inconsistent with a recent estimate of bottom water oxygenation in the Atlantic Sector of the Southern Ocean. Caveats, mainly associated with generic uncertainties in today's pelagic biogeochemical modules, remain.

Description: Climate change is especially strong in the region of the Arctic Ocean, and will have an important impact on its thermo-haline structure. We analyze the results of a hindcast simulation of a new 3D ocean model of the Arctic and North Atlantic oceans for the period 1970–2019. We compared the time period 1970–1999 with the time period 2010–2019. The comparison showed that there is a decrease of stratification between the two periods over most of the shallow Arctic shelf seas and in the core of the Transpolar Ice Drift. Fresh water inputs to the ocean surface decline, and inputs of momentum to the ocean increase, which can explain the decrease in stratification. The comparison also showed that the mixed layer becomes deeper during winter, in response to the weakened stratification owing to increased vertical mixing. The comparison of summer mixed layer depths between the two time periods follows a deepening pattern that is less evident. Regional exceptions include the Nansen Basin and the part of the Canadian Basin bordering the Canadian Archipelago, where the mixed layer shoals. Trends of freshwater fluxes imply that the changes of haline stratification in these regions are also influenced by other processes, for example, horizontal advection of fresh water, increased mixing and changes in the underlaying water masses. Runoff increase toward the Arctic Ocean can locally decrease but also increase salinity, and has an impact on stratification which can be explained by coastal dynamics. The results emphasize the non-linear nature of Arctic Ocean dynamics.

Description: We analyze an extensive set of global coupled biogeochemical ocean circulation models. The focus is on the equatorial Pacific. In all simulations, which are consistent with observed standing stocks of relevant biogeochemical species at the surface, we find spuriously enhanced (reduced) macronutrient (oxygen) concentrations in the deep eastern equatorial Pacific. This modeling problem, apparently endemic to global coupled biogeochemical ocean circulation models, was coined “nutrient trapping” by Najjar et al. (1992). In contrast to Aumont et al. (1999), we argue that “nutrient trapping” is still a persistent problem, even in eddy-permitting models and, further, that the scale of the problem retards model projections of nitrogen cycling. In line with previous work, our results indicate that a deficient circulation is at the core of the problem rather than an admittedly poor quantitative understanding of biogeochemical cycles. More specifically, we present indications that “nutrient trapping” in models is a result of a spuriously damped Equatorial Intermediate (zonal) Current System and Equatorial Deep Jets—phenomenon which await a comprehensive understanding and have, to date, not been successfully simulated.