Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm

Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.     

Application of a simple headcut advance model for gullies

Gully erosion begins in streambanks and uplands as a consequence of adjustments in driving forces on the landscape imposed by changes in land use or climate. The deleterious effects of gullies worldwide have led to many site‐specific studies of gully form and function. In the continental United States, gully erosion in agricultural land has destroyed valuable farmland yet, prediction of gully processes remains problematic on a national scale. This research has proposed a simple method to predict gully headcut advance. When combined with SWAT hydrologic flow routines, the model predicted gully headcut advance with reasonable accuracy on a daily time step for time periods exceeding two decades. The model was tested in two distinct land resource areas of the United States with differing climate, soils, cover and drainage. The inputs for the headcut model have been kept simple as the model will be applied over large areas. Model inputs consist of headcut height, headcut resistance (based on soil erodibility and a root‐cover factor), and daily flow. The model is compared with an annual time step model used in assessment of headcut advance and appears to offer a better way to assess gully headcut advance. Copyright © 2017 John Wiley & Sons, Ltd.     

Organic carbon fluxes of a glacier surface: A case study of Foxfonna,a small Arctic glacier

Arctic glaciers are rapidly responding to global warming by releasing organic carbon (OC) to downstream ecosystems. The glacier surface is arguably the most biologically active and biodiverse glacial habitat and therefore the site of important OC transformation and storage, although rates and magnitudes are poorly constrained. In this paper, we present measurements of OC fluxes associated with atmospheric deposition, ice melt, biological growth, fluvial transport and storage (in superimposed ice and cryoconite debris) for a supraglacial catchment on Foxfonna glacier, Svalbard (Norway), across two consecutive years. We found that in general atmospheric OC input (averaging 0.63 ± 0.25 Mg a^-1 total organic carbon, i.e. TOC, and 0.40 ± 0.22 Mg a^-1 dissolved organic carbon, i.e. DOC) exceeded fluvial OC export (0.46 ± 0.04 Mg a^-1 TOC and 0.36 ± 0.03 Mg a^-1 DOC). Early in the summer, OC was mobilised in snowmelt but its release was delayed by temporary storage in superimposed ice on the glacier surface. This delayed the export of 28.5% of the TOC in runoff. Biological production in cryoconite deposits was a negligible potential source of OC to runoff, while englacial ice melt was far more important on account of the glacier's negative ice mass balance (–0.89 and –0.42 m a^-1 in 2011 and 2012, respectively). However, construction of a detailed OC budget using these fluxes shows an excess of inputs over outputs, resulting in a net retention of OC on the glacier surface at a rate that would require c. 3 years to account for the OC stored as cryoconite debris. © 2018 John Wiley & Sons, Ltd.     

Numerical investigation of alternative fracture stiffness measures and their respective scaling behaviours

We study the mechanical deformation of fractures under normal stress, via tangent and specific fracture stiffnesses, for different length scales using numerical simulations and analytical insights. First, we revisit an equivalent elastic layer model that leads to two expressions: the tangent stiffness is the sum of an “intrinsic” stiffness and the normal stress, and the specific stiffness is the tangent stiffness divided by the fracture aperture at current stress. Second, we simulate the deformation of rough fractures using a boundary element method where fracture surfaces represented by elastic asperities on an elastic half‐space follow a self‐affine distribution. A large number of statistically identical “parent” fractures are generated, from which sub‐fractures of smaller dimensions are extracted. The self‐affine distribution implies that the stress‐free fracture aperture increases with fracture length with a power law in agreement with the chosen Hurst exponent. All simulated fractures exhibit an increase in the specific stiffness with stress and an average decrease with increase in length consistent with field observations. The simulated specific and tangent stiffnesses are well described by the equivalent layer model provided the “intrinsic” stiffness slightly decreases with fracture length following a power law. By combining numerical simulations and the analytical model, the effect of scale and stress on fracture stiffness measures can be easily separated using the concept of “intrinsic” stiffness. We learn that the primary reason for the variability in specific stiffness with length comes from the fact that the typical aperture of the self‐affine fractures itself scales with the length of the fractures.     

War landform mapping and classification on the Verdun battlefield (France) using airborne LiDAR and multivariate analysis

Acting as efficient earth-movers, soldiers can be viewed as significant geomorphological drivers of landscape change when replaced in the recent debates on Anthropocene Geomorphology. ‘Polemoforms’, generated by military activities, correspond with a set of human-made landforms of various sizes and geometries. They are particularly common on the World War One battlefield of Verdun (France) which ranks among the largest battles of attrition along the Western Front. The artillery bombardments and building of defensive positions in that battle significantly altered the landscape, resulting in thousands of shell craters, dugouts, and gun positions that have altered both the meso and microtopography. This paper proposes an innovative methodology to make an exhaustive inventory of these small-scale conflict-induced landforms (excluding linear features such as trenches) using a digital terrain model (DTM) acquired by airborne LiDAR on the whole battlefield. Morphometric analysis was conducted using Kohonen's self-organizing maps (SOMs) and hierarchical agglomerative clustering (HAC) in order to quantify and classify the high number of war landforms. This combined approach allowed for mapping more than one million landforms which can be classified into eight different shapes including shell craters and various soldier-made landforms (i.e. shelters, gun positions, etc.). Detection quality evaluation using field observations revealed the algorithm successfully classified 93% of shell craters and 74% of anthropologically constructed landforms. Finally, the iconographic database and map series produced will help archaeologists and foresters to better manage the historic site of Verdun, today covered by a large forest of ~10 000 ha. © 2019 John Wiley & Sons, Ltd.     

