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.

Milankovitch frequency peaks of ODP Site 113-693 (Table 3)

The first well logs collected below the Antarctic circle were obtained during Leg 113 at Site 693 on the Dronning Maud Land Margin (Antarctica) in the Weddell Sea. Gamma-ray, resistivity, and sonic logs were collected between 108.0 and 439.0 mbsf. The downhole logs show good agreement with the data collected from cores and provide a continuous measurement of the sedimentary record. These continuous log records show that the rather uniform Tertiary lithology seen in cores is characterized by high-frequency variability in the log data. Several thin hard streaks are identified, the largest of which coincides with a major Miocene hiatus. Associated with this hiatus is a change to lower illite content (and correspondingly lower gamma-ray counts) and to a significant increase in diatom content. Spectral analysis of the logs was performed on the lower Pliocene through upper Oligocene interval (108.0-343.0 mbsf). Between 108.0 and 245.0 mbsf, average sedimentation rates (50 and 26 m/m.y.) are high enough to show that variance is present in the orbital eccentricity (~95 k.y.) and obliquity (~41 k.y.) bands. Between 253.0 and 343.0 mbsf, the sedimentation rate (8 m/m.y.) is too low to resolve high frequency variations. The Milankovitch frequencies are best developed in the resistivity logs. Resistivity is responding to changes in porosity, which in these sediments is controlled by the abundance of biosiliceous sediments, particularly diatoms. The orbital forcing suggested by the Milankovitch frequencies may be influencing diatom productivity by inducing oscillations in upwelling, ice coverage, pack ice, and/or polynya. Although variations in diatom abundance were observed in the cores, they were not attributed to a Milankovitch signal, and therefore in this environment, downhole logs are an important contribution to the detection and understanding of orbitally influenced changes in sedimentation.

Sediment echosounding data files from the Amundsen Sea

The flow of ice streams, which account for most discharge from large ice sheets, is controlled by processes operating at their bed. Data from modern ice stream beds are difficult to obtain, but where ice advanced onto continental shelves during glacial periods extensive areas of the former bed can be imaged using modern swath sonar tools. We present new multibeam swath bathymetry data analyzed alongside sparse pre-existing data from the Amundsen Sea Embayment. The compilation is the most extensive, continuous area of multibeam data coverage yet obtained on the inner continental shelf of Antarctica. The data reveal streamlined subglacial bedforms that define a zone of paleo-ice stream convergence but, in contrast to previous models, do not show a simple down-flow progression of bedform types along paleo-ice stream troughs. We interpret high spatial variability of bedforms as indicating a complex mechanical and hydrodynamic regime at the former ice stream beds, consistent with observations from some modern ice streams. We conclude that care must be taken when using bedforms to infer paleo-ice stream velocities.

Effects of in situ CO2 enrichment on structural characteristics, photosynthesis, and growth of the Mediterranean seagrass Posidonia oceanica

Seagrass is expected to benefit from increased carbon availability under future ocean acidification. This hypothesis has been little tested by in situ manipulation. To test for ocean acidification effects on seagrass meadows under controlled CO2/pH conditions, we used a Free Ocean Carbon Dioxide Enrichment (FOCE) system which allows for the manipulation of pH as continuous offset from ambient. It was deployed in a Posidonia oceanica meadow at 11 m depth in the Northwestern Mediterranean Sea. It consisted of two benthic enclosures, an experimental and a control unit both 1.7 m**3, and an additional reference plot in the ambient environment (2 m**2) to account for structural artifacts. The meadow was monitored from April to November 2014. The pH of the experimental enclosure was lowered by 0.26 pH units for the second half of the 8-month study. The greatest magnitude of change in P. oceanica leaf biometrics, photosynthesis, and leaf growth accompanied seasonal changes recorded in the environment and values were similar between the two enclosures. Leaf thickness may change in response to lower pH but this requires further testing. Results are congruent with other short-term and natural studies that have investigated the response of P. oceanica over a wide range of pH. They suggest any benefit from ocean acidification, over the next century (at a pH of 7.7 on the total scale), on Posidonia physiology and growth may be minimal and difficult to detect without increased replication or longer experimental duration. The limited stimulation, which did not surpass any enclosure or seasonal effect, casts doubts on speculations that elevated CO2 would confer resistance to thermal stress and increase the buffering capacity of meadows.

