Mineralogy and petrography of C asteroid regolith: The Sutter's Mill CM meteorite

Based upon our characterization of three separate stones by electron and X‐ray beam analyses, computed X‐ray microtomography, Raman microspectrometry, and visible‐IR spectrometry, Sutter's Mill is a unique regolith breccia consisting mainly of various CM lithologies. Most samples resemble existing available CM2 chondrites, consisting of chondrules and calcium‐aluminum‐rich inclusion (CAI) set within phyllosilicate‐dominated matrix (mainly serpentine), pyrrhotite, pentlandite, tochilinite, and variable amounts of Ca‐Mg‐Fe carbonates. Some lithologies have witnessed sufficient thermal metamorphism to transform phyllosilicates into fine‐grained olivine, tochilinite into troilite, and destroy carbonates. One finely comminuted lithology contains xenolithic materials (enstatite, Fe‐Cr phosphides) suggesting impact of a reduced asteroid (E or M class) onto the main Sutter's Mill parent asteroid, which was probably a C class asteroid. One can use Sutter's Mill to help predict what will be found on the surfaces of C class asteroids such as Ceres and the target asteroids of the OSIRIS‐REx and Hayabusa 2 sample return missions (which will visit predominantly primitive asteroids). C class asteroid regolith may well contain a mixture of hydrated and thermally dehydrated indigenous materials as well as a significant admixture of exogenous material would be essential to the successful interpretation of mineralogical and bulk compositional data.     

Precise estimation of pressure–temperature paths from zoned minerals using Markov random field modeling: theory and synthetic inversion

The chemical zoning profile in metamorphic minerals is often used to deduce the pressure–temperature (P–T) history of rock. However, it remains difficult to restore detailed paths from zoned minerals because thermobarometric evaluation of metamorphic conditions involves several uncertainties, including measurement errors and geological noise. We propose a new stochastic framework for estimating precise P–T paths from a chemical zoning structure using the Markov random field (MRF) model, which is a type of Bayesian stochastic method that is often applied to image analysis. The continuity of pressure and temperature during mineral growth is incorporated by Gaussian Markov chains as prior probabilities in order to apply the MRF model to the P–T path inversion. The most probable P–T path can be obtained by maximizing the posterior probability of the sequential set of P and T given the observed compositions of zoned minerals. Synthetic P–T inversion tests were conducted in order to investigate the effectiveness and validity of the proposed model from zoned Mg–Fe–Ca garnet in the divariant KNCFMASH system. In the present study, the steepest descent method was implemented in order to maximize the posterior probability using the Markov chain Monte Carlo algorithm. The proposed method successfully reproduced the detailed shape of the synthetic P–T path by eliminating appropriately the statistical compositional noises without operator’s subjectivity and prior knowledge. It was also used to simultaneously evaluate the uncertainty of pressure, temperature, and mineral compositions for all measurement points. The MRF method may have potential to deal with several geological uncertainties, which cause cumbersome systematic errors, by its Bayesian approach and flexible formalism, so that it comprises potentially powerful tools for various inverse problems in petrology.     

Study on hydraulic characteristics of sabo dam with a flap structure for debris flow

The front part of the flow is very important and complex in the case of debris flow where there is an accumulation of large boulders. It is important to control or dampen the energy of the frontal part of a debris flow for the safety of the downstream area because the impact pressure of debris flow is much greater than that of clear fluid. The main objective of this study is to analyze the hydraulic characteristics of the proposed dam (i.e. closed-type dam with flap). The vertical pressure distribution of this type is compared with conventional dam types. In the experiments, the total pressure associated with major debris flows was recorded in real time by a system consisting of four dynamic pressure sensors installed on different types of dam. The results from experimental data clearly show that the dam with the flap has advantages of capturing the debris flow with large boulders and controls the total pressure by flow circulation due to presence of the flap structure compared to a closed-type dam without flap. Further-more, the empirical coefficients of hydrodynamic and solid collision models were proposed and com-pared with available coefficients.     

Tsunami generation of the 1993 Hokkaido Nansei-Oki earthquake

Heterogeneous fault motion of the 1993 Hokkaido Nansei-Oki earthquake is studied by using seismic, geodetic and tsunami data, and the tsunami generation from the fault model is examined. Seismological analyses indicate that the focal mechanism of the first 10 s, when about a third of the total moment was released, is different from the overall focal mechanism. A joint inversion of geodetic data on Okushiri Island and the tide gauge records in Japan and Korea indicates that the largest slip, about 6 m, occurred in a small area just south of the epicenter. This corresponds to the initial rupture on a fault plane dipping shallowly to the west. The slip on the northernmost subfault, which is dipping to the east, is about 2 m, while the slips on the southern subfaults, which are steeply dipping to the west, are more than 3 m. Tsunami heights around Okushiri Island are calculated from the heterogeneous fault model using different grid sizes. Computation on the smaller grids produces large tsunami height that are closer to the observed tsunami runup heights. Tsunami propagation in the nearly closed Japan Sea is examined as the free oscillation of the Japan Sea. The excitation of the free oscillation by this earthquake is smaller than that by the 1964 Niigata or 1983 Japan Sea earthquake.     

Theoretical investigation of the time variation of drainage density

