Pre-stack full waveform inversion of ultra-high-frequency marine seismic reflection data

The full waveform inversion (FWI) of seismic reflection data aims to reconstruct a detailed physical properties model of the subsurface, fitting both the amplitude and the traveltime of the reflections generated at physical discontinuities in the propagation medium. Unlike reservoir-scale seismic exploration, where seismic inversion is a widely adopted remote characterization tool, ultrahigh-frequency (UHF, 0.2–4.0 kHz) multichannel marine reflection seismology is still most often limited to a qualitative interpretation of the reflections’ architecture. Here we propose an elastic FWI methodology, custom-tailored for pre-stack UHF marine data in vertically heterogeneous media to obtain a decimetric-scale distribution of P-impedance, density and Poisson’s ratio within the shallow subseabed sediments. We address the deterministic multiparameter inversion in a sequential fashion. The complex trace instantaneous phase is first inverted for the P-wave velocity to make up for the lack of low frequency in the data and reduce the nonlinearity of the problem. This is followed by a short-offset P-impedance optimization and a further step of full offset range Poisson’s ratio inversion. Provided that the seismogram contains wide reflection angles (>40°), we show that it is possible to invert for density and decompose a posteriori the relative contribution of P-wave velocity and density to the P-impedance. A broad range of synthetic tests is used to prove the potential of the methodology and highlights sensitivity issues specific to UHF seismic. An example application to real data is also presented. In the real case, trace normalization is applied to minimize the systematic error deriving from an inaccurate source wavelet estimation. The inverted model for the top 15 m of the subseabed agrees with the local lithological information and core-log data. Thus, we can obtain a detailed remote characterization of the shallow sediments using a multichannel sub-bottom profiler within a reasonable computing cost and with minimal pre-processing. This has the potential to reduce the need of extensive geotechnical coring campaigns

Decimetric-resolution stochastic inversion of shallow marine seismic reflection data: dedicated strategy and application to a geohazard case study

Characterization of the top 10–50 m of the subseabed is key for landslide hazard assessment, offshore structure engineering design and underground gas-storage monitoring. In this paper, we present a methodology for the stochastic inversion of ultra-high-frequency (UHF, 0.2–4.0 kHz) pre-stack seismic reflection waveforms, designed to obtain a decimetric-resolution remote elastic characterization of the shallow sediments with minimal pre-processing and little a priori information. We use a genetic algorithm in which the space of possible solutions is sampled by explicitly decoupling the short and long wavelengths of the P-wave velocity model. This approach, combined with an objective function robust to cycle skipping, outperforms a conventional model parametrization when the ground-truth is offset from the centre of the search domain. The robust P-wave velocity model is used to precondition the width of the search range of the multiparameter elastic inversion, thereby improving the efficiency in high-dimensional parametrizations. Multiple independent runs provide a set of independent results from which the reproducibility of the solution can be estimated. In a real data set acquired in Finneidfjord, Norway, we also demonstrate the sensitivity of UHF seismic inversion to shallow subseabed anomalies that play a role in submarine slope stability. Thus, the methodology has the potential to become an important practical tool for marine ground model building in spatially heterogeneous areas, reducing the reliance on expensive and time-consuming coring campaigns for geohazard mitigation in marine areas

Gas migration pathways, controlling mechanisms and changes in sediment acoustic properties observed in a controlled sub-seabed CO2 release experiment

Carbon capture and storage (CCS) is a key technology to potentially mitigate global warming by reducing carbon dioxide (CO2) emissions from industrial facilities and power generation that escape into the atmosphere. To broaden the usage of geological storage as a viable climate mitigation option, it is vital to understand CO2 behaviour after its injection within a storage reservoir, including its potential migration through overlying sediments, as well as biogeochemical and ecological impacts in the event of leakage. The impacts of a CO2 release were investigated by a controlled release experiment that injected CO2 at a known flux into shallow, under-consolidated marine sediments for 37 days. Repeated high-resolution 2D seismic reflection surveying, both pre-release and syn-release, allows the detection of CO2-related anomalies, including: seismic chimneys; enhanced reflectors within the subsurface; and bubbles within the water column. In addition, reflection coefficient and seismic attenuation values calculated for each repeat survey, allow the impact of CO2 flux on sediment acoustic properties to be comparatively monitored throughout the gas release. CO2 migration is interpreted as being predominantly controlled by sediment stratigraphy in the early stages of the experiment. However, either the increasing flow rate, or the total injected volume become the dominant factors determining CO2 migration later in the experiment

3D seismic imaging of buried Younger Dryas mass movement flows: Lake Windermere, UK

