Three-dimensional poro-elasto-plastic model for wave-induced seabed response around submarine pipeline

The evaluation of the wave-induced excess pore pressure around a buried pipeline is particularly important for pipeline engineers involved in the design of offshore pipelines. Existing models for the wave-induced seabed response around submarine pipeline have been limited to poro-elastic soil behavior and de-coupled oscillatory and residual mechanisms for the rise in excess pore water pressure. To overcome the shortcoming of the existing models, in this study a three-dimensional poro-elasto-plastic soil model with submarine pipeline is established, in which both oscillatory and residual mechanisms can be simulated simultaneously. With the proposed model, a parametric study is conducted to investigate the relative differences of the predictions of the wave-induced pore pressure with poro-elasto-plastic model. Based on numerical examples, it can be concluded that the poro-elasto-plastic behaviors of soil have more significant influence on wave-induced pore pressure of seabed around submarine pipeline. As the seabed depth increases, the normalized pore pressures decrease rapidly at the upper part of seabed, and then change slightly at the lower part of the seabed. Soil permeability and wave period have obvious influence on the wave-induced normalized pore pressure.     

Precise Re–Os mineral isochron and Pb–Nd–Os isotope systematics of a mafic–ultramafic sill in the 2.0 Ga Onega plateau (Baltic Shield)

We present Re–Os, Sm–Nd and Pb–Pb isotope and trace element data for the Konchozero sill, a layered mafic–ultramafic intrusion in the Early Proterozoic Onega plateau, one of the oldest continental flood basalt provinces on Earth. The Sm–Nd and Pb–Pb combined mineral and whole-rock isochron ages of 1988±34 and 1985±57 Ma for the sill coincide with the age of ferropicrites from Pechenga (the Kola Peninsula). The lithostratigraphic, chemical and isotope evidence suggest the derivation of Pechenga lavas and the Onega plateau volcanics from a single mantle plume. Peridotite and gabbro whole-rock samples, and primary ulvospinel and ilmenite mineral separates from the sill yield a Re–Os isochron with a slope corresponding to an age of 1969±18 Ma, γ_Os(T) =−0.61±5.9. This age is consistent with the other isotope data, and indicates the closed-system behavior of Re and Os in the rocks. The peridotites and ulvospinel have high Os concentrations (2.5–14 ppb) and low ¹⁸⁷Re/¹⁸⁸Os ratios (0.35–1.1), thus allowing a more accurate determination of the weighted average initial ¹⁸⁷Os/¹⁸⁸Os of 0.1144±0.0019 (2σ_pop), γ_Os(T) =+0.77±1.7. This value is lower than that determined by Walker et al. (Geochim. Cosmochim. Acta 61 (1997) 3145–3160) for the Pechenga lavas (γ_Os(T) =+6.1±0.7), and implies a substantial Os-isotope heterogeneity in this ancient plume. Compared to the Onega plateau primary basalt magma, Pechenga ferropicrites are relatively enriched in iron and Ni, have lower (Nb/Th)_N ratios (2.1 vs 1.1) and less radiogenic Nd-isotope compositions (_Nd(T) = +3.1 and +1.4, respectively), but share similar low-radiogenic Pb-isotope characteristics (μ₁=8.57 and 8.60). Incorporation of small amounts (1.5%) of outer core material into the hotter central part of the plume and subsequent contamination of the Pechenga ferropicritic magmas with the 2.9 Ga Belomorian gneisses can explain the observed chemical and isotope variations in the two provinces provided that the core had <0.25 ppm of Pb.     

Dynamic properties and liquefaction resistance of two soil materials in an earthfill dam—Laboratory test results

