Holocene hydrological variability of Lake Ladoga,northwest Russia,as inferred from diatom oxygen isotopes

This article presents a new comprehensive assessment of the Holocene hydrological variability of Lake Ladoga, northwest Russia. The reconstruction is based on oxygen isotopes of lacustrine diatom silica (δ¹⁸O_diatom) preserved in sediment core Co 1309, and is complemented by a diatom assemblage analysis and a survey of modern isotope hydrology. The data indicate that Lake Ladoga has existed as a freshwater reservoir since at least 10.8 cal. ka BP. The δ¹⁸O_diatom values range from +29.8 to +35.0‰, and relatively higher δ¹⁸O_diatom values around +34.7‰ between c. 7.1 and 5.7 cal. ka BP are considered to reflect the Holocene Thermal Maximum. A continuous depletion in δ¹⁸O_diatom since c. 6.1 cal. ka BP accelerates after c. 4 cal. ka BP, indicating Middle to Late Holocene cooling that culminates during the interval 0.8–0.2 cal. ka BP, corresponding to the Little Ice Age. Lake‐level rises result in lower δ¹⁸O_diatom values, whereas lower lake levels cause higher δ¹⁸O_diatom values. The diatom isotope record gives an indication for a rather early opening of the Neva River outflow at c. 4.4–4.0 cal. ka BP. Generally, overall high δ¹⁸O_diatom values around +33.5‰ characterize a persistent evaporative lake system throughout the Holocene. As the Lake Ladoga δ¹⁸O_diatom record is roughly in line with the 60°N summer insolation, a linkage to broader‐scale climate change is likely.     

Another Year of Record Heat for the Oceans

Changes in ocean heat content(OHC), salinity, and stratification provide critical indicators for changes in Earth’s energy and water cycles. These cycles have been profoundly altered due to the emission of greenhouse gasses and other anthropogenic substances by human activities, driving pervasive changes in Earth’s climate system. In 2022, the world’s oceans, as given by OHC, were again the hottest in the historical record and exceeded the previous 2021 record maximum.According to IAP/CAS data, ...     

冰力学参数的超声波测试研究

利用非线性高能超声测试设备及超声波(纵波、横波)波速与物体力学参数的关系,对人造冰样进行了冰样力学参数(杨氏模量、泊松比、剪切模量、体积模量)随温度变化的研究。通过MATLAB进行所测数据曲线拟合,得到超声波波速在人造冰样中随温度的变化规律,进而由理论公式推导所测人造冰样力学参数随温度的变化规律。结果表明:冰样中超声波波速随温度降低而升高,冰样的杨氏模量、泊松比、剪切模量、体积模量也都随温度降低而升高。本研究有助于超声波检测法在冰样物理力学性质测量中的应用,为开展南极冰盖、海冰以及终年冻土等力学及流动特性研究提供理论模型和实验数据。     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

Reflected waves in finely layered tilted orthorhombic media

Upscaling in seismics is a homogenization of finely layered media in the zero-frequency limit. An upscaling technique for arbitrary anisotropic layers has been developed by Schoenberg and Muir. Applying this technique to a stack of layers of orthorhombic (ORT) symmetry whose vertical symmetry planes are aligned, results in an effective homogeneous layer with orthorhombic symmetry. If the symmetry planes in a horizontal orthorhombic layer are rotated with respect to vertical, the medium is referred to as tilted orthorhombic (TOR) medium, and the stack composed of TOR layers in zero-frequency limit will produce an effective medium of a lower symmetry than orthorhombic. We consider a P-wave that propagates through a stack of thin TOR layers, then it is reflected (preserving the mode) at some interface below the stack, and then propagates back through the same stack. We propose to use a special modified medium for the upscaling in case of this sequential down- and up-propagation: each TOR layer in the stack is replaced by two identical TOR layers whose tilt angles have the opposite algebraic sign. In this modified medium, one-way propagation of a seismic wave (any wave mode) is equivalent to propagation of a pure-mode reflection in the original medium. We apply this idea to study the contribution from an individual layer from the stack and show how the approach can be applied to a stack of TOR layers. To demonstrate the applicability of the model, we use well log data for the upscaling. The model we propose for the upscaling can be used in well-seismic ties to correct the effective parameters obtained from well log data for the presence of tilt, if latter is confirmed by additional measurements (for example, borehole imaging).     

Estimation of elemental abundances of the lunar regolith using clementine UVVIS+NIR data