Windermere is a glacially overdeepened lake located in the southeastern Lake District, UK. Using the threedimensional(3D) Chirp subbottom profiler, we image mass movement deposits related to the Younger Dryas(YD) within a decimetre-resolution 3D seismic volume, documenting their internal structure and interactionwith preexisting deposits in unprecedented detail. Three distinct flow events are identified and mappedthroughout the 3D survey area. Package structures and seismic attributes classify them as: a small (totalvolume of c. 1500 m3) debris flow containing deformed translated blocks; a large (inferred total volume ofc. 500,000 m3), homogeneous fine-grained mass flow deposit; and a debris flow (inferred total volume ofc. 60,000 m3) containing small (c. 8.0×2.0 m) deformed translated blocks. Geomorphological mapping oftheir distribution and interaction with preexisting sediments permit the reconstruction of a depositionalhistory for the stratigraphic units identified in the seismic volume.<br/

A novel sub-seabed CO\u3csub\u3e2\u3c/sub\u3e release experiment informing monitoring and impact assessment for geological carbon storage

© 2014 The Authors. Carbon capture and storage is a mitigation strategy that can be used to aid the reduction of anthropogenic CO2 emissions. This process aims to capture CO2 from large point-source emitters and transport it to a long-term storage site. For much of Europe, these deep storage sites are anticipated to be sited below the sea bed on continental shelves. A key operational requirement is an understanding of best practice of monitoring for potential leakage and of the environmental impact that could result from a diffusive leak from a storage complex. Here we describe a controlled CO2 release experiment beneath the seabed, which overcomes the limitations of laboratory simulations and natural analogues. The complex processes involved in setting up the experimental facility and ensuring its successful operation are discussed, including site selection, permissions, communications and facility construction. The experimental design and observational strategy are reviewed with respect to scientific outcomes along with lessons learnt in order to facilitate any similar future

A new varve sequence from Windermere, UK, records rapid ice retreat prior to the Lateglacial Interstadial (GI-1)

Annually laminated sediments (varves) provide excellent temporal resolution to study rapid environmental change, but are rare in the early part of the Last Termination (∼19–∼11.7 ka BP). We present a new >400 varve year (vyr) varve sequence in two floating parts from Windermere, a lake at the southern margin of the mountains of northwest England. This sequence records the final retreat of the Windermere glacier at the southern edge of the Lake District Ice Cap during the transition from Heinrich Stadial 1 (∼18–∼14.7 ka BP) into the Lateglacial Interstadial (∼14.7–∼12.9 ka BP). Laminated sediments from four lake cores from Windermere's northern and southern basins were investigated and shown to be varved. These sequences are integrated with seismic reflection evidence to reconstruct south-to-north deglaciation. Seismic and sedimentological evidence is consistent with gradual stepped ice retreat along the entire southern basin and into the northern basin between 255 and 700 vyr prior to the appearance of significant biota in the sediment that heralded the Lateglacial Interstadial, and had retreated past a recessional moraine (RM8) in the northern basin by 121 vyr prior to the interstadial. The Lateglacial Interstadial age of this biota-bearing unit was confirmed by 14C-dating, including one date from the northernmost core of ∼13.5 cal ka BP. A change in mineralogy in all four cores as the glacier retreated north of the Dent Group (the northernmost source of calcareous bedrock) and a decrease in coarse grains in the varves shows that the ice had retreated along the entire North Basin at ∼70 vyr prior to the Lateglacial Interstadial. The estimated retreat rate is 70–114 m yr−1 although buried De Geer moraines, if annual, may indicate retreat of 120 m yr−1 with a ≥3 year stillstand at a recessional moraine halfway along the basin. The glacier then retreated north of the lake basin, becoming land-terminating and retreating at 92.5–49 m yr−1. The northernmost core has a varve sequence ending at least 111 vyr after the other core chronologies, due to the increased proximity to remnant ice in the catchment uplands into the early Lateglacial Interstadial. We show that almost all of the glacier retreat in the Windermere catchment occurred before the abrupt warming at the onset of the Lateglacial Interstadial, in keeping with similar findings from around the Irish Sea Basin, and suggesting a similar retreat timescale for other radial valley glaciers of the Lake District Ice Cap. The seismic and core evidence also show the potential for a much longer varve chronology extending at least 400 and potentially over 1000 vyr further back into Heinrich Stadial 1 (18–14.7 ka BP), suggesting that glacier retreat in the Windermere valley initiated at least before 15.5 ka BP and perhaps 16 ka BP

The planform mobility of river channel confluences: Insights from analysis of remotely sensed imagery

