Back to the future: reflections and directions of South African marine bioinvasion research

Biological invasions continue to increase around the world, with impacts on many coastal marine systems. Here we review the South African marine invasion literature which, despite the field being relatively new, has grown to have significant presence in both the local and international arenas. Of the 79 papers reviewed, 70% focused on the establishment and spread of alien species, with modes of transport and introduction largely overlooked. An emphasis was also apparent towards field studies, in particular survey work, with few experimental studies. The overwhelming majority of papers focused on a single species, the Mediterranean mussel Mytilus galloprovincialis, reflecting the scale of this invasion and the tractable nature of rocky shores as study systems. With the exception of this one species, the impacts of marine alien species have rarely been quantified. We suggest that future research extends the taxonomic coverage of present work and develops a better understanding of the mechanisms of introduction, establishment and spread of marine alien species. Through an experimental approach, the drivers of altered ecological patterns and processes resulting from invasions should be addressed, providing insight into associated impacts. This approach will maintain the local applicability and international relevance of South African marine invasion research.     

Dynamics and role of the Durban cyclonic eddy in the KwaZulu-Natal Bight ecosystem

The semi-permanent Durban Eddy is a mesoscale, lee-trapped, cold-core cyclonic circulation that occurs off the east coast of South Africa between Durban in the north and Sezela, some 70 km to the south. When present, strong north-eastward countercurrents reaching 100 cm s^–1 are found inshore. It is hypothesised that the cyclone is driven by the strong south-westward flowing Agulhas Current offshore of the regressing shelf edge near Durban. Analysis of ADCP data and satellite imagery shows the eddy to be present off Durban approximately 55% of the time, with an average lifespan of 8.6 days, and inter-eddy periods of 4 to 8 days. After spin-up the eddy breaks loose from its lee position and propagates downstream on the inshore boundary of the Agulhas Current. The eddy is highly variable in occurrence, strength and downstream propagation speeds. There is no detectable seasonal cycle in eddy occurrence, with the Natal Pulse causing more variability than any seasonal signal. A thermistor array deployed in the eddy centre, together with ship CTD data, indicates upward doming of the thermal structure in the eddy core associated with cooler water and nutrients being moved higher in the water column, stimulating primary production. Together with the use of satellite imagery, our findings indicate a second mechanism of upwelling, viz. divergent upwelling in the northern limb of the eddy. Satellite-tracked surface drifters released in the eddy demonstrated the potential for nutrient-rich eddy water to be transported northwards along the inshore regions of the KwaZulu-Natal (KZN) Bight, thus contributing to the functioning of the bight ecosystem, as well as southwards along the KZN and Transkei coasts – both by the eddy migrating downstream and by eddy water being recirculated into the inshore boundary of the Agulhas Current itself.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

Defining the Potential of Nanoscale Re‐Os Isotope Systematics Using Atom Probe Microscopy

Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm³). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS ¹⁸⁷Os/¹⁸⁹Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM ¹⁸⁷Os/¹⁸⁹Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles.     

Textural and micro-analytical insights into mafic–felsic interactions during the Oruanui eruption,Taupo

Quenched juvenile mafic inclusions (enclaves) are an occasional but informative component in the deposits of large felsic eruptions. Typically, the groundmasses of these inclusions rapidly crystallize as the mafic magma is chilled against a more voluminous, cooler felsic host, providing a physical and chemical record of the nature and timing of mafic–felsic interactions. We examine mafic inclusions of two compositional lineages (tholeiitic and calc-alkaline) from deposits of the 25.4 ka Oruanui eruption (Taupo, New Zealand). 2-D quantitative textural data from analysis of back-scattered electron images reveal a marked diversity in the groundmass textures of the inclusions, including median crystal sizes (amphibole: 14–45 µm; plagioclase: 21–75 µm) and aspect ratios (amphibole: 1.7–4.2; plagioclase: 2.1–4.0), area number densities (amphibole: 122–2660 mm^?2; plagioclase: 117–2990 mm^?2), area fractions (?) of minerals (?_plag?=?23–45%, ?_amph?=?0–28%, ?_cpx?=?0–6%, ?_oxides?=?0.6–5.5%), and the relative abundance of plagioclase and amphibole (?_plag/?_amph?=?1.0–4.6). Textural parameters vary more significantly within, rather than between, the two compositional lineages, and in some cases show marked variations across individual clasts, implying that each inclusion’s cooling history, rather than bulk composition, was the dominant control on textural development. Groundmass mineral compositions are also diverse both within and between inclusions (e.g. plagioclase from An_34–92, with typical intra-clast variability of ~?20 mol%), and do not correlate with bulk chemistry. Diverse groundmass textures and mineral and glass chemistries are inferred to reflect complex interplay of a range of factors including the degree and rate of undercooling, bulk composition, water content and, possibly, intensive variables. Our data are inconsistent with breakup of a crystallizing ponded mafic layer at the base of the Oruanui melt-dominant body, instead implying that each inclusion partially crystallized as a discrete body with a unique cooling history. Extensive ingestion of mush-derived macro-crystals suggests that mechanical breakup of mafic feeder dikes occurred within a transition zone between the mush and melt-dominant magma body. In this zone, the mush lacked yield strength, as has been inferred from field studies of narrow (meters to few tens of meters) mush-melt transition zones preserved in composite intrusions. Evidence for plastic deformation of inclusions during eruption and the abundance of fresh residual glass in inclusions from all eruptive phases suggest that the inclusions formed syn-eruptively, and must have been formed recurrently at multiple stages throughout the eruption.     

Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

Analysis of stratigraphical sequences at Cocina Cave (Spain) using rare earth elements geochemistry

This study investigates the stratigraphical sequence of Cocina Cave (Spain) employing and testing for the first time the capability of rare earth elements as markers of human activities in caves. Located in Dos Aguas (Valencian Community, Spain), Cocina Cave is characterized by the presence of several Holocene archaeological deposits from the final Mesolithic to the present day and is a pivotal site for understanding the socio-ecological dynamics of the last hunter-gatherer inhabitants of the Iberian Peninsula and the transition to pastoral and agricultural economies in the Western Mediterranean. However, the identification of strata from particular time-periods in the cave is often difficult due to the homogeneity of layers, the poor archaeological record in some strata and the presence of severe disturbance phenomena. The methodological approach of this study consisted of cross-referencing rare earth elements and other chemical markers with the archaeological stratigraphical interpretation, in an attempt to not only support the identification of the anthropic contribution to the formation of Cocina Cave strata, but also to characterize and confirm different natural and occupational episodes, particularly those associated with hunter-gatherer, early agriculturalist, and shepherd activities. Sediments were collected from different excavation areas and analysed for major elements, trace elements, rare earth elements (REE), soil organic matter (SOM) amounts and pH. Multivariate statistics were employed to group samples according to their elemental profile, and these were then compared to the archaeological temporal interpretation. The obtained results showed that REE amount and fractionation geochemical processes were regulated by carbonates, phosphates and pH. The use of REE as markers was particularly useful as their concentrations and their calculated ratios and anomaly distributions were demonstrated to be highly consistent with the archaeological stratigraphical interpretation.