Does Poisson’s downward continuation give physically meaningful results?

The downward continuation (DWC) of the gravity anomalies from the Earth’s surface to the geoid is still probably the most problematic step in the precise geoid determination. It is this step that motivates the quasi-geoid users to opt for Molodenskij’s rather than Stokes’s theory. The reason for this is that the DWC is perceived as suffering from two major flaws: first, a physically meaningful DWC technique requires the knowledge of the irregular topographical density; second, the Poisson DWC, which is the only physically meaningful technique we know, presents itself mathematically in the form of Fredholm integral equation of the 1st kind. As Fredholm integral equations are often numerically ill-conditioned, this makes some people believe that the DWC problem is physically ill-posed. According to a revered French mathematician Hadamard, the DWC problem is physically well-posed and as such gives always a finite and unique solution. The necessity of knowing the topographical density is, of course, a real problem but one that is being solved with an ever increasing accuracy; so sooner or later it will allow us to determine the geoid with the centimetre accuracy.     

Metal speciation and attenuation in stream waters and sediments contaminated by landfill leachate

The degree of metal contamination (Zn, Pb, Cu, Ni, Cd) has been investigated in the vicinity of an old unmonitored municipal landfill in Prague, Czech Republic, where the leachate is directly drained into a surface stream. The water chemistry was coupled with investigation of the stream sediment (aqua regia extract, sequential extraction, voltammetry of microparticles) and newly formed products (SEM/EDS, XRD). The MINTEQA2 speciation-solubility calculation showed that the metals (Zn, Pb, Cu, Ni) are mainly present as carbonate complexes in leachate-polluted surface waters. These waters were oversaturated with respect to Fe(III) oxyhydroxides, calcite (CaCO₃) and other carbonate phases. Three metal attenuation mechanisms were identified in leachate-polluted surface waters: (i) spontaneous precipitation of metal-bearing calcite exhibiting significant concentrations of trace elements (Fe, Mn, Mg, Sr, Ba, Pb, Zn, Ni); (ii) binding to Fe(III) oxyhydroxides (mainly goethite, FeOOH) (Pb, Zn, Cu, Ni); and (iii) preferential bonding to sediment organic matter (Cu). These processes act as the key scavenging mechanisms and significantly decrease the metal concentrations in leachate-polluted water within 200 m from the direct leachate outflow into the stream. Under the near-neutral conditions governing the sediment/water interface in the landfill environment, metals are strongly bound in the stream sediment and remain relatively immobile.     

Study of the filamentary infrared dark cloud G192.76+00.10 in the S254–S258 OB complex

We present results of a high resolution study of the filamentary infrared dark cloud G192.76+00.10 in the S254-S258 OB complex in several molecular species tracing different physical conditions. These include three isotopologues of carbon monoxide(CO), ammonia(NH3) and carbon monosulfide(CS). The aim of this work is to study the general structure and kinematics of the filamentary cloud, and its fragmentation and physical parameters. The gas temperature is derived from the NH3(J, K) =(1,1),(2, 2) and ~(12)CO(2-1) lines, and the ~(13)CO(1-0), ~(13)CO(2-1) emission is used to investigate the overall gas distribution and kinematics. Several dense clumps are identified from the CS(2-1)data. Values of the gas temperature lie in the range 10-35 K, and column density N(H2) reaches the value 5.1 x 10~(22) cm~(-2). The width of the filament is of order 1 pc. The masses of the dense clumps range from ～ 30 M_☉ to ～ 160 M_☉. They appear to be gravitationally unstable. The molecular emission shows a gas dynamical coherence along the filament. The velocity pattern may indicate longitudinal collapse.     

Multistep inversion workflow for 3D long‐offset damped elastic waves in the Fourier domain

We present a new workflow for imaging damped three‐dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short‐offset and low‐frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine‐scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite‐difference frequency‐domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full‐waveform inversion workflow is applied to the 3D Marmousi‐2 and SEG/EAGE Salt models with long‐offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one‐third of the maximum offset length at a resolution consistent with the inverted frequencies.     

Ray tracing and geodesic deviation of the SH and SV reference rays in a heterogeneous generally anisotropic medium which is approximately uniaxial

The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. If we know that a medium is close to uniaxial (transversely isotropic), it may be advantageous to trace reference rays which resemble the SH–wave and SV–wave rays. This paper is devoted to defining and tracing these SH and SV reference rays of elastic S waves in a heterogeneous generally anisotropic medium which is approximately uniaxial (approximately transversely isotropic), and to the corresponding equations of geodesic deviation (dynamic ray tracing). All presented equations are simultaneously applicable to ordinary and extraordinary reference rays of electromagnetic waves in a generally bianisotropic medium which is approximately uniaxially anisotropic. The improvement of the coupling–ray–theory seismograms calculated along the proposed SH and SV reference rays, compared to the coupling–ray–theory seismograms calculated along the anisotropic common reference rays, has already been numerically demonstrated by the authors in four approximately uniaxial velocity models.     

