Description: We have developed a global biogeographic classification of the mesopelagic zone to reflect the regional scales over which the ocean interior varies in terms of biodiversity and function. An integrated approach was necessary, as global gaps in information and variable sampling methods preclude strictly statistical approaches. A panel combining expertise in oceanography, geospatial mapping, and deep-sea biology convened to collate expert opinion on the distributional patterns of pelagic fauna relative to environmental proxies (temperature, salinity, and dissolved oxygen at mesopelagic depths). An iterative Delphi Method integrating additional biological and physical data was used to classify biogeographic ecoregions and to identify the location of ecoregion boundaries or inter-regions gradients. We define 33 global mesopelagic ecoregions. Of these, 20 are oceanic while 13 are ‘distant neritic.’ While each is driven by a complex of controlling factors, the putative primary driver of each ecoregion was identified. While work remains to be done to produce a comprehensive and robust mesopelagic biogeography (i.e., reflecting temporal variation), we believe that the classification set forth in this study will prove to be a useful and timely input to policy planning and management for conservation of deep-pelagic marine resources. In particular, it gives an indication of the spatial scale at which faunal communities are expected to be broadly similar in composition, and hence can inform application of ecosystem-based management approaches, marine spatial planning and the distribution and spacing of networks of representative protected areas

Description: Microbial degradation of organic carbon in marine sediments is a key driver of global element cycles on multiple time scales. However, it is not known to what depth microorganisms alter organic carbon in marine sediments or how microbial rates of organic carbon processing change with depth, and thus time since burial, on a global scale. To better understand the connection between the dynamic carbon cycle and life’s limits in the deep subsurface, we have combined a number of global data sets with a reaction transport model (RTM) describing first, organic carbon degradation in marine sediments deposited throughout the Quaternary Period and second, a bioenergetic model for microbial activity. The RTM is applied globally, recognizing three distinct depositional environments – continental shelf, margin and abyssal zones. The results include the masses of particulate organic carbon, POC, stored in three sediment-depth layers: bioturbated Holocene (1.7 × 10^17 g C), non-bioturbated Holocene (2.6 × 10^18 g C) and Pleistocene (1.4 × 1020 g C) sediments. The global depth-integrated rates of POC degradation have been determined to be 6.8 × 10^13, 1.2 × 10^14 and 1.2 × 10^14 g C yr-1 for the same three layers, respectively. A number of maps depicting the distribution of POC, as well as the fraction that has been degraded have also been generated. Using POC degradation as a proxy for microbial catabolic activity, total heterotrophic processing of POC throughout the Quaternary is estimated to be between 10^-11 – 10^-6 g C cm-3 yr-1, depending on the time since deposition and location. Bioenergetic modeling reveals that laboratory-determined microbial maintenance powers are poor predictors of sediment biomass concentration, but that cell concentrations in marine sediments can be accurately predicted by combining bioenergetic models with the rates of POC degradation determined in this study. Our model can be used to quantitatively describe both the carbon cycle and microbial activity on a global scale for marine sediments less than 2.59 million years old.

Description: Riverine nutrient export is an important process in marine coastal biogeochemistry and also impacts global marine biology. The nitrogen cycle is a key player here. Internal feedbacks regulate not only nitrogen distribution, but also primary production and thereby oxygen concentrations. Phosphorus is another essential nutrient and interacts with the nitrogen cycle via different feedback mechanisms. After a previous study of the marine nitrogen cycle response to riverine nitrogen supply, we here additionally include phosphorus from river export with different phosphorus burial scenarios and study the impact of phosphorus alone and in combination with nitrogen in a global 3-D ocean biogeochemistry model. Again, we analyse the effects on near coastal and open ocean biogeochemistry. We find that the addition of bio-available riverine phosphorus alone or together with nitrogen affects marine biology on millennial timescales more than riverine nitrogen alone. Biogeochemical feedbacks in the marine nitrogen cycle are strongly influenced by the additional phosphorus. Where bio-available phosphorus is increased by river input, nitrogen concentrations increase as well, except for regions with high denitrification rates. High phosphorus burial rates decrease biological production significantly. Globally, riverine phosphorus leads to elevated primary production rates in the coastal and open oceans.

