ALBERT

Description: We present bathymetry and hydrological observations collected in the summer of 2014 from two fjord systems in southeastern Greenland with a multibeam sonar system. Our results provide a detailed bathymetric map of the fjord complex around the island of Skjoldungen in Skjoldungen Fjord and the outer part of Timmiarmiut Fjord and show far greater depths compared to the International Bathymetric Chart of the Arctic Ocean. The hydrography collected shows different properties in the fjords with the bottom water masses below 240 m in Timmiarmiut Fjord being 1–2 °C warmer than in the two fjords around Skjoldungen, but data also illustrate the influence of sills on the exchange of deeper water masses within fjords. Moreover, evidence of subglacial discharge in Timmiarmiut Fjord, which is consistent with satellite observations of ice mélange set into motion, adds to our increasing understanding of the distribution of subglacial meltwater.

Description: Ocean acidification resulting from the uptake of anthropogenic carbon dioxide (CO2) by the ocean is considered a major threat to marine ecosystems. Here we examined the effects of ocean acidification on microbial community dynamics in the eastern Baltic Sea during the summer of 2012 when inorganic nitrogen and phosphorus were strongly depleted. Large-volume in situ mesocosms were employed to mimic present, future and far future CO2 scenarios. All six groups of phytoplankton enumerated by flow cytometry ( 〈  20 µm cell diameter) showed distinct trends in net growth and abundance with CO2 enrichment. The picoeukaryotic phytoplankton groups Pico-I and Pico-II displayed enhanced abundances, whilst Pico-III, Synechococcus and the nanoeukaryotic phytoplankton groups were negatively affected by elevated fugacity of CO2 (fCO2). Specifically, the numerically dominant eukaryote, Pico-I, demonstrated increases in gross growth rate with increasing fCO2 sufficient to double its abundance. The dynamics of the prokaryote community closely followed trends in total algal biomass despite differential effects of fCO2 on algal groups. Similarly, viral abundances corresponded to prokaryotic host population dynamics. Viral lysis and grazing were both important in controlling microbial abundances. Overall our results point to a shift, with increasing fCO2, towards a more regenerative system with production dominated by small picoeukaryotic phytoplankton.

Description: Very short-lived substances (VSLS) contribute as source gases significantly to the tropospheric and stratospheric bromine loading. At present, an estimated 25% of stratospheric bromine is of oceanic origin. In this study, we investigate how climate change may impact the ocean- atmosphere flux of brominated VSLS, their atmospheric transport, and chemical transformations and evaluate how these changes will affect stratospheric ozone over the 21st century. Under the assumption of fixed ocean water concentrations and RCP6.0 scenario, we find an increase of the ocean- atmosphere flux of brominated VSLS of about 8-10% by the end of the 21st century compared to present day. A decrease in the tropospheric mixing ratios of VSLS and an increase in the lower stratosphere are attributed to changes in atmospheric chemistry and transport. Our model simulations reveal that this increase is counteracted by a corresponding reduction of inorganic bromine. Therefore the total amount of bromine from VSLS in the stratosphere will not be changed by an increase in upwelling. Part of the increase of VSLS in the tropical lower stratosphere results from an increase in the corresponding tropopause height. As the depletion of stratospheric ozone due to bromine depends also on the availability of chlorine, we find the impact of bromine on stratospheric ozone at the end of the 21st century reduced compared to present day. Thus, these studies highlight the different factors influencing the role of brominated VSLS in a future climate

