ALBERT

Description: PRIMAP2 is the next generation of the PRIMAP climate policy analysis suite.

Description: Characterizing the variability across timescales is important for understanding the underlying dynamics of the Earth system. It remains challenging to do so from palaeoclimate archives since they are more often than not irregular, and traditional methods for producing timescale-dependent estimates of variability, such as the classical periodogram and the multitaper spectrum, generally require regular time sampling. We have compared those traditional methods using interpolation with interpolation-free methods, namely the Lomb–Scargle periodogram and the first-order Haar structure function. The ability of those methods to produce timescale-dependent estimates of variability when applied to irregular data was evaluated in a comparative framework, using surrogate palaeo-proxy data generated with realistic sampling. The metric we chose to compare them is the scaling exponent, i.e. the linear slope in log-transformed coordinates, since it summarizes the behaviour of the variability across timescales. We found that, for scaling estimates in irregular time series, the interpolation-free methods are to be preferred over the methods requiring interpolation as they allow for the utilization of the information from shorter timescales which are particularly affected by the irregularity. In addition, our results suggest that the Haar structure function is the safer choice of interpolation-free method since the Lomb–Scargle periodogram is unreliable when the underlying process generating the time series is not stationary. Given that we cannot know a priori what kind of scaling behaviour is contained in a palaeoclimate time series, and that it is also possible that this changes as a function of timescale, it is a desirable characteristic for the method to handle both stationary and non-stationary cases alike.

Description: The surface mass balance scheme dEBM (diurnal Energy Balance Model) provides a novel interface between the atmosphere and land ice for Earth system modeling, which is based on the energy balance of glaciated surfaces. In contrast to empirical schemes, dEBM accounts for changes in the Earth’s orbit and atmospheric composition. The scheme only requires monthly atmospheric forcing (precipitation, temperature, shortwave and longwave radiation, and cloud cover). It is also computationally inexpensive, which makes it particularly suitable to investigate the ice sheets' response to long-term climate change. After calibration and validation, we analyze the surface mass balance of the Greenland Ice Sheet (GrIS) based on climate simulations representing two warm climate states: a simulation of the mid-Holocene (approximately 6000 years before present) and a climate projection based on an extreme emission scenario which extends to the year 2100. The former period features an intensified summer insolation while the 21st century is characterized by reduced outgoing longwave radiation. Specifically, we investigate whether the temperature–melt relationship, as used in empirical temperature-index methods, remains stable under changing insolation and atmospheric composition. Our results indicate that the temperature–melt relation is sensitive to changes in insolation on orbital timescales but remains mostly invariant under the projected warming climate of the 21st century.

Description: All relevant model outputs to reproduce main analyses and figures in Sakschewski et al. 2021: Variable tree rooting strategies are key to model distribution, productivity and evapotranspiration of tropical evergreen forests, Biogeosciences; https://doi.org/10.5194/bg-2020-97

Description: Model code of LPJmL4.0-VR as used in Sakschewski et al. 2021: Variable tree rooting strategies are key to model distribution, productivity and evapotranspiration of tropical evergreen forests, Biogeosciences; https://doi.org/10.5194/bg-2020-97

Description: This data archive contains model runs and data analysis files to the research article "Impact of declining renewable energy costs on electrification in low emission scenarios" by Gunnar Luderer, Silvia Madeddu, Leon Merfort, Falko Ueckerdt, Michaja Pehl, Robert Pietzcker, Marianna Rottoli, Felix Schreyer, Nico Bauer, Lavinia Baumstark, Christoph Bertram, Alois Dirnaichner, Florian Humpenöder, Antoine Levesque, Alexander Popp, Renato Rodrigues, Jessica Strefler, Elmar Kriegler

Description: Desalination remains in Saudi Arabia the supreme viable alternative to boost the future water supply requirements. The present research is focused on the environmental issue observations and modeling based on dissolved trace metals, hydrographic parameters, and nutrients at the discharging plume and nearby locations of Rabig desalination plant (RDP) and Yanbu desalination plant (YDP) located at the north of Jeddah city (~ 150 and 400 km, respectively). A multi-effect distillation (MED) technology is used in RDP whereas, in YDP, a combination of multistage flashing (MSF), reverse osmosis (RO), and MED technologies is in use for the production of potable waters. The concentration of dissolved Cu, Ni, and Zn was assessed by using cathodic stripping voltammetry. The observed pattern was reflecting an average value of dissolved Cu, Ni, and Zn at RDP that was 8.24, 5.28, and 12.69 nM, respectively. The correspondent changes at YDP were 17.53, 18.06, and 71.19 nM of Cu, Ni, and Zn, respectively. The statistical validation shows a significant positive correlation between metals and temperature with salinity and a significant negative correlation of dissolved oxygen and pH. The background concentration of nearby YDP is comparatively higher and which is ascribed by the negative effects of nearby industries. Detailed investigation of circulation pattern and brine discharge analysis near the Yanbu desalination plant also carried out with the help of 3d hydrodynamical numerical model Delft3d.

Description: Observations of the diversity, distribution and abundance of pelagic fauna are absent for many ocean regions in the Atlantic, but baseline data are required to detect changes in communities as a result of climate change. Gelatinous fauna are increasingly recognized as vital players in oceanic food webs, but sampling these delicate organisms in nets is challenging. Underwater (in situ) observations have provided unprecedented insights into mesopelagic communities in particular for abundance and distribution of gelatinous fauna. In September 2018, we performed horizontal video transects (50-1200 m) using the pelagic in situ observation system during a research cruise in the southern Norwegian Sea. Annotation of the video recordings resulted in 12 abundant and 7 rare taxa. Chaetognaths, the trachymedusa Aglantha digitale and appendicularians were the three most abundant taxa. The high numbers of fishes and crustaceans in the upper 100 m was likely the result of vertical migration. Gelatinous zooplankton included ctenophores (lobate ctenophores, Beroe spp., Euplokamis sp., and an undescribed cydippid) as well as calycophoran and physonect siphonophores. We discuss the distributions of these fauna, some of which represent the first record for the Norwegian Sea.

Description: Ectoenzymatic activity, prokaryotic heterotrophic abundances and production were determined in the Mediterranean Sea. Sampling was carried out in the sub-surface, the deep chlorophyll maximum layer (DCM), the core of the Levantine intermediate waters and in the deeper part of the mesopelagic layers. Michaelis–Menten kinetics were assessed using a large range of concentrations of fluorogenic substrates (0.025 to 50 µM). As a consequence, Km (Michaelis–Menten half-saturation constant) and Vm (maximum hydrolysis velocity) parameters were determined for both low- and high-affinity enzymes for alkaline phosphatase, aminopeptidase (LAP) and β-glucosidase (βGLU). Based on the constant derived from the high-LAP-affinity enzyme (0.025–1 µM substrate concentration range), in situ hydrolysis of N proteins contributed 48 % ± 30 % to the heterotrophic bacterial nitrogen demand within the epipelagic layers and 180 % ± 154 % in the Levantine intermediate waters and the upper part of the mesopelagic layers. The LAP hydrolysis rate was higher than bacterial N demand only within the deeper layer and only when considering the high-affinity enzyme. Based on a 10 % bacterial growth efficiency, the cumulative hydrolysis rates of C proteins and C polysaccharides contributed on average 2.5 % ± 1.3  % to the heterotrophic bacterial carbon demand in the epipelagic layers sampled (sub-surface and DCM). This study clearly reveals potential biases in current and past interpretations of the kinetic parameters for the three enzymes tested based on the fluorogenic-substrate concentration used. In particular, the LAP / βGLU enzymatic ratios and some of the depth-related trends differed between the use of high and low concentrations of fluorogenic substrates.

Description: An oxygen decrease of the intermediate-depth low-oxygen zones (300 to 700 m) is seen in time series for selected tropical areas for the period 1960 to 2008 in the eastern tropical Atlantic, the equatorial Pacific and the eastern tropical Indian Ocean. These nearly 5-decade time series were extended to 68 years by including rare historic data starting in 1950 and more recent data. For the extended time series between 1950 and 2018, the deoxygenation trend for the layer 300 to 700 m is similar to the deoxygenation trend seen in the shorter time series. Additionally, temperature, salinity, and nutrient time series in the upper-ocean layer (50 to 300 m) of these areas were investigated since this layer provides critical pelagic habitat for biological communities. Due to the low amount of data available, the results are often not statistically significant within the 95 % confidence interval but nevertheless indicate trends worth discussing. Generally, oxygen is decreasing in the 50 to 300 m layer, except for an area in the eastern tropical South Atlantic. Nutrients also showed long-term trends in the 50 to 300 m layer in all ocean basins and indicate overlying variability related to climate modes. Nitrate increased in all areas. Phosphate also increased in the Atlantic Ocean and Indian Ocean areas, while it decreased in the two areas of the equatorial Pacific Ocean. Silicate decreased in the Atlantic and Pacific areas but increased in the eastern Indian Ocean. Hence, oxygen and nutrients show trends in the tropical oceans, though nutrients trends are more variable between ocean areas than the oxygen trends; therefore, we conclude that those trends are more dependent on local drivers in addition to a global trend. Different positive and negative trends in temperature, salinity, oxygen and nutrients indicate that oxygen and nutrient trends cannot be completely explained by local warming.

Description: This study presents culture experiments of the cold water species Neogloboquadrina pachyderma (sinistral) and provides new insights into the incorporation of elements in foraminiferal calcite of common and newly established proxies for paleoenvironmental applications (shell Mg/Ca, Sr/Ca and Na/Ca). Specimens were collected from sea ice during the austral winter in the Antarctic Weddell Sea and subsequently cultured at different salinities and a constant temperature. Incorporation of the fluorescent dye calcein showed new chamber formation in the culture at salinities of 30, 31, and 69. Cultured foraminifers at salinities of 46 to 83 only revealed chamber wall thickening, indicated by the fluorescence of the whole shell. Signs of reproduction and the associated gametogenic calcite were not observed in any of the culture experiments. Trace element analyses were performed using an electron microprobe, which revealed increased shell Mg/Ca, Sr/Ca, and Na/Ca values at higher salinities, with Mg/Ca showing the lowest sensitivity to salinity changes. This study enhances the knowledge about unusually high element concentrations in foraminifera shells from high latitudes. Neogloboquadrina pachyderma appears to be able to calcify in the Antarctic sea ice within brine channels, which have low temperatures and exceptionally high salinities due to ongoing sea ice formation.

Description: We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic depth image crossing the 1994 M7.8 Java tsunami earthquake slip area. Seamount subduction is related to the uplift of the overriding plate by lateral shortening and vertical thickening, causing pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate size earthquake. The wrapping of the co-seismic slip contours around the seamount indicates that it diverted rupture propagation, documenting the control of forearc structures on seismic rupture.

Description: This paper uses two subsets of ensemble historical-Nat simulations and pi-Control simulations from CMIP5 as well as observational/reanalysis datasets to investigate responses of the tropical Pacific to the 11-yr solar cycle. A statistically significant 11-yr solar signal is found in the upper-ocean layers above the thermocline and tropospheric circulations. A warming response initially appears in the upper layers of the central equatorial Pacific in the solar maximum years in observations, then increases and shifts into the eastern Pacific at lagged 1–3 yr. Meanwhile, an anomalous updraft arises over the western equatorial Pacific and shifts eastwards in the following years with anomalous subsidence over the Maritime Continent. These lagged responses are confirmed by the historical-Nat simulations, except that the initial signal is located more to the west and all the responses are weaker than the observed. A simplified mixed-layer heat budget analysis based on the historical-Nat simulations suggests that the atmospheric forcing, especially the shortwave radiation, is the major contributor to the initial warming response, and the ocean heat transport effect is responsible for the eastward displacement of the lagged warming responses. In the solar maximum years, the zonal ocean temperature gradient in the western-central Pacific is reduced by the initial warming, and anomalous westerly winds appear over the western equatorial Pacific and extend into the eastern Pacific during the lagged years. These anomalous westerly winds reduce the wind-driven ocean dynamical transport, resulting in the initial warming in the central equatorial Pacific being amplified and the surface warming shifting eastward during the lagged 1–3 yr

Description: During the early part of the last glacial termination (17.2-15 ka) and coincident with a ∼ 35 ppm rise in atmospheric CO2, a sharp 0.3‰-0.4‰ decline in atmospheric δ13CO2 occurred, potentially constraining the key processes that account for the early deglacial CO2 rise. A comparable δ13C decline has also been documented in numerous marine proxy records from surface and thermocline-dwelling planktic foraminifera. The δ13C decline recorded in planktic foraminifera has previously been attributed to the release of respired carbon from the deep ocean that was subsequently transported within the upper ocean to sites where the signal was recorded (and then ultimately transferred to the atmosphere). Benthic δ13C records from the global upper ocean, including a new record presented here from the tropical Pacific, also document this distinct early deglacial δ13C decline. Here we present modeling evidence to show that rather than respired carbon from the deep ocean propagating directly to the upper ocean prior to reaching the atmosphere, the carbon would have first upwelled to the surface in the Southern Ocean where it would have entered the atmosphere. In this way the transmission of isotopically light carbon to the global upper ocean was analogous to the ongoing ocean invasion of fossil fuel CO2. The model results suggest that thermocline waters throughout the ocean and 500-2000m water depths were affected by this atmospheric bridge during the early deglaciation.

Description: The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray OBS network, consisting of 24 broadband ocean bottom seismometers deployed for eight months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ~10 km to ~16 km depth below sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene rift basin. The Basin inversion is suggested to be related to the Africa-Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift structures. Slightly different striking directions of faults in the basin centre compared to faults further east and hence away from the abandoned rift may mimic the counter-clockwise rotation of the Corsica-Sardinia block during ~20–16 Ma. The observed cluster events support the hypothesis of strengthening of crust and uppermost mantle during rifting related extension and thinning of continental crust.

Description: Environmental factors shape the structure and functioning of benthic communities. In coastal zones of the southwestern Baltic Sea, boulder fields represent one of the most productive habitats, supporting diverse benthic communities that provide many ecosystem services. In this study, the influence of the geological characteristics of boulder fields on the biodiversity of associated hard-bottom communities was investigated at two different spatial scales (few kilometers and tens of kilometers). The analyses on overall richness (taxonomic and functional) and community composition revealed how: (i) locally the size of boulders and (ii) regionally site-specific factors like the boulder density distribution and the sediment distribution can act as environmental driving forces. The overall richness of assemblages was shown to increase with increasing surface area of boulders, by up to 60% for species and up to 40% for functional richness. At both investigated scales, differences in compositional variability (β diversity) of the communities were detected. Locally, smallest boulders hosted more variable communities (β diversity up to 2 times higher), while at the regional level, indications of a larger habitat heterogeneity featuring the highest β diversity were observed. This study exemplifies how geological habitat characteristics shape the biodiversity of boulder field communities. The obtained information could be considered in assessment strategies, in order to avoid misclassifications of habitats naturally limited in biodiversity, making a step forward to the desired objective of protecting, conserving, and managing boulder field communities in the study area and at other comparable sites.

Description: Invasive alien species impacts might be mediated by environmental factors such as climatic warming. For invasive predators, multiple predator interactions could also exacerbate or dampen ecological impacts. These effects may be especially pronounced in highly diverse coastal ecosystems that are prone to profound and rapid regime shifts. We examine emergent effects of warming on the strength of intraspecific multiple predator effects from a highly successful invasive gammarid Gammarus tigrinus, using a functional response approach towards larval chironomids (feeding rates under different prey densities). Single predator maximum feeding rates were three-times higher at 24 °C compared to 18 °C overall, with potentially prey destabilising type II functional responses exhibited. However, pairs of gammarids exhibited intraspecific multiple predator effects that were in turn mediated by temperature regime, whereby synergisms were found at the lower temperature (i.e. positive non-trophic interactions) and antagonisms detected at the higher temperature (i.e. negative non-trophic interactions) under high prey densities. Accordingly, warming scenarios may worsen the impact of this invasive alien species, yet implications of temperature change are dependent on predator–predator interactions. Emergent effects between abiotic and biotic factors should be considered in ecological impact predictions across habitat types for invasive alien species.