Results are presented of laboratory resonant column and cyclic triaxial tests on specimens of two compacted soils (a sandy–silty clay and a sand–gravel mixture), planned to be used in the core and the shells, respectively, of a proposed earthfill dam. The values of low-amplitude shear modulus of the clayey material were found to increase with increasing confining pressure and decreasing water content, with deviations of ±20% from the predictions of the “Hardin equation”. On the other hand, the low-amplitude damping ratio was found to be at least four times higher than the values corresponding to natural undisturbed cohesive soils. The proposed G/G₀−γ_c curve for the compacted cohesive soil was found to be independent of confining pressure and small variations of the water content on either side of the optimum value and showed a remarkable agreement with recently published similar curves for natural cohesive soils. For the case of the sand–gravel mixture, normalized G/G₀–γ_c and D–γ_c curves are proposed, based on recently published results for gravelly soils and the limited data of the present study. The liquefaction resistance of the saturated sand–gravel mixture was found to be strongly dependent on its relative density, especially for high values of cyclic stress ratio. A relative density of at least 55% was found to be necessary to assure safety against earthquake-induced liquefaction of the material. The results presented herein may be used (directly or as guide) in the seismic analysis of (new or existing) earth dams constructed from similar soil materials and in addition they provide insight into the dynamic behavior of compacted soils.     

Dynamic response of a porous seabed around pipeline under three-dimensional wave loading

The evaluation of the wave-induced pore pressure around a buried pipeline is particularly important for pipeline engineers involved in the design of offshore pipelines. Most previous investigations of the wave-induced dynamic response around an offshore pipeline have limited to two-dimensional cases. In this paper, a three-dimensional model including buried pipeline is established, based on the existing DYNE3WAC models. Based on the proposed numerical model and poro-elastic soil material assumption, the effects of wave and soil characteristics, such as wave period, water depth, shear modulus and permeability, and configuration of pipelines, such as pipeline radius and pipeline buried depth, on the wave-induced excess pore pressure will be examined. Numerical results indicated that the normalized excess pore pressures versus z/h near the pipeline increase as the obliquity angle, wave period and water depth increase, and they decrease as the burial depth and radius of pipeline increase above the pipeline. Soil permeability has obvious influence on the wave-induced normalized excess pore pressure, and different soil material will result in distinct computation results.     

Determination of groundwater flow-through regimes of shallow lakes and wetlands from numerical analysis of stable isotope and chloride tracer distribution patterns

A numerical model previously developed to systematically examine groundwater flow in vertical section near shallow surface water bodies such as lakes, wetlands and ponds is further developed to include simulation of the distribution patterns of hydrogeochemical and stable isotopic tracers in relation to the surface water body and the geometry of distribution patterns of the tracers in the groundwater release zone of the lake. Many different possible flow regimes are identified, however, in this paper attention is focused on flow-through water bodies, as these are the flow regimes observed in field validation. Two shallow lakes on the Swan Coastal Plain of south-west Western Australia are the subject of field studies where hydrogeochemical and stable isotopes are used to validate the flow-through groundwater flow regime predicted by the modelling confirming the validity of the approach. The flow regime transition diagrams introduced in earlier papers are extended to include consideration of the hydrogeochemical and stable isotopic indices Cl_L/Cl₊ and (1000+δ_L)/(1000+δ₊). These ratios are introduced as an additional two of nine non-dimensional ratios that are necessary to analyse the problem. The ratios represent the chloride and isotopic composition (Cl_L and δ_L), respectively, in the groundwater release zones of the lakes, relative to these parameters in the groundwater capture zone (Cl₊ and δ₊) for the lake. Field data from the case studies plotted on appropriately configured transition diagrams demonstrates the overall validity of the modelling approach and its underlying assumptions. It is concluded that isotopic and hydrogeochemical data are invaluable in interpreting the interaction between lakes or wetlands and regional aquifers as it is very difficult to make physical or hydraulic measurements in the field that allow an understanding of lake–aquifer interaction. The tools and concepts developed that are summarized in the presented transition diagrams are invaluable starting points for the consideration and analysis of other case-specific examples of groundwater–surface water interaction and will improve the scientific basis of decision-making concerning lake and wetland management and groundwater interaction by water resource and environmental managers.     

Effects of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of pipeline buried in elasto-plastic seabed foundations

