High-resolution in situ oxygen microprofiles, porewater and solid phase geochemistry from the Crimean shelf (Black Sea) from Maria S. Merian cruise MSM15/1

The outer western Crimean shelf of the Black Sea is a natural laboratory to investigate effects of stable oxic versus varying hypoxic conditions on seafloor biogeochemical processes and benthic community structure. Bottom-water oxygen concentrations ranged from normoxic (175 µmol O2/L) and hypoxic (< 63 µmol O2/L) or even anoxic/sulfidic conditions within a few kilometers' distance. Variations in oxygen concentrations between 160 and 10 µmol/L even occurred within hours close to the chemocline at 134 m water depth. Total oxygen uptake, including diffusive as well as fauna-mediated oxygen consumption, decreased from 15 mmol/m**2/d on average in the oxic zone, to 7 mmol/m**2/d on average in the hypoxic zone, correlating with changes in macrobenthos composition. Benthic diffusive oxygen uptake rates, comprising respiration of microorganisms and small meiofauna, were similar in oxic and hypoxic zones (on average 4.5 mmol/m**2/d), but declined to 1.3 mmol/m**2/d in bottom waters with oxygen concentrations below 20 µmol/L. Measurements and modeling of porewater profiles indicated that reoxidation of reduced compounds played only a minor role in diffusive oxygen uptake under the different oxygen conditions, leaving the major fraction to aerobic degradation of organic carbon. Remineralization efficiency decreased from nearly 100 % in the oxic zone, to 50 % in the oxic-hypoxic zone, to 10 % in the hypoxic-anoxic zone. Overall, the faunal remineralization rate was more important, but also more influenced by fluctuating oxygen concentrations, than microbial and geochemical oxidation processes.

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.

Numerical modeling of marine radiocarbon reservoir ages during the last deglaciation

A critical problem in radiocarbon dating is the spatial and temporal variability of marine reservoir ages (MRAs). We assessed the MRA evolution during the last deglaciation by numerical modeling, applying a self-consistent iteration scheme in which an existing radiocarbon chronology (derived by Hughen et al., Quat. Sci. Rev., 25, pp. 3216-3227, 2006) was readjusted by transient, 3-D simulations of marine and atmospheric Delta14C. To estimate the uncertainties regarding the ocean ventilation during the last deglaciation, we considered various ocean overturning scenarios which are based on different climatic background states (PD: modern climate, GS: LGM climate conditions). Minimum and maximum MRAs are included in file 'MRAminmax_21-14kaBP.nc'. Three further files include MRAs according to equilibrium simulations of the preindustrial ocean (file 'C14age_preindustrial.nc'; this is an update of our results published in 2005) and of the glacial ocean (files 'C14age_spinupLGM_GS.nc' and 'C14age_spinupLGM_PD.nc').

Multicentennial resolution proxy records from splicings of MD96-2080 and MD02-2594, Agulhas Bank Splice (ABS)

[1] Planktonic d18O and Mg/Ca-derived sea surface temperature (SST) records from the Agulhas Corridor off South Africa display a progressive increase of SST during glacial periods of the last three climatic cycles. The SST increases of up to 4°C coincide with increased abundance of subtropical planktonic foraminiferal marker species which indicates a progressive warming due to an increased influence of subtropical waters at the core sites. Mg/Ca-derived SST maximizes during glacial maxima and glacial Terminations to values about 2.5°C above full-interglacial SST. The paired planktonic d18O and Mg/Ca-derived SST records yield glacial seawater d18O anomalies of up to 0.8 per mill, indicating measurably higher surface salinities during these periods. The SST pattern along our record is markedly different from a UK'37-derived SST record at a nearby core location in the Agulhas Corridor that displays SST maxima only during glacial Terminations. Possible explanations are lateral alkenone advection by the vigorous regional ocean currents or the development of SST contrasts during glacials in association with seasonal changes of Agulhas water transports and lateral shifts of the Agulhas retroflection. The different SST reconstructions derived from UK'37 and Mg/Ca pose a significant challenge to the interpretation of the proxy records and demonstrate that the reconstruction of the Agulhas Current and interocean salt leakage is not as straightforward as previously suggested.

High-resolution digital elevation model (DEM) of a thermokarst depression in Siberian ice-rich permafrost deposits

This data set provides a high-resolution digital elevation model (DEM) of a thermokarst depression (~7 km²) on ice-complex deposits in the Arctic Lena Delta, Siberia. The DEM based on a geodetic field survey and was used for quantitative land surface analyses and detailed description of the thermokarst depression morphology. Detailed morphometrical analyses, volume calculations, and solar radiation modeling were performed and statistically analyzed by Ulrich et al. (2010) to investigate the asymmetrical thermokarst depression development and directed lake migration previously proposed by Morgenstern et al. (2008). Furthermore, the high-resolution DEM in combination with satellite data allowed detailed analyses of spatial and temporal landscape changes due to thermokarst development (Günther, 2009).

