• European Marine Science
• Other research productsclose
• Open Access close
• European Commission close
• US close

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.

Here we present a merged and calibrated dataset of temperature, practical salinity and dissolved organic matter (DOM) fluorescence obtained from several Ice Tethered Profilers (ITPs) deployed across the central Arctic (2011-2016). The data offer a unique spatial coverage of the distribution of DOM in the surface 800 m below Arctic sea ice. A total of 5044 profiles are gathered. The ITP data are level 3 data products pressure-bin-averaged at 1-db vertical resolution with depth down to either 200 or approximately 750 m. Data (max 800m depth) from CTD casts made during two oceanographic cruises are also included. These were used as part of the calibration and validation of the ITP calibration routines. The cruises were PS94 (ARK-XXIX/3) with POLARSTERN in 2015 and NAACOS with DANA in 2012. The presented DOM fluorescence data are smoothed, corrected for instrument drift and calibrated to provide intercomparable data across the sensors. Fluorescence is reported in Raman Units (nm-1), and comparable to laboratory measurements conducted according to current community recommendations.

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.

The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting heat flux anomalies with various locations and intensities. In our model experiment, a minimum heat flux value of 970 mW m−2 located close to the East Greenland Ice-core Project (EGRIP) is required locally to reproduce the observed NEGIS velocities, giving basal melt rates consistent with previous estimates. The value cannot be attributed to geothermal heat flux alone and we suggest hydrothermal circulation as a potential explanation for the high local heat flux. By including high heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.

Ocean margin sediments have been considered as important sources of dissolved organic carbon (DOC) to the deep ocean, yet the contribution from advective settings has just started to be acknowledged. Here we present evidence showing that near-surface heating of sediment in the Guaymas Basin, a young extensional depression, causes mass production and discharge of reactive dissolved organic matter (DOM). In the sediment heated up to ~100 °C, we found unexpectedly low DOC concentrations in the pore waters, reflecting the combined effect of thermal desorption and advective fluid flow. Heating experiments suggested DOC production to be a rapid, abiotic process with the DOC concentration increasing exponentially with temperature. The high proportions of total hydrolyzable amino acids and presence of chemical species affiliated with activated hydrocarbons, carbohydrates and peptides indicate high reactivity of the DOM. Model simulation suggests that at the local scale, near-surface heating of sediment creates short and massive DOC discharge events that elevate the bottom-water DOC concentration. Because of the heterogeneous distribution of high heat flow areas, the expulsion of reactive DOM is spotty at any given time. We conclude that hydrothermal heating of young rift sediments alter deep-ocean budgets of bioavailable DOM, creating organic-rich habitats for benthic life.

Although it has been demonstrated that the speed and magnitude of the recent Arctic sea ice decline is unprecedented for the past 1450 years, few records are available to provide a paleoclimate context for Arctic sea ice extent. Bromine enrichment in ice cores has been suggested to indicate the extent of newly formed sea ice areas. Despite the similarities among sea ice indicators and ice core bromine enrichment records, uncertainties still exist regarding the quantitative linkages between bromine reactive chemistry and the first-year sea ice surfaces. Here we present a 120 000-year record of bromine enrichment from the RECAP (REnland ice CAP) ice core, coastal east Greenland, and interpret it as a record of first-year sea ice. We compare it to existing sea ice records from marine cores and tentatively reconstruct past sea ice conditions in the North Atlantic as far north as the Fram Strait (50–85∘ N). Our interpretation implies that during the last deglaciation, the transition from multi-year to first-year sea ice started at ∼17.5 ka, synchronously with sea ice reductions observed in the eastern Nordic Seas and with the increase in North Atlantic ocean temperature. First-year sea ice reached its maximum at 12.4–11.8 ka during the Younger Dryas, after which open-water conditions started to dominate, consistent with sea ice records from the eastern Nordic Seas and the North Icelandic shelf. Our results show that over the last 120 000 years, multi-year sea ice extent was greatest during Marine Isotope Stage (MIS) 2 and possibly during MIS 4, with more extended first-year sea ice during MIS 3 and MIS 5. Sea ice extent during the Holocene (MIS 1) has been less than at any time in the last 120 000 years.

