• European Marine Science
• Other research productsclose
• 2014-2023close
• European Commission close
• EC|FP7 close
• GB close
• RO close
• UA close

Tropospheric ozone (O3) concentrations depend on a combination of hemispheric, regional, and local-scale processes. Estimates of how much O3 is produced locally vs. transported from further afield are essential in air quality management and regulatory policies. Here, a tagged-ozone mechanism within the Weather Research and Forecasting model coupled with chemistry (WRF-Chem) is used to quantify the contributions to surface O3 in the UK from anthropogenic nitrogen oxide (NOx) emissions from inside and outside the UK during May–August 2015. The contribution of the different source regions to three regulatory O3 metrics is also examined. It is shown that model simulations predict the concentration and spatial distribution of surface O3 with a domain-wide mean bias of −3.7 ppbv. Anthropogenic NOx emissions from the UK and Europe account for 13 % and 16 %, respectively, of the monthly mean surface O3 in the UK, as the majority (71 %) of O3 originates from the hemispheric background. Hemispheric O3 contributes the most to concentrations in the north and the west of the UK with peaks in May, whereas European and UK contributions are most significant in the east, south-east, and London, i.e. the UK's most populated areas, intensifying towards June and July. Moreover, O3 from European sources is generally transported to the UK rather than produced in situ. It is demonstrated that more stringent emission controls over continental Europe, particularly in western Europe, would be necessary to improve the health-related metric MDA8 O3 above 50 and 60 ppbv. Emission controls over larger areas, such as the Northern Hemisphere, are instead required to lessen the impacts on ecosystems as quantified by the AOT40 metric.

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.

Coccolithophores are globally important marine calcifying phytoplankton that utilize a haplo-diplontic life cycle. The haplo-diplontic life cycle allows coccolithophores to divide in both life cycle phases and potentially expands coccolithophore niche volume. Research has, however, to date largely overlooked the life cycle of coccolithophores and has instead focused on the diploid life cycle phase of coccolithophores. Through the synthesis and analysis of global scanning electron microscopy (SEM) coccolithophore abundance data (n=2534), we find that calcified haploid coccolithophores generally constitute a minor component of the total coccolithophore abundance (≈ 2 %–15 % depending on season). However, using case studies in the Atlantic Ocean and Mediterranean Sea, we show that, depending on environmental conditions, calcifying haploid coccolithophores can be significant contributors to the coccolithophore standing stock (up to ≈30 %). Furthermore, using hypervolumes to quantify the niche of coccolithophores, we illustrate that the haploid and diploid life cycle phases inhabit contrasting niches and that on average this allows coccolithophores to expand their niche by ≈18.8 %, with a range of 3 %–76 % for individual species. Our results highlight that future coccolithophore research should consider both life cycle stages, as omission of the haploid life cycle phase in current research limits our understanding of coccolithophore ecology. Our results furthermore suggest a different response to nutrient limitation and stratification, which may be of relevance for further climate scenarios. Our compilation highlights the spatial and temporal sparsity of SEM measurements and the need for new molecular techniques to identify uncalcified haploid coccolithophores. Our work also emphasizes the need for further work on the carbonate chemistry niche of the coccolithophore life cycle.

In 2013, an ice core was recovered from Roosevelt Island in the Ross Sea, Antarctica, as part of the Roosevelt Island Climate Evolution (RICE) project. Roosevelt Island is located between two submarine troughs carved by paleo-ice-streams. The RICE ice core provides new important information about the past configuration of the West Antarctic Ice Sheet and its retreat during the most recent deglaciation. In this work, we present the RICE17 chronology and discuss preliminary observations from the new records of methane, the isotopic composition of atmospheric molecular oxygen (δ18O-Oatm), the isotopic composition of atmospheric molecular nitrogen (δ15N-N2) and total air content (TAC). RICE17 is a composite chronology combining annual layer interpretations, gas synchronization, and firn modeling strategies in different sections of the core. An automated matching algorithm is developed for synchronizing the high-resolution section of the RICE gas records (60–720 m, 1971 CE to 30 ka) to corresponding records from the WAIS Divide ice core, while deeper sections are manually matched. Ice age for the top 343 m (2635 yr BP, before 1950 C.E.) is derived from annual layer interpretations and described in the accompanying paper by Winstrup et al. (2017). For deeper sections, the RICE17 ice age scale is based on the gas age constraints and the ice age-gas age offset estimated by a firn densification model. Novel aspects of this work include: 1) stratigraphic matching of centennial-scale variations in methane for pre-anthropogenic time periods, a strategy which will be applicable for developing precise chronologies for future ice cores, 2) the observation of centennial-scale variability in methane throughout the Holocene which suggests that similar variations during the late preindustrial period need not be anthropogenic, and 3) the observation of continuous climate records dating back to ∼ 65 ka which provide evidence that the Roosevelt Island Ice Dome was a constant feature throughout the last glacial period.