In this paper, a two-dimensional integrated numerical model is developed to examine the influences of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of a pipeline buried in a porous seabed. In the wave model, the RANS (Reynolds Averaged Navier–Stokes) equation is used to govern the wave motion. In the seabed model, the residual soil response in the vicinity of an embedded pipeline is considered with the 2-D elasto-plastic solution, where the phase-resolved shear stress is used as a source for the build-up of residual pore pressure. Classical Biot׳s consolidation equation is used for linking the solid-pore fluid interaction. The validation of the proposed integrated numerical model is conducted by the comparisons with the previous experimental data. Numerical examples show that the pore pressures can accumulate to a large value, thus resulting in a larger area of liquefaction potential in the given anisotropic soil compared to that with isotropic solution. The influences of anisotropic parameters on the wave-induced residual soil response in the vicinity of pipeline are significant. A high rate of pore pressure accumulation and dissipation is observed and the liquefaction potential develops faster as the anisotropic parameters increase. Finally, a simplified approximation based on a detailed parametric investigations is proposed for the evaluation of maximum liquefaction depth (z_L) in engineering application.     

Finite element modelling for water waves-soil interaction

The soil permeability and shear modulus of many marine sediments vary with depth because of consolidation under overburden pressure. However, conventional theories for wave-induced soil response have assumed a homogeneous porous seabed, with constant soil permeability and shear modulus. This paper presents a finite element model for the wave-induced soil response in a porous seabed, with variable permeability and shear modulus as a function of burial depth. The soil matrix considered here is unsaturated and hydraulically anisotropic, and subjected to a three-dimensional short-crested wave system. The present finite element formulation is established by using a combination of semi-analytical techniques and the Galerkin method. The nodal effective stresses directly derived from the governing equations can be calculated accurately in the present model. Verification is available through the reduction to the simple case of homogeneous seabed. Three typical marine materials, course, fine sand and gravel, are considered in this study. The numerical results indicate that the soil permeability affects the wave-induced seabed response significantly especially for gravelled seabed, as does the soil shear modulus for sandy seabed.     

Numerical analysis of liquefaction of porous seabed around pipeline fixed in space under seismic loading

Liquefaction of seabed under seismic loading is one of the main points that govern the overall stability of submarine pipeline. However, most previous investigations concerned only with free seabed and searched for seismic accumulative excess pore pressure by solving Terzaghi's consolidation equation containing pore pressure source term. It is not able to introduce two-dimensional structures such as submarine pipelines in one-dimensional problem, and it is also not able to obtain the distribution of seismic accumulative excess pore pressure in seabed around submarine pipelines by such a way. In this study, a FEM numerical analysis method for determining the liquefaction of sandy seabed around a buried pipeline under seismic loading is presented. The empirical mode of dynamic increase of pore pressure under undrained shearing induced by seismic loading is incorporated with two-dimensional dynamic consolidation equation and a numerical procedure based on FEM is developed to assess the accumulative excess pore pressure. By numerical computations, the accumulative process of pore pressure and liquefaction potential of seabed soil during seismic loading is evaluated. From a series of numerical computations based on the presented model with various parameters, the effects of soil characteristic parameters and pipeline geometry on seismic accumulative excess pore pressure around submarine pipeline and along the depth of seabed are explored in detail.     

Strain localisation in bimineralic rocks: Experimental deformation of synthetic calcite–anhydrite aggregates

Deformation of synthetic calcite–anhydrite aggregates to large shear strains (up to γ = 12.4 at 600 °C, 300 MPa confining pressure and a constant angular displacement rate corresponding to a shear strain rate of 10^− 3 s^− 1) resulted in the first experimental observation of strain localisation from initially homogeneous rocks. In contrast to experiments on pure calcite and anhydrite, which deformed homogeneously to large strains (γ ≥ 5), all experiments on calcite–anhydrite mixtures resulted in heterogeneous deformation at γ > 1 and the formation of narrow localised bands in the microstructures at γ > 4. In these bands, the amount of strain is at least twice as large as in the rest of the sample and individual grains of the same phase cluster and align, thereby forming microstructural layering similar to planar fabrics in natural mylonites. A switch in deformation mechanism in anhydrite from dislocation creep to diffusion creep and/or grain boundary sliding occurs simultaneously with strain localisation. It is concluded that deformation-induced heterogeneous phase distributions cause local strength differences initiating strain localisation in the calcite–anhydrite mixtures. The study suggests that the presence of two phases in combination with a change in deformation mechanism may be responsible for strain localisation in natural poly-mineralic mylonites.     

