Sr-Nd isotopic analyses of sand-sized sediment from the San Francisco Bay coastal system

A diverse suite of geochemical tracers, including 87Sr/86Sr and 143Nd/144Nd isotope ratios, the rare earth elements (REEs), and select trace elements were used to determine sand-sized sediment provenance and transport pathways within the San Francisco Bay coastal system. This study complements a large interdisciplinary effort (Barnard et al., 2012) that seeks to better understand recent geomorphic change in a highly urbanized and dynamic estuarine-coastal setting. Sand-sized sediment provenance in this geologically complex system is important to estuarine resource managers and was assessed by examining the geographic distribution of this suite of geochemical tracers from the primary sources (fluvial and rock) throughout the bay, adjacent coast, and beaches. Due to their intrinsic geochemical nature, 143Nd/144Nd isotopic ratios provide the most resolved picture of where sediment in this system is likely sourced and how it moves through this estuarine system into the Pacific Ocean. For example, Nd isotopes confirm that the predominant source of sand-sized sediment to Suisun Bay, San Pablo Bay, and Central Bay is the Sierra Nevada Batholith via the Sacramento River, with lesser contributions from the Napa and San Joaquin Rivers. Isotopic ratios also reveal hot-spots of local sediment accumulation, such as the basalt and chert deposits around the Golden Gate Bridge and the high magnetite deposits of Ocean Beach. Sand-sized sediment that exits San Francisco Bay accumulates on the ebb-tidal delta and is in part conveyed southward by long-shore currents. Broadly, the geochemical tracers reveal a complex story of multiple sediment sources, dynamic intra-bay sediment mixing and reworking, and eventual dilution and transport by energetic marine processes. Combined geochemical results provide information on sediment movement into and through San Francisco Bay and further our understanding of how sustained anthropogenic activities which limit sediment inputs to the system (e.g., dike and dam construction) as well as those which directly remove sediments from within the Bay, such as aggregate mining and dredging, can have long-lasting effects.

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.

Mineralogy of sediment core GeoB3313-1

Geochemical and clay mineral parameters of a high accumulation marine sediment core from the Chilean continental slope (41°S) provide a 7700 yr record of rainfall variability in southern Chile related to the position of the Southern Westerlies. We especially use the iron content, measured with a time-resolution of ca. 10 yr on average, of 14C-accelerator mass spectrometry dated marine sediments as a proxy for the relative input of iron-poor Coastal Range and iron-rich Andean source rocks. Variations in this input are most likely induced by rainfall changes in the continental hinterland of the core position. Based on these interpretations, we find a pronounced rainfall variability on multi-centennial to millennial time-scales, superimposed on generally more arid conditions during the middle Holocene (7700 to 4000 cal yr B.P.) compared to the late Holocene (4000 to present). This variability and thus changes in the position of the Southern Westerlies are first compared to regional terrestrial paleoclimate data-sets from central and southern Chile. In order to derive possible wider implications and forcing mechanisms of the Holocene latitudinal shifts of the Southern Westerlies, we then compare our data to ice-core records from both tropical South America and coastal Antarctica. These records show similar bands of variability centered at ca. 900 and 1500 yr. Comparisons of band pass filters suggest a close connection of shifts of the Southern Westerlies to changes within the tropical climate system. The correlation to climate conditions in coastal Antarctica shows a more complicated picture with a phase shift at the beginning of the late Holocene coinciding with the onset of the modern state of El Niño-Southern Oscillation system. The presented data provide further evidence that the well known millennial-scale climate variability during the last glacial continued throughout the Holocene.

Plio-Pleistocene Mg/Ca and TEX86 sea surface temperature records of ODP Sites 184-1143 and 165-999,

The western warm pools of the Atlantic and Pacific oceans are a critical source of heat and moisture for the tropical climate system. Over the past five million years, global mean temperatures have cooled by 3-4 °C. Yet, current reconstructions of sea surface temperatures indicate that temperature in the warm pools has remained stable during this time. This stability has been used to suggest that tropical sea-surface temperatures are controlled by some sort of thermostat-like regulation. Here we reconstruct sea surface temperatures in the South China Sea, Caribbean Sea and western equatorial Pacific Ocean for the past five million years, using a combination of the Mg/Ca, TEXH86-and Uk'37 surface temperature proxies. Our data indicate that during the period of Pliocene warmth from about 5 to 2.6 million years ago, the western Pacific and western Atlantic warm pools were about 2 °C warmer than today. We suggest that the apparent lack of warming seen in the previous reconstructions was an artefact of low seawater Mg/Ca ratios in the Pliocene oceans. Taking this bias into account, our data indicate that tropical sea surface temperatures did change in conjunction with global mean temperatures. We therefore conclude that the temperature of the warm pools of the equatorial oceans during the Pliocene was not limited by a thermostat-like mechanism.