Much of our understanding of Earth's past climate states comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, major intervals in those records that lack the temporal resolution and/or age control required to identify climate forcing and feedback mechanisms. Here we document 66 million years of global climate by a new high-fidelity Cenozoic global reference benthic carbon and oxygen isotope dataset (CENOGRID). Using recurrence analysis, we find that on timescales of millions of years Earth's climate can be grouped into Hothouse, Warmhouse, Coolhouse and Icehouse states separated by transitions related to changing greenhouse gas levels and the growth of polar ice sheets. Each Cenozoic climate state is paced by orbital cycles, but the response to radiative forcing is state dependent.

The marine silicon cycle is intrinsically linked with carbon cycling in the oceans via biological production of silica by a wide range of organisms. The stable silicon isotopic composition (denoted by δ30Si) of siliceous microfossils extracted from sediment cores can be used as an archive of past oceanic silicon cycling. However, the silicon isotopic composition of biogenic silica has only been measured in diatoms, sponges and radiolarians, and isotopic fractionation relative to seawater is entirely unknown for many other silicifiers. Furthermore, the biochemical pathways and mechanisms that determine isotopic fractionation during biosilicification remain poorly understood. Here, we present the first measurements of the silicon isotopic fractionation during biosilicification by loricate choanoflagellates, a group of protists closely related to animals. We cultured two species of choanoflagellates, Diaphanoeca grandis and Stephanoeca diplocostata, which showed consistently greater isotopic fractionation (approximately −5 ‰ to −7 ‰) than cultured diatoms (−0.5 ‰ to −2.1 ‰). Instead, choanoflagellate silicon isotopic fractionation appears to be more similar to sponges grown under similar dissolved silica concentrations. Our results highlight that there is a taxonomic component to silicon isotope fractionation during biosilicification, possibly via a shared or related biochemical transport pathway. These findings have implications for the use of biogenic silica δ30Si produced by different silicifiers as proxies for past oceanic change.

Bottom trawling in the deep sea is one of the main drivers of sediment resuspension, eroding the seafloor and altering the content and composition of sedimentary organic matter (OM). The physical and biogeochemical impacts of bottom trawling were studied on the continental slope of the Gulf of Castellammare, Sicily (southwestern Mediterranean), through the analysis of two triplicate sediment cores collected at trawled and untrawled sites (∼550 m water depth) during the summer of 2016. Geochemical and sedimentological parameters (excess 210Pb, excess 234Th, 137Cs, dry bulk density, and grain size), elemental (organic carbon and nitrogen) and biochemical composition of sedimentary OM (proteins, carbohydrates, lipids), as well as its freshness (phytopigments) and degradation rates were determined in both coring locations. The untrawled site had a sedimentation rate of 0.15 cm yr−1 and presented a 6 cm thick surface mixed layer that contained siltier sediment with low excess 210Pb concentrations, possibly resulting from the resuspension, posterior advection, and eventual deposition of coarser and older sediment from adjacent trawling grounds. In contrast, the trawled site was eroded and presented compacted century-old sediment highly depleted in OM components, which were between 20 % and 60 % lower than those in the untrawled site. However, the upper 2 cm of the trawled site consisted of recently accumulated sediments enriched in excess 234Th, excess 210Pb, and phytopigments, while OM contents were similar to those from the untrawled core. This fresh sediment supported protein turnover rates of 0.025 d−1, which doubled those quantified in surface sediments of the untrawled site. The enhancement of remineralization rates in surface sediment of the trawled site was associated with the arrival of fresh particles on a chronically trawled deep-sea region that is generally deprived of OM. We conclude that the detrimental effects of bottom trawling can be temporarily and partially abated by the arrival of fresh and nutritionally rich OM, which stimulate the response of benthic communities. However, these ephemeral deposits are likely to be swiftly eroded due to the high trawling frequency over fishing grounds, highlighting the importance of establishing science-based management strategies to mitigate the impacts of bottom trawling.