Apparent stress scaling and statistical trends

Statistical relations between different earthquake parameters, such as M₀ (seismic moment), E_S (seismic energy), τ_a (apparent stress), A (rupture area), g (average slip acceleration), are investigated. For this purpose, a kinematic earthquake model representing averaged earthquake rupture process is formulated. The model implies a scaling relationship for τ_a as a function of three other parameters, related to kinematic (M₀), geometric (A) and material (g) source characteristics, which, according to the model, can change independently. This scaling relation is used to explain statistical trends that characterize different earthquake data sets (including micro-, small, moderate and large events) plotted in the log τ_a − log M₀ space, and to determine the area in this space, where typical earthquakes occur. The scaling relationship is interpreted in terms of the apparent stress minimum (i.e., the most uniform among the possible earthquake rupture patterns). It is concluded that, although the apparent stress increases on an average with increasing seismic moment, small and large earthquakes are essentially similar.     

Response spectra for differential motion of columns paper II: Out-of-plane response

It is shown that the common response spectrum method for synchronous ground motion can be extended to make it applicable for earthquake response analyses of extended structures experiencing differential out-of-plane ground motion. A relative displacement spectrum for design of first-story columns SDC (T, T_T, ζ, ζ_T, τ, δ) is defined. In addition to the natural period of the out-of-plane response, T, and the corresponding fraction of critical damping, ζ, this spectrum also depends on the fundamental period of torsional vibrations, T_T, and the corresponding fraction of critical damping, ζ_T, on the “travel time,” τ (of the waves in the soil over a distance of about one-half the length of the structure), and on a dimensionless factor δ, describing the relative response of the first floor. The new spectrum, SDC, can be estimated by using the empirical scaling equations for relative displacement spectra, SD, and for peak ground velocity, v_max. For recorded strong-motion acceleration, and for symmetric buildings, the new spectrum can be computed from Duhamel's integrals of two uncoupled equations for dynamics equilibrium describing translation and rotation of a two-degree-of-freedom system. This representation is accurate when the energy of the strong-motion is carried by waves in the ground the wavelengths of which are one order of magnitude or more longer than the characteristic length of the structure.     

Rock magnetic and geochemical proxies for iron mineral diagenesis in a tropical lake: Lago Verde, Los Tuxtlas, East–Central Mexico

Magnetic and non-magnetic mineral analyses were conducted on a lacustrine sequence from Lago Verde in the tropical coast along the Gulf of Mexico that covers the last 2000 years. The site witnessed the transformation of the environment since the early Olmec societies until forest clearance in the last century. Through these analyses we investigated the processes that affected the magnetic mineralogy in order to construct a model of past environmental changes, and compare this model with the archeological record and inferred climatic changes in the northern hemisphere of tropical America. Volcanic activity has played a major influence on sediment magnetic properties, as a purveyor of Ti-magnetites/Ti-maghemites, and as a factor of instability in the environment. Anoxic reductive conditions are evident in most of Lago Verde's sedimentary record.

