Curvature derived from LiDAR digital elevation models as simple indicators of debris-flow susceptibility

To mitigate the damage caused by debris flows resulting from heavy precipitation and to aid in evacuation plan preparation, areas at risk should be mapped on a scale appropriate for affected individuals and communities. We tested the effectiveness of simply identifying debris-flow hazards through automated derivation of surface curvatures using LiDAR digital elevation models. We achieved useful correspondence between plan curvatures and areas of existing debris-flow damage in two localities in Japan using the analysis of digital elevation models(DEMs). We found that plan curvatures derived from 10 m DEMs may be useful to indicate areas that are susceptible to debris flow in mountainous areas. In residential areas located on gentle sloping debris flow fans, the greatest damage to houses was found to be located in the elongated depressions that are connected to mountain stream valleys. Plan curvaturederived from 5 m DEM was the most sensitive indicators for susceptibility to debris flows.     

Evolution processes of Ordovician–Devonian arc system in the South‐Kitakami Massif and its relevance to the Ordovician ophiolite pulse

The South Kitakami Massif is one of the oldest geological domains in Japan having Silurian strata with acidic pyroclastic rocks and Ordovician–Silurian granodiorite–tonalite basement, suggesting that it was matured enough to develop acidic volcanisms in the Silurian period. On the northern and western margin of the South Kitakami Massif, an Ordovician arc ophiolite (Hayachine–Miyamori Ophiolite) and high‐pressure and low‐temperature metamorphic rocks (Motai metamorphic rocks) exhumed sometime in the Ordovician–Devonian periods are distributed. Chronological, geological, and petrochemical studies on the Hayachine–Miyamori Ophiolite, Motai metamorphic rocks, and other early Paleozoic geological units of the South Kitakami Massif are reviewed for reconstruction of the South Kitakami arc system during Ordovician to Devonian times with supplementary new data. The reconstruction suggests a change in the convergence polarity from eastward‐ to westward‐dipping subduction sometime before the Late Devonian period. The Hayachine–Miyamori Ophiolite was developed above the eastward‐dipping subduction through three distinctive stages. Two separate stages of overriding plate extension inducing decompressional melting with minor involvement of slab‐derived fluid occurred before and after a stage of melting under strong influence of slab‐derived fluids. The first overriding plate extension took place in the back‐arc side forming a back‐arc basin. The second one took place immediately before the ophiolite exhumation and near the fore‐arc region. We postulate that the second decompressional melting was triggered by slab breakoff, which was preceded by slab rollback inducing trench‐parallel wedge mantle flow and non‐steady fluid and heat transport leaving exceptionally hydrous residual mantle. The formation history of the Hayachine–Miyamori Ophiolite implies that weaker plate coupling may provide preferential conditions for exhumation of very hydrous mantle. Very hydrous peridotites involved in arc magmatism have not yet been discovered except for in the Cambrian–Ordovician periods, suggesting its implications for global geodynamics, such as the thermal state and water circulation in the mantle.     

Interdiffusion of Mg–Fe in olivine at 1,400–1,600 °C and 1 atm total pressure

The interdiffusion coefficient of Mg–Fe in olivine (D _Mg–Fe) was obtained at 1,400–1,600 °C at the atmospheric pressure with the oxygen fugacity of 10^?3.5–10^?2 Pa using a diffusion couple technique. The D _Mg–Fe shows the anisotropy (largest along the [001] direction and smallest along the [100] direction), and its activation energy (280–320 kJ/mol) is ~80–120 kJ/mol higher than that estimated at lower temperatures. The D _Mg–Fe at temperatures of >1,400 °C can be explained by the cation-vacancy chemistry determined both by the Fe³⁺/Fe²⁺ equilibrium and by the intrinsic point defect formation with the formation enthalpy of 220–270 kJ/mol depending on the thermodynamical model for the Fe³⁺/Fe²⁺ equilibrium in olivine. The formation enthalpy of 220–270 kJ/mol for the point defect (cation vacancy) in olivine is consistent with that estimated from the Mg self-diffusion in Fe-free forsterite. The increase in the activation energy of D _Mg–Fe at >1,400 °C is thus interpreted as the result of the transition of diffusion mechanism from the transition metal extrinsic domain to the intrinsic domain at the atmospheric pressure.     

Melting and Melt Segregation in the Mantle Wedge above a Subduction Zone: Evidence from the Chromite-bearing Peridotites of the Miyamori Ophiolite Complex, Northeastern Japan