Description: Background Edema is commonly seen after surgical fixation of ankle fractures. Rest, ice, compression, and elevation (RICE) is an established combination to prevent swelling but hardly able to stimulate lymphatic resorption. Recently, an epicutaneously applied negative pressure suction apparatus (LymphaTouch®) has been introduced to stimulate lymphatic flow. While postoperative recovery, soft tissue, and osseous healing as well as functional outcome are probably linked to the amount of postoperative swelling, estimates on this relative to prevention (RICE) or prevention + stimulated resorption (RICE + ) of fluid are scarce. Methods and analysis This is a single-center, evaluator-blinded randomized pilot trial to investigate postoperative swelling in adults requiring surgical fixation of a closed unilateral ankle fracture. A total of 50 patients will be recruited and randomly assigned to RICE or RICE + prior to surgery. All patients will undergo evaluator-blinded measurements of the ankle volume the day before surgery and subsequently from the evening of the 2nd postoperative day every 24 h until discharge. RICE will be initiated right after surgery and continued until discharge from the hospital in all patients. Additional application of negative pressure therapy (RICE + ) will be initiated on the morning of the 2nd postoperative day and repeated every 24 h until the time of discharge from the hospital. Outcome measures are (i) the relative amount and the time course of the postoperative swelling, (ii) the demand for analgesic therapy (type and amount) together with the perception of pain, (iii) the rate of complications, and (iv) mobility of the ankle joint and the recovery of walking abilities during a 12-weeks follow-up period. Serum and urine samples taken prior to sugery and during postoperative recovery will allow to evaluate the ratio of naturally occurring stable calcium isotopes (δ 44/42 Ca) as a marker of skeletal calcium accrual.

Description: In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D VP and VP/VS models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1–4.2 ML were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500 m in their epicenter and ∼ 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0–20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie–Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).

Description: The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons.

Description: The influence of climate change on the ecological impacts of invasive alien species (IAS) remains understudied, with deoxygenation of aquatic environments often-overlooked as a consequence of climate change. Here, we therefore assessed how oxygen saturation affects the ecological impact of a predatory invasive fish, the Ponto-Caspian round goby (Neogobius melanostomus), relative to a co-occurring endangered European native analogue, the bullhead (Cottus gobio) experiencing decline in the presence of the IAS. In individual trials and mesocosms, we assessed the effect of high, medium and low (90%, 60% and 30%) oxygen saturation on: (1) functional responses (FRs) of the IAS and native, i.e. per capita feeding rates; (2) the impact on prey populations exerted; and (3) how combined impacts of both fishes change over invasion stages (Pre-invasion, Arrival, Replacement, Proliferation). Both species showed Type II potentially destabilising FRs, but at low oxygen saturation, the invader had a significantly higher feeding rate than the native. Relative Impact Potential, combining fish per capita effects and population abundances, revealed that low oxygen saturation exacerbates the high relative impact of the invader. The Relative Total Impact Potential (RTIP), modelling both consumer species’ impacts on prey populations in a system, was consistently higher at low oxygen saturation and especially high during invader Proliferation. In the mesocosm experiment, low oxygen lowered RTIP where both species were present, but again the IAS retained high relative impact during Replacement and Proliferation stages at low oxygen. We also found evidence of multiple predator effects, principally antagonism. We highlight the threat posed to native communities by IAS alongside climate-related stressors, but note that solutions may be available to remedy hypoxia and potentially mitigate impacts across invasion stages.

Description: Direct comparison between paleo oceanic δ13C records and model results facilitates assessing simulated distributions and properties of water masses in the past. To accomplish this, we include a new representation of the stable carbon isotope 13C into the HAMburg Ocean Carbon Cycle model (HAMOCC), the ocean biogeochemical component of the Max Planck Institute Earth System Model (MPI-ESM). 13C is explicitly resolved for all existing oceanic carbon pools. We account for fractionation during air-sea gas exchange and for biological fractionation εp associated with photosynthetic carbon fixation during phytoplankton growth. We examine two εp parameterisations of different complexity: εpPopp varies with surface dissolved CO2 concentration (Popp et al., 1989), while εpLaws additionally depends on local phytoplankton growth rates (Laws et al., 1995). When compared to observations of δ13C in dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to δ13C in particulate organic carbon εpPopp shows a considerably improved performance than εpLaws, because the latter results in a too strong preference for 12C. The model also well reproduces the oceanic 13C Suess effect, i.e. the intrusion of the isotopically light anthropogenic CO2 into the ocean, based on comparison to other existing 13C models and to observation-based oceanic carbon uptake estimates over the industrial period. We further apply the approach of Eide et al. (2017a), who derived the first global oceanic 13C Suess effect estimate based on observations, to our model data that has ample spatial and temporal coverage. With this we are able to analyse in detail the underestimation of 13C Suess effect by this approach as it has been noted by Eide et al. (2017a). Based on our model we find underestimations of 13C Suess effect at 200 m by 0.24 ‰ in the Indian Ocean, 0.21 ‰ in the North Pacific, 0.26 ‰ in the South Pacific, 0.1 ‰ in the North Atlantic and 0.14 ‰ in the South Atlantic. We attribute the major sources of the underestimation to two assumptions in Eide et al. (2017a)'s approach: a spatially-constant preformed component of δ13CDIC in year 1940 and neglecting 13C Suess effect in CFC-12 free water.

Description: A hierarchy of global 1/4° (ORCA025) and Atlantic Ocean 1/20° nested (VIKING20X) ocean/sea-ice models is described. It is shown that the eddy-rich configurations performed in hindcasts of the past 50–60 years under CORE and JRA55-do atmospheric forcings realistically simulate the large-scale horizontal circulation, the distribution of the mesoscale, overflow and convective processes, and the representation of regional current systems in the North and South Atlantic. The representation, and in particular the long-term temporal evolution, of the Atlantic Meridional Overturning Circulation (AMOC) strongly depends on numerical choices for the application of freshwater fluxes. The interannual variability of the AMOC instead is highly correlated among the model experiments and also with observations, including the 2010 minimum observed by RAPID at 26.5° N pointing at a dominant role of the forcing. Regional observations in western boundary current systems at 53° N, 26.5° N and 11° S are explored in respect to their ability to represent the AMOC and to monitor the temporal evolution of the AMOC. Apart from the basin-scale measurements at 26.5° N, it is shown that in particular the outflow of North Atlantic Deepwater at 53° N is a good indicator of the subpolar AMOC trend during the recent decades, if the latter is provided in density coordinates. The good reproduction of observed AMOC and WBC trends in the most reasonable simulations indicate that the eddy-rich VIKING20X is capable in representing realistic forcing-related and ocean-intrinsic trends.

Description: This paper presents a modelling study on the fate of CHBr3 and its product gases in the troposphere within the context of tropical deep convection. A cloud-scale case study was conducted along the west coast of Borneo, where several deep convective systems were triggered on the afternoon and early evening of 19 November 2011. These systems were sampled by the Falcon aircraft during the field campaign of the SHIVA project and analysed using a simulation with the cloud-resolving meteorological model C-CATT-BRAMS at 2x2 km resolution that represents the emissions, transport by large-scale flow, convection, photochemistry, and washout of CHBr3 and its product gases (PGs). We find that simulated CHBr3 mixing ratios and the observed values in the boundary layer and the outflow of the convective systems agree. However, the model underestimates the background CHBr3 mixing ratios in the upper troposphere, which suggests a missing source at the regional scale. An analysis of the simulated chemical speciation of bromine within and around each simulated convective system during the mature convective stage reveals that 〉 85% of the bromine derived from CHBr3 and its PGs is transported vertically to the point of convective detrainment in the form of CHBr3 and that the remaining small fraction is in the form of organic PGs, principally insoluble brominated carbonyls produced from the photo-oxidation of CHBr3. The model simulates that within the boundary layer and free troposphere, the inorganic PGs are only present in soluble forms, i.e. HBr, HOBr, and BrONO2, and, consequently, within the convective clouds, the inorganic PGs are almost entirely removed by wet scavenging. We find that HBr is the most abundant PG in background lower-tropospheric air and that this prevalence of HBr is a result of the relatively low background tropospheric ozone levels at the regional scale. Contrary to a previous study in a different environment, for the conditions in the simulation, the insoluble Br-2 species is hardly formed within the convective systems and therefore plays no significant role in the vertical transport of bromine. This likely results from the relatively small quantities of simulated inorganic bromine involved, the presence of HBr in large excess compared to HOBr and BrO, and the relatively efficient removal of soluble compounds within the convective column.

Description: Earth system and climate modelling involves the simulation of processes on a wide range of scales and within and across various compartments of the Earth system. In practice, component models are often developed independently by different research groups, adapted by others to their special interests and then combined using a dedicated coupling software. This procedure not only leads to a strongly growing number of available versions of model components and coupled setups but also to model- and high-performance computing (HPC)-system-dependent ways of obtaining, configuring, building and operating them. Therefore, implementing these Earth system models (ESMs) can be challenging and extremely time consuming, especially for less experienced modellers or scientists aiming to use different ESMs as in the case of intercomparison projects. To assist researchers and modellers by reducing avoidable complexity, we developed the ESM-Tools software, which provides a standard way for downloading, configuring, compiling, running and monitoring different models on a variety of HPC systems. It should be noted that ESM-Tools is not a coupling software itself but a workflow and infrastructure management tool to provide access to increase usability of already existing components and coupled setups. As coupled ESMs are technically the more challenging tasks, we will focus on coupled setups, always implying that stand-alone models can benefit in the same way. With ESM-Tools, the user is only required to provide a short script consisting of only the experiment-specific definitions, while the software executes all the phases of a simulation in the correct order. The software, which is well documented and easy to install and use, currently supports four ocean models, three atmosphere models, two biogeochemistry models, an ice sheet model, an isostatic adjustment model, a hydrology model and a land-surface model. Compared to previous versions, ESM-Tools has lately been entirely recoded in a high-level programming language (Python) and provides researchers with an even more user-friendly interface for Earth system modelling. ESM-Tools was developed within the framework of the Advanced Earth System Model Capacity project, supported by the Helmholtz Association.

Description: The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

Description: An increasing number of dead zoning (hypoxia) has been reported as a consequence of declining levels of dissolved oxygen in coastal oceans all over the globe. Despite substantial efforts a quantitative description of hypoxia up to a level enabling reliable predictions has not been achieved yet for most regions of societal interest. This does also apply to Eckernförde Bight (EB) situated in the Baltic Sea, Germany. The aim of this study is to dissect underlying mechanisms of hypoxia in EB, to identify key sources of uncertainties, and to explore the potential of existing monitoring programs to predict hypoxia by developing and documenting a workflow that may be applicable to other regions facing similar challenges. Our main tool is an ultra-high spatially resolved general ocean circulation model based on a code framework of proven versatility in that it has been applied to various regional and even global simulations in the past. Our model configuration features a spacial horizontal resolution of 100 m (unprecedented in the underlying framework which is used in both global and regional applications) and includes an elementary representation of the biogeochemical dynamics of dissolved oxygen. In addition, we integrate artificial “clocks” that measure the residence time of the water in EB along with timescales of (surface) ventilation. Our approach relies on an ensemble of hindcast model simulations, covering the period from 2000 to 2018, designed to cover a range of poorly known model parameters for vertical background mixing (diffusivity) and local oxygen consumption within EB. Feed-forward artificial neural networks are used to identify predictors of hypoxia deep in EB based on data at a monitoring site at the entrance of EB. Our results consistently show that the dynamics of low (hypoxic) oxygen concentrations in bottom waters deep inside EB is, to first order, determined by the following antagonistic processes: (1) the inflow of low-oxygenated water from the Kiel Bight (KB) – especially from July to October – and (2) the local ventilation of bottom waters by local (within EB) subduction and vertical mixing. Biogeochemical processes that consume oxygen locally are apparently of minor importance for the development of hypoxic events. Reverse reasoning suggests that subduction and mixing processes in EB contribute, under certain environmental conditions, to the ventilation of the KB by exporting recently ventilated waters enriched in oxygen. A detailed analysis of the 2017 fish-kill incident highlights the interplay between westerly winds importing hypoxia from KB and ventilating easterly winds which subduct oxygenated water.

Description: We describe and test a new model of biological marine silicate cycling, implemented in the Kiel Marine Biogeochemical Model version 3 (KMBM3), embedded in the University of Victoria Earth System Climate Model (UVic ESCM) version 2.9. This new model adds diatoms, which are a key component of the biological carbon pump, to an existing ecosystem model. This new model combines previously published parameterizations of a diatom functional type, opal production and export with a novel, temperature-dependent dissolution scheme. Modelled steady-state biogeochemical rates, carbon and nutrient distributions are similar to those found in previous model versions. The new model performs well against independent ocean biogeochemical indicators and captures the large-scale features of the marine silica cycle to a degree comparable to similar Earth system models. Furthermore, it is computationally efficient, allowing both fully coupled, long-timescale transient simulations and “offline” transport matrix spinups. We assess the fully coupled model against modern ocean observations, the historical record starting from 1960 and a business-as-usual atmospheric CO2 forcing to the year 2300. The model simulates a global decline in net primary production (NPP) of 1.4 % having occurred since the 1960s, with the strongest declines in the tropics, northern midlatitudes and Southern Ocean. The simulated global decline in NPP reverses after the year 2100 (forced by the extended RCP8.5 CO2 concentration scenario), and NPP returns to 98 % of the pre-industrial rate by 2300. This recovery is dominated by increasing primary production in the Southern Ocean, mostly by calcifying phytoplankton. Large increases in calcifying phytoplankton in the Southern Ocean offset a decline in the low latitudes, producing a global net calcite export in 2300 that varies only slightly from pre-industrial rates. Diatom distribution moves southward in our simulations, following the receding Antarctic ice front, but diatoms are outcompeted by calcifiers across most of their pre-industrial Southern Ocean habitat. Global opal export production thus drops to 75 % of its pre-industrial value by 2300. Model nutrients such as phosphate, silicate and nitrate build up along the Southern Ocean particle export pathway, but dissolved iron (for which ocean sources are held constant) increases in the upper ocean. This different behaviour of iron is attributed to a reduction of low-latitude NPP (and consequently, a reduction in both uptake and export and particle, including calcite scavenging), an increase in seawater temperatures (raising the solubility of particulate iron) and stratification that “traps” the iron near the surface. These results are meant to serve as a baseline for sensitivity assessments to be undertaken with this model in the future.

Description: The region around the town Albstadt, SW Germany, was struck by four damaging earthquakes with magnitudes greater than 5 during the last century. These earthquakes occurred along the Albstadt Shear Zone (ASZ), which is characterized by more or less continuous microseismicity. As there are no visible surface ruptures that may be connected to the fault zone, we study its characteristics by its seismicity distribution and faulting pattern. We use the earthquake data of the state earthquake service of Baden-Württemberg from 2011 to 2018 and complement it with additional phase picks beginning in 2016 at the AlpArray and StressTransfer seismic networks in the vicinity of the ASZ. This extended data set is used to determine new minimum 1-D seismic vp and vs velocity models and corresponding station delay times for earthquake relocation. Fault plane solutions are determined for selected events, and the principal stress directions are derived. The minimum 1-D seismic velocity models have a simple and stable layering with increasing velocity with depth in the upper crust. The corresponding station delay times can be explained well by the lateral depth variation of the crystalline basement. The relocated events align about north–south with most of the seismic activity between the towns of Tübingen and Albstadt, east of the 9∘ E meridian. The events can be separated into several subclusters that indicate a segmentation of the ASZ. The majority of the 25 determined fault plane solutions feature an NNE–SSW strike but NNW–SSE-striking fault planes are also observed. The main fault plane associated with the ASZ dips steeply, and the rake indicates mainly sinistral strike-slip, but we also find minor components of normal and reverse faulting. The determined direction of the maximum horizontal stress of 140–149∘ is in good agreement with prior studies. Down to ca. 7–8 km depth SHmax is bigger than SV; below this depth, SV is the main stress component. The direction of SHmax indicates that the stress field in the area of the ASZ is mainly generated by the regional plate driving forces and the Alpine topography.

Description: Information on the feeding habits of species is essential to develop appropriate conservation actions. This study aimed to assess spatial and temporal variation in the diet of the Eurasian otter, Lutra lutra in the Anzali Wetland, through fecal and stable isotope analysis. Seven main prey items were observed in the analysis of 300 fresh spraints. The highest feeding index was observed for fish, followed by snakes. Among the fish species consumed, the index of preponderance of Prussian carp (Carassius gibelio) was the highest followed by pike (Esox lucius) and white bream (Blicca bjoerkna). Results of Shannon diversity index suggest spatial variation of species diversity within and between feeding items (p 〈 0.05); snakes, frogs, and oriental river prawn (Macrobrachium nipponense) showed a seasonal variation. The otter’s trophic level (TL) (3.79) was higher than the TLs of other Anzali Wetland predators, such as pike. Bayesian mixing model showed source proportion contributions of fish 49.5%, reptiles 16.7%, insects 14.8%, crustacean 10.5%, amphibians 4.3%, birds 4.1%, and mollusks 0.1%. When considering only fish species in the Bayesian mixing model, Prussian carp was the main fish prey in the otters’ diet accounting for 47%. Based on the results of this study, the Eurasian otter plays an important role in the ecology of the Anzali Wetland ecosystem even though it preys on exotic species, such as the Prussian carp and the oriental river prawn.