Quantifying long-term post-fire sediment delivery and erosion mitigation effectiveness

Large wildfires can have profound and lasting impacts not only from direct consumption of vegetation but also longer-term effects such as persistent soil erosion. The 2002 Hayman Fire burned in one of the watersheds supplying water to the Denver metropolitan area; thus there was concern regarding hillslope erosion and sedimentation in the reservoirs. The efficacy of various treatments for reducing erosion was tested, including hand scarification on contour, agricultural straw mulch, wood mulch, burned controls and unburned reference plots. Simulated rill erosion experiments were used both immediately after the fire and again 10 years post fire. To better understand untreated recovery, the same experiments were applied to control plots in post-fire years 1, 2, 3 and 4, and in unburned reference plots in years 4 and 10. Results indicate that control and scarified plots produced significantly greater sediment flux rates – 1.9 and 2.8 g s⁻¹ respectively – than the straw and wood mulch treatments – 0.9 and 1.1 g s⁻¹ – immediately after the fire. Mulch treatments reduced runoff rate, runoff velocity, and sediment concentration and flux rate. The straw mulch cover was no longer present, whereas the wood mulch was still there in year 10. Vegetation regrowth was slow and mulch treatments provided effective cover to reduce sediment right after the fire. In post-fire year 10, there were no significant differences in sediment flux rates across treatments; it is notable, however, that the wood mulch treatment (0.09 g s⁻¹) most closely approached the unburned condition (0.07 g s⁻¹). The burned control plots had high sediment flux rates until post-fire year 3, when flux rates significantly decreased and were statistically no longer higher than the unburned levels from year 4 and 10. These results will inform managers of the longer-term post-fire sediment delivery rates and of the ability of post-fire emergency hillslope treatments to mitigate erosion rates. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.     

Geomorphic process signatures reshaping sub-humid Mediterranean badlands: 2. Application to 5-year dataset

Badland landscapes exhibit high erosion rates and represent the main source of fine sediments in some catchments. Advances in high-resolution topographic methods allow analysis of topographic changes at high temporal and spatial scales. We apply the Mapping Geomorphic Processes in the Environment (MaGPiE) algorithm to infer the main geomorphic process signatures operating in two sub-humid badlands with contrasting morphometric attributes located in the Southern Pyrenees. By interrogating a 5-year dataset of seasonal and annual topographic changes, we examine the variability of geomorphic processes at multiple temporal scales. The magnitude of geomorphic processes is linked to landform attributes and meteorological variables. Morphometric differences between both adjacent badlands allow us to analyse the role of landform attributes in the main geomorphic process reshaping landscapes subjected to the same external forcing (i.e. rainfall and temperature). The dominant geomorphic process signatures observed in both badlands are different, despite their close proximity and the same rainfall and temperature regimes. Process signatures determining surface lowering in the gently sloping south-facing badland, characterized by lower connectivity and more vegetation cover, are driven by surface runoff-based processes, both diffuse (causing sheet washing) and concentrated (determining cutting and filling, rilling and gullying). The steeper, more connected north-facing slopes of the other badland are reshaped by means of gravitational processes, with mass wasting dominating topographic changes. In terms of processes determining surface raising, both mass wasting and cutting and filling are most frequently observed in both badlands. There is a clear near-balanced feedback between both surface-raising and -lowering processes that becomes unbalanced at larger temporal scales due to the thresholds overcome, as the volume associated with surface lowering becomes higher than that associated with raising-based processes. Rainfall variables control surface flow processes, while those variables associated with low temperature have a significant relation with mass movement-based processes and other localized processes such as regolith cohesion loss. Finally, our results point out that morphometry (slope and connectivity) together with vegetation cover are key factors determining geomorphic processes and associated topographic changes. © 2020 John Wiley & Sons, Ltd.     

From marine bands to hybrid flows: Sedimentology of a Mississippian black shale

