• European Marine Science
• Other research products close
• CH close

cti_cross_region_taxon This repository provides the data and codes used in the manuscript: Guillem Chust, Ernesto Villarino, Matthew McLean, Nova Mieszkowska, Lisandro Benedetti-Cecchi, Fabio Bulleri, Chiara Ravaglioli, Angel Borja, Iñigo Muxika, José A. Fernandes-Salvador, Leire Ibaibarriaga, Ainhize Uriarte, Marta Revilla, Fernando Villate, Arantza Iriarte, Ibon Uriarte, Soultana Zervoudaki, Jacob Carstensen, Paul J. Somerfield, Ana M. Queirós , Andrea J. McEvoy, Arnaud Auber, Manuel Hidalgo, Marta Coll, Joaquim Garrabou, Daniel Gómez-Gras, Cristina Linares, Francisco Ramírez, Núria Margarit, Mario Lepage, Chloé Dambrine, Jérémy Lobry, Myron A. Peck, Paula de la Barra, Anieke van Leeuwen, Gil Rilov, Erez Yeruham, Anik Brind'Amour, and Martin Lindegren. Cross-basin and cross-taxa patterns of marine community tropicalization and deborealization in warming European seas. Nature Communications Please check the github repo cti_cross_region_taxon for updates. Abstract Ocean warming and acidification, decreases in dissolved oxygen concentrations, and changes in primary production are causing an unprecedented global redistribution of marine life. The identification of underlying ecological processes underpinning marine species turnover, particularly the prevalence of tropicalization over deborealization, has been recently debated in the context of ocean warming. Here, we track changes in the mean thermal affinity of marine communities across European seas by calculating the Community Temperature Index for 65 biodiversity time series collected over four decades and containing 1,817 species from different communities (zooplankton, coastal benthos, pelagic and demersal invertebrates and fish). We show most communities and sites have clearly responded to ongoing ocean warming via abundance increases of warm-water species (tropicalization; 54%) and decreases of cold-water species (deborealization; 18%). Tropicalization dominated Atlantic sites compared to semi-enclosed basins such as the Mediterranean and Baltic Seas, probably due to physical barrier constraints to connectivity and species colonization. Semi-enclosed basins appeared to be particularly vulnerable to ocean warming, experiencing the fastest rates of warming and biodiversity loss through deborealization. Keywords: climate change, community temperature index, CTI, ocean connectivity. Acknowledgements This study has been supported by the European Union's Horizon 2020 research and innovation programme under grant agreement No 869300 (FutureMARES project) (G.C., E.V., J.F-S., M.P., M.Lin., C.D., D.G-G., G.R., L.B., C.R., M.C., F.R., G.R., P.B., M.P., M.Lep., J.L., A.B., N.M., J.G., F.B., L.B.C, A.Q., F.V., A.I., and I.U.), and by the Urban Klima 2050 -- LIFE 18 IPC 000001 project, which has been received funding from European Union's LIFE programme (G.C., E.V., L.I., A.B., A.U., and M.R.). Additional financial support was obtained from the Basque Government (PIBA2020-1-0028 & IT1723-22). We thank the NOAA Climate Prediction Center for providing Sea Temperature data through the NCEP Global Ocean Data Assimilation System (GODAS) www.cpc.ncep.noaa.gov/products/GODAS. We also thank Ocean Biodiversity Information System (OBIS) for providing global occurrences of the biological group studied here. Data from the Basque Country were obtained from the Basque Water Agency (URA) monitoring network, through a Convention with AZTI. M.C., J.G., D.G.-G. and F.R. acknowledge the 'Severo Ochoa Centre of Excellence' accreditation (CEX2019-000928-S). Authors M.H. and M.Lin. are grateful for the support from ICES Working Group on Comparative Ecosystem-based Analyses of Atlantic and Mediterranean marine systems (WGCOMEDA) for this research. This paper is contribution nº xxxx from AZTI, Marine Research, Basque Research and Technology Alliance (BRTA)

