Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

40Ar-39Ar plateau and isochron ages from the Ninetyeast and Broken Ridges (Table 1)

Concordant plateau and isochron ages are obtained from 40Ar-39Ar incremental heating experiments on volcanic rocks recovered by drilling at three Leg 121 sites along the Ninetyeast Ridge and two dredge locations on the southern scarp of the Broken Ridge, eastern Indian Ocean. The new data confirm a northerly increase in the age of volcanism along the Ninetyeast Ridge, from 38 to 82 Ma; this lineament links current hotspot volcanism near the Kerguelen islands with the Rajmahal flood basalt eruptions at M0 time (117 ± 1 Ma). The Broken Ridge was formed over the same hotspot at 88-89 Ma, but later experienced rift-related volcanism in Paleocene time (63 Ma). The geometry and distribution of ages along these prominent volcanic ridges and the Mascarene-Chagos-Laccadive-Maldive ridge system in the western Indian Ocean are most compatible with plate motions over fixed hotspots near Kerguelen and Reunion islands, respectively.

Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.

(Table 1) 40Ar/39Ar plateau and isochron ages for basalts from the Réunion hotspot track

Concordant plateau and isochron ages were calculated from 40Ar/39Ar incremental heating experiments on volcanic rocks recovered by drilling at four Leg 115 sites and two industry wells along the volcanic lineament connecting Reunion Island to the Deccan flood basalts, western Indian Ocean. The new ages provide unequivocal evidence that volcanic activity migrated southward along this sequence of linear ridges. The geometry and age distribution of volcanism are most compatible with origin above a stationary hotspot centered beneath Reunion. The hotspot became active with rapid eruption of the Deccan flood basalts, western India, and subsequent volcanic products record the northward motion of the Indian and African plates over the hotspot through Tertiary time. The radiometric ages are in general accord with basal biostratigraphic age estimates, although some adjustments in current magnetobiostratigraphic time scales may be required.

Time-series of seagrass and land cover maps from 1972 to 2010, in South East Queensland, Australia, with links to GeoTIFFs

Long term global archives of high-moderate spatial resolution, multi-spectral satellite imagery are now readily accessible, but are not being fully utilised by management agencies due to the lack of appropriate methods to consistently produce accurate and timely management ready information. This work developed an object-based remote sensing approach to map land cover and seagrass distribution in an Australian coastal environment for a 38 year Landsat image time-series archive (1972-2010). Landsat Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) imagery were used without in situ field data input (but still using field knowledge) to produce land and seagrass cover maps every year data were available, resulting in over 60 map products over the 38 year archive. Land cover was mapped annually using vegetation, bare ground, urban and agricultural classes. Seagrass distribution was also mapped annually, and in some years monthly, via horizontal projected foliage cover classes, sand and deep water. Land cover products were validated using aerial photography and seagrass maps were validated with field survey data, producing several measures of accuracy. An average overall accuracy of 65% and 80% was reported for seagrass and land cover products respectively, which is consistent with other studies in the area. This study is the first to show moderate spatial resolution, long term annual changes in land cover and seagrass in an Australian environment, created without the use of in situ data; and only one of a few similar studies globally. The land cover products identify several long term trends; such as significant increases in South East Queensland's urban density and extent, vegetation clearing in rural and rural-residential areas, and inter-annual variation in dry vegetation types in western South East Queensland. The seagrass cover products show that there has been a minimal overall change in seagrass extent, but that seagrass cover level distribution is extremely dynamic; evidenced by large scale migrations of higher seagrass cover levels and several sudden and significant changes in cover level. These mapping products will allow management agencies to build a baseline assessment of their resources, understand past changes and help inform implementation and planning of management policy to address potential future changes.

(Appendix) Recent surface areas and historical lengths of glaciers from Patagonia and Tierra del Fuego, South America

Fluctuations in the length of 72 glaciers in the Northern and Southern Patagonia Icefield (NPI and SPI, respectively) and the Cordillera Darwin Icefield (CDI) were estimated between 1945 and 2005. The information obtained from historical maps based on 1945 aerial photographs was compared to ASTER and Landsat satellite images and to information found in the literature. The majority of glaciers have retreated considerably, with maximum values of 12.2 km for Marinelli Glacier in the CDI, 11.6 km for O'Higgins Glacier in the SPI and 5.7 km for San Rafael Glacier in the NPI. Among the 20 glaciers that have retreated the most relative to their size, small (less than 50 km**2) and medium (between 50 and 200 km**2) glaciers are the most affected. However, no direct relation between glacier retreat and size was found for the 72 glaciers studied. The highest percentage retreat in the CDI was by the CDI-03 Glacier (37.9%) and Marinelli Glacier (37.6%). In the SPI, relative retreats were heterogeneous and fluctuated between 27.2% (Amelia Glacier) and 0.4% (Viedma Glacier). In the NPI, relative retreat was very high for Strindberg and Cachet glaciers (35.9% and 27.6%, respectively) but for the remaining glaciers in this icefield it ranged between 11.8% (Piscis Glacier) and 3.6% (San Quintin Glacier). In addition to surface area, the surface slope (calculated on the basis of the DEM SRTM) was also related to the relative retreat and no straightforward relation was found. From a global point of view, we suggest that glacier retreat in the region is controlled firstly by atmospheric warming, as it has been reported in this area. Besides the general increase in temperature observed, no signal of a geographical pattern for the fluctuations in glacier length was found. Consequently, glaciers appear to initially react to local conditions most probably induced by their exposition, geometry and hypsometry. The heterogeneity of rates of retreat suggests that differences in basin geometry, glacier dynamics and response time are key features to explain fluctuations of each glacier.

