Kinetic stability of a melted iron globule during chondrule formation. I. Non‐rotating model

Abstract— We have investigated the kinematics of the separation of iron globules from chondrules during chondrule formation. A simple model, which assumes that the system has no angular momentum, was used to calculate the energy of a system with an iron globule and a chondrule. The energies of three different states were calculated: 1) a melted iron globule fully embedded in a melted chondrule, 2) a melted iron globule on the surface of a melted chondrule, and 3) a melted iron globule being separated from a melted chondrule. We also calculated the lowest energy shape for a melted iron globule on the surface of a melted chondrule, and compared our result with the shapes of four natural samples of chondrules and iron globules in thin sections. The shapes were calculated using an assumed value for the interface energy between the four couples of melted chondrules and the iron globules, and agree well with the natural shapes of chondrules and iron globules. The results of our calculations show that the iron globules of these four samples would be strongly bound to the surface of the melted chondrule during chondrule formation, and separation would be difficult, if the iron globules had been on the surface of precursors of these chondrules. Our results also show that if these iron globules were initially inside and transported to the surface of the melted chondrule, most of them would be ejected from the inside to outside because of surface tension forces, as long as the energy losses due to viscous dissipation when the globules pass through the surface of melted chondrules were sufficiently small. Although further improvement of the model is required, our results demonstrate that this ejection process may be responsible for the depletion of siderophile elements in natural chondrules.     

Verifying optimality of rainfed agriculture using a stochastic model for drought occurrence

It may be paradoxical but subsistence rainfed agriculture is the predominant source of food in Sub-Saharan Africa where the production uncertainty is associated with the stochastic nature of rainfall. This paper attempts to comprehend the rationale of this situation by a mathematical approach. Considering the level of drought severity as the zero-reverting Ornstein–Uhlenbeck process, optimality of rainfed agriculture is investigated in the context of stochastic control theory. Occurrence of drought terminating growth of crops is modelled with the concept of first exit time. A stochastic control problem allowing for virtual cost of irrigation, water stress to crops, and benefits of farming is formulated with irrigation effort as the control variable. The Hamilton–Jacobi–Bellman equation governing the optimal control is studied to identify the set of cost functions optimizing rainfed agriculture in an inverse problem approach. Data and information were collected in the coastal savanna agro-ecological zone of Ghana, to identify model parameters, formulate the stochastic control problem, solve the inverse problem, and then verify optimality of rainfed agriculture. The results indicated that rainfed agriculture is not optimal when the crop is more tolerant to water stress.     

High-frequency polarity swings during the Gauss-Matuyama reversal from Baoji loess sediment

Paleomagnetic records of the Gauss-Matuyama reversal were obtained from two loess sections at Baoji on the Chinese Loess Plateau. Stepwise thermal demagnetization shows two obvious magnetization components. A low-temperature component isolated between 100 and 200–250°C is close to the present geomagnetic field direction, and a high-temperature component isolated above 200–250°C reveals clearly normal, reversed, and transitional polarities. Magnetostratigraphic results of both sections indicated that the Gauss-Matuyama reversal consists of a high-frequency polarity fluctuation zone, but the characteristic remanent magnetization directions during the reversal are clearly inconsistent. Rock magnetic experiments demonstrated that for all the specimens with normal, reversed, and transitional polarities magnetite and hematite are the main magnetic carriers. Anisotropy of magnetic susceptibility indicates that the studied loess sediments have a primary sedimentary fabric. Based on virtual geomagnetic pole latitudes, the Gauss-Matuyama reversal records in the two sections are accompanied by 14 short-lived geomagnetic episodes (15 rapid polarity swings) and 12 short-lived geomagnetic episodes (13 rapid polarity swings), respectively. Our new records, together with previous ones from lacustrine, marine, and aeolian deposits, suggest that high-frequency polarity swings coexist with the Gauss-Matuyama reversal, and that the Gauss-Matuyama reversal may have taken more than 11 kyr to complete. However, we need more detailed analyses of sections across polarity swings during reversals as well as more high-resolution reversal records to understand geomagnetic behavior and inconsistent characteristic remanent magnetization directions during polarity reversals.     

The Impact of Orographically-Induced Gravity Wave on the Diurnal Cycle of Rainfall over Southeast Kalimantan Island

Precipitation measurements from the Tropical Rainfall Measuring Mission (TRMM) satellite indicate that the southeastern area of Kalimantan (Borneo) Island receives much less rainfall than elsewhere on the island during the period from July to October.Results from sur-face meteorological observations show that the diurnal cycle of rainfall differs greatly between the eastern and western coasts of the island.Rainfall on the western coast of the island is frequent in the afternoon and evening,whereas almost all rainfall on the eastern coast occurs in the morning.Meanwhile,the Global Positioning System (GPS)-derived precipitable water (PW) on the eastern coast shows a substantial decrease in moisture in the af-ternoon and evening.Numerical experiments with a mesoscale model reveal that gravity waves driven by di-urnal heating of the elevated land surface of the moun-tains on Sulawesi Island,which lies approximately 300 kilometers to the east of Kalimantan Island,significantly affect the diurnal cycle of rainfall over southeast Kali-mantan Island.     

