Viscosity of peridotite liquid up to 13 GPa: Implications for magma ocean viscosities

The viscosity of synthetic peridotite liquid has been investigated at high pressures using in-situ falling sphere viscometry by combining a multi-anvil technique with synchrotron radiation. We used a newly designed capsule containing a small recessed reservoir outside of the hot spot of the heater, in which a viscosity marker sphere is embedded in a forsterite + enstatite mixture having a higher solidus temperature than the peridotite. This experimental setup prevents spheres from falling before a stable temperature above the liquidus is established and thus avoids difficulties in evaluating viscosities from velocities of spheres falling through a partially molten sample.

Experiments have been performed between 2.8 and 13 GPa at temperatures ranging from 2043 to 2523 K. Measured viscosities range from 0.019 (± 0.004) to 0.13 (± 0.02) Pa s. At constant temperature, viscosity increases with increasing pressure up to  8.5 GPa but then decreases between  8.5 and 13 GPa. The change in the pressure dependence of viscosity is likely associated with structural changes of the liquid that occur upon compression. By combining our results with recently published 0.1 MPa peridotite liquid viscosities [D.B. Dingwell, C. Courtial, D. Giordano, A. Nichols, Viscosity of peridotite liquid, Earth Planet. Sci. Lett. 226 (2004) 127–138.], the experimental data can be described by a non-Arrhenian, empirical Vogel-Fulcher-Tamman equation, which has been modified by adding a term to account for the observed pressure dependence of viscosity. This equation reproduces measured viscosities to within 0.08 log₁₀-units on average. We use this model to calculate viscosities of a peridotitic magma ocean along a liquid adiabat to a depth of  400 km and discuss possible effects on viscosity at greater pressures and temperatures than experimentally investigated.     

High-reliability zircon separation for hunting the oldest material on Earth: An automatic zircon separator with image-processing/microtweezers-manipulating system and double-step dating

Despite the recent development in radiometric dating of numerous zircons by LA-ICPMS, mineral separation still remains a major obstacle, particularly in the search for the oldest material on Earth. To improve the efficiency in zircon separation by an order of magnitude, we have designed/developed a new machine-an automatic zircon separator(AZS). This is designed particularly for automatic pick-up of100 μm-sized zircon grains out of a heavy mineral fraction after conventional separation procedures. The AZS operates in three modes:(1) image processing to choose targeted individual zircon grains out of all heavy minerals spread on a tray,(2) automatic capturing of the individual zircon grains with microtweezers, and(3) placing them one-by-one in a coordinated alignment on a receiving tray. The automatic capturing was designed/created for continuous mineral selecting without human presence for many hours. This software also enables the registration of each separated zircon grain for dating, by recording digital photo-image, optical(color) indices, and coordinates on a receiving tray. We developed two new approaches for the dating; i.e.(1) direct dating of zircons selected by LA-ICPMS without conventional resin-mounting/polishing,(2) high speed U-Pb dating, combined with conventional sample preparation procedures using the new equipment with multiple-ion counting detectors(LA-MIC-ICPMS).With the first approach, Pb-Pb ages obtained from the surface of a mineral were crosschecked with the interior of the same grain after resin-mounting/polishing. With the second approach, the amount of time required for dating one zircon grain is ca. 20 s, and a sample throughput of 150 grains per hour can be achieved with sufficient precision(ca. 0.5%).We tested the practical efficiency of the AZS, by analyzing an Archean Jack Hills conglomerate in Western Australia with the known oldest(4.3 Ga) zircon on Earth. Preliminary results are positive; we were able to obtain more than 194 zircons that are over 4.0 Ga out of ca. 3800 checked grains, and 9 grains were over 4300 Ma with the oldest at 4371 ± 7 Ma. This separation system by AZS, combined with the new approaches, guarantees much higher yield in the hunt for old zircons.     

Cause and countermeasure of long-period oscillations of moored ships and the quantification of surge and heave amplitudes

