In memory of G.P. Kalinin (1916–1975)

Chronicles

In memory of G.P. Kalinin (1916–1975)     

P–PINI: A cosmogenic nuclide burial dating method for landscapes undergoing non-steady erosion

Existing methods of cosmogenic nuclide burial dating perform well provided that sediment sources undergo steady rates of erosion and the samples experience continuous exposure to cosmic rays. These premises exert important limitations on the applicability of the methods. And yet, high mountain sediment sources are rife with transient processes, such as non-steady erosion by glacial quarrying and/or landsliding, or temporary cosmic-ray shielding beneath glaciers and/or sediment. As well as breaching the premises of existing burial dating methods, such processes yield samples with low nuclide abundances and variable ²⁶Al/¹⁰Be ratios that may foil both isochron and simple burial-age solutions. P–PINI (Particle-Pathway Inversion of Nuclide Inventories) is a new dating tool designed for dating the burial of sediments sourced from landscapes characterized by abrupt, non-steady erosion, discontinuous exposure, and catchments with elevation-dependent ²⁶Al/¹⁰Be production ratios. P–PINI merges a Monte Carlo simulator with established cosmogenic nuclide production equations to simulate millions of samples (¹⁰Be–²⁶Al inventories). The simulated samples are compared statistically with ¹⁰Be–²⁶Al measured in field samples to define the most probable burial age. Here, we target three published ¹⁰Be–²⁶Al datasets to demonstrate the versatility of the P–PINI model for dating fluvial and glacial sediments. (1) The first case serves as a robust validation of P–PINI. For the Pulu fluvial gravels (China), we obtain a burial age of 1.27 ± 0.10 Ma (1σ), which accords with the isochron burial age and two independent chronometers reported in Zhao et al. (2016) Quaternary Geochronology 34, 75–80. The second and third cases, however, reveal marked divergence between P–PINI and isochron-derived ages. (2) For the fluvial Nenana Gravel (USA), we obtain a minimum-limiting burial age of 4.5 ± 0.7 Ma (1σ), which is compatible with unroofing of the Alaska Range starting ∼ 6 Ma, while calling into question the Early Pleistocene isochron burial age presented in Sortor et al. (2021) Geology 49, 1473–1477. (3) For the Bünten Till (Switzerland), we obtain a limiting burial age of <204 ka (95th percentile range), which conforms with the classical notion of the most extensive glaciation in the northern Alpine Foreland assigned to the Riss glaciation (sensu marine isotope stage 6) contrary to the isochron burial age presented in Dieleman et al. (2022) Geosciences, 12, 39. Discrepancies between P–PINI and the isochron ages are rooted in the challenges posed by the diverse pre-burial ²⁶Al/¹⁰Be ratios produced under conditions characteristic of high mountain landscapes; i.e., non-steady erosion, discontinuous cosmic-ray exposure, and elevation-dependent ²⁶Al/¹⁰Be production ratios in the source region, which are incompatible with the isochron method, but easily accommodated by the stochastic design of P–PINI.     

Young eclogite from the Greater Himalayan Sequence,Arun Valley,eastern Nepal: P–T–t path and tectonic implications

