An empirical model for transient crater growth in granular targets based on direct observations

The present paper describes observations of crater growth up to the time of transient crater formation and presents a new empirical model for transient crater growth as a function of time. Polycarbonate projectiles were impacted vertically into soda-lime glass sphere targets using a single-stage light-gas gun. Using a new technique with a laser sheet illuminating the target [Barnouin-Jha, O.S., Yamamoto, S., Toriumi, T., Sugita, S., Matsui, T., 2007. Non-intrusive measurements of the crater growth. Icarus, 188, 506-521], we measured the temporal change in diameter of crater cavities (diameter growth). The rate of increase in diameter at early times follows a power law relation, but the data at later times (before the end of transient crater formation) deviates from the power law relation. In addition, the power law exponent at early times and the degree of deviation from a power law at later times depend on the target. In order to interpret these features, we proposed to modify Maxwell’s Z-model under the assumption that the strength of the excavation flow field decreases exponentially with time. We also derived a diameter growth model as: d(t)∝[1-exp(-βt)]^γ, where d(t) is the apparent diameter of the crater cavity at time t after impact, and β and γ are constants. We demonstrated that the diameter growth model could represent well the experimental data for various targets with different target material properties, such as porosity or angle of repose. We also investigated the diameter growth for a dry sand target, which has been used to formulate previous scaling relations. The obtained results showed that the dry sand target has larger degree of deviation from a power law, indicating that the target material properties of the dry sand target have a significant effect on diameter growth, especially at later times. This may suggest that the previously reported scaling relations should be reexamined in order to account for the late-stage behavior with the effect of target material properties.     

A 35 kyr record of organic matter composition and δ13C of n-alkanes in bog sediments close to Lake Baikal: Implications for paleoenvironmental studies

We characterized the compositions of organic compounds in a Cheremushka bog sediment core (deposited over the last 35 kyr), located at the eastern coast of Lake Baikal, to obtain basic information about the terrestrial organic matter (OM) which contributed to Lake Baikal sediments. The bog sediment was analyzed for the molecular composition of n-alkanes, lignin phenols and n-C₂₄ to C₃₀ alkanoic acids, as well as the carbon isotopic composition of plant wax derived n-C₂₇ to C₃₃ alkanes.Concentrations of lignin phenols [vanillyl (V) plus syringyl (S) phenols] normalized to total organic carbon (TOC) in the Holocene are twice those for the last glacial maximum (LGM), while concentrations of TOC-normalized n-C₂₄ to C₃₀ alkanoic acids do not change markedly in this period. Thus, the ratio of lignin phenols to n-C₂₄ to C₃₀ alkanoic acids increases from the LGM to the Holocene. This result is essentially consistent with pollen analysis indicating an expansion of woody plants in the Holocene and a prevailing herb-abundant environment for the LGM. The δ¹³C values of n-C₂₇ to C₃₃ alkanes (e.g. ?29‰ to ?33‰ for C₃₁) indicate the presence of C₃-dominant plants throughout the core.The contribution of terrestrial OM to Lake Baikal sediments was estimated using the biomarkers, on the assumption that the OM in the bog sediments is a representative of the terrestrial OM around the lake. Hence, the estimation using lignin phenol or n-C₂₄ to C₃₀ alkanoic acid parameters indicates that 11–24% of the TOC in the Academician Ridge sediments is land-derived for both the Holocene and the LGM, which is similar to the estimates from C/N values of bulk OM. However, the estimates for terrestrial OM using the n-C₂₇ to C₃₃ alkane parameter are generally higher than those using lignin phenol or n-C₂₄ to C₃₀ alkanoic acid parameters. The difference is thought to be associated with the difference in source and behavior of these biomarkers.     

Automated cloud tracking system for the Akatsuki Venus Climate Orbiter data

Japanese Venus Climate Orbiter, Akatsuki, is cruising to approach to Venus again although its first Venus orbital insertion (VOI) has been failed. At present, we focus on the next opportunity of VOI and the following scientific observations.We have constructed an automated cloud tracking system for processing data obtained by Akatsuki in the present study. In this system, correction of the pointing of the satellite is essentially important for improving accuracy of the cloud motion vectors derived using the cloud tracking. Attitude errors of the satellite are reduced by fitting an ellipse to limb of an imaged Venus disk. Next, longitude–latitude distributions of brightness (cloud patterns) are calculated to make it easy to derive the cloud motion vectors. The grid points are distributed at regular intervals in the longitude–latitude coordinate. After applying the solar zenith correction and a highpass filter to the derived longitude–latitude distributions of brightness, the cloud features are tracked using pairs of images. As a result, we obtain cloud motion vectors on longitude–latitude grid points equally spaced. These entire processes are pipelined and automated, and are applied to all data obtained by combinations of cameras and filters onboard Akatsuki. It is shown by several tests that the cloud motion vectors are determined with a sufficient accuracy. We expect that longitude–latitude data sets created by the automated cloud tracking system will contribute to the Venus meteorology.     

