A numerical simulation of the distribution of water temperature and salinity in the Seto Inland Sea

Constant flows, as well as oscillatory tidal flow, play an important role in the long-term dispersion of water in the Seto Inland Sea. Two kinds of numerical model (1-line and 2-line models) of the Seto Inland Sea have been developed to determine the role of density-induced currents, one type of the constant flow, in water dispersion in the Inland Sea. The seasonal variations of temperature, salinity and density fields are simulated and the density-induced current field is predicted at the same time. It is found that the most appropriate value of the longitudinal eddy diffusion coefficient,K _x, is 5×10⁶–7×10⁶ cm²sec^–1. The value of the overall mean dispersion coefficient is of the order of 10⁷cm²sec^–1 (Hayami and Unoki, 1970). Consequently, it is suggested that 50–70% of the total dispersion in the Seto Inland Sea can be attributed to currents other than density-induced currents,i.e., tidal currents, tide-induced currents and wind-driven currents.In winter, both density and velocity fields, calculated using the 1-line model, satisfy the conditions for the existence of a coastal front in Kii Channel and in the eastern Iyo-nada.     

The location and size of a solar hard X-ray burst on September 27, 1969

The location and size of a solar impulsive hard X-ray burst have been determined in one dimension to a considerable precision with a balloon-borne X-ray modulation collimator. The center of the X-ray source is on the line passing through the center of a big H flare region of 3 arc min. The size of the X-ray source is remarkably smaller and may be one arc min or less.     

Experimental calibration of sapphirine–spinel Fe–Mg exchange thermometer: Implication for constraints on P–T condition of Howard Hills, Napier Complex, East Antarctica

The Fe²⁺–Mg distribution coefficients between sapphirine and spinel:

were experimentally determined at pressures of 9–13 kbar and temperatures of 950–1150 °C using a natural ultrahigh-temperature (UHT) granulite with paragenesis of these minerals from the Napier Complex in East Antarctica [X_Mg = Mg / (Fe + Mg); X_Fe = Fe / (Fe + Mg)]. A new sapphirine–spinel geothermometer has been obtained as:

We applied the exchange thermometer to UHT or high-grade metamorphic rocks that were reported from various complexes in the world. If the K_D values of 2.63–4.34 obtained from low-Cr mineral pairs such as X_Cr^Spr < 0.016 and X_Cr^Spl < 0.047 were substituted into the equation, their temperature conditions would be estimated as 806–1050 °C at 11 kbar. The X_Cr means Cr / (Al + Cr(+ Fe³⁺)). These temperatures are reasonable retrograde or near peak metamorphic condition.     

Thermal metamorphism of CI and CM carbonaceous chondrites: An internal heating model

Abstract— Infrared diffuse reflectance spectra were measured for several thermally metamorphosed carbonaceous chondrites with CI-CM affinities which were recently found from Antarctica. Compared with other CI or CM carbonaceous chondrites, these Antarctic carbonaceous chondrites show weaker absorption bands near 3 μm due to hydrous minerals, and weaker absorption bands near 6.9 μm due to carbonates, interpreted as thermal metamorphic features. These absorption bands also disappear in the spectra of samples of the Murchison (CM) carbonaceous chondrite heated above 500 °C, implying that the metamorphic temperatures of the Antarctic carbonaceous chondrites considered here were higher than about 500 °C. Model calculations were performed to study thermal metamorphism of carbonaceous chondrites in a parent body internally heated by the decay of the extinct nuclide ²⁶Al. The maximum temperature of the interior of a body more than 20 km in radius is 500–700 °C for the bulk Al contents of CI and CM carbonaceous chondrites, assuming a ratio of ²⁶Al/²⁷Al = 5 × 10^?6 which has been previously proposed for an ordinary-chondrite parent body. The metamorphic temperatures experienced by the Antarctic carbonaceous chondrites considered here may be attainable by an internally heated body with an ²⁶Al/²⁷Al ratio similar to that inferred for an ordinary-chondrite parent body.     

Surface morphological features of boulders on Asteroid 25143 Itokawa

On the sub-kilometer S-type asteroid, 25143 Itokawa, some boulders on rough terrains seem to be exposed without any powdery material covering. Based on surface morphological features, there are two major types of boulders: one has rounded edges and corners (rounded boulders), while the other has angular edges and corners (angular boulders). The surface features of the rounded boulders suggest that they have hardness heterogeneity and that some may be breccias. The angular boulders appear to be more resistant to impact disruption than the rounded ones, which may be due to a difference in lithology. The major constituents of Itokawa may be LL chondrite-like brecciated lithology (rounded boulders) along with a remarkable number of boulders suggesting that lithology is atypical among LL chondrites (angular boulders). Some of both types of boulders contain intersecting and stepped planar foliations. Comparison with meteorite ALH76009 suggests that the planar foliations may be marks where rocks were torn apart. As lithified breccias cannot be formed on present-day sub-kilometer-sized Itokawa, it is reasonable that boulders with various lithologies on Itokawa were formed on its large ancestor(s). The rubble-pile structure of Itokawa suggested by its low density (∼1.9 g/cm³) indicates that boulders on Itokawa are reassembled fragments formed by catastrophic disruption of large ancestor(s).     

The spatial and seasonal distributions of air-transport origins to the Antarctic based on 5-day backward trajectory analysis

