Yamato 792947, 793408 and 82038: The most primitive H chondrites,with abundant refractory inclusions

Abstract— In this paper we report petrological and chemical data of the unusual chondritic meteorites Yamato (Y)‐792947, Y‐93408 and Y‐82038. The three meteorites are very similar in texture and chemical composition, suggesting that they are pieces of a single fall. The whole‐rock oxygen isotopes and the chemical compositions are indicative of H chondrites. In addition, the mineralogy, and the abundances of chondrule types, opaque minerals and matrices suggest that these meteorites are H3 chondrites. They were hardly affected by thermal and shock metamorphism. The degree of weathering is very low. We conclude that these are the most primitive H chondrites, H3.2–3.4 (S1), known to date. On the other hand, these chondrites contain extraordinarily high amounts of refractory inclusions, intermediate between those of ordinary and carbonaceous chondrites. The distribution of the inclusions may have been highly heterogeneous in the primitive solar nebula. The mineralogy, chemistry and oxygen isotopic compositions of inclusions studied here are similar to those in CO and E chondrites.     

Construction of Wide Field Cryogenic Telescope

We constructed a wide field cryogenic telescope (WFCT) containing a whole Ritchey-Chrétien system and a focal plane array in a cryostat for near infrared observations to cover the field of view of 0.4°. The telescope has a primary mirror of 220 mm and an engineering grade 256 × 256 InSb array. The optical components such as two mirrors, filters, spiders, and radiation shield tube are cooled down to 180 K as well as the InSb array to 35 K by a mechanical refrigerator. We show the results of the background surface brightnesses and the limiting magnitudes at 3.3 and 3.67 μm measured at Sutherland, South Africa. We describe the on-going upgrade of this instrument, equipped with a 1024 × 1024 ALADDIN InSb science grade array. This revised version was published online in July 2006 with corrections to the Cover Date.     

Tectonic setting and formation of the Setouchi volcanic rocks in the western Seto Inland Sea,Japan

The Setouchi volcanic rocks include high-Mg andesites (HMAs) and garnet-bearing dacite–rhyolite, and are sporadically distributed along the Median Tectonic Line, Japan. New U–Pb zircon ages and geological and geochemical data are presented for those rocks in the Western Setouchi region (W-Setouchi). Previous studies referred to the altered andesite in the W-Setouchi as “pre-Setouchi volcanic rocks.” However, on the basis of the new U–Pb age (14.4 Ma ± 0.3 Ma) and geochemical characteristics, we redefine it as the Jikamuro Formation, part of the Setouchi volcanic rocks. Incompatible elements are more enriched in the Jikamuro Formation rocks than in the Setouchi HMAs. The characteristic element compositions may be explained by mixing of compositionally different magmas, including subducted sediment melts, plus a contribution from crustal contamination. A stress-inversion technique with Bingham distribution method was applied to the orientations of felsic and mafic dikes within the Setouchi volcanic rocks, and indicates paleo-stress conditions during the period of Setouchi volcanism in the W-Setouchi. The analysis reveals NNW-extensional stresses and a strike-slip stress. We infer that the former represents extensional conditions during the main period of volcanism and the latter represents a stress transition during the most recent period of volcanism (after 12 Ma).     

High-resolution upper Pliocene to Pleistocene calcareous nannofossil biostratigraphy in Ocean Drilling Program Hole 1146A in the South China Sea

We established a high-resolution calcareous nannofossil biostratigraphy for the late Pliocene–Pleistocene by analyzing a 242 m-thick, continuous sedimentary succession from Ocean Drilling Program Site 1146, Hole A, in the South China Sea (SCS). A total of 14 calcareous nannofossil datums were detected in the SCS succession. They are, in descending order: first occurrence (FO) of Emiliania huxleyi, last occurrence (LO) of Pseudoemiliania lacunosa, LO of Reticulofenestra asanoi, FO of Gephyrocapsa parallela, FO of R. asanoi, LO of large Gephyrocapsa spp., FO of large G. spp., FO of Gephyrocapsa oceanica, FO of Gephyrocapsa caribbeanica, LO of Calcidiscus macintyrei, LO of Discoaster brouweri, LO of Discoaster pentaradiatus, LO of Discoaster surculus, and LO of Discoaster tamalis. The FO of E. huxleyi was not precisely detected due to poor preservation and dissolution of nannofossils in the underlying strata. We refined the previous calcareous nannofossil biostratigraphy in the SCS by identifying Gephyrocapsa species and four evolutionary extinction events of the genus Discoaster. The proposed calcareous nannofossil biostratigraphy correlates with those reported in other terrestrial and marine areas/sites and global benthic foraminiferal δ¹⁸O records. The age–depth curves based on nannofossil biostratigraphy indicate a significant increase in the sedimentation rates at the LO of R. asanoi (0.91–0.85 Ma). The timing of this increase corresponds to reef expansion in the Ryukyu Islands linked to a stepwise increase in Kuroshio Current intensity. This timing is broadly coeval with a sea surface temperature increase of ∼2 °C in the northwestern Pacific due to expansion of the Western Pacific Warm Pool towards the north and south subtropical regions. This can be explained by increased weathering and erosion of terrestrial areas in glacial periods and increased rainfall causing higher sediment transport in interglacial periods, which were both linked to Middle Pleistocene Transition-related climatic changes.     

