Modelling the decadal dynamics of reach-scale river channel evolution and floodplain turnover in CAESAR-Lisflood

Sedimentary deposits provide records of environmental change quantifying erosion fluxes conditioned by natural and anthropogenic disturbances. These fluxes are lagged by internal storage, particularly within floodplains, complicating reconstruction of environmental changes. The time sediment remains in storage underpins the interpretation of sedimentary records and accurate monitoring of pollutant fluxes. Turnover time is a measure of the timeframe to erode every floodplain surface. CAESAR-Lisflood is used to simulate fluvial evolution at reach scale, providing a basis for quantifying environmental changes on the timescales of sediment storage. We evaluate the accuracy of CAESAR-Lisflood simulations of channel changes and turnover times for alluvial floodplains using historical channel changes reconstructed for 10 reaches in northern England to quantify model accuracy in replicating mean annual erosion, deposition and channel lateral migration rates, alongside planform morphology. Here, a split-sample testing approach is adopted, whereby five of the reaches were calibrated and the resulting parameter values were applied to the other reaches to evaluate the transferability of parameter settings. The lowest overall integrated error identified the best-fit simulations and showed that modelled process rates were within ~25–50% of rates from historical reconstructions, generally. Calibrated parameters for some reaches are widely transferable, producing accurate geomorphic changes for some uncalibrated sites. However, large errors along some reaches indicate that reach-specific parameterization is recommended. Turnover times are underpinned by the assumption that areas of floodplain previously unvisited by the channel are reworked. This assumption has been challenged by studies that show floodplain (re)occupation rates vary spatially. However, this limitation is less important for the short-duration simulations presented here. The simulations reconstruct floodplain turnover times estimated by mapped rates mostly successfully, demonstrating the potential applicability of calibrated parameters over much longer timescales. Errors in the form of under-predicted erosion rates propagated, resulting in over-predicted turnover times by even greater magnitudes. © 2020 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Overwash-dominated stratigraphy of barriers with intermittent inlets

The morphodynamics and structure of barriers with persistent tidal inlets have been well studied. In contrast the stratigraphy and functioning of barrier systems with ephemeral inlets is poorly understood. This article examines the barrier-inlet systems of two intermittently closed open lagoons or temporarily open closed estuaries on the east coast of South Africa. Multiple geophysical surveys using ground-penetrating radar (GPR) were correlated with exposed sections of the barrier where inlet formation revealed the internal stratigraphy. Stratigraphic observations were placed in the context of the contemporary wave dynamics and mesoscale geomorphic evolution. The integrated databases reveal an absence of migrating channel features. Instead the stratigraphy is dominated by landward dipping sheets of alternating high- and low-amplitude reflectors. These correlate with gravel, shell debris and heavy mineral-lined beds formed by overwash processes. Where ephemeral inlet structures are preserved in the stratigraphy, their fills comprise aggrading, high-amplitude reflectors, linked to washover infilling of the inlet mouth. Multiple small channels in the more distal portions of the barrier in georadar stratigraphy are related to channelized washover flow. These barriers often breach during high swell and are subsequently sealed during fairweather wave conditions. Time series analysis of waves and satellite imagery shows a link between storms from the south and breach events. This is consistent with the truncations in subsurface images and inferred barrier lowering by overwash channelling. These barriers experience quasi-stable oscillations in their landward and seaward shore position, punctuated by periods of barrier rollover associated with the most intense storms. As overwash is responsible in part for both the constructive and destructive phases of the barrier, these barriers have low preservation potential. Persistent rollover driven by overwashing will terminate once accommodation space is eliminated and the barriers are eroded by storm activity. © 2019 John Wiley & Sons, Ltd.     

Edgar Hynes Thompson: pragmatism, ingenuity and mechanisms

In this contribution the author looks for an explanation of an uncharacteristically scathing review, by E. H. Thompson towards the end of his life, of a book, Theory of Errors and Generalized Matrix Inverses, by Arne Bjerhammar, a leading European academic.     

Terracing potential field data

Terracing is an operator that is applied to potential field data to produce regions of constant field amplitude that are separated by sharp boundaries, as an aid to interpretation. When applied to map data, the boundaries are defined by the zero contour of the 2D Laplacian derivative operator. An improved method is described here that defines the boundaries by the zero contour of the profile curvature. This approach gives superior results because the 2D Laplacian operator is composed only of derivatives in the EW and NS directions, while the profile curvature uses the curvature in the 'uphill' direction at each point, whatever that direction may be. The method is demonstrated on gravity data from South Africa. Source code in Matlab format is available from the authors on request.     

The simulation of SST, sea ice extents and ocean heat transports in a version of the Hadley Centre coupled model without flux adjustments

