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.     

Distribution of sea cucumbers in relation to sediment characteristics in coral reef lagoons and adjacent waters around Ishigaki Island,southern Japan

To elucidate which environmental factors affect lagoonal‐scale sea cucumber distributions in Ishigaki Island, Okinawa, Japan, intertidal and subtidal areas of three coral reef lagoons were classified into several ground divisions by bottom characteristics, and sea cucumber densities therein were compared with the composition of sediment cover, grain size and organic content, and coverage of macroalgae, seagrass, and massive corals. Holothuria atra, Holothuria leucospilota, Stichopus chloronotus, and Synapta maculata had highest densities in the nearshore areas but were rare in reef flats, probably because of wave disturbance and low areal cover of sand sediment as potential feeding environments. No relationship was observed between sea cucumber densities and sediment organic content or grain size. Thus, even if these sea cucumbers have selectivity for habitats with a high sediment organic content, the effect of such selectivity on the distribution seems to be limited to relatively small areas. The sea cucumber distributions can be classified by bottom sediment/biota composition into bedrock (H. leucospilota), sand (H. atra), and lagoonal types (St. chloronotus and Sy. maculata). These habitat selections were possibly related to various aspects of sea cucumber ecology such as refuge from predators or turbulence, or settlement and nursery place, which have implications for importance of the complexity of lagoonal‐scale topography and sediment/biota conditions for the coexistence of various holothurian species.     

Formation age of the lunar crater Giordano Bruno

Abstract— Using the Terrain Camera onboard the Japanese lunar explorer, SELENE (Kaguya), we obtained new high‐resolution images of the 22‐kilometer‐diameter lunar crater Giordano Bruno. Based on crater size‐frequency measurements of small craters (<200 m in diameter) superposed on its continuous ejecta, the formation age of Giordano Bruno is estimated to be 1 to 10 Ma. This is constructive evidence against the crater's medieval age formation hypothesis.     

Ion migration in MgSiO3-perovskite and olivine by molecular dynamics calculations

We applied molecular dynamics (MD) simulations to finding likely paths of atomic migration of the Mg ion in forsterite (Mg₂SiO₄) and the oxygen ion in MgSiO₃-perovskite to better understand atomic diffusion in minerals. Our simulations show that there exist two routes of Mg migration in the forsterite structure, that is, paths between the M1 and M2 sites and between the M1 and M1 sites. In the MgSiO₃-perovskite structure, some oxygen ions migrate to the next sites all together through the O1 vacant site showing co-operative movements. The O ions are relatively mobile mainly along the b axis in the perovskite structure. Meta-stable sites are often present between a stable site and another stable site on atomic migration. In spite of many assumptions, our MD simulations may show likely paths of atomic migration in crystal.     

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.     

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.     

Influence of Kuroshio Flow on the Horizontal Distribution of North Pacific Intermediate Water in the Shikoku Basin

The influence of the Kuroshio flow on the horizontal distribution of North Pacific Intermediate Water (NPIW) in the Shikoku Basin is examined based upon observational data collected by the training vessel “Seisui-maru” of Mie University together with oceanographic data compiled by the Japan Oceanographic Data Center (JODC). Although it has been stated that the NPIW with salinity less than 34.2 psu had been confined to the south of the Kuroshio main axis along the PT (KJ) Line on the eastern side of the Izu Ridge, a similar tendency can be detected on the western side of the Izu Ridge. Namely, the NPIW on the southern side of the Kuroshio main axis in the Shihoku Basin does not indicate a tendency to go northward across the Kuroshio main axis without an increase in salinity of more than 34.2 psu. However, the JODC data show that less saline water (<34.2 psu) was present on the northern side of the Kuroshio main axis south of the Kii Peninsula in May 1992. Satellite observed sea surface temperature (SST) data suggested that the Kuroshio approaches the Kii Peninsula after forming a small meander off Kyushu and some intrusions of the NPIW into the northern coastal side of the Kuroshio main axis occurred in this period. It is concluded that intrusion of the NPIW with salinity less than 34.2 psu to the northern coastal side through the Kuroshio main axis occurred during the decay period of the small meander path in May 1992. Based on these observational results, the source of the salinity minimum water on the northern coastal side of the Kuroshio main axis is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modeling the kinetics of silica nanocolloid formation and precipitation in geologically relevant aqueous solutions

The kinetics of the formation and precipitation of nanocolloidal silica from geologically relevant aqueous solutions is investigated. Changes in monomeric (SiO_2(mono)), nanocolloidal (SiO_2(nano)) and precipitated silica (SiO_2(ppt)) concentrations in aqueous solutions from pH 3 to 7, ionic strengths (IS) of 0.01 and 0.24 molal, and initial SiO₂ concentrations of 20.8, 12.5 and 4.2 mmolal (reported in [Icopini, G.A., Brantley, S.L., Heaney, P.J., 2005. Kinetics of silica oligomerization and nanocolloid formation as a function of pH and ionic strength at 25 °C. Geochim. Cosmochim. Acta69(2), 293-303.]) were fit using two kinetic models. The first model, termed the concentration model, is taken from Icopini et al. (2005) and assumes that the rate of change of SiO_2(mono) as a function of time has a fourth-order dependence on the concentration of SiO_2(mono) in solution. The second model, termed the supersaturation model, incorporates the equilibrium concentration of amorphous silica and predicts that polymerization will be a function of the degree of silica supersaturation in solution with respect to amorphous silica. While both models generally predicted similar rate constants for a given set of experimental conditions, the supersaturation model described the long-term equilibrium behavior of the SiO_2(mono) fraction more accurately, resulting in significantly better fits of the monomeric data. No difference was seen between the model fits of the nanocolloidal silica fraction. At lower pH values (3-4), a metastable equilibrium was observed between SiO_2(mono) and SiO_2(nano). This equilibrium SiO_2(mono) concentration was found to be 6 mmolal, or three times the reported solubility of bulk amorphous silica under the experimental conditions studied and corresponds to the predicted solubility of amorphous silica colloids approximately 3 nm in diameter. Atomic force microscopy was used to determine the average size of the primary nanocolloidal particles to be ∼3 nm, which is in direct agreement with the solubility calculations. Larger aggregates of the primary nanocolloids were also observed to range in size from 30 to 40 nm. This work provides the first kinetic models describing the formation and evolution of nanocolloidal silica in environmentally relevant aqueous solutions. Results indicate that nanocolloidal silica is an important species at low pH and neutral pH at low ionic strengths and may play a more important role in geochemical cycles in natural aqueous systems than previously considered.