Major components and major and trace elements from ODP Site 175-1085

We have analyzed 33 Pliocene bulk sediment samples from Ocean Drilling Program Site 1085 in the Cape Basin, located offshore of western Africa in the Angola-Benguela Current system, for 17 major and trace elements, and interpreted their associations and temporal variations in the context of an allied data set of CaCO3, opal, and Corg. We base our interpretations on elemental ratios, accumulation rates, inter-element correlations, and several multi-element statistical techniques. On the basis of qualitative assessment of downhole changes in the distributions of P and Ba, utilized as proxies of export production, we conclude that highs in bulk and biogenic accumulation that occur at 3.2 Ma, 3.0 Ma, 2.4 Ma, and 2.25 Ma were caused by increases in export production as well as terrigenous flux, and record a greater sequestering of organic matter during these time periods. Studies of refractory elements and other indicator proxies (SiO2, Al2O3, TiO2, Fe2O3, MgO, V, Cr, Sr, and Zr) strongly suggest that the terrigenous component of the bulk sediment is composed of two compositional end-members, one being 'basaltic' in composition and the other similar to an 'average shale'. The basaltic end-member comprises approximately 10-15% of the total bulk sediment and its presence is consistent with the local geology of source material in the drainage basin of the nearby Orange River. The increase in bulk accumulation at 2.4 Ma appears to reflect a greater relative increase in basaltic input than the relative increase in shale-type input. Although studies such as this cannot precisely identify the transport mechanisms of the different terrigenous components, these results are most consistent with variations in sea level (and associated changes in shelf geometry and fluvial input) being responsible for the changing depositional conditions along the Angolan Margin during this time period.

Grain-size distribution, sedimentation rate and x-ray fluorescence data of four sediment cores off the Senegal River estuary

Fine-grained sediment depocenters on continental shelves are of increased scientific interest since they record environmental changes sensitively. A north-south elongated mud depocenter extends along the Senegalese coast in mid-shelf position. Shallow-acoustic profiling was carried out to determine extent, geometry and internal structures of this sedimentary body. In addition, four sediment cores were retrieved with the main aim to identify how paleoclimatic signals and coastal changes have controlled the formation of this mud depocenter. A general paleoclimatic pattern in terms of fluvial input appears to be recorded in this depositional archive. Intervals characterized by high terrigenous input, high sedimentation rates and fine grain sizes occur roughly contemporaneously in all cores and are interpreted as corresponding to intensified river discharge related to more humid conditions in the hinterland. From 2750 to 1900 and from 1000 to 700 cal a BP, wetter conditions are recorded off Senegal, an observation which is in accordance with other records from NW-Africa. Nevertheless, the three employed proxies (sedimentation rate, grain size and elemental distribution) do not always display consistent inter-core patterns. Major differences between the individual core records are attributed to sediment remobilization which was linked to local hydrographic variations as well as reorganizations of the coastal system. The Senegal mud belt is a layered inhomogeneous sedimentary body deposited on an irregular erosive surface. Early Holocene deceleration in the rate of the sea-level rise could have enabled initial mud deposition on the shelf. These favorable conditions for mud deposition occur coevally with a humid period over NW-Africa, thus, high river discharge. Sedimentation started preferentially in the northern areas of the mud belt. During mid-Holocene, a marine incursion led to the formation of an embayment. Afterwards, sedimentation in the north was interrupted in association with a remarkable southward shift in the location of the active depocenter as it is reflected by the sedimentary architecture and confirmed by radiocarbon dates. These sub-recent shifts in depocenters location are caused by migrations of the Senegal River mouth. During late Holocene times, the weakening of river discharge allowed the longshore currents to build up a chain of beach barriers which have forced the river mouth to shift southwards.

GPS raw data (control points and ground control points) from the Liguria Region, Borghetto Santo Spirito, Italy

Monitoring the impact of sea storms on coastal areas is fundamental to study beach evolution and the vulnerability of low-lying coasts to erosion and flooding. Modelling wave runup on a beach is possible, but it requires accurate topographic data and model tuning, that can be done comparing observed and modeled runup. In this study we collected aerial photos using an Unmanned Aerial Vehicle after two different swells on the same study area. We merged the point cloud obtained with photogrammetry with multibeam data, in order to obtain a complete beach topography. Then, on each set of rectified and georeferenced UAV orthophotos, we identified the maximum wave runup for both events recognizing the wet area left by the waves. We then used our topography and numerical models to simulate the wave runup and compare the model results to observed values during the two events. Our results highlight the potential of the methodology presented, which integrates UAV platforms, photogrammetry and Geographic Information Systems to provide faster and cheaper information on beach topography and geomorphology compared with traditional techniques without losing in accuracy. We use the results obtained from this technique as a topographic base for a model that calculates runup for the two swells. The observed and modeled runups are consistent, and open new directions for future research.