Paleolimnological investigation of recent acidity changes in Sawanoike Pond, Kyoto, Japan

Sawanoike Pond, an artificial reservoir constructed ca. 150 yrs. ago, lies in an area geologically sensitive to acidification and contains water that is presently weakly acidic (pH 5.5; alkalinity 15 eq 1^-1). Sediments from a core obtained from the northern part and a core from the central part of the pond were analyzed for organic matter, total organic carbon, total nitrogen, natural and artificial radionuclides, pollen, diatoms, metals and spheroidal carbonaceous particles. A peaty silt layer at the bottom of the central core was replaced by a silt layer initially rich in humic substances, but became more mineral-rich near the surface. Although most of the dominant diatoms throughout the cores were acidophilous taxa, some acidobiontic species increased moderately in abundance after the middle 1960's. Diatom-inferred pH values for the northern core are relatively low before about 1930, then increase, reaching a maximum level between the 1930's and the late 1940's. They then decrease to the surface and reach minimum level during the middle 1980's. Diatom-inferred pH declined by 0.4-0.5 units in the northern core and 0.2-0.3 units in the central core since ca. 1910, when vegetation around the shore changed to a community similar to that of recent years. An increase in labile Pb, Zn, Cd and spheroidal carbonaceous particles, indicative of atmospheric deposition of fossil fuel combustion by-products, preceded the decline in inferred pH. The inferred pH decline of the pond could be an acidification due to acidic precipitation caused by fossil fuel combustion.     

日本南海海槽俯冲增生楔前缘的构造变形特征

对增生楔不同压力—温度条件下的构造变形、流体活动、沉积特征、岩石物性和化学组成等多方面的直接观测,可以帮助分析俯冲带地震的蕴育和发生的环境与机理。通过参加IODP的日本南海海槽发震带研究项目(NanTroSEIZE)第一阶段316航次所收集到的大量第一手数据和资料,分别在4个站位上(C0004,C0006,C0007,C0008)对日本南海海槽增生楔前缘岩芯尺度上的构造变形进行了详细分析,并且讨论了岩芯尺度上的构造变形与增生楔中大尺度的非序列分支逆冲断层和前缘逆冲断层的构造变形之间的关系。发现逆冲变形不是只在大尺度的逆冲断层面上进行,而是弥散分布在主逆冲断层面、次级逆冲断层面以及断层面之间的更小的尺度上。小尺度构造的倾向与大尺度断层的倾向有较好的一致性,表明它们是在相同的应力场下所形成的。在增生楔浅部高角度的正断层比较发育,显示张性应力场特征,同时所获得的岩芯尺度上的地层倾角较大并倾向与反射地震以及区域地质分析结果非常吻合,而在深部,特别是在大尺度逆冲断层发育带附近,各种类型的断层、滑移变形带、节理等非常普遍,同时层理与劈理的产状的复杂变化更多地受控于复杂的逆冲断层带的作用。     

Spectroscopic properties of Cr-doped and Cr,Li-doped synthetic forsterite crystals

Cr-droped and Cr,Li-doped forsterite crystals were grown and their optical properties were investigated. It was shown that when only Cr is doped, Cr³⁺ is substituted at the site of low crystal field, and the energy level ²E lie above the ⁴T₂ level, while ⁴T₂ is just above ²E when Cr and Li are codoped. The difference was rationalized by a deformation of the Cr substituted site with the introduction of Li.     

Advances in computational morphodynamics using the International River Interface Cooperative (iRIC) software

Results from computational morphodynamics modeling of coupled flow–bed–sediment systems are described for 10 applications as a review of recent advances in the field. Each of these applications is drawn from solvers included in the public-domain International River Interface Cooperative (iRIC) software package. For mesoscale river features such as bars, predictions of alternate and higher mode river bars are shown for flows with equilibrium sediment supply and for a single case of oversupplied sediment. For microscale bed features such as bedforms, computational results are shown for the development and evolution of two-dimensional bedforms using a simple closure-based two-dimensional model, for two- and three-dimensional ripples and dunes using a three-dimensional large-eddy simulation flow model coupled to a physics-based particle transport model, and for the development of bed streaks using a three-dimensional unsteady Reynolds-averaged Navier–Stokes solver with a simple sediment-transport treatment. Finally, macroscale or channel evolution treatments are used to examine the temporal development of meandering channels, a failure model for cantilevered banks, the effect of bank vegetation on channel width, the development of channel networks in tidal systems, and the evolution of bedrock channels. In all examples, computational morphodynamics results from iRIC solvers compare well to observations of natural bed morphology. For each of the three scales investigated here, brief suggestions for future work and potential research directions are offered. © 2019 The Authors Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd     

