Sulfate and nitrate firn concentrations on the Greenland ice sheet

Intercomparison of three new chemical ice core records from northern Greenland (covering the time span from approximately 1500 A.D. to present) with previously published records for southern and central Greenland reveals a uniform timing of anthropogenic changes in sulfate and nitrate firn concentrations over the entire ice sheet. The anthropogenic sulfate increase started around 1890, was interrupted by a transient decrease in the 1930s, and has resumed a major increase since 1950. Since the late 1970s though, a significant 30% decline in Greenland sulfate firn levels can be documented. The maximum anthropogenic increase in northern Greenland sulfate firn concentrations (up to 200-230 ppb) is 2-3 times larger than in southern and central Greenland. Nitrate records show an essentially steady increase since 1950 and, documented for the first time, a slight reduction during most recent years. Maximum nitrate firn levels of 100-130 ppb exceed the preindustrial background by 100% all over the Greenland ice sheet. Comparison with anthropogenic SO2 and NO x emission records indicates that the major increase in sulfate firn concentrations since 1950 can be attributed to Eurasian sources, while firn levels during the first half of this century appear to be dominated by North American emissions. A stronger North American source contribution is indicated over the entire 20th century in the case of nitrate. Applying a macroscopic deposition model separate time series for wet and dry deposition were derived which revealed a close correspondence of wet deposited sulfate with the timing of U.S. emissions, while the temporal evolution of Eurasian emissions is mainly reflected in the dry sulfate deposition record. During this century wet sulfate deposition increased by a factor of two while the total dry sulfate deposition flux increased by more than 500%. Wet and dry nitrate deposition both increased by 100% during the same period.

Sources and mixing state of size-resolved elemental carbon particles in a European megacity: Paris

An Aerosol Time-Of-Flight Mass Spectrometer (ATOFMS) was deployed to investigate the size-resolved chemical composition of single particles at an urban background site in Paris, France, as part of the MEGAPOLI winter campaign in January/February 2010. ATOFMS particle counts were scaled to match coincident Twin Differential Mobility Particle Sizer (TDMPS) data in order to generate hourly size-resolved mass concentrations for the single particle classes observed. The total scaled ATOFMS particle mass concentration in the size range 150–1067 nm was found to agree very well with the sum of concurrent High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS) and Multi-Angle Absorption Photometer (MAAP) mass concentration measurements of organic carbon (OC), inorganic ions and black carbon (BC) (R2 = 0.91). Clustering analysis of the ATOFMS single particle mass spectra allowed the separation of elemental carbon (EC) particles into four classes: (i) EC attributed to biomass burning (ECbiomass), (ii) EC attributed to traffic (ECtraffic), (iii) EC internally mixed with OC and ammonium sulfate (ECOCSOx), and (iv) EC internally mixed with OC and ammonium nitrate (ECOCNOx). Average hourly mass concentrations for EC-containing particles detected by the ATOFMS were found to agree reasonably well with semi-continuous quantitative thermal/optical EC and optical BC measurements (r2 = 0.61 and 0.65–0.68 respectively, n = 552). The EC particle mass assigned to fossil fuel and biomass burning sources also agreed reasonably well with BC mass fractions assigned to the same sources using seven-wavelength aethalometer data (r2 = 0.60 and 0.48, respectively, n = 568). Agreement between the ATOFMS and other instrumentation improved noticeably when a period influenced by significantly aged, internally mixed EC particles was removed from the intercomparison. 88% and 12% of EC particle mass was apportioned to fossil fuel and biomass burning respectively using the ATOFMS data compared with 85% and 15% respectively for BC estimated from the aethalometer model. On average, the mass size distribution for EC particles is bimodal; the smaller mode is attributed to locally emitted, mostly externally mixed EC particles, while the larger mode is dominated by aged, internally mixed ECOCNOx particles associated with continental transport events. Periods of continental influence were identified using the Lagrangian Particle Dispersion Model (LPDM) "FLEXPART". A consistent minimum between the two EC mass size modes was observed at approximately 400 nm for the measurement period. EC particles below this size are attributed to local emissions using chemical mixing state information and contribute 79% of the scaled ATOFMS EC particle mass, while particles above this size are attributed to continental transport events and contribute 21% of the EC particle mass. These results clearly demonstrate the potential benefit of monitoring size-resolved mass concentrations for the separation of local and continental EC emissions. Knowledge of the relative input of these emissions is essential for assessing the effectiveness of local abatement strategies.