Garnet geochronology was used to provide the first direct measurement of the timing of eclogitization in the central Himalaya. Lu–Hf dates from garnet separates in one relict eclogite from the Arun River Valley in eastern Nepal indicate an age of 20.7 ± 0.4 Ma, significantly younger than ultra-high pressure eclogites from the western Himalaya, reflecting either different origins or substantial time lags in tectonics along strike. Four proximal garnet amphibolites from structurally lower horizons are 14–15 Ma, similar to post-eclogitization ages published for rocks along strike in southern Tibet. P–T calculations indicate three metamorphic episodes for the eclogite: i) eclogite-facies metamorphism at ～ 670 °C and ≥ 15 kbar at 23–16 Ma; ii) a peak-T granulite event at ～ 780 °C and 12 kbar; and iii) late-stage amphibolite-facies metamorphism at ～ 675 °C and 6 kbar at ～ 14 Ma. The garnet amphibolites were metamorphosed at ～ 660 °C. Three models are considered to explain the observed P–T–t evolution. The first assumes that the Main Himalayan Thrust (basal thrust of the Himalayan thrust system) cuts deeper at Arun than elsewhere. While conceptually the simplest, this model has difficulty explaining both the granulite-facies overprint and the pulse of exhumation between 25 and 14 Ma. A second model assumes that (aborted) subduction, slab breakoff, and ascent of India's leading edge occurred diachronously: ～ 50 Ma in the western Himalaya, ～ 25 Ma in the central Himalaya of Nepal, and presumably later in the eastern Himalaya. This model explains the P–T–t path, particularly heating during initial exhumation, but implies significant along-strike diachroneity, which is generally lacking in other features of the Himalaya. A third model assumes repeated loss of mantle lithosphere, first by slab breakoff at ～ 50 Ma, and again by delamination at ～ 25 Ma; this model explains the P–T–t path, but requires geographically restricted tectonic behavior at Arun. The P–T–t history of the Arun eclogites may imply a change in the physical state of the Himalayan metamorphic wedge at 16–25 Ma, ultimately giving rise to the Main Central Thrust by 15–16 Ma.     

Cycles revealed in Jurassic paleoclimatic data (ca. 200–145 Ma B.P.)

Limestone varves from the Todilto Basin, New Mexico (United States), are studied. These deposits cover a time interval of 1592 years and belong to the Jurassic (ca. 200–145 Ma B.P.). Via the construction of a combined spectral periodogram, quasiharmonic components are revealed in the initial data. The periods of these components are close to the basic contemporary solar activity cycles: the periods of 13 and 20 years agree with the Schwabe and Hale cycles; that of 78 years is in accord with the Gleissberg cycle (with both branches of 69 and 105 years represented); and the components of 179 and 235 years can be probably related to the Suess or de Vries cycles. The results indicate the influence of solar activity on the terrestrial climate in the geological past.     

近距离观测的散射P′P′波

在震后2300～2450s的时间内,于震中距30°～50°之间发现了以前从未报道过的1Hz散射波。这些散射波可能是在上地幔和地壳发生的PKPbc到PKPbc的反方位角散射而产生的,为绘制地幔非均匀性的小尺度变化(10km)提供了一个新的手段。大孔径地震台阵(LASA)记录的阵列波束清晰地展示了散射波能量从噪声中逐步显露,在大约80s后达到最大振幅,并在150s后恢复到噪声水平。横向与径向慢度(ρt,ρr)的叠加显示在大约(2,-2)和(-2,-2)s/°存在两个峰值,表明这些波是沿着主弧路径(180°～360°)到达,而且明显是反方位角。基于上述观测资料,我们提出一种地幔和地表的PKPbc到PKPbc的散射机制,因为(1)这与散射波到时和散射波特有的慢度特征相吻合,(2)它的散射路径与之前观测到的深部地幔PK·KP散射路径类似(Chang and Cleary,1981)。观测到的上地幔散射波和PK·KP波均符合散射波的广义集合,我们称其为P′·d·P′,可以在地幔的任意深度d散射。     

P波偏振层析成像

论述了Ｐ波偏振层析成像方法。这是一个利用Ｐ波远震偏振资料去反演速度结构的方法,与走时层析成像相比,它有几个显著的优点：不受震源定位和发震时刻误差的影响;对深地幔速度结构不敏感而对接收器附近的速度结构和速度梯度最敏感,在这一意义上它与时反演是互补的。如果联合使用走时和偏振资料可以改善层析成像的结果。走时的变化对应于速度的变化,而偏振的变化则与速度梯度的变化相对应,因昆,要确定速度异常的边界,用偏振数     