River channel confluences are widely acknowledged as important geomorphological nodes that control the downstream routing of water and sediment, and which are locations for the preservation of thick fluvial deposits overlying a basal scour. Despite their importance, there has been little study of the stratigraphic characteristics of river junctions, or the role of confluence morphodynamics in influencing stratigraphic character and preservation potential. As a result, although it is known that confluences can migrate through time, models of confluence geomorphology and sedimentology are usually presented from the perspective that the confluence remains at a fixed location. This is problematic for a number of reasons, not least of which is the continuing debate over whether it is possible to discriminate between scour that has been generated by autocyclic processes (such as confluence scour) and that driven by allocyclic controls (such as sea-level change). This paper investigates the spatial mobility of river confluences by using the 40-year record of Landsat Imagery to elucidate the styles, rates of change and areal extent over which large river confluence scours may migrate. On the basis of these observations, a new classification of the types of confluence scour is proposed and applied to the Amazon and Ganges-Brahmaputra-Meghna (GBM) basins. This analysis demonstrates that the drivers of confluence mobility are broadly the same as those that drive channel change more generally. Thus in the GBM basin, a high sediment supply, large variability in monsoonal driven discharge and easily erodible bank materials result in a catchment where over 80% of large confluences are mobile over this 40-year window; conversely this figure is < 40% for the Amazon basin. These results highlight that: i) the potential areal extent of confluence scours is much greater than previously assumed, with the location of some confluences on the Jamuna (Brahmaputra) River migrating over a distance of 20 times the tributary channel width; ii) extensive migration in the confluence location is more common than currently assumed, and iii) confluence mobility is often tied to the lithological and hydrological characteristics of the drainage basins that determine sediment yield

The sedimentology of river confluences

Channel confluences are key nodes within large river networks, and yet surprisingly little is known about their spatial and temporal evolution. Moreover, because confluences are associated with vertical scour that typically extends to several times the mean channel depth, the deposits associated with such scours should have a high preservation potential within the rock record. Paradoxically, such scours are rarely observed, and their preservation and sedimentological interpretation are poorly understood. The present study details results from a physically‐based morphodynamic model that is applied to simulate the evolution and alluvial architecture of large river junctions. Boundary conditions within the model were defined to approximate the junction of the Ganges and Jamuna rivers, Bangladesh, with the model output being supplemented by geophysical datasets collected at this junction. The numerical simulations reveal several distinct styles of sedimentary fill that are related to the morphodynamic behaviour of bars, confluence scour downstream of braid bars, bend scour and major junction scour. Comparison with existing, largely qualitative, conceptual models reveals that none of these can be applied simply, although elements of each are evident in the deposits generated by the numerical simulation and observed in the geophysical data. The characteristics of the simulated scour deposits are found to vary according to the degree of reworking caused by channel migration, a factor not considered adequately in current conceptual models of confluence sedimentology. The alluvial architecture of major junction scours is thus characterized by the prevalence of erosion surfaces in conjunction with the thickest depositional sets. Confluence scour downstream of braid bar and bend scour sites may preserve some large individual sets, but these locations are typically characterized by lower average set thickness compared to major junction scour and by a lack of large‐scale erosional surfaces. Areas of deposition not related to any of the specific scour types highlighted above record the thinnest depositional sets. This variety in the alluvial architecture of scours may go some way towards explaining the paradox of ancient junction scours, that while abundant large scours are likely in the rock record, they have been reported rarely. The present results outline the likely range of confluence sedimentology and will serve as a new tool for recognizing and interpreting these deposits in the ancient fluvial record

Whole brain radiotherapy (WBRT) after local treatment of brain metastases in melanoma patients: Statistical Analysis Plan

Background: The WBRTMel trial is a multinational, open-label, phase III randomised controlled trial comparing whole brain radiotherapy (WBRT) to observation following local treatment of one to three melanoma brain metastases with surgery and/or stereotactic irradiation. The primary trial endpoint was to determine the effect of adding WBRT to local treatment on distant intracranial control, and the secondary endpoints were neurocognitive function, quality of life (QoL), performance status, overall survival, death from intracranial causes, death from melanoma and cost-effectiveness. Objective: The objective of this update is to outline and publish the pre-determined statistical analysis plan (SAP) before the database lock and the start of analysis. Methods: The SAP describes basic analysis principles, methods for dealing with a range of commonly encountered data analysis issues and the specific statistical procedures for analysing efficacy and safety outcomes. The SAP was approved after closure of recruitment and before completion of patient follow-up. It outlines the planned primary analyses and a range of subgroup and sensitivity analyses regarding the clinical and QoL outcomes. Health economic outcomes are not included in this plan but will be analysed separately. The SAP will be adhered to for the final data analysis of this trial to avoid analysis bias arising from knowledge of the data. Results: The resulting SAP is consistent with best practice and will allow open and transparent reporting. Conclusion: We have developed a SAP for the WBRTMel trial which will be followed to ensure high-quality standards of internal validity to minimise analysis bias