ALBERT

Description: Increasing atmospheric CO2 concentrations are potentially affecting marine ecosystems twofold, by warming and acidification. The rising amount of CO2 taken up by the ocean lowers the saturation state of calcium carbonate, complicating the formation of this key biomineral used by many marine organisms to build hard parts like skeletons or shells. Reliable time-series data of seawater pH are needed to evaluate the ongoing change and compare long-term trends and natural variability. For the high-latitude ocean, the region facing the strongest CO2 uptake, such time-series data are so far entirely lacking. Our study provides, to our knowledge, the first reconstruction of seasonal cycle and long-term trend in pH for a high-latitude ocean obtained from 2D images of stable boron isotopes from a coralline alga.

Description: Here we provide optimised vertical eddy diffusivity estimates for the PeECE III and KOSMOS 2013 mesocosm experiment, obtained from a model-based reanalysis. These diffusivities are derived from the observed temperature and salinity profiles that have been published in Schulz et al., 2008. Furthermore, we make our model code available, providing an adjustable tool to simulate vertical mixing in any other pelagic mesocosm. We also provide the interpolated and regridded temperature and salinity profiles of the PeECE III experiment as well as the density profiles which we calculated from the temperature and salinity profiles using the R package seacarb (Lavigne et al., 2011). These data files are required as input to run simulations of the PeECE III experiment with the 1D mesocosm mixing model. The columns of the environmental files (required input files for the model) from left to right are: Experiment year, month, day, Julian day, photosynthetically active radiation (PAR) [W/m^2], temperature [C], salinity [PSU], CO2 concentration [ppm], wind speed [m/s]. The rows list the respective value of each hour of the experiment. Temperature and salinity in this table are hourly interpolated values of the daily measurements published by the PeECE III team (2005). PAR has been calculated from global radiation data of Bergen provided by Olseth et al., 2005. In the temperature, salinity and density files, the rows indicate the depth (0.5 m resolution, the first row is the surface, the last row is the bottom), whereas the columns indicate the experiment time at an hourly resolution.

Description: No records exist to evaluate long-term pH dynamics in high-latitude oceans, which have the greatest probability of rapid acidification from anthropogenic CO2 emissions. We reconstructed both seasonal variability and anthropogenic change in seawater pH and temperature by using laser ablation high-resolution 2D images of stable boron isotopes (δ11B) on a long-lived coralline alga that grew continuously through the 20th century. Analyses focused on four multiannual growth segments. We show a long-term decline of 0.08 ± 0.01 pH units between the end of the 19th and 20th century, which is consistent with atmospheric CO2 records. Additionally, a strong seasonal cycle (∼0.22 pH units) is observed and interpreted as episodic annual pH increases caused by the consumption of CO2 during strong algal (kelp) growth in spring and summer. The rate of acidification intensifies from –0.006 ± 0.007 pH units per decade (between 1920s and 1960s) to –0.019 ± 0.009 pH units per decade (between 1960s and 1990s), and the episodic pH increases show a continuous shift to earlier times of the year throughout the centennial record. This is indicative of ecosystem shifts in shallow water algal productivity in this high-latitude habitat resulting from warming and acidification.

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 incorrect 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: 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 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 inhibits the overall vertical exchange between oxygenated surface waters and oxygen-deprived water at depth. At major upwelling sites, however (e.g. off the southern 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: 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 University of Victoria Earth System Climate Model (UVic ESCM). Transport matrices were extracted for 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: 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: Ocean eddies can both trigger mixing (during their formation and decay) and effectively shield water encompassed from being exchanged with ambient water (throughout their life time). These antagonistic effects of eddies complicate the interpretation of synoptic snapshots as 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. The water mass diversity responds to variations with depth of the circulation strength at the edge of the eddy.