Non‐equilibrium concepts lead to a unified explanation of the formation of chondrules and chondrites

Abstract— Calculations of the formation of seven types of chondrules in Semarkona from a gas of solar composition were performed with the FACT computer program to predict the chemistries of oxides (including silicates), developed by the authors and their colleagues. The constrained equilibrium theory was used in the calculations with two nucleation constraints suggested by nucleation theory. The first constraint was the blocking of Fe and other metal gaseous atoms from condensing to form solids or liquids because of very high surface free energies and high surface tensions of the solid and liquid metals, respectively. The second constraint was the blocking of the condensation of solids and the formation of metastable liquid oxides (including silicates) well below their liquidus temperatures. Our laboratory experiments suggested subcooling of type IIA chondrule compositions of 400 degrees or more below the liquidus temperature. The blocking of iron leads to a supersaturation of Fe atoms, so that the partial pressure of Fe (p_Fe) is larger than the partial pressure at equilibrium (p_Fe(eq)). The supersaturation ratio S = p_Fe/p_Fe(eq) becomes larger than 1 and increases rapidly with a decrease in temperature. This drives the reaction Fe + H₂O ? H₂ + FeO to the right. With S = 100, the activity of FeO in the liquid droplet is 100 times as large as the value at equilibrium. The FeO activities are a function of temperature and provide relative average temperatures of the crystallization of chondrules. Our calculations for the LL3.0 chondrite Semarkona and our study of some non‐equilibrium effects lead to accurate representations of the compositions of chondrules of types IA, IAB, IB, IIA, IIAB, IIB, and CC. Our concepts readily explain both the variety of FeO concentrations in the different chondrule types and the entire process of chondrule formation. Our theory is unified and could possibly explain the formation of chondrules in all chondritic meteorites as well as provide a simple explanation for the complex chemistries of chondrites, and especially for type 3 chondrites.     

East Asian monsoon variation and climate changes in Jeju Island, Korea, during the latest Pleistocene to early Holocene

A 4.96-m-long sediment core from the Hanon paleo-maar in Jeju Island, Korea was studied to investigate the paleoclimatic change and East Asian monsoon variations during the latest Pleistocene to early Holocene (23,000-9000 cal yr BP). High-resolution TOC content, magnetic susceptibility, and major element composition data indicate that Jeju Island experienced the coldest climate around 18,000 cal yr BP, which corresponds to the last glacial maximum (LGM). Further, these multi-proxy data show an abrupt shift in climatic regime from cold and arid to warm and humid conditions at around 14,000 cal yr BP, which represents the commencement of the last major deglaciation. After the last major deglaciation, the TOC content decreased from 13,300 to 12,000 cal yr BP and from 11,500 to 9800 cal yr BP, thereby reflecting the weakening of the summer monsoon. The LGM in Jeju Island occurred later in comparison with the Chinese Loess Plateau. Such a disparity in climatic change events between central China and Jeju Island appears to be caused by the asynchrony between the coldest temperature event and the minimum precipitation event in central China and by the buffering effect of the Pacific Ocean.     

Glaciomarine sedimentation and palaeo-glacial setting of Maxwell Bay and its tributary embayment, Marian Cove, South Shetland Islands, West Antarctica

High-resolution (3.5 kHz and multi-channel) seismic profiles and piston cores were collected from Maxwell Bay and its tributary embayment, Marian Cove, in the South Shetland Islands, Antarctica, during the Korea Antarctic Research Program (1992/93 and 1995/96) to elucidate the glaciomarine sedimentation processes and recent glacial history of the area. Seismic data from Maxwell Bay reveal a rugged bay margin and flattened basin floor covered with well-stratified hemipelagic muds. On the base-of-slope, acoustically transparent debris flows occur, indicating downslope resedimentation of glaciomarine sediments. Despite the subpolar and ice-proximal settings of Marian Cove, the seafloor is highly rugged with a thin sediment drape, suggesting that much of the area has been recently eroded by glaciers. Sediment cores from the cove penetrated three distinct fining-upward lithofacies: (1) basal till in the lower part of the core, accumulated just seaward of the grounding line of the tidewater glacier; (2) interlaminated sand and mud in the middle part, deposited in ice-proximal zone by a combination of episodic subglacial meltwater inflow and iceberg dumping; and (3) pebbly mud in the upper part, deposited in ice-distal zone by both surface meltwater plume and ice-rafting from the glacier front. A reconstruction of the glacial history of these areas since the late glacial maximum shows an ice sheet filling Maxwell Bay in late Wisconsin time and grounding of the tidewater glacier in Marian Cove until about 1300 yr BP.     

Estimation of Contamination of ERS-2 and POSEIDON Satellite Radar Altimetry Close to the Coasts of Australia

