Spatially explicit methane inventory for Switzerland with 500 m x 500 m resolution for 2011 (NetCDF 25 MB)

We present the first high-resolution (500 m × 500 m) gridded methane (CH4) emission inventory for Switzerland, which integrates the national emission totals reported to the United Nations Framework Convention on Climate Change (UNFCCC) and recent CH4 flux studies conducted by research groups across Switzerland. In addition to anthropogenic emissions, we also include natural and semi-natural CH4 fluxes, i.e., emissions from lakes and reservoirs, wetlands, wild animals as well as uptake by forest soils. National CH4 emissions were disaggregated using detailed geostatistical information on source locations and their spatial extent and process- or area-specific emission factors. In Switzerland, the highest CH4 emissions in 2011 originated from the agricultural sector (150 Gg CH4/yr), mainly produced by ruminants and manure management, followed by emissions from waste management (15 Gg CH4/yr) mainly from landfills and the energy sector (12 Gg CH4/yr), which was dominated by emissions from natural gas distribution. Compared to the anthropogenic sources, emissions from natural and semi-natural sources were relatively small (6 Gg CH4/yr), making up only 3 % of the total emissions in Switzerland. CH4 fluxes from agricultural soils were estimated to be not significantly different from zero (between -1.5 and 0 Gg CH4/yr), while forest soils are a CH4 sink (approx. -2.8 Gg CH4/yr), partially offsetting other natural emissions. Estimates of uncertainties are provided for the different sources, including an estimate of spatial disaggregation errors deduced from a comparison with a global (EDGAR v4.2) and a European CH4 inventory (TNO/MACC). This new spatially-explicit emission inventory for Switzerland will provide valuable input for regional scale atmospheric modeling and inverse source estimation.

Antarctic dataset in NetCDF format

The dataset described in this document has been put together for the purposes of numerical ice sheet modelling of the Antarctic Ice Sheet (AIS), containing data on the ice sheet configuration (e.g. ice surface and ice thickness) and boundary conditions, such as the surface air temperature and accumulation. It is now possible to download a community ice sheet model (e.g. Glimmer-CISM, Rutt et al., 2009 doi:10.1029/2008JF001015), but without adequate data it is difficult to utilise such models. More specifically, ice sheet models that are initialised and run forward from the present day ice sheet configuration, need input data to represent the present-day ice sheet configuration as closely as possible (unlike those spun-up from ice free conditions, which only require the bed/bathymetry). Whilst the BEDMAP dataset (Lythe et al., 2001) was a step forward when it was made, there are a number of inconsistencies within the dataset (see Section 3), and since its release, more data has become available. The dataset described here incorporates some major new datasets (e.g. AGASEA/BBAS ice thickness, Nitsche et al. (2006) bathymetry doi:10.1029/2007GC001694), but by no means incorporates all the new data available. This considerable task is left for a 'BEDMAP2', (an updated version of BEDMAP), however, the processing carried out in this document illustrates the requirements of a dataset for the purpose of high resolution ice sheet modelling, and bridges the gap until a BEDMAP2 is published. It is envisaged, however, that updated versions of the data set will be made available periodically when new regional data sets become available and can be readily incorporated.

The CMCC Eddy-permitting Global Ocean Physical Reanalysis (C-GLORS v5, 1980-2014), links to NetCDF files

The CMCC Global Ocean Physical Reanalysis System (C-GLORS) is used to simulate the state of the ocean in the last decades. It consists of a variational data assimilation system (OceanVar), capable of assimilating all in-situ observations along with altimetry data, and a forecast step performed by the ocean model NEMO coupled with the LIM2 sea-ice model. KEY STRENGTHS: - Data are available for a large number of ocean parameters - An extensive validation has been conducted and is freely available - The reanalysis is performed at high resolution (1/4 degree) and spans the last 30 years KEY LIMITATIONS: - Quality may be discontinuos and depend on observation coverage - Uncertainty estimates are simply derived through verification skill scores

Global distributions of microzooplankton abundance and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

Microzooplankton database. Originally published in: Buitenhuis, Erik, Richard Rivkin, Sévrine Sailley, Corinne Le Quéré (2010) Biogeochemical fluxes through microzooplankton. Global Biogeochemical Cycles Vol. 24, GB4015, doi:10.1029/2009GB003601 This new version has had some mistakes corrected.

Global Ocean Biomes: Mean and time-varying maps (NetCDF 7.8 MB)

Large-scale studies of ocean biogeochemistry and carbon cycling have often partitioned the ocean into regions along lines of latitude and longitude despite the fact that spatially more complex boundaries would be closer to the true biogeography of the ocean. Herein, we define 17 open-ocean biomes classified from four observational data sets: sea surface temperature (SST), spring/summer chlorophyll a concentrations (Chl a), ice fraction, and maximum mixed layer depth (maxMLD) on a 1° × 1° grid. By considering interannual variability for each input, we create dynamic ocean biome boundaries that shift annually between 1998 and 2010. Additionally we create a core biome map, which includes only the grid cells that do not change biome assignment across the 13 years of the time-varying biomes. These biomes can be used in future studies to distinguish large-scale ocean regions based on biogeochemical function.

World distribution of land cover changes, model in NetCDF format

