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.

Geochemistry of lavas from the Lau-Tonga arc and back-arc system

Detailed comparison of mineralogy, and major and trace geochemistry are presented for the modern Lau Basin spreading centers, the Sites 834-839 lavas, the modern Tonga-Kermadec arc volcanics, the northern Tongan boninites, and the Lau Ridge volcanics. The data clearly confirm the variations from near normal mid-ocean-ridge basalt (N-MORB) chemistries (e.g., Site 834, Central Lau Spreading Center) to strongly arc-like (e.g., Site 839, Valu Fa), the latter closely comparable to the modern arc volcanoes. Sites 835 and 836 and the East Lau Spreading Center represent transitional chemistries. Bulk compositions range from andesitic to basaltic, but lavas from Sites 834 and 836 and the Central Lau Spreading Center extend toward more silica-undersaturated compositions. The Valu Fa and modern Tonga-Kermadec arc lavas, in contrast, are dominated by basaltic andesites. The phenocryst and groundmass mineralogies show the strong arc-like affinities of the Site 839 lavas, which are also characterized by the existence of very magnesian olivines (up to Fo90-92) and Cr-rich spinels in Units 3 and 6, and highly anorthitic plagioclases in Units 2 and 9. The regional patterns of mineralogical and geochemical variations are interpreted in terms of two competing processes affecting the inferred magma sources: (1) mantle depletion processes, caused by previous melt extractions linked to backarc magmatism, and (2) enrichment in large-ion-lithophile elements, caused by a subduction contribution. A general trend of increasing depletion is inferred both eastward across the Lau Basin toward the modern arc, and northward along the Tongan (and Kermadec) Arc. Numerical modeling suggests that multistage magma extraction can explain the low abundances (relative to N-MORB) of elements such as Nb, Ta, and Ti, known to be characteristic of island arc magmas. It is further suggested that a subduction jump following prolonged slab rollback could account for the initiation of the Lau Basin opening, plausibly allowing a later influx of new mantle, as required by the recognition of a two-stage opening of the Lau Basin.

Age determination of sediment cores of the California Current System

Alkenone sea surface temperature (SST) records were generated from the Ocean Drilling Program's (ODP) Sites 1014 and 1016 to examine the response of the California Current System to global climate change during the last 136 ka. The temperature differences between these sites (Delta SST(NEP)=SST(ODP1014)-SST(ODP1016)) reflected the intensity of the California Current and varied between 0.4 and 6.1 °C. A high Delta SST(NEP) (weaker California Current) was found for late marine isotope stage (MIS) 2 and early MIS 5e, while a low Delta SST(NEP) (stronger California Current) was detected for mid-MIS 5e and MIS 1. Spectral analysis indicated that this variation pattern dominated 23- (precession) and 30-ka periods. Comparison of the Delta SST(NEP) and SST based on data from core MD01-2421 at the Japan margin revealed anti-phase variation; the high Delta SST(NEP) (weakening of the California Current) corresponded to the low SST at the Japan margin (the southward displacement of the NW Pacific subarctic boundary), and vice versa. This variation was synchronous with a model prediction of the tropical El Niño-Southern Oscillation behavior. These findings suggest that the intensity of the North Pacific High varied in response to precessional forcing, and also that the response has been linked with the changes of tropical ocean-atmosphere interactions.

(Table 2) Comparison of annual mass fluxes measured with shallower and deeper sediment traps off Cape Blanc (CB)

The ocean off NW Africa is the second most important coastal upwelling system with a total annual primary production of 0.33 Gt of carbon per year (Carr in Deep Sea Res II 49:59-80, 2002). Deep ocean organic carbon fluxes measured by sediment traps are also fairly high despite low biogenic opal fluxes. Due to a low supply of dissolved silicate from subsurface waters, the ocean off NW Africa is characterized by predominantly carbonate-secreting primary producers, i.e. coccolithophorids. These algae which are key primary producers since millions of years are found in organic- and chlorophyll-rich zooplankton fecal pellets, which sink rapidly through the water column within a few days. Particle flux studies in the Mauretanian upwelling area (Cape Blanc) confirm the hypothesis of Armstrong et al. (Deep Sea Res II 49:219-236, 2002) who proposed that ballast availability, e.g. of carbonate particles, is essential to predict deep ocean organic carbon fluxes. The role of dust as ballast mineral for organic carbon, however, must be also taken into consideration in the coastal settings off NW Africa. There, high settling rates of larger particles approach 400 m day**-1, which may be due to a particular composition of mineral ballast. An assessment of particle settling rates from opal-production systems in the Southern Ocean of the Atlantic Sector, in contrast, provides lower values, consistent with the assumptions of Francois et al. (Global Biogeochem Cycles 16(4):1087, 2002). Satellite chlorophyll distributions, particle distributions and fluxes in the water column off NW Africa as well as modelling studies suggest a significant lateral flux component and export of particles from coastal shelf waters into the open ocean. These transport processes have implications for paleo-reconstructions from sediment cores retrieved at continental margin settings.

High-resolution sediment record from a submarine meandering canyon system in the eastern Atlantic