It is broadly acknowledged that the precision of satellite-altimeter-measured instantaneous sea surface heights (SSH) is lower in coastal regions than in open oceans, due partly to contamination of the radar return from the coastal sea-surface state and from land topography. This study investigates the behavior of ERS-2 and POSEIDON altimeter waveform data in coastal regions and estimates a boundary around Australia's coasts in which the altimeter range may be poorly estimated by on-satellite tracking software. Over one million 20 Hz ERS-2 (March to April 1999) and POSEIDON (January 1998 to January 1999) radar altimeter waveform data were used over an area extending 350 km offshore Australia. The DS759.2 (5'resolution) ocean depth model and the GSHHS (0.2 km resolution) shoreline model were used together to define the coastal regions. Using the 50% threshold retracking points as the estimates of expected tracking gate, we determined that the sea surface height is contaminated out to maximum distance of between about 8 km and 22 km from the Australian shoreline for ERS-2, depending partly on coastal topography. Using the standard deviation of the mean waveforms as an indication of the general variability of the altimeter returns in the Australian coastal region shows obvious coastal contamination out to about 4 km for both altimeters, and less obvious contamination out to about 8 km for POSEIDON and 10 km for ERS-2. Therefore, ERS-2 and POSEIDON satellite altimeter data should be treated with some caution for distances less than about 22 km from the Australian coast and probably ignored altogether for distances less than 4 km.     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

Three‐dimensional simulation of tsunami run‐up around conical island

At the circular Babi Island in the Flores tsunami (1992) and pear shaped island in the Okushiri event (1993), unexpectedly large tsunami run‐up heights in the lee of conic islands were observed. The flume and basin physical model studies were conducted in the Coastal Hydraulic Laboratory, Engineering Research and Development Center, U.S. Army Corps of Engineers to provide a better understanding of the physical phenomena and verify numerical models used in predicting tsunami wave run‐up on beaches, islands, and vertical walls. Reasonably accurate comparison of run‐up height of solitary waves on a circular island has been obtained between laboratory experimental results and two‐dimensional computation model results. In this study we apply three‐dimensional RANS model to simulate wave run‐up on conical island. In the run‐up computation we obtain that 3D calculations are in very good comparison with laboratory and 2D numerical results. A close examination of the three‐dimensional velocity distribution around conical island to compare with depth‐integrated model is performed. It is shown that the velocity distribution along the vertical coordinate is not uniform: and velocity field is weaker in the bottom layer and higher on the sea surface. The maximum difference (about 40%) appears at the time when solitary wave reached the circular island.     

Source(s) and cycling of the nonhydrolyzable organic fraction of oceanic particles

A major fraction of particulate organic carbon (POC) in the deep ocean remains molecularly uncharacterized. In an effort to determine the chemical characteristics and source(s) of sinking POC, we studied a nonhydrolyzable fraction of sinking POC using ¹³C NMR (nuclear magnetic resonance) spectroscopy and analytical pyrolysis. ¹³C NMR spectra and products from analytical pyrolysis of the nonhydrolyzable fraction exhibit a strongly aliphatic character that is distinct from that of bulk POC. The aliphatic nature of this fraction is consistent with its low stable carbon isotope values. We hypothesize that the nonhydrolyzable fraction derives to a significant extent from a refractory component of organisms that selectively accumulates, resulting in its manifestation as a major part of POC sinking to the deep ocean and in underlying sediments.     

Study on the kinetics of iron oxide leaching by oxalic acid

The presence of iron oxides in clay or silica raw materials is detrimental to the manufacturing of high quality ceramics. Although iron has been traditionally removed by physical mineral processing, acid washing has been tested as it is more effective, especially for extremely low iron (of less than 0.1% w/w). However, inorganic acids such as sulphuric or hydrochloric acids easily contaminate the clay products with SO₄²⁻ and Cl⁻, and therefore should be avoided as much as possible. On the other hand, if oxalic acid is used, any acid left behind will be destroyed during the firing of the ceramic products. The characteristics of dissolution of iron oxides were therefore investigated in this study.The dissolution of iron oxides in oxalic acid was found to be very slow at temperatures within the range 25–60 °C, but its rate increases rapidly above 90 °C. The dissolution rate also increases with increasing oxalate concentration at the constant pH values set within the optimum range of pH2.5–3.0. At this optimum pH, the dissolution of fine pure hematite (Fe₂O₃) (105–140 μm) follows a diffusion-controlled shrinking core model. The rate expression expressed as 1 − (2 / 3)x − (1 − x)^2 / 3 where x is a fraction of iron dissolution was found to be proportional to [oxalate]^1.5.The addition of magnetite to the leach liquor at 10% w/w hematite was found to enhance the dissolution rate dramatically. Such addition of magnetite allows coarser hematite in the range 0.5–1.4 mm to be leached at a reasonable rate.