Data and Code archive for the interpolation of surface ocean carbon dioxide

We have developed a statistical gap-filling method adapted to the specific coverage and properties of observed fugacity of surface ocean CO2 (fCO2). We have used this method to interpolate the Surface Ocean CO2 Atlas (SOCAT) v2 database on a 2.5°×2.5° global grid (south of 70°N) for 1985-2011 at monthly resolution. The method combines a spatial interpolation based on a 'radius of influence' to determine nearby similar fCO2 values with temporal harmonic and cubic spline curve-fitting, and also fits long term trends and seasonal cycles. Interannual variability is established using deviations of observations from the fitted trends and seasonal cycles. An uncertainty is computed for all interpolated values based on the spatial and temporal range of the interpolation. Tests of the method using model data show that it performs as well as or better than previous regional interpolation methods, but in addition it provides a near-global and interannual coverage.

High-resolution dynamical downscaling products for the North Slope of Alaska and surrounding areas, links to model result files

Climatic changes are most pronounced in northern high latitude regions. Yet, there is a paucity of observational data, both spatially and temporally, such that regional-scale dynamics are not fully captured, limiting our ability to make reliable projections. In this study, a group of dynamical downscaling products were created for the period 1950 to 2100 to better understand climate change and its impacts on hydrology, permafrost, and ecosystems at a resolution suitable for northern Alaska. An ERA-interim reanalysis dataset and the Community Earth System Model (CESM) served as the forcing mechanisms in this dynamical downscaling framework, and the Weather Research & Forecast (WRF) model, embedded with an optimization for the Arctic (Polar WRF), served as the Regional Climate Model (RCM). This downscaled output consists of multiple climatic variables (precipitation, temperature, wind speed, dew point temperature, and surface air pressure) for a 10 km grid spacing at three-hour intervals. The modeling products were evaluated and calibrated using a bias-correction approach. The ERA-interim forced WRF (ERA-WRF) produced reasonable climatic variables as a result, yielding a more closely correlated temperature field than precipitation field when long-term monthly climatology was compared with its forcing and observational data. A linear scaling method then further corrected the bias, based on ERA-interim monthly climatology, and bias-corrected ERA-WRF fields were applied as a reference for calibration of both the historical and the projected CESM forced WRF (CESM-WRF) products. Biases, such as, a cold temperature bias during summer and a warm temperature bias during winter as well as a wet bias for annual precipitation that CESM holds over northern Alaska persisted in CESM-WRF runs. The linear scaling of CESM-WRF eventually produced high-resolution downscaling products for the Alaskan North Slope for hydrological and ecological research, together with the calibrated ERA-WRF run, and its capability extends far beyond that. Other climatic research has been proposed, including exploration of historical and projected climatic extreme events and their possible connections to low-frequency sea-atmospheric oscillations, as well as near-surface permafrost degradation and ice regime shifts of lakes. These dynamically downscaled, bias corrected climatic datasets provide improved spatial and temporal resolution data necessary for ongoing modeling efforts in northern Alaska focused on reconstructing and projecting hydrologic changes, ecosystem processes and responses, and permafrost thermal regimes. The dynamical downscaling methods presented in this study can also be used to create more suitable model input datasets for other sub-regions of the Arctic.

A suite of global, cross-scale topographic variables for environmental and biodiversity modeling, links to files in GeoTIFF format

Topographic variation, the spatial variation in elevation and terrain features, underpins a myriad of patterns and processes in geography and ecology and is key to understanding the variation of life on the planet. The characterization of this variation is scale-dependent, i.e. it varies with the distance over which features are assessed and with the spatial grain (grid cell resolution) of analysis. A fully standardized and global multivariate product of different terrain features has the potential to support many large-scale basic research and analytical applications, however to date, such technique is unavailable. Here we used the digital elevation model products of global 250 m GMTED and near-global 90 m SRTM to derive a suite of topographic variables: elevation, slope, aspect, eastness, northness, roughness, terrain roughness index, topographic position index, vector ruggedness measure, profile and tangential curvature, and 10 geomorphological landform classes. We aggregated each variable to 1, 5, 10, 50 and 100 km spatial grains using several aggregation approaches (median, average, minimum, maximum, standard deviation, percent cover, count, majority, Shannon Index, entropy, uniformity). While a global cross-correlation underlines the high similarity of many variables, a more detailed view in four mountain regions reveals local differences, as well as scale variations in the aggregated variables at different spatial grains. All newly-developed variables are available for download at http://www.earthenv.org and can serve as a basis for standardized hydrological, environmental and biodiversity modeling at a global extent.