In the framework of the EU-FP7 BACCHUS (impact of Biogenic versus Anthropogenic emissions on Clouds and Climate: towards a Holistic UnderStanding) project, an intensive field campaign was performed in Cyprus (March 2015). Remotely piloted aircraft system (RPAS), ground-based instruments, and remote-sensing observations were operating in parallel to provide an integrated characterization of aerosol–cloud interactions. Remotely piloted aircraft (RPA) were equipped with a five-hole probe, pyranometers, pressure, temperature and humidity sensors, and measured vertical wind at cloud base and cloud optical properties of a stratocumulus layer. Ground-based measurements of dry aerosol size distributions and cloud condensation nuclei spectra, and RPA observations of updraft and meteorological state parameters are used here to initialize an aerosol–cloud parcel model (ACPM) and compare the in situ observations of cloud optical properties measured by the RPA to those simulated in the ACPM. Two different cases are studied with the ACPM, including an adiabatic case and an entrainment case, in which the in-cloud temperature profile from RPA is taken into account. Adiabatic ACPM simulation yields cloud droplet number concentrations at cloud base (approximately 400 cm−3) that are similar to those derived from a Hoppel minimum analysis. Cloud optical properties have been inferred using the transmitted fraction of shortwave radiation profile measured by downwelling and upwelling pyranometers mounted on a RPA, and the observed transmitted fraction of solar radiation is then compared to simulations from the ACPM. ACPM simulations and RPA observations shows better agreement when associated with entrainment compared to that of an adiabatic case. The mean difference between observed and adiabatic profiles of transmitted fraction of solar radiation is 0.12, while this difference is only 0.03 between observed and entrainment profiles. A sensitivity calculation is then conducted to quantify the relative impacts of 2-fold changes in aerosol concentration, and updraft to highlight the importance of accounting for the impact of entrainment in deriving cloud optical properties, as well as the ability of RPAs to leverage ground-based observations for studying aerosol–cloud interactions.

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.

The mid-Pliocene is a valuable time interval for investigating equilibrium climate at current atmospheric CO2 concentrations because atmospheric CO2 concentrations are thought to have been comparable to the current day and yet the climate and distribution of ecosystems were quite different. One intriguing, but not fully understood, feature of the early to mid-Pliocene climate is the amplified Arctic temperature response and its impact on Arctic ecosystems. Only the most recent models appear to correctly estimate the degree of warming in the Pliocene Arctic and validation of the currently proposed feedbacks is limited by scarce terrestrial records of climate and environment. Here we reconstruct the summer temperature and fire regime from a subfossil fen-peat deposit on west–central Ellesmere Island, Canada, that has been chronologically constrained using cosmogenic nuclide burial dating to 3.9+1.5/-0.5 Ma. The estimate for average mean summer temperature is 15.4±0.8 ∘C using specific bacterial membrane lipids, i.e., branched glycerol dialkyl glycerol tetraethers. This is above the proposed threshold that predicts a substantial increase in wildfire in the modern high latitudes. Macro-charcoal was present in all samples from this Pliocene section with notably higher charcoal concentration in the upper part of the sequence. This change in charcoal was synchronous with a change in vegetation that included an increase in abundance of fire-promoting Pinus and Picea. Paleo-vegetation reconstructions are consistent with warm summer temperatures, relatively low summer precipitation and an incidence of fire comparable to fire-adapted boreal forests of North America and central Siberia. To our knowledge, this site provides the northernmost evidence of fire during the Pliocene. It suggests that ecosystem productivity was greater than in the present day, providing fuel for wildfires, and that the climate was conducive to the ignition of fire during this period. The results reveal that interactions between paleo-vegetation and paleoclimate were mediated by fire in the High Arctic during the Pliocene, even though CO2 concentrations were similar to modern values.