Description: Biogeochemical cycling of carbon (C) and nitrogen (N) in the ocean depends on both the composition and activity of underlying biological communities and on abiotic factors. The Southern Ocean is encircled by a series of strong currents and fronts, providing a barrier to microbial dispersion into adjacent oligotrophic gyres. Our study region straddles the boundary between the nutrient-rich Southern Ocean and the adjacent oligotrophic gyre of the South Indian Ocean, providing an ideal region to study changes in microbial productivity. Here, we measured the impact of C- and N- uptake on microbial community diversity, contextualized by hydrographic factors and local physico-chemical conditions across the Southern Ocean and South Indian Ocean. We observed that contrasting physico-chemical characteristics led to unique microbial diversity patterns, with significant correlations between microbial alpha diversity and primary productivity (PP). However, we detected no link between specific PP (PP normalized by chlorophyll a concentration) and microbial alpha and beta diversity. Prokaryotic alpha and beta diversity were correlated with biological N2 fixation, itself a prokaryotic process, and we detected measurable N2 fixation to 60° S. While regional water masses have distinct microbial genetic fingerprints in both the eukaryotic and prokaryotic fractions, PP and N2 fixation vary more gradually and regionally. This suggests that microbial phylogenetic diversity is more strongly bounded by physical oceanographic features, while microbial activity responds more to chemical factors. We conclude that concomitant assessments of microbial diversity and activity is central in understanding the dynamics and complex responses of microorganisms to a changing ocean environment.

Description: Cristamonadea is a large class of parabasalian protists that reside in the hindguts of wood-feeding insects, where they play an essential role in the digestion of lignocellulose. This group of symbionts boasts an impressive array of complex morphological characteristics, many of which have evolved multiple times independently. However, their diversity is understudied and molecular data remain scarce. Here we describe seven new species of cristamonad symbionts from Comatermes, Calcaritermes, and Rugitermes termites from Peru and Ecuador. To classify these new species, we examined cells by light and scanning electron microscopy, sequenced the symbiont small subunit ribosomal RNA (rRNA) genes, and carried out barcoding of the mitochondrial large subunit rRNA gene of the hosts to confirm host identification. Based on these data, five of the symbionts characterized here represent new species within described genera: Devescovina sapara n. sp., Devescovina aymara n. sp., Macrotrichomonas ashaninka n. sp., Macrotrichomonas secoya n. sp., and Macrotrichomonas yanesha n. sp. Additionally, two symbionts with overall morphological characteristics similar to the poorly-studied and probably polyphyletic ‘joeniid’ Parabasalia are classified in a new genus Runanympha n. gen.: Runanympha illapa n. sp., and Runanympha pacha n. sp.

Description: We present an ensemble of last glacial inception (LGI) simulations for the Northern Hemisphere that captures a significant fraction of inferred ice volume changes within proxy uncertainties. This ensemble was performed with LCice 1.0, a coupled ice sheet and climate model, varying parameters of both climate and ice sheet components, as well as the coupling between them. Certain characteristics of the spatiotemporal pattern of ice growth and subsequent retreat in both North America (NA) and Eurasia (EA) are sensitive to parameter changes while others are not. We find that the initial inception of ice over NA and EA is best characterized by the nucleation of ice at high-latitude and high-elevation sites. Subsequent spreading and merger along with large-scale conversion of snowfields dominate in different sectors. The latter plays an important role in the merging of eastern and western ice regions in NA. The inception peak ice volume in the ensemble occurs approximately at 111 ka and therefore lags the summer 60∘ N insolation minimum by more than 3 kyr. Ice volumes consistently peak earlier over EA than NA. The inception peak in North America is characterized by a merged Laurentide and Cordilleran ice sheet, with the Davis Strait covered in ice in ∼80 % of simulations. Ice also bridges Greenland and Iceland in all runs by 114 ka and therefore blocks the Denmark Strait. This latter feature would thereby divert the East Greenland Current and Denmark Strait overflow with a potentially significant impact on ocean circulation. The Eurasian ice sheet at its inception peak varies across ensemble runs between a continuous ice sheet and multiple smaller ice caps. In both continents, the colder high latitudes (i.e. Ellesmere and Svalbard) tend to grow ice through the entire simulation (to 102 ka), while lower latitudes lose ice after ∼110 ka. We find temperature decreases over the initial phases of the inception lead to the expansion of NA ice sheet area and that subsequent precipitation increases contribute to its thickening. EA ice sheet area also expands with decreasing temperatures, but sea ice limits any increases in precipitation, leading to an earlier retreat away from the EA maximum ice sheet volume. We also examine the extent to which the capture of both LGI ice growth and retreat constrains the coupled ice–climate model sensitivity to changing atmospheric pCO2. The 55-member sub-ensemble that meets our criteria for “acceptable” ice growth and retreat has an equilibrium climate sensitivity lower bound that is 0.3 ∘C higher than that of the full ensemble. This suggests some potential value of fully coupled ice–climate modelling of the last glacial inception to constrain future climate change.

Description: Open-ocean oxygen minimum zones (OMZs) occur in regions with high biological productivity and weak ventilation. They restrict marine habitats and alter biogeochemical cycles. Global models generally show a large model–data misfit with regard to oxygen. Reliable statements about the future development of OMZs and the quantification of their interaction with climate change are currently not possible. One of the most intense OMZs worldwide is located in the Arabian Sea (AS). We give an overview of the main model deficiencies with a detailed comparison of the historical state of 10 climate models from the 5th Coupled Model Intercomparison Project (CMIP5) that present our present-day understanding of physical and biogeochemical processes. Most of the models show a general underestimation of the OMZ volume in the AS compared to observations that is caused by an overly shallow layer of oxygen-poor water in the models. The deviation of oxygen values in the deep AS is the result of oxygen levels that are too high simulated in the Southern Ocean formation regions of Indian Ocean Deep Water in the models compared to observations and uncertainties in the deepwater mass transport from the Southern Ocean northward into the AS. Differences in simulated water mass properties and ventilation rates of Red Sea Water and Persian Gulf Water cause different mixing in the AS and thus influence the intensity of the OMZ. These differences in ventilation rates also point towards variations in the parameterizations of the overflow from the marginal seas among the models. The results of this study are intended to foster future model improvements regarding the OMZ in the AS.

Description: Marine particulate organic carbon-13 stable isotope ratios (δ13CPOC) provide insights in understanding carbon cycling through the atmosphere, ocean, and biosphere. They have been used to trace the input of anthropogenic carbon in the marine ecosystem due to the distinct isotopically light signature of anthropogenic emissions. However, δ13CPOC is also significantly altered during photosynthesis by phytoplankton, which complicates its interpretation. For such purposes, robust spatio-temporal coverage of δ13CP OC observations is essential. We collected all such available data sets, merged and homogenized them to provide the largest available marine δ13CPOC data set (Verwega et al., 2021). The data set consists of 4732 data points covering all major ocean basins beginning in the 1960s. We describe the compiled raw data, compare different observational methods, and provide key insights in the temporal and spatial distribution that is consistent with previously observed patterns. The main different sample collection methods (bottle, intake, net, trap) are generally consistent with each other when comparing within regions. An analysis of 1990s mean δ13CP OC values in an meridional section accross the Atlantic Ocean shows relatively high values (≥ −22 ‰) in the low latitudes (〈 30°) trending towards lower values in the Arctic Ocean (∼ −24 ‰) and Southern Ocean (≤ −28 ‰). The temporal trend since the 1960s shows a decrease of mean δ13CPOC by more than 3 ‰ in all basins except for the Southern Ocean which shows a weaker trend but contains relatively poor multi-decadal coverage.

Description: The Laptev Sea and the East Siberian Sea are remote areas of the Arctic region where detailed data on phytoplankton composition and spatial distribution remain limited. In the context of the ongoing environmental changes (increasing warming and ice melting) and prospective exploration activities (oil and gas production) on the Arctic shelves, understanding of the seasonal and interannual phytoplankton community dynamics is of critical importance. Our study provides new specifying data on species composition of phytoplankton over the vast area of the Laptev Sea shelf and the East Siberian Sea shelf. We found that the outer shelf of the Laptev and East Siberian seas was characterized by typical late spring diatom species (Chaetoceros furcellatus, Chaetoceros diadema, Chaetoceros debilis, Chaetoceros constrictus). On the inner shelf of the Laptev Sea, which is strongly affected by the Lena River water masses, the phytoplankton were characterized by the transition from the summer to an autumn stage of development. Local algal assemblages were composed by mixo- and heterotrophic dinoflagellates (Dinophysis and Protoperidinium genera) together with marine and brackish water-marine diatoms (Thalassiosira hyperborea, Thalassiosira baltica, Thalassiosira gravida, Thalassiosira nordenskioeldii) accompanied by sporadically occurring freshwater riverine planktonic diatom species (Aulacoseira granulata, Aulacoseira italica, Asterionella formosa). These variations in species composition over the Laptev Sea shelf were attributed to differences in the hydrography, marine chemical conditions, and the sea-ice regime.

Description: Permafrost degradation in the catchment of major Siberian rivers, combined with higher precipitation in a warming climate, could increase the flux of terrestrially derived dissolved organic matter (tDOM) into the Arctic Ocean (AO). Each year, ∼ 7.9 Tg of dissolved organic carbon (DOC) is discharged into the AO via the three largest rivers that flow into the Laptev Sea (LS) and East Siberian Sea (ESS). A significant proportion of this tDOM-rich river water undergoes at least one freeze–melt cycle in the land-fast ice that forms along the coast of the Laptev and East Siberian seas in winter. To better understand how growth and melting of land-fast ice affect dissolved organic matter (DOM) dynamics in the LS and ESS, we determined DOC concentrations and the optical properties of coloured dissolved organic matter (CDOM) in sea ice, river water and seawater. The data set, covering different seasons over a 9-year period (2010–2019), was complemented by oceanographic measurements (T, S) and determination of the oxygen isotope composition of the seawater. Although removal of tDOM cannot be ruled out, our study suggests that conservative mixing of high-tDOM river water and sea-ice meltwater with low-tDOM seawater is the major factor controlling the surface distribution of tDOM in the LS and ESS. A case study based on data from winter 2012 and spring 2014 reveals that the mixing of about 273 km3 of low-tDOM land-fast-ice meltwater (containing ∼ 0.3 Tg DOC) with more than 200 km3 of high-tDOM Lena River water discharged during the spring freshet (∼ 2.8 Tg DOC yr−1) plays a dominant role in this respect. The mixing of the two low-salinity surface water masses is possible because the meltwater and the river water of the spring freshet flow into the southeastern LS at the same time every year (May–July). In addition, budget calculations indicate that in the course of the growth of land-fast ice in the southeastern LS, ∼ 1.2 Tg DOC yr−1 (± 0.54 Tg) can be expelled from the growing ice in winter, together with brines. These DOC-rich brines can then be transported across the shelves into the Arctic halocline and the Transpolar Drift Current flowing from the Siberian Shelf towards Greenland. The study of dissolved organic matter dynamics in the AO is important not only to decipher the Arctic carbon cycle but also because CDOM regulates physical processes such as radiative forcing in the upper ocean, which has important effects on sea surface temperature, water column stratification, biological productivity and UV penetration.

Description: The organic mass fraction from sea spray aerosol (SSA) is currently a subject of intense research. The majority of this research is dedicated to measurements in ambient air. However a number of studies have recently started to focus on nascent sea spray aerosol. This work presents measurements collected during a 5-week cruise in May and June 2017 in the central and western Mediterranean Sea, an oligotrophic marine region with low phytoplankton biomass. Surface seawater was continuously pumped into a bubble-bursting apparatus to generate nascent sea spray aerosol. Size distributions were measured with a differential mobility particle sizer (DMPS). Chemical characterization of the submicron aerosol was performed with a time-of-flight aerosol chemical speciation monitor (ToF-ACSM) operating with 10 min time resolution and with filter-based chemical analysis on a daily basis. Using positive matrix factorization analysis, the ToF-ACSM non-refractory organic matter (OMNR) was separated into four different organic aerosol types, identified as primary OA (POANR), oxidized OA (OOANR), methanesulfonic acid type OA (MSA-OANR), and mixed OA (MOANR). In parallel, surface seawater biogeochemical properties were monitored providing information on phytoplankton cell abundance and seawater particulate organic carbon (1 h time resolution) and seawater surface microlayer (SML) dissolved organic carbon (DOC) (on a daily basis). Statistically robust correlations (for n〉500) were found between MOANR and nanophytoplankton cell abundance, as well as between POANR, OOANR, and particulate organic carbon (POC). Parameterizations of the contributions of different types of organics to the submicron nascent sea spray aerosol are proposed as a function of the seawater biogeochemical properties for use in models.

Description: The science guiding the EUREC4A campaign and its measurements is presented. EUREC4A comprised roughly 5 weeks of measurements in the downstream winter trades of the North Atlantic – eastward and southeastward of Barbados. Through its ability to characterize processes operating across a wide range of scales, EUREC4A marked a turning point in our ability to observationally study factors influencing clouds in the trades, how they will respond to warming, and their link to other components of the earth system, such as upper-ocean processes or the life cycle of particulate matter. This characterization was made possible by thousands (2500) of sondes distributed to measure circulations on meso- (200 km) and larger (500 km) scales, roughly 400 h of flight time by four heavily instrumented research aircraft; four global-class research vessels; an advanced ground-based cloud observatory; scores of autonomous observing platforms operating in the upper ocean (nearly 10 000 profiles), lower atmosphere (continuous profiling), and along the air–sea interface; a network of water stable isotopologue measurements; targeted tasking of satellite remote sensing; and modeling with a new generation of weather and climate models. In addition to providing an outline of the novel measurements and their composition into a unified and coordinated campaign, the six distinct scientific facets that EUREC4A explored – from North Brazil Current rings to turbulence-induced clustering of cloud droplets and its influence on warm-rain formation – are presented along with an overview of EUREC4A's outreach activities, environmental impact, and guidelines for scientific practice. Track data for all platforms are standardized and accessible at https://doi.org/10.25326/165 (Stevens, 2021), and a film documenting the campaign is provided as a video supplement.

Description: Supratidal sands are vitally important for coastal defence in the German Wadden Sea. They are less affected by human activities than other areas as they are located far off the mainland shore, touristical and commercial activities are generally prohibited. Therefore, supratidal sands are of high ecological interest. Nevertheless, the faunal inventory and distribution pattern of microorganisms on these sands were studied very little. The composition of living and dead foraminiferal assemblages was therefore investigated along a transect from the supratidal sand Japsand up to Hallig Hooge. Both assemblages were dominated by calcareous foraminifera of which Ammonia batava was the most abundant species. Elphidium selseyense and Elphidium williamsoni were also common in the living assemblage, but Elphidium williamsoni was comparably rare in the dead assemblage. The high proportions of Ammonia batava and Elphidium selseyense in the living assemblage arose from the reproduction season that differed between species. While Ammonia batava and Elphidium selseyense just finished their reproductive cycles, Elphidium williamsoni was just about to start. This was also confirmed by the size distribution patterns of the different species. The dead assemblage revealed 20 species that were not found in the living assemblage of which some were reworked from older sediments (e.g., Bucella frigida) and some were transported via tidal currents from other areas in the North Sea (e.g., Jadammina macrescens). The living foraminiferal faunas depicted close linkages between the open North Sea and the mainland. Key species revealing exchange between distant populations were Haynesina germanica, Ammonia batava and different Elphidium species. All these species share an opportunistic behaviour and are able to inhabit a variety of different environments; hence, they well may cope with changing environmental conditions. The benthic foraminiferal association from Japsand revealed that transport mechanisms via tides and currents play a major ecological role and strongly influence the faunal composition at this site.

Description: Nutrient depletion in the ecosystem mainly occurred via the base cations leach from the ecosystems. Isotope technology can provide new information on the nutrients’ behavior (cycle) in ecosystems. To better understand the strontium (Sr) biogeochemical cycle in the forest, radiogenic and stable Sr isotopes were determined in various environmental samples (rainwater, groundwater, calcareous soil, vegetation, and bedrocks) collected from a karst forest ecosystem, Southwest China. The Sr supplies are supported via the deposition of atmospheric materials and the weathering of carbonate rock in this area. These supplied sources observably influence the Sr isotope compositions (both stable and radiogenic isotope). Across the ecosystem, the varied ranges of 87Sr/86Sr ratios increased in the order: bedrock (0.708) ~ groundwater (0.708–0.709) ~ rainwater (0.709–0.711) ~ vegetation (0.709–0.710) 〈 surface calcareous soil (0.713–0.719) 〈 deep calcareous soil (0.714–0.724). Both rainwater and groundwater present similar 87Sr/86Sr ratios to the water-draining carbonate-dominated terrain. The 87Sr/86Sr ratios of vegetation are similar to that of rainwater but slightly heavier than groundwater. In contrast, the varied ranges of δ88/86Sr values followed a sequence of rainwater (0.08–0.21‰) 〈 groundwater (0.11–0.27‰) 〈 deep calcareous soil (0.02–0.29‰) 〈 tree leaves (0.12–0.17‰) 〈 surface calcareous soil (0.14–0.29‰) 〈 bedrock (0.40‰) 〈 shrub leaves (0.50‰). The δ88/86Sr values of tree leaves are between surface calcareous soil and deep calcareous soil, and the δ88/86Sr values of shrub leaves are higher than all of the compartments in the forest ecosystem, which indicated vegetation uptake might control the stable strontium fractionation. In combination with the radiogenic and stable Sr isotopes, our investigation presents that the 87Sr/86Sr provides sources’ information, whereas the δ88/86Sr records biogeochemical cycle information.

