U-Pb (zircon) geochronologic analyses for rocks of the Transantarctic Mountains (Table 4.3-1)

The Byrd Glacier discontinuity us a major boundary crossing the Ross Orogen, with crystalline rocks to the north and primarily sedimentary rocks to the south. Most models for the tectonic development of the Ross Orogen in the central Transantarctic Mountains consits of two-dimensional transects across the belt, but do not adress the major longitudinal contrast at Byrd Glacier. This paper presents a tectonic model centering on the Byrd Glacier discontinuity. Rifting in the Neoproterozoic producede a crustal promontory in the craton margin to the north of Byrd Glacier. Oblique convergence of the terrane (Beardmore microcontinent) during the latest Neroproterozoic and Early Cambrian was accompanied by subduction along the craton margin of East Antarctica. New data presented herein in the support of this hypothesis are U-Pb dates of 545.7 ± 6.8 Ma and 531.0 ± 7.5 Ma on plutonic rocks from the Britannia Range, subduction stepped out, and Byrd Glacier. After docking of the terrane, subduction stepped out, and Byrd Group was deposited during the Atdabanian-Botomian across the inner margin of the terrane. Beginning in the upper Botomian, reactivation of the sutured boundaries of the terrane resulted in an outpouring of clastic sediment and folding and faulting of the Byrd Group.

Tab. 1: U-Pb analyses of zircons from rock samples in MT. Newton and Cumpston Massif, East Antarctica

We present new U-Pb zircon (SHRIMP) data on rocks from Mt Newton and Cumpston Massif in the southern Prince Charles Mountains. Our data demonstrate that Mt Newton was affected by a newly proposed Palaeoproterozoic "Newton" Orogeny at c. 2100-2200 Ma. Sedimentation, felsic volcanism (c. 2200 Ma), metamorphism and folding, followed by granite intrusion (c. 2100 Ma), suggest development of a trough or aulacogene in the area during the early Palaeoproterozoic. An orthogneiss from Cumpston Massif yielded an age of c. 3180 Ma for granitic protolith emplacement, which is in good agreement with many U-Pb zircon ages from similar rocks in the southern Mawson Escarpment. A syn- to late-tectonic muscovite-bearing pegmatite from Cumpston Massif yielded a c. 2500 Ma date of emplacement, which indicates early Palaeoproterozoic activity in this block, probably in response to a tectono-magmatic episode in the Lambert Terrane bordering the Ruker Terrane in the northeast. The correlation of tectono-magmatic events in both the Ruker and Lambert terranes of the southern Prince Charles Mountains provides evidence for their common evolution during the Proterozoic.

Shipborne ocean color remote sensing

In this paper, modernized shipborne procedures are presented to collect and process above-water radiometry for remote sensing applications. A setup of five radiometers and a bidirectional camera system, which provides panoramic sea surface and sky images, is proposed for the collection of high-resolution radiometric quantities. Images from the camera system can be used to determine sky state and potential glint, whitecaps, or foam contamination. A peak in the observed remote sensing reflectance RRS spectra between 750-780 nm was typically found in spectra with relatively high surface reflected glint (SRG), which suggests this waveband could be a useful SRG indicator. Simplified steps for computing uncertainties in SRG corrected RRS are proposed and discussed. The potential of utilizing "unweighted multimodel averaging," which is the average of four or more common SRG correction models, is examined to determine the best approximation RRS. This best approximation RRS provides an estimate of RRS based on various SRG correction models established using radiative transfer simulations and field investigations. Applying the average RRS provides a measure of the inherent uncertainties or biases that result from a user subjectively choosing any one SRG correction model. Comparisons between inherent and apparent optical property derived observations were used to assess the robustness of the SRG multimodel averaging ap- proach. Correlations among the standard SRG models were completed to determine the degree of association or similarities between the SRG models. Results suggest that the choice of glint models strongly affects derived RRS values and can also influence the blue to green band ratios used for modeling biogeochemical parameters such as for chlorophyll a. The objective here is to present a uniform and traceable methodology for determining ship- borne RRS measurements and its associated errors due to glint correction and to ensure the direct comparability of these measurements in future investigations. We encourage the ocean color community to publish radiometric field measurements with matching and complete metadata in open access repositories.

Late Glacial and Holocene data from sediment cores MD99-2236 and MD99-2256