Geochemistry and descriptions of manganese nodules and crusts retrieved from the open ocean

The book is a compilation of all available data at the time of publication (1965) on the subject of marine minerals together with the author's original ideas regarding their exploitation. It is one of the most significant publications on ocean resources. It is particularly focused on manganese deposits, their description, sedimentary setting, formation and geochemistry.

Benthic and substrate cover data derived from photo-transect surveys in Lizard Island Reef conducted on December 10-15, 2011

Underwater georeferenced photo-transect surveys were conducted on December 10-15, 2011 at various sections of the reef at Lizard Island, Great Barrier Reef. For this survey a snorkeler or diver swam over the bottom while taking photos of the benthos at a set height using a standard digital camera and towing a GPS in a surface float which logged the track every five seconds. A standard digital compact camera was placed in an underwater housing and fitted with a 16 mm lens which provided a 1.0 m x 1.0 m footprint, at 0.5 m height above the benthos. Horizontal distance between photos was estimated by three fin kicks of the survey diver/snorkeler, which corresponded to a surface distance of approximately 2.0 - 4.0 m. The GPS was placed in a dry-bag and logged the position as it floated at the surface while being towed by the photographer. A total of 5,735 benthic photos were taken. A floating GPS setup connected to the swimmer/diver by a line enabled recording of coordinates of each benthic photo (Roelfsema 2009). Approximation of coordinates of each benthic photo was conducted based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software (www.geospatialexperts.com). Coordinates of each photo were interpolated by finding the GPS coordinates that were logged at a set time before and after the photo was captured. Benthic or substrate cover data was derived from each photo by randomly placing 24 points over each image using the Coral Point Count for Microsoft Excel program (Kohler and Gill, 2006). Each point was then assigned to 1 of 78 cover types, which represented the benthic feature beneath it. Benthic cover composition summary of each photo scores was generated automatically using CPCE program. The resulting benthic cover data of each photo was linked to GPS coordinates, saved as an ArcMap point shapefile, and projected to Universal Transverse Mercator WGS84 Zone 55 South.

Benthic and substrate cover data derived from photo-transect surveys in Lizard Island, Great Barrier Reef conducted on October 3-7, 2012

Underwater georeferenced photo-transect surveys were conducted on October 3-7, 2012 at various sections of the reef and lagoon at Lizard Island, Great Barrier Reef. For this survey a snorkeler swam while taking photos of the benthos at a set distance from the benthos using a standard digital camera and towing a GPS in a surface float which logged the track every five seconds. A Canon G12 digital camera was placed in a Canon underwater housing and photos were taken at 1 m height above the benthos. Horizontal distance between photos was estimated by three fin kicks of the survey snorkeler, which corresponded to a surface distance of approximately 2.0 - 4.0 m. The GPS was placed in a dry bag and logged the position at the surface while being towed by the photographer (Roelfsema, 2009). A total of 1,265 benthic photos were taken. Approximation of coordinates of each benthic photo was conducted based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software (www.geospatialexperts.com). Coordinates of each photo were interpolated by finding the GPS coordinates that were logged at a set time before and after the photo was captured. Benthic or substrate cover data was derived from each photo by randomly placing 24 points over each image using the Coral Point Count for Microsoft Excel program (Kohler and Gill, 2006). Each point was then assigned to 1 of 79 cover types, which represented the benthic feature beneath it. Benthic cover composition summary of each photo scores was generated automatically using CPCE program. The resulting benthic cover data of each photo was linked to GPS coordinates, saved as an ArcMap point shapefile, and projected to Universal Transverse Mercator WGS84 Zone 55 South.

Sedimentology of the northwestern Weddel Sea continental shelf and slope

