Petrologic and textural diversity among the PCA 02 howardite group,one of the largest pieces of the Vestan surface

Abstract— Nine howardites and two diogenites were recovered from the Pecora Escarpment Icefield (PCA) in 2002. Cosmogenic radionuclide abundances indicate that the samples are paired and that they constituted an approximately 1 m (diameter) meteoroid prior to atmospheric entry. At about 1 m in diameter, the PCA 02 HED group represents one of the largest single pre‐atmospheric pieces of the Vestan surface yet described. Mineral and textural variations were measured in six of the PCA 02 howardites to investigate meter‐scale diversity of the Vestan surface. Mineral compositions span the range of known eucrite and diogenite compositions. Additional non‐diogenitic groups of Mg‐ and Fe‐rich olivine are observed, and are interpreted to have been formed by exogenic contamination and impact melting, respectively. These howardites contain olivine‐rich impact melts that likely formed from dunite‐ and harzburgite‐rich target rocks. Containing the first recognized olivine‐rich HED impact melts, these samples provide meteoritic evidence that olivine‐rich lithologies have been exposed on the surface of Vesta. Finally, we present a new method for mapping distributions of lithologies in howardites using 8 elemental X‐ray maps. Proportions of diogenite and eucrite vary considerably among the PCA 02 howardites, suggesting they originated from a heterogeneous portion of the Vestan surface. While whole sample modes are dominated by diogenite, the finer grain size fractions are consistently more eucritic. This discrepancy has implications for near‐infrared spectral observations of portions of Vesta’s surface that are similar to the PCA 02 howardites, as the finer grained eucritic material will disproportionately dominate the spectra.     

Approximate evaluation of the distortion to the peak ocean wavelength and direction derived from SAR image spectra

A calculation procedure using the modulation transfer function approach is presented to estimate the distortions to the values of the peak ocean wavelength and direction derived from SAR images of these waves. The distortions are calculated for ranges of SAR parameters which are typical for aircraft and satellites, and ocean wave parameters which reflect sea states which can be imaged by these SARs. The calculated distortions are discussed in terms of overall trends, effect of particular SAR and wave parameters, and qualitatively compared with observed differences between SAR and in situ measurements.     

Cross-correlation beamforming

An areal distribution of sensors can be used for estimating the direction of incoming waves through beamforming. Beamforming may be implemented as a phase-shifting and stacking of data recorded on the different sensors (i.e., conventional beamforming). Alternatively, beamforming can be applied to cross-correlations between the waveforms on the different sensors. We derive a kernel for beamforming cross-correlated data and call it cross-correlation beamforming (CCBF). We point out that CCBF has slightly better resolution and aliasing characteristics than conventional beamforming. When auto-correlations are added to CCBF, the array response functions are the same as for conventional beamforming. We show numerically that CCBF is more resilient to non-coherent noise. Furthermore, we illustrate that with CCBF individual receiver-pairs can be removed to improve mapping to the slowness domain. An additional flexibility of CCBF is that cross-correlations can be time-windowed prior to beamforming, e.g., to remove the directionality of a scattered wavefield. The observations on synthetic data are confirmed with field data from the SPITS array (Svalbard). Both when beamforming an earthquake arrival and when beamforming ambient noise, CCBF focuses more of the energy to a central beam. Overall, the main advantage of CCBF is noise suppression and its flexibility to remove station pairs that deteriorate the signal-related beampower.     

A continuous tephrostratigraphic record from the Labrador Sea spanning the last 65 ka

Volcanic ash preserved in marine sediment sequences is key for independent synchronization of palaeoclimate records within and across different climate archives. Here we present a continuous tephrostratigraphic record from the Labrador Sea, spanning the last 65–5 ka, an area and time period that has not been investigated in detail within the established North Atlantic tephra framework. We investigated marine sediment core GS16-204-22CC for increased tephra occurrences and geochemically analysed the major element composition of tephra shards to identify their source volcano(es). In total we observed eight tephra zones, of which five concentration peaks show isochronous features that can be used as independent tie-points in future studies. The main transport mechanism of tephra shards to the site was near-instantaneous deposition by drifting of sea ice along the East Greenland Current. Our results show that the Icelandic Veidivötn volcanic system was the dominant source of tephra material, especially between late Marine Isotope Stage (MIS) 4 and early MIS 3. The Veidivötn system generated volcanic eruptions in cycles of ca. 3–5 ka. We speculate that the quantity of tephra delivered to the Labrador Sea was a result of variable Icelandic ice volume and/or changes in the transportation pathway towards the Labrador Sea.     

Fluvial dynamics and 14C-10Be disequilibrium on the Bolivian Altiplano