The Laurentide Ice Sheet (LIS) was a large, dynamic ice sheet in the early Holocene. The glacial events through Hudson Strait leading to its eventual demise are recorded in the well-dated Labrador shelf core, MD99-2236 from the Cartwright Saddle. We develop a detailed history of the timing of ice-sheet discharge events from the Hudson Strait outlet of the LIS during the Holocene using high-resolution detrital carbonate, ice rafted detritus (IRD), d18O, and sediment color data. Eight detrital carbonate peaks (DCPs) associated with IRD peaks and light oxygen isotope events punctuate the MD99-2236 record between 11.5 and 8.0 ka. We use the stratigraphy of the DCPs developed from MD99-2236 to select the appropriate DeltaR to calibrate the ages of recorded glacial events in Hudson Bay and Hudson Strait such that they match the DCPs in MD99-2236. We associate the eight DCPs with H0, Gold Cove advance, Noble Inlet advance, initial retreat of the Hudson Strait ice stream (HSIS) from Hudson Strait, opening of the Tyrrell Sea, and drainage of glacial lakes Agassiz and Ojibway. The opening of Foxe Channel and retreat of glacial ice from Foxe Basin are represented by a shoulder in the carbonate data. DeltaR of 350 years applied to the radiocarbon ages constraining glacial events H0 through the opening of the Tyrell Sea provided the best match with the MD99-2236 DCPs; DeltaR values and ages from the literature are used for the younger events. A very close age match was achieved between the 8.2 ka cold event in the Greenland ice cores, DCP7 (8.15 ka BP), and the drainage of glacial lakes Agassiz and Ojibway. Our stratigraphic comparison between the DCPs in MD99-2236 and the calibrated ages of Hudson Strait/Bay deglacial events shows that the retreat of the HSIS, the opening of the Tyrell Sea, and the catastrophic drainage of glacial lakes Agassiz and Ojibway at 8.2 ka are separate events that have been combined in previous estimates of the timing of the 8.2 ka event from marine records. SW Iceland shelf core MD99-2256 documents freshwater entrainment into the subpolar gyre from the Hudson Strait outlet via the Labrador, North Atlantic, and Irminger currents. The timing of freshwater release from the LIS Hudson Strait outlet in MD99-2236 matches evidence for freshwater forcing and LIS icebergs carrying foreign minerals to the SW Iceland shelf between 11.5 and 8.2 ka. The congruency of these records supports the conclusion of the entrainment of freshwater from the retreat of the LIS through Hudson Strait into the subpolar gyre and provides specific time periods when pulses of LIS freshwater were present to influence climate.

Age models for sediment cores off Peru

Planktonic foraminiferal faunas of the southeast Pacific indicate that sea surface temperatures (SST) have varied by as much as 8-10°C in the Peru Current, and by ?5-7°C along the equator, over the past 150,000 years. Changes in SST at times such as the Last Glacial Maximum reflect incursion of high-latitude species Globorotalia inflata and Neogloboquadrina pachyderma into the eastern boundary current and as far north as the equator. A simple heat budget model of the equatorial Pacific shows that observed changes in Peru Current advection can account for about half of the total variability in equatorial SSTs. The remaining changes in equatorial SST, which are likely related to local changes in upwelling or pycnocline depth, precede changes in polar climates as recorded by d18O. This partitioning of processes in eastern equatorial Pacific SST reveals that net ice-age cooling here reflects first a rapid response of equatorial upwelling to insolation, followed by a later response to changes in the eastern boundary current associated with high-latitude climate (which closely resembles variations in atmospheric CO2 as recorded in the Vostok ice core). Although precise mechanisms responsible for the equatorial upwelling component of climate change remain uncertain, one likely candidate that may operate independently of the ice sheets is insolation-driven changes in El Niño/Southern Oscillation (ENSO) frequency. Early responses of equatorial SST detected both here and elsewhere highlight the sensitivity of tropical systems to small changes in seasonal insolation. The scale of tropical changes we have observed are substantially greater than model predictions, suggesting a need for further quantitative assessment of processes associated with long-term climate change.

Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Turbulent flux and meteorological measurements during ARCTEX-2006 campaign

This dataset present result from the DFG- funded Arctic-Turbulence-Experiment (ARCTEX-2006) performed by the University of Bayreuth on the island of Svalbard, Norway, during the winter/spring transition 2006. From May 5 to May 19, 2006 turbulent flux and meteorological measurements were performed on the monitoring field near Ny-Ålesund, at 78°55'24'' N, 11°55'15'' E Kongsfjord, Svalbard (Spitsbergen), Norway. The ARCTEX-2006 campaign site was located about 200 m southeast of the settlement on flat snow covered tundra, 11 m to 14 m above sea level. The permanent sites used for this study consisted of the 10 m meteorological tower of the Alfred Wegener Institute for Polar- and Marine Research (AWI), the international standardized radiation measurement site of the Baseline Surface Radiation Network (BSRN), the radiosonde launch site and the AWI tethered balloon launch sites. The temporary sites - set up by the University of Bayreuth - were a 6 m meteorological gradient tower, an eddy-flux measurement complex (EF), and a laser-scintillometer section (SLS). A quality assessment and data correction was applied to detect and eliminate specific measurement errors common at a high arctic landscape. In addition, the quality checked sensible heat flux measurements are compared with bulk aerodynamic formulas that are widely used in atmosphere-ocean/land-ice models for polar regions as described in Ebert and Curry (1993, doi:10.1029/93JC00656) and Launiainen and Cheng (1995). These parameterization approaches easily allow estimation of the turbulent surface fluxes from routine meteorological measurements. The data show: - the role of the intermittency of the turbulent atmospheric fluctuation of momentum and scalars, - the existence of a disturbed vertical temperature profile (sharp inversion layer) close to the surface, - the relevance of possible free convection events for the snow or ice melt in the Arctic spring at Svalbard, and - the relevance of meso-scale atmospheric circulation pattern and air-mass advection for the near-surface turbulent heat exchange in the Arctic spring at Svalbard. Recommendations and improvements regarding the interpretation of eddy-flux and laser-scintillometer data as well as the arrangement of the instrumentation under polar distinct exchange conditions and (extreme) weather situations could be derived.