Description: The Last Glacial Maximum (LGM, ∼ 21 000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models has been used to generate LGM simulations as part of the Paleoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here, we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4, most of which are PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3, most of which are PMIP3-CMIP5). We show that the global averages of the PMIP4 simulations span a larger range in terms of mean annual surface air temperature and mean annual precipitation compared to the PMIP3-CMIP5 simulations, with some PMIP4 simulations reaching a globally colder and drier state. However, the multi-model global cooling average is similar for the PMIP4 and PMIP3 ensembles, while the multi-model PMIP4 mean annual precipitation average is drier than the PMIP3 one. There are important differences in both atmospheric and oceanic circulations between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large. Therefore, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land–sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the paleoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. These results point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance.

Description: We present a transient simulation of global vegetation and climate patterns of the mid- and late Holocene using the MPI-ESM (Max Planck Institute for Meteorology Earth System Model) at T63 resolution. The simulated vegetation trend is discussed in the context of the simulated Holocene climate change. Our model captures the main trends found in reconstructions. Most prominent are the southward retreat of the northern treeline that is combined with the strong decrease of forest in the high northern latitudes during the Holocene and the vast increase of the Saharan desert, embedded in a general decrease in precipitation and vegetation in the Northern Hemisphere monsoon margin regions. The Southern Hemisphere experiences weaker changes in total vegetation cover during the last 8000 years. However, the monsoon-related increase in precipitation and the insolation-induced cooling of the winter climate lead to shifts in the vegetation composition, mainly between the woody plant functional types (PFTs). The large-scale global patterns of vegetation almost linearly follow the subtle, approximately linear, orbital forcing. In some regions, however, non-linear, more rapid changes in vegetation are found in the simulation. The most striking region is the Sahel–Sahara domain with rapid vegetation transitions to a rather desertic state, despite a gradual insolation forcing. Rapid shifts in the simulated vegetation also occur in the high northern latitudes, in South Asia and in the monsoon margins of the Southern Hemisphere. These rapid changes are mainly triggered by changes in the winter temperatures, which go into, or move out of, the bioclimatic tolerance range of individual PFTs. The dynamics of the transitions are determined by dynamics of the net primary production (NPP) and the competition between PFTs. These changes mainly occur on timescales of centuries. More rapid changes in PFTs that occur within a few decades are mainly associated with the timescales of mortality and the bioclimatic thresholds implicit in the dynamic vegetation model, which have to be interpreted with caution. Most of the simulated Holocene vegetation changes outside the high northern latitudes are associated with modifications in the intensity of the global summer monsoon dynamics that also affect the circulation in the extra tropics via teleconnections. Based on our simulations, we thus identify the global monsoons as the key player in Holocene climate and vegetation change.

Description: This work introduces Waterburya agarophytonicola Bonthond and Shalygin gen. nov., sp. nov, a baeocyte producing cyanobacterium that was isolated from the rhodophyte Agarophyton vermiculophyllum (Ohmi) Gurgel et al., an invasive seaweed that has spread across the northern hemisphere. The new species genome reveals a diverse repertoire of chemotaxis and adhesion related genes, including genes coding for type IV pili assembly proteins and a high number of genes coding for filamentous hemagglutinin family (FHA) proteins. Among a genetic basis for the synthesis of siderophores, carotenoids and numerous vitamins, W. agarophytonicola is potentially capable of producing cobalamin (vitamin B-12), for which A. vermiculophyllum is an auxotroph. With a taxonomic description of the genus and species and a draft genome, this study provides as a basis for future research, to uncover the nature of this geographically independent association between seaweed and cyanobiont.

Description: We review the literature on the relationship between water temperature and size of freshwater phytoplankton, to examine the hypothesis that freshwater phytoplankton, like marine phytoplankton and many other groups of organisms, conform to Bergmann's Rule and become smaller with warming. We provide both experimental and field evidence in support of the above hypothesis, much of this evidence was hidden in studies focused on other issues, but presenting temperature and phytoplankton size data. Freshwater phytoplankton size shrinks with increasing temperature at both the species level (by cells or colonies becoming smaller) and at the community level (shift to smaller species). Exceptions to the Rule do occur but in most cases those exceptions can be explained by indirect effects of temperature on phytoplankton size, via processes such as grazing or nutrient availability. With global warming, freshwater phytoplankton are likely to be of smaller size. This article is dedicated to Colin S. Reynolds, who has had a leading role in our personal education and understanding of phytoplankton ecology.

Description: In order to investigate the impact of spatial resolution on the discrepancy between simulated δ18O and observed δ18O in Greenland ice cores, regional climate simulations are performed with the isotope-enabled regional climate model (RCM) COSMO_iso. For this purpose, isotope-enabled general circulation model (GCM) simulations with the ECHAM5-wiso general circulation model (GCM) under present-day conditions and the MPI-ESM-wiso GCM under mid-Holocene conditions are dynamically downscaled with COSMO_iso for the Arctic region. The capability of COSMO_iso to reproduce observed isotopic ratios in Greenland ice cores for these two periods is investigated by comparing the simulation results to measured δ18O ratios from snow pit samples, Global Network of Isotopes in Precipitation (GNIP) stations and ice cores. To our knowledge, this is the first time that a mid-Holocene isotope-enabled RCM simulation is performed for the Arctic region. Under present-day conditions, a dynamical downscaling of ECHAM5-wiso (1.1∘×1.1∘ ) with COSMO_iso to a spatial resolution of 50 km improves the agreement with the measured δ18O ratios for 14 of 19 observational data sets. A further increase in the spatial resolution to 7 km does not yield substantial improvements except for the coastal areas with its complex terrain. For the mid-Holocene, a fully coupled MPI-ESM-wiso time slice simulation is downscaled with COSMO_iso to a spatial resolution of 50 km. In the mid-Holocene, MPI-ESM-wiso already agrees well with observations in Greenland and a downscaling with COSMO_iso does not further improve the model–data agreement. Despite this lack of improvement in model biases, the study shows that in both periods, observed δ18O values at measurement sites constitute isotope ratios which are mainly within the subgrid-scale variability of the global ECHAM5-wiso and MPI-ESM-wiso simulation results. The correct δ18O ratios are consequently not resolved in the GCM simulation results and need to be extracted by a refinement with an RCM. In this context, the RCM simulations provide a spatial δ18O distribution by which the effects of local uncertainties can be taken into account in the comparison between point measurements and model outputs. Thus, an isotope-enabled GCM–RCM model chain with realistically implemented fractionating processes constitutes a useful supplement to reconstruct regional paleo-climate conditions during the mid-Holocene in Greenland. Such model chains might also be applied to reveal the full potential of GCMs in other regions and climate periods, in which large deviations relative to observed isotope ratios are simulated.

Description: The ocean and inland waters are two separate regimes, with concentrations in greenhouse gases differing on orders of magnitude between them. Together, they create the land–ocean aquatic continuum (LOAC), which comprises itself largely of areas with little to no data with regards to understanding the global carbon system. Reasons for this include remote and inaccessible sample locations, often tedious methods that require collection of water samples and subsequent analysis in the lab, and the complex interplay of biological, physical and chemical processes. This has led to large inconsistencies, increasing errors and has inevitably lead to potentially false upscaling. A set-up of multiple pre-existing oceanographic sensors allowing for highly detailed and accurate measurements was successfully deployed in oceanic to remote inland regions over extreme concentration ranges. The set-up consists of four sensors simultaneously measuring pCO2, pCH4 (both flow-through, membrane-based non-dispersive infrared (NDIR) or tunable diode laser absorption spectroscopy (TDLAS) sensors), O2 and a thermosalinograph at high resolution from the same water source. The flexibility of the system allowed for deployment from freshwater to open ocean conditions on varying vessel sizes, where we managed to capture day–night cycles, repeat transects and also delineate small-scale variability. Our work demonstrates the need for increased spatiotemporal monitoring and shows a way of homogenizing methods and data streams in the ocean and limnic realms.

Description: This study evaluates the potential usefulness of the halogenated compounds HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14, and PFC-116 as oceanographic transient tracers to better constrain ocean ventilation processes. We do this mainly in terms of four aspects of the characteristics of the potential tracers: input function (including atmospheric history and historical surface saturation), seawater solubility, feasibility of measurement, and stability in seawater; of these, atmospheric history and seawater solubility have been investigated in previous work. For the latter two aspects, we collected seawater samples and modified an established analytical technique for the Medusa–Aqua system to simultaneously measure these compounds. HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, and HFC-125 have been measured in depth profiles in the Mediterranean Sea for the first time and for reproducibility in the Baltic Sea. We found that the historical surface saturation of halogenated transient tracers in the Mediterranean Sea is estimated to have been nearly constant at 94 % based on historical data. Of the investigated compounds, HCFC-142b, HCFC-141b, and HFC-134a are found to currently be the most promising transient tracers in the ocean. The compounds that have the greatest potential as future tracers are PFC-14 and PFC-116, mainly hampered by the low solubility in seawater that creates challenging analytical conditions, i.e., low concentrations. HCFC-22 is found to be likely unstable in warm seawater, which compromises the potential as an oceanic transient tracer, although it is possibly useful in colder water. For the compounds HFC-125 and HFC-23, we were not able to fully evaluate their potential as tracers due to inconclusive results, especially regarding their solubility and stability in seawater, but also with regard to potential analytical challenges. On the other hand, HFC-125, HFC-23, and HCFC-22 might not need to be considered because there are alternative tracers with similar input histories that are better suited as transient tracers.

Description: Inferences about climate states and climate variability of the Holocene and the deglaciation rely on sparse paleo-observational proxy data. Combining these proxies with output from climate simulations is a means for increasing the understanding of the climate throughout the last tens of thousands of years. The analogue method is one approach to do this. The method takes a number of sparse proxy records and then searches within a pool of more complete information (e.g., model simulations) for analogues according to a similarity criterion. The analogue method is non-linear and allows considering the spatial covariance among proxy records. Beyond the last two millennia, we have to rely on proxies that are not only sparse in space but also irregular in time and with considerably uncertain dating. This poses additional challenges for the analogue method, which have seldom been addressed previously. The method has to address the uncertainty of the proxy-inferred variables as well as the uncertain dating. It has to cope with the irregular and non-synchronous sampling of different proxies. Here, we describe an implementation of the analogue method including a specific way of addressing these obstacles. We include the uncertainty in our proxy estimates by using “ellipses of tolerance” for tuples of individual proxy values and dates. These ellipses are central to our approach. They describe a region in the plane spanned by proxy dimension and time dimension for which a model analogue is considered to be acceptable. They allow us to consider the dating as well as the data uncertainty. They therefore form the basic criterion for selecting valid analogues. We discuss the benefits and limitations of this approach. The results highlight the potential of the analogue method to reconstruct the climate from the deglaciation up to the late Holocene. However, in the present case, the reconstructions show little variability of their central estimates but large uncertainty ranges. The reconstruction by analogue provides not only a regional average record but also allows assessing the spatial climate field compliant with the used proxy predictors. These fields reveal that uncertainties are also locally large. Our results emphasize the ambiguity of reconstructions from spatially sparse and temporally uncertain, irregularly sampled proxies.

Description: We report the distributions and stoichiometry ofdissolved zinc (dZn) and cobalt (dCo) in sub-tropical andsub-Antarctic waters of the south-eastern Atlantic Oceanduring austral spring 2010 and summer 2011/2012. In sub-tropical surface waters, mixed-layer dZn and dCo con-centrations during early spring were 1.60±2.58 nM and30±11 pM, respectively, compared with summer values of0.14±0.08 nM and 24±6 pM. The elevated spring dZn con-centrations resulted from an apparent offshore transport ofelevated dZn at depths between 20–55 m, derived from theAgulhas Bank. In contrast, open-ocean sub-Antarctic surfacewaters displayed largely consistent inter-seasonal mixed-layer dZn and dCo concentrations of 0.10±0.07 nM and11±5 pM, respectively. Trace metal stoichiometry, calcu-lated from concentration inventories, suggests a greater over-all removal for dZn relative to dCo in the upper water columnof the south-eastern Atlantic, with inter-seasonally decreas-ing dZn/dCo inventory ratios of 19–5 and 13–7 mol mol−1for sub-tropical surface water and sub-Antarctic surface wa-ter, respectively. In this paper, we investigate how the sea-sonal influences of external input and phytoplankton succes-sion may relate to the distribution of dZn and dCo and varia-tion in dZn/dCo stoichiometry across these two distinct eco-logical regimes in the south-eastern Atlantic.

Description: The ecosystem of the Baltic Sea is endangered by eutrophication. This has triggered expensive international management efforts. Some of these efforts are impeded by natural processes such as nitrogen-fixing cyanobacteria blooms that add bioavailable nitrogen to the already over-fertilized system and thereby enhance primary production, export of organic matter to depth, and associated oxygen consumption. Controls of cyanobacteria blooms are not comprehensively understood, and this adds to the uncertainty of model-based projections into the warming future of the Baltic Sea. Here we review our current understanding of cyanobacteria bloom dynamics. We summarize published field studies and laboratory experiments and dissect the basic principles ingrained in state-of-the-art coupled ocean–circulation biogeochemical models.

Description: A large number of water masses are presented in the Atlantic Ocean, and knowledge of their distributions and properties is important for understanding and monitoring of a range of oceanographic phenomena. The characteristics and distributions of water masses in biogeochemical space are useful for, in particular, chemical and biological oceanography to understand the origin and mixing history of water samples. Here, we define the characteristics of the major water masses in the Atlantic Ocean as source water types (SWTs) from their formation areas, and map out their distributions. The SWTs are described by six properties taken from the biased-adjusted Global Ocean Data Analysis Project version 2 (GLODAPv2) data product, including both conservative (conservative temperature and absolute salinity) and non-conservative (oxygen, silicate, phosphate and nitrate) properties. The distributions of these water masses are investigated with the use of the optimum multi-parameter (OMP) method and mapped out. The Atlantic Ocean is divided into four vertical layers by distinct neutral densities and four zonal layers to guide the identification and characterization. The water masses in the upper layer originate from wintertime subduction and are defined as central waters. Below the upper layer, the intermediate layer consists of three main water masses: Antarctic Intermediate Water (AAIW), Subarctic Intermediate Water (SAIW) and Mediterranean Water (MW). The North Atlantic Deep Water (NADW, divided into its upper and lower components) is the dominating water mass in the deep and overflow layer. The origin of both the upper and lower NADW is the Labrador Sea Water (LSW), the Iceland–Scotland Overflow Water (ISOW) and the Denmark Strait Overflow Water (DSOW). The Antarctic Bottom Water (AABW) is the only natural water mass in the bottom layer, and this water mass is redefined as Northeast Atlantic Bottom Water (NEABW) in the north of the Equator due to the change of key properties, especially silicate. Similar with NADW, two additional water masses, Circumpolar Deep Water (CDW) and Weddell Sea Bottom Water (WSBW), are defined in the Weddell Sea region in order to understand the origin of AABW.

Description: Proxy climate records are an invaluable source of information about the earth's climate prior to the instrumental record. The temporal and spatial coverage of records continues to increase; however, these records of past climate are associated with significant uncertainties due to non-climate processes that influence the recorded and measured proxy values. Generally, these uncertainties are timescale dependent and correlated in time. Accounting for structure in the errors is essential for providing realistic error estimates for smoothed or stacked records, detecting anomalies, and identifying trends, but this structure is seldom accounted for. In the first of these companion articles, we outlined a theoretical framework for handling proxy uncertainties by deriving the power spectrum of proxy error components from which it is possible to obtain timescale-dependent error estimates. Here in Part 2, we demonstrate the practical application of this theoretical framework using the example of marine sediment cores. We consider how to obtain estimates for the required parameters and give examples of the application of this approach for typical marine sediment proxy records. Our new approach of estimating and providing timescale-dependent proxy errors overcomes the limitations of simplistic single-value error estimates. We aim to provide the conceptual basis for a more quantitative use of paleo-records for applications such as model–data comparison, regional and global synthesis of past climate states, and data assimilation.

