Description: In a changing climate, marine pelagic biogeochemistry may modulate the atmospheric concentrations of climate-relevant species such as CO2 and N2O. To date, projections rely on earth system models, featuring simple pelagic biogeochemical model components, embedded into 3-D ocean circulation models. Most of these biogeochemical model components rely on the hyperbolic Michaelis–Menten (MM) formulation which specifies the limiting effect of light and nutrients on carbon assimilation by autotrophic phytoplankton. The respective MM constants, along with other model parameters, of 3-D coupled biogeochemical ocean-circulation models are usually tuned; the parameters are changed until a "reasonable" similarity to observed standing stocks is achieved. Here, we explore with twin experiments (or synthetic "observations") the demands on observations that allow for a more objective estimation of model parameters. We start with parameter retrieval experiments based on "perfect" (synthetic) observations which we distort, step by step, by low-frequency noise to approach realistic conditions. Finally, we confirm our findings with real-world observations. In summary, we find that MM constants are especially hard to constrain because even modest noise (10 %) inherent to observations may hinder the parameter retrieval already. This is of concern since the MM parameters are key to the model's sensitivity to anticipated changes in the external conditions. Furthermore, we illustrate problems caused by high-order parameter dependencies when parameter estimation is based on sparse observations of standing stocks. Somewhat counter to intuition, we find that more observational data can sometimes degrade the ability to constrain certain parameters.

Description: Deoxygenation in the Baltic Sea endangers fish yields and favours noxious algal blooms. Yet, vertical transport processes ventilating the oxygen-deprived waters at depth and replenishing nutrient-deprived surface waters (thereby fuelling export of organic matter to depth), are not comprehensively understood. Here, we investigate the effects of the interaction between surface currents and winds (also referred to as eddy/wind effects) on upwelling in an eddy-rich general ocean circulation model of the Baltic Sea. Contrary to expectations we find that accounting for current/wind effects does inhibit the overall vertical exchange between oxygenated surface waters and oxygen-deprived water at depth. At major upwelling sites, however, as e.g. off the south coast of Sweden and Finland, the reverse holds: the interaction between topographically steered surface currents with winds blowing over the sea results in a climatological sea surface temperature cooling of 0.5 K. This implies that current/wind effects drive substantial local upwelling of cold and nutrient-replete waters.

Description: Ocean eddies can both trigger mixing (during their formation and decay) and effectively shield water encompassed from being exchanged with ambient water (throughout their lifetimes). These antagonistic effects of eddies complicate the interpretation of synoptic snapshots typically obtained by ship-based oceanographic measurement campaigns. Here we use a coupled physical–biogeochemical model to explore biogeochemical dynamics within anticyclonic eddies in the eastern tropical South Pacific Ocean. The goal is to understand the diverse biogeochemical patterns that have been observed at the subsurface layers of the anticyclonic eddies in this region. Our model results suggest that the diverse subsurface nutrient patterns within eddies are associated with the presence of water masses of different origins at different depths.

Description: The Southern Ocean is a major sink for anthropogenic carbon. Yet, there is no quantitative consensus about how this sink will change when surface winds increase (as they are anticipated to do). Among the tools employed to quantify carbon uptake are global coupled ocean-circulation biogeochemical models. Because of computational limitations these models still fail to resolve potentially-important spatial scales. Instead, processes on these scales are parameterized. There is concern that deficiencies in these so-called eddy-parameterizations might imprint wrong sensitivities of projected oceanic carbon uptake. Here, we compare natural carbon uptake in the Southern Ocean simulated with contemporary eddy-parameterizations. We find that very differing parameterizations yield surprisingly similar oceanic carbon in response to strengthening winds. In contrast, we find (in an additional simulation) that the carbon uptake does differ substantially when the supply of bioavailable iron is altered within its envelope of uncertainty. We conclude that a more comprehensive understanding of bioavailable iron dynamics will substantially reduce the uncertainty of model-based projections of oceanic carbon uptake.

Description: Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. "Offline" numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the "transport matrix method" (TMM), which represents tracer transport as a sequence of sparse matrix–vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their "online" counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

Description: We present Nemo-Nordic, a Baltic and North Sea model based on the NEMO ocean engine. Surrounded by highly industrialized countries, the Baltic and North seas and their assets associated with shipping, fishing and tourism are vulnerable to anthropogenic pressure and climate change. Ocean models providing reliable forecasts and enabling climatic studies are important tools for the shipping infrastructure and to get a better understanding of the effects of climate change on the marine ecosystems. Nemo-Nordic is intended to be a tool for both short-term and long-term simulations and to be used for ocean forecasting as well as process and climatic studies. Here, the scientific and technical choices within Nemo-Nordic are introduced, and the reasons behind the design of the model and its domain and the inclusion of the two seas are explained. The model's ability to represent barotropic and baroclinic dynamics, as well as the vertical structure of the water column, is presented. Biases are shown and discussed. The short-term capabilities of the model are presented, especially its capabilities to represent sea level on an hourly timescale with a high degree of accuracy. We also show that the model can represent longer timescales, with a focus on the major Baltic inflows and the variability in deep-water salinity in the Baltic Sea.

Description: Anthropogenic emissions of greenhouse gases such as CO2 and N2O impinge on the Earth system, which in turn modulates atmospheric greenhouse gas concentrations. The underlying feedback mechanisms are complex and, at times, counterintuitive. So-called Earth system models have recently matured to standard tools tailored to assess these feedback mechanisms in a warming world. Applications for these models range from being targeted at basic process understanding to the assessment of geo-engineering options. A problem endemic to all these applications is the need to estimate poorly known model parameters, specifically for the biogeochemical component, based on observational data (e.g., nutrient fields). In the present study, we illustrate with an Earth system model that through such an approach biases and other model deficiencies in the physical ocean circulation model component can reciprocally compensate for biases in the pelagic biogeochemical model component (and vice versa). We present two model configurations that share a remarkably similar steady state (based on ad hoc measures) when driven by historical boundary conditions, even though they feature substantially different configurations (parameter sets) of ocean mixing and biogeochemical cycling. When projected into the future the similarity between the model responses breaks. Metrics such as changes in total oceanic carbon content and suboxic volume diverge between the model configurations as the Earth warms. Our results reiterate that advancing the understanding of oceanic mixing processes will reduce the uncertainty of future projections of oceanic biogeochemical cycles. Related to the latter, we suggest that an advanced understanding of oceanic biogeochemical cycles can be used for advancements in ocean circulation modules.