Elevation Models of Antarctica and Greenland derived from CryoSat-2 in the period 2011 to 2013

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be -503 ± 107 km**3/yr. Greenland contributes nearly 75% to the total volume change with -375 ± 24 km**3/yr.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Age models for sediment cores in the eastern tropical Pacific

We present new high-resolution N isotope records from the Gulf of Tehuantepec and the Nicaragua Basin spanning the last 50-70 ka. The Tehuantepec site is situated within the core of the north subtropical denitrification zone while the Nicaragua site is at the southern boundary. The d15N record from Nicaragua shows an 'Antarctic' timing similar to denitrification changes observed off Peru-Chile but is radically different from the northern records. We attribute this to the leakage of isotopically heavy nitrate from the South Pacific oxygen minimum zone (OMZ) into the Nicaragua Basin. The Nicaragua record leads the other eastern tropical North Pacific (ETNP) records by about 1000 years because denitrification peaks in the eastern tropical South Pacific (ETSP) before denitrification starts to increase in the Northern Hemisphere OMZ, i.e., during warming episodes in Antarctica. We find that the influence of the heavy nitrate leakage from the ETSP is still noticeable, although attenuated, in the Gulf of Tehuantepec record, particularly at the end of the Heinrich events, and tends to alter the recording of millennial timescale denitrification changes in the ETNP. This implies (1) that sedimentary d15N records from the southern parts of the ETNP cannot be used straightforwardly as a proxy for local denitrification and (2) that denitrification history in the ETNP, like in the Arabian Sea, is synchronous with Greenland temperature changes. These observations reinforce the conclusion that on millennial timescales during the last ice age, denitrification in the ETNP is strongly influenced by climatic variations that originated in the high-latitude North Atlantic region, while commensurate changes in Southern Ocean hydrography more directly, and slightly earlier, affected oxygen concentrations in the ETSP. Furthermore, the d15N records imply ongoing physical communication across the equator in the shallow subsurface continuously over the last 50-70 ka.

Age models, organic analyses and major-element analyses of sediment cores on a transect of the tropical African rainbelt

The distribution of rainfall in tropical Africa is controlled by the African rainbelt**1, which oscillates on a seasonal basis. The rainbelt has varied on centennial to millennial timescales along with changes in Northern Hemisphere high-latitude climate**2, 3, 4, 5, the Atlantic meridional overturning circulation**6 and low-latitude insolation**7 over the past glacial-interglacial cycle. However, the overall dynamics of the African rainbelt remain poorly constrained and are not always consistent with a latitudinal migration**2, 4, 5, 6, as has been proposed for other regions**8, 9. Here we use terrestrially derived organic and sedimentary markers from marine sediment cores to reconstruct the distribution of vegetation, and hence rainfall, in tropical Africa during extreme climate states over the past 23,000 years. Our data indicate that rather than migrating latitudinally, the rainbelt contracted and expanded symmetrically in both hemispheres in response to changes in climate. During the Last Glacial Maximum and Heinrich Stadial 1, the rainbelt contracted relative to the late Holocene, which we attribute to a latitudinal compression of atmospheric circulation associated with lower global mean temperatures**10. Conversely, during the mid-Holocene climatic optimum, the rainbelt expanded across tropical Africa. In light of our findings, it is not clear whether the tropical rainbelt has migrated latitudinally on a global scale, as has been suggested**8,9.

Late Pliocene age models of Agulhas Plateau sediment cores

An ensemble of new, high-resolution records of surface ocean hydrography from the Indian-Atlantic oceanic gateway, south of Africa, demonstrates recurrent and high-amplitude salinity oscillations in the Agulhas Leakage area during the penultimate glacial-interglacial cycle. A series of millennial-scale salinification events, indicating strengthened salt leakage into the South Atlantic, appear to correlate with abrupt changes in the North Atlantic climate and Atlantic Meridional Overturning Circulation (AMOC). This interhemispheric coupling, which plausibly involved changes in the Hadley Cell and midlatitude westerlies that impacted the interocean transport at the tip of Africa, suggests that the Agulhas Leakage acted as a source of negative buoyancy for the perturbed AMOC, possibly aiding its return to full strength. Our finding points to the Indian-to-Atlantic salt transport as a potentially important modulator of the AMOC during the abrupt climate changes of the Late Pleistocene.