Elevation adjustments of paired natural levees during flooding of the Saskatchewan River

Natural levees control the exchange of water between an alluvial channel and its floodplain, but little is known about the spatial distribution and evolution of levee heights. The summer 2005 flood of the Saskatchewan River (Cumberland Marshes, east‐central Saskatchewan) inundated large areas of floodplain for up to seven weeks, forming prominent new deposits on natural levees along main‐stem channels. Measurements of flood‐deposit thickness and crest heights of 61 levee pairs show that the thickest deposits occur on the lower pre‐flood levee in 80% of the sites, though no clear relationship exists between deposit thickness and magnitude of height difference. Only 16% of the pairs displayed thicker deposits on the higher levee, half of which occurred at sites where relatively clear floodbasin waters re‐entered turbid channels during general flooding. Difference in crest elevation (ΔE) between paired levees is approximately log‐normally distributed, both before and after the flood, though with different mean values. Supplemental observations from tank experiments indicate that during near‐bankfull flows, temporally and spatially variable deposition and erosion occur on levees due to backwater effects associated with nearby channel bars and irregular rises of the channel bed forced by channel extension. During floods, preferential deposition in lows tends to even out crest heights. Copyright © 2009 John Wiley & Sons, Ltd.     

Mass-dependent cadmium isotopic variations in nature with emphasis on the marine environment

We report a survey of natural mass-dependent cadmium isotope fractionation measured by thermal ionization mass spectrometry using a double-spike technique (DS-TIMS). Over sixty samples of natural terrestrial Cd from diverse environments, including MORB, OIB, continental loess, hydrogenic and hydrothermal ferromanganese deposits, and sphalerites (both oceanic and from major continental ore deposits) were analysed. Our results are expressed in terms of ε^112/110Cd, which are deviations in ¹¹²Cd/¹¹⁰Cd from our in-house JMC Cd standard in parts per 10⁴. The total ε^112/110Cd variation is relatively small, with a range of only 5 ε-units, and is one-to-two orders of magnitude smaller than that previously found in meteorites.The MORB, OIB and loess ε^112/110Cd values are similar and provide a good estimate for the bulk silicate Earth (BSE) value which is ? 0.95 ± 0.12 relative to our Cd standard (ε^112/110Cd = + 0.16 relative to Münster JMC Cd). Taken together, these data suggest little Cd isotope fractionation takes place during crust–mantle segregation. Cd isotopic compositions of continental sphalerite (ZnS) deposits worldwide and high-temperature oceanic hydrothermal sulphides show remarkably similar ε^112/110Cd values, consistent with our estimate for the BSE. In contrast, mid-temperature oceanic sulphides from a single extinct hydrothermal chimney display over 4 ε-units variation — along with the most negative values. These variations are most probably caused by precipitation/redissolution of sulphide phases en route within the hydrothermal system.The ε^112/110Cd variability found in worldwide marine Fe–Mn deposits reflects the seawater Cd isotope signal upon precipitation from ambient seawater. A decrease in ε^112/110Cd is observed in passing from shallow-water Fe–Mn deposits to those from deeper waters (> 2000 m depth). This shift is explained by biological fractionation related to the uptake of dissolved seawater Cd by phytoplankton in the upper water column. The relatively uniform ε^112/110Cd values close to zero at great depths are consistent with regeneration and remineralization of Cd at depth. Our data suggest that Cd isotopes – much like the Cd/Ca ratio in foraminifera – could potentially serve as a proxy for past changes in biological productivity. The temporal Cd isotope record in a Fe–Mn crust archive at 2000 m depth from the NE Atlantic suggests no gross long-term changes in Cd cycling took place over the past 8 Ma.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Frontal dynamics boost primary production in the summer stratified Mediterranean sea

