Common-reflection-surface stack — a real data example

The simulation of a zero-offset (ZO) stack section from multi-coverage reflection data is a standard imaging method in seismic processing. It significantly reduces the amount of data and increases the signal-to-noise ratio due to constructive interference of correlated events. Conventional imaging methods, e.g., normal moveout (NMO)/dip moveout (DMO)/stack or pre-stack migration, require a sufficiently accurate macro-velocity model to yield appropriate results, whereas the recently introduced common-reflection-surface stack does not depend on a macro-velocity model. For two-dimensional seismic acquisition, its stacking operator depends on three wavefield attributes and approximates the kinematic multi-coverage reflection response of curved interfaces in laterally inhomogeneous media. The common-reflection-surface stack moveout formula defines a stacking surface for each particular sample in the ZO section to be simulated. The stacking surfaces that fit best to actual events in the multi-coverage data set are determined by means of coherency analysis. In this way, we obtain a coherency section and a section of each of the three wavefield attributes defining the stacking operator. These wavefield attributes characterize the curved interfaces and, thus, can be used for a subsequent inversion. In this paper, we focus on an application to a real land data set acquired over a salt dome. We propose three separate one-parametric search and coherency analyses to determine initial common-reflection-surface stack parameters. Optionally, a subsequent optimization algorithm can be performed to refine these initial parameters. The simulated ZO section obtained by the common-reflection-surface stack is compared to the result of a conventional NMO/DMO/stack processing sequence. We observe an increased signal-to-noise ratio and an improved continuity along the events for our proposed method — without loss of lateral resolution.     

An ArcGIS® Approach to Include Tectonic Structures in Point Data Regionalization

Point data derived from drilling logs must often be regionalized. However, aquifers may show discontinuous surface structures, such as the offset of an aquitard caused by tectonic faults. One main challenge has been to incorporate these structures into the regionalization process of point data. We combined ordinary kriging and inverse distance weighted (IDW) interpolation to account for neotectonic structures in the regionalization process. The study area chosen to test this approach is the largest porous aquifer in Austria. It consists of three basins formed by neotectonic events and delimited by steep faults with a vertical offset of the aquitard up to 70 m within very short distances. First, ordinary kriging was used to incorporate the characteristic spatial variability of the aquitard location by means of a variogram. The tectonic faults could be included into the regionalization process by using breaklines with buffer zones. All data points inside the buffer were deleted. Last, IDW was performed, resulting in an aquitard map representing the discontinuous surface structures. This approach enables one to account for such surfaces using the standard software package ArcGIS^®; therefore, it could be adopted in many practical applications.     

Silica diagenesis in Cenozoic mudstones of the North Viking Graben: physical properties and basin modelling

Silica diagenesis can significantly change physical properties of the host strata and release large volumes of water. Predicting these changes and their timing is essential to understanding compaction, fluid flow and rock deformation in sedimentary basins. In this paper, the influence of silica diagenesis (opal‐A/CT transformation) on physical properties is determined, the sediment volume affected by these changes is mapped, and a new technique to model silica diagenesis is introduced. A petrophysical analysis of 16 exploration wells shows that the opal‐A/CT transformation leads to a porosity reduction of c.20% (from 49 to 29%) in Cenozoic mudstones of the North Viking Graben. Using three‐dimensional seismic reflection data, it is shown that the c.50 m thick opal‐A/CT transformation zone covers an area of >1500 km², equating to a minimum volume of 75 km³. The spatial and temporal evolution of opal‐A/CT transformation is simulated using an innovative basin modelling approach, the results of which indicate that the transformation started around Middle‐to‐Late Eocene times and then migrated upwards until it gradually fossilised between the Miocene and present. These findings are important, as they help understanding how these sediments compact and when fluids are released by diagenesis.     

Inventory and evolution of glacial lakes since the Little Ice Age: Lessons from the case of Switzerland