Chromite-bearing peridotites of the Ordovician Miyamori ophiolitecomplex exhibit spatial mineralogical variations on scales rangingfrom several centimeters to a few kilometers. The largest variationscorrespond to the entire structure of the complex, which featuresa layered zone of interstratified harzburgite, wehrlite, andvarious pyroxenites sandwiched between zones of unlayered harzburgiteand dunite containing only minor pyroxenite bands. All zonesexhibit the same deformation microstructures, tabular equigranularto porphyroclastic textures, and strong mineral aggregate lineation.Harzburgite from the unlayered zones is characterized by olivinevalues of 100Mg/(Mg+Fe)=91–93.5 and chromite values of100Cr/(Cr+Al+Fe³⁺)=40–75. These variables exhibit a positivecorrelation, which is typical of harzburgites and lherzolitesfrom the basal units of ophiolites and from xenoliths in alkalibasalts and kimberlites. The harzburgite is therefore interpretedas a residue from partial melting in the mantle. By contrast,harzburgite in the interlayered zone features olivine valuesof 100Mg/(Mg+Fe)=88–92 and chromite values of 100Cr/(Cr+Al+Fe3+)=40–60,and in this case the variables tend to show a negative correlationin any given locality and they partly overlap data from theintercalated wehrlite and dunite. The harzburgite of the layeredzone is interpreted as residual mantle that reacted with evolvedmelts that then crystallized as wehrlite and dunite. The harzburgitein the unlayered zones is more refractory than that in the layeredzone, even after removing effects of reaction. This differencecan be explained either by enhanced partial melting and meltextraction in the unlayered zones, possibly owing to the preferentialintroduction of a waterrich fluid, or by melt segregation fromthe unlayered zones and transfer to the layered zone in responseto a piezometric pressure gradient and compaction of a solidresidual matrix. Mineralogical evidence suggests that evolvedmelts migrated through conduits formed in the layered zone byfracturing or diapirism.     

Accelerated continental-scale snowmelt and ecohydrological impacts in the four largest Siberian river basins in response to spring warming

River runoff from the four largest Siberian river basins (the Ob, Yenisei, Lena, and Kolyma) considerably contributes to freshwater flux into the Arctic Ocean from the Eurasian continent. However, the effects of variation in snow cover fraction on the ecohydrological variations in these basins are not well understood. In this study, we analysed the spatiotemporal variability of the maximum snow cover fraction (SCF_max) in the four Siberian river basins. We compared the SCF_max from 2000 to 2016 with data in terms of monthly temperature and precipitation, night-time surface temperatures, the terrestrial water storage anomaly (TWSA), the normalised difference vegetation index (NDVI), and river runoff. Our results exhibit a decreasing trend in the April SCF_max values since 2000, largely in response to warming air temperatures in April. We identified snowmelt water as the dominant control on the observed increase in the runoff contribution in May across all four Siberian river basins. In addition, we detected that the interannual river runoff was predominantly controlled by interannual variations in the TWSA. The NDVI in June was strongly controlled by the timing of the snowmelt along with the surface air temperature and TWSA in June. The rate of increase in the freshwater flux from the four Siberian rivers decreased from 2000 to 2016, exhibiting large interannual variations corresponding to interannual variations in the TWSA. However, we identified a clear increase trend in the freshwater flux of ~4 km³/year when analysing the long-term 39-year historical record (1978–2016). Our results suggest that continued global warming will accelerate the transition towards the earlier timing of snowmelt and spring freshwater flux into the Arctic Ocean. Our findings also highlight the effects of earlier snowmelt on ecohydrological changes in the Northern Hemisphere.     

Thermal History of the Horoman Peridotite Complex: A Record of Thermal Perturbation in the Lithospheric Mantle

The ascent history of the Horoman peridotite complex, Hokkaido,northern Japan, is revised on the basis of a detailed studyof large ortho- and clinopyroxene grains 1 cm in size (megacrysts)in the Upper Zone of the complex. The orthopyroxene megacrystsexhibit distinctive M-shaped Al zoning patterns, which werenot observed in porphyroclastic grains less than 5 mm in sizedescribed in previous studies. Moreover, the Al and Ca contentsof the cores of the orthopyroxene megacrysts are lower thanthose of the porphyroclasts. The Upper Zone is inferred to haveresided not only at a higher temperature than previously suggestedbut also at a higher pressure (1070°C, 2·3 GPa) thanthe Lower Zone (950°C, 1·9 GPa), in the garnet stabilityfield, before the ascent of the two zones. The Horoman complexprobably represents a 12 ± 5 km thick section of lithosphericmantle with an 10 ± 8°C/km vertical thermal gradient.The current thickness of the Horoman complex is 3 km, whichis a result of shortening of the lithospheric mantle by 0·25± 0·1 during its ascent. The Upper Zone appearsto have experienced a heating event during its ascent throughthe spinel stability field, with a peak temperature as highas 1200°C. The effect of heating decreases continuouslytowards the base of the complex, and the lowermost part of theLower Zone underwent very minor heating at a pressure higherthan 0·5 GPa. The uplift and associated deformation,as well as heating, was probably driven by the ascent of a hotasthenospheric upper-mantle diapir into the Horoman lithosphere. KEY WORDS: Horoman; P–T trajectory; thermal history; Al diffusion in pyroxene; geothermobarometry     

Olivine-spinel geospeedometry: Analysis of diffusion-controlled Mg-Fe2+ exchange

