Enhancement of prestack diffraction data and attributes using a traveltime decomposition approach

Diffractions not only carry important information about small-scale subsurface structures, they also possess unique properties, which make them a powerful tool for seismic processing and imaging. Since a point diffractor scatters an incoming wave to all directions, a diffraction event implies better illumination than a reflection, because the rays travel through larger parts of the subsurface. Furthermore, unlike the reflection case, in which the emergence location of the reflected wave depends on the source position, in the case of non-Snell scattering, up-going and down-going raypaths are decoupled. Based on this decoupling, we introduce a diffraction traveltime decomposition principle, which establishes a direct connection between zero-offset and finite-offset diffraction wavefield attributes. By making use of this approach, we are able to enhance diffractions and obtain high-quality diffraction wavefield attributes at arbitrary offsets in the prestack domain solely based on zero-offset processing without any further optimization of attributes. We show the accuracy of the method by fitting diffraction traveltimes, and on simple waveform data. Application to complex synthetic data shows the ability of the proposed approach to enhance diffractions and provide high-quality wavefield attributes even in sparsely illuminated regions such as subsalt areas. The promising results reveal a high potential for improved prestack data enhancement and further applications such as efficient diffraction-based finite-offset tomography.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

Groundwater flow and contaminant transport modelling at an air weapons range

Numerical modelling was done at the Cold Lake Air Weapons Range, Canada, to test whether the dissolved RDX and nitrate detected in groundwater come from the same sources, and to predict whether contamination poses a threat to the surface water receptors near the site. Military live fire training activities may indeed pose a risk of contamination to groundwater resources, however field investigations on military bases are quite recent, and little information is available on the long-term behaviour of munition residues related contaminants. Very limited information was available about the contaminant source zones, which were assigned based on our knowledge of current training activities. The RDX plume was well represented with the model, but the heterogeneous distribution of nitrate concentrations was more difficult to reproduce. It was nonetheless determined that both contaminants originate from the same areas. According to the model, both contaminants should reach the nearby river, but concentrations in the river should remain very low if the source zone concentration does not change. Finally, the model allowed the recommendation of a new location for the main bombing target, which would offer added protection to the river and the lake into which it flows.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Applications of particle-tracking techniques to bank infiltration: a case study from El Paso,Texas, USA

This paper presents results of a small scale study that utilized particle-tracking techniques to evaluate transport of river water through an alluvial aquifer in a bank infiltration testing site in El Paso, Texas, USA. The particle-tracking survey was used to better define filtration parameters. Several simulations were generated to allow visualization of the effects of well placement and pumping rate on flow paths, travel time, the size of the pumping influence zone, and proportion of river-derived water and groundwater mixing in the pumping well. Simulations indicate that migration of river water into the aquifer is generally slow. Most water does not arrive at the well by the end of an 18-day pumping period at 0.54 m³/min pumping rate for a well located 18 m from the river. Forty-four percent of the water pumped from the well was river water. The models provided important information needed to design appropriate sampling schedules for bank filtration practices and ensured meeting adequate soil-retention times. The pumping rate has more effect on river water travel time than the location of the pumping well from the river. The examples presented in this paper indicate that operating the pumping well at a doubled distance from the river increased the time required for the water to travel to the well, but did not greatly change the capture zone.     

Predictability of large interannual Arctic sea-ice anomalies

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century.     

Scaling of oblique impacts in frictional targets: Implications for crater size and formation mechanisms

Almost every meteorite impact occurs at an oblique angle of incidence, yet the effect of impact angle on crater size or formation mechanism is only poorly understood. This is, in large part, due to the difficulty of inferring impactor properties, such as size, velocity and trajectory, from observations of natural craters, and the expense and complexity of simulating oblique impacts using numerical models. Laboratory oblique impact experiments and previous numerical models have shown that the portion of the projectile’s kinetic energy that is involved in crater excavation decreases significantly with impact angle. However, a thorough quantification of planetary-scale oblique impact cratering does not exist and the effect of impact angle on crater size is not considered by current scaling laws. To address this gap in understanding, we developed iSALE-3D, a three-dimensional multi-rheology hydrocode, which is efficient enough to perform a large number of well-resolved oblique impact simulations within a reasonable time. Here we present the results of a comprehensive numerical study containing more than 200 three-dimensional hydrocode-simulations covering a broad range of projectile sizes, impact angles and friction coefficients. We show that existing scaling laws in principle describe oblique planetary-scale impact events at angles greater than 30° measured from horizontal. The displaced mass of a crater decreases with impact angle in a sinusoidal manner. However, our results indicate that the assumption that crater size scales with the vertical component of the impact velocity does not hold for materials with a friction coefficient significantly lower than 0.7 (sand). We found that increasing coefficients of friction result in smaller craters and a formation process more controlled by impactor momentum than by energy.     

Correction to: Leveraging HPC accelerator architectures with modern techniques — hydrologic modeling on GPUs with ParFlow

Computational Geosciences - A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10065-y     

Experimental and petrological constraints on local-scale interaction of biotite-amphibole gneiss with H2O-CO2-（K, Na）Cl fluids at middle-crustal conditions： Example from the Limpopo Complex, South Africa

Reaction textures and fluid inclusions in the～2.0 Ga pyroxene-bearing dehydration zones within the Sand River biotite-hornblende orthogneisses（Central Zone of the Limpopo Complex） suggest that the formation of these zones is a result of close interplay between dehydration process along ductile shear zones triggered by H₂O-CO₂-salt fluids at 750—800℃and 5.5—6.2 kbar.partial melting,and later exsolution of residual brine and H₂O-CO₂ fluids during melt crystallization at 650—700℃.These processes caused local variations of water and alkali activity in the fluids,resulting in various mineral assemblages within the dehydration zone.The petrological observations are substantiated by experiments on the interaction of the Sand River gneiss with the H₂O-CO-2-（K,Na）Cl fluids at 750 and 800℃and 5.5 kbar.It follows that the interaction of biotite-amphibole gneiss with H₂O-CO₂-（K.Na）Cl fluids is accompanied by partial melting at 750—800℃.Orthopyroxene-bearing assemblages are characteristic for temperature 800℃and are stable in equilibrium with fluids with low salt concentrations,while salt-rich fluids produce clinopyroxene-bearing assemblages.These observations arc in good agreement with the petrological data on the dehydration zones within the Sand River orthoeneisses.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.