Gold Mineralization in Izu Peninsula,Central Japan,during Crustal Extension in Response to Double Subduction

Izu Peninsula in central Japan, the northern tip of the Izu‐Bonin arc, hosts numerous epithermal Au–Ag vein deposits of low‐sulfidation style. All have similar vein textures, mineralogy, and alteration. Geochemical data from fluid inclusions in vein quartz, the mineralogy and mineral chemistry of alteration, and stable isotope data indicate that auriferous hydrothermal activity occurred under subaerial conditions. The K–Ar ages of auriferous vein minerals are <1.5 Ma, indicating that the mineralization took place after extensive submarine volcanism for the host rocks. These observations suggest that Au–Ag mineralization was synchronous with the development of an extensional regime of the Izu block after its collision with the Honshu arc after 1.5 Ma. This collision resulted in the shifting of the Izu block far from the trench to the rear position, and the subduction of the Izu block along the Suruga trough to the west and along the Sagami trough to the east. The reararc position of the Izu block and double subduction resulted in crustal extension, upwelling of asthenospheric mantle, and tholeiitic magmatism reflected by mafic dyke swarms and subsequent monogenetic volcanic activity in the Izu peninsula. The timing of the Au mineralization in the Izu Peninsula during the beginning of lithospheric extension is similar to that of the Sado Au–Ag deposit on Sado island in the Japan Sea. Two mineralization events coincide with extensive tholeiitic mafic volcanism and injections of dyke swarms related to the back‐arc opening of the Japan Sea. The geological setting of the Au–Ag mineralization in Izu and Sado is also similar to that of the epithermal Au–Ag deposits in northern Nevada, where mineralization was contemporaneous with crustal extension and tholeiitic mafic magmatism derived from the asthenospheric mantle. This study suggests that epithermal Au mineralization at shallow crustal depths is a product of large‐scale lithospheric evolution.     

A deep section of accretionary complex: Susunai Complex in Sakhalin Island, Northwest Pacific Margin

Abstract A deep section of accretionary complex, the metamorphosed Susunai Complex, is observed on Sakhalin Is., Russia. High pressure part of pumpellyite-actinolite facies metavolcanics, metacherts and metapelites are well exposed and constitute a tectonic pile preserving primary structures related to underplating of the oceanic crust. Three stages of deformation, D₁ through D₃, suggest successive deformation during subduction, underplating and exhumation of the complex. Oceanic material in the complex is more abundant than other well documented ancient accretionary complexes (e.g. the Shimanto Belt in southwest Japan and the Ghost Rocks Formation in Alaska), which were shallowly underplated. At Susunai, deep down-stepping of a décollément has scraped off the upper part of the oceanic crust, primarily the pillowed basalt horizon. This down-stepping results from crustal weakening as overpressured water is released from the fractured oceanic crust during metamorphism.     

The mesozoic evolution of the mino terrane, central Japan: A geologic and paleomagnetic synthesis

The Mino tectono-stratigraphic terrane, central Japan, underlain by Permian to Jurassic sedimentary and volcanic rocks of various origins, was formed through accretion processes associated with the Mesozoic sea-floor spreading. This conclusion has been reached mainly from the following reasoning:

1. (1) the entire boundary of this terrane is defined by tectonic belts with high-pressure metamorphic rocks and serpentinized ultramafic rocks,

2. (2) the chemistry and petrology of the Permian greenstones demonstrate their affinity with abyssal tholeiitic and alkalic basalts,

3. (3) the widespread, but chaotic, occurrence of Permian greenstones, Triassic cherts, and Jurassic siliceous shales in the younger Jurassic clastic rocks of this terrane suggests extensive post-depositional mixing of strata,

4. (4) the sedimentology of the Jurassic sandstones strongly suggests that they are turbidity-current deposits supplied from cratonic lands,

5. (5) the South-Pacific type fossil assemblage in the Mino terrane shows strong contrast with the North-Pacific type fossil assemblage of the adjacent terranes,

6. (6) the paleomagnetism of the Permian and Jurassic greenstones, the Triassic cherts, and the Jurassic siliceous shales implies long-distance northward drift in Cretaceous time of these rocks from their original low latitudinal regions.

Along with this northward migration, the Mino terrane was accreted with extensive internal deformation to northeast Asia including the present Hida terrane. Recent accumulation of paleomagnetic and paleontologic data in the Pacific peripheral regions appears to support the existence of many allochthonous terranes which migrated from the equatorial regions. The Mino terane may be regarded as one example of these Circum-Pacific allochthons.     

A scaling model for quantification of earthquakes in and near Japan