Description: Bottom pressure observations on both sides of the Atlantic basin, combined with satellite measurements of sea level anomalies and wind stress data, are utilized to estimate variations of the Atlantic Meridional Overturning Circulation (AMOC) at 11∘ S. Over the period 2013–2018, the AMOC and its components are dominated by seasonal variability, with peak-to-peak amplitudes of 12 Sv for the upper-ocean geostrophic transport, 7 Sv for the Ekman and 14 Sv for the AMOC transport. The characteristics of the observed seasonal cycles of the AMOC and its components are compared to results from an ocean general circulation model, which is known to reproduce the variability of the Western Boundary Current on longer timescales. The observed seasonal variability of zonally integrated geostrophic velocity in the upper 300 m is controlled by pressure variations at the eastern boundary, while at 500 m depth contributions from the western and eastern boundaries are similar. The model tends to underestimate the seasonal pressure variability at 300 and 500 m depth, especially at the western boundary, which translates into the estimate of the upper-ocean geostrophic transport. In the model, seasonal AMOC variability at 11∘ S is governed, besides the Ekman transport, by the geostrophic transport variability in the eastern basin. The geostrophic contribution of the western basin to the seasonal cycle of the AMOC is instead comparably weak, as transport variability in the western basin interior related to local wind curl forcing is mainly compensated by the Western Boundary Current. Our analyses indicate that while some of the uncertainties of our estimates result from the technical aspects of the observational strategy or processes not being properly represented in the model, uncertainties in the wind forcing are particularly relevant for the resulting uncertainties of AMOC estimates at 11∘ S.

Description: Bromoform is the major by-product from chlorination of cooling water in coastal power plants. The number of power plants in East and Southeast Asian economies has increased rapidly, exceeding mean global growth. Bottom-up estimates of bromoform emissions based on few measurements appear to under-represent the industrial sources of bromoform from East Asia. Using oceanic Lagrangian analyses, we assess the amount of bromoform produced from power plant cooling-water treatment in East and Southeast Asia. The spread of bromoform is simulated as passive particles that are advected using the three-dimensional velocity fields over the years 2005/2006 from the high-resolution NEMO-ORCA0083 ocean general circulation model. Simulations are run for three scenarios with varying initial bromoform concentrations based on the range of bromoform measurements in cooling-water discharge. Comparing the modelled anthropogenic bromoform to in situ observations in the surface ocean and atmosphere, the two lower scenarios show the best agreement, suggesting initial bromoform concentrations in cooling water to be around 20–60 µg L−1. Based on these two scenarios, the model produces elevated bromoform in coastal waters of East Asia with average concentrations of 23 and 68 pmol L−1 and maximum values in the Yellow Sea, Sea of Japan and East China Sea. The industrially produced bromoform is quickly emitted into the atmosphere with average air–sea flux of 3.1 and 9.1 nmolm−2h−1 , respectively. Atmospheric abundances of anthropogenic bromoform are derived from simulations with the Lagrangian particle dispersion model FLEXPART based on ERA-Interim wind fields in 2016. In the marine boundary layer of East Asia, the FLEXPART simulations show mean anthropogenic bromoform mixing ratios of 0.4–1.3 ppt, which are 2–6 times larger compared to the climatological bromoform estimate. During boreal winter, the simulations show that some part of the anthropogenic bromoform is transported by the northeasterly winter monsoon towards the tropical regions, whereas during boreal summer anthropogenic bromoform is confined to the Northern Hemisphere subtropics. Convective events in the tropics entrain an additional 0.04–0.05 ppt of anthropogenic bromoform into the stratosphere, averaged over tropical Southeast Asia. In our simulations, only about 10 % of anthropogenic bromoform is outgassed from power plants located in the tropics south of 20∘ N, so that only a small fraction of the anthropogenic bromoform reaches the stratosphere. We conclude that bromoform from cooling-water treatment in East Asia is a significant source of atmospheric bromine and might be responsible for annual emissions of 100–300 Mmol of Br in this region. These anthropogenic bromoform sources from industrial water treatment might be a missing factor in global flux estimates of organic bromine. While the current emissions of industrial bromoform provide a significant contribution to regional tropospheric budgets, they provide only a minor contribution to the stratospheric bromine budget of 0.24–0.30 ppt of Br.

Description: Gully formation has been associated to groundwater seepage in unconsolidated sand- to gravel-sizedsediments. Our understanding of gully evolution by groundwater seepage mostly relies on experiments and nu-merical simulations, and these rarely take into consideration contrasts in lithology and permeability. In addition,process-based observations and detailed instrumental analyses are rare. As a result, we have a poor understandingof the temporal scale of gully formation by groundwater seepage and the influence of geological heterogeneityon their formation. This is particularly the case for coastal gullies, where the role of groundwater in their for-mation and evolution has rarely been assessed. We address these knowledge gaps along the Canterbury coastof the South Island (New Zealand) by integrating field observations, luminescence dating, multi-temporal un-occupied aerial vehicle and satellite data, time domain electromagnetic data and slope stability modelling. Weshow that gully formation is a key process shaping the sandy gravel cliffs of the Canterbury coastline. It is anepisodic process associated to groundwater flow that occurs once every 227 d on average, when rainfall intensi-ties exceed 40 mm d−1. The majority of the gullies in a study area southeast (SE) of Ashburton have undergoneerosion, predominantly by elongation, during the last 11 years, with the most recent episode occurring 3 yearsago. Gullies longer than 200 m are relict features formed by higher groundwater flow and surface erosion〉2 kaago. Gullies can form at rates of up to 30 m d−1via two processes, namely the formation of alcoves and tunnelsby groundwater seepage, followed by retrogressive slope failure due to undermining and a decrease in shearstrength driven by excess pore pressure development. The location of gullies is determined by the occurrenceof hydraulically conductive zones, such as relict braided river channels and possibly tunnels, and of sand lensesexposed across sandy gravel cliffs. We also show that the gully planform shape is generally geometrically similarat consecutive stages of evolution. These outcomes will facilitate the reconstruction and prediction of a prevalenterosive process and overlooked geohazard along the Canterbury coastline.

Description: Following the April 2015 eruption of the Calbuco volcano, an extensive ash plume spread across northern Patagonia and into the SE Pacific and SW Atlantic Ocean. Here we report the results of field surveys conducted in the marine region receiving the highest ash load following the eruption (Reloncaví Fjord). The fortuitous location of a long-term monitoring station in Reloncaví Fjord provided data to evaluate inshore phytoplankton bloom dynamics and carbonate chemistry during April–May 2015. Satellite derived chlorophyll-a measurements over the ocean regions affected by the ash plume in May 2015 were obtained to determine the spatial-temporal gradient in offshore phytoplankton response to ash. Additionally, leaching experiments were performed to quantify the release of total alkalinity, trace elements (Fe, Mn, Pb, Co, Cu, Ni and Cd) and major ions (Fl, Cl, SO4, NO3, Li, Na, NH4, K, Mg, Ca) from ash into solution. Within Reloncaví Fjord, integrated peak diatom abundances during the May 2015 austral bloom were higher than usual (up to 1.4 × 1011 cells m−2, integrated to 15 m depth), with the bloom intensity perhaps moderated due to high ash loadings in the two weeks following the eruption. In the offshore SE Pacific, a short duration phytoplankton bloom corresponded closely in space and time to the maximum observed ash plume, potentially in response to Fe-fertilization of a region where phytoplankton growth is typically Fe-limited at this time of year. Conversely, no clear fertilization was found in the area subject to an ash plume over the SW Atlantic where the availability of fixed nitrogen is thought to limit phytoplankton growth which was consistent with no significant release of fixed nitrogen from ash. In addition to release of nanomolar concentrations of dissolved Fe from ash suspended in seawater, it was observed that low loadings (〈 5 mg L−1) of freshly deposited ash were an unusually prolific source of Fe(II) into solution (up to 1.0 µmol Fe g−1), suggesting that the release of bioaccessible Fe from ash sources may generally be under-estimated when quantified from aged ash. This release of Fe(II) may make freshly deposited ash an unusually efficient dissolved Fe source with the 18–38 % fraction of dissolved Fe released as Fe(II) from Calbuco ash roughly comparable to literature values for Fe released into seawater from aerosols collected over the Pacific Ocean.

Description: In aquatic systems, invasive submerged macrophytes considerably alter the structure and functioning of communities, thus potentially compromising ecosystem services. The prolific spread of invasive macrophytes is often aided by vegetative fragment propagation, yet the contributions of various commonly occurring invertebrates to such fragmentation are often unquantified. In the present study, we examine fragmentary spread of invasive macrophytes by a group of shredder-herbivores, larval caddisflies. Through novel application of the comparative functional response (FR; resource use as a function of density) approach to the native case-building species Limnephilus lunatus, we compared utilisation of non-native waterweeds Elodea canadensis and E. nuttallii in mono- and polycultures. Furthermore, we quantified de-cased and cased caddisfly-induced fragment production and length changes among non-native E. canadensis, E. nuttallii, Crassula helmsii and Lagarosiphon major under two different plant orientations: horizontal (floating) versus vertical (upright) growth forms. Larval caddisflies exhibited Type II (hyperbolic) FRs towards both Elodea species, and utilised each plant at similar rates when plants were provided separately. When plant species were presented in combination horizontally, E. canadensis was significantly less utilised compared to E. nuttallii, corroborating observations in the field. De-cased larvae produced new plant fragments for all four aquatic macrophytes, whereas cased larvae fragmented plants significantly less. Elodea nuttalii and C. helmsii were fragmented the most overall. Crassula helmsii was utilised to the greatest extent when plants were horizontally orientated, and Elodea species when vertically orientated. This study identifies and quantifies a mechanism from a novel species group that may contribute to the spread of invasive macrophytes in aquatic systems. Whilst exploititative interactions are thought to impede invasion success, here we demonstrate how resource utilisation by a resident species may exacerbate propagule pressure from an invasive species.

Description: The Osborn Plateau is a large intraplate elevation in the eastern Indian Ocean, which has been poorly studied by geological and geophysical methods. On cruise SO258/1 aboard the R/V Sonne, new data were collected by Parasound seismic processing, multibeam survey, and, for the first time, dredging. Faults in the sedimentary cover, which extend to the seafloor surface, indicate high neotectonic activity in the area of the Osborn Plateau. It may have continued up to the present, as well as in the adjacent segment of the Ninetyeast Ridge, where strong earthquakes are recorded. Two reflectors A and C in the upper part of the sedimentary cover mark global regressive changes in the World Ocean level at the Pliocene–Pleistocene and Miocene–Pliocene and boundaries. Reflector B in sediments at the Lower–Upper Pliocene boundary is associated with a change in the regional hydrodynamic regime at the time in the eastern Indian Ocean. Reflector B at the Lower–Upper Pliocene boundary is associated with a change in the local hydrodynamic regime in the region of the plateau and was caused by uncompensated sedimentation. As a result of dredging, strongly altered vitroclastic tuffs were obtained, consisting of palagonized ash particles and lapilli. Analysis of geochemical data on the composition of palagonite made it possible to reconstruct the main geochemical features of the primary composition of glass in tuffs, in particular, high enrichment in incompatible elements (Nb, Zr, Ti, La, etc.). The results of studying the dredged rocks suggest that tuffs dredged on the Osborn Plateau were formed as a result of explosive volcanic eruptions of alkaline basalts under subaerial or relatively shallow conditions and represent the latest eruption products in the region.

Description: In a conserved culture of the purple sulfur bacterium Thiospirillum jenense DSM216T, cells of this species were easily recognized by cell morphology, large-size spirilla and visible flagellar tuft. The Tsp. jenense genome is 3.22 Mb in size and has a GC content of 48.7 mol%. It was readily identified as a member of the Chromatiaceae by the complement of proteins in its genome. A whole genome comparison clearly placed Tsp. jenense near Thiorhodovibrio and Rhabdochromatium species and somewhat more distant from Thiohalocapsa and Halochromatium species. This relationship was also found with the sequences of the photosynthetic reaction center protein PufM. The genome sequence supported important properties of this bacterium: the presence of ribulose-bisphosphate carboxylase and enzymes of the Calvin cycle of autotrophic carbon dioxide fixation but the absence of carboxysomes, an incomplete tricarboxylic acid cycle and the lack of malate dehydrogenase, the presence of a sulfur oxidation pathway including adenylylsulfate reductase (aprAB) but absence of assimilatory sulfate reduction, the presence of hydrogenase (hoxHMFYUFE), nitrogenase and a photosynthetic gene cluster (pufBALMC). The FixNOP type of cytochrome oxidase was notably lacking, which may be the reason that renders the cells highly sensitive to oxygen. Two minor phototrophic contaminants were found using metagenomic binning: one was identified as a strain of Rhodopseudomonas palustris and the second one has an average nucleotide identity of 82% to the nearest neighbor Rhodoferax antarcticus. It should be considered as a new species of this genus and Rhodoferax jenense is proposed as the name.

Description: For many coastal areas including the Baltic Sea, ambitious nutrient abatement goals have been set to curb eutrophication, but benefits of such measures were normally not studied in light of anticipated climate change. To project the likely responses of nutrient abatement on eelgrass (Zostera marina), we coupled a species distribution model with a biogeochemical model, obtaining future water turbidity, and a wave model for predicting the future hydrodynamics in the coastal area. Using this, eelgrass distribution was modeled for different combinations of nutrient scenarios and future wind fields. We are the first to demonstrate that while under a business as usual scenario overall eelgrass area will not recover, nutrient reductions that fulfill the Helsinki Commission’s Baltic Sea Action Plan (BSAP) are likely to lead to a substantial areal expansion of eelgrass coverage, primarily at the current distribution’s lower depth limits, thereby overcompensating losses in shallow areas caused by a stormier climate.

Description: Filaments and fronts play a crucial role for a net offshore and downward nutrient transport in Eastern Boundary Upwelling Regions (EBUS) and thereby reduce primary production. Often studies are either based on observations or model simulations but seldom both approaches are combined quantitatively to assess the importance of filaments for primary production and nutrient transport. Here we combine targeted interdisciplinary shipboard observations of a cold filament off Peru with submesoscale-permitting (1/45°) coupled physical (CROCO) and biogeochemical (PISCES) model simulations to (i) evaluate the model simulations in detail, including the timescales of biogeochemical modification of the newly upwelled water and (ii) quantify the net effect of submesoscale fronts and filaments on primary production of the Peruvian upwelling system. The observed filament contains relatively cold, fresh and nutrient-rich waters originating in the coastal upwelling. Enhanced nitrate concentrations and offshore velocities of up to 0.5 m s−1 within the filament suggest an offshore transport of nutrients. Surface chlorophyll in the filament is a factor 4 lower than at the upwelling front while surface primary production is a factor 2 higher, highlighting the additional value of direct rate measurements for model validation. The simulation exhibits filaments that are similar in horizontal and vertical scale compared to the observed filament. Nitrate concentrations and primary pro- duction within filaments in the model are comparable to observations as well, justifying further analysis of nitrate uptake and subduction using the model. Virtual Lagrangian floats were released in the subsurface waters along the shelf and biogeochemical variables tracked along the trajectories of floats upwelled near the coast. In the submesoscale-permitting (1/45°) simulation 43.0 % of upwelled floats and 40.6 % of upwelled nitrate is subducted within 20 days after upwelling, which corresponds to an increase of nitrate subduction compared to a mesoscale-resolving (1/9°) simulation by 13.9 %. This suggests that submesoscale processes further reduce primary production by amplifying the downward and offshore export of nutrients found in previous mesoscale studies, which are thus likely to underestimate the reduction in primary production due to eddy-fluxes. Moreover, this downward and offshore transport could also enhance the export of fresh organic matter below the photic zone and thereby potentially stimulate microbial activity in the upper offshore oxygen minimum zone.

Description: The sea surface microlayer (SML) represents the boundary layer at the air–sea interface. Microbial eukaryotes in the SML potentially influence air–sea gas exchange directly by taking up and producing gases and indirectly by excreting and degrading organic matter, which may modify the viscoelastic properties of the SML. However, little is known about the distribution of microbial eukaryotes in the SML. We studied the composition of the microbial community, transparent exopolymer particles and polysaccharides in the SML during the PEACETIME cruise along a west–east transect in the Mediterranean Sea, covering the western basin, Tyrrhenian Sea and Ionian Sea. At the stations located in the Ionian Sea, fungi – likely of continental origin and delivered by atmospheric deposition – were found in high relative abundances, making up a significant proportion of the sequences recovered. Concomitantly, bacterial and picophytoplankton counts decreased from west to east, while transparent exopolymer particle (TEP) abundance and total carbohydrate (TCHO) concentrations remained constant in all basins. Our results suggest that the presence of substrates for fungi, such as Cladosporium, known to take up phytoplankton-derived polysaccharides, in combination with decreased substrate competition by bacteria, might favor fungal dominance in the neuston of the Ionian Sea and other low-nutrient, low-chlorophyll (LNLC) regions.