Structure–soil–structure interaction: Literature review

The concept of structure–soil–structure dynamic interaction was introduced, and the research methods were discussed. Based on several documents, a systematic summary of the history and status of the structure–soil–structure dynamic interaction research that considers adjacent structures was proposed as a reference for researchers. This study is in the initial stage, given its complexity and excessive simplification of the model for soil and structures, and should be carried forward for its significance. An attempt was made to summarize the common major computer programs in this area of study. Furthermore, the advantages, disadvantages, and applicability of such programs were discussed. The existing problems and the future research trend in this field were also examined.     

Mixed fluid sources involved in diamond growth constrained by Sr–Nd–Pb–C–N isotopes and trace elements

Sub-micrometer inclusions in diamonds carry high-density fluids (HDF) from which the host diamonds have precipitated. The chemistry of these fluids is our best opportunity of characterizing the diamond-forming environment. The trace element patterns of diamond fluids vary within a limited range and are similar to those of carbonatitic/kimberlitic melts that originate from beneath the lithospheric mantle. A convecting mantle origin for the fluid is also implied by C isotopic compositions and by a preliminary Sr isotopic study (Akagi, T., Masuda, A., 1988. Isotopic and elemental evidence for a relationship between kimberlite and Zaire cubic diamonds. Nature 336, 665–667.). Nevertheless, the major element chemistry of HDFs is very different from that of kimberlites and carbonatites, varying widely and being characterized by extreme K enrichment (up to ～ 39 wt.% on a water and carbonate free basis) and high volatile contents. The broad spectrum of major element compositions in diamond-forming fluids has been related to fluid–rock interaction and to immiscibility processes.Elemental signatures can be easily modified by a variety of mantle processes whereas radiogenic isotopes give a clear fingerprint of the time-integrated evolution of the fluid source region. Here we present the results of the first multi radiogenic-isotope (Sr, Nd, Pb) and trace element study on fluid-rich diamonds, implemented using a newly developed off-line laser sampling technique. The data are combined with N and C isotope analysis of the diamond matrix to better understand the possible sources of fluid involved in the formation of these diamonds. Sr isotope ratios vary significantly within single diamonds. The highly varied but unsupported Sr isotope ratios cannot be explained by immiscibility processes or fluid-mineral elemental fractionations occurring at the time of diamond growth. Our results demonstrate the clear involvement of a mixed fluid, with one component originating from ancient incompatible element-enriched parts of the lithospheric mantle while the trigger for releasing this fluid source was probably carbonatitic/kimberlitic melts derived from greater depths. We suggest that phlogopite mica was an integral part of the enriched lithospheric fluid source and that breakdown of this mica releases K and radiogenic Sr into a fluid phase. The resulting fluids operate as a major metasomatic agent in the sub-continental lithospheric mantle as reflected by the isotopic composition and trace element patterns of G10 garnets.     

P波广角入射的透射P波特征研究（英文）

P波入射到地下反射界面时,如果入射角超过临界角,透射P波可能会出现极性反转。本文首先从Zoeppritz方程出发,研究了广角入射条件下的反射系数和透射系数特征,发现满足一定条件时,随着入射角的增加,透射P波系数的特征曲线从第一象限经过原点进入第三象限,透射P波系数发生极性反转,它的辐角发生突变。文中给出了透射P波系数为零时入射角的表达式及其近似公式,通过对单界面模型进行有限差分正演,证明了透射P波系数为零的存在性。储层从饱含水到饱含气变化时,储层的P波速度和密度明显降低,通过单界面储层模型的透射P波的合成记录,说明了储层饱含气时透射P波系数可能发生极性反转。经理论分析和数值模拟得到了一些结论,这些结论为井间地震勘探的应用奠定了理论基础。     

Laplace–Fourier-Domain Waveform Inversion for Fluid–Solid Media