Palaeo- and rock-magnetic investigations across Jurassic-Cretaceous boundary at St Bertrand’s Spring,Drôme,France: applications to magnetostratigraphy

Palaeo- and rock-magnetic investigations of the St Bertrand’s Spring (Le Ravin de Font de St Bertrand) locality in France were carried out in order to contribute to, and improve, the stratigraphy of the Jurassic-Cretaceous boundary interval. Magnetic susceptibility shows slightly diamagnetic behaviour in the lowermost part of the profile and an increase (paramagnetic) towards its middle and upper parts. Rock-magnetic measurements throughout the section show magnetite as the main magnetic fraction, together with traces of hematite. Additionally, thermal demagnetization indicates the presence of goethite. Our magnetostratigraphy indicates three normal/reversed polarity sequences; possibly encompassing the magnetozones M19r to the M17n. This suggests that the St Bertrand section straddles the Tithonian/Berriasian boundary and reaches the middle Berriasian sensu lato.     

Landslide-dammed lake sediment volume calculation using waterborne ERT and SONAR profiling

Lake sediment volume calculation is a challenging task, namely in cases when detailed drilling is complicated, expensive, or impossible, information on the pre-sedimentation surface unavailable, and record of siltation rate non-existent or too short. This study shows how waterborne, non-invasive geophysical survey, such as electrical resistivity tomography (ERT) can be very effective in acquiring the missing data, namely when combined with sound navigation ranging (SONAR) water depth measurements and supported by information from auxiliary sources. However, ERT surveying in water environment requires specific approaches, as we illustrate on the case of the Mladotice lake study. The lake was created after a landslide in May 1872, and since its formation, the depth has gradually decreased due to sedimentation. We have reconstructed the original surface, calculated the sediment volume, and compiled information on sedimentation to estimate its remaining life span. To achieve this, we measured nine waterborne ERT profiles across the lake. To reach the necessary depth, all ERT profiles were extended on land and crossed the lake using custom-built flotation pads. ERT profiling was combined with SONAR depth measurements, historical bathymetric surveys, borehole core analysis, sediment flux measurements, volumetric calculations, and water conductivity probing. The study has achieved three main results. First, practical applicability and advantages of stationary waterborne ERT profiling in combination with bathymetric sounding were demonstrated. Second, the original lake volume and accumulated sediment was calculated. We estimate that the volume of lake sediment is 187 000 m³, two-thirds of the original lake volume (over 275 000 m³). Finally, based on three volumetric data sets from 1972, 2003, and 2017, and recent monitoring of the sediment inflow, we propose scenarios of lake filling and its future development. Most interestingly, the sedimentation rate has decreased significantly in the last 20 years, suggesting that the lake may survive much longer than hitherto expected. © 2020 John Wiley & Sons, Ltd.     

The Lateglacial and Holocene in Central Europe: a multi‐proxy environmental record from the Bohemian Forest,Czech Republic

The Hercynian mountain ranges were islands of mountain glaciation and alpine tundra in a Central European ice‐free corridor during the Late Pleistocene. Today they are notable areas of glacial landforms, alpine‐forest free areas, peatlands and woodlands. However, our knowledge of the Lateglacial and early Holocene environmental changes in this region is limited. We present a new multi‐proxy reconstruction of a mid‐altitude environment in the Bohemian Forest spanning this period. A core (5.2 m length) in the ?erné Lake cirque (1028 m a.s.l.) was subjected to lithological, geochemical, pollen and macrofossil analysis supplemented by two optically stimulated luminescence (OSL) and 10 AMS radiocarbon dates. We determined the impact of regional and supraregional climate changes on the environment. The two most significant changes in sedimentation during the Lateglacial (17.6 and 15.8–15.5 cal. ka BP) were synchronous with regional glacial chronostratigraphy. Unlike Central European mountain ranges, in the Bohemian Forest the Younger Dryas was not coincident with glacier re‐advance, but was a dry, cold episode with low lake levels, which prevailed until the early Preboreal. Plant macrofossils indicate local establishment of Betula nana and Betula pendula/pubescens at 15.4–13.4 cal. ka BP. Comparison with Holocene records from Central Europe shows a similar immigration history of vegetation at mid and higher altitudes. The tree line exceeded an altitude of ~1000 m a.s.l. around 10.5 cal. ka BP and coincided with rapid geochemical changes in the sediment. The 8.2 ka BP event did not have any response in the sedimentary record, but corresponded to stabilization of the Picea abies population and expansion of Fagus. Fagus colonized the Bohemian Forest earlier than other Hercynian mid‐mountains, but never predominated in the composition of the forest at higher elevations. Abies alba was the last tree species that immigrated to the study area.