Late Holocene floodplain incision and alluvial fan formation in the central Grampian Highlands,Scotland: chronology,environment and implications

A section cut across an alluvial fan and the underlying floodplain terrace in the central Grampian Highlands provides an unusually complete record of late Holocene events. At ca. 2.7–2.4 cal kyr BP floodplain aggradation was replaced by net floodplain incision. Pollen evidence and charcoal counts provide no evidence for contemporaneous anthropogenic landscape change, and the timing of the transition suggests that it reflects an increase in high-magnitude erosive flood events following overall climatic deterioration. The overlying fan was deposited by torrential hyperconcentrated flows during three brief storm-generated depositional events at ca. 2.2–2.1, 1.9–1.8 and 0.9–0.7 cal kyr BP, separated and succeeded by prolonged periods of stability and peat accumulation. During these three events, a cumulative total of ca. 6750 m³ of sediment was deposited, probably in no more than a few hours over a timescale of two millennia. These findings imply that proposed links between human activity and the development of alluvial fans or debris cones require reassessment, and that different elements of the Holocene alluvial landscape have responded in different ways to the same climatic inputs. Aggregation of dating evidence relating to aggradation or incision of alluvial landforms at different scales therefore may produce misleading results. Copyright © 1999 John Wiley & Sons, Ltd.     

Probabilistic domain decomposition for the solution of the two-dimensional magnetotelluric problem

Probabilistic domain decomposition is proposed as a novel method for solving the two-dimensional Maxwell’s equations as used in the magnetotelluric method. The domain is split into non-overlapping sub-domains and the solution on the sub-domain boundaries is obtained by evaluating the stochastic form of the exact solution of Maxwell’s equations by a Monte-Carlo approach. These sub-domains can be naturally chosen by splitting the sub-surface domain into regions of constant (or at least continuous) conductivity. The solution over each sub-domain is obtained by solving Maxwell’s equations in the strong form. The sub-domain solver used for this purpose is a meshless method resting on radial basis function-based finite differences. The method is demonstrated by solving a number of classical magnetotelluric problems, including the quarter-space problem, the block-in-half-space problem and the triangle-in-half-space problem.     

Sensitivity of the elastic anisotropy and seismic reflection amplitude of the Eagle Ford Shale to the presence of kerogen

The Eagle Ford Shale of Central and South Texas is currently of great interest for oil and gas exploration and production. Laboratory studies show that the Eagle Ford Shale is anisotropic, with a correlation between anisotropy and total organic carbon. Organic materials are usually more compliant than other minerals present in organic‐rich shales, and their shapes and distribution are usually anisotropic. This makes organic materials an important source of anisotropy in organic‐rich shales. Neglecting shale anisotropy may lead to incorrect estimates of rock and fluid properties derived from inversion of amplitude versus offset seismic data. Organic materials have a significant effect on the PP and PS reflection amplitudes from the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface. The higher kerogen content of the Lower Eagle Ford compared with that of the Upper Eagle Ford leads to a negative PP reflection amplitude that dims with offset, whereas the PS reflection coefficient increases in magnitude with increasing offset. The PP and PS reflection coefficients at the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface all increase in magnitude with increasing volume fraction of kerogen.     

Development of ST:REAM: a reach‐based stream power balance approach for predicting alluvial river channel adjustment

River channel sediment dynamics are important in integrated catchment management because changes in channel morphology resulting from sediment transfer have important implications for many river functions. However, application of existing approaches that account for catchment‐scale sediment dynamics has been limited, largely due to the difficulty in obtaining data necessary to support them. It is within this context that this study develops a new, reach‐based, stream power balance approach for predicting river channel adjustment. The new approach, named ST:REAM (sediment transport: reach equilibrium assessment method), is based upon calculations of unit bed area stream power (ω) derived from remotely sensed slope, width and discharge datasets. ST:REAM applies a zonation algorithm to values of ω that are spaced every 50 m along the catchment network in order to divide the branches of the network up into relatively homogenous reaches. ST:REAM then compares each reach's ω value with the ω of its upstream neighbour in order to predict whether or not the reach is likely to be either erosion dominated or deposition dominated. The paper describes the application of ST:REAM to the River Taff in South Wales, UK. This test study demonstrated that ST:REAM can be rapidly applied using remotely sensed data that are available across many river catchments and that ST:REAM correctly predicted the status of 87.5% of sites within the Taff catchment that field observations had defined as being either erosion or deposition dominated. However, there are currently a number of factors that limit the usefulness of ST:REAM, including inconsistent performance and the need for additional, resource intensive, data to be collected to both calibrate the model and aid interpretation of its results. Copyright © 2014 John Wiley & Sons, Ltd.     