Global distributions of diazotrophs Gamma-A nifH genes abundance - Depth integrated values computed from a collection of source datasets - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. The present data set presents depth integrated values of diazotrophs Gamma-A nifH genes abundance, computed from a collection of source data sets.

Global distributions of diazotrophs abundance, biomass and nitrogen fixation rates - Collection of source datasets - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. The present collection presents the original data sets used to compile Global distributions of diazotrophs abundance, biomass and nitrogen fixation rates

Global distributions of diazotrophs abundance and biomass - Depth integrated values computed from a collection of source datasets - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. The present data set presents depth integrated values of diazotrophs abundance and biomass, computed from a collection of source data sets.

Global distributions of diazotrophs nitrogen fixation rates - Depth integrated values computed from a collection of source datasets - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. The present data set presents depth integrated values of diazotrophs nitrogen fixation rates, computed from a collection of source data sets.

Global distributions of diazotrophs abundance, biomass and nitrogen fixation rates - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. This is a gridded data product about diazotrophic organisms . There are 6 variables. Each variable is gridded on a dimension of 360 (longitude) * 180 (latitude) * 33 (depth) * 12 (month). The first group of 3 variables are: (1) number of biomass observations, (2) biomass, and (3) special nifH-gene-based biomass. The second group of 3 variables is same as the first group except that it only grids non-zero data. We have constructed a database on diazotrophic organisms in the global pelagic upper ocean by compiling more than 11,000 direct field measurements including 3 sub-databases: (1) nitrogen fixation rates, (2) cyanobacterial diazotroph abundances from cell counts and (3) cyanobacterial diazotroph abundances from qPCR assays targeting nifH genes. Biomass conversion factors are estimated based on cell sizes to convert abundance data to diazotrophic biomass. Data are assigned to 3 groups including Trichodesmium, unicellular diazotrophic cyanobacteria (group A, B and C when applicable) and heterocystous cyanobacteria (Richelia and Calothrix). Total nitrogen fixation rates and diazotrophic biomass are calculated by summing the values from all the groups. Some of nitrogen fixation rates are whole seawater measurements and are used as total nitrogen fixation rates. Both volumetric and depth-integrated values were reported. Depth-integrated values are also calculated for those vertical profiles with values at 3 or more depths.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Global impacts of the 1980s regime shift

Despite evidence from a number of Earth systems that abrupt temporal changes known as regime shifts are important, their nature, scale and mechanisms remain poorly documented and understood. Applying principal component analysis, change-point analysis and a sequential t-test analysis of regime shifts to 72 time series, we confirm that the 1980s regime shift represented a major change in the Earth's biophysical systems from the upper atmosphere to the depths of the ocean and from the Arctic to the Antarctic, and occurred at slightly different times around the world. Using historical climate model simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and statistical modelling of historical temperatures, we then demonstrate that this event was triggered by rapid global warming from anthropogenic plus natural forcing, the latter associated with the recovery from the El Chichón volcanic eruption. The shift in temperature that occurred at this time is hypothesized as the main forcing for a cascade of abrupt environmental changes. Within the context of the last century or more, the 1980s event was unique in terms of its global scope and scale; our observed consequences imply that if unavoidable natural events such as major volcanic eruptions interact with anthropogenic warming unforeseen multiplier effects may occur.

Global impacts of the 1980s regime shift

Despite evidence from a number of Earth systems that abrupt temporal changes known as regime shifts are important, their nature, scale and mechanisms remain poorly documented and understood. Applying principal component analysis, change-point analysis and a sequential t-test analysis of regime shifts to 72 time series, we confirm that the 1980s regime shift represented a major change in the Earth's biophysical systems from the upper atmosphere to the depths of the ocean and from the Arctic to the Antarctic, and occurred at slightly different times around the world. Using historical climate model simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and statistical modelling of historical temperatures, we then demonstrate that this event was triggered by rapid global warming from anthropogenic plus natural forcing, the latter associated with the recovery from the El Chichón volcanic eruption. The shift in temperature that occurred at this time is hypothesized as the main forcing for a cascade of abrupt environmental changes. Within the context of the last century or more, the 1980s event was unique in terms of its global scope and scale; our observed consequences imply that if unavoidable natural events such as major volcanic eruptions interact with anthropogenic warming unforeseen multiplier effects may occur.