Direct observations of magnetic minerals and ratios of geochemical (Fe, Ti), and ferrimagnetic (χf) and paramagnetic (χp) susceptibility (χ) data, are used as parameters for magnetite dissolution (χp/χ, Fe/χf), and precipitation (χf/Ti) of magnetic minerals. Intense volcanic activity and anoxia are recorded before A.D. 20, leading to the formation of framboidal pyrite. Increased erosion, higher evaporation rates, lower lake levels, anoxia and reductive diagenesis in non-sulphidic conditions are inferred for laminated sediments between A.D. 20–850. This deposit matches the period of historical crisis and multiyear droughts that contributed to the collapse of the Maya civilization. Dissolution of magnetite, a high organic content and framboidal pyrite point to anoxic, sulphidic conditions and higher lake levels after A.D. 850. Higher lake levels in Lago Verde broadly coincide with the increased precipitation documented during the Medieval Warm Period (A.D. 950–1350) in the northern tropical and subtropical regions of the American continent. For the Little Ice Age (A.D. 1400–1800), the relatively moist conditions inferred are in concordance with the glacial advances recorded in central Mexico. Higher erosion rates reflect destruction of the rainforest over the last 40 years.     

A 53 year seasonally resolved oxygen and carbon isotope record from a modern Gibraltar speleothem: Reconstructed drip water and relationship to local precipitation

The response of a climate proxy against measured temperature, rainfall and atmospheric circulation patterns at sub-annual resolution is the ultimate test of proxy fidelity but very few data exist showing the level of correspondence between speleothem climate proxies and the instrumental climate record. Cave sites on the Gibraltar peninsula provide a unique opportunity to calibrate speleothem climate proxies with the longest known available precipitation isotopes and instrumental records. An actively growing speleothem sampled from New St. Michaels Cave in 2004 is composed of paired laminae consisting of light columnar calcite and a darker microsparitic calcite. Stable isotope analysis of samples micromilled in 100 μm steps at the equivalent of bi-monthly intervals reveals fabric-correlated annual cycles in carbon isotopes, oxygen isotopes and trace elements responding to seasonal changes in cave microclimate, hydrology and ventilation patterns. Calcite δ¹³C values reach a minimum in the light columnar fabric and evidence from trace element behaviour and cave monitoring indicates that this grows under cave ‘winter’ conditions of highest pCO₂, whereas the dark microsparitic calcite, characterised by elevated δ¹³C and δ¹⁸O values grows under low ‘summer’ pCO₂ conditions. Drip water δ¹³C_DIC reaches a minimum in March–April, at which time the attenuated δ¹⁸O signal becomes most representative of winter precipitation. An age model based on cycle counting and the position of the ¹⁴C bomb carbon spike yields a precisely dated winter oxygen isotope proxy of cave seepage water for comparison with the GNIP and instrumental climate record for Gibraltar. The δ¹⁸O characteristics of calcite deposited from drip water representing winter precipitation for each year can be derived from the seasonally resolved record and allows reconstruction of the δ¹⁸O drip water representing winter precipitation for each year from 1951–2004. These data show an encouraging level of correspondence (r² = 0.47) with the δ¹⁸O of rainfall falling each year between October and March and on a decadal scale the δ¹⁸O of reconstructed winter drip water mirrors secular change in mean winter temperatures.     

水域隧道地震响应分析

本文基于Biot动力固结理论和弹性动力学理论,考虑海床（土壤）的两相性、黏弹性人工边界及流（水）-固耦合作用,建立了隧道-土-流体相互作用的力学模型,讨论了P波作用下有无水的情况以及水深、水域隧道埋深、海床土性质和地震波入射角等因素对隧道及其周围海床应力的影响。结果表明:隧道周围海床土的孔隙水压力和隧道内应力随着水深的增加而增加;地震波特性和海床土特性对隧道的内应力和海床土的孔隙水压力均有较大的影响;海床土的渗透性和隧道埋深对隧道的内应力影响较小,而对隧道周围海床土的孔隙水压力影响较大;地震动的入射角对隧道的内应力和隧道附近土层的孔隙水压力均有较大影响。      

Oxygen isotopic composition of low-temperature authigenic clinoptilolite

