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.     

Coastal upwelling along the north coast of Papua New Guinea and SST cooling over the pacific warm pool: A case study for the 2002/03 El Niño event

We investigate an overlooked mechanism—coastal upwelling—for sea surface temperature (SST) cooling in the western side of the mean location of the Pacific warm pool (WSWP: 5°S–5°N, 140°E–150°E) prior to El Niño onset. We analyze various observed data such as the TRIangle Trans-Ocean buoy Network (TRITON) moored buoy data, Conductivity-Temperature-Depth (CTD) data, satellite data and a hindcast experiment output by a high-resolution ocean general circulation model (OGCM). We focus on the precondition of the 2002/03 El Niño event, for which many datasets are available. Relatively cool water upwelled along the north coast of Papua New Guinea (PNG) during December 2001, prior to the onset of the 2002/03 El Niño event, and then spread out over a wider area to the northeast. Simultaneously, strong west-northerly surface winds occur along the north coast. Heat budget analysis of TRITON buoy data in the WSWP reveals that negative zonal heat advection due to eastward current is the main factor for cooling the mixed layer in the WSWP in contrast to the warming effect of the surface heat flux during the period. This cooling requires a source of colder water to the west. Similar analysis of OGCM outputs also suggests that the upwelled relatively cool water along the PNG north coast, and its northeastward extension to the equatorial region, contributes to cooling of the surface water over the WSWP mainly via negative zonal heat advection. Similar mechanisms are confirmed also for the 1982/83 and 1997/98 El Niño events by analyses of OGCM outputs and historical SST data. The low SST in the WSWP generated a positive zonal SST gradient together with high SST east of the WSWP. It may contribute to enhancement of the westerly surface wind in this region, leading to the onset of the 2002/03 El Niño event.     

Dislocation microstructures in naturally deformed silicate garnets

Dislocation microstructures of naturally deformed silicate garnets and olivines in garnet-peridotites and silicate garnets in eclogites from four localities have been observed with a transmission electron microscope (TEM) to clarify the dislocation characteristics of silicate garnets. We have obtained the following results: (1) dislocation densities of garnets in all the garnet-peridotites (ρ = 10⁵−10⁷ cm⁻²) are always nearly an order of magnitude lower than those of co-existing olivines; (2) dislocation densities of garnets in eclogites (ρ = 10⁵−10⁸ cm⁻²) which are embedded in garnet-peridotites are almost an order of magnitude higher than those of garnets in the surrounding garnet-peridotites; (3) the dominant Burgers vector, b, of mainly edge dislocations in garnet is 100 for specimens with dislocation density ρ = 10⁵−10⁶ cm⁻², while b=1/2111 for specimens with ρ = 10⁷−10⁸ cm⁻². Result (1) indicates that the observed dislocations in garnets were formed by plastic deformation under the same stresses as for co-existing olivines, and that there is a similar relationship between applied stress and dislocation density for garnets as for olivines. Result (2) suggests that the stress concentration occurred around eclogites embedded in garnet-peridotites, and the resulting differential stress in garnets in eclogites was further elevated by the interlocking of neighboring hard garnet grains. Finally, result (3) indicates that the dominant Burgers vector of mainly edge dislocations in garnet changes from 100 to 1/2111 with increasing applied differential stress.     

Monitoring of fault healing after the 1999 Kocaeli, Turkey, earthquake

The North Anatolian fault zone that ruptured during the mainshock of theM 7.4 Kocaeli (Izmit) earthquake of 17 August 1999 has beenmonitored using S wave splitting, in order to test a hypothesisproposed by Tadokoro et al. (1999). This idea is based on the observationof the M 7.2 1995 Hyogo-ken Nanbu (Kobe) earthquake, Japan.After the Hyogo-ken Nanbu earthquake, a temporal change was detectedin the direction of faster shear wave polarization in 2–3 years after the mainshock (Tadokoro, 1999). Four seismic stations were installed within andnear the fault zone at Kizanlik where the fault offset was 1.5 m, about80 km to the east of the epicenter of the Kocaeli earthquake. Theobservation period was from August 30 to October 27, 1999. Preliminaryresult shows that the average directions of faster shear wave polarization attwo stations were roughly parallel to the fault strike. We expect that thedirection of faster shear wave polarization will change to the same directionas the regional tectonic stress reflecting fault healing process. We havealready carried out a repeated aftershock observation at the same site in2000 for monitoring the fault healing process.     

ASCA Observation of the Quiescent X-Ray Counterpart to SGR 1627-41

We present a 2-10 keV ASCA observation of the field around the soft gamma repeater SGR 1627-41. A quiescent X-ray source, whose position is consistent both with that of a recently discovered BeppoSAX X-ray source and with the Interplanetary Network localization for this soft gamma repeater, was detected in this observation. In 2-10 keV X-rays, the spectrum of the X-ray source may be fit equally well by a power-law, blackbody, or bremsstrahlung function, with unabsorbed flux approximately 5x10-12 ergs cm-2 s-1. We do not confirm a continuation of a fading trend in the flux, and we find no evidence for periodicity, both of which were noted in the earlier BeppoSAX observations.     

Three-dimensional image of the anisotropic bodies beneath central Honshu, Japan