Hydrology dataset based on 111 CTD stations in the North Aegean Sea during June 1998 from the project INTERREGAEGEAN

The combination of two research projects offered us the opportunity to perform a comprehensive study of the seasonal evolution of the hydrological structure and the circulation of the North Aegean Sea, at the northern extremes of the eastern Mediterranean. The combination of brackish water inflow from the Dardanelles and the sea-bottom relief dictate the significant differences between the North and South Aegean water columns. The relatively warm and highly saline South Aegean waters enter the North Aegean through the dominant cyclonic circulation of the basin. In the North Aegean, three layers of distinct water masses of very different properties are observed: The 20-50 m thick surface layer is occupied mainly by Black Sea Water, modified on its way through the Bosphorus, the Sea of Marmara and the Dardanelles. Below the surface layer there is warm and highly saline water originating in the South Aegean and the Levantine, extending down to 350-400 m depth. Below this layer, the deeper-than-400 m basins of the North Aegean contain locally formed, very dense water with different i/S characteristics at each subbasin. The circulation is characterised by a series of permanent, semi-permanent and transient mesoscale features, overlaid on the general slow cyclonic circulation of the Aegean. The mesoscale activity, while not necessarily important in enhancing isopycnal mixing in the region, in combination with the very high stratification of the upper layers, however, increases the residence time of the water of the upper layers in the general area of the North Aegean. As a result, water having out-flowed from the Black Sea in the winter, forms a separate distinct layer in the region in spring (lying between "younger" BSW and the Levantine origin water), and is still traceable in the water column in late summer.

(Table 1) Age determination of sediment core SO139-74KL

A multiproxy record has been acquired from a piston core (SO139-74KL) taken offshore southern Sumatra, an area which is situated in the southwestern sector of the tropical Indo-Pacific Warm Pool. The high-resolution data sets (X-ray fluorescence, total organic carbon, and C37 alkenones) were used to track changes in paleoproductivity, freshwater budget, and sea surface temperature (SST) of the tropical climate system at orbital time scales over the past 300 ka. Our paleoclimatic data show that enhanced marine paleoproductivity was directly related to strengthening of coastal upwelling during periods of increased boreal summer insolation and associated SE monsoon strength with a precessional cyclicity. Changes in freshwater supply were primarily forced by precession-controlled changes in boreal NW winter monsoon rainfall enclosing an additional sea level component. SST variations of 2°-5°C occurred at eccentricity and precessional cyclicity. We suggest that the sea surface temperature variability off southern Sumatra is predominantly related to three major causes: (1) variations in upwelling intensity; (2) an elevated freshwater input into the southern Makassar Strait leading to reduced supply of warmer surface waters from the western Pacific and increased subsurface water transport via the Indonesian Throughflow into the Indian Ocean; and (3) long-term changes in the intensity or frequency of low-latitude climate phenomena, such as El Niño-Southern Oscillation.

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.

Bedform measurements from the San Francisco Bay Coastal System

The morphology of ~45,000 bedforms from 13 multibeam bathymetry surveys was used as a proxy for identifying net bedload sediment transport directions and pathways throughout the San Francisco Bay estuary and adjacent outer coast. The spatially-averaged shape asymmetry of the bedforms reveals distinct pathways of ebb and flood transport. Additionally, the region-wide, ebb-oriented asymmetry of 5% suggests net seaward-directed transport within the estuarine-coastal system, with significant seaward asymmetry at the mouth of San Francisco Bay (11%), through the northern reaches of the Bay (7-8%), and among the largest bedforms (21% for lambda > 50 m). This general indication for the net transport of sand to the open coast strongly suggests that anthropogenic removal of sediment from the estuary, particularly along clearly defined seaward transport pathways, will limit the supply of sand to chronically eroding, open-coast beaches. The bedform asymmetry measurements significantly agree (up to ~ 76%) with modeled annual residual transport directions derived from a hydrodynamically-calibrated numerical model, and the orientation of adjacent, flow-sculpted seafloor features such as mega-flute structures, providing a comprehensive validation of the technique. The methods described in this paper to determine well-defined, cross-validated sediment transport pathways can be applied to estuarine-coastal systems globally where bedforms are present. The results can inform and improve regional sediment management practices to more efficiently utilize often limited sediment resources and mitigate current and future sediment supply-related impacts.

Benthic foraminiferal stable isotope stratigraphy of ODP Site 138-846 in the tropical Pacific Ocean

