Pressure and temperature dependence of cation distribution in Mg-Mn olivine

Synthetic (Mg_0.51, Mn_0.49)₂SiO₄ olivine samples are heat-treated at three different pressures; 0, 8 and 12 GPa, all at the same temperature (～500° C). X-ray structure analyses on these single crystals are made in order to see the pressure effect on cation distribution. The intersite distribution coefficient of Mg and Mn in M1 and M2 sites, K _D = (Mn/Mg) _M1/(Mn/Mg) _M2, of these samples are 0.192 (0 GPa), 0.246 (8 GPa) and 0.281 (12 GPa), indicating cationic disordering with pressure. The small differences of cell dimensions between these samples are determined by powder X-ray diffraction. Cell dimensions b and c decrease, whereas a increases with pressure of equilibration. Cell volume decreases with pressure as a result of a large contraction of the b cell dimension. The effect of pressure on the free energy of the cation exchange reaction is evaluated by the observed relation between the cell volume and the site occupancy numbers. The magnitude of the pressure effect on cation distribution is only a fifth of that predicted from the observed change in volume combined with thermodynamic theory. This phenomenon is attributed to nonideality in this solid solution, and nonideal parameters are required to describe cation distribution determined in the present and previous experiments. We use a five-parameter equation to specify the cationic equilibrium on the basic of thermodynamic theory. It includes one energy parameter of ideal mixing, two parameters for nonideal effects, one volume parameter, and one thermal parameter originated from the lattice vibrational energy. The present data combined with some of the existing data are used to determine the five parameters, and the cation distribution in Mg-Mn olivine is described as a function of temperature, pressure, and composition. The basic framework of describing the cationic behavior in olivine-type mineral is worked out, although the result is preliminary: each of the determined parameters is not accurate enough to enable us to make a reliable prediction.     

Solid‐state 13C NMR characterization of insoluble organic matter from Antarctic CM2 chondrites: Evaluation of the meteoritic alteration level

Abstract— Chemical structures of the insoluble organic matter (IOM) from the Antarctic CM2 chondrites (Yamato [Y‐] 791198, 793321; Belgica [B‐] 7904; Asuka [A‐] 881280, 881334) and the Murchison meteorite were analyzed by solid‐state ¹³C nuclear magnetic resonance (NMR) spectroscopy. Different types of carbons were characterized, such as aliphatic carbon (Ali‐C), aliphatic carbon linked to hetero atom (Hetero‐Ali‐C), aromatic carbon (Aro‐C), carboxyls (COOR), and carbonyls (C=O). The spectra of the IOM from Murchison and Y‐791198 showed two major peaks: Ali‐C and Aro‐C, while the spectra from the other meteorites showed only one major peak of Aro‐C. Carbon distribution was determined both by manual integration and deconvolution. For most IOM, the Aro‐C was the most abundant (49.8–67.8%) of all carbon types. When the ratios of Ali‐C to Aro‐C (Ali/Aro) were plotted with the atomic hydrogen to carbon ratio (H/C), a correlation was observed. If we use the H/C as a parameter for the thermal alteration event on the meteorite parent body, this result shows a different extent of thermal alteration. In addition, IOM with a lower Ali/Aro showed a lower ratio of Ali‐C to COOR plus C=O (Ali / (COOR + C=O)). This result suggests that the ratio of CO moieties to aliphatic carbon in IOM might reflect chemical oxidation that was involved in hydrothermal alteration.     

Textural variations and impact history of the Millbillillie eucrite

Abstract— We have investigated 10 new specimens of the Millbillillie eucrite to study its textures and mineral compositions by electron probe microanalyser and scanning electron microscope. Although originally described as having fine-grained texture, the new specimens show diversity of texture. The compositions (Mg/Fe ratios) of the host pigeonites and augite lamellae are homogeneous, respectively, in spite of the textural variation. In addition to their chemical homogeneity, pyroxenes in coarse and fine-grained clasts are partly inverted to orthopyroxene. Chemical zoning of plagioclase during crystal growth is preserved. This eucrite includes areas of granulitic breccias and impact melts. Large scale textures show a subparallel layering suggesting incomplete mixing and deposition of impact melt and lithic fragments. An ³⁹Ar-⁴⁰Ar age determination for a coarse-grained clast indicates a strong degassing event at 3.55 ± 0.02 Ga. We conclude that Millbillillie is among the most equilibrated eucrites produced by thermal annealing after impact brecciation. According to the classification of impact breccias, Millbillillie can be classified as a mixture of granulitic breccias and impact melts. The last significant thermal event is characterized by network-like glassy veins that run through clasts and matrices. Consideration of textural observations and requirements for Ar-degassing suggests that the ³⁹Ar-⁴⁰Ar age could in principle date either the earilier brecciation and annealing event or the event which produced the veins.     