Description: Predicting the ecological impacts of invasive species is impeded by context-dependencies that can mediate the strength of trophic interactions. Reproductive status is a pervasive context across trophic and taxonomic groups, yet has mostly been neglected in studies of invader impact. The present study examines the influence of the reproductive cycle on predatory impacts of a known damaging invasive freshwater crustacean, Gammarus pulex (Amphipoda), using a comparative functional response approach. Across females that were non-ovigerous, or had immature- or mature-stage embryonic broods, all G. pulex exhibited potentially prey destabilising Type II functional responses towards chironomid larvae. Attack rates were highest by immature embryonic brood stages, and significantly lower by mature embryonic brood stages. Conversely, handling times were consistently lower, and hence maximum feeding rates higher, where broods were present as compared to absent. These predatory patterns may reflect changing resource demands for progeny development across the reproductive cycle, the influence of brood bulk, female moult stage and/or motivational aspects of feeding, such as avoidance of filial cannibalism. Accordingly, many aspects of reproduction may alter the intensity of predatory impacts by invasive species. Future studies considering invader impacts should consider the influence of sex and reproductive status to more holistically quantify and predict population-level invader impacts.

Description: Species invasions often occur on coasts and estuaries where abiotic conditions vary, e.g. salinity, temperature, runoff etc. Successful establishment and dispersal of non-indigenous species in many such systems are poorly understood, partially since the species tend to show genetic and ecological plasticity at population level towards many abiotic conditions, including salinity tolerance. Plasticity may be driven by shifting expression of heat shock proteins such as Hsp70, which is widely recognized as indicator of physical stress. In this study, we developed a qPCR assay for expression of the hsp70 gene in the invasive round goby (Neogobius melanostomus) and tested the expression response of fish collected from a brackish environment in the western Baltic Sea to three different salinities, 0, 10 and 30. hsp70 expression was highest in fresh water, indicating higher stress, and lower at brackish (ambient condition for the sampled population) and oceanic salinities, suggestive of low stress response to salinities above the population’s current distribution. The highest stress in fresh water was surprising since populations in fresh water exist, e.g. large European rivers and Laurentian Great Lakes. The results have implications to predictions for the species’ plasticity potential and possible range expansion of the species into other salinity regimes.

Description: Ocean deoxygenation is a rising threat to marine ecosystems and food resources under present climate warming conditions. Organic-rich sapropel layers deposited in the Mediterranean Sea provide a natural laboratory to study the processes that have controlled changes in seawater oxygen levels in the recent geological past. Our study is based on three sediment cores spanning the last 10 000 years and located on a bathymetric transect offshore from the western distributaries of the Nile delta. These cores are partly to continuously laminated in the sections recording sapropel S1, which is indicative of bottom-water anoxia above the western Nile deep-sea fan. We used a combination of microfacies analyses and inorganic and organic geochemical measurements to reconstruct changes in oxygenation conditions at seasonal to millennial timescales. Millimetre-thick laminations are composed of detrital, biogenic and chemogenic sublayers reflecting seasonal successions of sedimentation. Dark layers reflect the deposition of summer floods and two types of light layers correspond to autumn plankton blooms and authigenic carbonates formed in the water column during spring–early summer, respectively. The isotopic signature of the authigenic carbonates suggests permanent anoxic to euxinic bottom waters resulting in high levels of anaerobic remineralization of organic matter and highlights their potential to reconstruct seawater chemistry at times when benthic fauna was absent. Ratios of major elements combined with biomarkers of terrestrial and marine organic matter and redox-sensitive compounds allow changes in terrigenous input, primary productivity and past deoxygenation dynamics on millennial timescales to be tracked. Rapid fluctuations of oxygenation conditions in the upper 700 m water depth occurred above the Nile deep-sea fan between 10 and 6.5 ka BP, while deeper cores recorded more stable anoxic conditions. Synchronous changes in terrigenous input, primary productivity and past oxygenation dynamics after 6.5 ka BP show that runoff-driven eutrophication played a central role in rapid oxygenation changes in the south-eastern Levantine Basin. These findings are further supported by other regional records and reveal time-transgressive changes in oxygenation state driven by rapid changes in primary productivity during a period of long-term deep-water stagnation.

Description: Understanding mechanisms of tropical Pacific decadal variability (TPDV) is of high importance for differentiating between natural climate variability and human induced climate change as this region sustains strong global teleconnections. Here, we use an ocean general circulation model along with a Lagrangian tracer simulator to investigate the advection of density compensated temperature anomalies (“spiciness mechanism”) as a potential contributor to TPDV during the 1980–2016 period. Consistent with observations, we find the primary regions of spiciness generation in the eastern subtropics of each hemisphere. Our results indicate that 75% of the equatorial subsurface water originates in the subtropics, of which two thirds come from the Southern hemisphere. We further show two prominent cases where remotely generated spiciness anomalies are advected to the equatorial Pacific, impacting subsurface temperature. The relative contribution of Northern versus Southern Hemisphere prominence and/or interior versus western boundary pathways depends on the specific event. The anomalously warm case largely results from advection via the Southern hemisphere interior (65%), while the anomalously cold case largely results from advection via the Northern hemisphere western boundary (48%). The relatively slow travel times from the subtropics to the equator (〉 4 years) suggests that these spiciness anomalies underpin a potentially predictable contribution to TPDV. However, not all decadal peaks in equatorial spiciness can be explained by remotely generated spiciness anomalies. In those cases, we propose that spiciness anomalies are generated in the equatorial zone through changes in the proportion of Northern/Southern hemisphere source waters due to their different mean spiciness distribution.

Description: Numerical models have become indispensable tools for investigating submarine hydrothermal systems and for relating seafloor observations to physicochemical processes at depth. Particularly useful are multiphase models that account for phase separation phenomena, so that model predictions can be compared to observed variations in vent fluid salinity. Yet, the numerics of multiphase flow remain a challenge. Here we present a novel hydrothermal flow model for the system H 2 O–NaCl able to resolve multiphase flow over the full range of pressure, temperature, and salinity variations that are relevant to submarine hydrothermal systems. The method is based on a 2-D finite volume scheme that uses a Newton–Raphson algorithm to couple the governing conservation equations and to treat the non-linearity of the fluid properties. The method uses pressure, specific fluid enthalpy, and bulk fluid salt content as primary variables, is not bounded to the Courant time step size, and allows for a direct control of how accurately mass and energy conservation is ensured. In a first application of this new model, we investigate brine formation and mobilization in hydrothermal systems driven by a transient basal temperature boundary condition—analogue to seawater circulation systems found at mid-ocean ridges. We find that basal heating results in the rapid formation of a stable brine layer that thermally insulates the driving heat source. While this brine layer is stable under steady-state conditions, it can be mobilized as a consequence of variations in heat input leading to brine entrainment and the venting of highly saline fluids.

Description: Methane (CH4) is one of the substantial greenhouse gases in our atmosphere and its concentration has increased by ~ 4 % over the last decade. Although sources driving these increases are not well constrained, one potential contribution comes from wetlands, which are usually intertwined with rivers, channels and lakes, creating a considerable need to acquire higher resolution data to facilitate 5 modelling and predictions. Here we took a fully contained sensor set-up to obtain measurements of CO2, CH4, O2 and auxiliary parameters, installed on a houseboat for accessibility, to assess and analyse surface water concentrations within the Danube Delta, Romania. Over 3 seasons, we transected a ~ 400 km route with concentration mapping and additional stations for monitoring diel cycles. Overall, the delta was a source for CH4 throughout all seasons, with concentrations ranging between 0.113–15.6 μmol L−1. The dataset was split into three different subsystems; lakes, rivers and channels, with channels 10 showing the highest variability. We found large to extreme diel cycles in both the lakes and channels, with concentrations varying by an order of magnitude between these two systems. The observed strong diel cycle within the lake suggests daily vertical stratification allowing for macrophytes to create a temporal oxycline due to lack of light and movement between the stems as previously suggested. While throughout the day, there was a consistent overall surface concentration of CH4 at around 0.4 μmol L−1, there was a clear linear trend with an O2:CH4 molar ratio of −50:1 during the phase of nocturnal convection 15 with the two water stratified bodies mixing during the night. Daily spot sampling techniques and neglecting such diel cycles reducing the estimated average methane concentrations by 25 % and increase by 3.3 % for channels and lakes, respectively. On an individual lake basis, spot sampling can potentially incur an uncertainty range of a factor of 4.5. Analyses also included a ‘hot spot’, with a 10-fold stronger methane increase (4–15.6 μmol L−1) overnight compared to the lake, with an almost immediate and extreme decrease in CH4 following sunrise. Calculated diffusive CH4 fluxes for the overall delta yielded an 20 average of 49 ± 61 μmol m−2 h−1 corresponding to an extrapolated annual flux of 0.43 ± 0.53 μmol m−2 yr−1. Our data illustrate the importance of collecting information on diel cycles in different habitats to improve the emission estimates from wetland systems.

Description: The growth of El Niño/Southern Oscillation (ENSO) events is determined by the balance between ocean dynamics and thermodynamics. Here we quantify the contribution of the thermodynamic feedbacks to the sea surface temperature (SST) change during ENSO growth phase by integrating the atmospheric heat fluxes over the temporarily and spatially varying mixed layer to derive an offline “slab ocean” SST. The SST change due to ocean dynamics is estimated as the residual with respect to the total SST change. In observations, 1 K SST change in the Niño3.4 region is composed of an ocean dynamical SST forcing of + 2.6 K and a thermodynamic damping of − 1.6 K, the latter mainly by the shortwave-SST (− 0.9 K) and latent heat flux-SST feedback (− 0.7 K). Most climate models from the Coupled Model Intercomparison Project phase 5 (CMIP5) underestimate the SST change due to both ocean dynamics and net surface heat fluxes, revealing an error compensation between a too weak forcing by ocean dynamics and a too weak damping by atmospheric heat fluxes. In half of the CMIP5 models investigated in this study, the shortwave-SST feedback erroneously acts as an amplifying feedback over the eastern equatorial Pacific, resulting in a hybrid of ocean-driven and shortwave-driven ENSO dynamics. Further, the phase locking and asymmetry of ENSO is investigated in the CMIP5 model ensemble. The climate models with stronger atmospheric feedbacks tend to simulate a more realistic seasonality and asymmetry of the heat flux feedbacks, and they exhibit more realistic phase locking and asymmetry of ENSO. Moreover, the almost linear latent heat flux feedback contributes to ENSO asymmetry in the far eastern equatorial Pacific through an asymmetry in the mixed layer depth. This study suggests that the dynamic and thermodynamic ENSO feedbacks and their seasonality and asymmetries are important metrics to consider for improving ENSO representation in climate models.

Description: The Baltic Sea has a salinity gradient decreasing from fully marine (〉 25) in the west to below 7 in the central Baltic Proper. Habitat-forming and ecologically dominant mytilid mussels exhibit decreasing growth when salinity 〈 11; however, the mechanisms underlying reduced calcification rates in dilute seawater are not fully understood. Both [HCO−3] and [Ca2+] also decrease with salinity, challenging calcifying organisms through CaCO3 undersaturation (Ω≤1) and unfavourable ratios of calcification substrates ([Ca2+] and [HCO−3]) to the inhibitor (H+), expressed as the extended substrate–inhibitor ratio (ESIR). This study combined in situ monitoring of three southwest Baltic mussel reefs with two laboratory experiments to assess how various environmental conditions and isolated abiotic factors (salinity, [Ca2+], [HCO−3] and pH) impact calcification in mytilid mussels along the Baltic salinity gradient. Laboratory experiments rearing juvenile Baltic Mytilus at a range of salinities (6, 11 and 16), HCO−3 concentrations (300–2100 µmol kg−1) and Ca2+ concentrations (0.5–4 mmol kg−1) reveal that as individual factors, low [HCO−3], pH and salinity cannot explain low calcification rates in the Baltic Sea. Calcification rates are impeded when Ωaragonite ≤ 1 or ESIR ≤ 0.7 primarily due to [Ca2+] limitation which becomes relevant at a salinity of ca. 11 in the Baltic Sea. Field monitoring of carbonate chemistry and calcification rates suggest increased food availability may be able to mask the negative impacts of periodic sub-optimal carbonate chemistry, but not when seawater conditions are permanently adverse, as observed in two Baltic reefs at salinities 〈 11. Regional climate models predict a rapid desalination of the southwest and central Baltic over the next century and potentially a reduction in [Ca2+] which may shift the distribution of marine calcifiers westward. It is therefore vital to understand the mechanisms by which the ionic composition of seawater impacts bivalve calcification for better predicting the future of benthic Baltic ecosystems.

Description: Dimethyl sulfide (DMS), dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) were measured at the Boknis Eck Time Series Station (BE, Eckernförde Bay, SW Baltic Sea) during the period February 2009–December 2018. Our results show considerable interannual and seasonal variabilities in the mixed-layer concentrations of DMS, total DMSP (DMSPt) and total DMSO (DMSOt). Positive correlations were found between particulate DMSP (DMSPp) and particulate DMSO (DMSOp) as well as DMSPt and DMSOt in the mixed layer, suggesting a similar source for both compounds. The decreasing long-term trends, observed for DMSPt and DMS in the mixed layer, were linked to the concurrent trend of the sum of 19′-hexanoyloxyfucoxanthin and 19′-butanoyloxy-fucoxanthin, which are the marker pigments of prymnesiophytes and chrysophytes, respectively. Major Baltic inflow (MBI) events influenced the distribution of sulfur compounds due to phytoplankton community changes, and sediment might be a potential source for DMS in the bottom layer during seasonal hypoxia/anoxia at BE. A modified algorithm based on the phytoplankton pigments reproduces the DMSPp : Chl a ratios well during this study and could be used to estimate future surface (5 m) DMSPp concentrations at BE.

Description: As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated prestudy experiment with interactive stratospheric aerosol models simulating the volcanic aerosol cloud from an eruption resembling the 1815 Mt. Tambora eruption (VolMIP-Tambora ISA ensemble). The pre-study provided the ancillary ability to assess intermodel diversity in the radiative forcing for a large stratospheric-injecting equatorial eruption when the volcanic aerosol cloud is simulated interactively. An initial analysis of the VolMIP-Tambora ISA ensemble showed large disparities between models in the stratospheric global mean aerosol optical depth (AOD). In this study, we now show that stratospheric global mean AOD differences among the participating models are primarily due to differences in aerosol size, which we track here by effective radius. We identify specific physical and chemical processes that are missing in some models and/or parameterized differently between models, which are together causing the differences in effective radius. In particular, our analysis indicates that interactively tracking hydroxyl radical (OH) chemistry following a large volcanic injection of sulfur dioxide (SO2) is an important fac tor in allowing for the timescale for sulfate formation to be properly simulated. In addition, depending on the timescale of sulfate formation, there can be a large difference in effective radius and subsequently AOD that results from whether the SO2 is injected in a single model grid cell near the location of the volcanic eruption, or whether it is injected as a longitudinally averaged band around the Earth.

Description: Small pelagic fish off the coast of Peru in the Eastern Tropical South Pacific (ETSP) support around 10 % of the global fish catches. Their stocks fluctuate interannually due to environmental variability which can be exacerbated by fishing pressure. Because these fish are planktivorous, any change in fish abundance may directly affect the plankton and the biogeochemical system. To investigate the potential effects of variability in small pelagic fish populations on lower trophic levels, we used a coupled physical-biogeochemical model to build scenarios for the ETSP and compare these against an already published reference simulation. The scenarios mimic changes in fish predation by either increasing or decreasing mortality of the model's large and small zooplankton compartments. The results revealed that large zooplankton was the main driver of the response of the community. Its concentration increased under low mortality conditions and its prey, small zooplankton and large phytoplankton, decreased. The response was opposite, but weaker, in the high mortality scenarios. This asymmetric behaviour can be explained by the different ecological roles of large, omnivorous zooplankton, and small zooplankton, which in the model is strictly herbivorous. The response of small zooplankton depended on the antagonistic effects of mortality changes as well as the grazing pressure by large zooplankton. The results of this study provide a first insight on how the plankton ecosystem might respond if variations in fish populations were modelled explicitly.

