A new modern analogue reference dataset and its application to the 430-kyr pollen record from Lake Biwa

This study presents a newly compiled dataset of modern pollen and climate data from 798 sites across Japan and the Russian Far East. This comprehensive reference dataset combined with the modern analogue technique (MAT) provides a powerful tool for pollen-based reconstruction of the Quaternary Northwest Pacific climate. Pollen-derived reconstruction of the modern climate at the reference pollen-sampling sites matches well with the estimated modern climate values (R2 values vary between 0.79 and 0.95, and RMSEP values vary between 5.8 and 9.7% of the modern climatic range for all nine tested variables). The successful testing of the method encourages its application to the fossil pollen records. We used a coarse-resolution pollen record from Lake Biwa to reconstruct glacial-interglacial climate dynamics in central Japan since ~438 kyr and compared it to the earlier reconstruction based on a less representative reference dataset. The current and earlier results consistently demonstrate that the coldest glacial intervals experienced pronounced cooling in winter and moderate cooling in summer, supporting the growth of cool mixed forest (COMX) where warm mixed forest (WAMX) predominates today. During the last glacial, maximum (~24 kyr BP) mean temperatures of the coldest (MTCO) and warmest (MTWA) month were about -13 °C (RMSEP = 2.34 °C) and 21 °C (RMSEP = 1.66 °C) respectively, and annual precipitation (PANN) was about 800 mm (RMSEP = 158.06 mm). During the thermal optimums of the interglacial intervals, the temperatures of the coldest and warmest month were above 0 °C and 25 °C respectively, leading to the reconstruction of WAMX and temperate conifer forest (TECO). Although both these vegetation types grow in the southern part of Japan today, WAMX requires warmer space. The presence of WAMX during marine isotope stages (MIS) 11 and 1, and its absence during MIS 9 and MIS 5 contradict the marine isotope and Antarctic ice records, suggesting that the latter two interglacials were the warmest of the last 800 kyr. The apparent contradiction allows at least three different explanations including low temporal resolution of the pollen record; different trends in CO2 concentrations during 'short' and 'long' interglacials; and regional climate variability and non-linear response of different regions to the global forcing. More definitive conclusions will be possible on the basis of forthcoming high-resolution pollen records from central Japan.

Reference datasets and pollen record from sediment core YAK, Yakumo, SW Hokkaido

This study reconstructs middle and late Holocene vegetation and climate dynamics in the Oshima Peninsula, SW Hokkaido, using the published method of biome reconstruction and modern analogue technique applied to the Yakumo pollen record (42°17'03''N, 140°15'34''E) spanning the last 5500 years. Two previously published matrices assigning Japanese plant/pollen taxa to the major vegetation types (biomes) are tested using a newly compiled dataset of 78 surface pollen spectra from Hokkaido. With both matrices showing strengths and weaknesses in reconstructing cool mixed and temperate deciduous forests of Hokkaido, the results suggest the necessity to consider the whole list of identified terrestrial pollen taxa for generating robust vegetation reconstructions for northern Japan. Applied to the fossil pollen data, both biome-reconstruction approaches demonstrate consistently that oak-dominated cool mixed forest spread in the study region between 5.5 and 3.6 cal ka BP and was subsequently replaced by beech-dominated temperate deciduous forest. The pollen-based climate reconstruction suggests this change in the vegetation composition was caused by a shift from cooler and drier than present climate to warmer and wetter, similar to modern conditions about 3.6 cal ka BP. Comparing the pollen-based reconstruction results with the published marine records from the NW Pacific, the reconstructed vegetation and climate dynamics can be satisfactorily explained by the greater role played by the warm Tsushima Current in the Sea of Japan and in the Tsugaru Strait during the middle and late Holocene. An increase in sea surface temperatures west and south of the study site would favour air temperature rise and moisture uptake and cause an increase in precipitation and snow accumulation in the western part of Hokkaido during the late Holocene.

Pleistocene sedimentation in the Benguela Current, southeastern Atlantic

