Modeling boulder transport by coastal waves on cliff topography: Case study at Hachijo Island,Japan

Clifftop coastal boulders transported by storm waves or tsunamis have been reported around the world. Although numerical calculation of boulder transport is a strong tool for the identification of tsunami or storm boulders, and for estimation of the wave size emplacing boulders, models which can reasonably solve boulder transport from below a cliff or from a cliff-edge onto a cliff-top do not yet exist. In this study, we developed a new numerical formulation for cliff-top deposition of boulders from the cliff edge or below the cliff, with validation from laboratory tests. We then applied the model using storm and tsunami wave forcing to simulate the observed boulder deposits at the northwest coast of Hachijo Island, Japan. Using the model, the actual distribution of boulders was explained well using a reasonable storm wave height without assumption of anomalously high-water level by storm surge. Results show that boulder transport from the cliff edge or under the cliff onto the cliff-top was possible from a tsunami with periods of 5~10 min or storm waves with no storm surge. However, the actual distribution of boulders on the cliff was explained only by storm waves, but not by tsunami. Therefore, the boulders distributed at this site are likely of storm wave origin. Our developed model for the boulder transport calculation can be useful for identifying a boulder's origin and can reasonably calculate cliff-top deposition of boulders by tsunami and storm waves. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Physical properties of the shallow sediments in late Pleistocene formations,Ursa Basin,Gulf of Mexico,and their implications for generation and preservation of shallow overpressures

Understanding the evolution of abnormally high fluid pressures within sedimentary formations is critical for analysing hydrogeological processes and assessing drilling risks. We have constructed a two-dimensional basin model and have performed numerical simulations to increase the understanding of the history of fluid flow and shallow overpressures in the Pleistocene and Holocene formations in the Ursa basin, deepwater Gulf of Mexico. We measured physical properties of sediments, such as porosity and permeability, in the laboratory and estimated in situ pore pressures from preconsolidation pressures. We obtained porosity–effective stress relationships from measurements of bulk density, grain density and preconsolidation pressures in the laboratory. Porosity–effective stress relationships were also obtained from downhole density logs and measured pore pressures. The porosity–effective stress and porosity–permeability relationships obtained were applied in two-dimensional basin simulations. Results showed that high pore pressures developed shortly after sediment deposition. Peaks in pore pressure ratios were related to high sedimentation rates of mass transport deposits and the incision of the Ursa channel. Lateral flows from the area where the overburden is thick towards the area where it is thin have occurred at least since 30 ka. Present pore pressure and temperature distributions suggest that lateral flows play a role in re-distributing heat in the basin.     

Accurate ocean tide modeling in southeast Alaska and large tidal dissipation around Glacier Bay

An accurate prediction of ocean tides in southeast Alaska is developed using a regional, barotropic ocean model with a finite difference scheme. The model skill is verified by the observational tidal harmonics in southeast Alaska including Glacier Bay. The result is particularly improved in Glacier Bay compared to the previous model described by Foreman et al. (2000). The model bathymetry dominates the model skill. We re-estimate tidal energy dissipation in the Alaska Panhandle and suggest a value for tidal energy dissipation of 3.4 GW associated with the M₂ constituent which is 1.5 times the estimation of Foreman et al. (2000). A large portion of the M₂ energy budget entering through Chatham Strait is dissipated in the vicinity of Glacier Bay. Moreover, it is shown that the developed model has the potential to correct the ocean tide loading effect in geodetic data more efficiently than the model of Foreman et al. (2000), especially around Glacier Bay.     

Wind- and density-driven circulation in a bay on the Sanriku ria coast,Japan: study of Shizugawa Bay facing the Pacific Ocean

Flow fields in Shizugawa Bay on the Sanriku ria coast, which faces the Pacific Ocean, were investigated using hydrographic observations for the purpose of understanding oceanographic conditions and the process of water exchanges in the bay after the 2011 earthquake off the Pacific coast of Tohoku. In spring to summer, density-driven surface outflow is part of estuarine circulation and is induced by a pressure gradient force under larger longitudinal gradients in density along with lower salinity water in the innermost part of the bay, regardless of wind forcing. In winter to summer, another density-driven current with a thermal structure is induced by a pressure gradient force under the smaller longitudinal density gradients in calm wind conditions. Particularly in winter, Tsugaru Warm Current water can be transported in the surface layer inside the bay. Wind-driven bay-scale circulation with downwind and upwind currents in the surface and deeper layers, respectively, is induced by strong longitudinal wind forcing under the smaller longitudinal density gradients, irrespective of season. Particularly in fall to spring, this circulation can cause the intrusions of oceanic water associated with Oyashio water and Tsugaru Warm Current water in the deeper layer. These results suggest that wind- and density-driven currents can produce the active exchange of water from inside and outside the bay throughout the year.     

东海风应力旋度驱动的21天周期的海底压力变化及其传播

在2015年6月至2017年6月期间,跨越庆良间水道的由2台加载压力传感器的倒置式回声仪（PIES）和5台加载压力传感器和海流计的倒置式回声仪（CPIES）组成的观测断面获得了近2年的海底压力时间序列。该时间序列中存在着显著的21天周期的振荡（P_bot21）,该信号在2016年7月至10月期间尤为强烈。P_bot21与东海陆架上的风应力旋度存在较显著的3天延迟相关,其相关系数达到0.65。本文采用正压海洋模式解释了这一信号的产生、传播以及耗散过程,模式结果显示东海陆架上的风应力旋度驱动产生P_bot21并向琉球岛链传播,而深海上的风应力旋度不能驱动产生这一信号。在陆架上,P_bot21伴随21天周期的风应力旋度由海岸向东南方向传播,但由于摩擦作用,信号在离开风场后几天内即耗散。断面能否观测到P_bot21与陆架上21天周期风应力旋度场的分布相关,长江口东南方向风应力旋度驱动的P_bot21能被观测到,而长江口东北方向产生的P_bot21不能被观测到。     