Description: Trace elements play important roles as micronutrients in modulating marine productivity in the global ocean. The South Atlantic around 40° S is a prominent region of high productivity and a transition zone between the nitrate-depleted Subtropical Gyre and the iron-limited Southern Ocean. However, the sources and fluxes of trace elements to this region remain unclear. In this study, the distribution of the naturally occurring radioisotope 228Ra in the water column of the South Atlantic (Cape Basin and Argentine Basin) has been investigated along a 40° S zonal transect to estimate ocean mixing and trace element supply to the surface ocean. Ra-228 profiles have been used to determine the horizontal and vertical mixing rates in the near-surface open ocean. In the Argentine Basin, horizontal mixing from the continental shelf to the open ocean shows an eddy diffusion of Kx = 1.7 ± 1.4 (106 cm2 s−1) and an integrated advection velocity w = 0.6 ± 0.3 cm s−1. In the Cape Basin, horizontal mixing is Kx = 2.7 ± 0.8 (107 cm2 s−1) and vertical mixing Kz = 1.0–1.5 cm2 s−1 in the upper 600 m layer. Three different approaches (228Ra-diffusion, 228Ra-advection and 228Ra/TE-ratio) have been applied to estimate the dissolved trace-element fluxes from shelf to open ocean. These approaches bracket the possible range of off-shelf fluxes from the Argentine margin to be: 3.8–22 (× 103) nmol Co m−2 d−1, 7.9–20 (× 104) nmol Fe m−2 d−1 and 2.7–6.5 (× 104) nmol Zn m−2 d−1. Off-shelf fluxes from the Cape margin are: 4.3–6.2 (× 103) nmol Co m−2 d−1, 1.2–3.1 (× 104) nmol Fe m−2 d−1 and 0.9–1.2 (× 104) nmol Zn m−2 d−1. On average, at 40° S in the Atlantic, vertical mixing supplies 0.4–1.2 nmol Co m−2 d−1, 3.6–11 nmol Fe m−2 d−1, and 13–16 nmol Zn m−2 d−1 to the euphotic zone. Compared with atmospheric dust and continental shelf inputs, vertical mixing is a more important source for supplying dissolved trace elements to the surface 40° S Atlantic. It is insufficient, however, to provide the trace elements removed by biological uptake. Other inputs (e.g. particulate, or from winter deep-mixing) are required to balance the trace element budgets in this region.

Description: We present a climatology of the near-sea-surface temperature (NSST) anomaly and the sea-ice extent during the Last Glacial Maximum (LGM, 23 000–19 000 years before present) mapped on a global regular 1∘×1∘ grid. It is an extension of the Glacial Atlantic Ocean Mapping (GLAMAP) reconstruction of the Atlantic NSST based on the faunal and floral assemblage data of the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project and several recent estimates of the LGM sea-ice extent. Such a gridded climatology is highly useful for the visualization of the LGM climate, calculation of global and regional NSST averages, and estimation of the equilibrium climate sensitivity, as well as a boundary condition for atmospheric general circulation models. The gridding of the sparse NSST reconstruction was done in an optimal way using the Data-Interpolating Variational Analysis (DIVA) software, which takes into account the uncertainty in the reconstruction and includes the calculation of an error field. The resulting Glacial Ocean Map (GLOMAP) confirms the previous findings by the MARGO project regarding longitudinal and meridional NSST differences that were greater than today in all oceans. Taken at face value, the estimated global and tropical cooling would imply an equilibrium climate sensitivity at the lower end of the currently accepted range. However, because of anticipated changes in the seasonality and thermal structure of the upper ocean during the LGM as well as uneven spatial sampling, the estimated cooling and implied climate sensitivity are likely to be biased towards lower values.

Description: Strong regional sea-level trends, mainly related to basin-wide wind stress anomalies, have been observed in the western tropical Pacific over the last 3 decades. Analyses of regional sea level in the densely populated regions of the neighbouring Australasian Mediterranean Sea (AMS; also called tropical Asian seas) are hindered by its complex topography and respective studies are sparse. We used a series of global eddy-permitting ocean models, including a high-resolution configuration that resolves the AMS with 120∘ horizontal resolution, forced by a comprehensive atmospheric forcing product over 1958–2016 to characterize the patterns and magnitude of decadal sea-level variability in the AMS. The nature of this variability is elucidated further by sensitivity experiments with interannual variability restricted to either the momentum or buoyancy fluxes, building on an experiment employing a repeated-year forcing without interannual variability in all forcing components. Our results suggest that decadal fluctuations of the El Niño–Southern Oscillation (ENSO) account for over 80 % of the variability in all deep basins of the region, except for the central South China Sea (SCS). Changes related to the Pacific Decadal Oscillation (PDO) are most pronounced in the shallow Arafura and Timor seas and in the central SCS. On average, buoyancy fluxes account for less than 10 % of decadal SSH variability, but this ratio is highly variable over time and can reach values of up to 50 %. In particular, our results suggest that buoyancy flux forcing amplifies the dominant wind-stress-driven anomalies related to ENSO cycles. Intrinsic variability is mostly negligible except in the SCS, where it accounts for 25 % of the total decadal SSH variability.

Description: In the Comechingones pegmatitic field, central Argentina, leucogranite and pegmatite bodies crop out in a relatively narrow (25 × 10 km) belt, and were emplaced synkinematically with the main deformational event of the crustal-scale Guacha Corral shear zone during the Early Ordovician (~ 475 Ma). These leucogranites and pegmatites are geochemically evolved rocks with high silica and alkalis, low Fe2O3, MgO, TiO2 and CaO, and high ASI values. The leucogranites display quite variable Sr and Nd isotope compositions (initial 87Sr/86Sr ratios from 0.7048 to 0.7170, and εNd values from + 2.0 to − 3.1), some of which do not overlap with almost any other pre-Famatinian rock from the Sierras de Córdoba. The major and trace element geochemistry and the particular Sr and Nd isotope compositions of the leucogranites are here explained by the following processes: (1) water-fluxed partial melting of amphibolites at relatively low P–T conditions generating currently unexposed granodioritic melts with unradiogenic 87Sr/86Sr ratios and radiogenic εNd values; (2) fractionation of mostly plagioclase and monazite leading to compositions close to the leucogranite melts; and (3) assimilation of metasedimentary rocks with crustal isotopic signatures, modelled by assimilation and fractional crystallization processes. The major, trace and isotope compositions of the pegmatites suggest a derivation from partial melting of the same metasedimentary protoliths of the Sierras de Córdoba that were assimilated by leucogranite melts. We propose a feedback relationship among deformation, anatexis, magma evolution and mass transfer in the context of such a crustal-scale shear zone in the foreland of the Famatinian orogen.

