Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.

Model - data comparison of western African rainfall evolutions during the Last Interglacial (130-115 ka)

This dataset characterizes the evolution of western African precipitation indicated by marine sediment geochemical records in comparison to transient simulations using CCSM3 global climate model throughout the Last Interglacial (130-115 ka). It contains (1) defined tie-points (age models), newly published stable isotopes of benthic foraminifera and Al/Si log-ratios of eight marine sediment cores from the western African margin and (2) annual and seasonal rainfall anomalies (relative to pre-industrial values) for six characteristic latitudinal bands in western Africa simulated by CCSM3 (two transient simulations: one non-accelerated and one accelerated experiment).

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Sediment volume variations inferred from DGPS measurements within Bossons glacier proglacial area (Mont-Blanc massif, France)

This dataset presents Differential Global Positioning System data (DGPS) acquired within the Bossons glacier proglacial area. Bossons glacier is a rapidly retreating glacier and its proglacial area is deglaciated for ~30 years. Bossons stream is one of the outlets of the subglacial drainage system. It starts as a 800 m steep cascade reach, then flows through an area with gentler slope : the Plan des Eaux (PdE). PdE is a 300 m long, 50 m wide proglacial alluvial plain with an increasing channel mobility in the downstream direction but decreasing slope gradient and incision. As it may act a sediment trap, studying periglacial and proglacial erosion processes in the Bossons catchment requires to quantify PdE sediment volume evolution. A several meter-sized block located within Bossons proglacial area was set up as GPS base : its location was measured by one antenna (Topcon Hyper Pro) by performing 600 consecutive measurements throughout one day. A second antenna (Topcon Hyper Pro) was then used to measure XYZ location of points in the proglacial area with a ~2 m grid. Radio communication between the two antennas allowed differential calculations to be automatically carried out on field using the Topcon FC-250 hand controller. This methodology yields 3 cm XY and 1.5 cm Z uncertainties. DGPS data have been acquired through 10 campaigns from 2004 to 2014; campaigns from 2004 to 2008 cover a smaller area than those from 2010 to 2014. Digital Elevation Model (DEM) have been interpolated from DGPS data and difference between two DEMs yields deposited and eroded volume within PdE. Maps of PdE volume variation between two campaigns show that incision mainly occurs in the upper and lower sections where as deposition dominates in the middle section. Deposition, denudation and net rate (deposition rate - denudation rate) are calculated by normalizing volumes by DEM areas. Deposition dominates results with a mean net rate of 29 mm/yr. However, strong inter-annual variability exists and some years are dominated by denudation : -36 mm/yr and -100 mm/yr for 2006 and 2011, respectively. Nonetheless, oldest campaigns (2004 to 2008) were carried out on the lower part part of the alluvial plain and ruling them out to keep only complete DEM (2010 to 2014) yields a mean net rate of ~15 mm/yr. This results is coherent with field observations of both strong deposition (e.g. flood deposits) and strong erosion (e.g. 30 cm incision) evidences. Bossons glacier proglacial area is thus dynamic with year-to-year geormorphological changes but may leans toward increasing its mean elevation through a deposition dominated system.

Age model and color measurements of sediment cores from the eastern equatorial Pacific

High-resolution records of glacial-interglacial variations in biogenic carbonate, opal, and detritus (derived from non-destructive core log measurements of density, P-wave velocity and color; r >= 0.9) from 15 sediment sites in the eastern equatorial (sampling resolution is ~1 kyr) clear response to eccentricity and precession forcing. For the Peru Basin, we generate a high-resolution (21 kyr increment) orbitally-based chronology for the last 1.3 Ma. Spectral analysis indicates that the 100 kyr cycle became dominant at roughly 1.2 Ma, 200-300 kyr earlier than reported for other paleoclimatic records. The response to orbital forcing is weaker since the Mid-Brunhes Dissolution Event (at 400 ka). A west-east reconstruction of biogenic sedimentation in the Peru Basin (four cores; 91-85°W) distinguishes equatorial and coastal upwelling systems in the western and eastern sites, respectively. A north-south reconstruction perpendicular to the equatorial upwelling system (11 cores, 11°N-°3S) shows high carbonate contents (>= 50%) between 6°N and 4°S and highly variable opal contents between 2°N and 4°S. Carbonate cycles B-6, B-8, B-10, B-12, B-14, M-2, and M-6 are well developed with B-10 (430 ka) as the most prominent cycle. Carbonate highs during glacials and glacial-interglacial transitions extended up to 400 km north and south compared to interglacial or interglacial^glacial carbonate lows. Our reconstruction thus favors glacial-interglacial expansion and contraction of the equatorial upwelling system rather than shifting north or south. Elevated accumulation rates are documented near the equator from 6°N to 4°S and from 2°N to 4°S for carbonate and opal, respectively. Accumulation rates are higher during glacials and glacial-interglacial transitions in all cores, whereas increased dissolution is concentrated on Peru Basin sediments close to the carbonate compensation depth and occurred during interglacials or interglacial-glacial transitions.