Long-period oscillations of moored ships whose periods are about 1 or 2 min cause many troubles in many ports and harbours. It is necessary to investigate these phenomena and verify their causes and countermeasures in each case because they are strongly dependent on the environment of each port and harbour. From this point of view, long-period oscillations of moored ships in the Port of Shibushi in Japan were investigated by means of wave observations, the image processing of moored ship motions using the video camera and motion-capture software and numerical simulations. From observation results, the relationship between offshore long-period waves and long-period oscillations of moored ships was recognized and surge and heave amplitudes were quantified by using wave data in order to forecast moored ship motions. Furthermore, from observation and numerical results, it was revealed that long-period waves with the peak period of 120 s from the offshore typhoon kept or exaggerated the local harbour oscillation of 60–70 s and it caused long-period oscillations of moored ships. Numerical results in case of reducing the reflection coefficient of the target berth implied that it ceased the local harbour oscillation and it would give an effective countermeasure to reduce long-period oscillations of moored ships in the Port of Shibushi.     

Heterogeneous shock-induced thermal radiation in minerals

A 500 channel optical imaging intensifying and spectral digital recording system is used to record the shock-induced radiation emitted from 406 to 821 nm from transparent minerals during the time interval that a shock wave propagates through the sample. Initial results obtained for single crystals of gypsum, calcite and halite in the 30 to 40 GPa (300 to 400 kbar) pressure range demonstrate greybody emission spectra corresponding to temperatures in the 3,000 to 4,000 K range and emissivities ranging from 0.003 to 0.02. In the case of gypsum and calcite, distinctive line spectra, are superimposed on the thermal radiation. The observed color temperatures are a factor of 2 to 10 greater than the Hugoniot temperature, calculable on the basis of continuum thermodynamics and equation of state models for the shock states achieved in the three minerals. These observed high temperatures are believed to be real. We conclude that we are detecting a large number of closed spaced high temperature shear-band regions immediately behind the shock front. A shear instability model, such as proposed independently by Grady (1977, 1980), Ananin et al. (1974), and Horie (1980), in which small zones of highly deforming and melted material are produced which are the source of the detected thermal radiation and have a fractional effective area approximately numerically equal to the measured emissivity, can be used to predict an effective emissivity of 0.0065 directly behind the shock front. If shear band instability arises from viscous flow processes, Grady's model and mineral thermal properties yield for the shocked mineral viscosities values in the range of 10⁹ to 10¹⁵ P immediately behind the shock front.     

Diabatic heating and the low frequency dynamics in the tropics

Summary Using a General Circulation Model developed at FSU (FSUGCM), the role of the diabatic heating on the 30–60-day oscillation is investigated. To concentrate on the radiation and the moist convection processes, an aqua planet model is employed in this paper. We have obtained a 40-day oscillation with relatively lower frequency than other GCMs without strong heating in the lower troposphere. Unlike some GCMs and simple models, the convective area does not move eastward along with the oscillation. Adiabatic cooling due to the upward motion is mostly compensated by diabatic heating. This implies that Kelvin CISK theory might not explain our 40-day oscillation. We have also examined the impact of radiative heating on the low frequency oscillation. When we reduce the radiative cooling rate, our 40-day mode does not appear and a Kelvin CISK mode appears with a faster phase speed. The impact of the different convection schemes is also investigated. With an enhanced convection scheme, zonal wave number two with a 40-day period is generated.With 12 FiguresOn leave from Japan Meteorological Agency.     

Organic geochemistry of the Japan Sea sediments-1: Bulk organic matter and hydrocarbon analyses of Core KH-79-3, C-3 from the Oki Ridge for paleoenvironment assessments