Oxygen isotope ratios were obtained from authigenic clinoptilolites from Barbados Accretionary Complex, Yamato Basin, and Exmouth Plateau sediments (ODP Sites 672, 797, and 762) in order to investigate the isotopic fractionation between clinoptilolite and pore water at early diagenetic stages and low temperatures. Dehydrated clinoptilolites display isotopic ratios for the zeolite framework (δ¹⁸O_f) that extend from +18.7‰ to +32.8‰ (vs. SMOW). In combination with associated pore water isotope data, the oxygen isotopic fractionation between clinoptilolite and pore fluids could be assessed in the temperature range from 25°C to 40°C. The resulting fractionation factors of 1.032 at 25°C and 1.027 at 40°C are in good agreement with the theoretically determined oxygen isotope fractionation between clinoptilolite and water. Calculations of isotopic temperatures illustrate that clinoptilolite formation occurred at relatively low temperatures of 17°C to 29°C in Barbados Ridge sediments and at 33°C to 62°C in the Yamato Basin. These data support a low-temperature origin of clinoptilolite and contradict the assumption that elevated temperatures are the main controlling factor for authigenic clinoptilolite formation. Increasing clinoptilolite δ¹⁸O_f values with depth indicate that clinoptilolites which are now in the deeper parts of the zeolite-bearing intervals had either formed at lower temperatures (17–20°C) or under closed system conditions.     

Late Pleistocene to Holocene composite speleothem O and C chronologies from South Island, New Zealand—did a global Younger Dryas really exist?

Oxygen and carbon data from eight stalagmites from northwest South Island are combined to produce composite records of δ¹⁸O and δ¹³C from 23.4 ka to the present. The chronology is anchored by 43 thermal ionization mass spectrometry (TIMS) uranium series ages. Delta ¹⁸O values are interpreted as having a first order positive relationship to temperature, but also to be influenced by precipitation in a complex manner. Delta ¹³C is interpreted as responding negatively to increases in atmospheric CO₂ concentration, biological activity and precipitation amount.

Six climatic phases are recognized. After adjustment of 1.2‰ for the ice volume effect, the δ¹⁸O record between 23 and 18 ka varies around −3.72‰ compared to the Holocene average of −3.17‰. Late-glacial warming commenced between 18.2 and 17.8 ka and accelerated after 16.7 ka, culminating in a positive excursion between 14.70 and 13.53 ka. This was followed by a significant negative excursion between 13.53 and 11.14 ka of up to 0.55‰ depth that overlapped the Antarctic Cold Reversal (ACR) and spanned the Younger Dryas (YD). Positive δ¹⁸O excursions at 11.14 ka and 6.91–6.47 ka represent the warmest parts of the Holocene. The mid-Holocene from 6 to 2 ka was marked by negative excursions that coincide with increased glacial activity in the South Island. A short positive excursion from 0.71 to 0.57 ka was slightly later than the Medieval Warm Period of Europe.

Delta ¹³C values were high until 17.79 ka after which there was an abrupt decrease to 17.19 ka followed by a steady decline to a minimum at 10.97 ka. Then followed a general increase, suggesting a drying trend, to 3.23 ka followed by a further general decline. The abrupt decrease in δ-values after 17.79 ka probably corresponds to an increase in atmospheric CO₂ concentration, biological activity and wetness at the end of the Last Glaciation, but the reversal identified in the δ¹⁸O record from 13.53 to 11.14 ka was not reflected in δ¹³C changes. The lowest δ¹³C values coincided with the early Holocene climatic suboptimum when conditions were relatively wet as well as mild.

Major trends in the δ¹⁸O_c record are similar to the Northern Hemisphere, but second order detail is often distinctly different. Consequently, at the millennial scale, a more convincing case can be made for asymmetric climatic response between the two hemispheres rather than synchronicity.     

Preliminary stratigraphic and structural architecture of Bhutan: Implications for the along strike architecture of the Himalayan system

