Flow separation and roughness lengths over large bedforms in a tidal environment

This study characterises the shape of the flow separation zone (FSZ) and wake region over large asymmetric bedforms under tidal flow conditions. High resolution bathymetry, flow velocity and turbulence data were measured along two parallel transects in a tidal channel covered with bedforms. The field data are used to verify the applicability of a numerical model for a systematic study using the Delft3D modelling system and test the model sensitivity to roughness length. Three experiments are then conducted to investigate how the FSZ size and wake extent vary depending on tidally-varying flow conditions, water levels and bathymetry. During the ebb, a large FSZ occurs over the steep lee side of each bedform. During the flood, no flow separation develops over the bedforms having a flat crest; however, a small FSZ is observed over the steepest part of the crest of some bedforms, where the slope is locally up to 15°. Over a given bedform morphology and constant water levels, no FSZ occurs for velocity magnitudes smaller than 0.1 m s**-1; as the flow accelerates, the FSZ reaches a stable size for velocity magnitudes greater than 0.4 m s**-1. The shape of the FSZ is not influenced by changes in water levels. On the other hand, variations in bed morphology, as recorded from the high-resolution bathymetry collected during the tidal cycle, influence the size and position of the FSZ: a FSZ develops only when the maximum lee side slope over a horizontal distance of 5 m is greater than 10°. The height and length of the wake region are related to the length of the FSZ. The total roughness along the transect lines is an order of magnitude larger during the ebb than during the flood due to flow direction in relation to bedform asymmetry: during the ebb, roughness is created by the large bedforms because a FSZ and wake develops over the steep lee side. The results add to the understanding of hydrodynamics of natural bedforms in a tidal environment and may be used to better parameterise small-scale processes in large-scale studies.

Predicting bed form roughness: the influence of lee side angle

Flow transverse bedforms (ripples and dunes) are ubiquitous in rivers and coastal seas. Local hydrodynamics and transport conditions depend on the size and geometry of these bedforms, as they constitute roughness elements at the bed. Bedform influence on flow energy must be considered for the understanding of flow dynamics, and in the development and application of numerical models. Common estimations or predictors of form roughness (friction factors) are based mostly on data of steep bedforms (with angle-of-repose lee slopes), and described by highly simplified bedform dimensions (heights and lengths). However, natural bedforms often are not steep, and differ in form and hydraulic effect relative to idealised bedforms. Based on systematic numerical model experiments, this study shows how the hydraulic effect of bedforms depends on the flow structure behind bedforms, which is determined by the bedform lee side angle, aspect ratio and relative height. Simulations reveal that flow separation behind bedform crests and, thus, a hydraulic effect is induced at lee side angles steeper than 11 to 18° depending on relative height, and that a fully developed flow separation zone exists only over bedforms with a lee side angle steeper than 24°. Furthermore, the hydraulic effect of bedforms with varying lee side angle is evaluated and a reduction function to common friction factors is proposed. A function is also developed for the Nikuradse roughness (k s), and a new equation is proposed which directly relates k s to bedform relative height, aspect ratio and lee side angle.

Flow conditions, bed form dimensions and shear stress experiments

Large asymmetric bed forms commonly develop in rivers. The turbulence associated with flow separation that develops over their steep lee side is responsible for the form shear stress which can represent a substantial part of total shear stress in rivers. This paper uses the Delft3D modeling system to investigate the effects of bed form geometry and forcing conditions on flow separation length and associated turbulence, and bed form shear stress over angle-of-repose (30 lee side angle) bed forms. The model was validated with lab measurements that showed sufficient agreement to be used for a systematic analysis. The influence of flow velocity, bed roughness, relative height (bed form height/water depth), and aspect ratio (bed form height/length) on the variations of the normalized length of the flow separation zone, the extent of the wake region (where the turbulent kinetic energy (TKE) was more than 70% of the maximum TKE), the average TKE within the wake region and the form shear stress were investigated. Form shear stress was found not to scale with the size of the flow separation zone but to be related to the product of the normalized extent of the wake region (extent of the wake region/extent of water body above the bed form) and the average TKE within the wake region. The results add to understanding of the hydrodynamics of bed forms and may be used for the development of better parameterizations of smallscale processes for application in large-scale studies.

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.

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.

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.

Distribution of biological, anthropogenic, and volcanic constituents in the San Francisco Bay Coastal System