We present a Lagrangian convective transport scheme developed for global chemistry and transport models, which considers the variable residence time that an air parcel spends in convection. This is particularly important for accurately simulating the tropospheric chemistry of short-lived species, e.g., for determining the time available for heterogeneous chemical processes on the surface of cloud droplets. In current Lagrangian convective transport schemes air parcels are stochastically redistributed within a fixed time step according to estimated probabilities for convective entrainment as well as the altitude of detrainment. We introduce a new scheme that extends this approach by modeling the variable time that an air parcel spends in convection by estimating vertical updraft velocities. Vertical updraft velocities are obtained by combining convective mass fluxes from meteorological analysis data with a parameterization of convective area fraction profiles. We implement two different parameterizations: a parameterization using an observed constant convective area fraction profile and a parameterization that uses randomly drawn profiles to allow for variability. Our scheme is driven by convective mass fluxes and detrainment rates that originate from an external convective parameterization, which can be obtained from meteorological analysis data or from general circulation models. We study the effect of allowing for a variable time that an air parcel spends in convection by performing simulations in which our scheme is implemented into the trajectory module of the ATLAS chemistry and transport model and is driven by the ECMWF ERA-Interim reanalysis data. In particular, we show that the redistribution of air parcels in our scheme conserves the vertical mass distribution and that the scheme is able to reproduce the convective mass fluxes and detrainment rates of ERA-Interim. We further show that the estimated vertical updraft velocities of our scheme are able to reproduce wind profiler measurements performed in Darwin, Australia, for velocities larger than 0.6 m s−1. SO2 is used as an example to show that there is a significant effect on species mixing ratios when modeling the time spent in convective updrafts compared to a redistribution of air parcels in a fixed time step. Furthermore, we perform long-time global trajectory simulations of radon-222 and compare with aircraft measurements of radon activity.

The concentration of CO2 in the atmosphere is sensitive to changes in the depth at which sinking particulate organic matter is remineralized: often described as a change in the exponent “b” of the Martin curve. Sediment trap observations from deep and intermediate depths suggest there is a spatially heterogeneous pattern of b, particularly varying with latitude, but disagree over the exact spatial patterns. Here we use a biogeochemical model of the phosphorus cycle coupled with a steady-state representation of ocean circulation to explore the sensitivity of preformed phosphate and atmospheric CO2 to spatial variability in remineralization depths. A Latin hypercube sampling method is used to simultaneously vary the Martin curve independently within 15 different regions, as a basis for a regression-based analysis used to derive a quantitative measure of sensitivity. Approximately 30 % of the sensitivity of atmospheric CO2 to changes in remineralization depths is driven by changes in the subantarctic region (36 to 60∘ S) similar in magnitude to the Pacific basin despite the much smaller area and lower export production. Overall, the absolute magnitude of sensitivity is controlled by export production, but the relative spatial patterns in sensitivity are predominantly constrained by ocean circulation pathways. The high sensitivity in the subantarctic regions is driven by a combination of high export production and the high connectivity of these regions to regions important for the export of preformed nutrients such as the Southern Ocean and North Atlantic. Overall, regionally varying remineralization depths contribute to variability in CO2 of between around 5 and 15 ppm, relative to a global mean change in remineralization depth. Future changes in the environmental and ecological drivers of remineralization, such as temperature and ocean acidification, are expected to be most significant in the high latitudes where CO2 sensitivity to remineralization is also highest. The importance of ocean circulation pathways to the high sensitivity in subantarctic regions also has significance for past climates given the importance of circulation changes in the Southern Ocean.

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.