Gotland Deep (Baltic Sea) independent age model based on Palaeomagnetic Secular Variation (PSV) and lead (Pb) deposition profiles

Dating of sediment cores from the Baltic Sea has proven to be difficult due to uncertainties surrounding the 14C reservoir age and a scarcity of macrofossils suitable for dating. Here we present the results of multiple dating methods carried out on cores in the Gotland Deep area of the Baltic Sea. Particular emphasis is placed on the Littorina stage (8 ka ago to the present) of the Baltic Sea and possible changes in the 14C reservoir age of our dated samples. Three geochronological methods are used. Firstly, palaeomagnetic secular variations (PSV) are reconstructed, whereby ages are transferred to PSV features through comparison with varved lake sediment based PSV records. Secondly, lead (Pb) content and stable isotope analysis are used to identify past peaks in anthropogenic atmospheric Pb pollution. Lastly, 14C determinations were carried out on benthic foraminifera (Elphidium spec.) samples from the brackish Littorina stage of the Baltic Sea. Determinations carried out on smaller samples (as low as 4 µg C) employed an experimental, state-of-the-art method involving the direct measurement of CO2 from samples by a gas ion source without the need for a graphitisation step - the first time this method has been performed on foraminifera in an applied study. The PSV chronology, based on the uppermost Littorina stage sediments, produced ten age constraints between 6.29 and 1.29 cal ka BP, and the Pb depositional analysis produced two age constraints associated with the Medieval pollution peak. Analysis of PSV data shows that adequate directional data can be derived from both the present Littorina saline phase muds and Baltic Ice Lake stage varved glacial sediments. Ferrimagnetic iron sulphides, most likely authigenic greigite (Fe3S4), present in the intermediate Ancylus Lake freshwater stage sediments acquire a gyroremanent magnetisation during static alternating field (AF) demagnetisation, preventing the identification of a primary natural remanent magnetisation for these sediments. An inferred marine reservoir age offset (deltaR) is calculated by comparing the foraminifera 14C determinations to a PSV & Pb age model. This deltaR is found to trend towards younger values upwards in the core, possibly due to a gradual change in hydrographic conditions brought about by a reduction in marine water exchange from the open sea due to continued isostatic rebound.

Stable carbon and oxygen isotope ratios and age model of ODP Hole 182-1127B

The aim of this study was to evaluate the potential of constructing an oxygen and carbon isotope stratigraphy for the late Pleistocene succession from Hole 1127B drilled on the Great Australian Bight. Stable isotope analyses were performed on bulk- and fine-fraction (<38 µm) sediment samples. The oxygen isotope variations are generally smaller in magnitude than expected from global pelagic records. This is most likely due to the neriticly dominated sediment composition. Correlation of the oxygen isotope data with carbonate mineralogy and downhole logging data shows simultaneous variations and trends, which are particularly evident in the mid-Pleistocene sediments. Correlation of the oxygen isotope data with the classic SPECMAP curve is used to evaluate the stratigraphic potential of the Site 1127 sediments. This study indicates that an isotope stratigraphy based on planktonic and benthic foraminifers is needed to fully evaluate the response of cool-water carbonates deposited in a margin setting to global ice-volume fluctuations and, hence, the associated sea level variations.

Calcareous nannofossil datums and control points for the TCMS age model of ODP Site 175-1085