Full waveform inversion algorithms are widely used in the construction of subsurface velocity models. In the following study, we propose a Laplace–Fourier-domain waveform inversion algorithm that uses both Laplace-domain and Fourier-domain wavefields to achieve the reconstruction of subsurface velocity models. Although research on the Laplace–Fourier-domain waveform inversion has been published recently that study is limited to fluid media. Because the geophysical targets of marine seismic exploration are usually located within solid media, waveform inversion that is approximated to acoustic media is limited to the treatment of properly identified submarine geophysical features. In this study, we propose a full waveform inversion algorithm for isotropic fluid–solid media with irregular submarine topography comparable to a real marine environment. From the fluid–solid system, we obtained P and S wave velocity models from the pressure data alone. We also suggested strategies for choosing complex frequency bands constructed of frequencies and Laplace coefficients to improve the resolution of the restored velocity structures. For verification, we applied our Laplace–Fourier-domain waveform inversion for fluid–solid media to synthetic data that were reconstructed for fluid–solid media. Through this inversion test, we successfully restored reasonable velocity structures. Furthermore, we successfully extended our algorithm to a field data set.     

Challenges in operationalizing the water–energy–food nexus

Concerns about the water–energy–food (WEF) nexus have motivated many discussions regarding new approaches for managing water, energy and food resources. Despite the progress in recent years, there remain many challenges in scientific research on the WEF nexus, while implementation as a management tool is just beginning. The scientific challenges are primarily related to data, information and knowledge gaps in our understanding of the WEF inter-linkages. Our ability to untangle the WEF nexus is also limited by the lack of systematic tools that could address all the trade-offs involved in the nexus. Future research needs to strengthen the pool of information. It is also important to develop integrated software platforms and tools for systematic analysis of the WEF nexus. The experience made in integrated water resources management in the hydrological community, especially in the framework of Panta Rhei, is particularly well suited to take a lead in these advances.     

上地幔过渡带S—P和P—S转换的全球研究

本文研究了1984年至1986年3年期间全球数字地震台网(GDSN)记录的长周期数据,以讨论作为远震S前驱波的上地幔过渡带产生的S—P和P—S转换波的存在,这类转换波已在大量单台站记录上得到证实,许多记录的简单叠加加强了转换相位的出现,并证明不同板块区域之间过渡带性质没有重要的横向变化。400km间断面的S—P和P—S转换波在70°至85°之间距离上得到很好的观测,而670km间断面的转换波在87°以上距离观测得很好。已发表的震源机制分析和与合成地震图比较说明,转换震相的形态主要由地震辐射图型所控制,接收点下方经历S—P转换的震相可以从倾滑事件观测到。这些事件沿台站方向上辐射强SV和弱P波。在另一方面震源区域下方P—S转换经常可以从沿台站方向激发强P波和很小的SV波能量的事件中观测到,同时也从一些走滑事件中观测到。观测资料与合成地震图比较表明,Dziewonski和Anderson(1981)的PREM模型可以解释绝大部分观测数据。然而,从670km间断面上观测到的S—P和P—S转换波经常比合成地震图的振幅强。转换波和P波尾波的加强干涉,震源辐射的影响以及比在PREM模型中更高的670km间断面的速度反差都可能成为这一偏差的原因。     

P波地震图上的地面反射波与S—P转换波

辐射脉冲约１Ｓ或更短持续时间的震源理论研究表明：对于假定为实际地壳结构的双力偶震源，短周期Ｐ波地震图最突出的震相当是Ｐ，ＰＰ和ＳＰ波。已注意到显示Ｐ波与地面反射的某些地震图，并可以充分模拟，然而，由浅地震观测到的这种地震图数目很小，通常，观测到的地震图要私很简单，即它们显示出突出的初至，而很少有其他信息，要么地震图很复杂，即显示振幅大于或等于初至振幅的３０Ｓ以上的震相，简单和或复杂地震图上缺乏清晰     