A theoretical equation was developed to express the time variation of drainage density in a basin or geomorphic surface: D_i(t, T) is the drainage density at time T on the i-th basin or geomorphic surface, which was formed at time t; β(τ) is a factor related to the erosional force causing the development of the rivers of the basin or surface at time τ; δ_i is the maximum drainage density; and D_i is the initial drainage density on the i-th geomorphic surface or basin. The equation is based on the assumption that the drainage density increases with time until it reaches a specific upper limit δ_i(t)), the maximum drainage density, which is related to certain physical properties of the basin. The equations for various dated basins or geomorphic surfaces can be combined into one modified equation if the same relative erosional forces have acted on those basins or surfaces (β(t) = β(t) and if the basins or surfaces have the same physical properties δ_i(t) = δ_i(t), (D_i = D₀). The application of this equation to coastal terraces and glacial tills shows that the model is compatible with observed drainage densities on various dated basins or surfaces.     

Estimation of debris flow discharge coefficient considering sediment concentration

A sabo dam has a purpose to block the path of debris flow. However, when overflow occurs, a sabo dam works as a weir, a vertical obstruction, where the fluid must flow over. Many empirical formulas and discharge coefficients for weirs relating flow depth to discharge have been proposed to calculate overflow discharges. However, only a few studies about overflow discharge coefficients are available in the case of debris flow. In this paper, experiments and numerical simulations were done to estimate debris flow discharge coefficients by considering the sediment concentration. In the numerical simulation, a complete overflow equation and a free overfall equation were implemented to calculate debris overflow discharges at a sabo dam. To determine the discharge coefficients for each equation, single factor regression analysis was used. Laboratory experiments were done to calibrate and to compare with the simulation. Study results showed that the discharge coefficients increase as the sediment concentration increases. This finding suggests debris flow discharge coefficients are derived to calculate the debris overflow discharges at a sabo dam.     

Secular variations of carbon and helium isotopes at Izu-Oshima Volcano, Japan

Secular variations in ¹³C/¹²C ratios and chemical compositions of gas samples from October 1986 to July 1992 are reported from a 92–95 °C steam well located about 3 km north of Mt. Mihara, an active volcano on Izu-Oshima Island, Japan. The δ¹³C value steeply increased from −2.97‰ (relative to PDB carbonate) in December 1986 to −1.15‰ in March 1988 and then gradually decreased to −1.75‰ in July 1992. Over the same period, the CO₂ content changed similarly with time, even though the experimental error is relatively large. These variations are consistent with helium isotope changes. Initially rapid and then slow enhancements of ³He/⁴He ratio, δ¹³C value and CO₂ content are invoked by violent eruptions of Izu-Oshima volcano from 15 November to 18 December 1986. After the eruptive activity, depletion of magmatic gas emission and subsequent mixing with crustal fluids in the hydrothermal system may produce the gradual decreases of ³He/⁴He ratio, δ¹³C value and CO₂ content. Taking into account the rates of these decreases, we suggest that helium and carbon isotope ratios will return to the situation of before the magmatic eruption within 15 years.     

Middle Miocene to Quaternary primary basalt and high magnesian andesite magmas of North Hokkaido, Japan: Source mantle characteristics and degrees of partial melting

Abstract   Middle Miocene to Quaternary primitive basalts and high magnesian andesite (HMA) in North Hokkaido resulted from three periods of intense volcanism; early-stage (12–10 Ma), middle-stage (9–7 Ma) and late-stage (3–0 Ma). Based on the chemical compositions of olivines and chromian spinels and bulk chemistry of the primitive rocks, we examined depths of segregation of the calculated primary magmas and the degrees of partial melting of the source mantle. In the context of asthenospheric mantle upwelling, petrological data from the present study can be accounted for by the secular change in the depth of magma segregation from the upwelled asthenospheric mantle, which is composed of fertile peridotite. Thus, the early-stage primary magmas were generated by higher degrees of partial melting of the shallower part of hot asthenospheric mantle, whereas the middle- and late-stage primary magmas resulted from lower degrees of partial melting of a deeper part of the asthenospheric mantle. The early-stage HMA magma was generated by partial melting of the remnant subcontinental lithospheric mantle composed of refractory peridotite. This melting might have resulted from an increased geothermal gradient caused by upwelling of hot asthenosphere.     

What controls the lateral variation of large earthquake occurrence along the Japan Trench?

Abstract The lateral (along trench axis) variation in the mode of large earthquake occurrence near the northern Japan Trench is explained by the variation in surface roughness of the subducting plate. The surface roughness of the ocean bottom near the trench is well correlated with the large-earthquake occurrence. The region where the ocean bottom is smooth is correlated with'typical'large underthrust earthquakes (e.g. the 1968 Tokachioki event) in the deeper part of the seismogenic plate interface, and there are no earthquakes in the shallow part (aseismic zone). The region where the ocean bottom is rough (well-developed horst and graben structure) is correlated with large normal faulting earthquakes (e.g. the 1933 Sanriku event) in the outer-rise region, and large tsunami earthquakes (e.g. the 1896 Sanriku event) in the shallow region of the plate interface zone. In the smooth surface region, the coherent metamorphosed sediments form a homogeneous, large and strong contact zone between the plates. The rupture of this large strong contact causes great under-thrust earthquakes. In the rough surface region, large outer-rise earthquakes enhance the well-developed horst and grabens. As these structure are subducted with sediments in the graben part, the horsts create enough contact with the overriding block to cause an earthquake in the shallow part of the interface zone, and this earthquake is likely to be a tsunami earthquake. When the horst and graben structure is further subducted, many small strong contacts between the plates are formed, and they can cause only small underthrust earthquakes.