Spinifex-like textured metaperidotites from the Higo Metamorphic Rocks,Japan, a possible high-pressure dehydration product of antigorite serpentinite

Spinifex-like textured metaperidotites from the Higo Metamorphic Rocks (HMR), west-central Kyushu, Japan, may be formed by high-pressure dehydration of antigorite, and may indicate deep subduction of serpentinite reaching a pressure–temperature condition of 1.6 GPa and 740–750 °C. Three rock types have been identified based on mineral assemblage and rock texture: Type I (L) consisting of medium-grained (1–5 cm long) olivine + enstatite + chromite ±tremolite with secondary talc and anthophyllite that occurs in low-grade metamorphic rocks of the biotite zone, Type I (H) of coarse-grained (up to 10 cm long) olivine + enstatite (with clinoenstatite lamella) + chromite ±tremolite with secondary talc that occurs in high-grade metamorphic rocks of the garnet-cordierite zone, and Type II composed of Al-spinel + chlorite + olivine + apatite + ilmenite with minor sodic gedrite in the garnet-cordierite zone together with Type I (H). Olivines in all rock types are mostly serpentinized during exhumation. The chromite-olivine thermometer gives 560–690 °C for Type I (L) rocks, and the spinel-olivine thermometer gives 610–740 °C for Type II rocks. The peak metamorphic pressure will be higher than 1.6 GPa based on the location of the experimentally determined invariant point (P = 1.6 GPa and T = 670 °C) of antigorite + forsterite + enstatite + talc + H₂O. This estimate is consistent with the occurrence of chlorite in Type II rocks, which is stable up to 890 °C at 2.0 GPa. The spinifex-like textured metaperidotites occur as small bodies in the low P/T type gneisses, implying tectonic juxtaposition of them probably during exhumation of the HMR. Recent findings of medium pressure (0.9–1.2 GPa) granulites and gneisses from the HMR may indicate that the HMR has a deep root into the wedge mantle from which the spinifex-like textured metaperidotites have derived.     

Bayesian and Neural Network Approaches to Estimate Deep Temperature Distribution for Assessing a Supercritical Geothermal System: Evaluation Using a Numerical Model

The temperature distribution at depth is a key variable when assessing the potential of a supercritical geothermal resource as well as a conventional geothermal resource. Data-driven estimation by a machine-learning approach is a promising way to estimate temperature distributions at depth in geothermal fields. In this study, we developed two methodologies—one based on Bayesian estimation and the other on neural networks—to estimate temperature distributions in geothermal fields. These methodologies can be used to supplement existing temperature logs, by estimating temperature distributions in unexplored regions of the subsurface, based on electrical resistivity data, observed geological/mineralogical boundaries, and microseismic observations. We evaluated the accuracy and characteristics of these methodologies using a numerical model of the Kakkonda geothermal field, Japan, where a temperature above 500 °C was observed below a depth of about 3.7 km. When using geological and geophysical knowledge as prior information for the machine learning methods, the results demonstrate that the approaches can provide subsurface temperature estimates that are consistent with the temperature distribution given by the numerical model. Using a numerical model as a benchmark helps to understand the characteristics of the machine learning approaches and may help to identify ways of improving these methods.

     

The anisotropy of cosmic ray arrival directions around 10 eV

Anisotropy in the arrival directions of cosmic rays with energies above 10¹⁷ eV is studied using data from the Akeno 20 km² array and the Akeno Giant Air Shower Array (AGASA), using a total of about 114 000 showers observed over 11 years. In the first harmonic analysis, we have found a strong anisotropy of 4% around 10¹⁸ eV, corresponding to a chance probability of 0.2% after taking the number of independent trials into account. with two-dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.     

Calibration and validation of the ocean color version-3 product from ADEOS OCTS