Results are presented from a new version of the Hadley Centre coupled model (HadCM3) that does not require flux adjustments to prevent large climate drifts in the simulation. The model has both an improved atmosphere and ocean component. In particular, the ocean has a 1.25° × 1.25° degree horizontal resolution and leads to a considerably improved simulation of ocean heat transports compared to earlier versions with a coarser resolution ocean component. The model does not have any spin up procedure prior to coupling and the simulation has been run for over 400 years starting from observed initial conditions. The sea surface temperature (SST) and sea ice simulation are shown to be stable and realistic. The trend in global mean SST is less than 0.009 °C per century. In part, the improved simulation is a consequence of a greater compatibility of the atmosphere and ocean model heat budgets. The atmospheric model surface heat and momentum budget are evaluated by comparing with climatological ship-based estimates. Similarly the ocean model simulation of poleward heat transports is compared with direct ship-based observations for a number of sections across the globe. Despite the limitations of the observed datasets, it is shown that the coupled model is able to reproduce many aspects of the observed heat budget. Received: 1 October 1998 / Accepted: 20 July 1999     

Islands in a sea of sand: use of Acacia trees by tree skinks in the Kalahari Desert

In the Kalahari desert the tree skink Mabuya striata sparsa occupies Acacia trees separated by areas of open sand. On a single day in austral spring, the effects of tree size and structural complexity, presence of colonial nests of weavers, and fallen limbs on lizard abundance, and of time of day and air temperature on numbers of lizards active, location on trees, ground, or ground cover, and occupation of sunny, shaded, or partially shaded sites were examined. Multiple regression showed that the numbers of adult and total lizards increased with tree size, amount of ground cover, and air temperature. Neither number and size of sociable weaver nests nor complexity of tree shape affected lizard abundance. Proportions of lizards on trees and ground did not vary with time of day; neither did height above ground among individuals on trees. The proportions of lizards fully exposed to the sun was higher early and late in the day than at midday, when more lizards occupied shaded sites.     

Petrogenesis of carbonate-bearing nepheline syenites and carbonatites from Southern Victoria Land, Antarctica: origin of carbon and the effects of calcite-graphite equilibrium

Chemically-evolved carbonate-bearing nepheline syenites are intruded into basement metasediments of the Koettlitz Group on Dismal and Radian Ridges in the Pipecleaner Glacier region of Southern Victoria Land, Antarctica. Whole rock XRF data from the Dismal Nepheline Syenite defines a broad trend which is consistent with the removal of a cumulate fraction of approximate composition 70% hedenbergite, 15% nepheline, 10% titanite and 5% apatite. Stable isotope, major and trace element and mineralogical characteristics of the syenites are very similar to those of cross-cutting calcite-rich dykes indicating derivation from closely-related source magmas. The general association of carbonatites and nepheline syenites with extensional environments, suggests that the Dismal and Radian Ridge nepheline syenites and carbonatite dykes indicate a period of early Paleozoic (531 Ma) rifting or intrusion into localised tensional structures in an overall compressional regime. Assimilation of marble by the syenite magmas is evidenced by abundant rafts and xenoliths within the Dismal Nepheline Syenite, however, carbon and oxygen isotopic ratios from syenite and carbonatite calcites are distinctly lighter than values from the marble country rock indicating a magmatic source. Graphite and calcite commonly occur as aggregates in the Dismal Nepheline Syenite suggesting equilibrium between these two carbon-bearing phases. Isotopic fractionation between calcite and graphite, via the equilibrium: C + O₂ = CO₂ has resulted in enrichment of the ¹³C isotope in calcites from the Dismal Nepheline Syenite.     

Shipboard trials of an ozone-based ballast water treatment system

Legislation introduced by the United Nations International Maritime Organization (IMO) has focused primarily on standards defining successful treatments designed to remove invasive species entrained in ballast water. An earlier shipboard study found that ozone introduced into salt water ballast resulted in the formation of bromine compounds, measured as total residual oxidants (TRO) that were toxic to both bacteria and plankton. However, the diffuser system employed to deliver ozone to the ballast water tanks resulted in patchiness in TRO distribution and toxicity to entrained organisms. In this follow-up study, the shipboard diffuser system was replaced by a single Venturi-type injection system designed to deliver a more homogeneous biocide distribution. Within-tank variability in TRO levels and associated toxicity to zooplankton, phytoplankton and marine bacteria was measured via a matrix of tubes deployed to sample different locations in treated and untreated (control) tanks. Three trials were conducted aboard the oil tanker S/T Prince William Sound in the Strait of Juan de Fuca off Port Angeles, Washington State, USA, between June and December 2007. Mortalities of plankton and bacteria and oxidant concentrations were recorded for treated and untreated ballast water up to 3 days following treatment, and residual toxicity beyond this period was measured by bioassay of standard test organisms. Results indicated uniform compliance with current IMO standards, but only partial compliance with other existing and pending ballast water legislation.