• European Marine Science
• Other research productsclose
• Open Access close
• CH close

We present an age model for the 651 m deep Skytrain Ice Rise ice core (79°44.5'S, 78°32.7'W). The top 2000 years have previously been dated using age markers interpolated through annual layer counting. Below this, we align the Skytrain core to the AICC2012 age model using tie points in the ice and air phase, and apply the Paleochrono program to obtain the best fit to the tie points and glaciological constraints. In the gas phase, ties are made using methane and, in critical sections, δ18Oair; in the ice phase ties are through 10Be across the Laschamps Event, and through ice chemistry related to long-range dust transport and deposition. This strategy provides a good outcome to about 108 ka (~605 m). Beyond that there are signs of flow disturbance, with a section of ice probably repeated. Nonetheless values of CH4 and δ18Oair confirm that part of the last interglacial (LIG), from about 117-126 ka (617-628 m), is present and in chronological order. Below this there are clear signs of stratigraphic disturbance, with rapid oscillation of values in both the ice and gas phase at the base of the LIG section. Based on methane values, the warmest part of the LIG and the coldest part of the penultimate glacial are missing from our record. Ice below 631 m appears to be of age >150 ka.

Micropaleontological (coccolith, planktonic foraminifera) and geochemical (CaCO3, CaXRF, δ18O N. pachyderma, δ13C N. pachyderma) analyses from sediment core MD04-2718 retrieved from the Subantarctic zone (SAZ) of the Indian Southern Ocean covering the time interval from Marine Isotope Stage (MIS) 12 to MIS 10 (440,000– 360,000 years). Coccolith abundances and masses were performed using SYRACO software at Geoscience Paris-Saclay Laboratory (GEOPS). CaCO3 percentages were also measured at GEOPS. Foraminifera abundances and masses and stable oxygen and carbon isotopes were performed at Laboratoire des Sciences du Climat et de l'Environnement (LSCE). SST were calculated using two methods : the Modern Analogue Technique (Haddam et al., 2016) and the percentage of N. pachyderma s. (Govin et al., 2009). CaXRF measurements were performed at ETH, Zurich. These approaches allowed to decipher the variations in Carbonate Counter Pump and upwelling strength in the Indian Sector of the Southern Ocean between MIS 12 and MIS 10, including the Termination V and interglacial MIS 11.

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.

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.

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.

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.

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.

All data are from core ODP 1094 recovered from the Antarctic Zone of the Southern Ocean (Atlantic sector). Age model of ODP 1094 (1.5Ma) dδ18O, Mg/Ca, Mn/Ca and Mg/Ca-derived sea surface temperature and surface water d18O based on down core measurements of planktic Neogloboquadrina pachyderma (sinistral) from core ODP 1094 (downcore data and averaged for MIS). δ18O, Mg/Ca, Mn/Ca and Mg/Ca-derived bottom water temperature and bottom water d18O based on down core measurements of benthic Melonis pompilioides from core ODP 1094 (downcore data and averaged for MIS). δ18O of benthic Cibicidoides spp. from core ODP 1094.