ALBERT

Description: Biogeochemical models, capturing the major feedbacks of the pelagic ecosystem of the world ocean, are today often embedded into Earth system models which are increasingly used for decision making regarding climate policies. These models contain poorly constrained parameters (e.g., maximum phytoplankton growth rate), which are typically adjusted until the model shows reasonable behavior. Systematic approaches determine these parameters by minimizing the misfit between the model and observational data. In most common model approaches, however, the underlying functions mimicking the biogeochemical processes are nonlinear and non-convex. Thus, systematic optimization algorithms are likely to get trapped in local minima and might lead to non-optimal results. To judge the quality of an obtained parameter estimate, we propose determining a preferably large lower bound for the global optimum that is relatively easy to obtain and that will help to assess the quality of an optimum, generated by an optimization algorithm. Due to the unavoidable noise component in all observations, such a lower bound is typically larger than zero. We suggest deriving such lower bounds based on typical properties of biogeochemical models (e.g., a limited number of extremes and a bounded time derivative). We illustrate the applicability of the method with two real-world examples. The first example uses real-world observations of the Baltic Sea in a box model setup. The second example considers a three-dimensional coupled ocean circulation model in combination with satellite chlorophyll a.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: In this study ship- and AUV-based multibeam data from the German Mn-nodule license area in the Clarion-Clipperton Zone (CCZ; eastern Pacific) are linked to ground truth data from optical imaging. Photographs obtained by an AUV enable semi-quantitative assessments of nodule coverage at a spatial resolution in the range of meters. Together with high resolution AUV bathymetry this revealed a correlation of small-scale terrain variations (〈 5 m horizontally, 〈 1 m vertically) with nodule abundance. In the presented data set, increased nodule coverage could be correlated with slopes 〉 1.8° and concave terrain. On a more regional scale, factors such as the geological setting (existence of horst and graben structures, sediment thickness, outcropping basement) and influence of bottom currents seem to play an essential role for the spatial variation of nodule abundance and the related hard substrate habitat. AUV imagery was also successfully employed to map the distribution of re-settled sediment following a disturbance and sediment cloud generation during a sampling deployment of an Epibenthic Sledge. Data from before and after the "disturbance" allows a direct assessment of the impact. Automated image processing analyzed the nodule coverage at the seafloor, revealing nodule blanketing by resettling of suspended sediment within 16 hours after the disturbance. The visually detectable impact was spatially limited to a maximum of 100m distance from the disturbance track, downstream of the bottom water current. A correlation with high resolution AUV bathymetry reveals that the blanketing pattern varies in extent by tens of meters, strictly following the bathymetry, even in areas of only slightly undulating seafloor (〈 1 m vertical change). These results highlight the importance of detailed terrain knowledge when engaging in resource assessment studies for nodule abundance estimates and defining minable areas. At the same time, it shows the importance of high resolution mapping for detailed benthic habitat studies that show a heterogeneity at scales of 10 m to 100 m. Terrain knowledge is also needed to determine the scale of the impact by seafloor sediment blanketing during mining-operations.

Description: Climate-driven changes in environmental conditions have significant and complex effects on marine ecosystems. Variability in phytoplankton elements and biochemicals can be important for global ocean biogeochemistry and ecological functions, while there is currently limited understanding on how elements and biochemicals respond to the changing environments in key coccolithophore species such as Emiliania huxleyi. We investigated responses of elemental stoichiometry and fatty acids (FAs) in a strain of E. huxleyi under three temperatures (12, 18 and 24 °C), three N : P supply ratios (molar ratios 10:1, 24:1 and 63:1) and two pCO2 levels (560 and 2400 µatm). Overall, C : N : P stoichiometry showed the most pronounced response to N : P supply ratios, with high ratios of particulate organic carbon vs. particulate organic nitrogen (POC : PON) and low ratios of PON vs. particulate organic phosphorus (PON : POP) in low-N media, and high POC : POP and PON : POP in low-P media. The ratio of particulate inorganic carbon vs. POC (PIC : POC) and polyunsaturated fatty acid proportions strongly responded to temperature and pCO2, both being lower under high pCO2 and higher with warming. We observed synergistic interactions between warming and nutrient deficiency (and high pCO2) on elemental cellular contents and docosahexaenoic acid (DHA) proportion in most cases, indicating the enhanced effect of warming under nutrient deficiency (and high pCO2). Our results suggest differential sensitivity of elements and FAs to the changes in temperature, nutrient availability and pCO2 in E. huxleyi, which is to some extent unique compared to non-calcifying algal classes. Thus, simultaneous changes of elements and FAs should be considered when predicting future roles of E. huxleyi in the biotic-mediated connection between biogeochemical cycles, ecological functions and climate change.