The uploaded file contains three separate Python scripts aimed at modeling the sinking velocity of natural Trichodesmium colonies loaded with dust particles. Two scripts have been developed to calculate the dust influencing factor (K), and one script is dedicated to calculating the sinking velocity of particle-free Trichodesmium colonies (v0). All Python scripts are available on GitHub (https://github.com/Zhanzhu1110/Trichobuoyancy.git) and can be processed using the "Codespaces" function on GitHub. The file 'K given values' is a specialized script designed to simulate the dust factor (K), assuming seawater and dust densities of 1028 kg/m3 and 2500 kg/m3, respectively, with a colony radius of 1 mm.

This study presents simulations of Greenland surface melt for the Eemian interglacial period (∼130 000 to 115 000 years ago) derived from regional climate simulations with a coupled surface energy balance model. Surface melt is of high relevance due to its potential effect on ice core observations, e.g., lowering the preserved total air content (TAC) used to infer past surface elevation. An investigation of surface melt is particularly interesting for warm periods with high surface melt, such as the Eemian interglacial period. Furthermore, Eemian ice is the deepest and most compressed ice preserved on Greenland, resulting in our inability to identify melt layers visually. Therefore, simulating Eemian melt rates and associated melt layers is beneficial to improve the reconstruction of past surface elevation. Estimated TAC, based on simulated melt during the Eemian, could explain the lower TAC observations. The simulations show Eemian surface melt at all deep Greenland ice core locations and an average of up to ∼30 melt days per year at Dye-3, corresponding to more than 600 mm water equivalent (w.e.) of annual melt. For higher ice sheet locations, between 60 and 150 mmw.e.yr-1 on average are simulated. At the summit of Greenland, this yields a refreezing ratio of more than 25 % of the annual accumulation. As a consequence, high melt rates during warm periods should be considered when interpreting Greenland TAC fluctuations as surface elevation changes. In addition to estimating the influence of melt on past TAC in ice cores, the simulated surface melt could potentially be used to identify coring locations where Greenland ice is best preserved.

As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated pre-study 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 factor 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.

Coastal systems can act as filters for anthropogenic nutrient input into marine environments. Here, we assess the processes controlling the removal of phosphorus (P) and nitrogen (N) for four sites in the eutrophic Stockholm archipelago. Bottom water concentrations of oxygen (O2) and P are inversely correlated. This is attributed to the seasonal release of P from iron-oxide-bound (Fe-oxide-bound) P in surface sediments and from degrading organic matter. The abundant presence of sulfide in the pore water and its high upward flux towards the sediment surface (∼4 to 8 mmol m−2 d−1), linked to prior deposition of organic-rich sediments in a low-O2 setting (“legacy of hypoxia”), hinder the formation of a larger Fe-oxide-bound P pool in winter. This is most pronounced at sites where water column mixing is naturally relatively low and where low bottom water O2 concentrations prevail in summer. Burial rates of P are high at all sites (0.03–0.3 mol m−2 yr−1), a combined result of high sedimentation rates (0.5 to 3.5 cm yr−1) and high sedimentary P at depth (∼30 to 50 µmol g−1). Sedimentary P is dominated by Fe-bound P and organic P at the sediment surface and by organic P, authigenic Ca-P and detrital P at depth. Apart from one site in the inner archipelago, where a vivianite-type Fe(II)-P mineral is likely present at depth, there is little evidence for sink switching of organic or Fe-oxide-bound P to authigenic P minerals. Denitrification is the major benthic nitrate-reducing process at all sites (0.09 to 1.7 mmol m−2 d−1) with rates decreasing seaward from the inner to outer archipelago. Our results explain how sediments in this eutrophic coastal system can remove P through burial at a relatively high rate, regardless of whether the bottom waters are oxic or (frequently) hypoxic. Our results suggest that benthic N processes undergo annual cycles of removal and recycling in response to hypoxic conditions. Further nutrient load reductions are expected to contribute to the recovery of the eutrophic Stockholm archipelago from hypoxia. Based on the dominant pathways of P and N removal identified in this study, it is expected that the sediments will continue to remove part of the P and N loads.

Concentrations of dissolved 230Th in the ocean water column increase with depth due to scavenging and downward particle flux. Due to the 230Th scavenging process, any change in the calcium carbonate (CaCO3) fraction of the marine particle flux due to changes in biological CaCO3 hard-shell production as a consequence of progressing ocean acidification would be reflected in the dissolved 230Th activity. Our prognostic simulations with a biogeochemical ocean general circulation model using different scenarios for the reduction of CaCO3 production under ocean acidification and different greenhouse gas emission scenarios – the Representative Concentration Pathways (RCPs) 8.5 to 2.6 – reveal the potential for deep 230Th measurements to detect reduced CaCO3 production at the sea surface. The time of emergence of an acidification-induced signal on dissolved 230Th is of the same order of magnitude as for alkalinity measurements. Interannual and decadal variability in factors other than a reduction in CaCO3 hard-shell production may mask the ocean-acidification-induced signal in dissolved 230Th and make detection of the pure CaCO3-induced signal more difficult so that only really strong changes in marine CaCO3 export would be unambiguously identifiable soon. Nevertheless, the impacts of changes in CaCO3 export production on marine 230Th are stronger than those for changes in POC (particulate organic carbon) or clay fluxes.

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Version 2 of the Global Ocean Data Analysis Project (GLODAPv2) data product is composed of data from 724 scientific cruises covering the global ocean. It includes data assembled during the previous efforts GLODAPv1.1 (Global Ocean Data Analysis Project version 1.1) in 2004, CARINA (CARbon IN the Atlantic) in 2009/2010, and PACIFICA (PACIFic ocean Interior CArbon) in 2013, as well as data from an additional 168 cruises. Data for 12 core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have been subjected to extensive quality control, including systematic evaluation of bias. The data are available in two formats: (i) as submitted but updated to WOCE exchange format and (ii) as a merged and internally consistent data product. In the latter, adjustments have been applied to remove significant biases, respecting occurrences of any known or likely time trends or variations. Adjustments applied by previous efforts were re-evaluated. Hence, GLODAPv2 is not a simple merging of previous products with some new data added but a unique, internally consistent data product. This compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 6 µmol kg−1 in total alkalinity, 0.005 in pH, and 5 % for the halogenated transient tracers.The original data and their documentation and doi codes are available at the Carbon Dioxide Information Analysis Center (http://cdiac.ornl.gov/oceans/GLODAPv2/). This site also provides access to the calibrated data product, which is provided as a single global file or four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under the doi:10.3334/CDIAC/OTG.NDP093_GLODAPv2. The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2 methods and products and includes a broad overview of the secondary quality control results. The magnitude of and reasoning behind each adjustment is available on a per-cruise and per-variable basis in the online Adjustment Table.