Based on a high-resolution sediment record from a submarine meandering canyon system offshore the present-day hyperarid Saharan Africa, two phases of turbidity-current activity can be distinguished during the past 13,000 years. Frequent, siliciclastic turbidity currents can be related to deglacial sea-level history, whereas rhythmically recurring fine-grained and carbonate-rich turbidity currents with recurrence times of roughly 900 years are inferred for the Holocene. Various trigger mechanisms can be considered to initiate turbidity currents, but only a few can explain a periodic turbidite activity. A comparison of Holocene turbidite recurrence times and basic cycles of 900 and 1,800 years found in various Holocene paleoclimate studies suggests that a previously unrecognized climate-related coupling may be active.

(Tables 1-3) Pb and Sr isotopic data for sediments and rocks from DSDP holes in the Mariana island-arc system

Analyses of the isotopic composition of Pb in (1) western Pacific Ocean sediments [Jurassic(?) to Pleistocene in age, including clays and biogenic oozes], (2) Pacific Ocean basaltic rocks, (3) Mariana frontal arc volcanic rocks (Eocene to Miocene), and (4) Mariana active arc volcanic rocks [Pliocene (?) to Holocene] indicate that Pacific Ocean sediments could not have been a significant component of the source material for the Mariana arc volcanic rocks. Calculations involving the average concentrations and isotopic compositions of Pb in oceanic sediments, sea-floor basaltic rocks, and the Mariana arc volcanic rocks suggest that the sediment component must have been less than 1 percent of this source material. The Pb isotopic compositions of the Mariana arc volcanic rocks lie, within experimental error, along the trend of available Pacific Ocean basalt analyses in versus 207Pb/204Pb versus 206Pb/204Pb and 208Pb/204Pb versus 206Pb/204Pb diagrams. Isotopic analyses of Pb in Pacific Ocean sediments do not lie along this trend; they have higher 207Pb/204Pb and 208Pb/204Pb values for comparable 206Pb/204Pb ratios. Clayey sediments generally have higher 208Pb/204Pb and 207Pb/204Pb ratios than biogenic oozes regardless of the age of the sediment. Comparison of combined Sr and Pb isotopic analyses for (1) mantle-derived materials erupted through oceanic crust, (2) altered ocean-floor basaltic rocks, and (3) volcanic rocks from oceanic island arcs suggests that the Mariana arc volcanic rocks were derived, at least in part, from altered Pacific lithosphere subducted beneath the Mariana arc. Unaltered basalts from the Mariana inter-arc basin (Mariana Trough) have Pb and Sr isotopic compositions that are very similar to those reported for some Hawaiian volcanic rocks but distinct from Mariana active and frontal arc compositions. These observations, in addition to existing major-and trace-element data, support a mantle origin for the interarc basin volcanic rocks. Dacites dredged from the Mariana remnant arc (South Honshu Ridge) have Pb isotopic compositions that are within experimental error of the active-arc analyses, consistent with a genetic relation.

Tropical climate and vegetation simulations during the Heinrich event 1 using an Earth System Model of Intermediate Complexity (EMIC) - the University of Victoria Earth System-Climate Model (UVic ESCM)

Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.

Mineralogical compositions of sediment samples from the San Francisco Bay coastal system

The mineralogical compositions of 119 samples collected from throughout the San Francisco Bay coastal system, including bayfloor and seafloor, area beaches, cliff outcrops, and major drainages, were determined using X-ray diffraction (XRD). Comparison of the mineral concentrations and application of statistical cluster analysis of XRD spectra allowed for the determination of provenances and transport pathways. The use of XRD mineral identifications provides semi-quantitative compositions needed for comparisons of beach and offshore sands with potential cliff and river sources, but the innovative cluster analysis of XRD diffraction spectra provides a unique visualization of how groups of samples within the San Francisco Bay coastal system are related so that sand-sized sediment transport pathways can be inferred. The main vector for sediment transport as defined by the XRD analysis is from San Francisco Bay to the outer coast, where the sand then accumulates on the ebb tidal delta and also moves alongshore. This mineralogical link defines a critical pathway because large volumes of sediment have been removed from the Bay over the last century via channel dredging, aggregate mining, and borrow pit mining, with comparable volumes of erosion from the ebb tidal delta over the same period, in addition to high rates of shoreline retreat along the adjacent, open-coast beaches. Therefore, while previously only a temporal relationship was established, the transport pathway defined by mineralogical and geochemical tracers support the link between anthropogenic activities in the Bay and widespread erosion outside the Bay. The XRD results also establish the regional and local importance of sediment derived from cliff erosion, as well as both proximal and distal fluvial sources. This research is an important contribution to a broader provenance study aimed at identifying the driving forces for widespread geomorphic change in a heavily urbanized coastal-estuarine system.

Carbonate system during the May 2013 MedSeA cruise

Most global ocean models are based on the assumption of a "steady state" ocean. Here, we investigate the validation of this hypothesis for the anthropized Mediterranean Sea. In order to do so, we calculated the mixing coefficients of the water masses detected in this sea via an optimum multiparameter analysis referred to as the MIX approach, using data from the BOUM (2008) and MedSeA (2013) cruises. The comparison of the mixing coefficients of each water mass, between 2008 and 2013, indicates that some of their proportions have significantly changed. Surface water mass proportions did not change significantly (Delta0.05-0.1), while intermediate and deep water mass mixing coefficients of both Eastern and Western basins were significantly modified (~Delata 0.35). This study clearly shows that the Mediterranean seawater is not in a "steady state".

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.