Geologie cores on two profiles oriented normaly to the continental shelf and slope, have been investigated to reconstruct the Quaternary sedimentary history of the southeast continental border of South Orkney (NW Weddell Sea). The sediments were described macroscopically and their fabric investigated by use of X-radiographs. Laboratory work comprised detailed grain-size analysis, determination of the watercontent, carbonate, organic carbon and sand fraction.composition. Stable oxygen and carbon isotopes have been measured On planktonic foraminifera. Palaeomagnetism, analysis of 230Th-content and detailed comparison of the lithlogic Parameters with the oxygen isotope stages (Martinson curve) were used for stratigraphic classification of the sediments. The sediment cores from the continental slope comprise a maximum age of 300,000 years B. P.. Bottom currents, ice rafting and biogenic input are the main sources of sediment. Based on lithologic parameters a distinction between glacial and interglacial facies is possible. Silty clays without microfossils and few bioturbation characterise the sediments of the glacial facies. Only small amounts of icerafted debris can be recognized. This type of sediment was accumulated during times of lower sea-level and drastically reduced rate of bottom water production. Based on grain-size distribution, bottom current velocities of 0.01 cmls were calculated. Thick sea-ice coverage reduced biogenic production in the surface water, and as consequence benthic communities were depleted. Because of the reduced benthic life, sediments are only slithly bioturbated. At the beginning of the interglacial Stage, the sea-level rised rapidly, and calving rate of icebergs, combined with input of ice-rafted material, increased considerably. Sediments of this transition facies are silty cliiys with a high proportion of coarse ice-rafted debris, but without microfossils. With the onset of bottom water production in connection with shelf ice water, sediments of interglacial facies were formed. They consist of silty clays to clayey silts with considerable content of sand and gravel. Sediments are strongly bioturbated. Based On the sediment caracteristics, current velocities of the bottom water were calculated to be of 0.96 cmls for interglacials. At the southern slope of a NW/SE-striking ridge, bottom water current is channelized, resulting in a drastic increase of current velocities. Current velocities up to 7.5 cm/s lead to formation of residual sediments. While the continental slope has predominantly fine sediments, the South Orkney shelf are mainly sandy silts and silty sands with a high proportion of gravel. These sediments were formed dominantly by ice-rafting during Brunhes- and Matuyama-Epoch. Currents removed the fine fraction of the sediments. Based on microfossil contents it was not possible to differentiate sediments from glacial to interglacial. In the upper Parts of the cores graded sequences truncated by erosion were observed. These sequences were formed during Brunhes-Epoch by strong currents with velocities decreasing periodically from about 7.5 cm/s to about 1 cm/s. Sediments with a high proportion of siliceous microfossils but barren of foraminifera compose the lower part of the shelf cores. These sediments have formed during the warmer Matuyama-Epoch.

Investigations on sediment coring site PS1823 in the Riiser-Larsen Sea, Southern Ocean

Physical and sedimentological investigations were carried out on a 14 m long gravity core and a 0.5 m long box core from 4440 m water depth off Queen Maud Land, East-Antarctica. Strongly bioturbated hemipelagic muds of predominantly terrigenous origin and a very small biogenic part build up the 'Normal-Facies'. Several sandy to silty layers are inserted in the 'Normal-Facies'. These layers are seperated by lithology, structure and the investigated parameters of this study and are interpreted as turbidites. The source area for the turbidity currents is supposed to be at the uppermost continental margin, close to the shelf break and there is evidenee for this gravity transport within the erosive Ritscher-Canyon, which extends close to the core position. The distribution of biogenic components indicates an age of 1.3 million years or more, with an average sedimentation rate of about 1 cm/1000 years. Early diagenetic proeesses caused water loss by compaction, errosion and dissolution of biogenic components and precipitation and recrystallization of manganese micronodules. Cyclic fluctuations of the sediment-parameters within the 'Normal-Facies' enable the distinction of a 'Glazial'- and an 'Interglazial'-Facies. The 'Glazial'-Facies reflects glacial sedimentary conditions and shows a dark olive gray colour, high susceptibility, low silt/clay-ratios, only a few biogenic components and the regular occurence of interrelated turbidite layers. In contrast, the 'Interglazial'-Facies is dominated by a light olive or olive-brown colour, low susceptibility, high silt/clay-ratios and an increased number of biogenic components. This facies corresponds to interglacial conditions. Three main processes are supposed to have been responsible for the observed facies changes: (1) the bottom water mass circulation, (2) the gravity transport by turbidity currents and (3) the biogenic surface production. These processes are related to the quaternary climatic changes. The extension of the ice shelves directed the gravity transport to the deep sea and the formation of Antarctic Bottom Water, which in turn influenced the silt/clay-ratios in the sediment record. Fluctuations in sea ice coverage controlled the biogenic surface production.