Description: The Eddy Kinetic Energy (EKE) associated with the Subtropical Countercurrent (STCC) in the western subtropical South Pacific is known to exhibit substantial seasonal and decadal variability. Using an eddy-permitting ocean general circulation model, which is able to reproduce the observed, salient features of the seasonal cycles of shear, stratification, baroclinic production and the associated EKE, we investigate the decadal changes of EKE. We show that the STCC region exhibits, uniquely among the subtropical gyres of the world’s oceans, significant, atmospherically forced, decadal EKE variability. The decadal variations are driven by changing vertical shear between the STCC in the upper 300 m and the South Equatorial Current below, predominantly caused by variations in STCC strength associated with a changing meridional density gradient. In the 1970s, an increased meridional density gradient results in EKE twice as large as in later decades in the model. Utilizing sensitivity experiments, decadal variations in the wind field are shown to be the essential driver. Local wind stress curl anomalies associated with the Interdecadal Pacific Oscillation (IPO) lead to up- and downwelling of the thermocline, inducing strengthening or weakening of the STCC and the associated EKE. Additionally, remote wind stress curl anomalies in the eastern subtropical South Pacific, which are not related to the IPO, generate density anomalies that propagate westward as Rossby waves and can account for up to 30–40 % of the density anomalies in the investigated region.

Description: Although mesoscale ocean eddies are ubiquitous in the Southern Ocean, their average regional and seasonal association with phytoplankton has not been quantified systematically yet. To this end, we identify over 100000 mesoscale eddies with diameters of 50km and more in the Southern Ocean and determine the associated phytoplankton biomass anomalies using satellite-based chlorophyll-a (Chl) as a proxy. The mean Chl anomalies, δChl, associated with these eddies, comprising the upper echelon of the oceanic mesoscale, exceed ±10% over wide regions. The structure of these anomalies is largely zonal, with cyclonic, thermocline lifted, eddies having positive anomalies in the subtropical waters north of the Antarctic Circumpolar Current (ACC) and negative anomalies along its main flow path. The pattern is similar, but reversed for anticyclonic, thermocline deepened eddies. The seasonality of δChl is weak in subtropical waters, but pronounced along the ACC, featuring a seasonal sign switch. The spatial structure and seasonality of the mesoscale δChl can be explained largely by lateral advection, especially local eddy-stirring. A prominent exception is the ACC region in winter, where δChl is consistent with a modulation of phytoplankton light exposure caused by an eddy-induced modification of the mixed layer depth. The clear impact of mesoscale eddies on phytoplankton may implicate a downstream effect on Southern Ocean biogeochemical properties, such as mode water nutrient contents.

Description: Benthic foraminifera from Bottsand coastal lagoon, western Baltic Sea, have been studied since the mid-1960s. They were monitored annually in late autumn since 2003 at the terminal ditch of the lagoon. There were 12 different species recognised, of which three have not been recorded during earlier investigations. Dominant species showed strong interannual fluctuations and a steady increase in population densities over the last decade. Elphidium incertum, a stenohaline species of the Baltic deep water fauna, colonised the Bottsand lagoon in 2016, most likely during a period of salinities 〉19 units and water temperatures of 18 °C on average in early autumn. The high salinities probably triggered their germination from a propagule bank in the ditch bottom sediment. The new E. incertum population showed densities higher by an order of magnitude than those of the indigenous species. The latter did not decline, revealing that E. incertum used another food source or occupied a different microhabitat. Elphidium incertum survived transient periods of lower salinities in late autumn 2017, though with reduced abundances, and became a regular faunal constituent at the Bottsand lagoon.