We present calibration and validation results of the OCTS’s ocean color version-3 product, which mainly consists of the chlorophyll-a concentration (Chl-a) and the normalized water-leaving radiance (_nL_w). First, OCTS was calibrated for the inter-detector sensitivity difference, offset, and absolute sensitivity using external calibration source. It was also vicariously calibrated using in-situ measurements for water (Chl-a and_nL_w) and atmosphere (optical thickness), which were acquired synchronously with OCTS under cloud-free conditions. Second, the product was validated using selected 17 in-situ Chl-a and 11 in-situ_nL_w measurements. We confirmed that Chl-a was estimated with an accuracy of 68% for Chl-a less than 2 mg/m³, and_nL_w from 94% (band 2) to 128% (band 4). Geometric accuracy was improved to 1.3 km. Stripes were significantly reduced by modifying the detector normalization factor as a function of input radiance.     

Magma ascent beneath Unzen Volcano, SW Japan, deduced from the electrical resistivity structure

The resistivity structure of Unzen Volcano has been revealed by extensive magnetotelluric surveys since the first eruption on November 17, 1990. This structure comprises a highly resistive surface layer, a low-resistive second layer at several hundred meters depth, interpreted as a water-saturated layer, a resistive third layer, and a low-resistive fourth layer at 10 km depth, possibly related to the deep magmatic activity. The structure has influenced the volcanic activity of Unzen. This activity was characterized by a series of dramatic changes in eruption type: a minor phreatic eruption on November 17, 1990; phreatic eruptions after February 12, 1991, preceded by several weeks of volcanic tremor; phreatomagmatic eruptions after April 9, and dome effusion beginning May 19, 1991. This paper presents a hypothesis in which the top of the magma column rose about 20 m/day, reached the base of the water-saturated layer at the end of January, 1991, and approached the upper boundary of this layer on April 9. Thus, the temporal change of eruption type and associated phenomena are systematically explained by an interaction between magma and groundwater contained in the saturated layer.     

Removal of water‐surface reflection effects with a temporal minimum filter for UAV‐based shallow‐water photogrammetry

The recent development of structure‐from‐motion (SfM) and multi‐view stereo (MVS) photogrammetry techniques has enabled semi‐automatic high‐resolution bathymetry using aerial images taken by consumer‐grade digital cameras mounted on unmanned aerial vehicles (UAVs). However, the applicability of these techniques is sometimes limited by sun and sky reflections at the water surface, which render the point‐cloud density and accuracy insufficient. In this research, we present a new imaging technique to suppress the effect of these water‐surface reflections. In this technique, we order a drone to take a short video instead of a still picture at each waypoint. We then apply a temporal minimum filter to the video. This filter extracts the smallest RGB values in all the video frames for each pixel, and composes an image with greatly reduced reflection effects. To assess the performance of this technique, we applied it at three small shallow‐water sites. Specifically, we evaluated the effect of the technique on the point cloud density and the accuracy and precision of the photogrammetry. The results showed that the proposed technique achieved a far denser point cloud than the case in which a randomly chosen frame was used for each waypoint, and also showed better overall accuracy and precision in estimating water‐bottom elevation. The effectiveness of this new technique should depend on the surface wave state and sky radiance distribution, and this dependence, as well as the applicability to large areas, should be investigated in future research. Copyright © 2018 John Wiley & Sons, Ltd.     

Simulation of the effect of wind speedup in the formation of transverse dune fields

A computer simulation model for transverse‐dune‐field dynamics, corresponding to a uni‐directional wind regime, is developed. In a previous formulation, two distinct problems were found regarding the cross‐sectional dune shape, namely the erosion in the lee of dunes and the steepness of the windward slopes. The first problem is solved by introducing no erosion in shadow zones. The second issue is overcome by introducing a wind speedup (shear velocity increase) factor, which can be accounted for by adding a term to the original transport length, which is proportional to the surface height. By incorporating these features we are able to model dunes whose individual shape and collective patterns are similar to those observed in nature. Moreover we show how the introduction of a non‐linear shear‐velocity‐increase term leads to the reduction of dune height, and this may result in an equilibrium dune field configuration. This is thought to be because the non‐linear increase of the transport length makes the sand trapping efficiency lower than unity, even for higher dunes, so that the incoming and the outgoing sand flux are in balance. To fully describe the inter‐dune morphology more precise dynamics in the lee of the dune must be incorporated. Copyright © 2000 John Wiley & Sons, Ltd.