We infer variations in paleoproductivity and eolian input at ODP Site 1082 in the Walvis Basin from stable oxygen isotope compositions of the planktonic foraminifera Globorotalia inflata, total organic carbon mass accumulation rates (TOC MAR), and X-ray fluorescence analyses of Fe content. The most pronounced paleoclimatic changes correspond to the time at about 0.9 Ma, when glacial conditions in the northern hemisphere (NH) led to the onset of pronounced 100-kyr glacial-interglacial cycles. We used Fe intensity as a proxy for eolian terrigenous input, and TOC MAR as a paleoproductivity indicator. Paleoproductivity and eolian input show generally higher-amplitude variations of glacial-interglacial cyclicity from 1.5 to 0.58 Ma, indicating pronounced variations in upwellingfavorable winds in this area. At 0.58 Ma, paleoproductivity and eolian input shifted abruptly to lower-amplitude variations with a periodicity of 100 kyr while delta18O values show a trend toward more negative isotope values for the past 0.65 Myr. Especially during glacial periods, oxygen isotope values indicate increasingly warmer sea-surface temperatures toward the end of the Pleistocene. To evaluate the relative influences of NH glaciation and southern hemisphere (SH) insolation as potential forcing mechanisms for variations of eolian input and productivity in the northern Benguela system, we filtered our proxy records at orbital frequencies. The filtered records of Fe intensity and TOC MAR indicate a strong influence of the 100-kyr and 41-kyr frequency bands, supporting our assumption that strong ice buildup in the NH is the dominant trigger for climate changes on the continent and probably in trade-wind intensity. SH insolation and low-latitude precession-related insolation changes were important for paleoproductivity variations in the northern Benguela system, modifying the nutrient supply by southern ocean intermediate waters and the zonal direction of upwelling-inducing trades by the African monsoon system, respectively.

Respiration and ingestion rates of calanoid copepod in the northern Benguela Current System measured ex situ during the RSS Discovery D356 cruise in 2010

Respiration rates of 16 calanoid copepod species from the northern Benguela upwelling system were measured on board RRS Discovery in September/October 2010 to determine their energy requirements and assess their significance in the carbon cycle. Copepod species were sampled by different net types. Immediately after the hauls, samples were sorted to species and stages (16 species; females, males and C5 copepodids) according to Bradford-Grieve et al. (1999). Specimens were kept in temperature-controlled refrigerators for at least 12 h before they were used in experiments. Respiration rates of different copepod species were measured onboard by optode respirometry (for details see Köster et al., 2008) with a 10-channel optode respirometer (PreSens Precision Sensing Oxy-10 Mini, Regensburg, Germany) under simulated in situ conditions in temperature-controlled refrigerators. Experiments were run in gas-tight glass bottles (12-13 ml). For each set of experiments, two controls without animals were measured under exactly the same conditions to compensate for potential bias. The number of animals per bottle depended on the copepods size, stage and metabolic activity. Animals were not fed during the experiments but they showed natural species-specific movements. Immediately after the experiments, all specimens were deep-frozen at - 80 °C for later dry mass determination (after lyophilisation for 48 h) in the home lab. The carbon content (% of dry mass) of each species was measured by mass-spectrometry in association with stable isotope analysis and body dry mass was converted to units of carbon. For species without available carbon data, the mean value of all copepod species (44% dry mass) was applied. For the estimation of carbon requirements of copepod species, individual oxygen consumption rates were converted to carbon units, assuming that the expiration of 1 ml oxygen mobilises 0.44 mg of organic carbon by using a respiratory quotient (RQ) of 0.82 for a mixed diet consisting of proteins (RQ = 0.8-1.0), lipids (RQ = 0.7) and carbohydrates (RQ = 1.0) (Auel and Werner, 2003). The carbon ingestion rates were calculated using the energy budget and the potential maximum ingestion rate approach. To allow for physiological comparisons of respiration rates of deep- and shallow-living copepod species without the effects of ambient temperature and different individual body mass, individual respiration rates were temperature- (15°C, Q10=2) and size-adjusted. The scaling coefficient of 0.76 (R2=0.556) is used for the standardisation of body dry mass to 0.3 mg (mean dry mass of all analysed copepods), applying the allometric equation R= (R15°C/M0.76)×0.30.76, where R is respiration and M is individual dry mass in mg.

Variability of the Brazil Current during the late Holocene analysed on two marine sediment cores

