Digital surface model, hillshade and orthophoto of the world's largest fossil oyster reef, links to GeoTIFFs

The world's largest fossil oyster reef, formed by the giant oyster Crassostrea gryphoides and located in Stetten (north of Vienna, Austria) is studied by Harzhauser et al., 2015, 2016; Djuricic et al., 2016. Digital documentation of the unique geological site is provided by terrestrial laser scanning (TLS) at the millimeter scale. Obtaining meaningful results is not merely a matter of data acquisition with a suitable device; it requires proper planning, data management, and postprocessing. Terrestrial laser scanning technology has a high potential for providing precise 3D mapping that serves as the basis for automatic object detection in different scenarios; however, it faces challenges in the presence of large amounts of data and the irregular geometry of an oyster reef. We provide a detailed description of the techniques and strategy used for data collection and processing in Djuricic et al., 2016. The use of laser scanning provided the ability to measure surface points of 46,840 (estimated) shells. They are up to 60-cm-long oyster specimens, and their surfaces are modeled with a high accuracy of 1 mm. In addition to laser scanning measurements, more than 300 photographs were captured, and an orthophoto mosaic was generated with a ground sampling distance (GSD) of 0.5 mm. This high-resolution 3D information and the photographic texture serve as the basis for ongoing and future geological and paleontological analyses. Moreover, they provide unprecedented documentation for conservation issues at a unique natural heritage site.

Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Soil, snow, weather, and sub-surface storage data from a mountain catchment in the rain-snow transition zone

A comprehensive hydroclimatic data set is presented for the 2011 water year to improve understanding of hydrologic processes in the rain-snow transition zone. This type of dataset is extremely rare in scientific literature because of the quality and quantity of soil depth, soil texture, soil moisture, and soil temperature data. Standard meteorological and snow cover data for the entire 2011 water year are included, which include several rain-on-snow events. Surface soil textures and soil depths from 57 points are presented as well as soil texture profiles from 14 points. Meteorological data include continuous hourly shielded, unshielded, and wind corrected precipitation, wind speed, air temperature, relative humidity, dew point temperature, and incoming solar and thermal radiation data. Sub-surface data included are hourly soil moisture data from multiple depths from 7 soil profiles within the catchment, and soil temperatures from multiple depths from 2 soil profiles. Hydrologic response data include hourly stream discharge from the catchment outlet weir, continuous snow depths from one location, intermittent snow depths from 5 locations, and snow depth and density data from ten weekly snow surveys. Though it represents only a single water year, the presentation of both above and below ground hydrologic condition makes it one of the most detailed and complete hydro-climatic datasets from the climatically sensitive rain-snow transition zone for a wide range of modeling and descriptive studies.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Sea-surface temperature reconstruction of ODP Hole 175-1078C

A prominent feature in the Southeast Atlantic is the Angola-Benguela Front (ABF), the convergence between warm tropical and cold subtropical upwelled waters. At present, the sea-surface temperature (SST) gradient across the ABF and its position are influenced by the strength of southeasterly (SE) trade winds. Here, we present a record of changes in the ABF SST gradient over the last 25 kyr. Variations in this SST contrast indicate that periods of strengthened SE trade-wind intensity occurred during the Last Glacial Maximum, the Younger Dryas, and the Mid to Late Holocene, while Heinrich Event 1, the early part of the Bølling-Allerød, and the Early Holocene were periods of weakened SE trade-winds.

Investigation of the seasonality of N. pachyderma in the North Atlantic Ocean during Heinrich Stadial 1 using an ecosystem modeling approach