A statistical study of magnetosphere–ionosphere coupling in the Lyon–Fedder–Mobarry global MHD model

The statistics of magnetosphere–ionosphere (MI) coupling derived from a two-month long run of the Lyon–Fedder–Mobarry (LFM) global simulation model are investigated. MI coupling characteristics such as polar cap potential and field-aligned current (FAC), downward Poynting flux and vorticity of ionospheric convection are compared with observed statistical averages and with results from the Weimer 05 empirical model. The comparisons for eight different IMF clock-angle orientations show that the LFM model produces reasonably accurate average distributions of the Region I and Region II currents. Both current systems have average amplitudes similar to those observed by the Iridium satellite constellation; however, the average LFM amplitudes are smaller by a factor of two compared with the values from the Weimer 05 model. The comparisons of polar cap potential show that the LFM model produces reasonable patterns of ionospheric convection, but the average cross polar cap potential (CPCP) is greater than the observed results by a factor of approximately 2 and greater than Weimer 05 by a factor of 1.5. The differences in convection in LFM results relative to the Weimer 05 model accounts for much of the difference in the Poynting flux patterns and integrated power produced by the two models. The comparisons of average ionospheric field-aligned vorticity show good agreement on the dayside; however, the LFM model gives higher nightside vorticity which may imply that the ionospheric conductance on the nightside is too small in the simulation.     

Rainfall–interception–evaporation–runoff relationships in a semi-arid catchment,northern Limpopo basin,Zimbabwe

Abstract

Two types of monthly water balance models at basin scale are used: PE models use precipitation and potential evapotranspiration (PET) as their observed input data, whereas P models need only precipitation. Calibration proceeds by comparing model runoff and observed runoff. Calibration is entirely automatic with the exclusion of subjective elements. All models differ only by their actual evapotranspiration equations. PE models from previous papers are generalized essentially by replacing the constant evapotranspiration parameter by a periodic one, thus increasing the number of parameters by two (a “parameter” is an unknown constant to be estimated, and which is a characteristic of the river basin to be described). P models use a periodic “driving force”, which is intended to represent periodicity of hydrological phenomena, normally originating in the (unavailable) PET time series. These eight PE models and three P models are then applied to 55 river basins in 10 countries with widely diverging climates and soil conditions. A marked improvement of model performance in about one third of the basins is due to the introduction of the above mentioned periodic functions. Even when PET data are available it is sometimes useful to consider P models. P models scarcely perform less well than PE models. An engineer, wanting to try out as few models as possible on a given river basin, can restrict his attention to the optimization of two or three models. The paper is an extension of a long effort towards monthly water balance models, and is believed to give a solution in most circumstances.     

Bayesian comparison of different rainfall depth–duration–frequency relationships

Depth–duration–frequency curves estimate the rainfall intensity patterns for various return periods and rainfall durations. An empirical model based on the generalized extreme value distribution is presented for hourly maximum rainfall, and improved by the inclusion of daily maximum rainfall, through the extremal indexes of 24 hourly and daily rainfall data. The model is then divided into two sub-models for the short and long rainfall durations. Three likelihood formulations are proposed to model and compare independence or dependence hypotheses between the different durations. Dependence is modelled using the bivariate extreme logistic distribution. The results are calculated in a Bayesian framework with a Markov Chain Monte Carlo algorithm. The application to a data series from Marseille shows an improvement of the hourly estimations thanks to the combination between hourly and daily data in the model. Moreover, results are significantly different with or without dependence hypotheses: the dependence between 24 and 72 h durations is significant, and the quantile estimates are more severe in the dependence case.     

Regional scale rainfall–runoff modeling using VARX–MGARCH approach