Solar and thermal radiative effects during the 2011 extreme desert dust episode over Portugal

This paper analyses the influence of the extreme Saharan desert dust (DD) event on shortwave (SW) and longwave (LW) radiation at the EARLINET/AERONET Évora station (Southern Portugal) from 4 up to 7 April 2011. There was also some cloud occurrence in the period. In this context, it is essential to quantify the effect of cloud presence on aerosol radiative forcing. A radiative transfer model was initialized with aerosol optical properties, cloud vertical properties and meteorological atmospheric vertical profiles. The intercomparison between the instantaneous TOA shortwave and longwave fluxes derived using CERES and those calculated using SBDART, which was fed with aerosol extinction coefficients derived from the CALIPSO and lidar-PAOLI observations, varying OPAC dataset parameters, was reasonably acceptable within the standard deviations. The dust aerosol type that yields the best fit was found to be the mineral accumulation mode. Therefore, SBDART model constrained with the CERES observations can be used to reliably determine aerosol radiative forcing and heating rates. Aerosol radiative forcings and heating rates were derived in the SW (ARFSw, AHRSw) and LW (ARFLw, AHRLw) spectral ranges, considering a cloud-aerosol free reference atmosphere. We found that AOD at 440 nm increased by a factor of 5 on 6 April with respect to the lower dust load on 4 April. It was responsible by a strong cooling radiative effect pointed out by the ARFSw value (−99 W/m2 for a solar zenith angle of 60°) offset by a warming radiative effect according to ARFLw value (+21.9 W/m2) at the surface. Overall, about 24% and 12% of the dust solar radiative cooling effect is compensated by its longwave warming effect at the surface and at the top of the atmosphere, respectively. Hence, larger aerosol loads could enhance the response between the absorption and re-emission processes increasing the ARFLw with respect to those associated with moderate and low aerosol loads. The unprecedented results derived from this work complement the findings in other regions on the modifications of radiative energy budget by the dust aerosols, which could have relevant influences on the regional climate and will be topics for future investigations.

The role of the Rocky Mountains in shaping the atmospheric mean state and the response to tropical forcing in idealised simulations.

Understanding the natural and forced variability of the atmospheric general circulation and its drivers is one of the grand challenges in climate science. It is of paramount importance to understand to what extent the systematic error of climate models affects the processes driving such variability. This is done by performing a set of simulations (ROCK experiments) with an intermediate complexity atmospheric model (SPEEDY), in which the Rocky Mountains orography is increased or decreased to influence the structure of the North Pacific jet stream. For each of these modified-orography experiments, the climatic response to idealized sea surface temperature anomalies of varying intensity in the El Niño Southern Oscillation (ENSO) region is studied. ROCK experiments are characterized by variations in the Pacific jet stream intensity whose extension encompasses the spread of the systematic error found in Coupled Model Intercomparison Project (CMIP6) models. When forced with ENSO-like idealised anomalies, they exhibit a non-negligible sensitivity in the response pattern over the Pacific North American region, indicating that the model mean state can affect the model response to ENSO. It is found that the classical Rossby wave train response to ENSO is more meridionally oriented when the Pacific jet stream is weaker and more zonally oriented with a stronger jet. Rossby wave linear theory suggests that a stronger jet implies a stronger waveguide, which traps Rossby waves at a lower latitude, favouring a zonal propagation of Rossby waves. The shape of the dynamical response to ENSO affects the ENSO impacts on surface temperature and precipitation over Central and North America. A comparison of the SPEEDY results with CMIP6 models suggests a wider applicability of the results to more resources-demanding climate general circulation models (GCMs), opening up to future works focusing on the relationship between Pacific jet misrepresentation and response to external forcing in fully-fledged GCMs.