In this study we propose a regression model for the estimation of lunar elemental abundances from spectral features extracted from Clementine multispectral imagery in the visible and near-infrared domain. We extract a set of spectral features, including the continuum slope, the FWHM of the ferrous absorption trough near 1000 nm, and the wavelengths and relative depths of the absorption minima and inflection points present in the trough. As a “ground truth” for the elemental abundances we rely on the Lunar Prospector gamma ray spectrometer (LP GRS) data. With respect to the elemental abundances of the Apollo and Luna landing sites independently derived from returned samples, the best examined regression model is a second-order polynomial. The proposed regression-based approach allows an estimation of the elemental abundances of Ca, Al, Fe, Mg, and O at an accuracy of about 1 wt% and some tenths of a weight percent for Ti. We examine the influence of calibration of the Clementine UVVIS+NIR data and find that its effect on the results obtained with the regression approach is minor. Furthermore, we define a three-endmember model which allows the petrographic mapping of the lunar surface materials in terms of their Fe, Mg, and Al abundances. We examine the global distribution of Mg-rich rocks, the distribution of cryptomaria, and the occurrence of aluminous mare basalts in the Frigoris region. A possible regional compositional anomaly in northwestern Oceanus Procellarum is found, which corresponds to an extended area displaying spectral characteristics consistent with mare basalt containing significant amounts of olivine. On local scales, we examine in terms of our regression model the highland craters Proclus and Tycho, the compositionally anomalous central peaks of the craters Copernicus and Bullialdus, and the pyroclastic deposits on the floor of Alphonsus and on the northern rim of Petavius. As a general result, we show that the regression-based approach allows the detection of the main lunar terrain classes and rock types based on multispectral imagery in the visible and near-infrared domain.     

Thermal anomalies associated with shallow gas hydrates in the K-2 mud volcano, Lake Baikal

Thermal measurements and hydrate mapping in the vicinity of the K-2 mud volcano in Lake Baikal have revealed a particular type of association of thermal anomalies (29–121?mW?m^–2) near hydrate-forming layers. Detailed coring within K-2 showed that hydrates are restricted to two distinct zones at sub-bottom depths exceeding 70–300?cm. Temperature data from stations with hydrate recovery and degassing features all display low thermal gradients. Otherwise, the thermal gradients within the mud volcano are generally increased. These findings imply a more complicated thermal regime than often assumed for mud volcanoes, with important roles for both fluids and hydrates. The coexistence of neighbouring low and high thermal anomalies is interpreted to result from discharging and recharging fluid activity, rather than hydrate thermodynamics. It is suggested that hydrates play a key role in controlling the fluid circulation pattern at an early stage. At a later stage, the inflow of undersaturated lake water would favour the dissolution of structure I hydrates and the formation of structure II hydrates, the latter having been observed on top of structure I hydrates in the K-2 mud volcano.     

Stability of phlogopite in ultrapotassic kimberlite-like systems at 5.5–7.5 GPa

Hydrous K-rich kimberlite-like systems are studied experimentally at 5.5–7.5 GPa and 1200–1450?°C in terms of phase relations and conditions for formation and stability of phlogopite. The starting samples are phlogopite–carbonatite–phlogopite sandwiches and harzburgite–carbonatite mixtures consisting of Ol?+?Grt?+?Cpx?+?L (±Opx), according to the previous experimental results obtained at the same P–T parameters but in water-free systems. Carbonatite is represented by a K- and Ca-rich composition that may form at the top of a slab. In the presence of carbonatitic melt, phlogopite can partly melt in a peritectic reaction at 5.5 GPa and 1200–1350?°C, as well as at 6.3–7.0 GPa and 1200?°C: 2Phl?+?CaCO₃ (L)?Cpx?+?Ol?+?Grt?+?K₂CO₃ (L)?+?2H₂O (L). Synthesis of phlogopite at 5.5 GPa and 1200–1350?°C, with an initial mixture of H₂O-bearing harzburgite and carbonatite, demonstrates experimentally that equilibrium in this reaction can be shifted from right to left. Therefore, phlogopite can equilibrate with ultrapotassic carbonate–silicate melts in a?≥?150?°C region between 1200 and 1350?°C at 5.5 GPa. On the other hand, it can exist but cannot nucleate spontaneously and crystallize in the presence of such melts in quite a large pressure range in experiments at 6.3–7.0 GPa and 1200?°C. Thus, phlogopite can result from metasomatism of peridotite at the base of continental lithospheric mantle (CLM) by ultrapotassic carbonatite agents at depths shallower than 180–195 km, which creates a mechanism of water retaining in CLM. Kimberlite formation can begin at 5.5 GPa and 1350?°C in a phlogopite-bearing peridotite source generating a hydrous carbonate–silicate melt with 10–15 wt% SiO₂, Ca# from 45 to 60, and high K enrichment. Upon further heating to 1450?°C due to the effect of a mantle plume at the CLM base, phlogopite disappears and a kimberlite-like melt forms with SiO₂ to 20 wt% and Ca#?=?35–40.     

Kinematical characteristics of the factorized velocity model

The factorized velocity model that incorporates both vertical heterogeneity and constant anisotropy is one of the complicated analytical models used in seismic data processing and interpretation. In this paper, I derive the analytic equations for offset, traveltime and relative geometrical spreading for the quasi‐compressional (qP‐) waves that can be used for modelling and inversion of the traveltime parameters. I show that the presence of anelliptic anisotropy usually dominates over the vertical heterogeneity with respect to the non‐hyperbolicity of the factorized velocity model.