The variety of lithologies in the Yamato-86032 lunar meteorite: Implications for formation processes of the lunar crust

We performed a petrologic, mineralogical, geochemical, and isotopic study of several lithologies in the Y-86032 feldspathic breccia. This study leads us to conclude that Y-86032 likely originated on the lunar farside. Y-86032 is composed of several types of feldspathic clasts, granulitic breccias, and minor basaltic clasts set in a clastic matrix. We identify an “An97 anorthosite” that has An contents similar to those of nearside FANs. Mg′ (= molar Mg/(Mg + Fe) × 100) values vary significantly from ∼45 to ∼80 covering the ranges of both nearside FANs and the Mg′ gap between FANs and the Mg-suite. A light-gray feldspathic (LG) breccia making up ∼20% of the investigated slab (5.2 × 3.6 cm²) mainly consists of fragments of anorthosites (“An93 anorthosite”) more sodic than nearside FANs. LG also contains an augite-plagioclase clast which either could be genetically related to the An93 anorthosite or to slowly-cooled basaltic magma intruded into the precursor rock. The Na-rich nature of both An93 anorthosite and this clast indicates that the LG breccia was derived from a relatively Na-rich but incompatible-element-poor source. The Mg′ variation indicates that the “An97 anorthosite” is a genomict breccia of several types of primary anorthosites. Granulitic breccias in Y-86032 have relatively high Mg′ in mafic minerals. The highest Mg′ values in mafic minerals for the “An97 anorthosite” and granulitic breccias are similar to those of Mg-rich lithologies recently described in Dhofar 489. Basaltic clasts in the dark-gray matrix are aluminous, and the zoning trends of pyroxene are similar to those of VLT or LT basalts. The crystallization of these basaltic clasts pre-date the lithification age of the clastic matrix at ∼3.8 Ga. The low K contents of plagioclase in both the anorthositic and basaltic clasts and generally low incompatible element abundances in all the lithologies in Y-86032 indicate that KREEP was not involved during the formation of the precursor lithologies. This observation further suggests that urKREEP did not exist in the source regions of these igneous lithologies. All these facts support the idea that Y-86032 was derived from a region far distant from the PKT and that the lithic clasts and fragments are indigenous to that region. An An97 anorthositic clast studied here has distinct Sm-Nd isotopic systematics from those previously found for another An97 anorthositic clast and “An93 anorthosite”, and suggests either that An97 anorthosites come from isotopically diverse sources, or that the Sm-Nd isotopic systematics of this clast were reset ∼4.3 Ga ago. These lines of geochemical, isotopic, and petrologic evidence suggest that the lunar crust is geochemically more heterogeneous than previously thought.     

Summer Water Balance Characteristics of Koryto Glacier, Kamchatka Peninsula, Russia

The daily water balance for the drainage basin of Koryto Glacier, Kamchatka Peninsula, Russia, was calculated during the period from August to September 2000. The result shows that 14×10⁶ m³ of meltwater and 2×10⁶ m³ of rainwater entered the basin, while 26×10⁶ m³ of water drained from the basin through proglacial streams. Thus, about ?9×10⁶ m³ of water storage reduction occurred in the basin. Vertical displacements of the glacier surface showed that the volume change due to contraction of subglacial cavities was nearly 20% of the total storage change. The remaining fraction of water storage during the period is thought to be stored in englacial and supraglacial locations. The estimate of water balance components in the early ablation season in 2000 indicates that meltwater was already stored within the glacier before the spring, even during the previous year, and that the stored water drained through the ablation season.     