A stable-isotope stratigraphy at Site 846 (tropical Pacific, 3°06'S, 90°49'W, 3307 m water depth), based on the benthic foraminifers Cibicides wuellerstorfi and Uvigerina peregrina, yields a high-resolution record of deep-sea delta18O and delta13C over the past 1.8 Ma, with an average sampling interval of 3 k.y. Variance in the delta18O and delta13C records is concentrated in the well-known orbital periods of 100, 41, and 23 k.y. In the 100-k.y. band, both isotopic signals grow from relatively low amplitudes prior to 1.2 Ma, to high amplitudes in the late Quaternary since 0.7 Ma. The amplitude of delta18O and especially of delta13C decreases in the 41-k.y. band as it grows in the 100-k.y. band, consistent with a transfer of energy into an orbitally-paced internal oscillation. A weak 30-k.y. rhythm, present in both delta18O and delta13C, may reflect nonlinear interaction between the 41-k.y. and 100-k.y. bands in the evolving climate system. In the 23-k.y. and 19-k.y. bands associated with orbital precession, delta18O and delta13C are not coherent with each other on long time scales, and do not evolve like the 100-k.y. and 41-k.y. bands. This suggests that the source of the growing 100-k.y. oscillation is not a nonlinear response to precession, in contrast to predictions of some climate models. Sedimentation rates at this site also vary with a strong 100-k.y. cycle. Unlike the isotope records, the amplitude of 100-k.y. variations in sedimentation rate is relatively constant over the past 1.8 Ma, ranging from about 15 to 70 m/m.y. Prior to 0.9 Ma, sedimentation rates co-vary with orbital eccentricity, rather than with global climate as reflected by delta18O or delta13C. A source of this 100-k.y. cycle of sedimentation rate in the absence of similar ice volume fluctuations may be precessional heating of equatorial land masses, which in an energy balance climate model drives variations of monsoonal climates with a 100-k.y. rhythm. For the interval younger than 0.9 Ma, high sedimentation rates in the 100-k.y. band are consistently associated with glacial stages. This change of pattern suggests that when the amplitude of glacial cycles become large enough, their global effects overpower a local monsoon-driven variation in sedimentation rate at Site 846.

Images of a groundwater flow system in Nicaragua

Conceptualization of groundwater flow systems is necessary for water resources planning. Geophysical, hydrochemical and isotopic characterization methods were used to investigate the groundwater flow system of a multi-layer fractured sedimentary aquifer along the coastline in Southwestern Nicaragua. A geologic survey was performed along the 46 km2 catchment. Electrical resistivity tomography (ERT) was applied along a 4.4 km transect parallel to the main river channel to identify fractures and determine aquifer geometry. Additionally, three cross sections in the lower catchment and two in hillslopes of the upper part of the catchment were surveyed using ERT. Stable water isotopes, chloride and silica were analyzed for springs, river, wells and piezometers samples during the dry and wet season of 2012. Indication of moisture recycling was found although the identification of the source areas needs further investigation. The upper-middle catchment area is formed by fractured shale/limestone on top of compact sandstone. The lower catchment area is comprised of an alluvial unit of about 15 m thickness overlaying a fractured shale unit. Two major groundwater flow systems were identified: one deep in the shale unit, recharged in the upper-middle catchment area; and one shallow, flowing in the alluvium unit and recharged locally in the lower catchment area. Recharged precipitation displaces older groundwater along the catchment, in a piston flow mechanism. Geophysical methods in combination with hydrochemical and isotopic tracers provide information over different scales and resolutions, which allow an integrated analysis of groundwater flow systems. This approach provides integrated surface and subsurface information where remoteness, accessibility, and costs prohibit installation of groundwater monitoring networks.

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.

Calcareous green alga Halimeda tolerates ocean acidification conditions at tropical carbon dioxide seeps

We investigated ecological, physiological, and skeletal characteristics of the calcifying green alga Halimeda grown at CO2 seeps (pHtotal ? 7.8) and compared them to those at control reefs with ambient CO2 conditions (pHtotal ? 8.1). Six species of Halimeda were recorded at both the high CO2 and control sites. For the two most abundant species Halimeda digitata and Halimeda opuntia we determined in situ light and dark oxygen fluxes and calcification rates, carbon contents and stable isotope signatures. In both species, rates of calcification in the light increased at the high CO2 site compared to controls (131% and 41%, respectively). In the dark, calcification was not affected by elevated CO2 in H. digitata, whereas it was reduced by 167% in H. opuntia, suggesting nocturnal decalcification. Calculated net calcification of both species was similar between seep and control sites, i.e., the observed increased calcification in light compensated for reduced dark calcification. However, inorganic carbon content increased (22%) in H. digitata and decreased (-8%) in H. opuntia at the seep site compared to controls. Significantly, lighter carbon isotope signatures of H. digitata and H. opuntia phylloids at high CO2 (1.01 per mil [parts per thousand] and 1.94 per mil, respectively) indicate increased photosynthetic uptake of CO2 over HCO3- potentially reducing dissolved inorganic carbon limitation at the seep site. Moreover, H. digitata and H. opuntia specimens transplanted for 14 d from the control to the seep site exhibited similar delta13C signatures as specimens grown there. These results suggest that the Halimeda spp. investigated can acclimatize and will likely still be capable to grow and calcify in inline image conditions exceeding most pessimistic future CO2 projections.