Subsolidus Mg-Fe²⁺ exchange between olivine and spinel is governed by Mg-Fe²⁺ interdiffusion. Incomplete exchange results in Mg-Fe²⁺ heterogeneity in both olivine and spinel, which provides information on the thermal histories of the host rocks. A composite sphere model has been developed to obtain quantitative cooling rates or heating duration from the Mg-Fe²⁺ heterogeneity. The model assumes that a spherical core of spinel and a surrounding semi-infinite spherical shell of olivine interact by diffusion-controlled exchange of Mg and Fe²⁺. The differential equations describing the model are solved numerically by finite difference approximations. The numerical solution reveals that cooling rates or heating duration can be estimated from the relationship between the grain size of spinel and temperature calculated from the chemical compositions of the core of a spinel grain and of olivine far away from it. The calculated temperature is employed in place of $\text{Mg}\text{(Mg + Fe}{}^{\mathrm{2+}}\text{)}$ at the center of spinel to obtain the absolute temperature of thermal events.This olivine-spinel geospeedometer has been applied to peridotites. gabbro, and picrites from some ophiolite complexes in Japan to estimate their cooling rates. The estimated cooling rates for the peridotites range from 10^?4 to 10^?1 °C/yr, and those for the picrites from 10³ to 10⁴ °C/yr. The geospeedometer has been extended to estimate the heating duration of lherzolite xenoliths in basalt from the Ichinomegata crater, northeast Japan. The estimated heating duration of the xenoliths is less than one day.     

Evaluation of snow models in terrestrial biosphere models using ground observation and satellite data: impact on terrestrial ecosystem processes

Snow is important for water management, and an important component of the terrestrial biosphere and climate system. In this study, the snow models included in the Biome‐BGC and Terrestrial Observation and Prediction System (TOPS) terrestrial biosphere models are compared against ground and satellite observations over the Columbia River Basin in the US and Canada and the impacts of differences in snow models on simulated terrestrial ecosystem processes are analysed. First, a point‐based comparison of ground observations against model and satellite estimates of snow dynamics are conducted. Next, model and satellite snow estimates for the entire Columbia River Basin are compared. Then, using two different TOPS simulations, the default TOPS model (TOPS with TOPS snow model) and the TOPS model with the Biome‐BGC snow model, the impacts of snow model selection on runoff and gross primary production (GPP) are investigated. TOPS snow model predictions were consistent with ground and satellite estimates of seasonal and interannual variations in snow cover, snow water equivalent, and snow season length; however, in the Biome‐BGC snow model, the snow pack melted too early, leading to extensive underpredictions of snow season length and snow covered area. These biases led to earlier simulated peak runoff and reductions in summer GPP, underscoring the need for accurate snow models within terrestrial ecosystem models. Copyright © 2007 John Wiley & Sons, Ltd.     

Anisotropy of Mg isotopic fractionation during evaporation and Mg self-diffusion of forsterite in vacuum

Evaporation of solid materials under low-pressure conditions could play important roles in chemical and isotopic fractionations in the early solar system. We have studied anisotropy of isotopic fractionation of ²⁶Mg and ²⁵Mg during kinetic evaporation of forsterite (Mg₂SiO₄), which is potentially a powerful tool to understand thermal histories of crystals in the early solar system. Ion-microprobe depth profiling revealed that the Mg isotopic zoning profiles of forsterite evaporated at 1500-1700 °C are notably differing along the a-, b-, and c-axes, which can be attributed to anisotropy in self-diffusion coefficient of Mg (D) and an isotopic fractionation factor for evaporation of Mg (α). The D and α were obtained from zoning profiles by applying the diffusion-controlled isotopic fractionation model of Wang et al. [1999. Evaporation of single crystal forsterite: Evaporation kinetics, magnesium isotope fractionation, and implications of mass-dependent isotopic fractionation of a diffusion-controlled reservoir. Geochim. Cosmochim. Acta 63(6), 953-966.].The D is largest and smallest along the a- and c-axes, respectively. The activation energy of 560-670 kJ/mol indicates that Mg diffusion at 1500-1700 °C occurred in the intrinsic diffusion regime.The α seems to be larger along the a- or c-axes than along the b-axis. The α along the a- or c-axes show weak temperature dependence. The α along all the crystallographic orientations is closer to unity than that expected from the kinetic theory of gases. These lines of evidence suggest that surface processes such as breaking of bonds and surface diffusion are responsible for the isotopic fractionation.     

Fractal analysis of the fracture strength of lava dome material based on the grain size distribution of block-and-ash flow deposits at Unzen Volcano, Japan

A fractal theory of rock fragmentation is applied to block-and-ash flow deposits from the Fugendake dome, Unzen Volcano, Kyushu, Japan, in order to analyze the material strength and the energy required for size reduction of the source lava dome. Two fractal dimensions h and D_s, which are mutually interchangeable, represent the relative strength and energy for particles reduced to a given size. They can be theoretically estimated from the power relations of a reference grain size to the cumulative mass and number of fragments smaller than that size. The Unzen–Fugendake block-and-ash flow deposits have been further modified by size sorting and secondary fragmentation that occurred during flowage, so that the h value decreases (or D_s value increases) with increasing distance from the source. Coarse, reversely graded deposits are, however, found to retain the original size population relatively well. The D_s values estimated from deposits of this type are compatible with those previously reported from decompression–fragmentation experiments conducted on the same dome material. The employed fractal approach could thus give insights into the potential mode of dome collapse that generates block-and-ash flows.