Possible fluid-rock interactions on differentiated asteroids recorded in eucritic meteorites

The eucritic meteorites are basaltic rocks that originate from the upper part of the crust of some small bodies as exemplified possibly by asteroid 4-Vesta. A few eucrites appear to have been modified by different degrees of a late stage alteration process that caused significant variations in mineralogy. Three distinct alteration stages are identified: (1) Fe-enrichment along the cracks that cross cut the pyroxene crystals (“Fe-metasomatism”); secondary olivine and minute amounts of troilite are found only occasionally in cracks at this stage; (2) deposits of Fe-rich olivine (Fa_64-86) and minor amounts of troilite are frequent inside the cracks; sporadic secondary Ca-rich plagioclase (An_97-98) is associated with the fayalitic olivine; (3) at this stage, the Fe-enrichment of the pyroxene is accompanied by a marked Al-depletion; moreover, secondary Ca-rich plagioclase is more frequent and partly fills some cracks or rims of the primary plagioclase crystals. The composition of the secondary phases on one hand, the lack of incompatible trace element enrichment in the metasomatized pyroxenes on the other hand, rule out a silicate melt as the metasomatic agent. Although no hydrous phase has been yet identified in the studied samples, aqueous fluids are plausible candidates for explaining the deposits of ferroan olivine and anorthitic plagioclase inside the fractures of the studied unequilibrated eucrites.     

Tasks and problems involved in the handling of disaster waste upon April 2016 Kumamoto Earthquake,Japan

The potential impacts of tsunamis along the Catalan Coast (NW Mediterranean) are analysed using numerical modelling. The region is characterized by moderate to low seismic activity and by moderate- to low-magnitude earthquakes. However, the occurrence of historical strong earthquakes and the location of several active offshore faults in front of the coast suggest that the possibility of an earthquake-triggered tsunami is not negligible although of low probability. Up to five faults have been identified to generate tsunamis, being the highest associated possible seismic magnitudes of up to 7.6. Coastal flooding and port agitation are characterized using the Worst-case Credible Tsunami Scenario Analysis approach. The results show a multiple fault source contribution to tsunami hazard. The shelf dimensions and the existence of submerged canyons control the tsunami propagation. In wide shelves, waves travelling offshore may become trapped by refraction causing the wave energy to reach the coastline at some distance from the origin. The free surface water elevation increases at the head of the canyons due to the sharp depth gradients. The effects of potential tsunamis would be very harmful in low-lying coastal stretches, such as deltas, with a high population concentration, assets and infrastructures. The Ebro delta appears to be the most exposed coast, and about the 20% of the delta surface is prone to flooding due to its extremely low-lying nature. The activity at Barcelona port will be severely affected by inflow backflow current at the entrance of up to 2 m/s.     

Development of regional simulation of seismic ground‐motion and induced liquefaction enhanced by GPU computing

Much research has been conducted for physics‐based ground‐motion simulation to reproduce seismic response of soil and structures precisely and to mitigate damages caused by earthquakes. We aimed at enabling physics‐based ground‐motion simulations of complex three‐dimensional (3D) models with multiple materials, such as a digital twin (high‐fidelity 3D model of the physical world that is constructed in cyberspace). To perform one case of such simulation requires high computational cost and it is necessary to perform a number of simulations for the estimation of parameters or consideration of the uncertainty of underground soil structure data. To overcome this problem, we proposed a fast simulation method using graphics processing unit computing that enables a simulation with small computational resources. We developed a finite‐element‐based method for large‐scale 3D seismic response analysis with small programming effort and high maintainability by using OpenACC, a directive‐based parallel programming model. A lower precision variable format was introduced to achieve further speeding up of the simulation. For an example usage of the developed method, we applied the developed method to soil liquefaction analysis and conducted two sets of simulations that compared the effect of countermeasures against soil liquefaction: grid‐form ground improvement to strengthen the earthquake resistance of existing houses and replacement of liquefiable backfill soil of river wharves for seismic reinforcement of the wharf structure. The developed method accelerates the simulation and enables us to quantitatively estimate the effect of countermeasures using the high‐fidelity 3D soil‐structure models on a small cluster of computers.