Fate of carbonates within oceanic plates subducted to the lower mantle, and a possible mechanism of diamond formation

We report on high-pressure and high-temperature experiments involving carbonates and silicates at 30–80 GPa and 1,600–3,200 K, corresponding to depths within the Earth of approximately 800–2,200 km. The experiments are intended to represent the decomposition process of carbonates contained within oceanic plates subducted into the lower mantle. In basaltic composition, CaCO₃ (calcite and aragonite), the major carbonate phase in marine sediments, is altered into MgCO₃ (magnesite) via reactions with Mg-bearing silicates under conditions that are 200–300°C colder than the mantle geotherm. With increasing temperature and pressure, the magnesite decomposes into an assemblage of CO₂ + perovskite via reactions with SiO₂. Magnesite is not the only host phase for subducted carbon—solid CO₂ also carries carbon in the lower mantle. Furthermore, CO₂ itself breaks down to diamond and oxygen under geotherm conditions over 70 GPa, which might imply a possible mechanism for diamond formation in the lower mantle.     

Flux‐Free Fusion of Silicate Rock Preceding Acid Digestion for ICP‐MS Bulk Analysis

We present a new method for the decomposition of silicate rocks by flux‐free fusion in preparation for whole‐rock trace element determination (Sc, Rb, Sr, Y, Zr, Nb, Cs, Ba, rare earth elements and Hf) that is especially applicable to zircon‐bearing felsic rocks. The method was verified by analyses of RMs of mafic (JB‐1a, JB‐2, JGb‐1) and felsic rocks (JG‐3, JR‐3, JSd‐1, GSP‐2, G‐2). Pellets of powdered sample (up to 500 mg) without flux were weighed and placed in a clean platinum crucible. The samples were then fused in a Siliconit® tube furnace and quenched to room temperature. The optimum condition for the fusion of granitic rock was determined to be heating for 2–3 min at 1600 °C. The fused glass in the platinum crucible after heating was decomposed using HF and HClO₄ in a Teflon® beaker. Decomposed and diluted sample solutions were analysed using a quadrupole inductively coupled plasma‐mass spectrometer. Replicate analyses (n = 4 or 5) of the RMs revealed that analytical uncertainties were generally < 3% for all elements except Zr and Hf (～ 6%) in JG‐3. These higher uncertainties may be attributed to sample heterogeneity. Our analytical results for the RMs agreed well with recommended concentrations and recently published concentrations, indicating complete decomposition of our rock samples during fusion.     

Recycled crustal melt injection into lithospheric mantle: implication from cumulative composite and pyroxenite xenoliths

A rare composite xenolith and abundant cumulative pyroxenites obtained from the Mesozoic Fangcheng basalts on the eastern North China Craton record a complex history of melt percolation and circulation in the subcontinental lithospheric mantle. The composite xenolith has a dunite core and an olivine clinopyroxenite rim. The dunite is of cumulative origin and has a granular recrystallized texture and extremely low Mg# [100 Mg/(Mg + Fe) = 81–82] contents in olivines. The olivine clinopyroxenite contains larger clinopyroxene and/or orthopyroxene with a few fine-grained olivine and tiny phlogopite, feldspar, and/or carbonate minerals interstitial to clinopyroxene. The clinopyroxene has low Mg# (83–85). Compositional similarity between dunitic olivine and pyroxenitic one indicates a sequential crystallization of dunite and pyroxenite from a precursor melt. Pyroxenite xenoliths include olivine websterites and clinopyroxenites, both are of cumulative origin. Estimation of the melt from major oxides in olivines and REE concentrations in clinopyroxenes in these composite and pyroxenite xenoliths suggests a derivation from subducted crustal materials, consistent with the highly enriched EMII-like Sr and Nd isotopic ratios observed in the pyroxenites. Occurrence of phlogopite, feldspar and carbonate minerals in some xenoliths requires the melt rich in alkalis (K, Na), silica and volatiles (water and CO₂) at the latest stage as well, similar to highly silicic and potassic melts. Thus, the occurrence of these composite and pyroxenite xenoliths provides an evidence for voluminous injection of recycled crustal melts into the lithosphere beneath the southeastern North China Craton at the Late Mesozoic, a reason for the rapid lithospheric enrichment in both elemental and isotopic compositions.     

Trace-element characteristics of east–west mantle geochemical hemispheres

Recent statistical analyses on the isotopic compositions of oceanic, arc, and continental basalts have revealed that the Earth's mantle is broadly divided into eastern and western hemispheres. The present study aimed to characterize the isotopically defined east–west geochemical hemispheres using trace-element concentrations. Basalt data with Rb, Sr, Nd, Sm, Pb, Th, and U in addition to the isotopic ratios ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb were selected mostly from the GEOROC and PetDB databases. A total of 4787 samples were used to investigate the global geochemical variations. The results show that the wide trace-element variations are broadly explained by the melting of melt-metasomatized and fluid-metasomatized mantle sources. The larger amount of the fluid component derived from subducted plates in the eastern hemisphere than that in the western hemisphere is inferred from the basalts. These characteristics support the hypothesis that focused subduction towards the supercontinent created the mantle geochemical hemispheres.