Bayesian and Neural Network Approaches to Estimate Deep Temperature Distribution for Assessing a Supercritical Geothermal System: Evaluation Using a Numerical Model

The temperature distribution at depth is a key variable when assessing the potential of a supercritical geothermal resource as well as a conventional geothermal resource. Data-driven estimation by a machine-learning approach is a promising way to estimate temperature distributions at depth in geothermal fields. In this study, we developed two methodologies—one based on Bayesian estimation and the other on neural networks—to estimate temperature distributions in geothermal fields. These methodologies can be used to supplement existing temperature logs, by estimating temperature distributions in unexplored regions of the subsurface, based on electrical resistivity data, observed geological/mineralogical boundaries, and microseismic observations. We evaluated the accuracy and characteristics of these methodologies using a numerical model of the Kakkonda geothermal field, Japan, where a temperature above 500 °C was observed below a depth of about 3.7 km. When using geological and geophysical knowledge as prior information for the machine learning methods, the results demonstrate that the approaches can provide subsurface temperature estimates that are consistent with the temperature distribution given by the numerical model. Using a numerical model as a benchmark helps to understand the characteristics of the machine learning approaches and may help to identify ways of improving these methods.

     

Seismic response of steel structures with seesaw systems using viscoelastic dampers

Vibration control systems are being used increasingly worldwide to provide enhanced seismic protection for new and retrofitted buildings. This paper presents a new vibration control system on the basis of a seesaw mechanism with viscoelastic dampers. The proposed vibration control system comprises three parts: brace, seesaw member, and viscoelastic dampers. In this system, only tensile force appears in bracing members. Consequently, the brace buckling problem is negligible, which enables the use of steel rods for bracing members. By introducing pre‐tension in rods, long steel rods are applicable as bracing between the seesaw members and the moment frame connections over some stories. Seesaw mechanisms can magnify the damper deformation according to the damper system configuration. In this paper, first, the magnification factor, that is, the ratio of the damper deformation to the story drift, is delivered, which includes the rod deformation. Results of a case study demonstrate that the magnification factor of the proposed system is greater than unity for some cases. Seismic response analysis is conducted for steel moment frames with the proposed vibration control system. Energy dissipation characteristics are examined using the time‐history response results of energy. The maximum story drift angle distributions and time‐history response results of displacement show that the proposed system can reduce the seismic response of the frames effectively. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal variations in the nitrogen isotopic composition of settling particles at station K2 in the western subarctic North Pacific

Intensive observations using hydrographical cruises and moored sediment trap deployments during 2010 and 2012 at station K2 in the North Pacific Western Subarctic Gyre (WSG) revealed seasonal changes in δ ¹⁵N of both suspended and settling particles. Suspended particles (SUS) were collected from depths between the surface and 200 m; settling particles by drifting sediment traps (DST; 100–200 m) and moored sediment traps (MST; 200 and 500 m). All particles showed higher δ ¹⁵N values in winter and lower in summer, contrary to the expected by isotopic fractionation during phytoplankton nitrate consumption. We suggest that these observed isotopic patterns are due to ammonium consumption via light-controlled nitrification, which could induce variations in δ ¹⁵N(SUS) of 0.4–3.1 ‰ in the euphotic zone (EZ). The δ ¹⁵N(SUS) signature was reflected by δ ¹⁵N(DST) despite modifications during biogenic transformation from suspended particles in the EZ. δ ¹⁵N enrichment (average: 3.6 ‰) and the increase in C:N ratio (by 1.6) in settling particles suggests year-round contributions of metabolites from herbivorous zooplankton as well as TEPs produced by diatoms. Accordingly, seasonal δ ¹⁵N(DST) variations of 2.4–7.0 ‰ showed a significant correlation with primary productivity (PP) at K2. By applying the observed δ ¹⁵N(DST) vs. PP regression to δ ¹⁵N(MST) of 1.9–8.0 ‰, we constructed the first annual time-series of PP changes in the WSG. This new approach to estimate productivity can be a powerful tool for further understanding of the biological pump in the WSG, even though its validity needs to be examined carefully.     

Resonance capacity method for earthquake response analysis of hysteretic structures