Zircon U–Pb ages of Miocene granitic rocks in the Koshikijima Islands: Implications for Neogene tectonics in the Kyushu region,southwest Japan

The opening of the Japan Sea separated southwest Japan from the Eurasian continent during the Early to Middle Miocene. Since then, diverse igneous activities have occurred in relation to the subduction of the Philippine Sea Plate beneath southwest Japan. The Okinawa Trough formed in the back-arc region of the Ryukyu Arc since the Late Miocene. In the Koshikijima Islands, off the west coast of Kyushu and near the northern end of the Okinawa Trough, felsic to intermediate igneous rocks with Middle to Late Miocene radiometric ages occur as granitic intrusions and dikes. We obtained zircon U–Pb ages and whole-rock major- and trace-element compositions of Koshikijima granitic rocks to elucidate their magmagenesis. The U–Pb ages of granitic rocks in Kamikoshikijima and Shimokoshikijima and a dacite dike are about 10 Ma, suggesting that most magmatism on the Koshikijima Islands was coeval with early rifting in the Okinawa Trough. We infer that magmagenesis occurred via melting of lower crustal mafic rocks related to rifting in the Okinawa Trough based on the arc-like trace-element compositions of these I-type granites. Andesitic dikes preceded felsic igneous activity on the Koshikijima Islands, and their ages and petrochemistry will help elucidate the magmatism and tectonics in this area throughout the Miocene.     

Constituents of Vertical-component Coda Waves at Long Periods

Aki (1969) first modeled coda waves of a local earthquake as a superposition of scattered surface waves. This paper attempts to clarify the constituents of surface-wave coda at long periods at very long lapse times. For a large earthquake of magnitude 7 or larger, vertical component oscillation in periods from 90 to 180 s persists for more than 20 hours from the earthquake origin time. Although the early portion of the coda envelope is successfully modeled by assuming incoherent scattered Rayleigh waves by heterogeneities distributed all over the Earth, the later potion of the observed coda envelope (roughly later than 35,000 s) has systematically larger amplitude than theoretical prediction. To clarify the cause of this discrepancy, we studied the constituents of vertical-component seismograms of three large earthquakes recorded by the F-net in Japan using the f-k power spectral analysis. We found that the direct and scattered fundamental-mode Rayleigh waves of velocity about 3.7 km/s are dominant in the earlier part of each envelope. It justifies the use of a scattering model of the fundamental Rayleigh waves for synthesizing the envelope. At lapse times later than 20,000 s–35,000 s, higher modes with phase velocities around 20 km/s become dominant. The transition time to the dominance of higher modes is found to become earlier for a deeper focus earthquake. The small coda attenuation factor from (1.90±0.23) × 10⁻³ to (2.38±0.32) × 10⁻³ estimated from later coda envelopes recorded at IRIS stations distributed worldwide also agrees with the attenuation factor of spheroidal modes according to PREM. We may interpret that higher mode waves are uniformly distributed at large lapse time due to large velocity dispersion and/or scattering and they dominate over the fundamental mode waves because of smaller attenuation in the lower mantle. The coda attenuation measurement proposed by Aki is found to be useful even for long periods and at very large lapse times.     

Particle Image Velocimetry Measurements of Turbulent Flow Within Outdoor and Indoor Urban Scale Models and Flushing Motions in Urban Canopy Layers