Organic geochemical study of bulk organic matter (OM), hopanoid hydrocarbon and normal hydrocarbon (C₂₃C₃₅) was conducted for a 936-cm-long sediment core sample from the Oki Ridge of the Japan Sea (Core KH-79-3, C-3; 37°03.5 N, 134°42.6E, water depth 935 m). Stable carbon isotopic ratios were also measured for both bulk OM and individual hydrocarbons. The following results were obtained: (1) The weight ratios of total organic carbon to total nitrogen range from 6.2 to 9.4 in the core. The ¹³C values of bulk OM range from –25.1–20.7%.. The ¹³C values of OM in the sections of 140190 cm are lower (–25–24) than those in the other sections (–23–21). This result indicates that OM in the core except for the 140190 cm sections is essentially of marine origin. (2) The ¹³C value of diploptene (a hopanoid hydrocarbon) in the last glacial maximum (LGM), is –66.3 (vs. PDB), which indicates it originating in methanotrophic bacteria. This result provides evidence to support for the previous ideas (Oba et al., 1980, 1984; Masuzawa and Kitano, 1984) that the bottom waters in the Japan Sea were anoxic in LGM. (3) Long chain (C₂₃C₃₅) n-alkanes of higher-plant wax origin were found throughout the core. Their concentration is high in 140190 cm in depth, suggesting that eolian dust load was high in LGM. (4) The n-alkane/TOC ratio increases with decreasing ¹³C values of bulk OM. This result indicates that the load of terrestrial (probably eolian dust-derived) OM to the Japan Sea became higher in colder climates. (5) The CPI values of long-chain n-alkanes are different in different ¹³O stages of paleoclimate, probably reflecting variations in species of terrestrial higher plants as a result of climatological adaptations.     

Paleomagnetism of the middle Cretaceous Iritono granite in the Abukuma region, northeast Japan

We have studied the paleomagnetism of the middle Cretaceous Iritono granite of the Abukuma massif in northeast Japan together with ⁴⁰Ar–³⁹Ar dating. Paleomagnetic samples were collected from ten sites of the Iritono granite (102 Ma ⁴⁰Ar–³⁹Ar age) and two sites of its associated gabbroic dikes. The samples were carefully subjected to alternating field and thermal demagnetizations and to rock magnetic analyses. Most of natural remanent magnetizations show mixtures of two components: (1) H component, high coercivity (B_c > 50–90 mT) or high blocking temperature (T_b > 350–560 °C) component and (2) L component, relatively low B_c or low T_b component. H component was obtained from all the 12 sites to give a mean direction of shallow inclination and northwesterly declination (I = 29.9°, D = 311.0°, α₉₅ = 2.7°, N = 12). This direction is different from the geocentric axial dipole field at the present latitude (I = 56.5°) and the typical direction of the Cenozoic remagnetization in northeast Japan. Since rock magnetic properties indicate that the H component of the Iritono granite is carried mainly by magnetite inclusions in plagioclase, this component probably retains a primary one. Thus the shallow inclination indicates that the Abukuma massif was located at a low latitude (16.1 ± 1.6°N) about 100 Ma and then drifted northward by about 20° in latitude. The northwesterly deflection is attributed mostly to the counterclockwise rotation of northeast Japan due to Miocene opening of the Japan Sea. According to this model, the low-pressure and high-temperature (low-P/high-T) metamorphism of the Abukuma massif, which has been well known as a typical location, would have not occurred in the present location. On the other hand, the L component is carried mainly by pyrrhotite and its mean direction shows a moderate inclination and a northwesterly declination (I = 42.8°, D = 311.5°, α₉₅ = 3.3°, N = 9). Since this direction is intermediate between the H component and early Cenozoic remagnetization in northeast Japan, some thermal event would have occurred at lower temperature than pyrrhotite Curie point ( 320 °C) during the middle Cretaceous to early Cenozoic time to have resulted in partial remagnetization.     

Roughness length for heat over an urban canopy

The roughness length for heat z_T was evaluated over an urban canopy, using the measured sensible heat flux and radiometric temperature. To overcome thermal heterogeneity in the urban area, the measured radiometric temperature was transformed into the equivalent temperature of an upward longwave radiation flux. The equivalent temperature was found to provide an effective parameterization of the radiometric temperature. The daytime average of the resulting ln(z_T/z₀) was 10, where z₀ is the aerodynamic roughness length. This result generally agrees with previous studies; however, the anthropogenic heat is a large uncertainty, which could cause an error at least 240% in z_T.     

Thermal equation of state of majorite with MORB composition

In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V₀ = 1574.14(4) Å³) and the pressure derivative of the isothermal bulk modulus (K^′T = 4). This yielded an isothermal bulk modulus of K_T0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.022(5) GPa K⁻¹, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10⁻⁵ K⁻¹ and b = 1.0(5) × 10⁻⁸ K⁻². The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.