Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Middle Paleolithic occupation on a Marine Isotope Stage 5 lakeshore in the Nefud Desert,Saudi Arabia

Major hydrological variations associated with glacial and interglacial climates in North Africa and the Levant have been related to Middle Paleolithic occupations and dispersals, but suitable archaeological sites to explore such relationships are rare on the Arabian Peninsula. Here we report the discovery of Middle Paleolithic assemblages in the Nefud Desert of northern Arabia associated with stratified deposits dated to 75,000 years ago. The site is located in close proximity to a substantial relict lake and indicates that Middle Paleolithic hominins penetrated deeply into the Arabian Peninsula to inhabit landscapes vegetated by grasses and some trees. Our discovery supports the hypothesis of range expansion by Middle Paleolithic populations into Arabia during the final humid phase of Marine Isotope Stage 5, when environmental conditions were still favorable.     

Using monazite and zircon petrochronology to constrain the P–T–t evolution of the middle crust in the Bhutan Himalaya

The growth and dissolution behaviour of accessory phases (and especially those of geochronological interest) in metamorphosed pelites depends on, among others, the bulk composition, the prograde metamorphic evolution and the cooling path. Monazite and zircon are arguably the most commonly used geochronometers for dating felsic metamorphic rocks, yet crystal growth mechanisms as a function of rock composition, pressure and temperature are still incompletely understood. Ages of different growth zones in zircon and monazite in a garnet‐bearing anatectic metapelite from the Greater Himalayan Sequence in NW Bhutan were investigated via a combination of thermodynamic modelling, microtextural data and interpretation of trace‐element chemical ‘fingerprint’ indicators in order to link them to the metamorphic stage at which they crystallized. Differences in the trace‐element composition (HREE, Y, Eu_N/Eu*_N) of different phases were used to track the growth/dissolution of major (e.g. plagioclase, garnet) and accessory phases (e.g. monazite, zircon, xenotime, allanite). Taken together, these data constrain multiple pressure–temperature–time (P–T–t) points from low temperature (<550 °C) to upper amphibolite facies (partial melting, >700 °C) conditions. The results suggest that the metapelite experienced a cryptic early metamorphic stage at c. 38 Ma at <550 °C, ≥0.85 GPa during which plagioclase was probably absent. This was followed by a prolonged high‐T, medium‐pressure (~600 °C, 0.55 GPa) evolution at 35–29 Ma during which the garnet grew, and subsequent partial melting at >690 °C and >18 Ma. Our data confirm that both geochronometers can crystallize independently at different times along the same P–T path and that neither monazite nor zircon necessarily provides timing constraints on ‘peak’ metamorphism. Therefore, collecting monazite and zircon ages as well as major and trace‐element data from major and accessory phases in the same sample is essential for reconstructing the most coherent metamorphic P–T–t evolution and thus for robustly constraining the rates and timescales of metamorphic cycles.     

Modelling groundwater/surface water interaction in a managed riparian chalk valley wetland

Understanding hydrological processes in wetlands may be complicated by management practices and complex groundwater/surface water interactions. This is especially true for wetlands underlain by permeable geology, such as chalk. In this study, the physically based, distributed model MIKE SHE is used to simulate hydrological processes at the Centre for Ecology and Hydrology River Lambourn Observatory, Boxford, Berkshire, UK. This comprises a 10‐ha lowland, chalk valley bottom, riparian wetland designated for its conservation value and scientific interest. Channel management and a compound geology exert important, but to date not completely understood, influences upon hydrological conditions. Model calibration and validation were based upon comparisons of observed and simulated groundwater heads and channel stages over an equally split 20‐month period. Model results are generally consistent with field observations and include short‐term responses to events as well as longer‐term seasonal trends. An intrinsic difficulty in representing compressible, anisotropic soils limited otherwise excellent performance in some areas. Hydrological processes in the wetland are dominated by the interaction between groundwater and surface water. Channel stage provides head boundaries for broad water levels across the wetland, whilst areas of groundwater upwelling control discrete head elevations. A relic surface drainage network confines flooding extents and routes seepage to the main channels. In‐channel macrophyte growth and its management have an acute effect on water levels and the proportional contribution of groundwater and surface water. The implications of model results for management of conservation species and their associated habitats are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.