Description: The accumulation of gas hydrates in marine sediments is essentially controlled by the accumulation of particulate organic carbon (POC) which is microbially converted into methane, the thickness of the gas hydrate stability zone (GHSZ) where methane can be trapped, the sedimentation rate (SR) that controls the time that POC and the generated methane stays within the GHSZ, and the delivery of methane from deep-seated sediments by ascending pore fluids and gas into the GHSZ. Recently, Wallmann et al. (2012) presented transfer functions to predict the gas hydrate inventory in diffusion-controlled geological systems based on SR, POC and GHSZ thickness for two different scenarios: normal and full compacting sediments. We apply these functions to global data sets of bathymetry, heat flow, seafloor temperature, POC input and SR, estimating a global mass of carbon stored in marine methane hydrates from 3 to 455 Gt of carbon (GtC) depending on the sedimentation and compaction conditions. The global sediment volume of the GHSZ in continental margins is estimated to be 60–67 × 1015 m3, with a total of 7 × 1015 m3 of pore volume (available for GH accumulation). However, seepage of methane-rich fluids is known to have a pronounced effect on gas hydrate accumulation. Therefore, we carried out a set of systematic model runs with the transport-reaction code in order to derive an extended transfer function explicitly considering upward fluid advection. Using averaged fluid velocities for active margins, which were derived from mass balance considerations, this extended transfer function predicts the enhanced gas hydrate accumulation along the continental margins worldwide. Different scenarios were investigated resulting in a global mass of sub-seafloor gas hydrates of ~ 550 GtC. Overall, our systematic approach allows to clearly and quantitatively distinguish between the effect of biogenic methane generation from POC and fluid advection on the accumulation of gas hydrate, and hence, provides a simple prognostic tool for the estimation of large-scale and global gas hydrate inventories in marine sediments.

Description: Spatial predictions of total organic carbon (TOC) concentrations and stocks are crucial for understanding marine sediments’ role as a significant carbon sink in the global carbon cycle. In this study, we present a geospatial prediction of TOC concentrations and stocks at a 5 x 5 arc minute grid scale, using a deep learning model — a novel machine learning approach based on a new compilation of over 22,000 global TOC measurements and a new set of predictors, such as seafloor lithologies, grain size distribution, and an alpha-chlorophyll satellite data. In our study, we compared the predictions and discuss the limitations from various machine learning methods. Our findings reveal that the neural network approach outperforms methods such as k Nearest Neighbors and random forests, which tend to overfit to the training data, especially in highly heterogeneous and complex geological settings. We provide estimates of mean TOC concentrations and total carbon stock in both continental shelves and deep sea settings across various marine regions and oceans. Our model suggests that the upper 10 cm of oceanic sediments harbors approximately 171 Pg of TOC stock and has a mean TOC concentration of 0.68 %. Furthermore, we introduce a standardized methodology for quantifying predictive uncertainty using Monte Carlo dropout and present a map of information gain, that measures the expected increase in model knowledge achieved through in-situ sampling at specific locations which is pivotal for sampling strategy planning.

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: Highlights • Foraminifera responded faster to environmental perturbations than macroorganisms. • Population densities doubled and Asterigerinata mamilla bloomed under stronger flow. • Excess of threshold values and changes to new ecosystems induce acme occurrences. Abstract Benthic foraminifera are deemed sensitive indicators of environmental conditions. Triggers and magnitudes of faunal response to environmental changes are yet poorly constrained. Benthic foraminiferal faunas were monitored annually at Ria Formosa (Algarve, Portugal) coastal lagoon since 2013. Distinct environmental changes were recognised during the monitoring period. The relocation of a tidal inlet in winter 2015 effected faster flushing, higher tidal levels, and stronger currents in the Esteiro do Ancão tidal channel. The epibenthic foraminiferal species Asterigerinata mamilla increased in abundance and the population densities of the whole fauna were double as high as before inlet relocation. Enhanced sediment redeposition was recorded and extensive polychaete colonies successively replaced firmground patches with oysters. The standing stock of the foraminiferal fauna declined in the next year due to food impoverishment, while the high hydraulic energy levels and high percentages of Asterigerinata mamilla maintained. Benthic foraminifera responded much faster to environmental perturbations than macroorganisms identifying them as powerful proxies in environmental studies.