High-precision in-situ ion microprobe (SIMS) oxygen isotope analysis of zircons from two diorite intrusions associated with the late Caledonian Lochnagar pluton in Scotland has revealed large differences in the degree of heterogeneity in zircon δ¹⁸O between the diorites. Zircon crystals from the Cul nan Gad diorite (CnG) show a unimodal distribution of oxygen isotope values (δ¹⁸O = 6.0 ± 0.6‰ (2σ)) and no or only minor grain-scale variation. Those from the Allt Darrarie diorite (AD1) show a large range in δ¹⁸O and an apparent bimodal distribution with modes of 6.6 ± 0.4‰ and 7.3 ± 0.4‰. Variations of up to 1.2‰ occur between and within grains; both an increase and decrease in δ¹⁸O with zircon growth has been observed. The δ¹⁸O composition of growing zircon can only change if open-system processes affect the magma composition, i.e. if material of contrasting δ¹⁸O composition is added to the magma. The variability in AD1 is interpreted to represent a cryptic record of magma mixing. A ‘deep crustal hot zone’ is a likely site for generation of the dioritic magmas which developed by mixing of residual melts and crustal partial melts or by melting of mafic lower crustal rocks. The overall small number of zircons with mantle-like δ¹⁸O values (5.3 ± 0.6‰ (2σ)) in the Lochnagar diorites is largely the product of crustal differentiation rather than crustal growth.

The δ¹⁸O of quartz from the CnG and AD1 diorites shows only minor variation (CnG: 10.9 ± 0.5‰ (2σ), AD1: 11.7 ± 0.6‰ (2σ)) within single populations, with no evidence of mixing. Quartz–zircon isotopic disequilibrium is consistent with later crystallisation of quartz from late magmatic fluids, and in case of the AD1 diorite after the inferred magma mixing from a homogenised, higher δ¹⁸O melt.

High-precision SIMS oxygen isotope analysis of zircon provides a new approach to identifying and resolving previously undetected early-stage magma mixing and constraining the compositions and origins of the component magmas. A combination of zircon, quartz and whole-rock data has proven to be a powerful tool in reconstructing the petrogenetic evolution of diorite from early crystallisation to late alteration.     

Silicate perovskite analogue ScAlO3: temperature dependence of elastic moduli

The elastic moduli of ScAlO₃ perovskite, a very close structural analogue for MgSiO₃ perovskite, have been measured between 300 and 600 K using high precision ultrasonic interferometry in an internally heated gas-charged pressure vessel. This new capability for high temperature measurement of elastic wave speeds has been demonstrated on polycrystalline alumina. The temperature derivatives of elastic moduli of Al₂O₃ measured in this study agree within 15% with expectations based on published single-crystal data. For ScAlO₃ perovskite, the value of (∂K_S/∂T)_P is −0.033 GPa K⁻¹ and (∂G/∂T)_P is −0.015 GPa K⁻¹. The relative magnitudes of these derivatives agree with the observation in Duffy and Anderson [Duffy, T.S., Anderson, D.L., 1989. Seismic velocities in mantle minerals and the mineralogy of the upper mantle. J. Geophys. Res. 94, 1895–1912.] that |(∂K_S/∂T)_P| is typically about twice |(∂G/∂T)_P|. The value of (∂K_S/∂T)_P for ScAlO₃ is intermediate between those inferred less directly from V(P,T) studies of Fe-free and Fe- and Al-bearing MgSiO₃ perovskites [Wang, Y., Weidner, D.J., Liebermann, R.C., Zhao, Y., 1994. P–V–T equation of state of (Mg,Fe)SiO₃ perovskite: constraints on composition of the lower mantle. Phys. Earth Planet. Inter. 83, 13–40; Mao, H.K., Hemley, R.J., Shu, J., Chen, L., Jephcoat, A.P., Wu, Y., Bassett, W.A., 1991. Effect of pressure, temperature and composition on the lattice parameters and density of (Mg,Fe) SiO₃ perovskite to 30 GPa. J. Geophys. Res. 91, 8069–8079; Zhang, Weidner, D., 1999. Thermal equation of state of aluminum-enriched silicate perovskite. Science 284, 782–784]. The value of |(∂G/∂T)|_P for ScAlO₃ is similar to those of most other mantle silicate phases but lower than the recent determination for MgSiO₃ perovskite [Sinelnikov, Y., Chen, G., Neuville, D.R., Vaughan, M.T., Liebermann, R.C., 1998. Ultrasonic shear wave velocities of MgSiO₃ perovskite at 8 GPa and 800K and lower mantle composition. Science 281, 677–679].