Bio-physical glider measurements from a unique process-oriented experiment in the Eastern Alboran Sea (AlborEx) allowed us to observe the distribution of the deep chlorophyll maximum (DCM) across an intense density front, with a resolution (～ 400 m) suitable for investigating sub-mesoscale dynamics. This front, at the interface between Atlantic and Mediterranean waters, had a sharp density gradient (Δρ ～ 1 kg/m³ in ～ 10 km) and showed imprints of (sub-)mesoscale phenomena on tracer distributions. Specifically, the chlorophyll-a concentration within the DCM showed a disrupted pattern along isopycnal surfaces, with patches bearing a relationship to the stratification (buoyancy frequency) at depths between 30 and 60 m. In order to estimate the primary production (PP) rate within the chlorophyll patches observed at the sub-surface, we applied the Morel and Andrè (J Geophys Res 96:685–698 1991) bio-optical model using the photosynthetic active radiation (PAR) from Argo profiles collected simultaneously with glider data. The highest production was located concurrently with domed isopycnals on the fresh side of the front, suggestive that (sub-)mesoscale upwelling is carrying phytoplankton patches from less to more illuminated levels, with a contemporaneous delivering of nutrients. Integrated estimations of PP (1.3 g C m^?2d^?1) along the glider path are two to four times larger than the estimations obtained from satellite-based algorithms, i.e., derived from the 8-day composite fields extracted over the glider trip path. Despite the differences in spatial and temporal sampling between instruments, the differences in PP estimations are mainly due to the inability of the satellite to measure DCM patches responsible for the high production. The deepest (depth > 60 m) chlorophyll patches are almost unproductive and probably transported passively (subducted) from upper productive layers. Finally, the relationship between primary production and oxygen is also investigated. The logarithm of the primary production in the DCM interior (chlorophyll (Chl) > 0.5 mg/m³) shows a linear negative relationship with the apparent oxygen utilization, confirming that high chlorophyll patches are productive. The slope of this relationship is different for Atlantic, mixed interface waters and Mediterranean waters, suggesting the presence of differences in planktonic communities (whether physiological, population, or community level should be object of further investigation) on the different sides of the front. In addition, the ratio of optical backscatter to Chl is high within the intermediate (mixed) waters, which is suggestive of large phytoplankton cells, and lower within the core of the Atlantic and Mediterranean waters. These observations highlight the relevance of fronts in triggering primary production at DCM level and shaping the characteristic patchiness of the pelagic domain. This gains further relevance considering the inadequacy of optical satellite sensors to observe DCM concentrations at such fine scales.     

Occurrence and distribution of bacteria indicators,chemical tracers and pathogenic vibrios in Singapore coastal waters

Water quality in Singapore's coastal area was evaluated with microbial indicators, pathogenic vibrios, chemical tracers and physico-chemical parameters. Sampling sites were grouped into two clusters (coastal sites at (i) northern and (ii) southern part of Singapore). The coastal sites located at northern part of Singapore along the Johor Straits exhibited greater pollution. Principal component analysis revealed that sampling sites at Johor Straits have greater loading on carbamazepine, while turbidity poses greater influence on sampling sites at Singapore Straits. Detection of pathogenic vibrios was also more prominent at Johor Straits than the Singapore Straits. This study examined the spatial variations in Singapore's coastal water quality and provided the baseline information for health risk assessment and future pollution management.     

Evaluation of recent hydro-climatic changes in four tributaries of the Niger River Basin (West Africa)

West Africa experienced severe drought during the 1970s and 1980s, posing a threat to water resources. A wetter climate more recently suggests recovery from the drought. The Mann-Kendall trend and Theil-Sen’s slope estimator were applied to detect probable trends in weather elements in four sub-basins of the Niger River Basin between 1970 and 2010. The cross-entropy method was used to detect breakpoints in rainfall and runoff, Spearman’s rank test for correlation between the two, and cross-correlation analysis for possible lags. Results showed an overall increase in rainfall and runoff and a decrease in sunshine duration. Spearman’s coefficients suggest significant (5%) moderate to strong rainfall–runoff correlation for three sub-basins. A significant lower runoff was observed around 1979, with a rainfall break around 1992, indicating possible cessation of the drought. Temperatures increased significantly, at 0.02–0.05°C year^-1, with a negative wind speed trend for most stations. Half of the stations exhibited an increase in potential evapotranspiration.

EDITOR M.C. Acreman

ASSOCIATE EDITOR Not assigned     