Retreating glaciers give way to new landscapes with lakes as an important element. In this study, we combined available data on lake outlines with historical orthoimagery and glacier outlines for six time periods since the end of the Little Ice Age (LIA; ~1850). We generated a glacial lake inventory for modern times (2016) and traced the evolution of glacial lakes that formed in the deglaciated area since the LIA. In this deglaciated area, a total of 1192 lakes formed over the period of almost 170 years, 987 of them still in existence in 2016. Their total water surface in 2016 was 6.22 ± 0.25 km². The largest lakes are > 0.4 km² (40 ha) in size, while the majority (> 90%) are smaller than 0.01 km². Annual increase rates in area and number peaked in 1946–1973, decreased towards the end of the 20th century, and reached a new high in the latest period 2006–2016. For a period of 43 years (1973–2016), we compared modelled overdeepenings from previous studies to actual lake genesis. For a better prioritization of formation probability, we included glacier-morphological criteria such as glacier width and visible crevassing. About 40% of the modelled overdeepened area actually got covered by lakes. The inclusion of morphological aspects clearly aided in defining a lake formation probability to be linked to each modelled overdeepening. Additional morphological variables, namely dam material and type, surface runoff, and freeboard, were compiled for a subset of larger and ice-contact lakes in 2016, constituting a basis for future hazard assessment.     

Catchments as reactors: a comprehensive approach for water fluxes and solute turnover