An astronomically calibrated age model for the Pliocene section of Ocean Drilling Program Leg 175 Cape Basin Site 1085 based on magnetic susceptibility data was developed using shipboard biostratigraphic datums. The composite core magnetic susceptibility record was compiled using shipboard correlations between Holes 1085A and 1085B and then tuned to the record of orbital variations in eccentricity to generate an orbitally tuned age model. Magnetic susceptibility apparently records climate variations in the Cape Basin. Strong power spectra values at the 100- and 400-k.y. frequency suggest an orbital control on the beat of Pliocene climate change in the Cape Basin.

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.

Downhole variations in grain sizes at DSDP/ODP Hole 504B

The maximum grain sizes of plagioclase and magnetite in the groundmass of the sheeted dike complex drilled at Hole 504B have been measured. Downhole variations through a 440-m-long section show a crude zig-zag pattern consisting of a gradual decrease or increase followed by an abrupt jump. The gradual decrease or increase in grain size extends over many lithologic units, and hence, does not reflect variations in grain size within a single dike. Such a zig-zag pattern is well explained by grain-size variations through multiple dikes. By using the observed inclination of sheeted dikes of 81° ± 2.5°, thickness of the multiple dikes varies from 0.7 to 8.5 m and averages to 4 ± 1 m. The average thickness of individual dikes forming multiple dikes is 0.8 m. We expect such multiple dikes to be formed during rifting events beneath mid-oceanic spreading ridges. If the average expansion at rifting episodes is twice as wide as the average width of the multiple dike units, the full spreading rate of 7.2 cm/yr of Cocos Ridge gives 112 ± 33 yr for a time interval of the rifting. A simple one-dimensional conductive cooling model is applied to solidification of multiple dikes. Numerical simulations show that the grain-size variations observed through the drill hole are more consistent with a model where a new injection of a dike occurs periodically with a constant time interval rather than one where the next dike intrudes just after the solidification of the previous one. Grain-size variations within simple dikes from Iritono, Japan, and those for Makaopuhi lava lake, Hawaii, show that square root of crystallization time is linearly correlated with the logarithm of plagioclase size. By using an empirically derived relationship between these two variables, the variations of plagioclase size through Hole 504B are directly compared with the calculated times for crystallization. Each rifting episode at the Costa Rica Rift lasts for several years, and periodic injection of a new dike occurs into the center of a previously solidified multiple dike at time intervals varying from 1 to 12 months.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Age model and grain size distribution of the fine fraction of ODP Site 114-704

The physical properties of sediments beneath an upwelling area in the southern part of the Atlantic Ocean (ODP Hole 704A) were investigated. Highly significant correlations characterize the relationship of carbonate content to bulk density (R = 0.85), carbonate content to porosity (R = 0.84), and carbonate content to impedance (R = 0.84). No relationship exists between carbonate content and compressional-wave velocity (R = 0.24), indicating that amplitude variations in impedance are primarily controlled by variations in bulk density, which, in turn, are controlled by climatically driven biogenic opal and carbonate deposition. In general, maxima in impedance correspond to maxima in carbonate content (minima in opal content). The impedance record exhibits its most drastic change at about 2.4 Ma, marking dramatic increases in the average content of biogenic opal and the beginning of large-amplitude fluctuations. Between 0.7 and 0.4 Ma carbonate content, bulk density, and grain density decrease while opal content drastically increases. Similar changes have been observed in sediments beneath an upwelling cell off northwest Africa, indicating an oceanwide enhancement in upwelling or in the calcite corrosiveness of bottom water that appears to be isochronous.

Trace metal concentration and age model of ODP Hole 165-1002C sediments

Titanium and iron concentration data from the anoxic Cariaco Basin, off the Venezuelan coast, can be used to infer variations in the hydrological cycle over northern South America during the past 14,000 years with subdecadal resolution. Following a dry Younger Dryas, a period of increased precipitation and riverine discharge occurred during the Holocene 'thermal maximum'. Since ~5400 years ago, a trend toward drier conditions is evident from the data, with high-amplitude fluctuations and precipitation minima during the time interval 3800 to 2800 years ago and during the 'Little Ice Age'. These regional changes in precipitation are best explained by shifts in the mean latitude of the Atlantic Intertropical Convergence Zone (ITCZ), potentially driven by Pacific-based climate variability. The Cariaco Basin record exhibits strong correlations with climate records from distant regions, including the high-latitude Northern Hemisphere, providing evidence for global teleconnections among regional climates.