The role of Pre- and Protohistoric anthropogenic land cover changes needs to be quantified i) to establish a baseline for comparison with current human impact on the environment and ii) to separate it from naturally occurring changes in our environment. Results are presented from the simple, adaptation-driven, spatially explicit Global Land Use and technological Evolution Simulator (GLUES) for pre-Bronze age demographic, technological and economic change. Using scaling parameters from the History Database of the Global Environment as well as GLUES-simulated population density and subsistence style, the land requirement for growing crops is estimated. The intrusion of cropland into potentially forested areas is translated into carbon loss due to deforestation with the dynamic global vegetation model VECODE. The land demand in important Prehistoric growth areas - converted from mostly forested areas - led to large-scale regional (country size) deforestation of up to 11% of the potential forest. In total, 29 Gt carbon were lost from global forests between 10 000 BC and 2000 BC and were replaced by crops; this value is consistent with other estimates of Prehistoric deforestation. The generation of realistic (agri-)cultural development trajectories at a regional resolution is a major strength of GLUES. Most of the pre-Bronze age deforestation is simulated in a broad farming belt from Central Europe via India to China. Regional carbon loss is, e.g., 5 Gt in Europe and the Mediterranean, 6 Gt on the Indian subcontinent, 18 Gt in East and Southeast Asia, or 2.3 Gt in subsaharan Africa.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

Larcform 1 - single column model intercomparison of Arctic air formation, link to model results in NetCDF format

Weather and climate models struggle to represent lower tropospheric temperature and moisture profiles and surface fluxes in Arctic winter, not least because they lack or misrepresent physical processes that are specific to high latitudes. The Arctic boundary layer in winter has been observed to be in either a radiatively clear or cloudy state: The radiatively clear state is characterized by strong surface radiative cooling leading to the build-up of surface-based temperature inversions, whereas the cloudy state occurs when cloud liquid water is present in the atmospheric column, allowing little or no surface radiative cooling and leading to weaker and typically elevated temperature inversions. Many large-scale models have been shown to lack the cloudy state, and some do substantially underestimate stability in the clear state. We here present results from the first Lagrangian ARCtic air FORMation experiment (Larcform 1), a GASS (Global atmospheric system studies) single-column model intercomparison which reproduces these biases of large-scale models in an idealised setup.

Objective Mapping of Argo data in the Weddell Gyre: a gridded dataset of upper ocean water properties, link to data files in NetCDF format

The Weddell Gyre plays a crucial role in the regulation of climate by transferring heat into the deep ocean through deep and bottom water mass formation. However, our understanding of Weddell Gyre water mass properties is limited to regions of data availability, primarily along the Prime Meridian. The aim is to provide a dataset of the upper water column properties of the entire Weddell Gyre. Objective mapping was applied to Argo float data in order to produce spatially gridded, time composite maps of temperature and salinity for fixed pressure levels ranging from 50 to 2000 dbar, as well as temperature, salinity and pressure at the level of the sub-surface temperature maximum. While the data are currently too limited to incorporate time into the gridded structure, the data are extensive enough to produce maps of the entire region across three time composite periods (2002-2005, 2006-2009 and 2010-2013), which can be used to determine how representative conclusions drawn from data collected along general RV transect lines are on a gyre scale perspective. The time composite data sets are provided as netCDF files; one for each time period. Mapped fields of conservative temperature, absolute salinity and potential density are provided for 41 vertical pressure levels. The above variables as well as pressure are provided at the level of the sub-surface temperature maximum. Corresponding mapping errors are also included in the netCDF files. Further details are provided in the global attributes, such as the unit variables and structure of the corresponding data array (i.e. latitude x longitude x vertical pressure level). In addition, all files ending in "_potTpSal" provide mapped fields of potential temperature and practical salinity.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Optical depth of Atlantic dust as netCDF files

African dust outbreaks are the result of complex interactions between the land, atmosphere, and oceans, and only recently has a large body of work begun to emerge that aims to understand the controls on-and impacts of-African dust. At the same time, long-term records of dust outbreaks are either inferred from visibility data from weather stations or confined to a few in situ observational sites. Satellites provide the best opportunity for studying the large-scale characteristics of dust storms, but reliable records of dust are generally on the scale of a decade or less. Here the authors develop a simple model for using modern and historical data from meteorological satellites, in conjunction with a proxy record for atmospheric dust, to extend satellite-retrieved dust optical depth over the northern tropical Atlantic Ocean from 1955 to 2008. The resultant 54-yr record of dust has a spatial resolution of 1° and a monthly temporal resolution. From analysis of the historical dust data, monthly tropical northern Atlantic dust cover is bimodal, has a strong annual cycle, peaked in the early 1980s, and shows minimums in dustiness during the beginning and end of the record. These dust optical depth estimates are used to calculate radiative forcing and heating rates from the surface through the top of the atmosphere over the last half century. Radiative transfer simulations show a large net negative dust forcing from the surface through the top of the atmosphere, also with a distinct annual cycle, and mean tropical Atlantic monthly values of the surface forcing range from -3 to -9 W/m**2. Since the surface forcing is roughly a factor of 3 larger in magnitude than the top-of-the-atmosphere forcing, there is also a positive heating rate of the midtroposphere by dust.