Organic-rich mudstones have long been of interest as conventional and unconventional source rocks and are an important organic carbon sink. Yet the processes that deposited organic-rich muds in epicontinental seaways are poorly understood, partly because few modern analogues exist. This study investigates the processes that transported and deposited sediment and organic matter through part of the Bowland Shale Formation, from the Mississippian Rheic–Tethys seaway. Field to micron-scale sedimentological analysis reveals a heterogeneous succession of carbonate-rich, siliceous, and siliciclastic, argillaceous muds. Deposition of these facies at basinal and slope locations was moderated by progradation of the nearby Pendle delta system, fourth-order eustatic sea-level fluctuation and localized block and basin tectonism. Marine transgressions deposited bioclastic ‘marine band’ (hemi)pelagic packages. These include abundant euhaline macrofaunal tests, and phosphatic concretions of organic matter and radiolarian tests interpreted as faecal pellets sourced from a productive water column. Lens-rich (lenticular) mudstones, hybrid, debrite and turbidite beds successively overlie marine band packages and suggest reducing basin accommodation promoted sediment deposition via laminar and hybrid flows sourced from the basin margins. Mud lenses in lenticular mudstones lack organic linings and bioclasts and are equant in early-cemented lenses and in plan-view, and are largest and most abundant in mudstones overlying marine band packages. Thus, lenses likely represent partially consolidated mud clasts that were scoured and transported in bedload from the shelf or proximal slope, as a ‘shelf to basin’ conveyor, during periods of reduced basin accommodation. Candidate in situ microbial mats in strongly lenticular mudstones, and as rip-up fragments in the down-dip hybrid beds, suggest that these were potentially key biostabilizers of mud. Deltaic mud export was fast, despite the intrabasinal complexity, likely an order of magnitude higher than similar successions deposited in North America. Epicontinental basins remotely linked to delta systems were therefore capable of rapidly accumulating both sediment and organic matter.     

A comparison of biological trends from four marine ecosystems: Synchronies, differences, and commonalities

Major features of four marine ecosystems were analyzed based on a broad range of fisheries-associated datasets and a suite of oceanographic surveys. The ecosystems analyzed included the Gulf of Maine/Georges Bank in the Northwest Atlantic Ocean, the Norwegian/Barents Seas in the Northeast Atlantic Ocean, and the eastern Bering Sea and the Gulf of Alaska in the Northeast Pacific Ocean. We examined survey trends in major fish abundances, total system fish biomass, and zooplankton biomasses. We standardized each time series and examined trends and anomalies over time, using both time series and cross-correlational statistical methods. We compared dynamics of functionally analogous species from each of these four ecosystems. Major commonalities among ecosystems included a relatively stable amount of total fish biomass and the importance of large calanoid copepods, small pelagic fishes and gadids. Some of the changes in these components were synchronous across ecosystems. Major differences between ecosystems included gradients in the magnitude of total fish biomass, commercial fish biomass, and the timing of major detected events. This work demonstrates the value of comparative analysis across a wide range of marine ecosystems, suggestive of very few but none-the-less detectable common features across all northern hemisphere ocean systems.     

The release of dissolved actinium to the ocean: A global comparison of different end-members

The measurement of short-lived ²²³Ra often involves a second measurement for supported activities, which represents ²²⁷Ac in the sample. Here we exploit this fact, presenting a set of 284 values on the oceanic distribution of ²²⁷Ac, which was collected when analyzing water samples for short-lived radium isotopes by the radium delayed coincidence counting system. The present work compiles ²²⁷Ac data from coastal regions all over the northern hemisphere, including values from ground water, from estuaries and lagoons, and from marine end-members. Deep-sea samples from a continental slope off Puerto Rico and from an active vent site near Hawaii complete the overview of ²²⁷Ac near its potential sources.The average ²²⁷Ac activities of nearshore marine end-members range from 0.4 dpm m^− 3 at the Gulf of Mexico to 3.0 dpm m^− 3 in the coastal waters of the Korean Strait. In analogy to ²²⁸Ra, we find the extension of adjacent shelf regions to play a substantial role for ²²⁷Ac activities, although less pronounced than for radium, due to its weaker shelf source. Based on previously published values, we calculate an open ocean ²²⁷Ac inventory of 1.35 * 10¹⁸ dpm ²²⁷Ac_ex in the ocean, which corresponds to 37 moles, or 8.4 kg. This implies a flux of 127 dpm m⁻² y^− 1 from the deep-sea floor. For the shelf regions, we obtain a global inventory of ²²⁷Ac of 4.5 * 10¹⁵ dpm, which cannot be converted directly into a flux value, as the regional loss term of ²²⁷Ac to the open ocean would have to be included.Ac has so far been considered to behave similarly to Ra in the marine environment, with the exception of a strong Ac source in the deep-sea due to ²³¹Pa_ex. Here, we present evidence of geochemical differences between Ac, which is retained in a warm vent system, and Ra, which is readily released [Moore, W.S., Ussler, W. and Paull, C.K., 2008-this issue. Short-lived radium isotopes in the Hawaiian margin: Evidence for large fluid fluxes through the Puna Ridge. Marine Chemistry]. Another potential mechanism of producing deviations in ²²⁷Ac/²²⁸Ra and daughter isotope ratios from the expected production value of lithogenic material is observed at reducing environments, where enrichment in uranium may occur. The presented data here may serve as a reference for including ²²⁷Ac in circulation models, and the overview provides values for some end-members that contribute to the global Ac distribution.