Shear-wave splitting from local deep earthquakes is investigated to clarify the volume and the location of two anisotropic bodies in the mantle wedge beneath central Honshu, Japan. We observe a spatial variation in splitting parameters depending on the combination of sources and receivers, nearly N–S fast in the northern region, nearly E–W fast in the southern region and small time delays in the eastern region. Using forward modelling, two models with 30 and 10 per cent anisotropy are tested by means of a global search for the locations of anisotropic bodies with various volumes. The optimum model is obtained for 30 per cent anisotropy, which means a 5 per cent velocity difference between fast and slow polarized waves. The northern anisotropic body has a volume of 1.00° (longitude) × 0.5° (latitude) × 75 km (depth), with the orientation of the symmetry axis being N20°E. The southern anisotropic body has a volume of 1.25° × 1.25° × 100 km with the symmetry axis along N95°E. Our results show that the anisotropic bodies are located in low-velocity and low- Q regions of the mantle. This, together with petrological data and the location of volcanoes in the arc, suggests that the possible cause of the anisotropy is the preferred alignment of cracks filled with melt.     

Source mechanisms and tectonic significance of historical earthquakes along the nankai trough,Japan

Two recent and three historical earthquakes which occurred along the Nankai trough, marking the northern plate boundary between the Philippine Sea and the Asian Plate, are studied mainly on the basis of the data of crustal deformations and tsunami waves. These earthquakes are the 1946 Nankaido, the 1944 Tonankai, the 1854 Ansei I, II and the 1707 Hoei earthquakes. They are all interpreted as low-angle thrust faults at the plate boundary, with the oceanic side underthrusting northwestward against southwestern Japan. The fault parameters of the historical earthquakes are assumed here to be common to those of the recent two earthquakes, except for the magnitude of dislocation.The entire fault region, which extends for 530 km from western Shikoku Island in the west to the Tokai district in the east, is divided into four fault planes, which are denoted the planes A, B, C and D, from west to east, respectively. Then, the five earthquakes may be attributed to the planes A, B, C and D, in the following manner: the Nankaido earthquake, A + B; the Tonankai earthquake, C; the Ansei II earthquake, A + B; the Ansei I earthquake, C + D; and the Hoei earthquake, A + B + C + D.The latest cycle of earthquake migration seems incomplete as proved by the recent inactivity in D. Consequently, the future major earthquake next to occur is expected there, off the Tokai district. Eight further ancient earthquakes from A.D. 684 to 1605 are also discussed. Taking the results of the foregoing studies into consideration, their sequence is well interpreted by the four migration cycles. Topographical data, tilt of coastal terraces and location of hinge lines, prove that the thrusting has continued all along the extension of the Nankai trough for at least 300,000 years.     

Simulation of groundwater dynamics in the North China Plain by coupled hydrology and agricultural models

We simulated the effects of irrigation on groundwater flow dynamics in the North China Plain by coupling the NIES Integrated Catchment‐based Ecohydrology (NICE) model with DSSAT‐wheat and DSSAT‐maize, two agricultural models. This combined model (NICE‐AGR) was applied to the Hai River catchment and the lower reach of the Yellow River (530 km wide by 840 km long) at a resolution of 5 km. It reproduced excellently the soil moisture, evapotranspiration and crop production of summer maize and winter wheat, correctly estimating crop water use. So, the spatial distribution of crop water use was reasonably estimated at daily steps in the simulation area. In particular, NICE‐AGR reproduced groundwater levels better than the use of statistical water use data. This indicates that NICE‐AGR does not need detailed statistical data on water use, making it very powerful for evaluating and estimating the water dynamics of catchments with little statistical data on seasonal water use. Furthermore, the simulation reproduced the spatial distribution of groundwater level in 1987 and 1988 in the Hebei Plain, showing a major reduction of groundwater level due mainly to overpumping for irrigation. Copyright © 2006 John Wiley & Sons, Ltd.     

Resonant excitation of the solarg-modes through coupling of 5-min oscillations

It is proposed that the interacting solarg-modes are driven by the beat forcing ofg-modes by two 5-min oscillations, excited by ak-mechanism. The observational implications of the proposed mechanism are discussed.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Estimation of groundwater use by crop production simulated by DSSAT‐wheat and DSSAT‐maize models in the piedmont region of the North China Plain

Sustainable water use is in serious crisis in the piedmont region of the Taihang Mountains in the North China Plain, owing to rapid groundwater drawdown. Estimating the water requirement for agriculture, the biggest user of groundwater, will be helpful in understanding groundwater decline. Through the use of DSSAT‐3·5 wheat and maize models, we assessed water use in winter wheat and maize, two staple crops in the region, in 1987–2001. Trends between groundwater change and simulated agricultural water use were compared. The results showed that groundwater decline was sensitive to simulated crop water requirement and irrigation requirement. According to regression analysis, 100 mm of water requirement by cultivated land (mainly wheat and maize) resulted in about 0·64 m of groundwater decline. This relationship might be useful in understanding the regional water balance and to help decision‐makers control groundwater decline through controlling crop water use or through long‐distance water transfer. The study demonstrated the usefulness of using the DSSAT model for estimating crop water use and the effectiveness of clarifying the reason for groundwater decline using the simulation results of water use. Copyright © 2006 John Wiley & Sons, Ltd.