The Resonance Capacity Method is proposed for the earthquake response analysis of hysteretic structures. Resonance Capacity is a physical quantity of structures which is related to the hysteretic energy absorbed by structures in one cycle and is equated to the acceleration, velocity and displacement amplitudes α₀, d₀ and d₀ of earthquake ground motions at resonance.¹ According to the idealized trapezoidal approximation of earthquake ground motions in the logarithmic period–velocity plane as proposed by Veletsos and Newmark,⁸ the Resonance Capacity property applies in each period range, short, medium and long, where α₀, v₀ and d₀ respectively are approximately constant. In the medium range of periods, the energy dissipated in hysteretic loops and the deformation amplitudes of a single-degree system with elasto–plastic force–deformation relationships are calculated for the case of El Centro 1940, 18 May earthquake, by this Resonance Capacity Method. The result is compared with results from conventional numerical response analyses obtained by Berg and Thomaides,¹⁴ Kato and Akiyama¹² and Veletsos and Newmark,⁸ and the general agreement is seen to be good. Therefore, it may be possible to apply this Resonance Capacity Method over the entire range of periods. By means of this method the earthquake response analysis of hysteretic systems can be performed easily, and the hysteretic energy and fatigue characteristics of structures may be taken into account directly, up to the point of fracture.     

Near inertial motion excited by wind change in a margin of the Typhoon 9019

An excitation of inertial oscillation in the upper layer east of course of Typhoon 9019 was fortuitously observed at three surface buoys deployed during the Ocean Mixed Layer Experiment (OMLET). The observed inertial oscillation was compared with wind fluctuation measured at Ocean Weather Station T (29°N, 135°E) which was placed at the center of a triangle with three vertexes occupied by the respective surface buoys. Inertial oscillation is effectively excited in the mixed layer at the eastern margin of the typhoon by a rapid decrease of wind rather than by prevailing strong wind. It is shown by means of a least square deviation that the inertial oscillation observed in the mixed layer has a period of 23.9 hours shorter than the local inertial period of 24.7 hours. This shorter period suggests that the inertial oscillation has the finite velocities of phase and group as an inertial internal wave. A theoretically obtained ratio of vertical component of group velocity to that of phase velocity, approximately agrees with observed value. The inertial internal wave is excited by fluctuation of divergence with near inertial period in the mixed layer.     

SAR-imaged spiral eddies in Mutsu Bay and their dynamic and kinematic models

A pattern of slick streaks winding into a spiral, known as a spiral eddy, was identified in 5 images taken by the ERS-1/2 synthetic aperture radar (SAR) in Mutsu Bay (Japan); dynamic and kinematic models of these spiral eddies have been proposed. Common characteristics of the five spiral eddies are: 1) an eddy diameter of about 15 km; 2) their location in the western part of the bay; and 3) their cyclonic direction of rotation. Moreover, the wind conditions over the bay were common: prior to acquiring the images, a strong easterly wind continued blowing for more than one day. The wind field on the bay is known to be orographically steered and has strong windstress vorticity, which generates cyclonic circulation. The diameter and location of the circulation simulated with a numerical ocean model corresponded well to those of the identified spiral eddies. Based on these facts, we propose a dynamic model for the movement of a slick streak, and a kinematic model for the formation of a spiral eddy. We have assumed calm air, a microlayer and seawater with a cyclonic circulation in the dynamic model. The balance of forces is established in the microlayer among the frictional force from the seawater, the frictional force from the calm air, the gravitational force, and the Coriolis force. As a result, the velocity vector of the microlayer deflects slightly towards the center of the cyclonic circulation. We have assumed a point source of the microlayer in the kinematic model. The shapes of a slick streak simulated with the models agree well with the identified patterns in the SAR images.     

Possibility of recent changes in vertical distribution and size composition of chlorophyll-a in the western North Pacific region

Our analysis of the last three decades of retrospective data of vertical distributions and size composition of chlorophyll-a (Chl-a) over the western North Pacific has revealed significant changes of three indices related to Chl-a during summer season, as follows: (1) decreasing linear trend of the proportion of Chl-a in surface layer to that of the whole water column by 0.4 and 2.3% year⁻¹ in the subtropical area along 137°E (STA₁₃₇) during 1972 to 1997 and in the Kuroshio Extension area along 175°E (KEA₁₇₅) during 1990 to 2001; (2) increasing linear trend of the depth of subsurface Chl-a maximum (DCM) by 0.4 and 2.6 m year⁻¹ in STA₁₃₇ and KEA₁₇₅; and (3) decreasing linear trend of larger-size Chl-a (>3 μm) by 0.1 and 2.5% year⁻¹ in STA₁₃₇ and KEA₁₇₅, respectively. Water density (σ _θ ) at 75 m depth had also decreased by 0.006 and 0.05 year⁻¹ in STA₁₃₇ and KEA₁₇₅, respectively. The ratio of biogenic opal to biogenic CaCO₃ in the sinking flux decreased by 0.015 year⁻¹ in the subtropical region from 1997 to 2005. These findings may indicate that the subsurface chlorophyll maximum is deepening and larger phytoplankton such as diatoms has been decreasing during the past decade, associated with the decreasing density of surface water caused by warming in the western North Pacific, especially in the summer.