Sustainable water quality management requires a profound understanding of water fluxes (precipitation, run-off, recharge, etc.) and solute turnover such as retention, reaction, transformation, etc. at the catchment or landscape scale. The Water and Earth System Science competence cluster (WESS, http://www.wess.info/) aims at a holistic analysis of the water cycle coupled to reactive solute transport, including soil–plant–atmosphere and groundwater–surface water interactions. To facilitate exploring the impact of land-use and climate changes on water cycling and water quality, special emphasis is placed on feedbacks between the atmosphere, the land surface, and the subsurface. A major challenge lies in bridging the scales in monitoring and modeling of surface/subsurface versus atmospheric processes. The field work follows the approach of contrasting catchments, i.e. neighboring watersheds with different land use or similar watersheds with different climate. This paper introduces the featured catchments and explains methodologies of WESS by selected examples.     

Exploring spiral galaxy potentials with hydrodynamical simulations

We study how well the complex gas velocity fields induced by massive spiral arms are modelled by the hydrodynamical simulations that we used recently to constrain the dark matter fraction in nearby spiral galaxies. More specifically, we explore the dependence of the positions and amplitudes of features in the gas flow on the temperature of the interstellar medium (assumed to behave as a one-component isothermal fluid), the non-axisymmetric disc contribution to the galactic potential, the pattern speed  Ω_p  , and finally the numerical resolution of the simulation. We argue that, after constraining the pattern speed reasonably well by matching the simulations to the observed spiral arm morphology, the amplitude of the non-axisymmetric perturbation (the disc fraction) is left as the primary parameter determining the gas dynamics. However, owing to the sensitivity of the positions of the shocks to modelling parameters, one has to be cautious when quantitatively comparing the simulations to observations. In particular, we show that a global least-squares analysis is not the optimal method for distinguishing different models, as it tends to slightly favour low disc fraction models. Nevertheless, we conclude that, given observational data of reasonably high spatial resolution and an accurate shock-resolving hydro-code, this method tightly constrains the dark matter content within spiral galaxies. We further argue that, even if the perturbations induced by spiral arms are weaker than those of strong bars, they are better suited for this kind of analysis because the spiral arms extend to larger radii where effects like inflows due to numerical viscosity and morphological dependence on gas sound speed are less of a concern than they are in the centres of discs.     

P-T-X data on P21/c-clinopyroxenes and their displacive phase transitions

The P2₁/c clinopyroxene kanoite (ideally MnMgSi₂O₆) was studied as a function of pressure and temperature using powder X-ray diffraction, differential scanning calorimetry (DSC) and optical methods. The temperature of the P2₁/c to high-temperature (HT) C2/c transition ranges from 425?°C in endmember MnMgSi₂O₆ to 125?°C in natural samples with an aegirine component. Compiling pigeonite and clinoenstatite–clinoferrosilite literature data, the temperature of the transformation was found to decrease linearly with M2 cation size. A synchrotron powder diffraction study in a heated diamond-anvil cell (DAC) yielded compression and thermal expansion data for low kanoite of composition Mn_1.2Mg_0.4Fe_0.4Si₂O₆. The high-pressure (HP) phase transition from P2₁/c to HPC2/c was reversed at 5.8 GPa at 417?°C. The high-temperature phase transition from P2₁/c to HTC2/c was rather indistinct and occurred at approximately 530?°C and 0.76 GPa. In a separate experiment, the HT transition was observed optically in a hydrothermal DAC between 0.0 and 0.4 GPa. The in-situP-T data of both experiments yielded an increase in transition temperature with increasing pressure (approx. 149?°C/GPa) and suggest a change in character of the transition from first order to continuous with increasing pressure. The data indicate that the HTC2/c and HPC2/c polymorphs are distinct phases with different stability fields. Since the high-temperature and the high-pressure polymorphs of kanoite were shown to be isotypic with other low-Ca clinopyroxenes such as the (Mg,Fe)SiO₃ series, the conclusions we draw from this study are valid for all clinopyroxenes with small (<0.88 Å) M1 and M2 cation sizes. The petrologic implications of these conclusions for the occurrence of “clinoenstatite” in the Alpe Arami peridotite are discussed.     

Phosphorus Cycling and Burial in Sediments of a Seasonally Hypoxic Marine Basin

Recycling of phosphorus (P) from sediments contributes to the development of bottom-water hypoxia in many coastal systems. Here, we present results of a year-long assessment of P dynamics in sediments of a seasonally hypoxic coastal marine basin (Lake Grevelingen, the Netherlands) in 2012. Sequential phosphorus extractions (SEDEX) and X-ray absorption spectroscopy (XAS) indicate that P was adsorbed to Fe-(III)-(oxyhydr)oxides when cable bacteria were active in the surface sediments in spring. With the onset of summer hypoxia, sulphide-induced dissolution of the Fe-(III)-(oxyhydr)oxides led to P release to the pore water and overlying water. The similarity in authigenic Ca-P concentrations in the sediment and suspended matter suggest that Ca-P is not formed in situ. The P burial efficiency was ≤ 32%. Hypoxia-driven sedimentary P recycling had a major impact on the water-column chemistry in the basin in 2012. Water-column monitoring data indicate up to ninefold higher surface water concentrations of phosphate in the basin in the late 1970s and a stronger hypoxia-driven seasonal P release from the sediment. The amplified release of P from the sediment in the past is attributed to the presence of a larger pool of Fe-bound P in the basin prior to the first onset of hypoxia. Given that P is not limiting, primary production in the basin has not been affected by the decadal changes in P availability and recycling over the past 40 years. The changes in P dynamics on decadal time scales were not recorded in sediment profiles of total P or organic C/total P.     

High precision triple oxygen isotope composition of small size urban micrometeorites indicating constant influx composition in the early geologic past

In this study, we present a method for high precision Δ′¹⁷O (Δ′¹⁷O_RL = ln(δ¹⁷O + 1) – λ_RL ln(δ¹⁸O + 1)) analysis of small mass silicate and oxide materials. The analyses were conducted by laser fluorination in combination with gas chromatography and continuous flow isotope ratio monitoring gas spectrometry. We could analyze the oxygen isotope composition of samples down to 1 μg, which corresponded to about 13 nmol O₂. The analytical error (we report the 1σ external reproducibility of a single analysis) in δ¹⁸O increases with decreasing sample sizes from ~0.2‰ for ~20 μg samples to ~0.9‰ for 1 μg samples. For Δ′¹⁷O, we achieved an external reproducibility of 0.04‰ for a sample mass range between 1 and 27 μg. The uncertainty in Δ′¹⁷O is smaller than the uncertainty in δ¹⁸O due to the correlated errors in δ¹⁷O and δ¹⁸O. We applied the method to urban micrometeorites, that is, small meteorites (<2 mm) that were sampled from a rooftop in Berlin, Germany. A total of 10 melted micrometeorites (S-type cosmic spherules, masses between 11 and 22 μg) were analyzed. The oxygen isotope compositions are comparable to that of modern Antarctic collections, indicating that the urban micrometeorites sample the same population. No indication for terrestrial weathering had been identified in the studied set of urban micrometeorites making them suitable materials for the study of micrometeorite origins.