Our understanding of the centennial-scale variability of the Brazil Current (BC) during the late Holocene is elusive because of the lack of appropriate records. Here we used the Mg/Ca and oxygen isotopic composition of planktonic foraminifera from two marine sediment cores collected at 27° S and 33° S off southeastern South America to assess the late Holocene variability in the upper water column of the BC. Our results show in phase fluctuations of up to 3 °C in sea surface temperatures (SST), and 0.8 per mil in oxygen isotopic composition of surface sea water, a proxy for relative sea surface salinity (SSS). Time-series analyses of our records indicate a cyclicity with a period of ca. 730 yr. We suggest that the observed cyclicity reflects variability in the strength of the BC associated to changes in the Atlantic meridional overturning circulation (AMOC). Positive (negative) SST and SSS anomalies are related to a strong (weak) BC and a weak (strong) AMOC. Moreover, periods of peak strength in the BC occur synchronously to a weak North Brazil Current, negative SST anomalies in the high latitudes of the North Atlantic, and positive (negative) precipitation anomalies over southeastern South America (equatorial Africa), further corroborating our hypothesis. This study shows a tight coupling between the variability of the BC and the high latitudes of the North Atlantic mediated by the AMOC even under late Holocene boundary conditions.

Respiration rates and abundances of dominant copepods of the northern Benguela Current System

Respiration rates of 16 calanoid copepod species from the northern Benguela upwelling system were measured on board RRS Discovery in September/October 2010 to determine their energy requirements and assess their significance in the carbon cycle. Individual respiration rates were standardised to a mean copepod body mass and a temperature regime typical of the northern Benguela Current. These adjusted respiration rates revealed two different activity levels (active and resting) in copepodids C5 of Calanoides carinatus and females of Rhincalanus nasutus, which reduced their metabolism during dormancy by 82% and 62%, respectively. An allometric function (Imax) and an energy budget approach were performed to calculate ingestion rates. Imax generally overestimated the ingestion rates derived from the energy budget approach by >75%. We suggest that the energy budget approach is the more reliable approximation with a total calanoid copepod (mainly females) consumption of 78 mg C m-2 d-1 in neritic regions and 21 mg C m-2 d-1 in oceanic regions. The two primarily herbivorous copepods C. carinatus (neritic) and Nannocalanus minor (oceanic) contributed 83% and 5%, respectively, to total consumption by calanoid copepods. Locally, C. carinatus can remove up to 90% of the diatom biomass daily. In contrast, the maximum daily removal of dinoflagellate biomass by N. minor was 9%. These estimates imply that C. carinatus is an important primary consumers in the neritic province of the northern Benguela system, while N. minor has little grazing impact on phytoplankton populations further offshore. Data on energy requirements and total consumption rates of dominant calanoid copepods of this study are essential for the development of realistic carbon budgets and food-web models for the northern Benguela upwelling system.

Processed 2 minutes-averaged continuous VM-ADCP (vessel-mounted Acoustic Doppler Current Profiler) profiles during POLARSTERN cruise ANT-XXVIII/3

This data set was obtained during the R. V. POLARSTERN cruise ANT-XXVIII/3. Current velocities were measured nearly continuously when outside territorial waters along the ship's track with a vessel-mounted TRD Instruments' 153.6-kHz Ocean Surveyor ADCP. The transducers were located 11 m below the water line and were protected against ice floes by an acoustically transparent plastic window. The current measurements were made using a pulse of 2s and vertical bin length of 4 m. The ship's velocity was calculated from position fixes obtained by the Global Positioning System (GPS). Heading, roll and pitch data from the ship's gyro platforms and the navigation data were used to convert the ADCP velocities into earth coordinates. Accuracy of the ADCP velocities mainly depends on the quality of the position fixes and the ship's heading data. Further errors stem from a misalignment of the transducer with the ship's centerline. The ADCP data were processed using the Ocean Surveyor Sputum Interpreter (OSSI) software developed by GEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel. The averaging interval was set to 120 seconds. The reference layer was set to bins 5 to 16 avoiding near surface effects and biases near bin 1. Sampling interval setting: 2s; Number of bins: 80; Bin length: 4m; Pulse length: 4m; Blank beyond transmit length: 4m. Data processing setting: Top reference bin: 5; Bottom reference bin: 16; Average: 120s; Misalignment amplitude: 1.0276 +/- 0.1611, phase: 0.8100 +/- 0.7190. The precision for single ping and 4m cell size reported by TRDI is 0.30m/s. Resulting from the single ping precision and the number of pings (most of the time 36) during 120seconds the velocity accuracy is nearly 0.05m/s. (Velocity accuracy = single ping precision divided by square root of the number of pings).