We investigate the spatial characteristics of urban-like canopy flow by applying particle image velocimetry (PIV) to atmospheric turbulence. The study site was a Comprehensive Outdoor Scale MOdel (COSMO) experiment for urban climate in Japan. The PIV system captured the two-dimensional flow field within the canopy layer continuously for an hour with a sampling frequency of 30 Hz, thereby providing reliable outdoor turbulence statistics. PIV measurements in a wind-tunnel facility using similar roughness geometry, but with a lower sampling frequency of 4 Hz, were also done for comparison. The turbulent momentum flux from COSMO, and the wind tunnel showed similar values and distributions when scaled using friction velocity. Some different characteristics between outdoor and indoor flow fields were mainly caused by the larger fluctuations in wind direction for the atmospheric turbulence. The focus of the analysis is on a variety of instantaneous turbulent flow structures. One remarkable flow structure is termed ‘flushing’, that is, a large-scale upward motion prevailing across the whole vertical cross-section of a building gap. This is observed intermittently, whereby tracer particles are flushed vertically out from the canopy layer. Flushing phenomena are also observed in the wind tunnel where there is neither thermal stratification nor outer-layer turbulence. It is suggested that flushing phenomena are correlated with the passing of large-scale low-momentum regions above the canopy.     

Wind-Tunnel Experiments for Gas Dispersion in an Atmospheric Boundary Layer with Large-Scale Turbulent Motion

Large-scale turbulent motions enhancing horizontal gas spread in an atmospheric boundary layer are simulated in a wind-tunnel experiment. The large-scale turbulent motions can be generated using an active grid installed at the front of the test section in the wind tunnel, when appropriate parameters for the angular deflection and the rotation speed are chosen. The power spectra of vertical velocity fluctuations are unchanged with and without the active grid because they are strongly affected by the surface. The power spectra of both streamwise and lateral velocity fluctuations with the active grid increase in the low frequency region, and are closer to the empirical relations inferred from field observations. The large-scale turbulent motions do not affect the Reynolds shear stress, but change the balance of the processes involved. The relative contributions of ejections to sweeps are suppressed by large-scale turbulent motions, indicating that the motions behave as sweep events. The lateral gas spread is enhanced by the lateral large-scale turbulent motions generated by the active grid. The large-scale motions, however, do not affect the vertical velocity fluctuations near the surface, resulting in their having a minimal effect on the vertical gas spread. The peak concentration normalized using the root-mean-squared value of concentration fluctuation is remarkably constant over most regions of the plume irrespective of the operation of the active grid.     

Effects of Snowfall on Drifting Snow and Wind Structure Near a Surface

Wind-tunnel and numerical experiments were performed to investigate the effects of snowfall on the wind profile and the development of drifting snow. Wind profiles and mass-flux profiles of drifting snow were measured with and without artificial snowfall over a snow surface within the tunnel. Wind and shear-stress profiles and the impact speeds of the snowflakes during snowfall were also investigated numerically. During snowfall, snowflakes transfer part of their horizontal momentum to the air, which increases the stress close to the snow surface; however, the resultant modifications of the wind profiles are small. Because snowflakes have large momentum, the decomposed snow crystals that result from their collision with a surface can form a saltation layer, even over a hard snow surface where entrainment of the grains from the surface does not occur. Additionally, during snowfall, the threshold friction velocity can be lower than the impact threshold because snowflake fragmentation facilitates snow drifting. The broken crystals contribute to the increase in the number of drifting snow grains, even below the impact threshold.     

Estimation of sub-annual inter-catchment groundwater flow using short-term water balance method

Predicting inter-catchment groundwater flow (IGF) is essential because IGF greatly affects stream water discharge and water chemistry. However, methods for estimating sub-annual IGF and clarifying its mechanisms using minimal data are limited. Thus, we quantified the sub-annual IGF and elucidated its driving factors using the short-term water balance method (STWB) for three forest headwater catchments in Japan (named here catchment A, B and As). Our previous study using the chloride mass balance indicated that annual IGF of catchment A (49.0 ha) can be negligible. Therefore, we calculated the daily evapotranspiration (ET) rate using the Priestley–Taylor expression and the 5-year water balance in catchment A (2010–2014). The sub-annual IGF of the three catchments was then calculated by subtracting the ET rate from the difference between rainfall and stream discharge during the sub-annual water balance periods selected using the STWB. The IGF rates of catchment B (7.0 ha), which is adjacent to catchment A, were positive in most cases, indicating that more groundwater flowed out of the catchment than into it, and exhibited positive linear relationships with rainfall and stream discharge. This suggested that as the catchments became wetter, more groundwater flowed out of catchment B. Conversely, the IGF rates of catchment As (5.3 ha), included in catchment A, were negative in most cases, indicating that more groundwater flowed into the catchment than out from it, and exhibited negative linear relationships with rainfall and stream discharge. Given the topography of the catchments studied, infiltration into the bedrock was the probable reason for the IGF outflow from catchment B. We hypothesized that in catchment As, the discrepancy between the actual hydrological boundary and the surface topographic boundary could have caused an IGF inflow. This study provides a useful tool for determining an IGF model structure to be incorporated into rainfall-runoff models.