Description: Reanalysis data sets are widely used to understand atmospheric processes and past variability, and are often used to stand in as "observations" for comparisons with climate model output. Because of the central role of water vapor (WV) and ozone (O3) in climate change, it is important to understand how accurately and consistently these species are represented in existing global reanalyses. In this paper, we present the results of WV and O3 intercomparisons that have been performed as part of the SPARC (Stratosphere–troposphere Processes and their Role in Climate) Reanalysis Intercomparison Project (S-RIP). The comparisons cover a range of timescales and evaluate both inter-reanalysis and observation-reanalysis differences. We also provide a systematic documentation of the treatment of WV and O3 in current reanalyses to aid future research and guide the interpretation of differences amongst reanalysis fields. The assimilation of total column ozone (TCO) observations in newer reanalyses results in realistic representations of TCO in reanalyses except when data coverage is lacking, such as during polar night. The vertical distribution of ozone is also relatively well represented in the stratosphere in reanalyses, particularly given the relatively weak constraints on ozone vertical structure provided by most assimilated observations and the simplistic representations of ozone photochemical processes in most of the reanalysis forecast models. However, significant biases in the vertical distribution of ozone are found in the upper troposphere and lower stratosphere in all reanalyses. In contrast to O3, reanalysis estimates of stratospheric WV are not directly constrained by assimilated data. Observations of atmospheric humidity are typically used only in the troposphere, below a specified vertical level at or near the tropopause. The fidelity of reanalysis stratospheric WV products is therefore mainly dependent on the reanalyses' representation of the physical drivers that influence stratospheric WV, such as temperatures in the tropical tropopause layer, methane oxidation, and the stratospheric overturning circulation. The lack of assimilated observations and known deficiencies in the representation of stratospheric transport in reanalyses result in much poorer agreement amongst observational and reanalysis estimates of stratospheric WV. Hence, stratospheric WV products from the current generation of reanalyses should generally not be used in scientific studies.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: Oxygen optode measurements on floats and gliders suffer from a slow time response and various sources of drift in the calibration coefficients. Based on two dual-O2 Argo floats, we show how to post-correct for the effect of the optode's time response and give an update on optode in situ drift stability and in-air calibration. Both floats are equipped with an unpumped Aanderaa 4330 optode and a pumped Sea-Bird SBE63 optode. Response times for the pumped SBE63 were derived following Bittig et al. (2014) and the same methods were used to correct the time response bias. Using both optodes on each float, the time response regime of the unpumped Aanderaa optode was characterized more accurately than previously possible. Response times for the pumped SBE63 on profiling floats are in the range of 25–40 s, while they are between 60 and 95 s for the unpumped 4330 optode. Our parameterization can be employed to post-correct the slow optode time response on floats and gliders. After correction, both sensors agree to within 2–3 µmol kg−1 (median difference) in the strongest gradients (120 µmol kg−1 change over 8 min or 20 dbar) and better elsewhere. However, time response correction is only possible if measurement times are known, i.e., provided by the platform as well as transmitted and stored with the data. The O2 in-air measurements show a significant in situ optode drift of −0.40 and −0.27 % yr−1 over the available 2 and 3 years of deployment, respectively. Optode in-air measurements are systematically biased high during midday surfacings compared to dusk, dawn, and nighttime. While preference can be given to nighttime surfacings to avoid this in-air calibration bias, we suggest a parameterization of the daytime effect as a function of the Sun's elevation to be able to use all data and to better constrain the result. Taking all effects into account, calibration factors have an uncertainty of 0.1 %. In addition, in-air calibration factors vary by 0.1–0.2 % when using different reanalysis models as a reference. The overall accuracy that can be achieved following the proposed correction routines is better than 1 µmol kg−1.

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: The assessment of the ocean biota's role in climate climate change is often carried out with global biogeochemical ocean models that contain many components, and involve a high level of parametric uncertainty. Examination the models' fit to climatologies of inorganic tracers, after the models have been spun up to steady state, is a common, but computationally expensive procedure to assess model performance and reliability. Using new tools that have become available for global model assessment and calibration in steady state, this paper examines two different model types – a complex seven-component model (MOPS), and a very simple two-component model (RetroMOPS) – for their fit to dissolved quantities. Before comparing the models, a subset of their biogeochemical parameters has been optimised against annual mean nutrients and oxygen. Both model types fit the observations almost equally well. The simple model, which contains only nutrients and dissolved organic phosphorus (DOP), is sensitive to the parameterisation of DOP production and decay. The spatio-temporal decoupling of nitrogen and oxygen, and processes involved in their uptake and release, renders oxygen and nitrate valuable tracers for model calibration. In addition, the non-conservative nature of these tracers (with respect to their upper boundary condition) introduces the global bias as a useful additional constraint on model parameters. Dissolved organic phosphorous at the surface behaves antagonistically to phosphate, and suggests that observations of this tracer – although difficult to measure – may be an important asset for model calibration

Description: The balance in microbial net consumption of nitrogen and phosphorus was investigated in samples collected in two mesotrophic coastal environments: the Baltic Sea (Tvärminne field station) and the North Sea (Espegrend field station). For this, we have refined a bioassay based on the response in alkaline phosphatase activity (APA) over a matrix of combinations in nitrogen and phosphorus additions. This assay not only provides information on which element (N or P) is the primary limiting nutrient, but also gives a quantitative estimate for the excess of the secondary limiting element (P+ or N+, respectively), as well as the ratio of balanced net consumption of added N and P over short timescales (days). As expected for a Baltic Sea late spring–early summer situation, the Tvärminne assays (n =  5) indicated N limitation with an average P+ =  0.30 ± 0.10 µM-P, when incubated for 4 days. For short incubations (1–2 days), the Espegrend assays indicated P limitation, but the shape of the response surface changed with incubation time, resulting in a drift in parameter estimates toward N limitation. Extrapolating back to zero incubation time gave P limitation with N+ ≈  0.9 µM-N. The N : P ratio (molar) of nutrient net consumption varied considerably between investigated locations: from 2.3 ± 0.4 in the Tvärminne samples to 13 ± 5 and 32 ± 3 in two samples from Espegrend. Our assays included samples from mesocosm acidification experiments, but statistically significant effects of ocean acidification were not found by this method.

Description: The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.