Description: We present results from an ensemble of eight climate models, each of which has carried out simulations of the early Eocene climate optimum (EECO, ∼ 50 million years ago). These simulations have been carried out in the framework of the Deep-Time Model Intercomparison Project (DeepMIP; http://www.deepmip.org, last access: 10 January 2021); thus, all models have been configured with the same paleogeographic and vegetation boundary conditions. The results indicate that these non-CO2 boundary conditions contribute between 3 and 5 ∘C to Eocene warmth. Compared with results from previous studies, the DeepMIP simulations generally show a reduced spread of the global mean surface temperature response across the ensemble for a given atmospheric CO2 concentration as well as an increased climate sensitivity on average. An energy balance analysis of the model ensemble indicates that global mean warming in the Eocene compared with the preindustrial period mostly arises from decreases in emissivity due to the elevated CO2 concentration (and associated water vapour and long-wave cloud feedbacks), whereas the reduction in the Eocene in terms of the meridional temperature gradient is primarily due to emissivity and albedo changes owing to the non-CO2 boundary conditions (i.e. the removal of the Antarctic ice sheet and changes in vegetation). Three of the models (the Community Earth System Model, CESM; the Geophysical Fluid Dynamics Laboratory, GFDL, model; and the Norwegian Earth System Model, NorESM) show results that are consistent with the proxies in terms of the global mean temperature, meridional SST gradient, and CO2, without prescribing changes to model parameters. In addition, many of the models agree well with the first-order spatial patterns in the SST proxies. However, at a more regional scale, the models lack skill. In particular, the modelled anomalies are substantially lower than those indicated by the proxies in the southwest Pacific; here, modelled continental surface air temperature anomalies are more consistent with surface air temperature proxies, implying a possible inconsistency between marine and terrestrial temperatures in either the proxies or models in this region. Our aim is that the documentation of the large-scale features and model–data comparison presented herein will pave the way to further studies that explore aspects of the model simulations in more detail, for example the ocean circulation, hydrological cycle, and modes of variability, and encourage sensitivity studies to aspects such as paleogeography, orbital configuration, and aerosols.

Description: Ocean warming is affecting marine ectothermic herbivores as well as the macroalgal species they consume and this has the potential to alter their trophic interaction. However, it is currently still unknown how these two important components of benthic food webs will react to a warming environment. Consumption rates of grazers change with increasing temperatures, but it is unclear whether this is also true for feeding preferences. In this study, multiple-choice feeding assays with the tropical abalone Haliotis squamata from Western Indonesia were conducted in August 2018. After brief acclimation of either the grazer or the macroalgae to moderately elevated water temperatures (maximum 2 °C above the long-term average) in the laboratory, three species of living macroalgae were simultaneously offered to the abalone in feeding assays. Consumption rates of H. squamata were lower under elevated water temperatures, while its feeding preference switched: At 27 °C (2 °C below long-term average), abalone preferred non-acclimated Gracilaria salicornia, but switched to non-acclimated Amphiroa spp. at 31 °C. Interestingly, no such switch in preference occurred when the macroalgae, but not the grazers were acclimated. This indicates that the grazer will presumably be the driver of this potential change in interactions between H. squamata and its macroalgal food. Ocean warming may result in changes in the structure of benthic communities, mediated by changes in the feeding behaviour of herbivorous invertebrates.

Description: In 2019 the RV Meteor cruised from Guadalupe in April/May to Cape Verde in June/July and to Namibia in August/September. The distance is about 10,000 km. The ship has a moon pool for installation of instruments. In Cape Verde we had a first glimpse of the already sparsely populated moon pool. We reached Namibian waters in mid-August. In mid-September, 47 days later and 6000 km south, the ship's moon pool was sampled in the port of Walvis Bay. 13 different taxa could be identified belonging to two phyla, four classes, six orders and 10 families. Most of these species have not yet been observed in the port or in the adjacent areas and are new records for the entire Namibian coast. The goose barnacles Conchoderma auritum (Linnaeus, 1767), Conchoderma virgatum Spengler, 1789 and Lepas anatifera Linnaeus, 1758 were particularly noticeable. They were only surpassed by the large number of amphipods. The species Ericthonius brasiliensis (Dana, 1853), Jassa marmorata Holmes, 1905, Stenothoe senegalensis Krapp-Schickel, 2015 and Paracaprella pusilla Mayer, 1890 are particularly noteworthy here. In addition, the pycnogonid species Endeis straughani Clark, 1970 and the titan acorn barnacle Megabalanus coccopoma (Darwin, 1854) should be mentioned, which occurred very frequently as well. The present study shows, on the one hand, an example of the transatlantic spread of bioinvasive species by ships as vectors and, on the other hand, a convenient method for sampling ship hulls.

Description: The extensive use of silver nanoparticles (AgNPs) in consumer products, medicine, and industry leads to their release into the environment. Thus, a characterization of the concentration, size, fate, and toxicity of AgNPs under environmental conditions is required. In this study, we present the characterization and optimization of an asymmetric flow field-flow fractionation (AF4) system coupled with UV/Vis spectrophotometer and dynamic light scattering (DLS) detector as a powerful tool for the size separation and multi-parameter characterization of AgNPs in complex matrices. The hyphenated AF4-UV/Vis-DLS system was first characterized using individual injections of the different size fractions. We used electrostatically stabilized AgNPs of 20-, 50-, and 80-nm nominal diameters coated with lipoic acid. We investigated the effect of applied cross-flows, carrier solutions, focus times, and quantity of injected particles on the nature of the AF4 fractograms and on the integrity of the AgNPs. Best size separation of a 1:1 mixture of 20- and 80-nm AgNPs was achieved using cross-flows of 0.5 and 0.7 mL/min with 1 mM NaCl and 0.05% v/v Mucasol as carrier solutions. We also researched the behavior of AgNPs in natural waters using the hyphenated AF4-UV/Vis-DLS system, under determined optimal conditions.

Description: This encyclopedia is a reference for aquatic physical, biological and biogeochemical sciences, collecting and connecting a number of topics, concepts and facts about aquatic systems and their scientific investigation. The scope of the book comprises the aquatic physiosphere-biosphere transition zone, an entity that encompasses both inanimate matter and collectives (the physiosphere) as well as living organisms and collectives (the biosphere). This combined approach is meaningful because both realms are intimately linked and because available methods of investigation are often similar. Much can be gained from considering both spheres at and across their interface jointly, and while there is a strong focus on marine systems, most concepts presented are also applicable to freshwater systems. This presented snapshot of knowledge of the transition zone between the aquatic physiosphere and biosphere is taken from a very specific angle: the point of view of a modeler. Modeling is not only a state-of-the-art mode of scientific investigation, but also requires the explicit specification of all assumptions (helping to avoid fallacies), and offers the advantage of being quantitative and allows for theoretical "what if" scenarios. As in any reference work, equal emphasis is given to fundamental facts, the definition of terms and the explanation of concepts, in an attempt to establish a joint language for physicists, biologists and biogeochemists. Although originating from a modeler's approach to nature, the resulting suite of compatible concepts may also be useful beyond modeling purposes. Furthermore, the material is presented in a condensed, straightforward way. Hence, the length of each entry is limited to one (occasionally two) pages, thus offering a quick introductory overview. This excludes lengthy derivations and very specialized details. The book is geared towards researchers, teachers and advanced students in the field of aquatic (marine and limnic) sciences, in particular those interested or involved in interdisciplinary work.

Description: There is a long-standing debate on how the El Niño/Southern Oscillation (ENSO) amplitude may change during the twenty-first century in response to global warming. Here we identify the sources of uncertainty in the ENSO amplitude projections in models participating in the Coupled Model Intercomparison Phase 5 (CMIP5) and Phase 6 (CMIP6), and quantify scenario uncertainty, model uncertainty and uncertainty due to internal variability. The model projections exhibit a large spread, ranging from increasing standard deviation of up to 0.6 °C to diminishing standard deviation of up to − 0.4 °C by the end of the twenty-first century. The ensemble-mean ENSO amplitude change is close to zero. Internal variability is the main contributor to the uncertainty during the first three decades; model uncertainty dominates thereafter, while scenario uncertainty is relatively small throughout the twenty-first century. The total uncertainty increases from CMIP5 to CMIP6: while model uncertainty is reduced, scenario uncertainty is considerably increased. The models with “realistic” ENSO dynamics have been analyzed separately and categorized into models with too small, moderate and too large ENSO amplitude in comparison to instrumental observations. The smallest uncertainties are observed in the sub-ensemble exhibiting realistic ENSO dynamics and moderate ENSO amplitude. However, the global warming signal in ENSO-amplitude change is undetectable in all sub-ensembles. The zonal wind-SST feedback is identified as an important factor determining ENSO amplitude change: global warming signal in ENSO amplitude and zonal wind-SST feedback strength are highly correlated across the CMIP5 and CMIP6 models.

Description: Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multimillion-year decrease (i.e., Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating hydrographic restriction via South Atlantic–Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological timescales (〉 1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (〉 6 %) and relatively low (∼ 3.5 %) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive intervals before and after OAE 1a. In all cases, these short-term intervals of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we show that fluctuations in weathering-derived nutrient input from the southern African continent, linked to changes in orbitally driven humidity and aridity, were the underlying drivers of repetitive episodes of enhanced organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short-term fluctuations in riverine nutrient fluxes. We explain this relationship using the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.

Description: Over recent decades Antarctic sea-ice extent has increased, alongside widespread ice shelf thinning and freshening of waters along the Antarctic margin. In contrast, Earth system models generally simulate a decrease in sea ice. Circulation of water masses beneath large-cavity ice shelves is not included in current Earth System models and may be a driver of this phenomena. We examine a Holocene sediment core off East Antarctica that records the Neoglacial transition, the last major baseline shift of Antarctic sea ice, and part of a late-Holocene global cooling trend. We provide a multi-proxy record of Holocene glacial meltwater input, sediment transport, and sea-ice variability. Our record, supported by high-resolution ocean modelling, shows that a rapid Antarctic sea-ice increase during the mid-Holocene (∼ 4.5 ka) occurred against a backdrop of increasing glacial meltwater input and gradual climate warming. We suggest that mid-Holocene ice shelf cavity expansion led to cooling of surface waters and sea-ice growth that slowed basal ice shelf melting. Incorporating this feedback mechanism into global climate models will be important for future projections of Antarctic changes.

Description: Biodiversity is declining on a global scale and the spread of invasive alien species (IAS) is a major driver, particularly through predatory impacts. Thus, effective means of assessing and predicting the consequences of IAS predation on native prey population stability remains a vital goal for conservation. Here, we applied two classic ecological concepts, consumer functional response (FR) and prey switching, to predict and understand the ecological impacts of juveniles of the lionfish (Pterois volitans), a notorious and widespread marine invader. Functional responses and prey switching propensities were quantified towards three representative prey species: Artemia salina, Palaemonetes varians, and Gammarus oceanicus. Lionfish exhibited potentially destabilising Type II FRs towards individual prey species, owing to high consumption rates at low prey densities, whilst FR magnitudes differed among prey species. Functional response attack rates (a) were highest, and handling times (h) lowest, towards A. salina, followed by P. varians and then G. oceanicus. Maximum feeding rates (1/h) and functional response ratios (FRR; a/h) also followed this impact gradient for the three prey species. Lionfish, however, displayed a potentially population stabilising prey switching propensity (i.e. frequency-dependent predation) when multiple prey species were presented simultaneously, where disproportionately less of rare prey, and more of abundant prey, were consumed. Whilst FR and FRR magnitudes indicate marked per capita lionfish predatory impacts towards prey species, a strong prey switching propensity may reduce in-field impacts by offering low density prey refuge in biodiverse communities. Our results thus corroborate field patterns documenting variable impacts of lionfish, with prey extirpations less likely in diverse communities owing to frequency-dependent predation.

Description: Global freshwaters are increasingly threatened by pollutants emanating from human activities around watersheds. Microplastic pollution is an increasing problem for rivers worldwide, potentially threatening ecological integrity, ecosystem services and human health. We present quantifications and characterisations of sediment microplastic pollution in a subtropical river system in southern Africa, and relate distributions to wastewater treatment works, abiotic variables and urban environments. We additionally apply several diversity indices to decipher how microplastic types differ across the river system seasonally. Over two thousand microplastic particles were found across five sites and three seasons in the river system, comprising microbeads of various colours and microfibres. Microplastic concentrations were highest and most diverse in the hot–wet (mean range 76.0 ± 10.0–285.5 ± 44.5 microplastic kg−1) season as compared to the cool–dry (16.5 ± 4.5–27.0 ± 5.0 microplastic kg−1) and hot–dry (13.0 ± 4.0–29.0 ± 10.0 microplastic kg−1) seasons, and were mostly dominated by microfibres. However, no clear patterns were found in relation to wastewater treatment operations spatially, or in relation to abiotic variables in the river system. This study therefore finds a diverse range of microplastic types widely distributed in the river system that differ across seasons. Our results provide important, novel insights into plastic pollution in an understudied area of the Global South, and point to extensive pollution from sources outside of wastewater treatment works.

Description: A novel Actinobacterium strain YIM 131861 T, was isolated from lichen collected from the South Bank Forest of the Baltic Sea, Germany. It was Gram-stain-positive, strictly aerobic, catalase positive and oxidase negative, yellow pigmented. Cells were motile with a polar flagellum, irregular rod shaped and did not display spore formation. The strain grew at 15 − 30 °C (optimum 25 °C), at pH 6.0 − 10.0 (optimum pH 7.0) and in the presence of 0 − 1.5% (w/v) NaCl (optimum 1%). Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain YIM 131861 T belonged to the genus Glaciibacter, and exhibited a high sequence similarity (96.4%) with Glaciibacter superstes NBRC 104264 T. The genomic DNA G + C content of strain YIM 131861 T was 68.2 mol%. The average nucleotide identity (ANI) and digital DNA–DNA hybridization (dDDH) values between strain YIM 131861 T and Glaciibacter superstes NBRC 104264 T were 73.2 and 19.9% based on the draft genome sequence. The cell-wall peptidoglycan type was B2γ and contained the 2, 4-diaminobutyric acid as the diagnostic amino acid. Whole cell sugars were galactose, rhamnose, ribose and glucose. It contained MK-12 and MK-13 as the predominant menaquinones. The major cellular fatty acids (〉 10%) were identified as anteiso-C15:0, iso-C16:0 and anteiso-C17:0. The polar lipids were diphosphatidylglycerol, phosphatidylglycerol and two unknown glycolipids. Based on the results of phenotypic, chemotaxonomic and phylogenetic analyses, strain YIM 131861 T should belong to the genus Glaciibacter and represents a novel species of the genus Glaciibacter, for which the name Glaciibacter flavus sp. nov. is proposed. The type strain is YIM 131861 T (= CGMCC 1.16588 T = NBRC 113572 T).

Description: The deep chlorophyll maximum (DCM) is a ubiquitous feature of phytoplankton vertical distribution in stratified waters that is relevant to our understanding of the mechanisms that underpin the variability in photoautotroph ecophysiology across environmental gradients and has implications for remote sensing of aquatic productivity. During the PEACETIME (Process studies at the air-sea interface after dust deposition in the Mediterranean Sea) cruise, carried out from 10 May to 11 June 2017, we obtained 23 concurrent vertical profiles of phytoplankton chlorophyll a, carbon biomass and primary production, as well as heterotrophic prokaryotic production, in the western and central Mediterranean basins. Our main aims were to quantify the relative role of photoacclimation and enhanced growth as underlying mechanisms of the DCM and to assess the trophic coupling between phytoplankton and heterotrophic prokaryotic production. We found that the DCM coincided with a maximum in both the biomass and primary production but not in the growth rate of phytoplankton, which averaged 0.3 d−1 and was relatively constant across the euphotic layer. Photoacclimation explained most of the increased chlorophyll a at the DCM, as the ratio of carbon to chlorophyll a (C:Chl a) decreased from ca. 90–100 (g:g) at the surface to 20–30 at the base of the euphotic layer, while phytoplankton carbon biomass increased from ca. 6 mg C m−3 at the surface to 10–15 mg C m−3 at the DCM. As a result of photoacclimation, there was an uncoupling between chlorophyll a-specific and carbon-specific productivity across the euphotic layer. The ratio of fucoxanthin to total chlorophyll a increased markedly with depth, suggesting an increased contribution of diatoms at the DCM. The increased biomass and carbon fixation at the base of the euphotic zone was associated with enhanced rates of heterotrophic prokaryotic activity, which also showed a surface peak linked with warmer temperatures. Considering the phytoplankton biomass and turnover rates measured at the DCM, nutrient diffusive fluxes across the nutricline were able to supply only a minor fraction of the photoautotroph nitrogen and phosphorus requirements. Thus the deep maxima in biomass and primary production were not fuelled by new nutrients but likely resulted from cell sinking from the upper layers in combination with the high photosynthetic efficiency of a diatom-rich, low-light acclimated community largely sustained by regenerated nutrients. Further studies with increased temporal and spatial resolution will be required to ascertain if the peaks of deep primary production associated with the DCM persist across the western and central Mediterranean Sea throughout the stratification season.

Description: The late-tectonic 511.4 ± 0.6 Ma-old Nomatsaus intrusion (Donkerhoek batholith, Damara orogen, Namibia) consists of moderately peraluminous, magnesian, calc-alkalic to calcic granites similar to I-type granites worldwide. Major and trace-element variations and LREE and HREE concentrations in evolved rocks imply that the fractionated mineral assemblage includes biotite, Fe–Ti oxides, zircon, plagioclase and monazite. Increasing K2O abundance with increasing SiO2 suggests accumulation of K-feldspar; compatible with a small positive Eu anomaly in the most evolved rocks. In comparison with experimental data, the Nomatsaus granite was likely generated from meta-igneous sources of possibly dacitic composition that melted under water-undersaturated conditions (X H2O: 0.25–0.50) and at temperatures between 800 and 850 °C, compatible with the zircon and monazite saturation temperatures of 812 and 852 °C, respectively. The Nomatsaus granite has moderately radiogenic initial 87Sr/86Sr ratios (0.7067–0.7082), relatively radiogenic initial εNd values (− 2.9 to − 4.8) and moderately evolved Pb isotope ratios. Although initial Sr and Nd isotopic compositions of the granite do not vary with SiO2 or MgO contents, fSm/Nd and initial εNd values are negatively correlated indicating limited assimilation of crustal components during monazite-dominated fractional crystallization. The preferred petrogenetic model for the generation of the Nomatsaus granite involves a continent–continent collisional setting with stacking of crustal slices that in combination with high radioactive heat production rates heated the thickened crust, leading to the medium-P/high-T environment characteristic of the southern Central Zone of the Damara orogen. Such a setting promoted partial melting of metasedimentary sources during the initial stages of crustal heating, followed by the partial melting of meta-igneous rocks at mid-crustal levels at higher P–T conditions and relatively late in the orogenic evolution.

Description: Regional anomalies of steric sea level are either due to redistribution of heat and freshwater anomalies or due to ocean-atmosphere buoyancy fluxes. Interannual to decadal variability in sea level across the tropical Pacific is mainly due to steric variations driven by wind stress anomalies. The importance of air--sea buoyancy fluxes is less clear. We use a global, eddy permitting ocean model and a series of sensitivity experiments with quasi-climatological momentum and buoyancy fluxes to identify the contribution of buoyancy fluxes for interannual to decadal sea level variability in the tropical Pacific. We find their contribution on interannual timescales to be strongest in the central tropical Pacific at around 10° latitude in both hemispheres and also relevant in the very east of the tropical domain. Buoyancy flux forced anomalies are in phase with variations driven by wind stress changes but their effect on the prevailing anomalies and the importance of heat and fresh water fluxes vary locally. In the eastern tropical basin interannual sea level variability is amplified by anomalous heat fluxes, while the importance of fresh water fluxes is small and neither has any impact on decadal timescales. In the western tropical Pacific the variability on interannual and decadal timescales is dampened by both, heat and freshwater fluxes. The mechanism involves westward propagating Rossby waves that are triggered during ENSO events by anomalous buoyancy fluxes in the central tropical Pacific and counteract the prevailing sea level anomalies once they reach the western part of the basin.

Description: The inflow of relatively warm and salty water from the Indian Ocean into the South Atlantic via Agulhas leakage is important for the global overturning circulation and the global climate. In this study, we analyse the robustness of Agulhas leakage estimates as well as the thermohaline property modifications of Agulhas leakage south of Africa. Lagrangian experiments with both the newly developed tool Parcels and the well established tool Ariane were performed to simulate Agulhas leakage in the eddy-rich ocean–sea-ice model INALT20 (1/20∘ horizontal resolution) forced by the JRA55-do atmospheric boundary conditions. The average transport, its variability, trend and the transit time from the Agulhas Current to the Cape Basin of Agulhas leakage is simulated comparably with both Lagrangian tools, emphasizing the robustness of our method. Different designs of the Lagrangian experiment alter in particular the total transport of Agulhas leakage by up to 2 Sv, but the variability and trend of the transport are similar across these estimates. During the transit from the Agulhas Current at 32∘ S to the Cape Basin, a cooling and freshening of Agulhas leakage waters occurs especially at the location of the Agulhas Retroflection, resulting in a density increase as the thermal effect dominates. Beyond the strong air–sea exchange around South Africa, Agulhas leakage warms and salinifies the water masses below the thermocline in the South Atlantic.

Description: LivingAmmoniaspecies and an inventory of dead assemblages from Adriatic subtidal, nearshoreenvironments were investigated at four stations off Bellaria, Italy.Ammonia falsobeccarii, Ammonia parkinso-niana, Ammonia tepida, andAmmonia venetawere recognized in the living (rose-bengal-stained) fauna, andAmmonia bellarian. sp. is described herein for the first time.Ammonia beccariiwas only found in the deadassemblage. The biometry of 368 living individuals was analysed by using light microscopic and scanning elec-tron microscopic images of three aspects. A total of 15 numerical and 8 qualitative parameters were measuredand assessed, 5 of which were recognized to be prone to a certain subjectivity of the observer. The accuracy ofnumerical data as revealed by the mean residuals of parallel measurements by different observers ranged from0.5 % to 5.5 %. The results indicated a high degree of intraspecific variability. The test sizes of the individualspecies were log-normally distributed and varied among the stations. Parameters not related to the growth ofthe individuals, i.e. flatness of the tests, dimensions of the second-youngest chamber, proloculus, umbilical andpore diameter, sinistral–dextral coiling, and umbilical boss size, were recognized as being species-distinctive incombination. They may well supplement qualitative criteria that were commonly used for species discriminationsuch as a lobate outline, a subacute or rounded peripheral margin, or the degree of ornamentation on the spiraland umbilical sides. The averages of the measured parameters were often lower than the range of previouslypublished values, mainly because the latter were retrieved from a few adult specimens and not from the wholeassemblage as in the present approach. We conclude that the unprecedented high proportions ofAmmonia bec-cariiin the northern Adriatic may well be artificial. A robust species identification without genetic analyses ispossible by considering designated biometric parameters. This approach is also applicable to earlier literaturedata, and their re-assessment is critical for a correct denomination of recent genotypes.

Description: Improving the understanding of changes in the mean and variability of climate variables as well as their interrelation is crucial for reliable climate change projections. Comparisons between general circulation models and paleoclimate archives using indirect proxies for temperature or precipitation have been used to test and validate the capability of climate models to represent climate changes. The oxygen isotopic ratio δ18O, a proxy for many different climate variables, is routinely measured in speleothem samples at decadal or higher resolution, and single specimens can cover full glacial–interglacial cycles. The calcium carbonate cave deposits are precisely dateable and provide well preserved (semi-)continuous albeit multivariate climate signals in the lower and mid-latitudes, where the measured δ18O in the mineral does not directly represent temperature or precipitation. Therefore, speleothems represent suitable archives to assess climate model abilities to simulate climate variability beyond the timescales covered by meteorological observations (101–102 years). Here, we present three transient isotope-enabled simulations from the Hadley Center Climate Model version 3 (iHadCM3) covering the last millennium (850–1850 CE) and compare them to a large global dataset of speleothem δ18O records from the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database version 2 (Comas-Bru et al., 2020b). We systematically evaluate offsets in mean and variance of simulated δ18O and test for the main climate drivers recorded in δ18O for individual records or regions. The time-mean spatial offsets between the simulated δ18O and the speleothem data are fairly small. However, using robust filters and spectral analysis, we show that the observed archive-based variability of δ18O is lower than simulated by iHadCM3 on decadal and higher on centennial timescales. Most of this difference can likely be attributed to the records' lower temporal resolution and averaging or smoothing processes affecting the δ18O signal, e.g., through soil water residence times. Using cross-correlation analyses at site level and modeled grid-box level, we find evidence for highly variable but generally low signal-to-noise ratios in the proxy data. This points to a high influence of cave-internal processes and regional climate particularities and could suggest low regional representativity of individual sites. Long-range strong positive correlations dominate the speleothem correlation network but are much weaker in the simulation. One reason for this could lie in a lack of long-term internal climate variability in these model simulations, which could be tested by repeating similar comparisons with other isotope-enabled climate models and paleoclimate databases.

Description: Karst groundwater systems are characterized by the presence of multiple porosity types. Of these, subsurface conduits that facilitate concentrated, heterogeneous flow are challenging to resolve geologically and geophysically. This is especially the case in evaporite karst systems, such as those present on the shores of the Dead Sea, where rapid geomorphological changes are linked to a fall in base level by over 35 m since 1967. Here we combine field observations, remote-sensing analysis, and multiple geophysical surveying methods (shear wave reflection seismics, electrical resistivity tomography, ERT, self-potential, SP, and ground-penetrating radar, GPR) to investigate the nature of subsurface groundwater flow and its interaction with hypersaline Dead Sea water on the rapidly retreating eastern shoreline, near Ghor Al-Haditha in Jordan. Remote-sensing data highlight links between the evolution of surface stream channels fed by groundwater springs and the development of surface subsidence patterns over a 25-year period. ERT and SP data from the head of one groundwater-fed channel adjacent to the former lakeshore show anomalies that point to concentrated, multidirectional water flow in conduits located in the shallow subsurface (〈 25 m depth). ERT surveys further inland show anomalies that are coincident with the axis of a major depression and that we interpret as representing subsurface water flow. Low-frequency GPR surveys reveal the limit between unsaturated and saturated zones (〈 30 m depth) surrounding the main depression area. Shear wave seismic reflection data nearly 1 km further inland reveal buried paleochannels within alluvial fan deposits, which we interpret as pathways for groundwater flow from the main wadi in the area towards the springs feeding the surface streams. Finally, simulations of density-driven flow of hypersaline and undersaturated groundwaters in response to base-level fall perform realistically if they include the generation of karst conduits near the shoreline. The combined approaches lead to a refined conceptual model of the hydrological and geomorphological processes developed at this part of the Dead Sea, whereby matrix flow through the superficial aquifer inland transitions to conduit flow nearer the shore where evaporite deposits are encountered. These conduits play a key role in the development of springs, stream channels and subsidence across the study area.