Multivariate time series modeling approaches are known as useful tools for describing, simulating, and forecasting hydrologic variables as well as their changes over the time. These approaches also have temporal and cross-sectional spatial dependence in multiple measurements. Although the application of multivariate linear and nonlinear time series approaches such as vector autoregressive with eXogenous variables (VARX) and multivariate generalized autoregressive conditional heteroscedasticity (MGARCH) models are commonly used in financial and economic sciences, these approaches have not been extensively used in hydrology and water resources engineering. This study employed VARX and VARX–MGARCH approaches in modeling mean and conditional heteroscedasticity of daily rainfall and runoff records in the basin of Zarrineh Rood Dam, Iran. Bivariate diagonal VECH (DVECH) model, as a main type of MGARCH, shows how the conditional variance–covariance and conditional correlation structure vary over the time between residuals series of the fitted VARX. For this purpose, five model fits, which consider different combinations of twofold rainfall and runoff, including both upstream and downstream stations, have been investigated in the present study. The VARX model, with different orders, was applied to the daily rainfall–runoff process of the study area in each of these model fits. The Portmanteau test revealed the existence of conditional heteroscedasticity in the twofold residuals of fitted VARX models. Therefore, the VARX–DVECH model is proposed to capture the heteroscedasticity existing in the daily rainfall–runoff process. The bivariate DVECH model indicated both short-run and long-run persistency in the conditional variance–covariance matrix related to the twofold innovations of rainfall–runoff processes. Furthermore, the evaluation criteria for the VARX–DVECH model revealed the improvement of VARX model performance.     

Incremental dynamic analysis and FEMA P695 seismic performance evaluation of a cold-formed steel–framed building with gravity framing and architectural sheathing

The objective of this paper is to present incremental dynamic analysis (IDA) and seismic performance evaluation results for a two-story cold-formed steel (CFS)–framed building. The archetype building was designed to current U.S. standards and then subjected to full-scale shake table tests under the U.S. National Science Foundation Network for Earthquake Engineering Simulation (NEES) program. Test results showed that the building's stiffness and capacity were considerably higher than expected and the building suffered only nonstructural damage even at excitations in excess of Maximum Considered Earthquake levels for a high seismic zone. For the archetype building, three-dimensional finite element models at different modeling fidelity levels were created using OpenSees. The models are subjected to IDA using the far-field ground motion records prescribed in Federal Emergency Management Agency (FEMA) P695. Seismic performance quantification following the FEMA P695 procedure shows that if the modeling fidelity only follows the state-of-the-practice, ie, only includes shear walls, unsafe collapse margin ratios are predicted. State-of-the-art models that account for participation from CFS gravity walls and architectural sheathing have overall performance that are consistent with testing, and IDA results indicate acceptable collapse margin ratios, predicated primarily on large system overstrength. Neglecting the lateral force resistance of the gravity system and nonstructural components, as done in current design, renders a safe design in the studied archetype, but largely divorced from actual system behavior. The modeling protocols established here provide a means to analyze a future suite of CFS-framed archetype buildings for developing further insight on the seismic response modification coefficients for CFS-framed buildings.     

Intensity–duration–frequency curves exploiting neighbouring extreme precipitation data

ABSTRACT

Records of precipitation extremes are essential for hydrological design. In urban hydrology, intensity–duration–frequency curves are typically estimated from observation records. However, conventional approaches seldom consider the areal extent of events. If they do, duration-dependent area reduction factors are used, but precipitation is measured at only a few locations. Due to the high spatial variability of precipitation, it is relatively unlikely that a gauged observation network will capture the extremes that occur during a precipitation event. Therefore, the area reduction approach cannot be regarded as the reduction of an observed maximum. To investigate precipitation extremes, spatial aspects need to be considered using different approaches. Here, we both address the conventional practice of area reduction and consider a within-area chance of increased precipitation, defined as the maximum precipitation intensity observed in a cluster within a selected domain. The results show that (1) the risk of urban flooding is routinely underestimated in current design practice, and (2) traditional calculations underestimate extremes by as much as 30–50%. We show how they can be revised sensibly.