Determining sediment transfer times is key to understanding source-to-sink dynamics and the transmission of environmental signals through the fluvial system. Previous work on the Bolivian Altiplano applied the in situ cosmogenic ¹⁴C-¹⁰Be-chronometer to river sands and proposed sediment storage times of ~10–20 kyr in four catchments southeast of Lake Titicaca. However, the fidelity of those results hinges upon isotopic steady-state within sediment supplied from the source area. With the aim of independently quantifying sediment storage times and testing the ¹⁴C-¹⁰Be steady-state assumption, we dated sediment storage units within one of the previously investigated catchments using radiocarbon dating, cosmogenic ¹⁰Be-²⁶Al isochron burial dating, and ¹⁰Be-²⁶Al depth-profile dating. Palaeosurfaces appear to preserve remnants of a former fluvial system, which has undergone drainage reversal, reduction in catchment area, and local isostatic uplift since ~2.8 Ma. From alluvium mantling the palaeosurfaces we gained a deposition age of ~580 ka, while lower down fluvial terraces yielded ≤34 ka, and floodplains ~3–1 ka. Owing to restricted channel connectivity with the terraces and palaeosurfaces, the main source of channel sediment is via reworking of the late Holocene floodplain. Yet modelling a set of feasible scenarios reveals that floodplain storage and burial depth are incompatible with the ¹⁴C-¹⁰Be disequilibrium measured in the channel. Instead we propose that the ¹⁴C-¹⁰Be offset results from: (i) non-uniform erosion whereby deep gullies supply hillslope-derived debris; and/or (ii) holocene landscape transience associated with climate or human impact. The reliability of the ¹⁴C-¹⁰Be chronometer vitally depends upon careful evaluation of sources of isotopic disequilibrium in a wide range of depositional and erosional landforms in the landscape. © 2018 John Wiley & Sons, Ltd.     

Point‐spread functions for interferometric imaging

Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.     

A Simple Analytical Formula for the Leakage Flux Through a Perforated Aquitard

In this methods note, we present a simple analytical formula to quantify the steady‐state leakage flux (Q) over a perforated aquitard. The flux depends on the aquitard thickness (D), the radius of the perforation (R), the hydraulic conductivity of the material inside the perforation (k_fill), the conductivities of the overlying and underlying aquifers (k₁ and k₂, respectively), and the head difference between the two aquifers (ΔH): This equation assumes an aquitard separating two homogeneous and infinite aquifers (R ? aquifer thickness) in which radial flow to and from the perforation occurs, with no other recharge or discharge boundaries near the perforation. The flux through a perforation in a hypothetical case study with D = 10 m, k₁ = 10 m/d, k₂ = 20 m/d, R = 0.072 m, and ΔH = 1 m ranges between less than 1 mL/d if the hole is backfilled with bentonite (k_fill = 10^?4 m/d), to several liters per day if the perforation is backfilled with sand from the overlying aquifer (k_fill = 10 m/d), to several m³/d if the perforation forms an open conduit (k_fill = 10⁵ m/d). The leakage fluxes calculated with this model agree well with those calculated using a numerical model (MODFLOW).     

Two‐dimensional controlled‐source electromagnetic interferometry by multidimensional deconvolution: spatial sampling aspects

We use numerically modelled data sets to investigate the sensitivity of electromagnetic interferometry by multidimensional deconvolution to spatial receiver sampling. Interferometry by multidimensional deconvolution retrieves the reflection response below the receivers after decomposition of the fields into upward and downward decaying fields and deconvolving the upward decaying field by the downward decaying field. Thereby the medium above the receiver level is replaced with a homogeneous half‐space, the sources are redatumed to the receiver level and the direct field is removed. Consequently, in a marine setting the retrieved reflection response is independent of any effect of the water layer and the air above. A drawback of interferometry by multidimensional deconvolution is a possibly unstable matrix inversion, which is necessary to retrieve the reflection response. Additionally, in order to correctly separate the upward and the downward decaying fields, the electromagnetic fields need to be sampled properly. We show that the largest possible receiver spacing depends on two parameters: the vertical distance between the source and the receivers and the length of the source. The receiver spacing should not exceed the larger of these two parameters. Besides these two parameters, the presence of inhomogeneities close to the receivers may also require a dense receiver sampling. We show that by using the synthetic aperture concept, an elongated source can be created from conventionally acquired data in order to overcome these strict sampling criteria. Finally, we show that interferometry may work under real‐world conditions with random noise and receiver orientation and positioning errors.     

Bifurcation modelling in a meandering gravel–sand bed river

The Rhine bifurcation at Pannerden forms the major distribution point for water supply in the Netherlands, distributing not only water and sediment but also flooding risks and navigability. Its morphological stability has been a concern for centuries. We present experiences from more than two decades of numerical morphological modelling of this bifurcation with a gravel–sand bed and a meandering planform. Successive computations have shown the importance of upstream approach conditions, the necessity to include physical mechanisms for grain sorting and alluvial roughness, and the need to assume a thicker active layer of the river bed than is suggested by laboratory flume experiments using a constant discharge. The active layer must be thicker in the model to account for river bed variations due to higher‐frequency discharge variations that are filtered out in morphological modelling. We discuss limitations in calibration and verification, but argue that, notwithstanding these limitations, 2D and 3D morphological models are valuable tools, not only for pragmatic applications to engineering problems, but also for revealing the limitations of established knowledge and understanding of the relevant physical processes. The application of numerical models to the Pannerden bifurcation appeared to reveal shortcomings in established model formulations that do not pose particular problems in other cases. This application is therefore particularly useful for setting the agenda for further research. Copyright © 2012 John Wiley & Sons, Ltd.