Description: Highlights • More diverse non-native taxa generally include more economically costly species. • Chordates, nematodes and pathogens are among significantly over-represented taxa. • Monetary cost magnitude links positively to numbers of costly invasive species. • Costs are biased towards a few ‘hyper-costly’ invasive species groups. • Future invasion rates will continue to harbour new economically costly species. Abstract A dominant syndrome of the Anthropocene is the rapid worldwide spread of invasive species with devastating environmental and socio-economic impacts. However, the dynamics underlying the impacts of biological invasions remain contested. A hypothesis posits that the richness of impactful invasive species increases proportionally with the richness of non-native species more generally. A competing hypothesis suggests that certain species features disproportionately enhance the chances of non-native species becoming impactful, causing invasive species to arise disproportionately relative to the numbers of non-native species. We test whether invasive species with reported monetary costs reflect global numbers of established non-native species among phyla, classes, and families. Our results reveal that numbers of invasive species with economic costs largely reflect non-native species richness among taxa (i.e., in 96 % of families). However, a few costly taxa were over- and under-represented, and their composition differed among environments and regions. Chordates, nematodes, and pathogenic groups tended to be the most over-represented phyla with reported monetary costs, with mammals, insects, fungi, roundworms, and medically-important microorganisms being over-represented classes. Numbers of costly invasive species increased significantly with non-native richness per taxon, while monetary cost magnitudes at the family level were also significantly related to costly invasive species richness. Costs were biased towards a few ‘hyper-costly’ taxa (such as termites, mosquitoes, cats, weevils, rodents, ants, and asters). Ordination analysis revealed significant dissimilarity between non-native and costly invasive taxon assemblages. These results highlight taxonomic groups which harbour disproportionately high numbers of costly invasive species and monetary cost magnitudes. Collectively, our findings support prevention of arrival and containment of spread of non-native species as a whole through effective strategies for mitigation of the rapidly amplifying impacts of invasive species. Yet, the hyper- costly taxa identified here should receive greater focus from managers to reduce impacts of current invasive species.

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr−1, ELUC 0.9 ± 0.5 GtC yr−1, GATM 4.3 ± 0.1 GtC yr−1, SOCEAN 2.5 ± 0.5 GtC yr−1, and SLAND 2.8 ± 0.8 GtC yr−1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr−1, 2.2% above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and assuming an ELUC of 1.0 ± 0.5 GtC yr−1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr−1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70% from EFF (390 ± 20 GtC) and 30% from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3).

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2, and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), EFF was 9.0 ± 0.5 GtC yr−1, ELUC was 0.9 ± 0.5 GtC yr−1, GATM was 4.4 ± 0.1 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 3.0 ± 0.8 GtC yr−1. For the year 2014 alone, EFF grew to 9.8 ± 0.5 GtC yr−1, 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr−1 that took place during 2005–2014. Also, for 2014, ELUC was 1.1 ± 0.5 GtC yr−1, GATM was 3.9 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and SLAND was 4.1 ± 0.9 GtC yr−1. GATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger SLAND for that year. The global atmospheric CO2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in EFF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of EFF and assumed constant ELUC for 2015, cumulative emissions of CO2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO2) for 1870–2015, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2015).