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.     

Detecting,sourcing, and age-dating dredged sediments on the open shelf,southern California,using dead mollusk shells

Molluscan shell debris is an under-exploited means of detecting, sourcing, and age-dating dredged sediments in open-shelf settings. Backscatter features on the Southern California shelf are suggestive of dredged sediment hauled from San Diego Bay but deposited significantly inshore of the EPA-designated ocean disposal site. We find that 36% of all identifiable bivalve shells > 2 mm (44% of shells > 4 mm) in sediment samples from this 'short dump' area are from species known to live exclusively in the Bay; such shells are absent at reference sites of comparable water depth, indicating that their presence in the short-dump area signals non-compliant disposal rather than natural offshore transport or sea level rise. These sediments lack the shells of species that invaded California bays in the 1970s, suggesting that disposal preceded federal regulations. This inexpensive, low-tech method, with its protocol for rejecting alternative hypotheses, will be easy to adapt in other settings.     

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.     

Can suspended sediment concentrations be estimated from multispectral imagery using only image-derived information?

Sampling for suspended sediment concentrations (SSC) in inland waters is traditionally based on collecting samples at sparse locations and in limited intervals. A number of investigators explored the utility of earth-observing satellites and air-borne sensors for monitoring of SSC over vast areas. Two approaches are commonly deployed: (1) empirical relationships between a chosen remotely sensed quantity and the actual in-situ SSC; and (2) bio-optical models founded on radiative transfer modeling. Unfortunately, in-situ measurements are often unavailable for direct image calibration, and inherent optical properties of optically active constituents (specific scattering and absorption coefficients) are usually unknown. This paper examines the possibility to retrieve SSC from multispectral satellite imagery without any in-situ data, i.e. using only image-derived information. The fundamental principle of image selfcalibration relies on the fact that in the visual domain of wavelengths (∼400–700 nm) the at-sensor reflectance becomes “saturated“ at high SSC, whereas the near-infrared domain (∼700–900 nm) remains almost perfectly linearly related to sediment concentrations. The core idea of the self-calibrating procedure is rather simple and is based on fitting an exponential function between reflectance and SSC, with SSC replaced by a linear relationship between SSC and reflectance in the near-infrared domain. As a first approximation of the non-linearity between reflectance and SSC levels in the 400–700 nm range, we used the equation proposed by Schiebe et al. (1992), although other equations, especially those arising from optical theory could be used as well. The technique is illustrated on a moderately sediment-laden reservoir and two scenes acquired from Landsat ETM+. The standard error of the estimated SSC was below 15 mg/L (i.e. ∼25 % relative error for the observed range of SSC). Although the proposed algorithm does not yield better results than other models mentioned in the literature, the primary advantage of the outlined methodology is that no in-situ measurements (water sampling nor spectral profiling) are needed — i.e. only image-derived information is used.     

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.     

Effects of social and climatic factors on building activity in the Czech lands between 1450 and 1950: a dendrochronological analysis

The development of settlement and building activity is the result of socioeconomic, political and demographic changes in the past. However, accurate information on temporal variation in building activity is rather limited. Dendrochronological databases containing dated historical wooden constructions provide an important resource. We used 6514 tree-felling dates to reconstruct building activity in the Czech lands for the period 1450–1950. Comparing felling dates with historical events demonstrated that building activity was negatively associated with intense wars, particularly during the Thirty Years' War (1618–1648). After the Peace of Westphalia (1648), socioeconomic renewal and demographic growth were reflected in an upsurge of building activity, especially ecclesiastical buildings. While the construction of ecclesiastical and noble buildings culminated around the 1720s, rural buildings peaked in the 1780s and the 1820s. Although no direct effect of climate was demonstrated, adverse climatic conditions leading to harvest failures and subsequent famines (e.g. the ‘Hunger Years’ 1770–1772) significantly contributed to declines in building activity. In contrast, a higher number of felling dates were detected when strong and/or frequent windstorms occurred. This study provides a comprehensive understanding of building activity in Central Europe and advocates the use of dendrochronological databases for the investigation of human activities in history.