Quantifying potential propagule pressure of aquatic invasive species from the commercial shipping industry in Canada

We quantify and compare different measures of potential propagule pressure (PPP) of aquatic invasive species (AIS) from commercial vessels in Canada. We used ship arrivals and ballast water discharge volumes as proxies for PPP from ballast water organisms, and wetted surface area (WSA) as a proxy for hull fouling PPP, to determine their relative contributions to total PPP. For three regions studied, PPP proxies correlated significantly across ports and some vessel categories. Relative contributions of ship arrivals, ballast discharge, and WSAs to PPP, evidenced by non-significant correlations across these measures, varied across regions, ports, vessel types, and seasons. Flow-through (dominant on east and west coasts) and empty-refill (in Great Lakes-St. Lawrence region) were the major ballast water exchange methods employed by the vessels surveyed. These methods have different biological efficacy for AIS removal, influencing PPP. Our study illustrates benefits and limitations of using different PPP proxies to estimate invasion risk.     

Setting new thermal standards for transitional and coastal (TraC) waters

The construction of a new generation of coastal power stations in the UK and other western European countries has highlighted the absence of robust standards for thermal discharges to transitional and coastal waters. The proposed discharge volumes are greater than hitherto, yet there has been little independent critical examination of their potential impact, whilst much of the existing guidance has been adapted from freshwater practice. This review considers the available knowledge on the tolerance and behaviour of fish and other marine biota to heated effluents. Appropriate ways are suggested of grouping fish species to reflect their sensitivity to thermal effects. The plethora of existing standards are considered and their validity assessed in a framework of predicted seawater temperature rise. Those species or groups of organisms most likely to be affected are identified and finally specific recommendations for thermal standards consistent with long term sustainability are proposed.     

Modelling converted seismic waveforms in isotropic and anisotropic 1-D gradients: discontinuous versus continuous gradient representations

Over the past decade, there have been numerous receiver function studies directed at imaging the lithosphere-asthenosphere boundary (LAB). Although it is generally accepted that receiver function phases observed in these studies are derived from physical mode conversions at depth within the lithosphere-asthenosphere transition, it is still debatable as to whether these phases are directly indicative of the LAB. This is because interpretation of receiver function LAB signals relies on understanding the elastic characteristics of the Earth??s outer thermal boundary layer. The main issues for receiver function imaging are the sharpness of the elastic material property transition and, more importantly, what specifically are the material gradients. To test the various transition models, a forward modelling approach is required that allows accurate waveform synthetics for a range of discontinuous and continuous gradients in anisotropic, elastic media. We present a derivation of the reflection and transmission response for continuous one-dimensional (1-D) gradients in generally anisotropic elastic media. We evaluate the influence of 1-D isotropic and anisotropic elastic gradients on the seismic waveform by comparing numerical results of models for discontinuous and continuous transitions. The results indicate that discontinuous representations using layers each with uniform parameters and with thicknesses on the order of approximately 1/3 to 1/8 of the dominant seismic wavelength can be used to accurately model P-to-S and S-to-P mode conversions due to continuous transitions of both isotropic and anisotropic elastic properties. From a practical point of view, when comparing synthetic modelling with observation, this constraint can be relaxed further. The presence of signal noise and/or the result of receiver function stacking techniques will likely obscure these subtle waveform e ff ects. Hence this study suggests that accurate synthetic waveforms for LAB transitions can be modelled with discontinuous gradient representations using a reasonable number of discrete transition layers with layer thicknesses no greater than 1/2 to 1/3 the dominant seismic wavelength.     