Combining the results from the previous studies and current measurements on ScAlO₃ perovskite, we extracted the parameters (q and γ₀) needed to fully specify its Mie–Grüneisen–Debye equation-of-state. In this study, we have demonstrated that acoustic measurements of K_S(T), unlike V(P,T) data, tightly constrain the value of q. It is concluded that ScAlO₃ has ‘normal’ γ₀ (1.3) and high q (3.6). The high value of q indicates that ScAlO₃ has very strong intrinsic temperature dependence of the bulk modulus; similar behaviour has been observed in measurements on Fe- and Al-bearing silicate perovskites (Mao et al., 1991; Zhang and Weidner, 1999).     

Pipe–soil interaction for segmented buried pipelines subjected to dip faults

This paper describes an investigation of pipe–soil interaction equations suggested by currently used pipeline seismic design codes and the applicability of these equations to segmented pipelines. The results of computer‐aided analyses were compared to results obtained in full‐scale experiments on a segmented ductile iron pipeline 93 mm in diameter and 15 m in length. The pipeline was installed 600 mm below the ground surface in a sandy soil compacted to two different subgrade reaction values. The type of fault considered was a reverse fault with an intersection angle of 60° with the pipeline, and the fault movement was a total of 350 mm in three same steps in the fault trace direction. The findings of this study demonstrate the necessity of considering the nature of soil behavior in pipe–soil interaction equations and the effects of connection joints on the integrated response of pipelines to fault‐induced ground deformations. A new combination of equations constituting a direction‐wise selection from among the equations proposed by currently used guidelines is introduced as a new series to describe pipe–soil interaction for segmented pipelines and is verified using the results of full‐scale experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

A critical evaluation of young (near-zero) K–Ar ages

An evaluation of the precision and resolution of the unspiked K–Ar dating method is presented with particular regard to the statistical significance of ages that are measured near or at the detection limit of the technique. Near-zero (historical) ages can be measured by the unspiked K–Ar technique with a precision that is essentially controlled by the precision with which the ⁴⁰Ar/³⁶Ar of the sample can be resolved from the present-day atmospheric value of 295.5. The best analytical precision on the isotopic ratio is ±0.05% (1σ) by this technique, which currently limits the lower detection limit of unspiked K–Ar ages to samples featuring at least 0.14% of radiogenic ⁴⁰Ar. The corresponding youngest resolvable K–Ar age depends on the K content and atmospheric contamination of the sample. Total-fusion analysis of high-K refractory minerals like sanidine is not practicable via K–Ar, and the lowest resolvable age for medium-K samples more amenable to complete fusion is around 1.5 ka (on a single-run basis). It is argued that near-zero age measured with a probability density straddling or narrowing the time-origin cannot be handled without accounting for the non-negativity constraint imposed by the physical requirement of a positive age. The pertinent equations are derived both for the single-run case and for the case of independent replicates made on a single sample. We show that pooled K–Ar replicates can theoretically reduce the nominal uncertainty of individual unspiked ages (typically ±1.5 ka, 2σ) to a value that is close to the smallest ⁴⁰Ar/³⁹Ar isochron age uncertainty achievable on sanidine in the 0–2 ka range (±0.2 ka, 2σ). However, this performance is obtained at the cost of prohibitively large-sample statistics (n≥15) for medium-K feldspars datable via K–Ar. Coupled with the inability of the K–Ar approach to obviate the problems of excess/fractionated ⁴⁰Ar and/or xenocrystic contamination, this makes the ⁴⁰Ar/³⁹Ar technique the method of choice for dating historical events by the K–Ar scheme.