Description: The quantitative reconstruction of past seawater salinity has yet to be achieved and the search for a direct and independent salinity proxy is ongoing. Recent culture and field studies show a significant positive correlation of Na/Ca with salinity in benthic and planktonic foraminiferal calcite. For accurate paleoceanographic reconstructions, consistent and reliable calibrations are necessary, which are still missing. In order to assess the reliability of foraminiferal Na/Ca as a direct proxy for seawater salinity, this study presents electron microprobe Na/Ca data, measured on cultured specimens of Trilobatus sacculifer. The culture experiments were conducted over a wide salinity range of 26 to 45, while temperature was kept constant. To further understand potential controlling factors of Na incorporation, measurements were also performed on foraminifera cultured at various temperatures in the range of 19.5 °C to 29.5 °C under constant salinity conditions. Foraminiferal Na/Ca ratios positively correlate with seawater salinity (Na/Caforam = 0.97 + 0.115 ⋅ Salinity, R = 0.97, p 〈 0.005). Temperature on the other hand exhibits no statistically significant relationship with Na/Ca ratios indicating salinity to be the dominant factor controlling Na incorporation. The culturing results are corroborated by measurements on T. sacculifer from Caribbean and Gulf of Guinea surface sediments. In conclusion, planktonic foraminiferal Na/Ca can be applied as a reliable proxy for reconstructing sea surface salinities, albeit species-specific calibrations might be necessary.

Description: We report here the results of total mercury (HgT) determinations along the 2014 Geotraces Geovide cruise (GA01 transect) in the North Atlantic Ocean (NA) from Lisbon (Portugal) to the coast of Labrador (Canada). HgT concentrations in unfiltered samples (HgTUNF) were log-normally distributed and ranged between 0.16 and 1.54 pmol L−1, with a geometric mean of 0.51 pmol L−1 for the 535 samples analysed. The dissolved fraction (〈 0.45 µm) of HgT (HgTF), determined on 141 samples, averaged 78 % of the HgTUNF for the entire data set, 84 % for open seawaters (below 100 m) and 91 % if the Labrador Sea data are excluded, where the primary production was high (with a winter convection down to 1400 m). HgTUNF concentrations increased eastwards and with depth from Greenland to Europe and from subsurface to bottom waters. The HgTUNF concentrations were similarly low in the subpolar gyre waters ( ∼  0.45 pmol L−1), whereas they exceeded 0.60 pmol L−1 in the subtropical gyre waters. The HgTUNF distribution mirrored that of dissolved oxygen concentration, with highest concentration levels associated with oxygen-depleted zones. The relationship between HgTF and the apparent oxygen utilization confirms the nutrient-like behaviour of Hg in the NA. An extended optimum multiparameter analysis allowed us to characterize HgTUNF concentrations in the different source water types (SWTs) present along the transect. The distribution pattern of HgTUNF, modelled by the mixing of SWTs, show Hg enrichment in Mediterranean waters and North East Atlantic Deep Water and low concentrations in young waters formed in the subpolar gyre and Nordic seas. The change in anthropogenic Hg concentrations in the Labrador Sea Water during its eastward journey suggests a continuous decrease in Hg content in this water mass over the last decades. Calculation of the water transport driven by the Atlantic Meridional Overturning Circulation across the Portugal–Greenland transect indicates northward Hg transport within the upper limb and southward Hg transport within the lower limb, with resulting net northward transport of about 97.2 kmol yr−1.

Description: The chemistry–climate model ECHAM-HAMMOZ contains a detailed representation of tropospheric and stratospheric reactive chemistry and state-of-the-art parameterizations of aerosols using either a modal scheme (M7) or a bin scheme (SALSA). This article describes and evaluates the model version ECHAM6.3-HAM2.3-MOZ1.0 with a focus on the tropospheric gas-phase chemistry. A 10-year model simulation was performed to test the stability of the model and provide data for its evaluation. The comparison to observations concentrates on the year 2008 and includes total column observations of ozone and CO from IASI and OMI, Aura MLS observations of temperature, HNO3, ClO, and O3 for the evaluation of polar stratospheric processes, an ozonesonde climatology, surface ozone observations from the TOAR database, and surface CO data from the Global Atmosphere Watch network. Global budgets of ozone, OH, NOx, aerosols, clouds, and radiation are analyzed and compared to the literature. ECHAM-HAMMOZ performs well in many aspects. However, in the base simulation, lightning NOx emissions are very low, and the impact of the heterogeneous reaction of HNO3 on dust and sea salt aerosol is too strong. Sensitivity simulations with increased lightning NOx or modified heterogeneous chemistry deteriorate the comparison with observations and yield excessively large ozone budget terms and too much OH. We hypothesize that this is an impact of potential issues with tropical convection in the ECHAM model.