Aminostratigraphy and thermoluminescence dating of coastal aeolianites and the later Quaternary history of a failed delta: The River Murray mouth region,South Australia

A geochronological framework for the sequential development of coastal barrier aeolianite complexes in the mouth region of the River Murray, Australia's largest river system is presented based on amino acid racemization and thermoluminescence dating. The sedimentary successions represent a foreshortened and condensed sequence of coastal barriers compared with those of the Coorong Coastal Plain in southern South Australia where the barrier complexes are more widely separated in response to tectonic uplift. The barriers have formed during interglacial sea-level highstands and are correlatives of genetically equivalent landforms of the Coorong Coastal Plain. Thermoluminescence dating and the extent of amino acid racemization in aeolianite ‘whole-rock’ sediment samples, reveal a general increase in age of the barriers landwards from the modern coastline. In detail, however, the individual barriers represent composite structures having formed in more than one interglaciation, due to the reoccupation of Pleistocene shoreline positions during sea-level highstands of similar amplitude, in a zone of gradual basin subsidence. The most seaward Pleistocene aeolianite at Surfer Beach is of interstadial age (Marine Isotope Stage 5c, 105 ± 5 ka; MIS 5c), and correlates with the Robe Range of the Coorong Coastal Plain. The last interglacial shoreline (130 ± 15 ka; MIS 5e) is particularly well-defined in the River Murray mouth region. It is represented by a complex association of coastal parabolic dunes superimposed on a transverse dune system, which runs parallel with the former coastline, and also includes associated estuarine, lagoonal and open ocean beach facies. Landward of the last interglacial succession are distinct barriers relating to the penultimate interglaciation (215 ± 35 ka; MIS 7), as well as earlier interglaciations (350 ± 65 ka; MIS 9 or 11 and 470 ± 70 ka; MIS 11 or 13). The coastal barriers have been successively breached by the ancestral River Murray at times of lower sea level during glacial cycles. Former mouths of the River Murray during interglacial sea-level highstands are likely to have existed near Tauwitchere Island during MIS 7, and between Goolwa and Hindmarsh Island and near the southern-most part of Lake Albert during the last interglacial (MIS 5e). The River Murray mouth region represents a failed delta as the limited sediment brought to this area since late middle Pleistocene time has been either rapidly incorporated within aeolian deposits during sea-level highstands, or transported to the edge of the Lacepede Shelf during glacial maxima. The Holocene and modern River Murray has not established a marine delta, but deposits its load in the settling basins of the terminal lakes. Only a small digitate delta has formed where the river enters Lake Alexandrina.     

Shaking‐table tests of flat‐bottom circular silos containing grain‐like material

According to Eurocode 8, the seismic design of flat‐bottom circular silos containing grain‐like material is based on a rough estimate of the inertial force imposed on the structure by the ensiled content during an earthquake: 80% of the mass of the content multiplied by the peak ground acceleration. A recent analytical consideration of the horizontal shear force mobilised within the ensiled material during an earthquake proposed by some of the authors has resulted in a radically reduced estimate of this load suggesting that, in practice, the effective mass of the content is significantly less than that specified. This paper describes a series of laboratory tests that featured shaking table and a silo model, which were conducted in order to obtain some experimental data to verify the proposed theoretical formulations and to compare with the established code provisions. Several tests have been performed with different heights of ensiled material – about 0.5 mm diameter Ballotini glass – and different magnitudes of grain–wall friction. The results indicate that in all cases, the effective mass is indeed lower than the Eurocode specification, suggesting that the specification is overly conservative, and that the wall–grain friction coefficient strongly affects the overturning moment at the silo base. At peak ground accelerations up to around 0.35 g, the proposed analytical formulation provides an improved estimate of the inertial force imposed on such structures by their contents. Copyright © 2015 John Wiley & Sons, Ltd.