Transport of moisture-bearing air to the Antarctic is one of the important factors that control the mass balance of the ice sheet. Here, we investigate the distribution of air-parcel transport using a backward trajectory analysis over the entire Antarctic ice sheet, based on whether the air parcel was located inside or outside of Antarctica at 5 days before arrival. At this time, we considered the air from outside Antarctica to be moisture rich. Oceanic air was found to dominate in West Antarctica throughout the year, whereas air from inland was more prominently distributed around East Antarctica, especially in summer. In East Antarctica, there was a significant seasonal variation: air from inland dominated in summer, while air of oceanic origin dominated in winter. The distribution of air parcels that came from oceanic/inland sources was similar to the accumulation map (based on satellite data), which indicates that oceanic air parcels could be a substitute for moisture transport to the Antarctic. To determine the future impacts of climate change (e.g., sea level rise), more precise predictions of the variations in the surface mass balance will be required. Our results contribute towards the improved understanding of the spatial distributions of accumulation and aerosols found in Antarctic snow and ice cores.     

The reduction effects of mangrove forest on a tsunami based on field surveys at Pakarang Cape,Thailand and numerical analysis

Using an integrated approach including satellite imagery analysis, field measurements, and numerical modeling, we investigated the damage to mangroves caused by the 2004 Indian Ocean tsunami at Pakarang Cape in Pang Nga Province, Thailand. Comparing pre- and post-tsunami satellite imagery of the study area, we found that approximately 70% of the mangrove forest was destroyed by the tsunami. Based on field observations, we found that the survival rate of mangroves increased with increasing stem diameter. Specifically, we found that 72% of Rhizophora trees with a 25–30 cm stem diameter survived the tsunami impact, whereas only 19% with a 15–20 cm stem diameter survived. We simulated the 2004 Indian Ocean tsunami using the nonlinear shallow-water wave theory to reproduce the tsunami inundation flow and investigated the bending moment acting on the mangrove trees. Results of the numerical model showed that the tsunami inundated areas along the mangrove creeks, and its current velocity reached 5.0 m s⁻¹. Based on the field measurements and numerical results, we proposed a fragility function for mangroves, which is the relationship between the probability of damage and the bending stress caused by the maximum bending moment. We refined the numerical model to include the damage probability of mangrove forests using the obtained fragility function to investigate the tsunami reduction effect of mangrove forest. Under simple numerical conditions related to the mangrove forest, ground level, and incident wave, the model showed that a mangrove forest of Rhizophora sp. with a density of 0.2 trees m⁻² and a stem diameter of 15 cm in a 400 m wide area can reduce the tsunami inundation depth by 30% when the incident wave is assumed to have a 3.0 m inundation depth and a wave period of 30 min at the shoreline. However, 50% of the mangrove forest is destroyed by a 4.5 m tsunami inundation depth, and most of the mangrove forest is destroyed by a tsunami inundation depth greater than 6 m. The reduction effect of tsunami inundation depth decreased when the tsunami inundation depth exceeded 3 m, and was mostly lost when the tsunami inundation depth exceeded 6 m.     

A multi‐geochronological study of the Hakusan volcano,central Japan

The Hakusan volcano, central Japan, is located in a region where two subducting plates (the Pacific Plate and the Philippine Sea Plate) overlap near the junction of four plates adjacent to the Japanese Islands (the Pacific Plate, the Philippine Sea Plate, the Eurasia Plate, and the North American Plate). The Hakusan volcano consists of products from four major volcanic episodes: Kagamuro, Ko‐hakusan, and Shin‐Hakusan I and II. To date the eruption events of the Hakusan volcano we applied thermoluminescence and fission track methods. ²³⁸U(²³⁴U)–²³⁰Th disequilibrium and ²⁰⁶Pb/²³⁸U methods were applied to date the zircon crystallization ages for estimating the magma residence time before the eruptions. The eruption ages we obtained are ca 250 ka for Kagamuro, ca 100 ka and ca 60 ka for Ko‐Hakusan, ca 50 ka for Shin‐Hakusan I, and <10 ka for Shin‐Hakusan II. They are concordant with previous reports based on K–Ar dating. Some of the pyroclastic rocks, possibly originating from Shin‐Hakusan II activities, are dated to be ca 36 ka or 50 ka, and belong to the Shin‐Hakusan I activity. The zircon crystallization ages show several clusters prior to eruption. The magma residence time was estimated for each volcanic activity by comparing the major crystallization events and eruption ages, and we found a gradual decrease from ca. 500 ky for the Kagamuro activity to ca. 5 ky for the Shin‐Hakusan II activity. This decrease in residence time may be responsible for the decrease in volume of erupted material estimated from the current topography of the region. The scale of volcanic activity, which was deduced from the number of crystallized zircons, is more or less constant throughout the Hakusan volcanic activity. Therefore, the decrease in magma residence time is most likely the result of stress field change.     

Measurements of the Turbulent Energy Dissipation Rate around the Shelf Break in the East China Sea

Using the micro-structure profiler, TurboMAP, large values for the turbulent energy dissipation rate ε were found just above the bottom of the shelf and around the thermocline near the continental shelf break in the East China Sea. The values found above the bottom are produced by the bottom stress due to tidal currents, resulting in a distinct bottom mixed layer where the vertical eddy diffusivity Kz is also large. Distinct maxima in the values of ε detected around the thermocline are located at the depth of the fine-scale shear maxima detected with the moored ADCP. The vertical profiles of ε were compared with those of the current velocity, and it was found that the maxima in ε appear to correspond to those of the shear with fine scale. The magnitude of the observed ε coincided approximately with the ε calculated from the fine-scale shear and the buoyancy frequency according to the parameterization proposed by Gregg (1989), if the large-scale mean shear caused by the Kuroshio is subtracted. However, it is not clear whether the parameterization for the internal wave fields in the open ocean is applicable to the estimation of ε in the shelf break. Whereas the most predominant value of ε was found just above the bottom and around the thermocline, the maxima of ε could be found in the internal area. They could have been caused by the propagation of the vertically high wave number internal tides along the characteristic ray.