Fossil shells of planktonic foraminifera serve as the prime source of information on past changes in surface ocean conditions. Because the population size of planktonic foraminifera species changes throughout the year, the signal preserved in fossil shells is biased towards the conditions when species production was at its maximum. The amplitude of the potential seasonal bias is a function of the magnitude of the seasonal cycle in production. Here we use a planktonic foraminifera model coupled to an ecosystem model to investigate to what degree seasonal variations in production of the species Neogloboquadrina pachyderma may affect paleoceanographic reconstructions during Heinrich Stadial 1 (~18-15 cal. ka B.P.) in the North Atlantic Ocean. The model implies that during Heinrich Stadial 1 the maximum seasonal production occurred later in the year compared to the Last Glacial Maximum (~21-19 cal. ka B.P.) and the pre-industrial era north of 30 ºN. A diagnosis of the model output indicates that this change reflects the sensitivity of the species to the seasonal cycle of sea-ice cover and food supply, which collectively lead to shifts in the modeled maximum production from the Last Glacial Maximum to Heinrich Stadial 1 by up to six months. Assuming equilibrium oxygen isotopic incorporation in the shells of N. pachyderma, the modeled changes in seasonality would result in an underestimation of the actual magnitude of the meltwater isotopic signal recorded by fossil assemblages of N. pachyderma wherever calcification is likely to take place.

(Tables 2, 3) Element ratios of planktonic foraminifera from western Pacific surface sediments

Accurately quantifying deep-sea calcite dissolution is crucial for understanding the role of the marine carbonate system in regulating atmospheric pCO2 over millennia. We compare a foraminifer-fragmentation-based calcite dissolution proxy (Globorotalia menardii fragmentation index (MFI)) to Mg/Ca, Sr/Ca, and Mg/Sr in several species of deep dwelling planktonic foraminifers. We conducted microfossil and geochemical analyses on the same core top samples taken at different depths on the Ontong Java Plateau to maximize the dissolution signal and minimize the temperature overprint on our data. We also compare elemental ratios from planktonic foraminifer tests to modern bottom water [CO3]2- undersaturation and model-derived estimates of percent calcite dissolved in deep-sea sediments. We find clear linear decreases in Mg/Ca or Mg/Sr in G. menardii and Pulleniatina obliquiloculata with increasing (1) bottom water [CO3]2- undersaturation, (2) percent calcite dissolved in sediments calculated with biogeochemical modeling, (3) MFI, and (4) percent calcite dissolved derived from MFI. These findings lend further support to MFI as a calcite dissolution proxy for deep-sea sediments. In contrast, we find no significant correlation between Sr/Ca and independent dissolution indicators. Our results suggest that Mg/Ca and Mg/Sr from deep dwelling foraminifers could potentially be used as calcite dissolution proxies in combination with independent water temperature estimates. Likewise, establishing the relationship between MFI and dissolution-induced changes in the Mg/Ca of surface-dwelling foraminifers could provide a tool to correct Mg/Ca-derived sea surface temperature reconstructions for calcite dissolution.

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.

Basic measurements of radiation from the Baseline Surface Radiation Network (BSRN) in the years 1992 to 2012, reference list of 6378 datasets

In the framework of the global energy balance, the radiative energy exchanges between Sun, Earth and space are now accurately quantified from new satellite missions. Much less is known about the magnitude of the energy flows within the climate system and at the Earth surface, which cannot be directly measured by satellites. In addition to satellite observations, here we make extensive use of the growing number of surface observations to constrain the global energy balance not only from space, but also from the surface. We combine these observations with the latest modeling efforts performed for the 5th IPCC assessment report to infer best estimates for the global mean surface radiative components. Our analyses favor global mean downward surface solar and thermal radiation values near 185 and 342 Wm**-2, respectively, which are most compatible with surface observations. Combined with an estimated surface absorbed solar radiation and thermal emission of 161 Wm**-2 and 397 Wm**-2, respectively, this leaves 106 Wm**-2 of surface net radiation available for distribution amongst the non-radiative surface energy balance components. The climate models overestimate the downward solar and underestimate the downward thermal radiation, thereby simulating nevertheless an adequate global mean surface net radiation by error compensation. This also suggests that, globally, the simulated surface sensible and latent heat fluxes, around 20 and 85 Wm**-2 on average, state realistic values. The findings of this study are compiled into a new global energy balance diagram, which may be able to reconcile currently disputed inconsistencies between energy and water cycle estimates.

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.