Hydrological system analysis and modelling of the Kara River basin (West Africa) using a lumped metric conceptual model

This paper discusses the analysis and modelling of the hydrological system of the basin of the Kara River, a transboundary river in Togo and Benin, as a necessary step towards sustainable water resources management. The methodological approach integrates the use of discharge parameters, flow duration curves and the lumped conceptual model IHACRES. A Sobol sensitivity analysis is performed and the model is calibrated by applying the shuffled complex evolution algorithm. Results show that discharge generation in three nested catchments of the basin is affected by landscape physical characteristics. The IHACRES model adequately simulates the rainfall–runoff dynamics in the basin with a mean modified Nash-Sutcliffe efficiency measure of 0.6. Modelling results indicate that parameters controlling rainfall transformation to effective rainfall are more sensitive than those routing the streamflow. This study provides insights into understanding the catchment’s hydrological system. Nevertheless, further investigations are required to better understand detailed runoff generation processes.

EDITOR M.C. Acreman; ASSOCIATE EDITOR N Verhoest     

Geophysical flows under location uncertainty,Part I Random transport and general models

A stochastic flow representation is considered with the Eulerian velocity decomposed between a smooth large scale component and a rough small-scale turbulent component. The latter is specified as a random field uncorrelated in time. Subsequently, the material derivative is modified and leads to a stochastic version of the material derivative to include a drift correction, an inhomogeneous and anisotropic diffusion, and a multiplicative noise. As derived, this stochastic transport exhibits a remarkable energy conservation property for any realizations. As demonstrated, this pivotal operator further provides elegant means to derive stochastic formulations of classical representations of geophysical flow dynamics.     

Geophysical flows under location uncertainty,Part II Quasi-geostrophy and efficient ensemble spreading

Models under location uncertainty are derived assuming that a component of the velocity is uncorrelated in time. The material derivative is accordingly modified to include an advection correction, inhomogeneous and anisotropic diffusion terms and a multiplicative noise contribution. In this paper, simplified geophysical dynamics are derived from a Boussinesq model under location uncertainty. Invoking usual scaling approximations and a moderate influence of the subgrid terms, stochastic formulations are obtained for the stratified Quasi-Geostrophy and the Surface Quasi-Geostrophy models. Based on numerical simulations, benefits of the proposed stochastic formalism are demonstrated. A single realization of models under location uncertainty can restore small-scale structures. An ensemble of realizations further helps to assess model error prediction and outperforms perturbed deterministic models by one order of magnitude. Such a high uncertainty quantification skill is of primary interests for assimilation ensemble methods. MATLAB^® code examples are available online.     

A wind tunnel investigation of particle segregation,ripple formation and armouring within sand beds of systematically varied texture

This paper reports on a wind tunnel investigation of particle segregation, ripple formation and surface armouring within sand beds of systematically varied particle size distribution, from coarsely skewed to bimodal. By design, the system was closed with no external inputs of mass from an external particle feed. Particles too coarse to travel in saltation for the given range in wind speed were dyed red in order to distinguish them in optical images from finer sand particles, which could be entrained into the unidirectional airflow. A 3D laser scanner measured the changing bed topography at regular time intervals during 18 experiments involving varied combinations of wind speed and bed texture. Image classification techniques were used to investigate the coincident self‐organization of the two populations of particles, as distinguished by their colour. As soon as saltation commenced, some of the red particles segregated into thin discontinuous patches. Particle trapping and sheltering on these rough patches was strongly favoured, causing them to grow preferentially. During the earliest stages of formation, bedform growth coincided with: (i) rapid coarsening of the surface texture; and (ii) the merging of proto‐ripple ‘crests’ to generate larger rhythmic bedforms of lower frequency. Consistent with previous work, ripple size was observed to increase under stronger winds when not exceeding the threshold for entrainment of the coarse‐mode or red particles from the crest. With declining rates of mass transport and particle segregation as the bed surface armoured, and the consequent deceleration of ripple propagation through to the end of each experiment, all surfaces eventually attained a steady‐state morphometry. At saturation, the largest ripples developed on beds having the lowest initial concentration of red particles. Copyright © 2016 John Wiley & Sons, Ltd.