A simple method is developed to determine seismic moments of earthquakes. The method is qualified through criteria such as simplicity of calculations, coverage of wide magnitude range, and insensitivity to detailed instrumental response. The method is applied to 163 major earthquakes which occurred underneath Japan and the Japan Sea in the time from 1926 to 1977. Magnitudes of these earthquakes, which have been determined by the Japan Meteorological Agency, (M_JMA) cover the range from 4.3 to 7.5. At first, source spectra are analyzed through a very simple way introducing two new parameters: characteristic period T_c and seismic-moment factor M_c. The former is defined as an average value of apparent periods of seismic waves with the maximum trace amplitude at many stations. The latter is an average of products of maximum trace amplitude and its apparent period multiplied by epicentral distance. It is shown that T_c corresponds to the period of the corner frequency of an earthquake and M_c to the seismic-moment density at the period of T_c. A scaling model of earthquake source spectra is presented which satisfies the empirical relations between the surface-wave magnitude M_s and M_JMA, and M_JMA and the body-wave magnitude m_b. Those relations are independent of the Gutenberg and Richter relation between M_s and m_b, because M_JMA is determined from maximum amplitudes of seismic waves with a period of about 4 sec. The static seismic moment of each earthquake can be estimated from calculated M_c using the source spectra of the scaling model. Seismic moments of 18 earthquakes determined by conventional analyses from near- and/or far-field observations are consistent with static seismic moments thus estimated over the range from 2 × 10²³ to 3 × 10²⁷ dyne cm. This shows the potential in practice of the present method, especially in the routine processing of seismic data.     

Stochastic analysis of bed-thickness distribution of sediments

On formation of a bed and distribution of bed thickness, A. N. Kolmogorov presented a mathematical explanation that if repetitive alternations of material accumulation and erosion form a sequence of beds, the resultant bed-thickness distribution curve takes a shape truncated by the ordinate at zero thickness. In this truncated distribution curve, its continuation and extension from positive to negative thickness represents the distribution of beds with negative thickness, that is, the depth of erosion. When a distribution curve, including both positive and negative parts, is expressed by a function f(x),the ratio \(\int_0^\infty {f(x)dx to} \int_{ - \infty }^\infty {f(x)dx} \) ,called Kolmogorov's coefficient and designated as p,is a parameter representing the degree of accumulation in the depositional environment. On the assumption that f(x)is described by the Gaussian distribution function, the coefficient pfor Permian and Pliocene sequences in central Japan was calculated. The coefficients also were obtained from published data for different types of sediments from other areas. It was determined that they are more or less different depending on their depositional environments. The calculated results are summarized as follows: $$\begin{gathered} p = 0.80 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ p = 0.65 - 0.95for{\text{ }}nearshore{\text{ }}sediments \hfill \\ p = 0.55 - 0.95for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ p = 0.90 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ In addition, a ratio \(q = \int_0^\infty {xf(x)dx/} \int_{ - \infty }^\infty {|x|f(x)dx} \) ,called Kolmogorov's ratio in this paper, is introduced for estimating a degree of total thickness actually observed in the field relative to total thickness once present in a basin. The calculated results of Kolmogorov's ratio are as follows: $$\begin{gathered} q = 0.88 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ q = 0.68 - 0.98for{\text{ }}nearshore{\text{ }}sediments \hfill \\ q = 0.55 - 0.96for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ q = 0.92 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ The sedimentological significance of these values is discussed.     

Contributions from mafic alkaline magmas to the Bingham porphyry Cu–Au–Mo deposit, Utah, USA

The Bingham porphyry Cu-Au-Mo deposit, Utah, may only be world-class because of substantial contributions of sulfur and metals from mafic alkaline magma to an otherwise unremarkable calc-alkaline system. Volcanic mafic alkaline rocks in the district are enriched in Cr, Ni, and Ba as well as Cu, Au, platinum group elements (PGE), and S. The bulk of the volcanic section that is co-magmatic with ore-related porphyries is dacitic to trachytic in composition, but has inherited the geochemical signature of high Cr, Ni, and Ba from magma mixing with the mafic alkaline rocks. The volcanic section that most closely correlates in time with ore-related porphyries is very heterogeneous containing clasts of scoriaceous latite, latitic, and minette, and flows of melanephelinite, shoshonite, and olivine latite in addition to volumetrically dominant dacite/trachyte. Bingham ore-related porphyries show ample evidence of prior mixing with mafic alkaline magmas. Intrusive porphyries that have not been previously well-studied have several chemical and mineralogical indications of magma mixing. These "mixed" lithologies include the hybrid quartz monzonite porphyry, biotite porphyry, and minette dikes. Even some of the more silicic latite and monzonite porphyries retain high Cr and Ba contents indicative of mixing and contain trace amounts of sapphire (<1 mm). The heterogeneous block and ash flow deposits also contain sapphire and are permissively correlated with the intrusions based on chemical, mineralogical, and isotopic data. Magma mixing calculations suggest about 10% of the monzonitic/latitic ore-related magma may have been derived from mafic alkaline magma similar to the melanephelinite. If the original S content of the mafic magma was about 2,000-4,000 ppm, comparable with similar magmas, then the mafic magma may have been responsible for contributing more than half of the S and a significant portion of the Cu, Au, and PGE in the Bingham deposit.