Hide, survive and hitch a ride, the dispersal of Antarctic krill into the nursery habitats

Antarctic krill (Euphausia superba), a key species of Southern Ocean food webs plays a central role in ecosystem processes, community dynamics of apex predators and as a commercial fishery target. A decline in krill abundance during the late 20th century in the SW Atlantic sector has been linked to a concomitant decrease in sea ice, based on the hypothesis that sea ice acts as a feeding ground for overwintering larvae. However, evidence supporting this hypothesis has been scarce due to logistical challenges of collecting data in austral winter. Here we report on a winter study that involved diver observations of larval krill in their under-ice environment, ship-based studies of krill, sea ice physical characteristics, and biophysical model analyses of krill-ocean-ice interactions. We present evidence that complex under-ice topography is vital for larval krill in terms of dispersal and advection into high productive nursery habitats, rather than the provision by the ice environment of food. Further, ongoing changes in sea ice will lead to increases in sea-ice regimes favourable for overwintering larval krill but shifting southwards. This will result in ice-free conditions in the SW Atlantic, which will be conducive for enhancing food supplies due to sufficient light and iron availability, thus enhancing larvae development and growth. However, the associated impact on dispersal and advection may lead to a net shift in krill from the SW Atlantic to regions further east by the eastward flowing ACC and the northern branch of the Weddell Gyre, with profound consequences for the Southern Ocean pelagic ecosystem.

Seawater carbonate chemistry, and Lithophyllum sp. and Feldmannia spp. coverage and dry weight during experiments, 2009

Climate-driven change represents the cumulative effect of global through local-scale conditions, and understanding their manifestation at local scales can empower local management. Change in the dominance of habitats is often the product of local nutrient pollution that occurs at relatively local scales (i.e. catchment scale), a critical scale of management at which global impacts will manifest. We tested whether forecasted global-scale change [elevated carbon dioxide (CO2) and subsequent ocean acidification] and local stressors (elevated nutrients) can combine to accelerate the expansion of filamentous turfs at the expense of calcifying algae (kelp understorey). Our results not only support this model of future change, but also highlight the synergistic effects of future CO2 and nutrient concentrations on the abundance of turfs. These results suggest that global and local stressors need to be assessed in meaningful combinations so that the anticipated effects of climate change do not create the false impression that, however complex, climate change will produce smaller effects than reality. These findings empower local managers because they show that policies of reducing local stressors (e.g. nutrient pollution) can reduce the effects of global stressors not under their governance (e.g. ocean acidification). The connection between research and government policy provides an example whereby knowledge (and decision making) across local through global scales provides solutions to some of the most vexing challenges for attaining social goals of sustainability, biological conservation and economic development.

Meteorological observations and eddy covariance raw data from polygonal tundra in the Lena River Delta, Siberia

We present the first ecosystem-scale methane flux data from a northern Siberian tundra ecosystem covering the entire snow-free period from spring thaw until initial freeze-back. Eddy covariance measurements of methane emission were carried out from the beginning of June until the end of September in the southern central part of the Lena River Delta (72°22' N, 126°30' E). The study site is located in the zone of continuous permafrost and is characterized by Arctic continental climate with very low precipitation and a mean annual temperature of -14.7°C. We found relatively low fluxes of on average 18.7 mg/m**2/d, which we consider to be because of (1) extremely cold permafrost, (2) substrate limitation of the methanogenic archaea, and (3) a relatively high surface coverage of noninundated, moderately moist areas. Near-surface turbulence as measured by the eddy covariance system in 4 m above the ground surface was identified as the most important control on ecosystem-scale methane emission and explained about 60% of the variance in emissions, while soil temperature explained only 8%. In addition, atmospheric pressure was found to significantly improve an exponential model based on turbulence and soil temperature. Ebullition from waterlogged areas triggered by decreasing atmospheric pressure and near-surface turbulence is thought to be an important pathway that warrants more attention in future studies. The close coupling of methane fluxes and atmospheric parameters demonstrated here raises questions regarding the reliability of enclosure-based measurements, which inherently exclude these parameters.