Although conventional sediment parameters (mean grain size, sorting, and skewness) and provenance have typically been used to infer sediment transport pathways, most freshwater, brackish, and marine environments are also characterized by abundant sediment constituents of biological, and possibly anthropogenic and volcanic, origin that can provide additional insight into local sedimentary processes. The biota will be spatially distributed according to its response to environmental parameters such as water temperature, salinity, dissolved oxygen, organic carbon content, grain size, and intensity of currents and tidal flow, whereas the presence of anthropogenic and volcanic constituents will reflect proximity to source areas and whether they are fluvially- or aerially-transported. Because each of these constituents have a unique environmental signature, they are a more precise proxy for that source area than the conventional sedimentary process indicators. This San Francisco Bay Coastal System study demonstrates that by applying a multi-proxy approach, the primary sites of sediment transport can be identified. Many of these sites are far from where the constituents originated, showing that sediment transport is widespread in the region. Although not often used, identifying and interpreting the distribution of naturally-occurring and allochthonous biologic, anthropogenic, and volcanic sediment constituents is a powerful tool to aid in the investigation of sediment transport pathways in other coastal systems.

Coral-algae metabolism and diurnal changes in the CO2-carbonate system of bulk sea water

Precise measurements were conducted in continuous flow seawater mesocosms located in full sunlight that compared metabolic response of coral, coral-macroalgae and macroalgae systems over a diurnal cycle. Irradiance controlled net photosynthesis (Pnet), which in turn drove net calcification (Gnet), and altered pH. Pnet exerted the dominant control on [CO3]2- and aragonite saturation state (Omega arag) over the diel cycle. Dark calcification rate decreased after sunset, reaching zero near midnight followed by an increasing rate that peaked at 03:00 h. Changes in Omega arag and pH lagged behind Gnet throughout the daily cycle by two or more hours. The flux rate Pnet was the primary driver of calcification. Daytime coral metabolism rapidly removes dissolved inorganic carbon (DIC) from the bulk seawater and photosynthesis provides the energy that drives Gnet while increasing the bulk water pH. These relationships result in a correlation between Gnet and Omega arag, with Omega arag as the dependent variable. High rates of H+ efflux continued for several hours following mid-day peak Gnet suggesting that corals have difficulty in shedding waste protons as described by the Proton Flux Hypothesis. DIC flux (uptake) followed Pnet and Gnet and dropped off rapidly following peak Pnet and peak Gnet indicating that corals can cope more effectively with the problem of limited DIC supply compared to the problem of eliminating H+. Over a 24 h period the plot of total alkalinity (AT) versus DIC as well as the plot of Gnet versus Omega arag revealed a circular hysteresis pattern over the diel cycle in the coral and coral-algae mesocosms, but not the macroalgae mesocosm. Presence of macroalgae did not change Gnet of the corals, but altered the relationship between Omega arag and Gnet. Predictive models of how future global changes will effect coral growth that are based on oceanic Omega arag must include the influence of future localized Pnet on Gnet and changes in rate of reef carbonate dissolution. The correlation between Omega arag and Gnet over the diel cycle is simply the response of the CO2-carbonate system to increased pH as photosynthesis shifts the equilibria and increases the [CO3]2- relative to the other DIC components of [HCO3]- and [CO2]. Therefore Omega arag closely tracked pH as an effect of changes in Pnet, which also drove changes in Gnet. Measurements of DIC flux and H+ flux are far more useful than concentrations in describing coral metabolism dynamics. Coral reefs are systems that exist in constant disequilibrium with the water column.

Heavy mineral analyses of sediment samples from the San Francisco Bay coastal system

Heavy or high-specific gravity minerals make up a small but diagnostic component of sediment that is well suited for determining the provenance and distribution of sediment transported through estuarine and coastal systems worldwide. By this means, we see that surficial sand-sized sediment in the San Francisco Bay Coastal System comes primarily from the Sierra Nevada and associated terranes by way of the Sacramento and San Joaquin Rivers and is transported with little dilution through the San Francisco Bay and out the Golden Gate. Heavy minerals document a slight change from the strictly Sierran-Sacramento mineralogy at the confluence of the two rivers to a composition that includes minor amounts of chert and other Franciscan Complex components west of Carquinez Strait. Between Carquinez Strait and the San Francisco Bar, Sierran sediment is intermingled with Franciscan-modified Sierran sediment. The latter continues out the Gate and turns southward towards beaches of the San Francisco Peninsula. The Sierran sediment also fans out from the San Francisco Bar to merge with a Sierran province on the shelf in the Gulf of the Farallones. Beach-sand sized sediment from the Russian River is transported southward to Point Reyes where it spreads out to define a Franciscan sediment province on the shelf, but does not continue southward to contribute to the sediment in the Golden Gate area.